An OpenFOAM pressure-based coupled CFD solver for turbulent and compressible flows in turbomachinery applications

ABSTRACT In this article, a recently developed pressure-based, fully coupled solver capable of predicting fluid flow at all speeds is extended to deal with turbulent flows in a rotating frame of reference and emphasizing turbomachinery applications. The pressure–velocity coupling at the heart of the Navier-Stokes equations is resolved by deriving a pressure equation in a similar fashion to a segregated SIMPLE algorithm but with implicit treatment of the velocity and pressure fields. The resulting system of coupled equations is solved using an algebraic multigrid solver. The above numerical procedures have been implemented within OpenFOAM®, which is an open-source code framework capable of dealing with industrial-scale flow problems. The OpenFOAM-based coupled solver is validated using experimental and numerical data available from reference literature test cases as well as with a segregated solver based on the SIMPLE algorithm. This is done in addition to evaluating its performance by solving an industrial problem. In comparison with the segregated solver, the coupled solver results indicate substantial reduction in computational cost with increased robustness.

[1]  Luca Mangani,et al.  Assessment of Various Turbulence Models in a High Pressure Ratio Centrifugal Compressor With an Object Oriented CFD Code , 2011 .

[2]  Feng Lin,et al.  Numerical Investigation on the Originating Mechanism of Unsteadiness in Tip Leakage Flow for a Transonic Fan Rotor , 2008 .

[3]  Sophie Papst,et al.  Computational Methods For Fluid Dynamics , 2016 .

[4]  M. Darwish,et al.  A high-resolution pressure-based algorithm for fluid flow at all speeds , 2001 .

[5]  Antony Jameson,et al.  Multigrid solution of the Euler equations using implicit schemes , 1985 .

[6]  D. Brandt,et al.  Multi-level adaptive solutions to boundary-value problems math comptr , 1977 .

[7]  William P. Timlake,et al.  An Accelerated Relaxation Algorithm for Iterative Solution of Elliptic Equations , 1968 .

[8]  M. Darwish,et al.  Pressure‐Based Algorithms for Single‐Fluid and Multifluid Flows , 2009 .

[9]  Luca Mangani,et al.  Conjugate Heat Transfer Analysis of NASA C3X Film Cooled Vane With an Object-Oriented CFD Code , 2010 .

[10]  M. Darwish,et al.  New family of adaptive very high resolution schemes , 1998 .

[11]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[12]  G. Raithby,et al.  A multigrid method based on the additive correction strategy , 1986 .

[13]  Susie Cristine Keller.  ADDITIVE CORRECTION MULTIGRID METHOD APPLIED TO DIFUSION PROBLEMS WITH UNSTRUCTURED GRIDS , 2004 .

[14]  S. Patankar,et al.  Pressure based calculation procedure for viscous flows at all speeds in arbitrary configurations , 1988 .

[15]  J. Spurk.  Boundary Layer Theory , 2019, Fluid Mechanics.

[16]  S. Muzaferija,et al.  Adaptive finite volume method for flow prediction using unstructured meshes and multigrid approach , 1994 .

[17]  C. Rhie,et al.  A numerical study of the flow past an isolated airfoil with separation , 1981 .

[18]  M. Darwish,et al.  A Fully Coupled Navier-Stokes Solver for Fluid Flow at All Speeds , 2014 .

[19]  M. Darwish, F. Moukalled, B. Sekar.  A UNIFIED FORMULATION OF THE SEGREGATED CLASS OF ALGORITHMS FOR MULTIFLUID FLOW AT ALL SPEEDS , 2001 .

[20]  M. Darwish,et al.  TVD schemes for unstructured grids , 2003 .

[21]  Michael D. Hathaway,et al.  Laser anemometer measurements in a transonic axial-flow fan rotor , 1989 .

[22]  V. Guinot.  Approximate Riemann Solvers , 2010 .

[23]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[24]  M. Darwish,et al.  A UNIFIED FORMULATION OF THE SEGREGATED CLASS OF ALGORITHMS FOR FLUID FLOW AT ALL SPEEDS , 2000 .

[25]  A. Gosman,et al.  High resolution NVD differencing scheme for arbitrarily unstructured meshes , 1999 .

[26]  Luca Mangani,et al.  DEVELOPMENT AND VALIDATION OF A C++ OBJECT ORIENTED CFD CODE FOR HEAT TRANSFER ANALYSIS , 2007 .

[27]  M. Darwish,et al.  NORMALIZED VARIABLE AND SPACE FORMULATION METHODOLOGY FOR HIGH-RESOLUTION SCHEMES , 1994 .

[28]  J. P. V. Doormaal,et al.  ENHANCEMENTS OF THE SIMPLE METHOD FOR PREDICTING INCOMPRESSIBLE FLUID FLOWS , 1984 .

[29]  I. Demirdzic,et al.  A COLLOCATED FINITE VOLUME METHOD FOR PREDICTING FLOWS AT ALL SPEEDS , 2005 .

[30]  M. Darwish,et al.  A COMPARATIVE ASSESSMENT OF THE PERFORMANCE OF MASS CONSERVATION-BASED ALGORITHMS FOR INCOMPRESSIBLE MULTIPHASE FLOWS , 2002 .

[31]  Luca Mangani,et al.  Investigation of Innovative Trailing Edge Cooling Configurations With Enlarged Pedestals and Square or Semicircular Ribs: Part 1—Experimental Results , 2008 .

[32]  M. Darwish,et al.  AN EFFICIENT VERY-HIGH-RESOLUTION SCHEME BASED ON AN ADAPTIVE-SCHEME STRATEGY , 1998 .

[33]  Luca Mangani,et al.  Application of an Object-Oriented CFD Code to Heat Transfer Analysis , 2008 .

[34]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[35]  Gene H. Golub,et al.  Numerical solution of saddle point problems , 2005, Acta Numerica.

[36]  A. Gosman,et al.  PREDICTION OF IMPELLER- INDUCED FLOW IN MIXING VESSELS USING MULTIPLE FRAMES OF REFERENCE , 1994 .

[37]  R. P. Fedorenko.  A relaxation method for solving elliptic difference equations , 1962 .

[38]  Luca Mangani,et al.  AN OBJECT-ORIENTED CFD CODE FOR OPTIMIZATION OF HIGH PRESSURE RATIO COMPRESSORS , 2012 .

[39]  M. Darwish,et al.  Development of a Novel Fully Coupled Solver in OpenFOAM: Steady-State Incompressible Turbulent Flows , 2014 .

[40]  R. Pletcher,et al.  On solving the compressible Navier-Stokes equations for unsteady flows at very low Mach numbers , 1993 .

[41]  C. Rhie,et al.  A numerical study of the turbulent flow past an isolated airfoil with trailing edge separation , 1982 .

[42]  Juan J. Alonso,et al.  Fluid/structure coupled aeroelastic computations for transonic flows in turbomachinery , 2002 .

[43]  Allan D. Grosvenor.  RANS Prediction of Transonic Compressive Rotor Performance Near Stall , 2007 .

[44]  E. Turkel,et al.  PRECONDITIONING TECHNIQUES IN COMPUTATIONAL FLUID DYNAMICS , 1999 .

[45]  George D. Raithby,et al.  Numerical solution of problems in incompressible fluid flow: treatment of the température-velocity coupling , 1986 .

[46]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[47]  M. S. Mihelc,et al.  Analytical and Experimental Evaluation of the Heat Transfer Distribution over the Surfaces of Turbine Vanes , 1983 .

[48]  M. Darwish,et al.  A Coupled Incompressible Flow Solver on Structured Grids , 2007 .

[49]  Katja Bachmeier,et al.  Numerical Heat Transfer And Fluid Flow , 2016 .

[50]  Luca Mangani,et al.  CONJUGATE HEAT TRANSFER ANALYSIS OF AN INTERNALLY COOLED TURBINE BLADES WITH AN OBJECT ORIENTED CFD CODE , 2009 .

[51]  M. Darwish,et al.  The χ-schemes: a new consistent high-resolution formulation based on the normalized variable methodology , 2003 .

[52]  F. Moukalled,et al.  A coupled finite volume solver for the solution of incompressible flows on unstructured grids , 2009, J. Comput. Phys..

[53]  Bruno Facchini,et al.  Conjugate Heat Transfer Simulation of a Radially Cooled Gas Turbine Vane , 2004 .

[54]  M. Darwish,et al.  Development of a Novel Fully Coupled Solver in OpenFOAM: Steady-State Incompressible Turbulent Flows in Rotational Reference Frames , 2014 .