[Amit Mathur]

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;
}
}

)

[kalyangoparaju]

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

[mech01]

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

[Akshay]

Hi mech01

It is the latter.

[mech01]

Originally Posted by Akshay

Hi mech01

It is the latter.

Latter means its one block

[Akshay]

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

[mech01]

Originally Posted by Akshay

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

Thank you Akshay