CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > 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: 57
Rep Power: 5
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, 12 views)
AmirBaqa1987 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 05:53.