CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Porous media model in OpenFOAM 14

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   July 2, 2007, 02:47
Default The release notes for v. 1.4 r
  #1
gac
New Member
 
gac
Join Date: Mar 2009
Location: Pennsylvania, USA
Posts: 2
Rep Power: 0
gac is on a distinguished road
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.
gac is offline   Reply With Quote

Old   July 2, 2007, 03:46
Default In the tutorials directory, yo
  #2
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: http://olesenm.github.io/
Posts: 777
Rep Power: 18
olesen will become famous soon enough
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.
olesen is offline   Reply With Quote

Old   July 2, 2007, 05:12
Default You can find class files in:
  #3
New Member
 
Masato Otsuki
Join Date: Mar 2009
Location: Tokyo, Japan
Posts: 26
Rep Power: 8
otsuki is on a distinguished road
You can find class files in:
OpenFOAM-1.4/src/finiteVolume/cfdTools/general/porousMedia
You can find brief comments in "porousZones.H".
otsuki is offline   Reply With Quote

Old   July 4, 2007, 08:03
Default Hi When I try to run the a
  #4
New Member
 
abhishek k n
Join Date: Mar 2009
Location: Gothenburg, Sweden
Posts: 16
Rep Power: 8
knabhishek is on a distinguished road
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
knabhishek is offline   Reply With Quote

Old   July 4, 2007, 09:13
Default Do you get it to work without
  #5
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: http://olesenm.github.io/
Posts: 777
Rep Power: 18
olesen will become famous soon enough
Do you get it to work without FoamX?
olesen is offline   Reply With Quote

Old   July 4, 2007, 09:41
Default No I can't get it to run witho
  #6
New Member
 
abhishek k n
Join Date: Mar 2009
Location: Gothenburg, Sweden
Posts: 16
Rep Power: 8
knabhishek is on a distinguished road
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??
knabhishek is offline   Reply With Quote

Old   July 4, 2007, 10:35
Default The points are missing because
  #7
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,912
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
The points are missing because you didn't run blockMesh yet
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   August 1, 2007, 14:24
Default Hi, i've just downloaded and
  #8
stchouan
Guest
 
Posts: n/a
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.
  Reply With Quote

Old   August 1, 2007, 17:34
Default Hi Stephane, Compressible t
  #9
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
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.
eugene is offline   Reply With Quote

Old   September 12, 2007, 07:35
Default Hi When I try to run the an
  #10
Member
 
Paul Mauk
Join Date: Mar 2009
Posts: 39
Rep Power: 8
plmauk is on a distinguished road
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 is offline   Reply With Quote

Old   September 12, 2007, 08:23
Default I've tried to run this case wi
  #11
Member
 
Paul Mauk
Join Date: Mar 2009
Posts: 39
Rep Power: 8
plmauk is on a distinguished road
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.
plmauk is offline   Reply With Quote

Old   September 12, 2007, 21:35
Default Hi Paul, The swich pressure
  #12
New Member
 
Masato Otsuki
Join Date: Mar 2009
Location: Tokyo, Japan
Posts: 26
Rep Power: 8
otsuki is on a distinguished road
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
otsuki is offline   Reply With Quote

Old   September 13, 2007, 08:22
Default Hi, Masato! It seems to work.
  #13
Member
 
Paul Mauk
Join Date: Mar 2009
Posts: 39
Rep Power: 8
plmauk is on a distinguished road
Hi, Masato!
It seems to work.
Thank you very mutch!
plmauk is offline   Reply With Quote

Old   October 30, 2007, 07:01
Default Hi, everyone! The case don't
  #14
Member
 
Paul Mauk
Join Date: Mar 2009
Posts: 39
Rep Power: 8
plmauk is on a distinguished road
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?
plmauk is offline   Reply With Quote

Old   October 30, 2007, 07:55
Default To Paul Mauk: you should run b
  #15
Senior Member
 
mkraposhin's Avatar
 
Matvej Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 172
Rep Power: 8
mkraposhin is on a distinguished road
To Paul Mauk: you should run blockMesh before running case
mkraposhin is offline   Reply With Quote

Old   October 31, 2007, 07:19
Default To Matvej Karposhin: Thank you
  #16
Member
 
Paul Mauk
Join Date: Mar 2009
Posts: 39
Rep Power: 8
plmauk is on a distinguished road
To Matvej Karposhin: Thank you very match for helping. Would you like to share some experience in Open Foam with me? In russian?
plmauk is offline   Reply With Quote

Old   October 31, 2007, 07:48
Default Paul Mauk, ii would be nice, i
  #17
Senior Member
 
mkraposhin's Avatar
 
Matvej Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 172
Rep Power: 8
mkraposhin is on a distinguished road
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?
mkraposhin is offline   Reply With Quote

Old   October 31, 2007, 08:07
Default Of course, I can speak russian
  #18
Member
 
Paul Mauk
Join Date: Mar 2009
Posts: 39
Rep Power: 8
plmauk is on a distinguished road
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.
plmauk is offline   Reply With Quote

Old   November 1, 2007, 06:57
Default Paul Mauk, here is my e-mail:
  #19
Senior Member
 
mkraposhin's Avatar
 
Matvej Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 172
Rep Power: 8
mkraposhin is on a distinguished road
Paul Mauk, here is my e-mail: mkraposhin@inbox.ru, same as in profile.
mkraposhin is offline   Reply With Quote

Old   September 7, 2009, 04:29
Default LHS and RHS of + have different dimensions
  #20
New Member
 
srikara's Avatar
 
Srikara Mahishi
Join Date: Mar 2009
Location: Bangalore
Posts: 22
Rep Power: 8
srikara is on a distinguished road
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: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
srikara is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
CFX: what is "A true volume-porous media model?" whiz CFX 17 June 6, 2013 00:19
Porous media in openfoam arturo OpenFOAM Running, Solving & CFD 31 March 16, 2013 17:07
Power law model for porous media Karl FLUENT 0 December 13, 2004 16:14
turbulence model in porous media zhzhguo FLUENT 2 November 1, 2004 08:10
Question about porous media model Glenn Price FLUENT 5 June 28, 2000 11:48


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