CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

InterFoam with boussinesq approximation

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   March 28, 2007, 10:35
Default Hello, I've been trying to
Posts: n/a

I've been trying to implement the boussinesq approximation into two phase solver interFoam. I have implemented density as a temperature dependent function in the twoPhaseMixture class and added
+ fvm::div(phi,T)
- fvm::laplacian(DT, T)

as the "energy" equation to be solved within the PISO-loop.

As the test case I use a right circular cylinder half filled with the liquid as a axis symmetric simulation (same setup has been used for isothermal simulation and working perfect). Gravity is acting along the symmetry axis and the contact between the liquid and wall is 90°, so the liquid is at rest and no motion is induced by capillary forces.

I'm encountering a increase of the liquid fraction during the simulation when the wall that touches the contact line heats up the liquid.
To verify the cause of the phenomena I did some tests. When gravity is turned off and only diffusive terms play a role the liquid/gas ratio stays constant and one can observe the diffusive temperature transport in the liquid and gas. With gravity turned on and only heating the part of the wall that is in contact with the gas produces nice convection/diffusion transport in the gas phase and no change in the liquid/gas ratio is observed.

I'm guessing that the additional density gradient near the surface accounts for additional buoyancy with is not balanced by sigma*k (surface tension * curvature) which causes some flux normal to the liquid/gas interface which is then adjusted by the gamma equation, causing extra liquid volume.

I've tried to adjust the energy equation to account for no convective fluxes across the liquid/gas interface in same manner as in the gamma equation, but without any luck.

Unfortunately I'm out of ideas how to overcome this problem. Did anyone already implement two phase flow with boussinesq approximation?

I would real appreciate if anyone could give me a hint where to look at or point me to some papers that discuss the implementation of the boussinesq approximation in a two phase VOF-code.


  Reply With Quote

Old   November 24, 2016, 21:07
Join Date: Jan 2016
Location: china
Posts: 38
Rep Power: 7
wangchengan2003 is on a distinguished road
Send a message via Skype™ to wangchengan2003
Dear Malte,

I meet the same problem like you. If you have found the solution?


wangchengan2003 is offline   Reply With Quote

Old   October 25, 2020, 13:21
New Member
Join Date: Nov 2019
Location: Brazil
Posts: 5
Rep Power: 3
chymalaia is on a distinguished road
Hi, Chengan

I also have to implement the boussinesq approximation in a two phase VOF-code. Because my project consists in modeling a Stratified Rayleigh-Bénard-Poiseuille flow and its interface... Has anyone found a solution to this type of problem yet?


chymalaia is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Abdul Aziz Jaafar Main CFD Forum 19 May 3, 2019 12:00
Boussinesq approximation Joshua Siemens 0 July 9, 2008 12:13
boussinesq approximation jamal FLUENT 2 March 25, 2008 09:57
Boussinesq approximation Volker Phoenics 0 September 5, 2003 16:48
Boussinesq approximation again Gabriel Main CFD Forum 3 May 11, 2000 10:24

All times are GMT -4. The time now is 17:51.