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