Search found 16 matches
- 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...
- 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...
- 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 ...
- 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...
- 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...
- 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
I used ' cs_user_source_terms.f90' under ustsnv for Boussinesq approximation.
Regards
Kunal
- 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...
- 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...
- 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.
Regards
Kunal
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)
Kunal
- 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...