code_saturne User's Forum

Hello Yvan, Thanks for your quick reply. I have tried your suggestion but I get the same behavior with some iterations in the velocity field in the solid part. In the setup.log appears the text "Fluid-solid mode (disable dynamics in the solid part)" Regards, Luciano setup.log run_solver.log

Hello, I'm running in CS 6.0.4 a simple example using the internal coupling feature, but when setting different densities (rho_air=1, rho_solid=100) spurious velocities appears at the solid region as seen in the picture. The densities and temperature are constant and the gravity is turned on. If gra...

Hello,

I did a test in a small case using CS 6.0.4 and the temperature field in celsius. I used LES dyn model restarting from K-omega SST and I don't get this behaviour in the temperature iterations. Maybe it was a bug that was fixed in the subsequent versions.


Regards,
Luciano

Hello Ok, but setting different blencv in your domain you are mixing first and second order schemes. May be the ischcv = 3 ( Blending between Second Order Linear Upwind and Centered scheme) is what you need, because your will be using two second order schemes and the unphysical oscillations can be r...

Hello,

Did you try with another convective scheme like SOLU (ischcv). Also, you can set different limiters (isstpc)

Regards,
Luciano

Hello, I think that before call cs_parall_allgather_r() you have to do a parallel sum (parcpt) in order to sum "nlelt1" to get the total "n_g_elts". After you get "n_g_elts" you can call to cs_parall_allgather_r() Let me know if I'm right, because I never use it in CS, ...

Hello, After partitioning the mesh each process will have a portion or not of the boundary WALL1, so each process has to create a local array and then you need to join all these array with a mpi gather. You can take a look to the file cs_user_extra_operations-parallel_operations.f90 in the example d...

Hello,

You can check the Doxygen documentation

https://www.code-saturne.org/cms/sites/ ... an_8c.html

and also this old post

viewtopic.php?f=2&t=1831&start=10

Regards,
Luciano

Hello, ibeta it is a (integer) pointer for thermal expansion coefficient. To access to the thermal expansion coeff you can use this example used inphyvar.f90, I think. double precision, dimension(:), pointer :: cpro_beta .... .... if (ibeta.ge.0) then iok1 = 0 call field_get_val_s(ibeta, cpro_beta) ...

Hello, In the file /turb/divrit.f90 are the same lines that you mention, ! AFM model ! "-C_theta*k/eps*( xi* uT'.Grad u + eta*beta*g_i*T'^2)" if (iturt(iscal).eq.20) then if (itt.gt.0.and.ibeta.ge.0) then temp(ii) = temp(ii) - ctheta(iscal)*xtt* & etaafm*cpro_beta(iel)*grav(ii)*cvara_t...