00001 { 00002 solve 00003 ( 00004 fvm::ddt(rho, e) 00005 + fvm::div(phi, e) 00006 - fvm::laplacian(turbulence->alphaEff(), e) 00007 == 00008 - p*fvc::div(phi/fvc::interpolate(rho)) 00009 ); 00010 00011 thermo.correct(); 00012 }