applications/solvers/compressible/sonicFoam/eEqn.H

   1 {
   2     solve
   3     (
   4         fvm::ddt(rho, e)
   5       + fvm::div(phi, e)
   6       - fvm::laplacian(turbulence->alphaEff(), e)
   7       ==
   8       - p*fvc::div(phi/fvc::interpolate(rho))
   9     );
  10
  11     thermo.correct();
  12 }