[tanveerfathima]

Hai. I am new to OpenFOAM and facing difficulty in creating a blockMesh. I have written down a blockMeshDict file for a tapered cylinder and was trying to mesh it in cavity case just to see how it looks like in paraview.
This is the code:

vertices
(
      (1 0 0) //0
      (0.707 0.707 0) //1
      (0 1 0) //2
      (-0.707 0.707 0) //3
      (-1 0 0) //4
      (-0.707 -0.707 0) //5
      (0 -1 0) //6
      (0.707 -0.707 0) //7
      (2 0 5) //8
      (1.414  1.414 5) //9
      (0 2 0) //10
      (-1.414 1.414 5) //11
      (-2 0 5) //12
      (-1.414 -1.414 5) //13
      (0 -2 5) //14
      (1.414 -1.414 5) //15
);

blocks
(
      hex (0 1 2 3 11 10 9 8) (10 10 10) simpleGrading(1 1 1)
      hex (4 5 6 7 15 14 13 12) (10 10 10) simpleGrading(1 1 1)
);

edges
(
      arc 0 2 (0.707 0.707 0)
      arc 2 4 (-0.707 0.707 0)
      arc 4 6 (-0.707 -0.707 0)
      arc 6 0 (0.707 -0.707 0)
      arc 8 10 (1.414 1.414 5)
      arc 10 12 (-1.414 1.414 5)
      arc 12 14 (-1.414 -1.414 5)
      arc 14 8 (1.414 -1.414 5)
);

boundary
(
 fixedWalls
 {
    type wall;
    faces
    (
      (0 1 9 8)
      (1 2 10 9)
      (2 3 11 10)
      (3 4 12 11)
      (4 5 13 12)
      (5 6 14 13)
      (6 7 15 14)
      (7 0 8 15)   
);
}
  frontAndBack
  {
    type wall;
    faces
    (
      (0 1 2 3)
      (3 4 5 6)
      (6 7 0 1)
      (8 9 10 11)
      (11 12 13 14)
      (14 15 8 9)
    );
   }
);

mergePatchPairs
(
);

I am getting this error while meshing it:

Create time

Creating block mesh from
    "/home/tanu/OpenFOAM/OpenFOAM-v1806/cylinder/cavity/system/blockMeshDict"
Creating block edges
No non-planar block faces defined
Creating topology blocks
Creating topology patches

Creating block mesh topology


--> FOAM FATAL ERROR:
face 3 in patch 0 does not have neighbour cell face: 4(3 4 12 11)

    From function Foam::labelList Foam::polyMesh::facePatchFaceCells(const faceList&, const labelListList&, const faceListList&, Foam::label) const
    in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 116.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&)sh: 1: addr2line: not found
 addr2line failed
#1  Foam::error::abort()sh: 1: addr2line: not found
 addr2line failed
#2  Foam::polyMesh::facePatchFaceCells(Foam::List<Foam::face> const&, Foam::List<Foam::List<int> > const&, Foam::List<Foam::List<Foam::face> > const&, int) constsh: 1: addr2line: not found
 addr2line failed
#3  Foam::polyMesh::setTopology(Foam::List<Foam::cellShape> 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>&)sh: 1: addr2line: not found
 addr2line failed
#4  Foam::polyMesh::polyMesh(Foam::IOobject const&, Foam::Field<Foam::Vector<double> >&&, 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)sh: 1: addr2line: not found
 addr2line failed
#5  Foam::blockMesh::createTopology(Foam::IOdictionary const&, Foam::word const&)sh: 1: addr2line: not found
 addr2line failed
#6  Foam::blockMesh::blockMesh(Foam::IOdictionary const&, Foam::word const&)sh: 1: addr2line: not found
 addr2line failed
#7  ?sh: 1: addr2line: not found
 addr2line failed
#8  __libc_start_mainsh: 1: addr2line: not found
 addr2line failed
#9  ?sh: 1: addr2line: not found
 addr2line failed
Aborted (core dumped)

What changes should I make in the file?

[yambanshee]

Please make use of the code tags when attaching code.


Also, with a question like this: it is very helpful to attach an image showing how the verticies/blocks look. I have attached something of the sorts here (all be it a terrible quaility one: I'm running out of paper).



In doing this I have found serveral issues:



Firstly, the two blocks you have defined do not connect to each other. This is your largest issue. This can be 'fixed' by making another block that goes along the lines of (4 7 0 3 12 15 8 11)


Secondly, your blocks will all be twisted due to incorrect ordering of verticies. your first block should be along the lines of (0 1 2 3 8 9 10 11) and second (4 5 6 7 12 13 14 15)


Thirdly, your definition of edges aren't right. As an example, you define an arc between points 0 and 2, but there is no direct connection between them. To accomplish this, you will need to define an arc between 0 and 1, and then 1 and 2.


My last point would be a concern over mesh quality given your block construction. I suspect that the cells by the verticies 1, 2, 5, 6, 9, 10, 13 and 14 will be quite low quality. The standard block structure for a pipe is attached in the second figure