|
[Sponsors] |
January 14, 2014, 03:58 |
interFoam with boussinesq buoyant flow
|
#1 |
New Member
edward jr shin
Join Date: Sep 2013
Posts: 2
Rep Power: 0 |
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 TEqn.H fvScalarMatrix TEqn ( fvm::ddt(T) + fvm::div(phi, T) - fvm::laplacian(DTEff, T) ); UEqn.H fvc::reconstruct ( 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() pEqn.H surfaceScalarField phig ( ( fvc::interpolate(interface.sigmaK())*fvc::snGrad(a lpha1) // - ghf*fvc::snGrad(rho) + ghf*fvc::interpolate(rho)*twoPhaseProperties.TEf() *fvc::snGrad(T) )*rAUf*mesh.magSf() ); 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? |
|
January 14, 2014, 06:03 |
|
#2 |
Senior Member
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 21 |
Because rho is discontinuous.
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
mass flow in is not equal to mass flow out | saii | CFX | 12 | March 19, 2018 05:21 |
Question from tutorial: Buoyant flow in a partitioned cavity | noppawit | CFX | 13 | March 19, 2014 02:18 |
Slug Flow, interFoam, problems with Contact Angle | PrzemekPL | OpenFOAM Running, Solving & CFD | 13 | February 18, 2014 22:10 |
channel flow using interfoam in OpenFoam2.0 | bojiezhang | OpenFOAM | 2 | March 12, 2012 21:29 |
interFoam - stratified flow - problem with shear stress at interface | AnjaMiehe | OpenFOAM Running, Solving & CFD | 8 | June 14, 2010 06:49 |