code_saturne User's Forum

Hello, I am a little bit confused about the correct way of setting Code_Saturne to use a Boussinesq approach or a low-Mach approach. I use v6.0.0. In the GUI, one has the possibility of selecting Incompressible or Compressible flow. I guess that in both cases the Incompressible option has tu be chos...

Hello, In code_Saturne v5, what I understand is that the variable time step option cannot be used with LES, since second order time scheme in automatically activated and does not support it. My question: is there a way to use a CFL driven time step with LES simulations? For example can I somehow swi...

Dear Luciano,

Thanks to your help, I managed to calculate the forces on the faces I wanted to.

Thank you very much!

Best regards,

Daniele

Hello, I agree with you about the example. Nevertheless, in code_saturne doxygen "iforbr" is described to be "field id of the stresses at boundary". That's why I am confused. It looks as if iforbr was just mu*du_i/dx_i. Moreover, are pressure forces taken into account inside ifor...

Dear Luciano, Thanks a lot, you are right, the correct field is iforbr. I am now able to calculate the Stress on the cell surfaces. I still have one problem. In order to obtain the force, I have to multiply the Stress by the boundary surface of the cell. The surface value is given by surfbn(ifac) if...

Dear Luciano,

Thank you for your suggestion. Adding the option does not prevent the problem to occur
I attached also the listing showing some details of the error.

Best regards,

Daniele

Dear Yvan,

Thank you for your help. Unfortunately, even changing the format I get the same error.
I attach the cs_user_extra_operations.f90 I use, just in case it could help more.

Thanks a lot.

Best regards,

Daniele

Hello, I am still trying to get the forces on selected faces, but I have some issues. I simply write the following, inside the cs_user_extra_operations.f90: call getfbr('STRUCT_WALL and y<0.5', nlelt, lstelt) do ilelt = 1, nlelt ifac = lstelt(ilelt) fx_1=fx_1+forstr(1,ifac) enddo This seems to work ...

Hello Yvan, Attached you will find the xml and the ALE routines. I am simulating a vibrating tube in a tube bundle. I set the mesh viscosity constant (with a high value) on subdomains corresponding to the tubes boundary layer, and mesh viscosity as a function of the distance in the mesh around the m...

Hello, I have set up the ALE simulation through the subroutines: it seems to work correctly, I don't have the boundary condition problem anymore. Is there a way to double check the ALE boundary conditions (fixed wall, sliding, ...) defined through the cs_user_boundary_conditions-ale.f90 (I mean, be ...