
[Sponsors] 
buoyantPimpleFoam with boussinesq is incompressible? 

LinkBack  Thread Tools  Search this Thread  Display Modes 
September 15, 2020, 06:48 
buoyantPimpleFoam with boussinesq is incompressible?

#1 
New Member
Visakh
Join Date: Dec 2010
Location: Bangalore, India
Posts: 7
Rep Power: 12 
Hello everyone,
I am trying to solve 2D RayleighBenard convection with OpenFOAM 7. I am using buoyantPimpleFoam with boussinesq model for equation of state. Does this make the solver incompressible? Does the solution differ in any way with that of the (deprecated) buoyantBoussinesqPimpleFoam? Also, does the value of molWeight play any role in the solution? Because I'm trying to get a Nu vs Ra, Pr correlation, my only input parameters are Ra and Pr, or indirectly, mu, Cp, rho0, beta, Pr and deltaT, L from the boundary conditions. How does sensibleInternalEnergy instead of sensibleEnthalpy make a difference in this context? And finally, is there a better way to solve nondimensional version of equations (than solving th usual dimensional equations by adjusting the fluid properties to get the desired Ra, Pr)? Attaching my thermophysicalProperties dictionary: Code:
thermoType { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState Boussinesq; specie specie; energy sensibleEnthalpy; } mixture { specie { molWeight 18; } equationOfState { rho0 996.5; T0 300; beta 3.03e04; } thermodynamics { Cp 4195; Hf 0; } transport { mu 0.8e2; Pr 6.0; } } 

September 15, 2020, 17:49 

#2  
Member
Join Date: Oct 2017
Posts: 74
Rep Power: 5 
Some answers:
Quote:
It is used to calculate the specific gas constant. 

September 18, 2020, 08:04 

#3 
New Member
Visakh
Join Date: Dec 2010
Location: Bangalore, India
Posts: 7
Rep Power: 12 

September 18, 2020, 11:55 

#4 
Senior Member
Join Date: Sep 2013
Posts: 269
Rep Power: 16 
nMoles and molWeight is not used for pureMixture. It is however used for the other options and therefore a parameter that needs to be present although it is not used in your case. Yeah c++ inheritance thingy. They got rid of some parts of it over the years though.
There was a debate about buoyantPimpleFoam and buoyantBoussinesqPimpleFoam and their differences on here a while ago. Since the Boussinesq approximation was incorporated into the equation of state it is not really incompressible if i remember correctly. However i am not certain if this is how it was implemented. You might find the thread though. What is true though is that the density is not a function of pressure. 

September 18, 2020, 18:21 

#5  
Member
Join Date: Oct 2017
Posts: 74
Rep Power: 5 
Quote:
Quote:
The German Wikipedia article on incompressible fluid contains the following statement: Quote:
Quote:


September 24, 2020, 03:24 

#6 
Senior Member
Join Date: Sep 2013
Posts: 269
Rep Power: 16 
Yes that is correct, sorry for the confusion. for rho(p,T) or rho(p) you get compressible for rho(T) incompressible behaviour. In other words, without pressure dependence you can't compress your volume by pushing on it. You can however change it's volume by heating it.
The topic i meant was the following. Wrong implementation of the Boussinesq assumption in heat transfer solvers the main argument here is probably that the density update shouldn't happen in all terms. I am however not familiar with the approximation and implementation to know which is the case. 

September 24, 2020, 03:44 

#7 
Senior Member
Agustín Villa
Join Date: Apr 2013
Location: Fuenlabrada
Posts: 280
Rep Power: 12 
This is a nice thread as I see, I am wondering the same and I built my own Boussinesq solver based on older OF versions.
Let's assume that buoyantSimple/PimpleFoam with Boussinesq is incompressible. What are the turbulence models employed on these simulations? Boussinesq solvers use incompressible turbulence models, but I assume that with this modification you're obliged use compressible models. Most of the incompressible models for Reynolds stresses can be "translated into" compressible models, but if I remember correctly, only the turbulent Prandtl number is used for the turbulent heat fluxes. If you are using custom models for heat fluxes (THFM) , then some coding is needed to pass your models from the incompressible to the compressible framework. This is why I'd build my own Boussinesq solvers until I find out a good way to add my custom THFM models 

Tags 
boussinesq approximation, buoyanboussinesqpimple, rayleighbenard 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
[IHFOAM] The IHFOAM Thread  Phicau  OpenFOAM Community Contributions  360  Yesterday 12:37 
[swak4Foam] swakExpression not writing to log  alexfells  OpenFOAM Community Contributions  3  March 16, 2020 18:19 
buoyantPimpleFoam and incompressible fluid (temperature dependent properties)  agustinvo  OpenFOAM Running, Solving & CFD  2  May 26, 2016 17:33 
1.7.x > buoyantPimpleFoam > hRhoThermo > incompressible and icoPoly3ThermoPhysics?  will.logie  OpenFOAM Programming & Development  1  February 16, 2011 20:52 
1.7.x > buoyantPimpleFoam > hRhoThermo > incompressible and icoPoly3ThermoPhysics?  will.logie  OpenFOAM  0  December 16, 2010 07:08 