CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Meshing & Mesh Conversion

internal boundary

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 19, 2014, 14:32
Default internal boundary
  #1
Member
 
maryam morta
Join Date: Sep 2013
Posts: 54
Rep Power: 2
mary mor is on a distinguished road
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
Attached Images
File Type: png inltet outlet.png (7.9 KB, 9 views)
mary mor 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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Foam::error::PrintStack almir OpenFOAM 45 July 8, 2014 04:12
Moving boundary destroy other internal meshing yhaomin2007 OpenFOAM 1 November 7, 2013 02:34
BlockMesh number internal faces and boundary ynos OpenFOAM Native Meshers: blockMesh 6 December 13, 2011 05:36
CFX13 Post Periodic interface EtaEta CFX 7 December 8, 2011 17:15
Internal Faces to boundary patch conversion robingilbert OpenFOAM 4 March 24, 2010 11:40


All times are GMT -4. The time now is 20:28.