|
[Sponsors] |
buoyantBoussinesqSimpleFoam with porous zone (source term) problem |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
May 26, 2015, 08:05 |
buoyantBoussinesqSimpleFoam with porous zone (source term) problem
|
#1 |
New Member
Join Date: Feb 2014
Posts: 4
Rep Power: 12 |
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. Geometry: 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); coordinateSystem { type cartesian; origin (0 0 0); coordinateRotation { 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. |
|
Tags |
openfoam 2.1.x, porous boundary, porous domain |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] swak4foam building problem | GGerber | OpenFOAM Community Contributions | 54 | April 24, 2015 16:02 |
Problem compiling a custom Lagrangian library | brbbhatti | OpenFOAM Programming & Development | 2 | July 7, 2014 11:32 |
[swak4Foam] Error bulding swak4Foam | sfigato | OpenFOAM Community Contributions | 18 | August 22, 2013 12:41 |
problem with pump as source term | hari | FLUENT | 1 | September 19, 2005 11:38 |
UDF Scalar Code: HT 1 | Greg Perkins | FLUENT | 8 | October 20, 2000 12:40 |