Code_Saturne User's Forum

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 functions cs_math_3_cross_product()...

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

Hello, As Yvan says, you can write your expression (Boussinesq-type term) through the MEI as #air expansion coefficient (1/K) b=0.00343; #Ref temperature tref=300; # Density of air density = rho0*(1-b*(temperature-tref)); Then in "Gravity" you have to define the acceleration direction. I attach an e...

Hello, The case is small, so I uploaded the mesh, the xml and the subroutines. In order to find where is the problem (I think in my subroutines) I have solved the same problem but now the inner sphere is held fixed and the sourrounding fluid is initialized with a rotating velocity throught the GUI m...

Hello Yvan and thanks for answering. The first step was check the if the velocity field it is what I expected and then the efforts. The velocity was ok, but the efforts looks strange, so I decide to simplify the problem and solve a rotating sphere with the fluid at rest. For this case the inner sphe...

Hi everybody, I'm solving a problem of a rotating object (cylinder) in a fluid, initially at rest, and I want to compute the forces and moments acting on it. The object is placed at the origen and a tangencial (sliding) velocity is prescribed at the surface through cs_user_boundary_conditions.f90. T...