code_saturne User's Forum

Hello, As a convention, we take negative mass flux for inlets (because as a convention, the normal to the boundary is outwarding) so no matter of the direction of z. Concerning your crash, what I advice is to extrude your mesh next to the outlet (because the pressure BC is supposing that the mesh is...

Hello, You can specify the initial water elevation by modifying this vectorial field in cs_user_initialisation.f90. Then, by default, a imposed pressure is imposed at the free surface (with p0 value) and a homogeneous Neumann on the velocity. Concerning the mesh velocity, the normal component to the...

Hello, If you want to specified an initial mesh_displacement (for an initial water elevation), you can do it in cs_user_initialization. The boundary condition imposed in case of free surface is an imposed pressure (at p0) and an homogeneous Neumann on the fluid velocity (so no additional shear stres...

Hello,
Could you provide some screenshots to understand your problem?

Best regards
Martin

Hello, I cannot debug from my phone, but that it may happen that the cylinder is moving too much, and therefore the algorithm which is computing the mesh displacement (which is not really robust) make nodes crossing one to the other (and you obtain a negative cell volume). What I would advice is: - ...

Hello,
You can export the stress of the fluid on the structure directly in the GUI. The normal component is therefore the pressure (more or less).
By the way, The coupling between Code_Saturne and aster is already working in version 4.0.3.
Best regards
Martin

Hello Vilas, Ok now I understand what you try to do. You can already do that specifying idften(iu)=6 in cs_user_parameters.f90 (let say in usipsu subroutine), and by computing your mu_T in phyvar (see in src/base). Then mu_T is a 6*ncelet array (the components are [mu_T(11), mu_T(22), mu_T(33), mu_T...

Hello, Do you want to add "div(R)" explicitly (if yes the routine divrij already do that) or do you have a constitutive law stating "R =-2 nu Grad u +..." With nu a tensor (if yes you can use directly bilscv subroutine and specify the option on the variable "idften(ivar)=6&q...

Hello, Which version do you use? Cfdivs is used only for the compressible module to add explicitly div(sigma.u) in the total energy equation. Which equation would you like to solve? If you want to solve momentum, look at predvv.f90, and the convection diffusion step is done in cs_convection_diffusio...

Hello, At least for the term "mu_T grad u" you can have an idea of how to do so by greping iporos.eq.2 in the code. Unfortunately the brick to code "mu_T transpose(grad u)" does not exist yet for a tensor mu_T, even if it is not that difficult to adapt it from the scalar brick. I...