//CMC mixture fraction variance equation 
//cf) mfVar should not be ZERO

        SDR =  turbulence->epsilon() * mfVar / turbulence->k();
        volVectorField Gradmf = fvc::grad(mf);

	//Info<<"CMC mixture fraction variance equation start"<<endl;
        solve
        (
                fvm::ddt(rho,mfVar) 
		+ mvConvection->fvmDiv(phi, mfVar)
                - fvm::laplacian(1.47*turbulence->mut(), mfVar) //let 1/Sc = 1.47
                //- fvm::laplacian(1.47*turbulence->muEff(), mfVar) //let 1/Sc = 1.47
                ==
                  2 * ( 1.47 * turbulence->mut() ) * (Gradmf & Gradmf)
                //  2 * ( 1.47 * turbulence->muEff() ) * (Gradmf & Gradmf)
                - 2 * rho * SDR
                //+ 2*(0.5)*sqrt(mag(mfVar))*Sevap*(1-mf)
        );
	//Info<<"CMC mixture fraction variance equation end"<<endl;