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

interFoam with boussinesq buoyant flow

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

LinkBack Thread Tools Display Modes
Old   January 14, 2014, 04:58
Default interFoam with boussinesq buoyant flow
New Member
edward jr shin
Join Date: Sep 2013
Posts: 2
Rep Power: 0
join4us is on a distinguished road
Hi. there.

I want to solve a hotpot simulation.

The pot has water and air in half on pot.

Base solver is interFoam.

I modified twoPhaseMixture.C and twoPhaseMixture.H for thermal expansion and thermal diffusive.
and add
fvScalarMatrix TEqn
+ fvm::div(phi, T)
- fvm::laplacian(DTEff, T)

fvc::interpolate(interface.sigmaK())*fvc::snGrad(a lpha1)
// - ghf*fvc::snGrad(rho)
+ ghf*fvc::interpolate(rho)*twoPhaseProperties.TEf() *fvc::snGrad(T)
- fvc::snGrad(p_rgh)
) * mesh.magSf()

surfaceScalarField phig
fvc::interpolate(interface.sigmaK())*fvc::snGrad(a lpha1)
// - ghf*fvc::snGrad(rho)
+ ghf*fvc::interpolate(rho)*twoPhaseProperties.TEf() *fvc::snGrad(T)

The simulation gives good result except pressure field.
p=p_rgh + rho*gh

I attached flow field pictures.

Why pressure field look like step?

Could you help me to solve the problem well?
Attached Images
File Type: png alpha.png (13.6 KB, 30 views)
File Type: png p_rgh.png (41.1 KB, 33 views)
File Type: png p.png (19.2 KB, 36 views)
File Type: png T.png (13.6 KB, 27 views)
File Type: png U.png (72.7 KB, 36 views)
join4us is offline   Reply With Quote

Old   January 14, 2014, 07:03
Senior Member
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 15
Bernhard is on a distinguished road
Because rho is discontinuous.
Bernhard is offline   Reply With Quote


Thread Tools
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 On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Question from tutorial: Buoyant flow in a partitioned cavity noppawit CFX 13 March 19, 2014 03:18
Slug Flow, interFoam, problems with Contact Angle PrzemekPL OpenFOAM Running, Solving & CFD 13 February 18, 2014 23:10
channel flow using interfoam in OpenFoam2.0 bojiezhang OpenFOAM 2 March 12, 2012 22:29
interFoam - stratified flow - problem with shear stress at interface AnjaMiehe OpenFOAM Running, Solving & CFD 8 June 14, 2010 06:49
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07

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