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

[Other] internal boundary

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 19, 2014, 15:32
Default internal boundary
  #1
Member
 
maryam morta
Join Date: Sep 2013
Posts: 54
Rep Power: 12
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, 16 views)
mary mor is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Foam::error::PrintStack almir OpenFOAM Running, Solving & CFD 91 December 21, 2022 05:50
Radiation in semi-transparent media with surface-to-surface model? mpeppels CFX 11 August 22, 2019 08:30
Question about adaptive timestepping Guille1811 CFX 25 November 12, 2017 18:38
Error - Solar absorber - Solar Thermal Radiation MichaelK CFX 12 September 1, 2016 06:15
CFX13 Post Periodic interface EtaEta CFX 7 December 8, 2011 18:15


All times are GMT -4. The time now is 02:59.