|February 27, 2012, 02:16||
snappyHexMesh decomposition error
Join Date: Feb 2012
Posts: 4Rep Power: 6
Hi, I'm pretty new at OpenFOAM and I'm trying to so a dsmcFoam simuation using a mesh created with snappyHexMesh. The mesh looks fine in ParaView but I get warnings when I run dsmcInitialise, and the same warnings when I run the simulation which cause it to crash. The error says:
--> FOAM Warning :
From function Foam::List<Foam::FixedList<Foam::label, 4> >Foam::Cloud<ParticleType>::faceTetIndices(label fI, label cI) const
in file meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C at line 561
No base point for face 25085, 3(3298 2596 2567), produces a valid tet decomposition.
Another warning says:
No base point for face xxxxx, (xxxx)produces a decomposition that has a minimum volume greater than tolerance
I've looked all over the place trying to figure out what's causing this error but with no luck. Does anyone know what this error means or how I might fix it?
|June 25, 2013, 11:33||
Join Date: Mar 2013
Posts: 12Rep Power: 5
hi, Esmaeil and John,
have you found the way to solve this problem above? Could you give me a hint for it ? thank you!
|July 29, 2013, 09:57||
Join Date: Jul 2013
Posts: 4Rep Power: 5
Looks like the problem is with face tets.
If your mesh is generated with snappyHexMesh, try the following:
In snappyhexMeshDict file, look for "minTetQuality".
//- Minimum quality of the tet formed by the face-centre
// and variable base point minimum decomposition triangles and
// the cell centre. Set to very negative number (e.g. -1E30) to
// <0 = inside out tet,
// 0 = flat tet
// 1 = regular tet
For a regular tet, it should be positive. Try setting it to 1. If it does not work, set it to a very low positive number. Say 1e-10.
One way to check whether the issue with face tet is resolved or not before the start of the simulation is to do an "checkMesh -allGeometry". It will report if there are any issues with face tets.
Hope this helps.
|Thread||Thread Starter||Forum||Replies||Last Post|
|GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh||gschaider||OpenFOAM||300||October 29, 2014 19:00|
|Saving ParaFoam views and case||sail||OpenFOAM Paraview & paraFoam||9||November 25, 2011 16:46|
|ParaView for OF-1.6-ext||Chrisi1984||OpenFOAM Installation||0||December 31, 2010 07:42|
|DecomposePar links against liblamso0 with OpenMPI||jens_klostermann||OpenFOAM Bugs||11||June 28, 2007 17:51|
|user defined function||cfduser||CFX||0||April 29, 2006 10:58|