Mark A. Christona,*, Ranjeet S. Patilb
aCCS-2 Continuum Dynamics, Los Alamos National Laboratory, Los Alamos, NM 87545, USA bMechanical Engineering Department, University of New Mexico, Albuquerque, NM 87131, USA
ABSTRACT
We extend the optimal second-order finite element projection method developed by Gresho et al.[1] in order to treat low-Mach number reacting flows using simple/reduced-order chemistry. The reacting flow formulation considers the zero Mach-number limit for combustion allowing for large heat release and the concomitant excursions in temperature and density while avoiding the need for a fully-compressible numerical treatment. In order to avoid spurious oscillations due to non-smooth data, a flux-corrected advection algorithm based on an explicit predictor/semi-implicit corrector and characteristic-based limiting is used. Results are presented for variable density non-reacting/reacting wakes and illustrate the process of bluff-body flame stabilization.
Keywords:
Low-Mach number; Reacting flow; Finite element; Projection method; CFD