Direct numerical simulation of a high-pressure turbulent reacting temporal mixing layer

Direct Numerical Simulation realizations were created of a temporal mixing layer in which combustion occurs under high-pressure (high-p) turbulent conditions. The model combines the formulation of Masi, et al. (2013) for describing multi-species mixing under high-p conditions and a single-step chemical reaction of rate consistent with ignition prediction (Borghesi, and Bellan, 2015). In each simulation the computations are pursued past a time at which a maximum average-volumetric p is attained; most analysis is performed at this time, t^*_(pp). The ensemble of realizations explores the effect of the initial Reynolds number, Re_0, of the initial pressure, p_0, and of the initial composition of the two mixing-layer streams. The results show that the thermodynamic energy added by the reaction at the small scales is partially dissipated and partially backscattered. The formation of turbulent small scales is initiated by the morphological changes in the flow through stretching and twisting rather than vice versa. The reaction establishes primarily in the oxidizer stream and is preponderantly of diffusion type. Overwhelmingly, the higher reaction rates occur in the diffusion flame, particularly in regions of high density-gradient magnitude. At higher p_0 the reaction rate reaches higher values and occurs in regions of higher density gradients. The range of reaction rates is independent of the Re_0 value but the magnitude of the density gradients increases with Re_0. When H_2O and CO_2 are initially present, uphill diffusion dominates over regular diffusion and occurs in regions of smaller density-gradient magnitude whereas regular diffusion occurs in regions of larger density-gradient magnitude where the reaction is more vigorous. H_2O is more prone than CO_2 to regular diffusion in the larger density-gradient magnitude regions. When H_2O and CO_2 only form in the flame, both H_2O and CO_2 are subject to regular diffusion over the entire range of density-gradient values. The dissipation probability density function is a log normal distribution at large dissipation values.

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

[2]  Ghislain Lartigue,et al.  Large-eddy simulation of supercritical fluid injection , 2013 .

[3]  R. Moser,et al.  Mixing transition and the cascade to small scales in a plane mixing layer , 1991 .

[4]  S. Sreedhara,et al.  Direct numerical simulation of autoignition in a non-premixed, turbulent medium , 2000 .

[5]  Imad A. Hannoun,et al.  Turbulence structure near a sharp density interface , 1988, Journal of Fluid Mechanics.

[6]  T. Takeno,et al.  A numerical study on flame stability at the transition point of jet diffusion flames , 1996 .

[7]  C. Law,et al.  Ignition of counterflowing methane versus heated air under reduced and elevated pressures , 1997 .

[8]  V. Giovangigli,et al.  Detailed modeling of planar transcritical H2–O2–N2 flames , 2011 .

[9]  Giulio Borghesi,et al.  Complex chemistry DNS of n-heptane spray autoignition at high pressure and intermediate temperature conditions , 2013 .

[10]  M. Pfitzner,et al.  LES of Low to High Turbulent Combustion in an Elevated Pressure Environment , 2014 .

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

[12]  D. Parsons,et al.  Self-consistent theory of the transition of an unstable thermodynamic system from spinodal to heterophase kinetics , 1991 .

[13]  T. Poinsot,et al.  Numerical simulations of autoignition in turbulent mixing flows , 1997 .

[14]  Stefan Hickel,et al.  Large-eddy simulation of nitrogen injection at trans- and supercritical conditions , 2016 .

[15]  Silvia M. Müller,et al.  A method to parallelize tridiagonal solvers , 1991, Parallel Comput..

[16]  F. Egolfopoulos,et al.  Non-premixed hydrocarbon ignition at high strain rates , 1998 .

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

[18]  G. Borghesi,et al.  Irreversible entropy production rate in high-pressure turbulent reactive flows , 2015 .

[19]  J. Bellan,et al.  Mixing rules for multicomponent mixture mass diffusion coefficients and thermal diffusion factors. , 2004, The Journal of chemical physics.

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

[21]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

[22]  S. Pope,et al.  Direct numerical simulation of a statistically stationary, turbulent reacting flow , 1999 .

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

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

[25]  Josette Bellan,et al.  Direct Numerical Simulations of O/H Temporal Mixing Layers Under Supercritical Conditions , 2002 .

[26]  J. Bellan,et al.  High-Pressure Binary Mass Diffusion Coefficients for Combustion Applications , 2004 .

[27]  P. Moin,et al.  DIRECT NUMERICAL SIMULATION: A Tool in Turbulence Research , 1998 .

[28]  E. Bruce Nauman,et al.  Nonlinear diffusion and phase separation , 2001 .

[29]  Alexandre Ern,et al.  Thermal diffusion effects in hydrogen-air and methane-air flames , 1998 .

[30]  Joel Keizer,et al.  Statistical Thermodynamics of Nonequilibrium Processes , 1987 .

[31]  Josette Bellan,et al.  Multi-species turbulent mixing under supercritical-pressure conditions: modelling, direct numerical simulation and analysis revealing species spinodal decomposition , 2013, Journal of Fluid Mechanics.

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

[33]  Jung Yul Yoo,et al.  The Turbulent Flows of Supercritical Fluids with Heat Transfer , 2013 .

[34]  E. Mastorakos Ignition of turbulent non-premixed flames , 2009 .

[35]  Alexandre Ern,et al.  Multicomponent transport algorithms , 1994 .

[36]  M. Carpenter,et al.  Several new numerical methods for compressible shear-layer simulations , 1994 .

[37]  G. Borghesi,et al.  A priori and a posteriori investigations for developing large eddy simulations of multi-species turbulent mixing under high-pressure conditions , 2015 .

[38]  C. Westbrook,et al.  Simplified Reaction Mechanisms for the Oxidation of Hydrocarbon Fuels in Flames , 1981 .

[39]  C. Cha,et al.  Direct numerical simulations of the double scalar mixing layer: Part II: Reactive scalars , 2007 .

[40]  H. Knapp,et al.  Vapor-liquid equilibria for mixtures of low-boiling substances , 1982 .

[41]  N. Easwar,et al.  Turbulent suppression of spinodal decomposition , 1984 .

[42]  Tim Edwards,et al.  CRACKING AND DEPOSITION BEHAVIOR OF SUPERCRITICAL HYDROCARBON AVIATION FUELS , 2006 .

[43]  Igor Pioro,et al.  Heat transfer to supercritical fluids flowing in channels—empirical correlations (survey) , 2004 .