[roenby]

IsoAdvector and interIsoFoam have been extended to work in regions including a porous medium.

This extension has now been integrated into the latest OpenFOAM-v2112 release. Please see release notes for further details:

https://www.openfoam.com/news/main-news/openfoam-v2112/solver-and-physics#interisofoam

Any feedback will be much appreciated.

Cheers,
Johan

[Bloerb]

Awesome. Have you ever thought about adding functionality which basically removes one phase from the calculation? So basically a zeroGradient at the interface.

In nearly all multi-phase applications only one phase is relevant for evaluation. And one major disadvantage of interfoam like solvers is, that they are failing for highly viscous fluids. Like a droplet of high viscosity silicon oil falling in air. You get extremly unrealistic high velocities in the air phase.

Or that one needs to use a compressible solver if one phase is incompressible (e.g water) and the other one is pretty much irrelevant for the user, but compressible. Like water waves, where the air has no effect on the solution, or like injection molding. Flow3D for example uses this and calls it TrueVOF. It's basically interFoam but the movement at the free surface is basically a zeroGradient boundary. And all the air phase cells are simply set to zero.

Since you are calculating normals at the interface and the cells surounding the interface in your isoAdvector aproach this would be an awesome and probably easily implemented addition. You'd simply need to see which cells are forming the interface. Set the velocity and pressure in those cells / extrapolate them, and afterwards zero out the solution in the other phase / set the matrix solution for those.