00001 tmp<fv::convectionScheme<scalar> > mvConvection 00002 ( 00003 fv::convectionScheme<scalar>::New 00004 ( 00005 mesh, 00006 fields, 00007 phi, 00008 mesh.divScheme("div(phi,ft_b_h_hu)") 00009 ) 00010 ); 00011 00012 if (composition.contains("ft")) 00013 { 00014 volScalarField& ft = composition.Y("ft"); 00015 00016 solve 00017 ( 00018 fvm::ddt(rho, ft) 00019 + mvConvection->fvmDiv(phi, ft) 00020 - fvm::laplacian(turbulence->alphaEff(), ft) 00021 ); 00022 } 00023 00024 // ************************ vim: set sw=4 sts=4 et: ************************ //