TEqn.H
Go to the documentation of this file.
00001 { 00002 volScalarField kByCv 00003 ( 00004 "kByCv", 00005 (alpha1*k1/Cv1 + alpha2*k2/Cv2) 00006 + (alpha1*rho1 + alpha2*rho2)*turbulence->nut() 00007 ); 00008 00009 solve 00010 ( 00011 fvm::ddt(rho, T) 00012 + fvm::div(rhoPhi, T) 00013 - fvm::laplacian(kByCv, T) 00014 + p*fvc::div(phi)*(alpha1/Cv1 + alpha2/Cv2) 00015 ); 00016 00017 // Update compressibilities 00018 psi1 = eos1->psi(p, T); 00019 psi2 = eos2->psi(p, T); 00020 }