Search found 16 matches

by kunalck
Fri Oct 25, 2019 3:05 pm
Forum: code_saturne usage
Topic: Low value of Turbulent heat fluxes in Natural convection
Replies: 9
Views: 5075

Re: Low value of Turbulent heat fluxes in Natural convection

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 plo...
by kunalck
Tue May 07, 2019 9:50 am
Forum: code_saturne usage
Topic: Several uniform heat source
Replies: 10
Views: 5298

Re: Several uniform heat source

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/vol...
by kunalck
Mon May 06, 2019 8:07 am
Forum: code_saturne usage
Topic: Several uniform heat source
Replies: 10
Views: 5298

Re: Several uniform heat source

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' also ...
by kunalck
Mon Apr 29, 2019 2:26 pm
Forum: code_saturne usage
Topic: Low value of Turbulent heat fluxes in Natural convection
Replies: 9
Views: 5075

Re: Low value of Turbulent heat fluxes in Natural convection

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 in GUI...
by kunalck
Mon Apr 29, 2019 2:21 pm
Forum: code_saturne usage
Topic: Low value of Turbulent heat fluxes in Natural convection
Replies: 9
Views: 5075

Re: Low value of Turbulent heat fluxes in Natural convection

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 pro...
by kunalck
Tue Apr 09, 2019 9:00 am
Forum: code_saturne usage
Topic: Low value of Turbulent heat fluxes in Natural convection
Replies: 9
Views: 5075

Re: Low value of Turbulent heat fluxes in Natural convection

Hello,

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

Regards
Kunal
by kunalck
Mon Apr 08, 2019 3:23 pm
Forum: code_saturne usage
Topic: Low value of Turbulent heat fluxes in Natural convection
Replies: 9
Views: 5075

Low value of Turbulent heat fluxes in Natural convection

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 prob...
by kunalck
Mon Feb 18, 2019 10:31 am
Forum: code_saturne usage
Topic: RMS temperature / Variance calculation
Replies: 1
Views: 1235

RMS temperature / Variance calculation

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 RES...
by kunalck
Thu Jan 31, 2019 1:48 pm
Forum: code_saturne usage
Topic: Accessing Reynolds Stress Components
Replies: 4
Views: 2022

Re: Accessing Reynolds Stress Components

Hello Yvan,

Thank you so much.
Following worked fine.

Code: Select all

double precision, dimension(:,:), pointer :: cvara_rij
call field_get_val_prev_v(ivarfl(irij), cvara_rij)
Regards
Kunal
by kunalck
Wed Jan 30, 2019 1:57 pm
Forum: code_saturne usage
Topic: Accessing Reynolds Stress Components
Replies: 4
Views: 2022

Re: Accessing Reynolds Stress Components

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 alt...