CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[snappyHexMesh] snappyHexMesh deterministic?

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

LinkBack Thread Tools Search this Thread Display Modes
Old   September 16, 2021, 05:14
Default snappyHexMesh deterministic?
Tom Lauriks
Join Date: Apr 2020
Posts: 34
Rep Power: 5
Tom Lauriks is on a distinguished road
Dear all,

I need to know whether the snappyHexMesh process is completely deterministic.

I知 running into errors that do not seem deterministic. There are 2 possible causes for this, which makes it hard for me to identify the actual cause.

Firstly, I知 having problems with our clusters (I have two clusters available). Intel compilers were used. I知 pretty sure that I saw openfoam errors on our clusters, that were not deterministic. By this I mean that for one and the same model, sometimes an error occurred and sometimes not, when running it multiple times.

Secondly, I知 trying to mesh a 3D drawing of buildings (part of a city, Antwerp, Belgium), using snappyHexMesh. I知 well aware that my geometry is challening for SHM to mesh. I致e already removed several artifacts from my drawing, that can cause problems for SHM. It is however impossible to remove all such artifacts, because several hundreds of vertices and edges in the drawing could possibly have to be removed. Therefore, my strategy is, to try to mesh the partly cleaned geometry. If this fails, I値l divide it into smaller parts and so on. In the end, I should be able to identify which parts of the geometry are causing problems for SHM.
The problem however is, that I run into errors, that occur after different times spent for meshing for one and the same geometry part. Also, sometimes one of the geometry parts meshes on one of my clusters and not on another cluster. Can anyone confirm whether SHM is supposed to be completely deterministic?

I値l paste the error here that I知 getting during the meshing:
[52] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[52] #1  Foam::sigFpe::sigHandler(int) at ??:?
[52] #2  ? in "/lib64/"
[52] #3  Foam::tetrahedron<Foam::Vector<double>, Foam::Vector<double> const&>::circumRadius() const at ??:?
[52] #4  Foam::polyMeshTetDecomposition::findSharedBasePoint(Foam::polyMesh const&, int, Foam::Vector<double> const&, double, bool) at ??:?
[52] #5  Foam::polyMeshGeometry::checkFaceTets(bool, double, Foam::polyMesh const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::List<int> const&, Foam::List<Foam::Pair<int> > const&, Foam::HashSet<int, Foam::Hash<int> >*) at ??:?
[52] #6  Foam::motionSmootherAlgo::checkMesh(bool, Foam::polyMesh const&, Foam::dictionary const&, Foam::List<int> const&, Foam::List<Foam::Pair<int> > const&, Foam::HashSet<int, Foam::Hash<int> >&) at ??:?
[52] #7  Foam::motionSmootherAlgo::scaleMesh(Foam::List<int>&, Foam::List<Foam::Pair<int> > const&, Foam::dictionary const&, Foam::dictionary const&, bool, int) at ??:?
[52] #8  Foam::motionSmootherAlgo::scaleMesh(Foam::List<int>&, Foam::List<Foam::Pair<int> > const&, bool, int) at ??:?
[52] #9  Foam::snappySnapDriver::scaleMesh(Foam::snapParameters const&, int, Foam::List<Foam::Pair<int> > const&, Foam::motionSmoother&) at ??:?
[52] #10  Foam::snappySnapDriver::doSnap(Foam::dictionary const&, Foam::dictionary const&, double, double, Foam::snapParameters const&) at ??:?
[52] #11  ? at ??:?
[52] #12  __libc_start_main in "/lib64/"
[52] #13  ? at ??:?
[r1c03cn3:31527:0:31527] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2))
@font-face{ font-family:"Times New Roman"; }@font-face{ font-family:"宋体"; }@font-face{ font-family:"Calibri"; }@font-face{ font-family:"SimSun"; }p.MsoNormal{ mso-style-name:Normal; mso-style-parent:""; margin:0pt; margin-bottom:.0001pt; mso-pagination:none; text-align:justify; text-justify:inter-ideograph; font-family:Calibri; mso-fareast-font-family:SimSun; mso-bidi-font-family:'Times New Roman'; font-size:10,5000pt; mso-font-kerning:1,0000pt; }span.msoIns{ mso-style-type:export-only; mso-style-name:""; text-decoration:underline; text-underline:single; color:blue; }span.msoDel{ mso-style-type:export-only; mso-style-name:""; text-decoration:line-through; color:red; }div.Section0{page:Section0;}
Tom Lauriks is offline   Reply With Quote


cluster, deterministic, intel compiler, meshing, snappyhexmesh

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
[CAD formats] Creating waterproof STL using snappyHexMesh or salome Tobi OpenFOAM Meshing & Mesh Conversion 58 May 13, 2020 07:01
[snappyHexMesh] Running snappyHexMesh in parallel - optimizing peterhess OpenFOAM Meshing & Mesh Conversion 2 January 3, 2018 03:54
[snappyHexMesh] Tutorial crashes: snappyHexMesh floating point exception. jasv OpenFOAM Meshing & Mesh Conversion 4 May 10, 2016 03:55
Strange Results With snappyHexMesh calebamiles OpenFOAM Running, Solving & CFD 0 August 14, 2011 17:02
[snappyHexMesh] stitchMesh and snappyHexMesh gdbaldw OpenFOAM Meshing & Mesh Conversion 0 December 23, 2009 03:09

All times are GMT -4. The time now is 10:00.