code_saturne User's Forum

Hello Everyone,

I am currently running Rayleigh Benard Natural Convection simulations. Ra 1e5 and Pr 0.025
I have encountered a strange error in temperature profile. I am noticing these rings along the plumes of upward flow i.e. Mushroom like structures in temperature contours. I have attached a ...

Hello Yvan,

Thank you for your reply. Currently I am using following code which adds heat source in each cell and then removes from the boundary cells. But I am not entirely sure if this the correct way of writing it. Boundary faces have colors 1to5.

do iel=1,ncel
crvexp(iel) = volume(iel)*pwatt ...

Hello Yvan,

I am trying to impose a volumetric heat source thoughout the domain(not on boundary cells- they have different BCs) in an non regular geometry. I am trying to use following
do iel=1,ncel
crvexp(iel) = volume(iel)*pwatt/voltf;
enddo
for imposing the heat source term. Do these 'ncel ...

Hi,

Thanks for the clarification. Another question if you do not mind. If you use Boussinesq through cs_user_source_terms.f90 under ustsnv what did you do with the gravity vector in the generation rate of turbulence (Gk)?

Cheers,
Costas

Hello,
I gave the value of gravity as a negative vector ...

Hello Everyone,

I am running Rayleigh Benard natural convection case with Rayleigh number to be 630000 and Pr as 0.71. I am using AFM as turbulent heat flux model. For boussinesq implementation, I am adding boussiesq(rho*beta*(T-t_ref)*g) source term in momentum equation.
I am facing an unusual ...

Hello,

I used ' cs_user_source_terms.f90' under ustsnv for Boussinesq approximation.

Regards
Kunal

Hello Everyone,

I am running Rayleigh Benard natural convection case with Rayleigh number to be 630000 and Pr as 0.71. I am using AFM as turbulent heat flux model. For boussinesq implementation, I am adding boussiesq(rho*beta*(T-t_ref)*g) source term in momentum equation.
I am facing an unusual ...

Hello Everyone,

I am running a steady state RANS simulation of Rayleigh Benard natural convection where I would like to compare temperature variance across the domain with DNS data. I wanted to know is there any way to extract RMS temperature or temperature variance (\vec{t^2}) as extracting the ...

Hello Yvan,

Thank you so much.
Following worked fine. Regards
Kunal

Hello Yvan,
I used following piece of code as suggested by you and encountered this error : "Rank mismatch in argument 'p' at (1) (rank-1 and rank-2)" . Compile log attached below. Does this mean there's a dimensional mismatch with Rij and defined pointer ? Could you please suggest any alternate ...