[dduque]

Hi everybody.

We were greatly pleased by the modifications in compressibleInterFoam introduced in the new OpenFOAM v2.2.0 release. We were precisely working in extending the EOS to include non-isoT simulations, and the task is now basically done for us.

However, I notice the equation for temperatura in the new Teqn.H reads:

       fvm::ddt(rho, T)
      + fvm::div(rhoPhi, T)
      - fvm::laplacian(kByCv, T)
      + p*fvc::div(phi)*(alpha1/Cv1 + alpha2/Cv2)

Don't you have a term

rho g·u / c_v

on the right hand side of this eq?

This is the power production due to gravity (well, rho g·u is the power, divided by c_v gives a temperature variation).

Thanks. Best,

Daniel