System of mix mass-fraction beta SDEs.
This file implements the time integration of a system of stochastic differential equations (SDEs) with linear drift and quadratic diagonal diffusion, whose invariant is the joint beta distribution. There are two differences compared to the plain beta SDE (see DiffEq/Beta.h):
- First, the parameters, b, and kappa are specified via functions that constrain the beta SDE to be consistent with the turbulent mixing process. In particular, the SDE is made consistent with the no-mix and fully mixed limits. See, e.g., MixMassFractionBetaCoeffConst::update().
- Second, there two additional random variables computed, the same as also computed by the mass-fraction beta equation, see also DiffEq/MassFractionBeta.h.
In a nutshell, the equation integrated governs a set of scalars, , , as
with parameter vectors , , and . This is similar to DiffEq/Beta.h, but the parameters, and constrained. Here . The fluctuation about the mean, , is defined as usual: , and and are user-specified constants. Also, is an isotropic vector-valued Wiener process with independent increments. The invariant distribution is the joint beta distribution. This system of SDEs consists of N independent equations. For more on the beta SDE, see https:/
In addition to integrating the above SDE, there are two additional functions of are computed as
These equations compute the instantaneous mixture density, , and instantaneous specific volume, , for equation in the system. These quantities are used in binary mixing of variable-density turbulence between two fluids with constant densities, and . The additional parameters, and are user input parameters and kept constant during integration. Since we compute the above variables, and , and call them mixture density and specific volume, respectively, , governed by the beta SDE is a mass fraction.
All of this is unpublished, but will be linked in here once published.
- namespace walker
- Walker declarations and definitions.