Scalar flux consistent models

In transported PDF methods (Pope 2000), the closure model for A, V, ip) will be a known function26 ofV. Thus, (U,Aj) will be closed and will depend on the moments of U and their spatial derivatives.27 Moreover, Reynolds-stress models derived from the PDF transport equation are guaranteed to be realizable (Pope 1994b), and the corresponding consistent scalar flux model can easily be found. We shall return to this subject after looking at typical conditional acceleration and conditional diffusion models. 

We shall see that a conditional acceleration model in the form of (6.48) is equivalent to a stochastic Lagrangian model for the velocity fluctuations whose characteristic correlation time is proportional to e/k. As discussed below, this implies that the scalar flux (u,

joint velocity, composition PDF level, and thus that a consistent scalar-flux transport equation can be derived from the PDF transport equation. 

Unlike Lagrangian composition codes that use two-equation turbulence models, closure at the level of second-order RANS turbulence models is achieved. In particular, the scalar fluxes are treated in a consistent manner with respect to the turbulence model, and the effect of chemical reactions on the scalar fluxes is treated exactly. 

One expects from lowest-order QCD that the short range part of the potential (which arises mainly from one-gluon exchange) is vector in character. On the other hand, lattice QCD calculations and some flux-tube models indicate that the confining potential is a scalar. The experimental evidence is not conclusive, but is consistent with this picture. If these ideas are correct, then there should not be any appreciable spin-spin splitting in quarkonium states with L > 0 because in such states the wave function is small at short distances because of the centrifugal barrier. Measurements of the spin-spin splitting in L > 0 states therefore provide important tests of QCD-motivated potentials. 

In these equations fi is the coluirm mass of dry air, V is the velocity (u, v, w), and (jf) is a scalar mixing ratio. These equations are discretized in a finite volume formulation, and as a result the model exactly (to machine roundoff) conserves mass and scalar mass. The discrete model transport is also consistent (the discrete scalar conservation equation collapses to the mass conservation equation when = 1) and preserves tracer correlations (c.f. Lin and Rood (1996)). The ARW model uses a spatially 5th order evaluation of the horizontal flux divergence (advection) in the scalar conservation equation and a 3rd order evaluation of the vertical flux divergence coupled with the 3rd order Runge-Kutta time integration scheme. The time integration scheme and the advection scheme is described in Wicker and Skamarock (2002). Skamarock et al. (2005) also modified the advection to allow for positive definite transport. 

---