Code_Saturne User's Forum

Hello, Maybe the error is that the file does not exist and can not be appended the data. Try adding open(impout,file='force_1_x.dat', status="unknown" , position="append") When I write some data to a file, in the first iteration I create the file, write some header with info abou...

Hello,

Yes, set a prefixed (Dirichlet) value for the scalars at the fluid inlets (in the boundary condition tab). Additionally you can set different scalar for the different inlets.

Regards,

Luciano

Hello,

This error is strange because in the subroutine uitssc the lambda coefficient is used for groundwater flows.

Another way to track the scalars is set an inlet value for each scalar and the track the evolution. Using this way you avoid the use of the source term.

Regards,

Luciano

Hello, cs_glob_rank_id will give you the rank of a MPI process in case of parallelism and it will take the values between 0<=cs_glob_rank_id< number of process. In case of a serial run cs_glob_rank_id =-1 , so only process of rank=0 (parallel) or -1 (serial) will write to this file. In you code, the...

Hello, But these "points" should be associated to some surface or area in order to compute forces, or chose the boundary face in the rigid body which contain or is placed the point. After you have selected the face or surface(faces) you can computed the forces using the cs_user_extra_opera...

Hello, I think that the answer is yes (I'm not totally sure because I never have done), but if you have created a new variable (scalar, vector, etc) and you store some data, then you can access from both, fortran or C. Maybe someone from the development team could answer more precisely to your quest...

Hello, You can call the fortran subrutine in a C function doing some BIND(C) to be compatible with C. May be this link could help you. https://gcc.gnu.org/onlinedocs/gfortran/Interoperable-Subroutines-and-Functions.html Maybe you can check in the post.f90 file the subrutine cs_f_post_write_var , whi...

Hello, I really don't understand very well the question....but the forces computed are normal and tangencial to specific faces and also you can get the X,Y,Z with respect to the global reference. If the object is moving in the Z and X direction I don't understand why do you get forces in the Y direc...

Hello Yvan, Thanks for the help, maybe my case setup is wrong. Additionally, I get the same behavior when I use "anisotropic_turbulent_viscosity" in order to set anisotropic conductivities. I have added a quite similar block of code for the contribution of the boundary faces to the cs_anis...

Hello Yvan, Hello Luciano, Is the top left picture using the modified or unmodified version ? The top left picture is the unmodified version from the master branch. On the right, I do not see the orange curve. Is it superimposed with another Yes, the orange is superimposed with the another curves. T...