code_saturne User's Forum

Does anyone have a suggestion for solving this problem, I am very keen on understanding how CS can be used to compute these budget terms!

Thank you for your response. Main terms of interest that I want to calculate are for example for the velocity field u_i = U_i + u'_i the terms: turbulent dissipation: \epsilon = \nu \left\langle \frac{\partial u'_i}{\partial x_j} \left[ \frac{\partial u'_i}{\partial x_j}+\frac{\partial u'_j}{\partia...

I am attempting to compute turbulent kinetic energy budget terms in an LES in CS V5. For example, in order to calculate the dissipation term I require to compute the average of the product of the velocity gradient tensor with itself. Calculation of time averages is currently done in cs_user_paramete...

Hi Yvan, I found the problem and it was (as expected) a simple mistake on my part. I was not setting the intended velocity at the inlet when specifying the entry velocity in the 'advanced option'. On this point I have a follow up question regarding this issue: what is the difference between the foll...

Ok so the problem does not seem to be with the type of turbulence method because when I run without the advanced subroutine in the cs_user_les_inflow.f90 I now have correct looking mass fluxes at the inlet and outlet! Could it be that I am simply inadvertently assigning zero inlet velocity when usin...

Thanks for the quick reply. I am not sure how to activate the debug build as I am running on a cluster and I don't think I can pick and choose at the moment. However I will try the other Les inflow options and see if that solves the NaN problem. If you would like I can also send you the actual mesh ...

Hello, You can see the listing file inside the Same directory of SRC. I can't post my mesh because it's too big but in general it is a good quality hexahedron mesh. The CFL number is as you can see in the listing file. I have also tried with much smaller time steps but the result is the same. I find...

Can someone tell me why I get the strange output in listing (zero residuals and the message "Incoming flow detained on faces...") This happened when I started attempting to use the les inflow module with the SEM enabled. Please refer to my SRC and listing (both contained inside the zip file)

Brian - are you saying that in the case of two blade rows (say a rotor and stator row) it is possible to have a single moving rotor mesh passage which slides past a single stator mesh passage (provided there is an integer ratio between the number of blades) in code_saturne!? If so I would be very in...

Apologies - I think I see now that the rotation vector specified in rotation_define has elements that are the omega in each direction! However I am not that the above is the correct way to set up the simulation in a rotating reference frame. Could someone confirm?