00001 { 00002 fvScalarMatrix hsEqn 00003 ( 00004 fvm::ddt(rho, hs) 00005 + mvConvection->fvmDiv(phi, hs) 00006 - fvm::laplacian(turbulence->alphaEff(), hs) 00007 // - fvm::laplacian(turbulence->muEff(), hs) // unit lewis no. 00008 == 00009 DpDt 00010 + chemistrySh 00011 ); 00012 00013 hsEqn.relax(); 00014 hsEqn.solve(); 00015 00016 thermo.correct(); 00017 00018 Info<< "T gas min/max = " << min(T).value() << ", " 00019 << max(T).value() << endl; 00020 }