Large Eddy Simulation of Shear Coaxial Rocket Injector: Real Fluid Effects

The implementation and verification of real-fluid effects towards the high-fidelity large eddy simulation of rocket combustors is reported. The non-ideal fluid behavior is modeled using a cubic Peng-Robinson equation of state; a thermodynamically consistent approach is used to convert conserved into primitive variables. The viscosity is estimated by Chung et al.’s method in the supercritical gas phase. In the transcritical liquid phase, a simple, accurate and efficient method to estimate the viscosity as a function of temperature and pressure is proposed. The highly non-linear coupling of the primitive thermodynamic variables requires special consideration in regions of high-density gradients to avoid spurious numerical oscillations. The characterization of the non-linearity of the equation of state identifies the regions of high sensitivity. In these regions, small relative changes in the pressure lead to significant changes in density and/or temperature, therefore, numerical instabilities tend to be amplified in these regions. To avoid non-physical oscillations, a first-order and second-order essentially non-oscillatory (ENO) schemes are locally applied to the flux computation on the faces identified with a dual-threshold relative density sensor. The evaluation of the sensor and capabilities of the non-oscillatory schemes on canonical test cases are presented. Finally, these schemes are used to model two canonical cases.

[1]  C. Segal,et al.  Disintegrating supercritical jets in a subcritical environment , 2013, Journal of Fluid Mechanics.

[2]  Hiroshi Terashima,et al.  Approach for simulating gas-liquid-like flows under supercritical pressures using a high-order central differencing scheme , 2012, J. Comput. Phys..

[3]  Guilhem Lacaze,et al.  A non-premixed combustion model based on flame structure analysis at supercritical pressures , 2012 .

[4]  Sébastien Candel,et al.  Experiments and numerical simulation of mixing under supercritical conditions , 2012 .

[5]  Nobuhiro Yamanishi,et al.  High-Resolution Numerical Method for Supercritical Flows with Large Density Variations , 2011 .

[6]  Anthony Ruiz,et al.  NUMERICAL INVESTIGATION OF DESTABILIZATION OF SUPERCRITICAL ROUND TURBULENT JETS USING LARGE EDDY SIMULATION , 2011 .

[7]  Sébastien Candel,et al.  Large-Eddy Simulation of oxygen/methane flames under transcritical conditions , 2011 .

[8]  Oskar Haidn,et al.  Effect of Recess in High-Pressure Liquid Oxygen/Methane Coaxial Injection and Combustion , 2009 .

[9]  Laurent Selle,et al.  Modelling of subgrid-scale phenomena in supercritical transitional mixing layers: an a priori study , 2007, Journal of Fluid Mechanics.

[10]  M. Oschwald,et al.  Atomization and Flames in LOX/H2- and LOx/CH4- Spray Combustion , 2007 .

[11]  M. Habiballah,et al.  EXPERIMENTAL STUDIES OF HIGH-PRESSURE CRYOGENIC FLAMES ON THE MASCOTTE FACILITY , 2006 .

[12]  J. Bellan,et al.  Real-Gas Effects on Mean Flow and Temporal Stability of Binary-Species Mixing Layers , 2003 .

[13]  Claudio Bruno,et al.  Supercritical Combustion Properties , 2003 .

[14]  R. Branam,et al.  Length scales in cryogenic injection at supercritical pressure , 2002 .

[15]  Josette Bellan,et al.  Direct numerical simulation of a transitional supercritical binary mixing layer: heptane and nitrogen , 2002, Journal of Fluid Mechanics.

[16]  Josette Bellan,et al.  Consistent Boundary Conditions for Multicomponent Real Gas Mixtures Based on Characteristic Waves , 2002 .

[17]  Douglas G Talley,et al.  Visual characteristics and initial growth rates of round cryogenic jets at subcritical and supercritical pressures , 2002 .

[18]  Josette Bellan,et al.  Direct numerical simulations of supercritical fluid mixing layers applied to heptane–nitrogen , 2001, Journal of Fluid Mechanics.

[19]  Roland Span,et al.  A Reference Equation of State for the Thermodynamic Properties of Nitrogen for Temperatures from 63.151 to 1000 K and Pressures to 2200 MPa , 2000 .

[20]  J. Vera Molecular thermodynamics of fluid‐phase equilibria by John M. Prausnitz, Rüdiger N. Lichtenthaler and Edmundo Comes de Azevedo, Third Edition, 1999; Prentice Hall PTR, Upper‐Saddle River, New Jersey 07458, xxiii + 860 pages; price $131.95 CND; ISBN 0‐13‐977745‐8 , 2000 .

[21]  Michael Oschwald,et al.  Supercritical nitrogen free jet investigated by spontaneous Raman scattering , 1999 .

[22]  Joseph C. Oefelein,et al.  Modeling High-Pressure Mixing and Combustion Processes in Liquid Rocket Engines , 1998 .

[23]  W. Mayer,et al.  Atomization and Breakup of Cryogenic Propellants Under High-Pressure Subcritical and Supercritical Conditions , 1998 .

[24]  Sébastien Candel,et al.  Experimental investigation of shear coaxial cryogenic jet flames , 1998 .

[25]  Josette Bellan,et al.  Efficient High-Pressure State Equations , 1997 .

[26]  Hiroshi Tamura,et al.  Flow Visualization of Supercritical Propellant Injection in a Firing LOX/GH2 Rocket Engine , 1995 .

[27]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[28]  K. E. Starling,et al.  Generalized multiparameter correlation for nonpolar and polar fluid transport properties , 1988 .

[29]  K. E. Starling,et al.  Applications of kinetic gas theories and multiparameter correlation for prediction of dilute gas viscosity and thermal conductivity , 1984 .

[30]  D. Peng,et al.  A New Two-Constant Equation of State , 1976 .

[31]  E. Andrade,et al.  LVIII. A theory of the viscosity of liquids.—Part II , 1934 .

[32]  C. Raman,et al.  A Theory of the Viscosity of Liquids , 1923, Nature.