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