00001 { 00002 fvScalarMatrix hEqn 00003 ( 00004 fvm::ddt(rho, h) 00005 + mvConvection->fvmDiv(phi, h) 00006 - fvm::laplacian(turbulence->alphaEff(), h) 00007 == 00008 DpDt 00009 ); 00010 00011 hEqn.relax(); 00012 hEqn.solve(); 00013 00014 thermo.correct(); 00015 } 00016 00017 // ************************ vim: set sw=4 sts=4 et: ************************ //