code_saturne User's Forum

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

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