code_saturne User's Forum

Hi,
I am Running 2D convergence nozzle with Enthalpy equation in version 7.1.2. the area ratio is 4 so the pressure and velocity change along x direction is significant.
The Enthalpy equation described in theory guide page 17 is as follows: here dp/dt is the lagrangian derivative of pressure.

Since pressure and velocities are changing I should expect change in Enthalpy (should be decreasing ) but I see not change at all.
I have attached the setup and subroutines I am using.

Subroutines used

Hello,

I am not sure what is really coded here, since we do not compute Lagrangian derivatives (to my knowledge). If this is a simple unsteady term, that would be more consistent with your results. I'll be happy if someone who is more familiar with this part of the theory guides would answer in my place here, but I have the impression that something might be lost when moving from the Lagrangian theory to Eulerian operators (i.e. things might be simplified and some terms dropped).

Checking the actual code is of course the safest solution, even if it takes a bit of extra time. That would be in covofi.f90 for versions 8.0 and older, and cs_solve_equation.c in the master (dev) branch.

Best regards,

Yvan