Abstract
Preserving scalar boundedness is an important prerequisite to performing large-eddy simulations of turbulent reacting flows. A number of popular combustion models use a conserved-scalar, mixture-fraction to parameterize reactions that, by definition, is bound between zero and one. To avoid unphysical clipping, the numerical scheme solving the conserved-scalar transport equation must preserve these bounds, while minimizing the amount of numerical diffusivity. To this end, a flux correction method is presented and applied to the quadratic-upwind biased interpolative convective scheme that ensures preservation of the scalar's physical bounds while retaining the low numerical diffusivity of the original quadratic-upwind biased interpolative convective scheme. It is demonstrated that this bounded quadratic-upwind biased interpolative convective scheme outperforms the third-order weighted essentially nonoscillatory scheme in maintaining spatial accuracy and reducing numerical dissipation errors both in generic test cases as well as direct numerical simulation of canonical flows.
Original language | English (US) |
---|---|
Pages (from-to) | 2879-2886 |
Number of pages | 8 |
Journal | AIAA journal |
Volume | 44 |
Issue number | 12 |
DOIs | |
State | Published - Dec 2006 |
Externally published | Yes |
ASJC Scopus subject areas
- Aerospace Engineering