VODE: a variable-coefficient ODE solver
SIAM Journal on Scientific and Statistical Computing
Efficient implementation of weighted ENO schemes
Journal of Computational Physics
Direct numerical simulation of a turbulent reactive plume on a parallel computer
Journal of Computational Physics
A semi-implicit numerical scheme for reacting flow: I. stiff chemistry
Journal of Computational Physics
A semi-implicit numerical scheme for reacting flow: II. stiff, operator-split formulation
Journal of Computational Physics
Conservative high-order finite-difference schemes for low-Mach number flows
Journal of Computational Physics
An analysis of operator splitting techniques in the stiff case
Journal of Computational Physics
Additive Runge-Kutta schemes for convection-diffusion-reaction equations
Applied Numerical Mathematics
Adaptive low Mach number simulations of nuclear flame microphysics
Journal of Computational Physics
Fourth-Order Time-Stepping for Stiff PDEs
SIAM Journal on Scientific Computing
Modeling Low Mach Number Reacting Flow with Detailed Chemistry and Transport
Journal of Scientific Computing
Time-accurate calculation of variable density flows with strong temperature gradients and combustion
Journal of Computational Physics
A high-order public domain code for direct numerical simulations of turbulent combustion
Journal of Computational Physics
Hi-index | 31.45 |
We present an improved numerical scheme for numerical simulations of low Mach number turbulent reacting flows with detailed chemistry and transport. The method is based on a semi-implicit operator-splitting scheme with a stiff solver for integration of the chemical kinetic rates, developed by Knio et al. [O.M. Knio, H.N. Najm, P.S. Wyckoff, A semi-implicit numerical scheme for reacting flow II. Stiff, operator-split formulation, Journal of Computational Physics 154 (2) (1999) 428-467]. Using the material derivative form of continuity equation, we enhance the scheme to allow for large density ratio in the flow field. The scheme is developed for direct numerical simulation of turbulent reacting flow by employing high-order discretization for the spatial terms. The accuracy of the scheme in space and time is verified by examining the grid/time-step dependency on one-dimensional benchmark cases: a freely propagating premixed flame in an open environment and in an enclosure related to spark-ignition engines. The scheme is then examined in simulations of a two-dimensional laminar flame/vortex-pair interaction. Furthermore, we apply the scheme to direct numerical simulation of a homogeneous charge compression ignition (HCCI) process in an enclosure studied previously in the literature. Satisfactory agreement is found in terms of the overall ignition behavior, local reaction zone structures and statistical quantities. Finally, the scheme is used to study the development of intrinsic flame instabilities in a lean H"2/air premixed flame, where it is shown that the spatial and temporary accuracies of numerical schemes can have great impact on the prediction of the sensitive nonlinear evolution process of flame instability.