|
[Sponsors] |
Boussinesq approximation And variable density model |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 8, 2020, 08:03 |
Boussinesq approximation And variable density model
|
#1 |
New Member
saad
Join Date: Jan 2019
Posts: 7
Rep Power: 7 |
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 |
|
March 8, 2020, 11:02 |
|
#2 |
Senior Member
|
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.
|
|
March 8, 2020, 13:52 |
|
#3 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,850
Rep Power: 73 |
Quote:
Other than the classical Bousinnesq model, have also a look to the low Mach formulation, classically used in cases of combustion. |
||
March 8, 2020, 20:03 |
|
#4 | ||
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,747
Rep Power: 66 |
Quote:
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:
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? |
|||
March 9, 2020, 06:45 |
|
#5 |
Member
Mingming Zhang
Join Date: Dec 2019
Posts: 35
Rep Power: 6 |
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 |
|
|
|
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 |