code_saturne User's Forum

Dear Benoît, Could you please help me to confirm another question? When I did the test case, I made a small change to your source term subroutine. I changed it from the Boussinesq approximation to ideal gas law. I defined the density as P/(R*T) in another subroutine, and call the density in the sour...

Dear Benoît, Thank you very much for your codes. It's very helpful. Using your source term subroutine, I did a test case that is described in the paper in my previous post. I compared Code_Saturne result with other solvers' results, which is good. The result, mesh and setup file are attached. Hope i...

Dear developers, Please could you help me with this restarting problem? I am running version 7.2 alpha, compiled from GitHub source code, on a cluster. When I restart, I always found an error, "Invalid (not-a-number) values detected for a field." The listing, xml and subroutine I used are ...

Dear Yvan, Thanks a lot for your help. I coded the ideal gas law into cs_user_physical_properties.f90 and got good results. The code is attached below. Many thanks, Ruonan call field_get_val_s(icrom, cpro_crom) call field_get_val_s(iprtot, cpro_prtot) call field_get_val_s(ivarfl(isca(iscalt)), cvar_...

Hello Benoît and Yvan, Thank you very much! I really appreciate the help from you both. I now can understand the icorio subroutine, I can see the rotating frame is working. Please let me explain my object and problem further: What I want to simulate is Rayleigh-Benard convection, but it's not in gra...

Hello Yvan, I think I may found another bug in the turbomachinery model. Could you please help to confirm? Following your last suggestion, I tried "setting as a transient turbomachinery case, with only a rotor". I found in RANS it's ok, but in LES, it doesn't work. The rotating boundary co...

Dear developers, Could you please help me with this question? Following the heated square cavity tutorial, I used Boussinesq approximation, make density change with temperature. Now I want to make density change with both pressure and temperature, that is the ideal gas law, rho=P/(R*T). I think this...

Hello Yvan, I really appreciate your time and help! Please find some idea based on my tests: 1. When I define rotation vector either in GUI or subroutine, the case can run with no error, thanks to your debugging. When I post-process the result, I can find two additional variables, "Abs_Pressure...

Hello Yvan,

May I confirm if the plus or minus of omega component in the "icorio" follows the right-hand law? For example, if I look along the x axis and the model is rotating anticlockwise, should I set omega_x=-Ω rad/s?

Best regards,
Ruonan

Hello Yvan, I can't thank you enough. Your bug fixing helped me a lot. Now the "icorio" subroutine can work with no error. Another question please: Now I want the annuli to rotate along x axis with 8.159 rad/s, and I want the whole domain to be solved in a rotating reference frame. I set o...