An easily implemented task-based parallel scheme for the Fourier pseudospectral solver applied to 2D Navier–Stokes turbulence

Abstract An efficient parallel scheme is proposed for performing direct numerical simulation (DNS) of two-dimensional Navier–Stokes turbulence at high Reynolds numbers. We illustrate the resulting numerical code by displaying relaxation to states close to those that have been predicted by statistical–mechanical methods which start from ideal (Euler) fluid mechanics. The validation of these predictions by DNS requires unusually long computation times on single-cpu workstations, and suggests the use of parallel computation. The performance of our MPI Fortran 90 code on the SGI Origin 3800 is reported, together with its comparison with another parallel method. A few computational results that illustrate tests of the statistical–mechanical predictions are presented.

[1]  Glenn Joyce,et al.  Negative temperature states for the two-dimensional guiding-centre plasma , 1973, Journal of Plasma Physics.

[2]  Shiyi Chen,et al.  On statistical correlations between velocity increments and locally averaged dissipation in homogeneous turbulence , 1993 .

[3]  James C. McWilliams,et al.  The emergence of isolated coherent vortices in turbulent flow , 1984, Journal of Fluid Mechanics.

[4]  G. Eyink,et al.  Negative-temperature states and large-scale, long-lived vortices in two-dimensional turbulence , 1993 .

[5]  C. Y. Chu,et al.  Comparison of two-dimensional FFT methods on the hypercube , 1989, C3P.

[6]  T. A. Zang,et al.  Spectral methods for fluid dynamics , 1987 .

[7]  R. Robert,et al.  Statistical equilibrium states for two-dimensional flows , 1991, Journal of Fluid Mechanics.

[8]  T. Lundgren,et al.  Statistical mechanics of two‐dimensional vortices in a bounded container , 1976 .

[9]  Javier Jiménez,et al.  The structure of intense vorticity in isotropic turbulence , 1993, Journal of Fluid Mechanics.

[10]  Alternative statistical-mechanical descriptions of decaying two-dimensional turbulence in terms of ``patches'' and ``points'' , 2002, physics/0211024.

[11]  Steven A. Orszag,et al.  A case study in parallel computing: I. Homogeneous turbulence on a hypercube , 1991 .

[12]  M. Cross,et al.  Statistical Mechanics, Euler’s Equation, and Jupiter’s Red Spot , 1992 .

[13]  Evangelos A. Coutsias,et al.  Pseudospectral Solution of the Two-Dimensional Navier-Stokes Equations in a Disk , 1999, SIAM J. Sci. Comput..

[14]  R. Pelz.  The parallel Fourier pseudospectral method , 1991 .

[15]  Robert B. Ross,et al.  Using MPI-2: Advanced Features of the Message Passing Interface , 2003, CLUSTER.

[16]  J. Chasnov.  On the decay of two‐dimensional homogeneous turbulence , 1997 .

[17]  Sean Oughton,et al.  Relaxation in two dimensions and the "sinh-Poisson" equation , 1992 .

[18]  R. Peyret.  Spectral Methods for Incompressible Viscous Flow , 2002 .

[19]  S. Gauthier,et al.  2D pseudo-spectral parallel Navier–Stokes simulations of compressible Rayleigh–Taylor instability , 2002 .

[20]  Anthony Skjellum,et al.  Using MPI - portable parallel programming with the message-parsing interface , 1994 .

[21]  S. Kida,et al.  Late states of incompressible 2D decaying vorticity fields , 1997, chao-dyn/9709020.

[22]  Anthony Skjellum,et al.  A High-Performance, Portable Implementation of the MPI Message Passing Interface Standard , 1996, Parallel Comput..

[23]  D. Martínez,et al.  Decaying, two-dimensional, Navier-Stokes turbulence at very long times , 1991 .