00001 { 00002 fvScalarMatrix hEqn 00003 ( 00004 fvm::ddt(rho, h) 00005 + fvm::div(phi, h) 00006 - fvm::laplacian(turbulence->alphaEff(), h) 00007 == 00008 DpDt 00009 ); 00010 00011 hEqn.relax(); 00012 hEqn.solve(mesh.solver(h.select(finalIter))); 00013 00014 thermo.correct(); 00015 }