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.
ASJC Scopus subject areas
- Aerospace Engineering