Approach for simulating gas-liquid-like flows under supercritical pressures using a high-order central differencing scheme

This study proposes an approach for simulations of cryogenic fluid mixing under supercritical pressures using high-order schemes. In this approach, we introduce a pressure evolution equation and consistently construct numerical diffusion terms to maintain the velocity and pressure equilibriums at fluid interfaces. The interfaces with high density and temperature ratio are successfully captured without the generation of spurious oscillations, while a high-order central differencing scheme resolves the flow fields. The present method preserves the mass and momentum conservation properties, while the poor energy conservation property is recognized. The one-dimensional single and multi-species interface advection and two-dimensional cryogenic jet mixing problems demonstrate the superiority and robustness of the present method over a conventional fully conservative method.

[1]  Smadar Karni,et al.  Compressible bubbles with surface tension , 1998 .

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

[3]  Keh-Ming Shyue,et al.  Regular Article: A Fluid-Mixture Type Algorithm for Compressible Multicomponent Flow with van der Waals Equation of State , 1999 .

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

[5]  Vigor Yang,et al.  Liquid rocket thrust chambers : aspects of modeling, analysis, and design , 2004 .

[6]  Ronald Fedkiw,et al.  A General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations , 2000 .

[7]  Smadar Karni,et al.  Hybrid Multifluid Algorithms , 1996, SIAM J. Sci. Comput..

[8]  D. Gaitonde,et al.  Pade-Type Higher-Order Boundary Filters for the Navier-Stokes Equations , 2000 .

[9]  Nicolas Guezennec,et al.  Numerical Studies of Mixing and Flame-Turbulence Interactions in Shear Coaxial Injector Flows under Trans-Critical Conditions , 2012 .

[10]  Andrew W. Cook,et al.  Artificial Fluid Properties for Large-Eddy Simulation of Compressible Turbulent Mixing , 2007 .

[11]  Rémi Abgrall,et al.  Computations of compressible multifluids , 2001 .

[12]  R. Abgrall.  How to Prevent Pressure Oscillations in Multicomponent Flow Calculations , 1996 .

[13]  Vigor Yang,et al.  A numerical study of cryogenic fluid injection and mixing under supercritical conditions , 2004 .

[14]  M. McLinden,et al.  NIST Standard Reference Database 23: Reference Fluid Thermodynamic and Transport Properties-REFPROP, Version 8.0 , 2007 .

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

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

[17]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[18]  Smadar Karni,et al.  Multicomponent Flow Calculations by a Consistent Primitive Algorithm , 1994 .

[19]  Keh-Ming Shyue,et al.  A fluid-mixture type algorithm for compressible multicomponent flow with Mie-Grüneisen equation of state , 2001 .

[20]  K. Thompson.  Time-dependent boundary conditions for hyperbolic systems, II , 1990 .

[21]  G. Soave.  Equilibrium constants from a modified Redlich-Kwong equation of state , 1972 .

[22]  Soshi Kawai,et al.  A high‐resolution scheme for compressible multicomponent flows with shock waves , 2011 .

[23]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[24]  Rémi Abgrall,et al.  An adaptive shock-capturing algorithm for solving unsteady reactive flows , 2003 .

[25]  Vigor Yang,et al.  A unified treatment of general fluid thermodynamics and its application to a preconditioning scheme , 2003 .

[26]  Soshi Kawai,et al.  Localized artificial diffusivity scheme for discontinuity capturing on curvilinear meshes , 2008, J. Comput. Phys..

[27]  Takashi Yabe,et al.  Unified Numerical Procedure for Compressible and Incompressible Fluid , 1991 .

[28]  Anthony Ruiz,et al.  Large-Eddy Simulation of Supercritical-Pressure Round Jets , 2010 .

[29]  W.,et al.  Time Dependent Boundary Conditions for Hyperbolic Systems , 2003 .

[30]  W. Cabot,et al.  A high-wavenumber viscosity for high-resolution numerical methods , 2004 .

[31]  C. Whitson,et al.  Estimating Diffusion Coefficients of Dense Fluids , 1993 .

[32]  Josette Bellan,et al.  Supercritical (and subcritical) fluid behavior and modeling: drops, streams, shear and mixing layers, jets and sprays , 2000 .

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