[physicosm]

Dear foamers,


recently, I had to run some multiphase porous simulations and I realized that there isn't any tutorial case regarding the porousInterFoam. Therefore, I've decided to create one, at least for foam-extend 4.0, in the case that somebody needs it in the future.


Before creating the case, I recompiled the solver so, it will be independent from interFoam. However, a minor change also should be made. In the UEqn.H, the volScalarField mu is defined such as

....

    // Calculate and cache mu for the porous media
    volScalarField mu(twoPhaseProperties.mu());
....

This is not a global definition and thus, it cannot be called later when the solver adds resistance in the UEqn due to porosity:

...

    pZones.addResistance(UEqn);
...

Thus, mu should be defined globally in the createFields.H as,

    // Need to store mu
    volScalarField mu
    (
        IOobject
        (
            "mu",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        twoPhaseProperties.mu()
    );

I have attached the solver named porosityInterFoam and it runs in foam-extend 4.0.
porosityInterFoam.zip


Further, a case has been also attached and this describes the two-phase flow in a channel with a porous media.
porosityInterFoam-case.zip


Here are some results: results.zip



Enjoy freely!


Thomas N.