Development and validation of a massively parallel flow solver for turbomachinery flows

The development and validation of the unsteady, three-dimensional, multiblock, parallel turbomachinery  ow solver TFLO is presented. The unsteady Reynolds-averaged Navier–Stokes equations are solved using a cellcentered discretization on arbitrary multiblock meshes. The solution procedure is based on efŽ cient explicit Runge– Kutta methods with several convergenceacceleration techniques such as multigrid, implicit residual smoothing, and local time stepping. The solver is parallelized using domain decomposition, a single program multiple data strategy, and the message passing interface standard. Details of the communication scheme and load balancing algorithms are discussed. A general and efŽ cient procedure for parallel interblade row interfacing is developed. The dual-time stepping technique is used to advance unsteady computations in time. The focus is on improving the parallel efŽ ciency and scalability of the  ow solver, as well as on its initial validation of steady-state calculations in multiblade row environment. The result of this careful implementation is a solver with demonstrated scalability up to 1024 processors. For validation and veriŽ cation purposes, results from TFLO are compared with both existing experimental data and computational results from other computational  uid dynamics codes used in aircraft engine industry.

[1]  R. A. Delaney,et al.  Numerical Prediction of Turbine Vane-Blade Aerodynamic Interaction , 1989 .

[2]  Mark G. Turner,et al.  An Average Passage Closure Model for General Meshes , 1999 .

[3]  Juan J. Alonso,et al.  Multigrid unsteady Navier-Stokes calculations with aeroelastic applications , 1995 .

[4]  D. A. Johnson,et al.  Transonic, turbulent boundary-layer separation generated on an axisymmetric flow model , 1986 .

[5]  U. K. Singh,et al.  Time marching methods for turbomachinery flow calculation , 1979 .

[6]  Ron-Ho Ni,et al.  Prediction of 3D multi-stage turbine flow field using a multiple-grid Euler solver , 1989 .

[7]  Antony Jameson,et al.  A new high resolution scheme for compressible viscous flows with shocks , 1995 .

[8]  Antony Jameson,et al.  Three-dimensional computations of time-dependent incompressible flows with an implicit multigrid-driven algorithm on parallel computers , 1997 .

[9]  Feng Liu,et al.  A Strongly Coupled Time-Marching Method for Solving the Navier—Stokes andk-ω Turbulence Model Equations with Multigrid , 1996 .

[10]  Edward J. Hall.  Aerodynamic Modeling of Multistage Compressor Flowfields: Part 2 — Modeling Deterministic Stresses , 1997 .

[11]  John J. Adamczyk,et al.  Simulation of three-dimensional viscous flow within a multistage turbine , 1990 .

[12]  Luigi Martinelli,et al.  Large eddy simulation of compressible flow using a parallel, multigrid driven algorithm , 1996 .

[13]  Antony Jameson,et al.  A new implicit algorithm with multigrid for unsteady incompressible flow calculations , 1995 .

[14]  Man Mohan Rai,et al.  Navier-Stokes Simulations of Rotor/Stator Interaction Using Patched and Overlaid Grids , 1987 .

[15]  R. M. Zacharias,et al.  1995 ASME Gas Turbine Award Paper: Development and Application of a Multistage Navier–Stokes Solver: Part I—Multistage Modeling Using Bodyforces and Deterministic Stresses , 1998 .

[16]  M. Giles,et al.  UNSFLO : a numerical method for unsteady inviscid flow in turbomachinery , 1988 .

[17]  Mark D. Sanetrik,et al.  Time-accurate Navier-Stokes calculations with multigrid acceleration , 1993 .

[18]  J. Adamczyk.  Model equation for simulating flows in multistage turbomachinery , 1996 .

[19]  Juan J. Alonso,et al.  A Two-Dimensional Multigrid-Driven Navier-Stokes Solver for Multiprocessor Architectures , 2000 .

[20]  Andrea Arnone,et al.  Rotor-Stator Interaction Analysis Using the Navier-Stokes Equations and a Multigrid Method , 1995 .

[21]  Daniel J. Dorney,et al.  Unsteady analysis of hot streak migration in a turbine stage , 1990 .

[22]  Craig J. Fischberg,et al.  Development and Application of a Multistage Navier-Stokes Solver: Part I — Multistage Modeling Using Bodyforces and Deterministic Stresses , 1995 .

[23]  Juan Jose Alonso.  Parallel computation of unsteady and aeroelastic flows using an implicit multigrid-driven algorithm , 1997 .

[24]  Michael B. Giles,et al.  Stator/rotor interaction in a transonic turbine , 1988 .

[25]  C. R. Lejambre,et al.  Development and Application of a Multistage Navier-Stokes Flow Solver: Part II — Application to a High Pressure Compressor Design , 1995 .

[26]  Juan J. Alonso,et al.  Development and Validation of a Massively Parallel Flow Solver for Turbomachinery Flows , 2001 .

[27]  William N. Dawes.  Simulation of unsteady blade row interaction with CFD: applications , 1998 .

[28]  Karl Engel,et al.  Numerical Investigation of the Clocking Effects in a Multistage Turbine , 1996 .

[29]  Olivier Sgarzi,et al.  A Parallel Computation of 3D Unsteady Flows in a Full Stage of Transonic Turbine. The Tumult Project , 1999 .

[30]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[31]  A. Jameson.  Time dependent calculations using multigrid, with applications to unsteady flows past airfoils and wings , 1991 .

[32]  John C. Vassberg,et al.  An Efficient Multiblock Method for Aerodynamic Analysis and Design on Distributed Memory Systems , 1997 .

[33]  Ron-Ho Ni,et al.  Using 3D Euler flow simulations to assess effects of periodic unsteady flow through turbines , 1990 .

[34]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[35]  Philip C. E. Jorgenson,et al.  Explicit Runge-Kutta method for unsteady rotor/stator interaction , 1988 .

[36]  Paul G. A. Cizmas,et al.  Parallel Computation of Turbine Blade Clocking , 1999 .

[37]  John J. Adamczyk,et al.  A Model for Closing the Inviscid Form of the Average-Passage Equation System , 1986 .

[38]  Paul G. A. Cizmas,et al.  Parallel Computation of Rotor-Stator Interaction , 1998 .

[39]  E J Hall.  Aerodynamic modelling of multistage compressor flow fields Part 2: Modelling deterministic stresses , 1998 .

[40]  R. Delaney,et al.  Investigation of unsteady flow through transonic turbine stage. I - Analysis , 1990 .

[41]  Feng Liu,et al.  Computation of turbomachinery flow by a convective-upwind-split-pressure (CUSP) scheme , 1998 .

[42]  Roger L. Davis,et al.  Prediction of 3-D unsteady flow in multi-stage turbomachinery using an implicit dual time-step approach , 1996 .

[43]  A. Jameson.  ANALYSIS AND DESIGN OF NUMERICAL SCHEMES FOR GAS DYNAMICS, 1: ARTIFICIAL DIFFUSION, UPWIND BIASING, LIMITERS AND THEIR EFFECT ON ACCURACY AND MULTIGRID CONVERGENCE , 1995 .

[44]  M. L. Celestina,et al.  A model for closing the inviscid form of the average−passage equation system , 1986 .

[45]  Karen Gundy-Burlet.  Computations of unsteady multistage compressor flows in a workstation environment , 1991 .

[46]  A. Jameson.  ANALYSIS AND DESIGN OF NUMERICAL SCHEMES FOR GAS DYNAMICS, 2: ARTIFICIAL DIFFUSION AND DISCRETE SHOCK STRUCTURE , 1994 .

[47]  E J Hall.  Aerodynamic modelling of multistage compressor flow fields Part 1: Analysis of rotor-stator-rotor aerodynamic interaction , 1998 .

[48]  Meng-Sing Liou,et al.  Multigrid Time-Accurate Integrations of Navier-Stokes Equations , 1993 .

[49]  Edward J. Hall.  Aerodynamic Modeling of Multistage Compressor Flowfields: Part 1 — Analysis of Rotor/Stator/Rotor Aerodynamic Interaction , 1997 .

[50]  Juan J. Alonso,et al.  Automatic aerodynamic optimization on distributed memory architectures , 1996 .