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.
Copyright © 2006 by the American Institute of Aeronautics and Astronautics, Inc. Presented as Paper 1282 at the 43rd AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, 10–13 January 2005; received 15 June 2005; revision received 21 August 2006; accepted for publication 28 August 2006. This work was carried out while V. Raman was a Research Associate at the Center for Turbulence Research. The authors would like to thank Frank Ham, Magnus Svard, and Heinz Pitsch for fruitful discussions.
Published - HERaiaaj06.pdf