[ajucgnr]

While creating a mesh for an axisymmetric geometry of a cylinder I get the following error while mentioning the face for the axis.


ERROR:

Trying to specify a boundary face 4(0 4 4 0) 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 2 named axis.

From function void Foam::polyMesh::setTopology(const cellShapeList&, const faceListList&, const wordList&, Foam::labelList&, Foam::labelList&, Foam::label&, Foam::label&, Foam::cellList&)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 324.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::polyMesh::setTopology(Foam::List<Foam::cellS hape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::List<int>&, Foam::List<int>&, int&, int&, Foam::List<Foam::cell>&) at ??:?
#3 Foam::polyMesh::polyMesh(Foam::IOobject const&, Foam::Xfer<Foam::Field<Foam::Vector<double> > > const&, Foam::List<Foam::cellShape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::PtrList<Foam::dictionary> const&, Foam::word const&, Foam::word const&, bool) at ??:?
#4 Foam::blockMesh::createTopology(Foam::IOdictionary const&, Foam::word const&) at ??:?
#5 Foam::blockMesh::blockMesh(Foam::IOdictionary const&, Foam::word const&) at ??:?
#6 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/blockMesh"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/blockMesh"
Aborted (core dumped)



The following is the blockMeshDict entry:




/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.001;

vertices
(
(0 0 0) //0
(-1 -1.732050807 0) //1
(-1 0 0) //2
(-1 1.732050807 0) //3
(0 0 -100) //4
(-1 -1.732050807 -100) //5
(-1 0 -100) //6
(-1 1.732050807 -100) //7

);

blocks
(
hex (0 1 2 0 4 5 6 4) (10 10 20) simpleGrading (1 1 1) //0
hex (0 2 3 0 4 6 7 4) (10 10 20) simpleGrading (1 1 1) //1


);

edges
(

);

boundary
(
front
{
type cyclic;
neighbourPatch back;
faces
( (0 4 7 3)

);
}
back
{
type cyclic;
neighbourPatch front;
faces
(
(1 5 4 0)


);
}

axis
{
type empty;
faces
(
(0 4 4 0)

);
}
);



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


I am trying to create two blocks for the entire wedge section of the cylinder (for making the geometry easy I have replaced the circular arc edge with just a straight line). It seems the error is because the face (0 4 4 0) is already used for making the Internal face connecting the two blocks. Thus a repetition of the face is not allowed. This is my guess. What could possibly be done to define the axis (0 4 4 0)?