CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (http://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   internal boundary (http://www.cfd-online.com/Forums/openfoam-meshing/130104-internal-boundary.html)

mary mor February 19, 2014 15:32

internal boundary
 
1 Attachment(s)
Hi dear all
I want to make block of meshes like the image attached in the post.

The nine partitions below are introduced in blockMeshDict. I'd like to introduce the internal boundary inlet as a patch. The blockMeshDict file is this:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  1.4                                  |
|  \\  /    A nd          | Web:      http://www.openfoam.org              |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/

FoamFile
{
    version        2.0;
    format          ascii;
    class          dictionary;
    object          blockMeshDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


convertToMeters 1;

vertices
(
    (0 -0.7 0)//0
    (1 -0.7 0)//1
    (2 -0.7 0)//2
    (4 -0.7 0)//3

    (0 -0.04 0)//4
    (1 -0.04 0)//5
    (2 -0.04 0)//6
    (4 -0.04 0)//7

    (0 0.04 0)//8
    (1 0.04 0)//9
    (2 0.04 0)//10
    (4 0.04 0)//11

    (0 0.08  0)//12
    (1 0.08  0)//13
    (2 0.08  0)//14
    (4 0.08  0)//15

    (0 -0.7 0.01)//16
    (1 -0.7 0.01)//17
    (2 -0.7 0.01)//18
    (4 -0.7 0.01)//19

    (0 -0.04 0.01)//20
    (1 -0.04 0.01)//21
    (2 -0.04 0.01)//22
    (4 -0.04 0.01)//23

    (0 0.04 0.01)//24
    (1 0.04 0.01)//25
    (2 0.04 0.01)//26
    (4 0.04 0.01)//27

    (0 0.08  0.01)//28
    (1 0.08  0.01)//29
    (2 0.08  0.01)//30
    (4 0.08  0.01)//31
);


blocks
(
    hex (0 1 5 4 16 17 21 20) (60 40 1) simpleGrading (1 0.25 1)
    hex (1 2 6 5 17 18 22 21) (60 40 1) simpleGrading (1 0.25 1)
    hex (2 3 7 6 18 19 23 22) (125 40 1) simpleGrading (1 0.25 1)

    hex (4 5 9 8 20 21 25 24) (60 10 1) simpleGrading (1 1 1)
    hex (5 6 10 9 21 22 26 25) (60 10 1) simpleGrading (1 1 1)
    hex (6 7 11 10 22 23 27 26) (125 10 1) simpleGrading (1 1 1)

    hex (8 9 13 12 24 25 29 28) (60 5 1) simpleGrading (1 3 1)
    hex (9 10 14 13 25 26 30 29) (60 5 1) simpleGrading (1 3 1)
    hex (10 11 15 14 26 27 31 30) (125 5 1) simpleGrading (1 3 1)
);

edges
(
);

patches
(
    patch left
    (
        (0 4 20 16)
        (4 8 24 20)
        (8 12 28 24)
    )

    patche inlet
    (
        (1 5 21 17)
        (5 9 25 21)
        (9 13 29 25)
    )

    patch right
    (
        (3 7 23 19)
        (7 11 27 23)
        (11 15 31 27)
    )

    patch bottom
    (
        (0 1 17 16)
        (1 2 18 17)
        (2 3 19 18)
    )

    patch top
    (
        (12 13 29 28)
        (13 14 30 29)
        (14 15 31 30)
    )

    empty frontAndBack
    (
        // Back
        (0 1 5 4)
        (1 2 6 5)
        (2 3 7 6)
        (4 5 9 8)
        (5 6 10 9)
        (6 7 11 10)
        (8 9 13 12)
        (9 10 14 13)
        (10 11 15 14)

        // Front
        (16 17 21 20)
        (17 18 22 21)
        (18 19 23 22)
        (20 21 25 24)
        (21 22 26 25)
        (22 23 27 26)
        (24 25 29 28)
        (25 26 30 29)
        (26 27 31 30)
    )
);

mergePatchPairs
();

// ************************************************************************* //


But after running, it gives me the error below:
Code:

--> FOAM FATAL ERROR:
Trying to specify a boundary face 4(1 5 21 17) on the face on cell 0 which is either an internal face or already belongs to some other patch.  This is face 0 of patch 1 named inlet.

    From function polyMesh::polyMesh
(
    const IOobject&,
    const Xfer<pointField>&,
    const cellShapeList& cellsAsShapes,
    const faceListList& boundaryFaces,
    const wordList& boundaryPatchTypes,
    const wordList& boundaryPatchNames,
    const word& defaultBoundaryPatchType
)
    in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 489.

FOAM aborting

Aborted

Is it possible to define an internal boundary as a patch? Or I'm doing anything wrong?

Thanks
Best regards


All times are GMT -4. The time now is 00:19.