Search found 245 matches

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

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: 1539

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: 3739

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: 3739

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: 7549

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 The group "joint5" in the selection criteria: "joint5 or join6" does not correspond to any boundary face. May be "join t 5" 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: 4219

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
by Luciano Garelli
Mon May 09, 2016 2:57 am
Forum: code_saturne usage
Topic: Axail fan modelling usin cs_fan.h
Replies: 16
Views: 4219

Re: Axail fan modelling usin cs_fan.h

Hello,

I suspected that the axial torque has the influence that you mention.

Wth respect to the pressure curve, is a function of volumetric flow rate or velocity, because I dont understand what means with flow in th doxygen doc.

Regards

Luciano
by Luciano Garelli
Fri May 06, 2016 7:17 pm
Forum: code_saturne usage
Topic: Axail fan modelling usin cs_fan.h
Replies: 16
Views: 4219

Re: Axail fan modelling usin cs_fan.h

Hello, Some time ago I was working with other CFD software in a simulation with fan, so it is very interesting for me too. I was following the comments of Yvan, I added the "cs_fan.c" to my 4.2.1 version and also have to add "cs_math.c" and "cs_math.h" because the funct...
by Luciano Garelli
Mon Apr 25, 2016 3:11 pm
Forum: code_saturne usage
Topic: Calculation of forces and moments
Replies: 7
Views: 2218

Re: Calculation of forces and moments

Thanks Yvan, Now the results looks good and the values are as expected. One additional comment/suggestion for future versions, when choosing a slinding wall in the Boundary Conditions options through the GUI, only can be applied a traslational velocity. A rotational velocity will be very helpfull an...
by Luciano Garelli
Thu Apr 21, 2016 10:01 pm
Forum: code_saturne usage
Topic: Adding a source term in the momentum equation
Replies: 10
Views: 3107

Re: Adding a source term in the momentum equation

In the subroutine that you post, you use crvexp(3,ncelet) ....... do iel = 1, ncel crvexp(iel) = volume(iel)*beta*rho*(Tref-Temp)*gi enddo And in the documetation crvexp has two argument crvexp(i,iel) for vectorial fields, http://code-saturne.org/doxygen/src/cs__user__source__terms_8f90.html where &...