|
[Sponsors] |
[blockMesh] problem with mesh generation in blockMesh |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 22, 2013, 13:26 |
problem with mesh generation in blockMesh
|
#1 |
New Member
Luis Miguel
Join Date: Apr 2013
Location: Colombia
Posts: 13
Rep Power: 13 |
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 ( ); |
|
April 22, 2013, 14:12 |
|
#2 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37 |
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) Niels |
|
April 23, 2013, 14:34 |
|
#3 |
New Member
Luis Miguel
Join Date: Apr 2013
Location: Colombia
Posts: 13
Rep Power: 13 |
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?....
|
|
April 23, 2013, 14:45 |
|
#4 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37 |
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 |
|
April 24, 2013, 11:20 |
|
#5 |
New Member
Luis Miguel
Join Date: Apr 2013
Location: Colombia
Posts: 13
Rep Power: 13 |
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 |
|
April 24, 2013, 12:57 |
|
#6 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37 |
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 Code:
waveFoam Kind regards Niels |
|
April 24, 2013, 13:54 |
|
#7 |
New Member
Luis Miguel
Join Date: Apr 2013
Location: Colombia
Posts: 13
Rep Power: 13 |
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 |
|
April 24, 2013, 14:01 |
|
#8 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37 |
Too poor resolution?
- Niels |
|
April 24, 2013, 15:51 |
|
#9 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37 |
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 |
|
April 30, 2013, 15:52 |
|
#10 |
New Member
Luis Miguel
Join Date: Apr 2013
Location: Colombia
Posts: 13
Rep Power: 13 |
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?
|
|
June 3, 2013, 12:58 |
problem with blockMesh
|
#11 |
New Member
Bhupesh
Join Date: Jun 2012
Location: Germany
Posts: 14
Rep Power: 13 |
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. |
|
June 3, 2013, 13:01 |
problem with blockMesh
|
#12 |
New Member
Bhupesh
Join Date: Jun 2012
Location: Germany
Posts: 14
Rep Power: 13 |
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 ( ); // ************************************************** *********************** // |
|
June 4, 2013, 02:29 |
|
#13 |
Member
Yosmcer Mocktai
Join Date: Apr 2013
Location: Behind a computer
Posts: 50
Rep Power: 17 |
Hello, the error log gives you one warning and one error:
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 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. |
|
June 4, 2013, 11:07 |
problem with blockMesh
|
#14 |
New Member
Bhupesh
Join Date: Jun 2012
Location: Germany
Posts: 14
Rep Power: 13 |
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. |
|
June 4, 2013, 14:39 |
some new warnings
|
#15 |
New Member
Bhupesh
Join Date: Jun 2012
Location: Germany
Posts: 14
Rep Power: 13 |
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 ( ); // ************************************************** *********************** // |
|
June 6, 2013, 04:45 |
|
#16 |
Member
Yosmcer Mocktai
Join Date: Apr 2013
Location: Behind a computer
Posts: 50
Rep Power: 17 |
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. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Mesh generation problem | psapkota | Main CFD Forum | 3 | June 30, 2016 13:42 |
[ICEM] surface mesh merging problem | everest | ANSYS Meshing & Geometry | 44 | April 14, 2016 06:41 |
[mesh manipulation] Importing Multiple Meshes | thomasnwalshiii | OpenFOAM Meshing & Mesh Conversion | 18 | December 19, 2015 18:57 |
[snappyHexMesh] No layers in a small gap | bobburnquist | OpenFOAM Meshing & Mesh Conversion | 6 | August 26, 2015 09:38 |
Problem with mesh generation | Nestor | FLUENT | 3 | November 1, 2006 02:54 |