CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Porous media model in OpenFOAM 14 (http://www.cfd-online.com/Forums/openfoam-solving/59343-porous-media-model-openfoam-14-a.html)

gac July 2, 2007 02:47

The release notes for v. 1.4 r
 
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 July 2, 2007 03:46

In the tutorials directory, yo
 
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 July 2, 2007 05:12

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

knabhishek July 4, 2007 08:03

Hi When I try to run the a
 
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 July 4, 2007 09:13

Do you get it to work without
 
Do you get it to work without FoamX?

knabhishek July 4, 2007 09:41

No I can't get it to run witho
 
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 July 4, 2007 10:35

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

stchouan August 1, 2007 14:24

Hi, i've just downloaded and
 
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 August 1, 2007 17:34

Hi Stephane, Compressible t
 
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 September 12, 2007 07:35

Hi When I try to run the an
 
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 September 12, 2007 08:23

I've tried to run this case wi
 
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 September 12, 2007 21:35

Hi Paul, The swich pressure
 
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 September 13, 2007 08:22

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

plmauk October 30, 2007 07:01

Hi, everyone! The case don't
 
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 October 30, 2007 07:55

To Paul Mauk: you should run b
 
To Paul Mauk: you should run blockMesh before running case

plmauk October 31, 2007 07:19

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

mkraposhin October 31, 2007 07:48

Paul Mauk, ii would be nice, i
 
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 October 31, 2007 08:07

Of course, I can speak russian
 
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 November 1, 2007 06:57

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

srikara September 7, 2009 04:29

LHS and RHS of + have different dimensions
 
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:

Quote:

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::printStack(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::operator+(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::operator+<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


All times are GMT -4. The time now is 15:17.