code_saturne User's Forum

Thank you for your response. In fact, when the atmospheric model (dry atmosphere) is activated, with the wall_temperature chosen in GUI, it does not output the real boundary condition imposed at the rough walls. It gives some other strange large values for the temperature. By the way, for the other ...

Version 3.0.6 I am using the thinwall to add internal boundary condition. In cs_user_mesh.c, I select the name of the internal faces, for example, 'wall0' In cs_user_boundary_condition.f90, I imposed the boundary condition on the thinwall as call getfbr('wall0', nlelt, lstelt) do ilelt = 1, nlelt if...

Version 3.0.6 We would like to check the boundary condition, for example, the temperature at the wall, is set correctly, by visualizing it in Ensight-10. Apparently, there was a change since version 3.0 that only the cells attached to the wall can be visualized by default. When I tried to visualize ...

Hello, I tried to set rtpa, it gives the following errors during the compilation Boussinesq4/SRC/cs_user_initialization.f90(130): error #6917: This array with adjustable declarators must be declared a dummy argument. [RTPA] double precision dt(ncelet), rtp(ncelet,*), propce(ncelet,*), rtpa(ncelet,*)...

Dear Yvan, I have a further question about the temporal second order scheme. For the Rayleigh-Benard convection, we have the base flow, an analytical solution: u = 0; v = 0; w =0; linear profile for temperature. I use the cs_user_intialisation.f90 to initialize this field. The table rtp can be set, ...

Thank you for your reply. Here I give some feedbacks and pose some new questions. I indeed observed the Rayleigh-Benard cells by runing code_saturne. In fact, we need the width D between the plates is set to about 10mm for air, if the \beta is set to 1/T_0 as suggested in the tutorial. However, if w...

I am not answering your question, but asking you a question. How do you run adimensional simulation for Navier-Stokes equation? The properties of fluid always has some unities when we input it by GUI or user subroutines. It is need to do some mathematical manipulation to find the right Renoylds numb...

Thank you for your reply. It is time to give some feedback. I did some search in the code version 3.0.6. It seems that we can declare a new "propce", but the number of variables is not automatically incremented for the global number of variables, unless they are already declared in the spe...

I am using the version 3.0.6 I added a user defined parameter served to output a scalar, which is based on the flow properties: velocity amplitude, temperature and radiative intensity. Since it is unnecessary to be solved by the transport equation, I think it is better to declare this variable as a ...

Dear Yvan, Thank you for your reply. In SRC/REFERENCE/cs_user_parameters.f90, I found the subroutine usipsu, but there is no option for ischtp. Could you tell me in which file it is? I am using the version 3.0.5 Best regards Zhenlan Hello, In any case, for DNS-type computations, you will need a CFL ...