April 23, 2012
Abstract
         Equation of mixture fraction and variance of mixture fraction in the
     flamelet solver
   // Solving equation of variance of mixture fraction
   fvScalarMatrix csiEqn
   (
        (
             fvm::div(phi, csi)
             - fvm::Sp(fvc::div(phi), csi)
             - fvm::laplacian(turbulence->mut()/sigmat, csi)
        )
   )
                                          e − ∇ • µt ∇Z
                     ∇ • φZe − Sp((∇ • φ)Z)              e                 (1)
                                                   σ
  // Solving equation of variance of mixture fraction
  fvScalarMatrix csiv2Eqn,
  (
         (
            fvm::div(phi, csiv2)
            - fvm::Sp(fvc::div(phi), csiv2)
            - fvm::laplacian(turbulence->mut()/sigmat, csi)
            - Cg*turbulence->mut()*magSqr(fvc::grad(csi))
            + Cd*rho*turbulence->epsilon()/(turbulence->k()+k small)*csiv2
         )
  )
                   f00 )−∇• µt ∇Z                                  
                                              p
∇•φZf00 −Sp((∇•φ)Z               f00 +Cg∗µt ∗| ∇Z  f00 |+Cd∗ρ∗             f00
                                                                          ∗Z
                             σ                                 k + ksmall
                                                                           (2)