[gac]

The release notes for v. 1.4 refer to a "Porous media model including power-law inertial and viscous models for either explicit or implicit implementation in any of the pressure-velocity solvers (demonstration example included)."

However, there is no mention of the porous media model in the User Guide, and I can't find the "demonstration example."

Any assistance is appreciated.

[olesen]

In the tutorials directory, you should find both
rhoExplicitPorousSimpleFoam and rhoImplicitPorousSimpleFoam with an angleDuct example.

For higher flow resistances, the implicit formulation should be more robust. I don't think that it costs much more either.

[otsuki]

You can find class files in:
OpenFOAM-1.4/src/finiteVolume/cfdTools/general/porousMedia
You can find brief comments in "porousZones.H".

[knabhishek]

Hi

When I try to run the angleDuct tutorial the following error is reported.

FoamXError "getApplication::Invalid application class name 'rhoExplicitPorousSimpleFoam'."
In function "IPropertiesImpl"
in file "IPropertiesImpl.C" at line 910

[olesen]

Do you get it to work without FoamX?

[knabhishek]

No I can't get it to run without FoamX: in fact I don't know how to create a case from the terminal using solvers that don't have a tutorial associated with them (like dieselEngineFoam). I have so far been using the drop down menu in FoamX having only access to the default solvers upon creating a case. Maybe knowing how to do this will solve my problem.

When i run the case angleDuct on terminal I am using this command:

rhoPorousSimpleFoam ~/ angledDuctImplicit

the following error is reported:

FOAM FATAL ERROR : Cannot find file "points" in directory "constant/polyMesh"

From function Time::findInstance(const word& dir, const word& name) in file db/Time/findInstance.C at line 133.

I have aslo noticed there is no solver with name rhoImplicitPorousSimpleFoam and rhoExplicitPorousSimpleFoam

But there is a rhoPorousSimpleFoam. Maybe I need to use rhoExplicitPorousSimpleFoam in conjunction with something??

[gschaider]

The points are missing because you didn't run blockMesh yet

Hi,
i've just downloaded and installed the last version of OpenFoam and i was wondering if anybody had any experience with it on solving coupled and non linear convection-diffusion problems like those in reservoir or bassin simulation in the oil industry. Particularly with heterogenous and anisotropic porous media, has anyone ever used OpenFoam to simulate compressible two-phase flows in porous media?
Thanks for your answers.

[eugene]

Hi Stephane,

Compressible two-phase flow has been done for transonic sprays, but I am not aware of any two-phase flow calculations in conjunction with porous media. To get it working would require some coding and depending on your working fluids, the possible implementation of additional equations of state. You can mail me directly if you require additional info.

[plmauk]

Hi

When I try to run the angleDuct tutorial the following error is reported:

FoamXError "getApplication::Invalid application class name 'rhoExplicitPorousSimpleFoam'."

I know, the same question was asked by "abhishek kn"
but I couldnt find any answer...

[plmauk]

I've tried to run this case without FoamX,
but it did not work. I've got following message:

FOAM FATAL IO ERROR : keyword pressureImplicitPorousity is undefined in dictionary "/home/cfduser1/OpenFOAM/cfduser1-1.4.1/run/tutorials/rhoExplicitPorousSimpleFoa m/angledDuctExplicit/system/fvSolution::SIMPLE"

file: /home/cfduser1/OpenFOAM/cfduser1-1.4.1/run/tutorials/rhoExplicitPorousSimpleFoam /angledDuctExplicit/system/fvSolution::SIMPLE from line 75 to line 76.

From function dictionary::lookupEntry(const word& keyword) const
in file db/dictionary/dictionary.C at line 146.

[otsuki]

Hi Paul,

The swich pressureImplicitPorousity has no effect
for rhoExplicitPorousSimpleFoam.
Modify angledDuctExplicit/system/fvSolution as:

SIMPLE
{
nNonOrthogonalCoccectors 0;
pMin pMin [1 -1 -2 0 0 0 0] 100;
pressureImplicitPorousity No;
}

and try again.
Masato

[plmauk]

Hi, Masato!
It seems to work.
Thank you very mutch!

[plmauk]

Hi, everyone!
The case don't work once again.

FOAM FATAL ERROR : Cannot find file "points" in directory "constant/polyMesh"

From function Time::findInstance(const word& dir, const word& name)
in file db/Time/findInstance.C at line 133.

FOAM exiting

Can anybody help me please?

[mkraposhin]

To Paul Mauk: you should run blockMesh before running case

[plmauk]

To Matvej Karposhin: Thank you very match for helping. Would you like to share some experience in Open Foam with me? In russian?

[mkraposhin]

Paul Mauk, ii would be nice, if you can speak russian language (english is still very difficult to me).
What kind of problems are you intersted in OpenFoam?

[plmauk]

Of course, I can speak russian much better than english. I think it would be better, you tell me your e-mail adress and I will write an e-mail, to discribe my problems.

[mkraposhin]

Paul Mauk, here is my e-mail: mkraposhin@inbox.ru, same as in profile.

[srikara]

Hi,
I am running the rhoPorousSimpleFoam solver to solve a porous media simulation in OF 1.5. I get the following error message when I enter the command:

Starting time loop

Time = 0.1



LHS and RHS of + have different dimensions
dimensions : [0 2 -1 0 0 0 0] + [1 -1 -1 0 0 0 0]
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/smahishi/OpenFOAM/OpenFOAM
-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/smahishi/OpenFOAM/OpenFOAM-1.5/lib/linux64Gcc
DPOpt/libOpenFOAM.so"
#2 Foam:perator+(Foam::dimensionSet const&, Foam::dimensionSet const&) in "/h
ome/smahishi/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#3 Foam::tmp<Foam::GeometricField<Foam::typeOfSum<dou ble, double>::type, Foam::
fvPatchField, Foam::volMesh> > Foam:perator+<double, double, Foam::fvPatchFiel
d, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foa m::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::vol Mesh> const&) in "/home/smahishi/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libco mpressibleRASModels.so"
#4 Foam::compressible::RASModel::muEff() const in "/home/smahishi/OpenFOAM/Open FOAM-1.5/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#5 Foam::compressible::RASModels::kEpsilon::divDevRho Reff(Foam::GeometricField< Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/home/smahi shi/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#6 main in "/home/smahishi/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPO pt/rhoPorousSimpleFoam"
#7 __libc_start_main in "/lib64/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/home/smahishi/OpenFOAM/OpenFOAM-1.5 /applications/bin/linux64GccDPOpt/rhoPorousSimpleFoam"


From function operator+(const dimensionSet& ds1, const dimensionSet& ds2)
in file dimensionSet/dimensionSet.C at line 403.

FOAM aborting

Aborted

Does anybody know where this addition is taking place? It has something to do with muEff, that much I understand. Could anybody please let me know what I need to modify here?

Thank you in advance,
Srikara