Reaction analogy based forcing for incompressible scalar turbulence

Abstract

We present a novel reaction analogy (RA) based forcing method for generating statistically stationary scalar fields in incompressible turbulence. The new method can produce more general scalar probability density functions (PDFs), for example, quasi-double-delta PDF, than current methods, while ensuring that scalar fields remain bounded, unlike existent forcing methodologies that can potentially violate naturally existing bounds. Such features are useful for generating initial fields in nonpremixed combustion, inlet conditions for spatially developing flows, or for studying non-Gaussian scalar turbulence. The RA method mathematically models hypothetical chemical reactions that convert reactants in a mixed state back into its pure unmixed components. Various types of chemical reactions are formulated and the corresponding mathematical expressions derived such that the reaction term is smooth in scalar space and is consistent with mass conservation. For large values of the scalar forcing rate, the method produces statistically stationary quasi-double-d scalar PDFs. Quasiuniform, Gaussian, and stretched exponential scalar statistics are recovered for smaller values of the scalar forcing rate. The shape of the scalar PDF can be further controlled by changing the stoichiometric coefficients of the reaction. The ability of the new method to produce fully developed passive scalar fields with quasi-Gaussian PDFs is also investigated, by exploring the convergence of the scalar variance spectrum to the Obukhov-Corrsin scaling and of the third-order mixed structure function to the "four-thirds" Yaglom's law.