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.
Enthalpy equation in code_saturne
Forum rules
Please read the forum usage recommendations before posting.
Please read the forum usage recommendations before posting.
Enthalpy equation in code_saturne
 Attachments

 setup.xml
 (11.17 KiB) Downloaded 150 times
Re: Enthalpy equation in code sauterne
Subroutines used
 Attachments

 cs_user_initializationbase.c
 (3.68 KiB) Downloaded 151 times

 cs_user_physical_properties.f90
 (5.76 KiB) Downloaded 145 times

 cs_user_parameters.c
 (3.29 KiB) Downloaded 150 times

 Posts: 4121
 Joined: Mon Feb 20, 2012 3:25 pm
Re: Enthalpy equation in code sauterne
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
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