CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   porousInterFoam Solver (

sprobst76 July 15, 2010 10:29

porousInterFoam Solver
Dear All,

I am a new user of OpenFOAM and not really experienced yet. I am interested in the simulation the flow through porous media. For my cases the flow is in general unsaturated and therefore a solver with two phases would be very useful. In the latest version of OpenFOAM there is a new solver called porousInterFoam which can perhaps help me with my problems. But there is no example and I have also not found a description for this solver.

Can anybody give me more informations about the solver.

Thank you in advance!

Best regards

andrea.pasquali July 23, 2010 10:20

In the new OpenFOAM 1.7 there is the porousInterFoam Solver.
There is not a tutoria but you can take the interFoam tutorial adding a porous medium inside


tabaer February 25, 2011 09:43

The details of porousInterFoam are also of interest to me.

From a cursory look at the source, it seems that this solver is applicable to problems in which a fully saturated fluid wave is being forced through a porous material. Accounting for capillary wetting in the porous material via a saturation-capillary pressure relation or the influence of saturation level on permeability do not seem to be included in this solver.

Has anyone had any experience with this solver that could comment on these points ?


julienreveillon April 29, 2014 07:11

For those who are still looking for a tutorial as I was, I gave a porous project to some of my students and after looking (and not finding) for a while, I finally wrote a small document. Only small problem a small variable name (mu) recognition problem. A working solution is detailed here:

John Degenkolb September 3, 2014 09:56

Dear all,

I'd like to use the porousInterFoam Solver for a laminar, creeping flow though porous media i.e. simulating Darcy's law only.

It works quite fine, but for bigger meshes, the calculation time becomes very long. So my question is how I can tune the solver to make it quicker and more tailored to my problem.

I already outcommented (dumb as I am) the line

// + fvm::div(rhoPhi, U)

in the UEqn.H file.
The computation still works and the results are the same, but it only saves arround 13 % of time.
It seems, that the guys who do what I do on a professional basis, use a much simpler algorithm:


Best regards,


All times are GMT -4. The time now is 23:39.