code_saturne User's Forum

Than you for the precisions. Do you plot the norm of the shear stress or the x component (x being in the fluid direction)? Pressure being defined up to a constant for an incompressible flow, we can shift it from a constant (but the integral of a constant over a closed solid is 0). Another question: ...

hello,

Which version do you use ?

Best Regards,
Martin

Hello, variables are solved in increments, so if a source term write S^{exp} + S^{imp} var then it gives with \delta var = var^{n+1} - var^n : \rho Vol \delta var / \Delta t - S^{imp} \delta var + ... = S^{exp} + S^{imp} var^n Note implicite par is put on the left hand side if - S^{imp} > 0 . Martin

Thank you for your answer. I have no time to check the resu tonight but I can try to answer some of your questions : - z ground is used to imposed initial and boundary profile for the velocity and the (potential) temperature with relief. In Monin Obukhov profiles we take the distance to the ground i...

Dear SimBu,
I suspect there is something wrong in you user subroutines because a ground has already been computed (from 0 to 50m) and should not be computed again.

Do you have any user subroutines which would explain the problem?
Best regards
Martin

Hello, Just to answer to your first question: in the incompressible module the friction source term is neglected and not in the code. You can add it as a user source terme. For the compressible module, the solved equation is the total energy and the terme is handelled. Concerning the wall function, ...

Hello,

I think you should impose a flux rather a Dirichlet BC because otherwise the ongoing value of scalar is only from diffusion flux.

Could test that?
There is also a possibility to perform scalar balance in the GUI.

Best Regards
Martin

Hello, Indeed eq I.5.25 has a typo (no Clog). I will correct that. Concerning the two velocity scale, yes there are other references you can read (Craft, and also PHD of Jean-Francois Wald). The derivation of the BCs of k and epsilon is detailed in appendix (which help reading the code). Indeed the ...

Dear user, Thank you for the feedback. We are aware of this problem which is due to the fact, if I understand well, that the additional force you have on the mean fluid (with the head loss) is missing in the lagrangian module (all terms on the momentum are not automatically taken as they should). We...

Dear Daniele, The CDO solver used for ALE is not compatible yet with FSI (and INTERNAL coupling) because the way boundary conditions are treated are really different (using the GUI it is transparent, but you cannot use cs_user_boundary_conditions_ale.f90 for the CDO solver). We check if it is possib...