Problem with rhoPorousMRFSimpleFoam and how to define porous media
 

Dear All,
I am facing problem while running rhoPorousSimpleFoam solver in OpenFoam 2.0.x. I am trying to simulate pressure driven flows (Pressure is fixed boundary condition at inlet and outlet) in a simple pipe with porous zone. The K-Epsilon turbulence model is used. The solver is oscillating with higher values of epsilon on boundary. I am also facing problem in defining Porous media. If some one have some document regarding to that please send to me.
In some cases it runs by making so many changes in relaxation parameters etc. But the results are not validated, especially the pressure drop across the porous media both for isotropic and anisotropic cases. I am copying below the porous zone file defined for isotropic(homogeneous) porous media.


1
(
cellTable_8
{
coordinateSystem
{
e1 (0.90709 0.295279 0.299993);
e2 (-0.329827 0.033477 0.943447);
}

Darcy
{

d d [0 -2 0 0 0 0 0] (1.14e8 -1 -1);

f f [0 -1 0 0 0 0 0] (422.955 422.955 422.955);
}
thermalModel
{
type none; // fixedTemperature;

// fixedTemperature coefficients
T 350;
}
}

)

Amit,

Did you solve your problem.

I am assuming that you have a high relaxation factor for epsilon. Decreasing it might help. If not, try to change the initial values of epsilon.

Also, to crenate a porousZone from a mesh not built in foam, you need to create cellSet and then convert it to a cellZone using setSet function.

Kalyan

Meaning of value in porous zone file
 

Hi can any one tell me what is ment by 1 in porouse zone file. Is it default value of porosity. or it mean one cell zone or one block.

Thanks

Hi mech01

It is the latter.

Latter means its one block

Absolutely...it means you are defining only one zone i.e cellTable_8

Thank you Akshay