[AmirBaqa1987]

Hi everyone

I want to set up a simple case where I have flow through a pipe and the pipe is filled with porous media(metal foam).

my case is incompressible,laminar and steady state so I've selected porousSimpleFoam solver.
I modify blockMeshDict like below:

convertToMeters 1;

vertices
(
    (0 0 0)      //0
    (0 0 0.001)
    (0 0.001 0.001)      //2
    (0 0.001 0)
    (0 0 0.003)   //4
    (0 0.0021213 0.0021213)
    (0 0.003 0 )      //6

    (1 0 0)      //7
    (1 0 0.001)
    (1 0.001 0.001)      //9
    (1 0.001 0)
    (1 0 0.003)   //11
    (1 0.0021213 0.0021213)
    (1 0.003 0 )      //13
    

);

blocks
(
    // inlet block
    hex (7 8 9 10 0 1 2 3) porosity (9 9 300) simpleGrading (1 1 1)
    hex (8 11 12 9 1 4 5 2) porosity (51 9 300) simpleGrading (0.03075 1 1)
    hex (10 9 12 13 3 2 5 6) porosity (9 51 300) simpleGrading (1 0.03075 1)
);

edges
(
    arc 4 5 (0 0.001148 0.00277)
    arc 5 6 (0 0.00277 0.001148)
    arc 11 12 (1 0.001148 0.00277)
    arc 12 13 (1 0.00277 0.001148)
);
boundary
(
    inlet
    {
        type patch;
        faces
        (
            (0 1 2 3)
            (1 4 5 2)
            (2 5 6 3)
        );
    }

    outlet
    {
        type patch;
        faces
        (
            (7 10 9 8)
            (8 9 12 11)
            (9 10 13 12)
        );
    }

    side1
    {
        type symmetryPlane;
        neighbourPatch side2;
        faces
        (
            (0 7 8 1)
            (1 8 11 4) 
        );

       transform rotational;
        rotationAxis (1 0 0);
        rotationCentre (0 0 0);
    }

    side2
    {
        type symmetryPlane;
        neighbourPatch side1;
        faces
        (
            (0 3 10 7)
            (3 6 13 10)
        );

        transform rotational;
        rotationAxis (1 0 0);
        rotationCentre (0 0 0);
    }

    walls
    {
        type wall;
        faces
        (
            (4 11 12 5)
           (5 12 13 6)
        );
    }

   
);

I have calculated the D and F coefficients,the porous media is homogeneous and isotropic , the porousZones is like below:

1
(
    porosity
    {
        coordinateSystem
        {
            e1  (1 0 0);
            e2  (0 1 0);
        }

        Darcy
        {
            d   d [0 -2 0 0 0 0 0] (9.208e8 9.208e8 9.2208e8);
            f   f [0 -1 0 0 0 0 0] (1.375e4 1.375e4 1.375e4);
        }
    }
)

and U file is like below:

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform (0.09934 0 0);                                 
    }
    outlet
    {
        type            zeroGradient;
    }
    side1
    {
        type            symmetryPlane;
       
    }
    side2
    {
        type           symmetryPlane;
       
    }
    walls
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
 
}

my questions:

1.is every thing correct ? specially the Red lines
2.should I use slip boundary condition for walls in U file or not?


I've attached pipe picture.

Thanks in advance
Arjang

[deepbandivadekar]

Sorry, bumping up a thread too old.

But I'd like to know this too. Did you work this out? Is it correct?
Anyone else more experienced has any feedback?

[clapointe]

An alternative is creating cellSets with topoSet and converting them to zones. Then you can use fvOptions to apply the momentum sink as mentioned previously.

Caelan

[deepbandivadekar]

Originally Posted by clapointe

An alternative is creating cellSets with topoSet and converting them to zones. Then you can use fvOptions to apply the momentum sink as mentioned previously.

Caelan

So that'd mean what he has posted is correct?

[clapointe]

Yes. I'm sure about defining the porous regions within blockMeshDict but think it can be done -- I'd mentioned topoSet to create the zones as a workaround. I'm sure I can dig up an old example case if that'd help.

Caelan

[deepbandivadekar]

Originally Posted by clapointe

Yes. I'm sure about defining the porous regions within blockMeshDict but think it can be done -- I'd mentioned topoSet to create the zones as a workaround. I'm sure I can dig up an old example case if that'd help.

Caelan

Yes, please. That'd help.

[clapointe]

Here's an old case I made a while ago -- for OpenFoam 2.2.x I think. Hope it helps.

Caelan