CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Native Meshers: blockMesh

problem with mesh generation in blockMesh

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 22, 2013, 13:26
Default problem with mesh generation in blockMesh
  #1
New Member
 
Luis Miguel
Join Date: Apr 2013
Location: Colombia
Posts: 13
Rep Power: 5
luigi21 is on a distinguished road
I'm trying to compile the mesh below in waveFoam, but I always get the same error:

FOAM FATAL ERROR:
Inconsistent number of faces between block pair 1 and 2

From function blockMesh::calcMergeInfo()
in file blockMesh/blockMeshMerge.C at line 221.

if somebody knows which the error is? i'll be very gratefull....


vertices
(
(0 0 0)
(8 0 0)
(12 0 0)
(20 0 0)
(8.35 0.35 0)
(11.65 0.35 0)
(0 1 0)
(8 1 0)
(8.35 1 0)
(11.65 1 0)
(12 1 0)
(20 1 0)
(0 0 1)
(8 0 1)
(12 0 1)
(20 0 1)
(8.35 0.35 1)
(11.65 0.35 1)
(0 1 1)
(8 1 1)
(8.35 1 1)
(11.65 1 1)
(12 1 1)
(20 1 1)
);

blocks
(
hex (0 1 7 6 12 13 19 18) (40 5 5) simpleGrading (1 1 1)
hex (1 4 8 7 13 16 20 19) (2 5 5) simpleGrading (1 1 1)
hex (4 5 9 8 16 17 21 20) (16 3 5) simpleGrading (1 1 1)
hex (5 2 10 9 17 14 22 21) (2 5 5) simpleGrading (1 1 1)
hex (2 3 11 10 14 15 23 22) (40 5 5) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
leftWall
{
type wall;
faces
(
(0 12 18 6)
);
}
rightWall
{
type wall;
faces
(
(3 15 23 11)
);
}
lowerWall
{
type wall;
faces
(
(0 1 13 12)
(1 4 16 13)
(4 5 17 16)
(2 14 17 5)
(2 3 15 14)
);
}
atmosphere
{
type patch;
faces
(
(6 18 19 7)
(7 19 20 8)
(8 20 21 9)
(9 21 22 10)
(10 22 23 11)
);
}
);

mergePatchPairs
(
);
luigi21 is offline   Reply With Quote

Old   April 22, 2013, 14:12
Default
  #2
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,702
Rep Power: 27
ngj will become famous soon enoughngj will become famous soon enough
Hi Luis,

Your third hex-definition has the wrong number of points in the vertical direction. It should read like this:

Code:
hex (4 5 9 8 16 17 21 20) (16 5 5) simpleGrading (1 1 1)
Kind regards

Niels
ngj is offline   Reply With Quote

Old   April 23, 2013, 14:34
Default
  #3
New Member
 
Luis Miguel
Join Date: Apr 2013
Location: Colombia
Posts: 13
Rep Power: 5
luigi21 is on a distinguished road
thanks a lot for the help Niels....the mesh work well but it was rotated, the waves are propagating in the bottom of the mesh. how can I rotate the mesh or waves field for the propagation of waves are in the top of the mesh and not in the bottom?....
luigi21 is offline   Reply With Quote

Old   April 23, 2013, 14:45
Default
  #4
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,702
Rep Power: 27
ngj will become famous soon enoughngj will become famous soon enough
Hi Luis,

Could you please post some graphics, because I am unsure of what you mean. Though, the easy answer is for you to re-define your blockMeshDict, so it reflects the geometry, which you want to model.

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   April 24, 2013, 11:20
Default
  #5
New Member
 
Luis Miguel
Join Date: Apr 2013
Location: Colombia
Posts: 13
Rep Power: 5
luigi21 is on a distinguished road
Hi niels,

The problem is a wave channel with a submerged structure, In the attached images you can watch that the waves are propagating in the wrong face of the mesh (y,z) when the right direction is (x,y)....so how can i rotate either mesh or wave direction of propagation?

link images
http://imgur.com/FDHYCVf
http://imgur.com/TqljqoN
http://imgur.com/TMsWaxT
luigi21 is offline   Reply With Quote

Old   April 24, 2013, 12:57
Default
  #6
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,702
Rep Power: 27
ngj will become famous soon enoughngj will become famous soon enough
OK. All you will need is to set the direction of the gravity vector in the vertical direction. Right now you probably have it parallel with your horizontal plane.

Furthermore, I would suggest that you use

Code:
setWaveField
before executing

Code:
waveFoam
It looks as if your domain is empty, and that will create problems for you. This is resolved using setWaveField.

Kind regards

Niels
ngj is offline   Reply With Quote

Old   April 24, 2013, 13:54
Default
  #7
New Member
 
Luis Miguel
Join Date: Apr 2013
Location: Colombia
Posts: 13
Rep Power: 5
luigi21 is on a distinguished road
hi niels,

Thanks a lot for all the help...all the simulation worked fine, but now I have a little doubt, why the water looks like if this were flying, the bottom of the channel have 0 flow and above of this points the flow grow up? look at the pic after the structure...

http://imgur.com/5oBzRa4
luigi21 is offline   Reply With Quote

Old   April 24, 2013, 14:01
Default
  #8
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,702
Rep Power: 27
ngj will become famous soon enoughngj will become famous soon enough
Too poor resolution?

- Niels
ngj is offline   Reply With Quote

Old   April 24, 2013, 15:51
Default
  #9
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,702
Rep Power: 27
ngj will become famous soon enoughngj will become famous soon enough
What I mean is that 5 computational cells over the entire width/height of the computational domain is by far too little.

Kind regards

Niels
ngj is offline   Reply With Quote

Old   April 30, 2013, 15:52
Default
  #10
New Member
 
Luis Miguel
Join Date: Apr 2013
Location: Colombia
Posts: 13
Rep Power: 5
luigi21 is on a distinguished road
I've had some problems with boundary conditions, is the wavefoam's article free (A Wave Generation Toolbox for the Open-Source CFD Library)?....where can I get it?
luigi21 is offline   Reply With Quote

Old   June 3, 2013, 12:58
Default problem with blockMesh
  #11
New Member
 
Bhupesh
Join Date: Jun 2012
Location: Germany
Posts: 14
Rep Power: 6
bhups45 is on a distinguished road
I am unable to find out the cause of the following error while executing blockMesh. The error is:

blockMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.0-5be49240882f
Exec : blockMesh
Date : Jun 03 2013
Time : 21:06:31
Host : "ubuntu"
PID : 2459
Case : /home/bhupesh/OpenFOAM/bhupesh-2.2.0/run/tutorials/incompressible/icoFoam/pistonCylinder
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Creating block mesh from
"/home/bhupesh/OpenFOAM/bhupesh-2.2.0/run/tutorials/incompressible/icoFoam/pistonCylinder/constant/polyMesh/blockMeshDict"
Creating curved edges
Creating topology blocks
Creating topology patches

Creating block mesh topology
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 128
zero or negative pyramid volume: -9.5e-05 for face 2
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 128
zero or negative pyramid volume: -9.5e-05 for face 2
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 128
zero or negative pyramid volume: -9.5e-05 for face 3
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 128
zero or negative pyramid volume: -9.5e-05 for face 3


--> FOAM FATAL ERROR:
face 10 in patch 4 does not have neighbour cell face: 4(30 31 35 34)

From function polyMesh::facePatchFaceCells(const faceList& patchFaces,const labelListList& pointCells,const faceListList& cellsFaceShapes,const label patchID)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 127.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam:olyMesh::facePatchFaceCells(Foam::List<Foam ::face> const&, Foam::List<Foam::List<int> > const&, Foam::List<Foam::List<Foam::face> > const&, int) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#3 Foam:olyMesh::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>&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam:olyMesh:olyMesh(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) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::blockMesh::createTopology(Foam::IOdictionary const&, Foam::word const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libblockMesh.so"
#6 Foam::blockMesh::blockMesh(Foam::IOdictionary const&, Foam::word const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libblockMesh.so"
#7
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/blockMesh"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/blockMesh"
Aborted (core dumped)

The blockMeshDict file is :

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

convertToMeters 1;

vertices
(
(0.5 0.25 -0.25)//0
(0.006 0.25 -0.25)//1
(-0.006 0.25 -0.25)//2
(-1.5 0.25 -0.25)//3
(-1.5 -0.25 -0.25)//4
(-0.006 -0.25 -0.25)//5
(0.006 -0.25 -0.25)//6
(0.5 -0.25 -0.25)//7

(0.5 0.25 0.25)//8
(0.006 0.25 0.25)//9
(-0.006 0.25 0.25)//10
(-1.5 0.25 0.25)//11
(-1.5 -0.25 0.25)//12
(-0.006 -0.25 0.25)//13
(0.006 -0.25 0.25)//14
(0.5 -0.25 0.25)//15

(0.006 0.25 0)//16
(0.006 0.06 -0.25)//17
(0.006 0.06 0)//18
(0.006 0.06 0.25)//19
(0.006 0 -0.25)//20
(0.006 0 -0.06)//21
(0.006 0 0.06)//22
(0.006 0 0.25)//23
(0.006 -0.06 -0.25)//24
(0.006 -0.06 0)//25
(0.006 -0.06 0.25)//26
(0.006 -0.25 0)//27

(-0.006 0.25 0)//28
(-0.006 0.06 -0.25)//29
(-0.006 0.06 0)//30
(-0.006 0.06 0.25)//31
(-0.006 0 -0.25)//32
(-0.006 0 -0.06)//33
(-0.006 0 0.06)//34
(-0.006 0 0.25)//35
(-0.006 -0.06 -0.25)//36
(-0.006 -0.06 0)//37
(-0.006 -0.06 0.25)//38
(-0.006 -0.25 0)//39

(-0.354 0 -0.06)//40
(-0.354 0.06 0)//41
(-0.354 0 0.06)//42
(-0.354 -0.06 0)//43

);

blocks
(
hex (0 1 6 7 8 9 14 15) (20 20 20) simpleGrading (1 1 1)
hex (1 2 5 6 9 10 13 14) (20 20 20) simpleGrading (1 1 1)
hex (2 3 4 5 10 11 12 13) (20 20 20) simpleGrading (1 1 1)
hex (1 16 18 19 2 28 30 31) (20 20 20) simpleGrading (1 1 1)
hex (16 9 17 18 28 10 29 30) (20 20 20) simpleGrading (1 1 1)
hex (18 17 20 21 30 29 32 33) (20 20 20) simpleGrading (1 1 1)
hex (21 20 24 25 33 32 36 37) (20 20 20) simpleGrading (1 1 1)
hex (25 24 14 27 37 36 13 39) (20 20 20) simpleGrading (1 1 1)
hex (26 25 27 6 38 37 39 5) (20 20 20) simpleGrading (1 1 1)
hex (23 22 25 26 35 34 37 38) (20 20 20) simpleGrading (1 1 1)
hex (19 18 22 23 31 30 34 33) (20 20 20) simpleGrading (1 1 1)
hex (22 18 21 25 42 41 40 43) (20 20 20) simpleGrading (1 1 1)
);

edges
(
arc 25 21 (0.006 -0.00424 -0.00424)
arc 21 18 (0.006 0.00424 -0.00424)
arc 18 22 (0.006 0.00424 0.00424)
arc 22 25 (0.006 -0.00424 0.00424)
arc 37 33 (-0.006 -0.00424 -0.00424)
arc 33 30 (-0.006 0.00424 -0.00424)
arc 30 34 (-0.006 0.00424 0.00424)
arc 34 37 (-0.006 -0.00424 0.00424)

);

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

middleWall
{
type wall;
faces
(
(1 16 18 19)
(16 9 17 18)
(18 17 20 21)
(21 20 24 25)
(25 24 14 27)
(26 25 27 6)
(23 22 25 26)
(19 18 22 23)

(10 28 30 29)
(28 2 31 30)
(30 31 35 34)
(34 35 38 37)
(37 38 5 39)
(36 37 39 13)
(32 33 37 36)
(29 30 33 32)
);
}

cylinderInlet
{
type patch;
faces
(
(22 18 21 25)
);
}

cylinderWall
{
type wall;
faces
(
(21 18 30 33)
(18 22 34 30)
(22 25 37 34)
(25 21 33 37)
);
}
);

mergePatchPairs
(
);

// ************************************************** *********************** //
I am trying to make a geometry of a cylinder with a wall in it and there is a circular hole in the middle wall through which fluid passes.
bhups45 is offline   Reply With Quote

Old   June 3, 2013, 13:01
Default problem with blockMesh
  #12
New Member
 
Bhupesh
Join Date: Jun 2012
Location: Germany
Posts: 14
Rep Power: 6
bhups45 is on a distinguished road
I am unable to find out the cause of the following error while executing blockMesh. Someone please help me out with this error. The error is:

blockMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.0-5be49240882f
Exec : blockMesh
Date : Jun 03 2013
Time : 21:06:31
Host : "ubuntu"
PID : 2459
Case : /home/bhupesh/OpenFOAM/bhupesh-2.2.0/run/tutorials/incompressible/icoFoam/pistonCylinder
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Creating block mesh from
"/home/bhupesh/OpenFOAM/bhupesh-2.2.0/run/tutorials/incompressible/icoFoam/pistonCylinder/constant/polyMesh/blockMeshDict"
Creating curved edges
Creating topology blocks
Creating topology patches

Creating block mesh topology
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 128
zero or negative pyramid volume: -9.5e-05 for face 2
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 128
zero or negative pyramid volume: -9.5e-05 for face 2
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 128
zero or negative pyramid volume: -9.5e-05 for face 3
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 128
zero or negative pyramid volume: -9.5e-05 for face 3


--> FOAM FATAL ERROR:
face 10 in patch 4 does not have neighbour cell face: 4(30 31 35 34)

From function polyMesh::facePatchFaceCells(const faceList& patchFaces,const labelListList& pointCells,const faceListList& cellsFaceShapes,const label patchID)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 127.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam:olyMesh::facePatchFaceCells(Foam::List<Foam ::face> const&, Foam::List<Foam::List<int> > const&, Foam::List<Foam::List<Foam::face> > const&, int) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#3 Foam:olyMesh::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>&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam:olyMesh:olyMesh(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) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::blockMesh::createTopology(Foam::IOdictionary const&, Foam::word const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libblockMesh.so"
#6 Foam::blockMesh::blockMesh(Foam::IOdictionary const&, Foam::word const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libblockMesh.so"
#7
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/blockMesh"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/blockMesh"
Aborted (core dumped)

I am trying to make a geometry of a cube with a wall inside it and there is a circular hole in the middle wall through which fluid passes. The blockMeshDict file is :

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

convertToMeters 1;

vertices
(
(0.5 0.25 -0.25)//0
(0.006 0.25 -0.25)//1
(-0.006 0.25 -0.25)//2
(-1.5 0.25 -0.25)//3
(-1.5 -0.25 -0.25)//4
(-0.006 -0.25 -0.25)//5
(0.006 -0.25 -0.25)//6
(0.5 -0.25 -0.25)//7

(0.5 0.25 0.25)//8
(0.006 0.25 0.25)//9
(-0.006 0.25 0.25)//10
(-1.5 0.25 0.25)//11
(-1.5 -0.25 0.25)//12
(-0.006 -0.25 0.25)//13
(0.006 -0.25 0.25)//14
(0.5 -0.25 0.25)//15

(0.006 0.25 0)//16
(0.006 0.06 -0.25)//17
(0.006 0.06 0)//18
(0.006 0.06 0.25)//19
(0.006 0 -0.25)//20
(0.006 0 -0.06)//21
(0.006 0 0.06)//22
(0.006 0 0.25)//23
(0.006 -0.06 -0.25)//24
(0.006 -0.06 0)//25
(0.006 -0.06 0.25)//26
(0.006 -0.25 0)//27

(-0.006 0.25 0)//28
(-0.006 0.06 -0.25)//29
(-0.006 0.06 0)//30
(-0.006 0.06 0.25)//31
(-0.006 0 -0.25)//32
(-0.006 0 -0.06)//33
(-0.006 0 0.06)//34
(-0.006 0 0.25)//35
(-0.006 -0.06 -0.25)//36
(-0.006 -0.06 0)//37
(-0.006 -0.06 0.25)//38
(-0.006 -0.25 0)//39

(-0.354 0 -0.06)//40
(-0.354 0.06 0)//41
(-0.354 0 0.06)//42
(-0.354 -0.06 0)//43

);

blocks
(
hex (0 1 6 7 8 9 14 15) (20 20 20) simpleGrading (1 1 1)
hex (1 2 5 6 9 10 13 14) (20 20 20) simpleGrading (1 1 1)
hex (2 3 4 5 10 11 12 13) (20 20 20) simpleGrading (1 1 1)
hex (1 16 18 19 2 28 30 31) (20 20 20) simpleGrading (1 1 1)
hex (16 9 17 18 28 10 29 30) (20 20 20) simpleGrading (1 1 1)
hex (18 17 20 21 30 29 32 33) (20 20 20) simpleGrading (1 1 1)
hex (21 20 24 25 33 32 36 37) (20 20 20) simpleGrading (1 1 1)
hex (25 24 14 27 37 36 13 39) (20 20 20) simpleGrading (1 1 1)
hex (26 25 27 6 38 37 39 5) (20 20 20) simpleGrading (1 1 1)
hex (23 22 25 26 35 34 37 38) (20 20 20) simpleGrading (1 1 1)
hex (19 18 22 23 31 30 34 33) (20 20 20) simpleGrading (1 1 1)
hex (22 18 21 25 42 41 40 43) (20 20 20) simpleGrading (1 1 1)
);

edges
(
arc 25 21 (0.006 -0.00424 -0.00424)
arc 21 18 (0.006 0.00424 -0.00424)
arc 18 22 (0.006 0.00424 0.00424)
arc 22 25 (0.006 -0.00424 0.00424)
arc 37 33 (-0.006 -0.00424 -0.00424)
arc 33 30 (-0.006 0.00424 -0.00424)
arc 30 34 (-0.006 0.00424 0.00424)
arc 34 37 (-0.006 -0.00424 0.00424)

);

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

middleWall
{
type wall;
faces
(
(1 16 18 19)
(16 9 17 18)
(18 17 20 21)
(21 20 24 25)
(25 24 14 27)
(26 25 27 6)
(23 22 25 26)
(19 18 22 23)

(10 28 30 29)
(28 2 31 30)
(30 31 35 34)
(34 35 38 37)
(37 38 5 39)
(36 37 39 13)
(32 33 37 36)
(29 30 33 32)
);
}

cylinderInlet
{
type patch;
faces
(
(22 18 21 25)
);
}

cylinderWall
{
type wall;
faces
(
(21 18 30 33)
(18 22 34 30)
(22 25 37 34)
(25 21 33 37)
);
}
);

mergePatchPairs
(
);

// ************************************************** *********************** //
bhups45 is offline   Reply With Quote

Old   June 4, 2013, 02:29
Default
  #13
Member
 
Yosmcer Mocktai
Join Date: Apr 2013
Location: Behind a computer
Posts: 48
Rep Power: 9
Yosmcer will become famous soon enough
Hello, the error log gives you one warning and one error:

  1. zero or negative pyramid volume
  2. face 10 in patch 4 does not have neighbour cell face: 4(30 31 35 34)
I think there is several problems in your mesh.
First, I see that your circular junction is not really circular, arcs are rounded inside instead of outside (see picture attached).
Secondly, block 1 (remeber we start counting with 0) and block 11 intersect.
Thirdly, check all the blocks definitions from 3 to 10 to be sure all perform as you wish.

To quickly check the blocs configuration, remove* (don't forget to save your blockMeshDict file before) all the boundary patches definitions, and just run:
Code:
paraFoam -block
Without the need to make any blockMesh before.

I use the "surface with edges" in paraVieuw to see the blocks and the edges. The blocks are represented with no curves, that's why it is better to also display the edges to see how the curves (arcs) are defined.

Good corrections

* I advice to remove them because "paraFoam -block" will not run as it consider that there is the face error. Once you corrected the blocks, you can add the faces and check with "paraFoam -block" if it still work, and correct errors if some.
Attached Images
File Type: jpg blocks1.jpg (15.3 KB, 28 views)
Yosmcer is offline   Reply With Quote

Old   June 4, 2013, 11:07
Default problem with blockMesh
  #14
New Member
 
Bhupesh
Join Date: Jun 2012
Location: Germany
Posts: 14
Rep Power: 6
bhups45 is on a distinguished road
Thanks Yosmcer for your reply.

By seeing the co-ordinate axes in your figure,I think that I have interchanged the y & z axes co-ordinates. Do you think this might be the problem?
Also I have to ask you, in case of this geometry, which face will be defined as upper face while defining blocks for the middle wall?
Thnx in advace.
bhups45 is offline   Reply With Quote

Old   June 4, 2013, 14:39
Default some new warnings
  #15
New Member
 
Bhupesh
Join Date: Jun 2012
Location: Germany
Posts: 14
Rep Power: 6
bhups45 is on a distinguished road
Hello Yosmcer,

I changed the orientation of axes and followed your advice of using 'paraFoam -block' first for viewing blocks.
I am getting the blocks as my requirements but I am again getting the same warning of negative volume and it is showing some blocks maybe defined inside out which I am unable to figure out why because I have read about the rules we have to follow while defining blocks in blockMesh from User Guide and I think I am following it correctly.
Can you please help me with this numbering of blocks?

My blockMesh file is:

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

convertToMeters 1;

vertices
(
(0.5 0.25 -0.25)//0
(0.006 0.25 -0.25)//1
(-0.006 0.25 -0.25)//2
(-1.5 0.25 -0.25)//3
(-1.5 -0.25 -0.25)//4
(-0.006 -0.25 -0.25)//5
(0.006 -0.25 -0.25)//6
(0.5 -0.25 -0.25)//7

(0.5 0.25 0.25)//8
(0.006 0.25 0.25)//9
(-0.006 0.25 0.25)//10
(-1.5 0.25 0.25)//11
(-1.5 -0.25 0.25)//12
(-0.006 -0.25 0.25)//13
(0.006 -0.25 0.25)//14
(0.5 -0.25 0.25)//15

(0.006 0 0.25)//16
(0.006 -0.25 0.06)//17
(0.006 0 0.06)//18
(0.006 0.25 0.06)//19
(0.006 -0.25 0)//20
(0.006 -0.06 0)//21
(0.006 0.06 0)//22
(0.006 0 0.25)//23
(0.006 -0.25 -0.06)//24
(0.006 0 -0.06)//25
(0.006 0.25 -0.06)//26
(0.006 0 -0.25)//27

(-0.006 0 0.25)//28
(-0.006 -0.25 0.06)//29
(-0.006 0 0.06)//30
(-0.006 0.25 0.06)//31
(-0.006 -0.25 0)//32
(-0.006 -0.06 0)//33
(-0.006 0.06 0)//34
(-0.006 0.25 0)//35
(-0.006 -0.25 -0.06)//36
(-0.006 0 -0.06)//37
(-0.006 0.25 -0.06)//38
(-0.006 0 -0.25)//39

(-0.354 -0.06 0)//40
(-0.354 0 0.06)//41
(-0.354 0.06 0)//42
(-0.354 0 -0.06)//43

);

blocks
(
hex (8 9 14 15 0 1 6 7) (20 20 20) simpleGrading (1 1 1)
hex (9 10 13 14 1 2 5 6) (20 20 20) simpleGrading (1 1 1)
hex (10 11 12 13 2 3 4 5) (20 20 20) simpleGrading (1 1 1)
hex (1 16 18 19 2 28 30 31) (20 20 20) simpleGrading (1 1 1)
hex (16 9 17 18 28 10 29 30) (20 20 20) simpleGrading (1 1 1)
hex (18 17 20 21 30 29 32 33) (20 20 20) simpleGrading (1 1 1)
hex (21 20 24 25 33 32 36 37) (20 20 20) simpleGrading (1 1 1)
hex (25 24 14 27 37 36 13 39) (20 20 20) simpleGrading (1 1 1)
hex (26 25 27 6 38 37 39 5) (20 20 20) simpleGrading (1 1 1)
hex (23 22 25 26 35 34 37 38) (20 20 20) simpleGrading (1 1 1)
hex (19 18 22 23 31 30 34 33) (20 20 20) simpleGrading (1 1 1)
hex (22 18 21 25 42 41 40 43) (20 20 20) simpleGrading (1 1 1)
);

edges
(
arc 25 21 (0.006 -0.00424 -0.00424)
arc 21 18 (0.006 -0.00424 0.00424)
arc 18 22 (0.006 0.00424 0.00424)
arc 22 25 (0.006 0.00424 -0.00424)
arc 37 33 (-0.006 -0.00424 -0.00424)
arc 33 30 (-0.006 -0.00424 0.00424)
arc 30 34 (-0.006 0.00424 0.00424)
arc 34 37 (-0.006 0.00424 -0.00424)
arc 43 40 (0.006 -0.00424 -0.00424)
arc 40 41 (0.006 -0.00424 0.00424)
arc 41 42 (0.006 0.00424 0.00424)
arc 42 43 (0.006 0.00424 -0.00424)
);

boundary
(

);

mergePatchPairs
(
);

// ************************************************** *********************** //
bhups45 is offline   Reply With Quote

Old   June 6, 2013, 04:45
Default
  #16
Member
 
Yosmcer Mocktai
Join Date: Apr 2013
Location: Behind a computer
Posts: 48
Rep Power: 9
Yosmcer will become famous soon enough
Sorry for the delay.

With "paraFoam" block, you can select wich block you want to see (for example, only display block 2 and 11). that allow you to see each block independently and how they are jointed.

In the first attached picture, I selected surface with edges, and no block. So I see only the edges. It seems there is still a problem.
In the second one, I displayed only the block 11. You quickly see that this block is ill.
In the third picture, you can see that the 11'th block is define in the 2nd one. The volume of the 11th is a part of the volume of the 2nd, one same volume cant be share between two blocks.
Volume 1 contain volumes 2-10, that's the same problem.

Coordinate system is defined as (x, y, z).

Can you make a drawing of what you would like to have? I'm not sure as there is several places where blocs cover the sames volume.
Attached Images
File Type: jpg edges1.jpg (11.8 KB, 18 views)
File Type: jpg blocks2.jpg (22.3 KB, 23 views)
File Type: jpg blocks3.jpg (11.9 KB, 15 views)
Yosmcer 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
[ICEM] surface mesh merging problem everest ANSYS Meshing & Geometry 44 April 14, 2016 06:41
No layers in a small gap bobburnquist OpenFOAM Native Meshers: snappyHexMesh and Others 6 August 26, 2015 09:38
[ICEM] problem of Non-concidence mesh jeevankumarb ANSYS Meshing & Geometry 1 March 4, 2011 09:47
Problem with mesh generation Nestor FLUENT 3 November 1, 2006 03:54
Mesh generation problem Lam FLUENT 2 December 25, 2003 13:37


All times are GMT -4. The time now is 14:17.