Code_Saturne User's Forum

Hello, As Benoît mention the problem was in the pressure. In your case if you set a P_{ref}=0 and the outlet pressure according to the cavitation number \sigma you will get the cavitation of the fluid at the throat. I have run you case with a dt=5e-6 [s] and the results are reasonable according to t...

Hello, Inside of the boundary faces loop, you can access to the coordinates of the centers of the boundary faces thought the pointer cdgfbo ! Boundary faces of group '3' call getfbr('3', nlelt, lstelt) do ilelt = 1, nlelt ifac = lstelt(ilelt) ! y-coordinate Face center ycoord = cdgfbo(2, ifac) enddo...

Hello Amir, I was looking the cavitation.f90 file and the theory documentation, because I never work with cavitation. In the theory documentation of the CS 4.0.0 the vaporization source terms are m^+ = -\dfrac{C_{prod} \rho_l \min \left( p-p_V,0 \right)}{0.5\rho_lu_\infty^2t_\infty}, m^- = -\dfrac{C...

Hello amir,

Could you also post the case and the mesh files?

Regards,

Luciano

Hello Amir, In the listing file I have noted that the Uinf for the Merkel model is set to Uinf=-0.1e+14 , wherewith the evaporation and condensation rates will be almost zero. You can check this in the therory documentation (http://code-saturne.org/cms/sites/default/files/docs/4.0/theory.pdf) and in...

Hello,

I think that the subroutine that you are looking for is in /REFERENCE/cs_user_particle_tracking.f90

Also you can check /EXAMPLES/cs_user_lagrangian_particle_tracking_module.f90

Regards,

Luciano

Ok, Did you check cs_user_boundary_conditions-base.f90 , cs_user_boundary_conditions-advanced.f90 , cs_user_boundary_conditions-mapped_inlet.f90 ? The file with the velocity comes from another simulation with a different mesh or same mesh? If the profile comes from a different mesh you will have to ...

Hello,

Do you want to use this velocity field to initialize the simulation?. If so, please check the cs_user_initialization-base.f90 and cs_user_initialization-compressible.f90 files at /SRC/EXAMPLES

Regards,

Luciano

Hello,

I think that you have a typo at a boundary condition. I have checked your listing and find this
May be "joint5" is "join5"

Regards,

Luciano

Hello,

In my test case I have set inlet and outlet.

I attach the case file, the mesh and the fan function. If this case run without problem post your case in order to help you to debug the problem.

Regards,

Luciano