CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Boussinesq approximation And variable density model

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 2 Post By sbaffini
  • 1 Post By LuckyTran

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 8, 2020, 08:03
Default Boussinesq approximation And variable density model
  #1
New Member
 
saad
Join Date: Jan 2019
Posts: 7
Rep Power: 7
highend is on a distinguished road
I have a question about the Boussinesq approximation.
I know that the approximation is rho = rho0* ( 1+beta*(T-T0))
where rho0 is the reference density at T0 (reference temperature)
And beta is the thermal expansion coefficient.

My question is do T0 change with the system or is it constant?
Imagine you have a natural circulation system single-phase water system where the pressure is atmospheric pressure and you have a heat source and a sink(cooling).
If you give initial condition T0 = 300 ,rho0 = 1000, T0=T , rho0 = rho, beta =-0.0007
Beta will obviously be updated based on beta= - (1/rho0 )* delta rho/delta T

So the rho will change based on the equation but do I change T0 based on the previous time step value or keep it constant.

Also , can any body provide a variable density model to calculate the density based on temperature and pressure change other than the Boussinesq approximation? With reference literature please.

I read in Relap7 manual some equation where they used the equation of state to get the density change. There was a linear equation of state that seemed to be for water but they did not provide any reference to the equation source to get the correct parameters.

Any help will be much appreciated
highend is offline   Reply With Quote

Old   March 8, 2020, 11:02
Default
  #2
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,192
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
The very idea of the Boussinesq model is to keep T0, beta and rho0 constant, as it is a linearization around the 0 state. If you change stuff during iterations it is not Boussinesq anymore, you are actually using the true model in a certain sense.
FMDenaro and LuckyTran like this.
sbaffini is offline   Reply With Quote

Old   March 8, 2020, 13:52
Default
  #3
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,882
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by highend View Post
I have a question about the Boussinesq approximation.
I know that the approximation is rho = rho0* ( 1+beta*(T-T0))
where rho0 is the reference density at T0 (reference temperature)
And beta is the thermal expansion coefficient.

My question is do T0 change with the system or is it constant?
Imagine you have a natural circulation system single-phase water system where the pressure is atmospheric pressure and you have a heat source and a sink(cooling).
If you give initial condition T0 = 300 ,rho0 = 1000, T0=T , rho0 = rho, beta =-0.0007
Beta will obviously be updated based on beta= - (1/rho0 )* delta rho/delta T

So the rho will change based on the equation but do I change T0 based on the previous time step value or keep it constant.

Also , can any body provide a variable density model to calculate the density based on temperature and pressure change other than the Boussinesq approximation? With reference literature please.

I read in Relap7 manual some equation where they used the equation of state to get the density change. There was a linear equation of state that seemed to be for water but they did not provide any reference to the equation source to get the correct parameters.

Any help will be much appreciated



Other than the classical Bousinnesq model, have also a look to the low Mach formulation, classically used in cases of combustion.
FMDenaro is offline   Reply With Quote

Old   March 8, 2020, 20:03
Default
  #4
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,754
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
Quote:
Originally Posted by sbaffini View Post
The very idea of the Boussinesq model is to keep T0, beta and rho0 constant, as it is a linearization around the 0 state. If you change stuff during iterations it is not Boussinesq anymore, you are actually using the true model in a certain sense.

I just want to echo this statement. Yes, the point of Boussinesq is to keep all those those parameters constant. If you do any more (i.e. allow any of those to vary), you very quickly approach solving a variable density model. At that point you are no longer following the Boussinesq philosophy, you're just picking and choosing what dependencies & complexities you want to include in your model.



Quote:
Originally Posted by highend View Post
Also , can any body provide a variable density model to calculate the density based on temperature and pressure change other than the Boussinesq approximation? With reference literature please.

I'm confused by the question because these are really common. People use all the time the full ideal gas law, incompressible ideal gas law, polynomial functions of temperature and so on. What kind of novelties are you looking for that isn't readily available?
sbaffini likes this.
LuckyTran is offline   Reply With Quote

Old   March 9, 2020, 06:45
Default
  #5
Member
 
optimux's Avatar
 
Mingming Zhang
Join Date: Dec 2019
Posts: 35
Rep Power: 6
optimux is on a distinguished road
Yes, two papers I read a couple of months ago dealing with variable density flow, actually metal alloy. The densities of both solid and liquid change with temperature and composition. The density dependence on composition is very strong. Check out the follwoing 2 coupled papers: Numerical method for solution of strongly coupled binary alloy solidification problems, Daming Xu & Qingchun Li, 1991; Gravity-and solidification-shrinkage-induced liquid flow in a horizontally solidified alloy ingot, Daming Xu & Qingchun Li, 1991. Hope it helps.

Mingming
optimux is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Boussinesq approximation (temperature gradient range) azizi0407 FLUENT 3 December 11, 2023 13:05
Adding Boussinesq Approximation to multiphaseInterFoam? Alona OpenFOAM Programming & Development 2 November 20, 2019 19:43
Fire and smoke modeling CFDST CFX 84 April 30, 2019 08:01
Using rhok for a modified kepsilon model. How to access to the value? be_inspired OpenFOAM Programming & Development 1 October 19, 2017 13:46
Natural convection through fin in Ansys FLUENT hifjur FLUENT 9 March 16, 2016 12:52


All times are GMT -4. The time now is 21:48.