CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (https://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   [blockMesh] problem with mesh generation in blockMesh (https://www.cfd-online.com/Forums/openfoam-meshing/116576-problem-mesh-generation-blockmesh.html)

luigi21 April 22, 2013 13:26

problem with mesh generation in blockMesh
 
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
(
);

ngj April 22, 2013 14:12

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

luigi21 April 23, 2013 14:34

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?....

ngj April 23, 2013 14:45

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

luigi21 April 24, 2013 11:20

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

ngj April 24, 2013 12:57

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

luigi21 April 24, 2013 13:54

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

ngj April 24, 2013 14:01

Too poor resolution?

- Niels

ngj April 24, 2013 15:51

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

luigi21 April 30, 2013 15:52

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?

bhups45 June 3, 2013 12:58

problem with blockMesh
 
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::printStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::polyMesh::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::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>&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 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) 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 June 3, 2013 13:01

problem with blockMesh
 
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::printStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::polyMesh::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::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>&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 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) 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
(
);

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

Yosmcer June 4, 2013 02:29

1 Attachment(s)
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.

bhups45 June 4, 2013 11:07

problem with blockMesh
 
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 June 4, 2013 14:39

some new warnings
 
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
(
);

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

Yosmcer June 6, 2013 04:45

3 Attachment(s)
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.


All times are GMT -4. The time now is 19:29.