code_saturne User's Forum

Hello, since version 6.0 there is an option to activate the Boussinesq approximation (indeed rho varies only in the buoyancy term and is supposed to be constant everywhere else).
Look for idilat = 0 in cs_user_parameter file.

Best Regards
Martin

Hello, If you want to model a free surface, you can also test the ALE module and then “free surface” boundary condition will be used. It will impose p0 at the freesurface and let it move with the flow. So if you have breaking waves this method does not work. We have also experimental options for tha...

Hello, It should be compatible but we don’t have any validation test case for that so that would be nice if you could set up a toy model test case to check if it is working correctly. If not, that will not be a big deal to fix (ALE is mainly impacting the other fields updating convective mass flux, ...

Hello, Concerning the initialisation you can read paragraph 6.3 of the user guide of version 5.0 for instance. Concerning the inlet BC for turbulence, you can either chose established profile with an hydraulic diameter given by the user, or hydraulic diameter and a turbulence intensity. See for inst...

Hello, Is it a valve you want to model? If yes, this boundary should be considered as moving walls. If it is inlet outlet boundary conditions, as mention in the user and the theory guide, isolib BC type retain the flow (imposing a 0 mass flux) when entering (it is a free outlet, but not inlet is all...

Hello Daniele,
That is not normal. What is the version you are using?
Best regards
Martin

Dear Odion,

a call function was missing for the vector coupling. I added that to the trunk version. The computation is not crashing anymore but the solvers are not converging properly.

I look at it in the coming days.

Best Regards

Martin

Hello,
You can activate a variable diffusivity and give the same value for all the cells. We have to correct this combination.

Best regards

Martin

Hello again,
Just a small additional remark: is "ni" a component of a vector? If so, I would suggest you to move to the trunk version because we add the capability to transport vector species (or if not it will be in the coming days).
Best regards
Martin

Hello Paul, Reading your problem, I think you can do it using DRIFT modelling. Your scalar is nothing else that transported by the velocity u plus a drift (proportional to grad phi). This is possible in version 4.0 to activate drift to a scalar (look at DRIFT in cs_user_parameters.f90) but usually t...