CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

multiple porous blocks and correct setting in blockMesh and porous Zone Dictionaries

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   November 27, 2013, 06:53
Default multiple porous blocks and correct setting in blockMesh and porous Zone Dictionaries
  #1
Member
 
Arjang Behnoud
Join Date: Oct 2012
Posts: 61
Rep Power: 7
AmirBaqa1987 is on a distinguished road
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:

Code:
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:
Code:
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:
Code:
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
Attached Images
File Type: jpg pipe.jpg (55.5 KB, 20 views)
AmirBaqa1987 is offline   Reply With Quote

Old   January 17, 2018, 11:09
Default
  #2
Member
 
Deep
Join Date: Oct 2017
Posts: 44
Rep Power: 2
deepbandivadekar is on a distinguished road
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?
deepbandivadekar is offline   Reply With Quote

Old   January 19, 2018, 14:12
Default
  #3
Member
 
Join Date: Aug 2015
Posts: 31
Rep Power: 4
clapointe is on a distinguished road
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
clapointe is offline   Reply With Quote

Old   Yesterday, 11:06
Default
  #4
Member
 
Deep
Join Date: Oct 2017
Posts: 44
Rep Power: 2
deepbandivadekar is on a distinguished road
Quote:
Originally Posted by clapointe View Post
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?
deepbandivadekar is offline   Reply With Quote

Old   Yesterday, 14:05
Default
  #5
Member
 
Join Date: Aug 2015
Posts: 31
Rep Power: 4
clapointe is on a distinguished road
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
clapointe is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 04:10.