A methodology for high performance computation of fully inhomogeneous turbulent flows

A large-eddy simulation methodology for high performance parallel computation of statistically fully inhomogeneous turbulent flows on structured grids is presented. Strategies and algorithms to improve the memory efficiency as well as the parallel performance of the subgrid-scale model, the factored scheme, and the Poisson solver on shared-memory parallel platforms are proposed and evaluated. A novel combination of one-dimensional red–black/line Gauss–Seidel and two-dimensional red–black/line Gauss–Seidel methods is shown to provide high efficiency and performance for multigrid relaxation of the Poisson equation. Parallel speedups are measured on various shared-distributed memory systems. Validations of the code are performed in large-eddy simulations of turbulent flows through a straight channel and a square duct. Results obtained from the present solver employing a Lagrangian dynamic subgrid-scale model show good agreements with other available data. The capability of the code for more complex flows is assessed by performing a large-eddy simulation of the tip-leakage flow in a linear cascade. Copyright © 2006 John Wiley & Sons, Ltd.

[1]  Trong T. Bui,et al.  A parallel finite volume algorithm for large-eddy simulation of turbulent flows , 1998 .

[2]  F. Durst,et al.  Computational fluid dynamics applications on parallel-vector computers: computations of stirred vessel flows , 2002 .

[3]  P. Moin,et al.  A dynamic localization model for large-eddy simulation of turbulent flows , 1995, Journal of Fluid Mechanics.

[4]  William L. Briggs,et al.  A multigrid tutorial, Second Edition , 2000 .

[5]  William J. Devenport,et al.  Wake of a Compressor Cascade with Tip Gap, Part 1: Mean Flow and Turbulence Structure , 2004 .

[6]  Donghyun You.  Study of tip clearance flow in a turbomachinery cascade using large eddy simulation , 2004 .

[7]  Barbara M. Chapman,et al.  Performance Oriented Programming for NUMA Architectures , 2001, WOMPAT.

[8]  Matthias S. Müller.  Some Simple OpenMP Optimization Techniques , 2001, WOMPAT.

[9]  Haecheon Choi,et al.  Turbulent Drag Reduction: Studies of Feedback Control and Flow Over Riblets , 1993 .

[10]  Meng Wang,et al.  A Computational Methodology for Large-Eddy Simulation of Tip-Clearance Flows , 2003 .

[11]  Jun Zhang.  Multigrid Method and Fourth-Order Compact Scheme for 2D Poisson Equation with Unequal Mesh-Size Discretization , 2002 .

[12]  S. Menon,et al.  AN UNSTEADY INCOMPRESSIBLE NAVIER-STOKES SOLVER FOR LARGE EDDY SIMULATION OF TURBULENT FLOWS , 1999 .

[13]  Pierluigi Amodio,et al.  A Parallel Gauss-Seidel Method for Block Tridiagonal Linear Systems , 1995, SIAM J. Sci. Comput..

[14]  S. McCormick,et al.  A multigrid tutorial (2nd ed.) , 2000 .

[15]  Peter S. Pacheco.  Parallel programming with MPI , 1996 .

[16]  P. Moin,et al.  Suitability of upwind-biased finite difference schemes for large-eddy simulation of turbulent flows , 1997 .

[17]  Glenn R. Luecke,et al.  Scalability and performance of OpenMP and MPI on a 128‐processor SGI Origin 2000 , 2001, Concurr. Comput. Pract. Exp..

[18]  H. F. Jordan,et al.  Is SOR Color-Blind? , 1986 .

[19]  Jun Zhangy.  Multigrid Method and Fourth Order Compact Diierence Scheme for 2d Poisson Equation with Unequal Meshsize Discretization , 2001 .

[20]  Surya Pratap Vanka,et al.  Large eddy simulation of turbulence‐driven secondary flow in a square duct , 1991 .

[21]  Pierluigi Amodio,et al.  The conditioning of Toeplitz band matrices , 1996 .

[22]  Asmund Huser,et al.  Direct numerical simulation of turbulent flow in a square duct , 1992 .

[23]  T. Lund,et al.  Generation of Turbulent Inflow Data for Spatially-Developing Boundary Layer Simulations , 1998 .

[24]  P. Moin,et al.  Turbulence statistics in fully developed channel flow at low Reynolds number , 1987, Journal of Fluid Mechanics.

[25]  P. Moin,et al.  A dynamic subgrid‐scale eddy viscosity model , 1990 .

[26]  Shanhong Ji,et al.  Flutter Computation of Turbomachinery Cascades Using a Parallel Unsteady Navier-Stokes Code , 1999 .

[27]  Pierluigi Amodio,et al.  Parallel iterative solvers for boundary value methods , 1996 .

[28]  Bengt Sundén,et al.  Parallelization of a multi-blocked CFD code via three strategies for fluid flow and heat transfer analysis , 2004 .

[29]  Raymond T. Pierrehumbert,et al.  Surface quasi-geostrophic dynamics , 1995, Journal of Fluid Mechanics.

[30]  Yu Wang,et al.  Wake of a Compressor Cascade with Tip Gap, Part 2: Effects of Endwall Motion , 2004 .

[31]  David J. Evans,et al.  Parallel S.O.R. iterative methods , 1984, Parallel Comput..

[32]  William H. Press,et al.  Numerical recipes , 1990 .

[33]  S. Gavrilakis,et al.  Numerical simulation of low-Reynolds-number turbulent flow through a straight square duct , 1992, Journal of Fluid Mechanics.

[34]  P. Moin,et al.  Numerical experiments on the flow past A circular cylinder at sub-critical reynolds number , 1994 .

[35]  Jay Hoeflinger,et al.  Producing scalable performance with OpenMP: Experiments with two CFD applications , 2001, Parallel Comput..

[36]  Ugo Piomelli,et al.  Large-eddy simulation of rotating channel flows using a localized dynamic model , 1995 .

[37]  C. Meneveau,et al.  A Lagrangian dynamic subgrid-scale model of turbulence , 1994, Journal of Fluid Mechanics.

[38]  Marco Zagha,et al.  OriginTM 2000 and Onyx2® Performance Tuning and Optimization Guide , 1993 .