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!
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
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 ?
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: http://www.cfdandco.com/page13/page20/index.html
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
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:
|All times are GMT -4. The time now is 09:05.|