code_saturne User's Forum

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

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

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

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...

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...

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 &...

Hello, In the predefined symbols the temperature doesn't appear, but is used in the example tab at least in CS 4.0. In the example that I have posted I only use the GUI and it seems to works. Additionally, you mention that the density is a constant, but you are using an expansion coefficient, so the...