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

buoyantBoussinesqSimpleFoam with porous zone (source term) problem

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

LinkBack Thread Tools Search this Thread Display Modes
Old   May 26, 2015, 09:05
Exclamation buoyantBoussinesqSimpleFoam with porous zone (source term) problem
New Member
Join Date: Feb 2014
Posts: 4
Rep Power: 12
cissystef is on a distinguished road
Hi everyone,

Apart from previous post using porousBafflePressure for a zero-thickness porous plate, I also wanted to test it out using DarcyForchheimer source term for a constructed porous zone.


So the geometry is the same, a duct (0.6*0.6*3) with a porous zone inside (0.6*0.6*0.1).

0 folder - initial condition:

the most important 3 parameters for setup:

U: walls (fixedValue; uniform (0 0 0));
inlet (fixedValue; uniform (0.88 0 0));
outlet (zeroGradient);

T: walls and outlet (zeroGradient)
inlet (fixedValue, uniform 300);

p_rgh: walls and inlet (fixedFluxPressure, uniform 0)
outlet (fixedValue, uniform 0)

System folder:

However, I found out that the results are not consistent with the experiment The DarcyForchheimerCoeffs I set in fvOptions are obtained from Delta p = -(D*mu*U + 0.5*I*rho*|U|^2 )L:

d d [0 -2 0 0 0 0 0] (6.3e4 0 0);
f f [0 -1 0 0 0 0 0] (17.73 0 0);

type cartesian;
origin (0 0 0);
type axesRotation;
e1 (1 0 0);
e2 (0 1 0);

My flow direction is along x-axis, and my porosity directions (not quite sure about this) set for e1 e2 are as indicated above. I did change e1 e2 to a different orientation, the results are still quite similar.

I used snappyHexMesh to refine the porous zone, and topoSet to create cell zones for porous block.

My problem is again, results (pressure drop across the porous zone) are not matching with the experiment data. The pressure drop created by using source term is 1/10 of the experiment value, no matter what D value you choose. The convergence plot is good, although T residuals is above 0.1.

Can someone explain what might cause this problem?

Thank you very much.
cissystef is offline   Reply With Quote


openfoam 2.1.x, porous boundary, porous domain

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
[swak4Foam] swak4foam building problem GGerber OpenFOAM Community Contributions 54 April 24, 2015 17:02
Problem compiling a custom Lagrangian library brbbhatti OpenFOAM Programming & Development 2 July 7, 2014 12:32
[swak4Foam] Error bulding swak4Foam sfigato OpenFOAM Community Contributions 18 August 22, 2013 13:41
problem with pump as source term hari FLUENT 1 September 19, 2005 12:38
UDF Scalar Code: HT 1 Greg Perkins FLUENT 8 October 20, 2000 13:40

All times are GMT -4. The time now is 03:59.