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

[blockMesh] Creating a tapered cylinder

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 2, 2018, 00:05
Default Creating a tapered cylinder
  #1
New Member
 
Tanveer
Join Date: Sep 2018
Posts: 6
Rep Power: 4
tanveerfathima is on a distinguished road
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:
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:

Code:
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?

Last edited by wyldckat; March 6, 2019 at 20:30. Reason: Added [CODE][/CODE] markers
tanveerfathima is offline   Reply With Quote

Old   December 4, 2018, 07:22
Default
  #2
Senior Member
 
Zander Meiring
Join Date: Jul 2018
Posts: 125
Rep Power: 4
yambanshee is on a distinguished road
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
Attached Images
File Type: jpeg 1.jpeg (20.9 KB, 13 views)
File Type: jpg 2.jpg (16.9 KB, 8 views)
yambanshee is offline   Reply With Quote

Reply

Tags
blockmesh, tapered_cylinder

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
[Commercial meshers] Fluent3DMeshToFoam simvun OpenFOAM Meshing & Mesh Conversion 50 January 19, 2020 16:33
Flow past a cylinder at Re 1e05 using LES, drag force coefficient is to low Scabbard Main CFD Forum 21 June 19, 2018 14:58
Forces Acting on a Rotating Cylinder (Moving Mesh) dreamchaser CFX 5 April 25, 2015 07:01
Problem in running ICEM grid in Openfoam Tarak OpenFOAM 6 September 9, 2011 18:51
Creating cylinder in STAR-CD Sachin Siemens 2 March 6, 2008 03:53


All times are GMT -4. The time now is 10:44.