Search found 239 matches

by Luciano Garelli
Fri May 27, 2016 8:37 pm
Forum: code_saturne usage
Topic: Jacobi: error (divergence) solving for Velocity
Replies: 17
Views: 6907

Re: Jacobi: error (divergence) solving for Velocity

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...
by Luciano Garelli
Thu May 26, 2016 3:50 pm
Forum: code_saturne usage
Topic: Velocity profile from file
Replies: 6
Views: 3450

Re: Velocity profile from file

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...
by Luciano Garelli
Tue May 24, 2016 9:02 pm
Forum: code_saturne usage
Topic: Jacobi: error (divergence) solving for Velocity
Replies: 17
Views: 6907

Re: Jacobi: error (divergence) solving for Velocity

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...
by Luciano Garelli
Fri May 20, 2016 2:16 pm
Forum: code_saturne usage
Topic: Jacobi: error (divergence) solving for Velocity
Replies: 17
Views: 6907

Re: Jacobi: error (divergence) solving for Velocity

Hello amir,

Could you also post the case and the mesh files?

Regards,

Luciano
by Luciano Garelli
Tue May 17, 2016 3:28 pm
Forum: code_saturne usage
Topic: Jacobi: error (divergence) solving for Velocity
Replies: 17
Views: 6907

Re: Jacobi: error (divergence) solving for Velocity

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...
by Luciano Garelli
Mon May 16, 2016 4:36 pm
Forum: code_saturne usage
Topic: Particle drag law
Replies: 2
Views: 1400

Re: Particle drag law

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
by Luciano Garelli
Fri May 13, 2016 1:41 pm
Forum: code_saturne usage
Topic: Velocity profile from file
Replies: 6
Views: 3450

Re: Velocity profile from file

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 ...
by Luciano Garelli
Thu May 12, 2016 4:23 pm
Forum: code_saturne usage
Topic: Velocity profile from file
Replies: 6
Views: 3450

Re: Velocity profile from file

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
by Luciano Garelli
Tue May 10, 2016 7:37 pm
Forum: code_saturne usage
Topic: Jacobi: error (divergence) solving for Velocity
Replies: 17
Views: 6907

Re: Jacobi: error (divergence) solving for Velocity

Hello,

I think that you have a typo at a boundary condition. I have checked your listing and find this

Code: Select all

The group "joint5" in the selection criteria:
"joint5 or join6"
 does not correspond to any boundary face.
May be "joint5" is "join5"

Regards,

Luciano
by Luciano Garelli
Mon May 09, 2016 1:47 pm
Forum: code_saturne usage
Topic: Axail fan modelling usin cs_fan.h
Replies: 16
Views: 3572

Re: Axail fan modelling usin cs_fan.h

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