Third M.I.T. Conference on Computational Fluid and Solid Mechanics June 14–17, 2005  

A finite element projection method for low-Mach number reacting flows

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

  Full Text
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

* Corresponding author. Tel.: +1 (505) 667 1540; Fax: +1 (505) 665 4972; E-mail: