Flux corrected finite volume scheme for preserving scalar boundedness in reacting large-eddy simulations

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.