TY - JOUR

T1 - A new parallel strategy for two-dimensional incompressible flow simulations using pseudo-spectral methods

AU - Yin, Z.

AU - Yuan, Li

AU - TANG, Tao

N1 - Funding Information:
We thank Linbo Zhang, Zhongze Li, and Ying Bai for the support of using local parallel computers. Z.Y. also thanks Professor W.H. Matthaeus who supplied the original serial FORTRAN 77 Navier–Stokes pseudo-spectral codes. This work is supported by National Natural Science Foundation of China (G10476032). T.T. thanks the supports from International Research Team on Complex System, Chinese Academy of Sciences, and from Hong Kong Research Grant Council.

PY - 2005/11/20

Y1 - 2005/11/20

N2 - A novel parallel technique for Fourier-Galerkin pseudo-spectral methods with applications to two-dimensional Navier-Stokes equations and inviscid Boussinesq approximation equations is presented. It takes advantage of the programming structure of the phase-shift de-aliased scheme for pseudo-spectral codes, and combines the task-distribution strategy [Z. Yin, H.J.H. Clercx, D.C. Montgomery, An easily implemented task-based parallel scheme for the Fourier pseudo-spectral solver applied to 2D Navier-Stokes turbulence, Comput. Fluid 33 (2004) 509] and parallelized Fast Fourier Transform scheme. The performances of the resulting MPI Fortran90 codes with the new procedure on SGI 3800 are reported. For fixed resolution of the same problem, the peak speed of the new scheme can be twice as fast as the old parallel methods. The parallelized codes are used to solve some challenging numerical problems governed by the Navier-Stokes equations and the Boussinesq equations. Two interesting physical problems, namely, the double-valued ω-ψ structure in two-dimensional decaying turbulence and the collapse of the bubble cap in the Boussinesq simulation, are solved by using the proposed parallel algorithms.

AB - A novel parallel technique for Fourier-Galerkin pseudo-spectral methods with applications to two-dimensional Navier-Stokes equations and inviscid Boussinesq approximation equations is presented. It takes advantage of the programming structure of the phase-shift de-aliased scheme for pseudo-spectral codes, and combines the task-distribution strategy [Z. Yin, H.J.H. Clercx, D.C. Montgomery, An easily implemented task-based parallel scheme for the Fourier pseudo-spectral solver applied to 2D Navier-Stokes turbulence, Comput. Fluid 33 (2004) 509] and parallelized Fast Fourier Transform scheme. The performances of the resulting MPI Fortran90 codes with the new procedure on SGI 3800 are reported. For fixed resolution of the same problem, the peak speed of the new scheme can be twice as fast as the old parallel methods. The parallelized codes are used to solve some challenging numerical problems governed by the Navier-Stokes equations and the Boussinesq equations. Two interesting physical problems, namely, the double-valued ω-ψ structure in two-dimensional decaying turbulence and the collapse of the bubble cap in the Boussinesq simulation, are solved by using the proposed parallel algorithms.

KW - Boussinesq equations

KW - Navier-Stokes equations

KW - Parallel computing

KW - Pseudo-spectral methods

KW - Task distribution

UR - http://www.scopus.com/inward/record.url?scp=23044461426&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2005.04.010

DO - 10.1016/j.jcp.2005.04.010

M3 - Journal article

AN - SCOPUS:23044461426

SN - 0021-9991

VL - 210

SP - 325

EP - 341

JO - Journal of Computational Physics

JF - Journal of Computational Physics

IS - 1

ER -