Low value of Turbulent heat fluxes in Natural convection
Forum rules
Please read the forum usage recommendations before posting.
Please read the forum usage recommendations before posting.
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 problem regarding the turbulent heat flux magnitude. The turbulent heat fluxes for this case are incredibly low value(1e-38 to 1e-36). Whereas the temperature profiles are predicted well when compared to DNS results.
Also the turbulent kinetic energy (i.e. rij values) is also incredibly low.
Code_saturne v 5.1.3.
Turbulence model : RSM-EB
turbulent heat flux model : AFM
Kindly suggest few changes.
Regards
Kunal
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 problem regarding the turbulent heat flux magnitude. The turbulent heat fluxes for this case are incredibly low value(1e-38 to 1e-36). Whereas the temperature profiles are predicted well when compared to DNS results.
Also the turbulent kinetic energy (i.e. rij values) is also incredibly low.
Code_saturne v 5.1.3.
Turbulence model : RSM-EB
turbulent heat flux model : AFM
Kindly suggest few changes.
Regards
Kunal
- Attachments
-
- 63e4_afm_data0.csv
- (2.15 MiB) Downloaded 235 times
Re: Low value of Turbulent heat fluxes in Natural convection
Hi,
May I asked how did you implement Boussinessq?
i.e cs_user_source_terms.f90 under ustsnv or cs_user_physical_properties.f90 under usphyv?
Cheers,
Costas
May I asked how did you implement Boussinessq?
i.e cs_user_source_terms.f90 under ustsnv or cs_user_physical_properties.f90 under usphyv?
Cheers,
Costas
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
Re: Low value of Turbulent heat fluxes in Natural convection
Hello,
I have started a case with natural convection under 6.0-beta, Rij-SSG + AFM, the temperature / turbulent heat fluxes have a very high amplitude, 10e50. Switch back to 5.0 if you need results quickly. Otherwise, wait for fixes.
P.S.: tests performed with the temperature flagged as buoyant in cs_user_parameters.c and the SIMPLEC time scheme. Switching to PISO and nterup = 3 did not improve the situation.
Best regards,
Cédric
I have started a case with natural convection under 6.0-beta, Rij-SSG + AFM, the temperature / turbulent heat fluxes have a very high amplitude, 10e50. Switch back to 5.0 if you need results quickly. Otherwise, wait for fixes.
P.S.: tests performed with the temperature flagged as buoyant in cs_user_parameters.c and the SIMPLEC time scheme. Switching to PISO and nterup = 3 did not improve the situation.
Code: Select all
void
cs_user_parameters(cs_domain_t *domain)
{
/* Declare the temperature as a buoyant scalar */
cs_field_set_key_int(cs_field_by_name("temperature"),
cs_field_key_id("is_buoyant"),
1);
}
Cédric
-
- Posts: 9
- Joined: Fri Apr 27, 2018 2:33 pm
Re: Low value of Turbulent heat fluxes in Natural convection
Thank you for your feedback on the AFM model. Your temperature profile looks good probably because the GGDH part of the AFM works well whereas the additional terms ( bouyancy, velocity gradient,...) might be bugged.kunalck wrote: ↑Mon Apr 08, 2019 3:23 pm 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 problem regarding the turbulent heat flux magnitude. The turbulent heat fluxes for this case are incredibly low value(1e-38 to 1e-36). Whereas the temperature profiles are predicted well when compared to DNS results.
Also the turbulent kinetic energy (i.e. rij values) is also incredibly low.
Code_saturne v 5.1.3.
Turbulence model : RSM-EB
turbulent heat flux model : AFM
Kindly suggest few changes.
Regards
Kunal
Can you post your test case so that I can have a look ?
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
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
Re: Low value of Turbulent heat fluxes in Natural convection
Hello Jean-Francois,Jean-Francois Wald wrote: ↑Tue Apr 09, 2019 12:54 pmThank you for your feedback on the AFM model. Your temperature profile looks good probably because the GGDH part of the AFM works well whereas the additional terms ( bouyancy, velocity gradient,...) might be bugged.kunalck wrote: ↑Mon Apr 08, 2019 3:23 pm 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 problem regarding the turbulent heat flux magnitude. The turbulent heat fluxes for this case are incredibly low value(1e-38 to 1e-36). Whereas the temperature profiles are predicted well when compared to DNS results.
Also the turbulent kinetic energy (i.e. rij values) is also incredibly low.
Code_saturne v 5.1.3.
Turbulence model : RSM-EB
turbulent heat flux model : AFM
afm_dns.png
63e4_afm_data0.csv
Kindly suggest few changes.
Regards
Kunal
Can you post your test case so that I can have a look ?
Thank you for your feedback. I have also run the same case with unsteady solver(CS_v4.0) and I could not prevent the laminarization. Although, temperature profile was like previous case(predicting well compared with DNS).
Re: Low value of Turbulent heat fluxes in Natural convection
Hello,
I gave the value of gravity as a negative vector in GUI. And when using Boussinesq through cs_user_source_terms.f90, then the I would use gy for gravity.
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 plot of the temperature across the observed ring.
Also I would like to point out that the structures/ perturbations in temperature are observed only for Pr 0.025 cases. The rings/perturbations are no observed for Pr 0.7 or 0.006. Also I would like to add that the overall profile of the temperature on averaging across the domain along x-axis is slightly affected and can be observed only at a closer look. I just wanted to know what could be causing it and what can be done to avoid these in future ?
Heat flux model - AFM
Turbulence model - RSM:EB
unsteady RANS
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 plot of the temperature across the observed ring.
Also I would like to point out that the structures/ perturbations in temperature are observed only for Pr 0.025 cases. The rings/perturbations are no observed for Pr 0.7 or 0.006. Also I would like to add that the overall profile of the temperature on averaging across the domain along x-axis is slightly affected and can be observed only at a closer look. I just wanted to know what could be causing it and what can be done to avoid these in future ?
Heat flux model - AFM
Turbulence model - RSM:EB
unsteady RANS
-
- Posts: 4082
- Joined: Mon Feb 20, 2012 3:25 pm
Re: Low value of Turbulent heat fluxes in Natural convection
Hello,
Which version of the code are you using ? If you are not using version 6.0, I recommend upgrading, as you may
- either find a bug has been fixed
- help get a bug fixed (assuming it is a bug) in one of the next 6.0.z releases
while with older versions, only major bugs may be fixed (but if you have a regression relative to an older version, we do need to check that also).
Could you post your setup (without the mesh, unless it is small).
Best regards,
Yvan
Which version of the code are you using ? If you are not using version 6.0, I recommend upgrading, as you may
- either find a bug has been fixed
- help get a bug fixed (assuming it is a bug) in one of the next 6.0.z releases
while with older versions, only major bugs may be fixed (but if you have a regression relative to an older version, we do need to check that also).
Could you post your setup (without the mesh, unless it is small).
Best regards,
Yvan