|
[Sponsors] |
[snappyHexMesh] SnappyHexMesh Point is not inside Mesh but it is |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 18, 2009, 10:19 |
SnappyHexMesh Point is not inside Mesh but it is
|
#1 |
Senior Member
Wolfgang Heydlauff
Join Date: Mar 2009
Location: Germany
Posts: 136
Rep Power: 21 |
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 mesh-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. |
|
August 4, 2009, 06:56 |
|
#2 |
New Member
Join Date: Aug 2009
Posts: 1
Rep Power: 0 |
hi,
did you find out why this happens? I have the same problem. |
|
August 4, 2009, 23:24 |
|
#3 |
Member
Andrew King
Join Date: Mar 2009
Location: Perth, Western Australia, Australia
Posts: 82
Rep Power: 17 |
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
__________________
Dr Andrew King Fluid Dynamics Research Group Curtin University |
|
September 4, 2009, 10:18 |
|
#4 |
Senior Member
Anonymous
Join Date: Mar 2009
Posts: 110
Rep Power: 17 |
I'm also having this problem and increasing the number of deciaml places is not helping. Has anyone found a proper workaround for this?
|
|
September 4, 2009, 10:27 |
|
#5 |
Senior Member
Anonymous
Join Date: Mar 2009
Posts: 110
Rep Power: 17 |
...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 ); |
|
December 16, 2009, 15:33 |
|
#6 |
Member
James Baker
Join Date: Dec 2009
Posts: 35
Rep Power: 17 |
Worked perfect, thanks!
|
|
September 3, 2010, 09:59 |
|
#7 |
Senior Member
Join Date: Mar 2009
Location: My oyster
Posts: 124
Rep Power: 17 |
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 |
|
July 12, 2016, 11:46 |
|
#8 | |
New Member
Sandip Wadekar
Join Date: Oct 2014
Posts: 17
Rep Power: 12 |
Quote:
Thank you |
||
June 19, 2017, 08:37 |
|
#9 |
New Member
ABHINAVMOHAN
Join Date: Jun 2017
Posts: 16
Rep Power: 9 |
vertices
( (48 660 -2015) (2000 660 -2015) (2000 3000 -2015) (48 3000 -2015) (48 660 4000) (2000 660 4000) (2000 3000 4000) (48 3000 4000) ); blocks ( hex (3 2 6 7 0 1 5 4) (20 20 20) simpleGrading (1 1 1) Same problem I guess tried everything still can't figure it out. |
|
October 4, 2017, 08:38 |
|
#10 |
New Member
Filippo
Join Date: Jun 2017
Posts: 1
Rep Power: 0 |
Hi abhinavmohan!
I think that the correct order of the points in blockMeshDict should be: blocks ( hex (0 1 2 3 4 5 6 7) (20 20 20) simpleGrading (1 1 1) ); Hope it helps! |
|
October 4, 2017, 14:25 |
|
#11 |
New Member
Thomas Evans
Join Date: Dec 2015
Posts: 21
Rep Power: 11 |
abhinavmohan:
What happened is, you did not form the volume you thought you were forming. Which is why you should run checkMesh on the grid when any tools complain about the mesh. Or even beforehand. The order of points matters, as illustrated in the user guide, section 5.2, table 5.1, "Vertex, face and edge numbering for cellShapes." But I guess "trying everything" doesn't include reading the manual, huh? |
|
December 1, 2019, 15:11 |
|
#12 |
New Member
Ashwin
Join Date: Mar 2018
Location: Houston Texas
Posts: 10
Rep Power: 8 |
Sometimes it is better to clear the polymesh folder and remesh.
For example, all my definitions were correct, (location in mesh, ordering of blockmesh etc) but was still failing (with the message at the begining of this thread). I then deleted polymesh and remeshed (with exactly the same definitions as before) and it worked correctly. It appears that the old polymesh definitions were retained and not cleared. Nevertheless, if everything fails, one thing to try is to clear out the polymesh folder and restart the meshing process. |
|
November 15, 2020, 16:06 |
|
#13 |
Member
Deutschland
Join Date: Sep 2020
Posts: 69
Rep Power: 6 |
hey,
I have been trying to do https://holzmann-cfd.com/community/t.../suzannes-head example in order to learn how creating geometry using blender / Salome and using it in the OpenFOAM simulation. I created my own background mesh using Salome by following the tutorial shown in the link and when I try to run the command "mpirun -np 4 snappyHexMesh -parallel -overwrite" I got the following error Code:
-------------------------------------------------------------------------- WARNING: Linux kernel CMA support was requested via the btl_vader_single_copy_mechanism MCA variable, but CMA support is not available due to restrictive ptrace settings. The vader shared memory BTL will fall back on another single-copy mechanism if one is available. This may result in lower performance. Local host: LAPTOP-BKEON11O -------------------------------------------------------------------------- /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2006 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : _b45f8f6f58-20200629 OPENFOAM=2006 Arch : "LSB;label=32;scalar=64" Exec : snappyHexMesh -parallel -overwrite Date : Nov 15 2020 Time : 20:39:29 Host : LAPTOP-BKEON11O PID : 644 I/O : uncollated Case : /home/sam/OpenFOAM/OpenFOAM-v2006/tutorials/susan nProcs : 4 Hosts : ( (LAPTOP-BKEON11O 4) ) Pstream initialized with: floatTransfer : 0 nProcsSimpleSum : 0 commsType : nonBlocking polling iterations : 0 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Read mesh in = 0.02 s Overall mesh bounding box : (-1 -0.827093 -0.601021) (1 2.36662 0.671601) Relative tolerance : 1e-06 Absolute matching distance : 3.97736e-06 Reading refinement surfaces. Read refinement surfaces in = 0.01 s Reading refinement shells. Read refinement shells in = 0 s Setting refinement level of surface to be consistent with shells. For geometry suzanne.stl detected 0 uncached triangles out of 12100 For geometry leftears.stl detected 0 uncached triangles out of 1872 For geometry rightears.stl detected 0 uncached triangles out of 1872 Checked shell refinement in = 0.06 s Reading features. Read features in = 0 s Determining initial surface intersections ----------------------------------------- Edge intersection testing: Number of edges : 78797 Number of edges to retest : 0 Number of intersected edges : 0 Edge intersection testing: Number of edges : 78797 Number of edges to retest : 78797 Number of intersected edges : 26 Calculated surface intersections in = 0.06 s Initial mesh : cells:25346 faces:78797 points:28200 Cells per refinement level: 0 25346 Adding patches for surface regions ---------------------------------- Patch Type Region ----- ---- ------ monkey: 3 wall monkey monkeyLBack: 4 wall monkeyLBack monkeyRBack: 5 wall monkeyRBack Added patches in = 0 s Edge intersection testing: Number of edges : 78797 Number of edges to retest : 0 Number of intersected edges : 26 Selecting decompositionMethod simple [4] Refinement phase ---------------- [0] [0] [1] [1] [1] --> FOAM FATAL ERROR: [1] Point (1.2 0.1 0) is not inside the mesh or on a face or edge. Bounding box of the mesh:(-1 -0.827093 -0.601021) (1 2.36662 0.671601) [2] [2] [2] --> FOAM FATAL ERROR: [2] Point (1.2 0.1 0) is not inside the mesh or on a face or edge. Bounding box of the mesh:(-1 -0.827093 -0.601021) (1 2.36662 0.671601) [2] [2] From static Foam::labelList Foam::refinementParameters::findCells(bool, const Foam::polyMesh&, const pointField&) [2] in file [3] [3] [3] --> FOAM FATAL ERROR: [3] Point (1.2 0.1 0) is not inside the mesh or on a face or edge. Bounding box of the mesh:(-1 -0.827093 -0.601021) (1 2.36662 0.671601) [3] [3] From static Foam::labelList Foam::refinementParameters::findCells(bool, const Foam::polyMesh&, const pointField&) [3] in file snappyHexMeshDriver/refinementParameters/refinementParameters.C at line 241.[0] --> FOAM FATAL ERROR: [0] Point (1.2 0.1 0) is not inside the mesh or on a face or edge. Bounding box of the mesh:(-1 -0.827093 -0.601021) (1 2.36662 0.671601) [0] [0] From static Foam::labelList Foam::refinementParameters::findCells(bool, const Foam::polyMesh&, const pointField&) [0] in file snappyHexMeshDriver/refinementParameters/refinementParameters.C at line 241. [0] FOAM parallel run exiting [0] [1] [1] From static Foam::labelList Foam::refinementParameters::findCells(bool, const Foam::polyMesh&, const pointField&) [1] in file snappyHexMeshDriver/refinementParameters/refinementParameters.C at line 241. [1] FOAM parallel run exiting [1] snappyHexMeshDriver/refinementParameters/refinementParameters.C at line 241. [2] FOAM parallel run exiting [2] -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 2 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see the output from other processes, depending on exactly when Open MPI kills them. -------------------------------------------------------------------------- [3] FOAM parallel run exiting [3] [LAPTOP-BKEON11O:00640] 3 more processes have sent help message help-btl-vader.txt / cma-permission-denied [LAPTOP-BKEON11O:00640] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages [LAPTOP-BKEON11O:00640] 3 more processes have sent help message help-mpi-api.txt / mpi-abort and How do I fix it? Kind regards vava10 |
|
November 15, 2020, 22:05 |
|
#14 | |
Senior Member
Join Date: Aug 2013
Posts: 407
Rep Power: 16 |
Hi,
Quote:
A quick fix would be to change your "locationInMesh" point to (0.2 0.1 0) and see if the resulting mesh is the one you want. Hope this helps. Cheers, Antimony |
||
March 25, 2024, 06:46 |
--> FOAM FATAL ERROR: Neither insidePoint/insidePoints nor outsidePoint/outsidePoin
|
#15 |
New Member
Join Date: Feb 2024
Posts: 3
Rep Power: 2 |
im getting the above error, please help me rectify it. thank you https://drive.google.com/file/d/1y5m...usp=drive_link
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[snappyHexMesh] snappyHexMesh can't create refimenet region | clktp | OpenFOAM Meshing & Mesh Conversion | 3 | February 27, 2022 10:26 |
Gambit problems | Althea | FLUENT | 22 | January 4, 2017 04:19 |
[snappyHexMesh] Error with a turbine mesh with SnappyHexMesh | DOliveira | OpenFOAM Meshing & Mesh Conversion | 1 | July 16, 2015 09:02 |
[snappyHexMesh] snappyHexMesh & Mesh around a missile | sasanghomi | OpenFOAM Meshing & Mesh Conversion | 2 | October 15, 2014 20:54 |
Random point inside the mesh | filipe | OpenFOAM Programming & Development | 0 | October 16, 2009 12:24 |