[Tom Lauriks]

Dear all,

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

I’m 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’m having problems with our clusters (I have two clusters available). Intel compilers were used. I’m 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’m trying to mesh a 3D drawing of buildings (part of a city, Antwerp, Belgium), using snappyHexMesh. I’m well aware that my geometry is challening for SHM to mesh. I’ve 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’ll 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’ll paste the error here that I’m getting during the meshing:

[52] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[52] #1  Foam::sigFpe::sigHandler(int) at ??:?
[52] #2  ? in "/lib64/libpthread.so.0"
[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/libc.so.6"
[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;}