CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Mesh Utilities (http://www.cfd-online.com/Forums/openfoam-meshing-utilities/)
-   -   SnappyHexMesh Point is not inside Mesh but it is (http://www.cfd-online.com/Forums/openfoam-meshing-utilities/61484-snappyhexmesh-point-not-inside-mesh-but.html)

wolle1982 February 18, 2009 10:19

Hi everyOne, maybe this is
 
Hi everyOne,

maybe this is an easy one. Meshing with snappyHexMesh I sometimes get the error message:

Point (51 35 45) is not inside the mesh or on a face or edge.
Bounding box of the meshhttp://www.cfd-online.com/OpenFOAM_D...lipart/sad.gif-1500 -300 -500) (500 1200 1200)

From function refinementParameters::findCells(const polyMesh&) const
in file autoHexMesh/autoHexMeshDriver/refinementParameters/refinementParameters.C at line 104.

FOAM exiting


it tells me, that the point, that defines the keeping mesh zone, is not inside the initial mesh....but it is. i can try with a whole bunch of other coordinates with the same result (i don't think they all accidentely match a cell-face).

maybe someone of you had the same problem sometime and knows the answer. the case itself is a very simple one. i've already managed to mesh with sHM but this one is bugging me, because there seems no plausible explaination.

beefcake883 August 4, 2009 05:56

hi,

did you find out why this happens? I have the same problem.

andersking August 4, 2009 22:24

They key is in "... or on a face or edge", snappyHexMesh can't have the deciding point lie on a face or edge. Try something like (51.3141592654 35.3141592654 45.3141592654), instead of (51 35 45) - ie. something that is unlikely to come up during the cell divisions.

Cheers
Andrew

madad2005 September 4, 2009 09:18

I'm also having this problem and increasing the number of deciaml places is not helping. Has anyone found a proper workaround for this?

madad2005 September 4, 2009 09:27

...and I've just found the problem! Make sure you have defined your bounding box properly, starting with the face at minX then the face at maxX. Don't know if this is properly explaining what the issue was, but this works for me. See below:

vertices
(
( 750 -750 750)
( 750 750 750)
( 750 750 -750)
( 750 -750 -750)
(-750 -750 750)
(-750 750 750)
(-750 750 -750)
(-750 -750 -750)
);
blocks
(
//hex(0 3 2 1 4 7 6 5) (20 20 20) simpleGrading (1 1 1) <--this is wrong
hex (4 7 6 5 0 3 2 1) (20 20 20) simpleGrading (1 1 1) // this is correct
);

fijinx December 16, 2009 15:33

Worked perfect, thanks!

ziad September 3, 2010 08:59

The proper way to do it is to set up your blockMesh domain as per section 5.3 in the user guide. You should use the setup in figure 5.5 on page U-142.

right-handed coordinate system x1x2x3
x1 direction described by moving from vertex 0 to vertex 1
x2 direction described by moving from vertex 1 to vertex 2
vertices 0, 1, 2, 3 define the plane x3 = 0
vertex 4 is found by moving by moving from vertex 0 in the x3 direction
vertices 5,6 and 7 are similarly found by moving in the x3 direction from vertices 1,2
and 3 respectively


All times are GMT -4. The time now is 03:02.