[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

[dduque]

No answers? Well, I'll try to implement it myself, while listening to this tune: Gravity's Gone.

[dduque]

OK, I have just modified Teqn.H to solve this equation, including the
gravity power term, and it seems to be working.

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

Best.

[markusrehm]

Hi Daniel,

consider posting a bug report to the Mantis:

http://www.openfoam.org/mantisbt

I am not sure if the developers follow the forum.

Markus

[dduque]

Actually, I was wrong: there is NO gravity term in the energy equation. At least, in this energy equation. There is one in the conservation version of the energy equation (for E=rho e + rho u^2/2), but dissapears in this version. Sorry!

[alundilong]

Originally Posted by dduque

Actually, I was wrong: there is NO gravity term in the energy equation. At least, in this energy equation. There is one in the conservation version of the energy equation (for E=rho e + rho u^2/2), but dissapears in this version. Sorry!

Yes. However, I did not see any gravity term, rho g&U, in the "total energy form".