zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Adding check of mesh after initial point merge (which gets rid of zero

Browse Source 
area faces).  The intention is that if checkMesh returns a non zero
number of bad faces, then no amount of collapsing and backtracking can
satisfy the criteria (as long as the value of mergeClosenessCoeff is
appropriate - a more rigorous zero area face collapse may be
appropriate) and the mesher will halt (or at least not filter at
all). The mesh quality criteria and the surface conformation controls
can be altered and the internalDelaunayVertices read in to try again,
or more cell size resolution can be applied in areas where problems
are occuring.
This commit is contained in:
graham
2010-01-12 14:11:06 +00:00
parent 9f9e97a97f
commit 4e98d8b99e
1 changed files with 33 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

48
src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshCalcDualMesh.C
View File

@ -109,29 +109,48 @@ void Foam::conformalVoronoiMesh::calcDualMesh
points.setSize(dualVertI);
// Merge close points
// No-risk face filtering to get rid of zero area faces and
// establish if the mesh can be produced at all to the
// specified criteria
Info<< nl << " Merging close points" << endl;
mergeCloseDualVertices(points);
// Smooth the surface of the mesh
label nInitialBadQualityFaces = checkPolyMeshQuality(points);
Info<< nl << " Smoothing surface" << endl;
smoothSurface(points);
// Collapse faces throughout the mesh
Info<< nl << " Collapsing unnecessary faces" << endl;
Info<< "Initial check before face collapse, found "
<< nInitialBadQualityFaces << " bad quality faces"
<< endl;
HashSet<labelPair, labelPair::Hash<> > deferredCollapseFaces;
collapseFaces(points, deferredCollapseFaces);
if (nInitialBadQualityFaces == 0)
{
// Risky and undo-able face filtering to reduce the face count
// as much as possible staying within the specified criteria
label nBadQualityFaces = checkPolyMeshQuality(points);
Info<< nl << " Smoothing surface" << endl;
Info<< "Found " << nBadQualityFaces << " bad quality faces" << endl;
smoothSurface(points);
Info<< nl << " Collapsing unnecessary faces" << endl;
collapseFaces(points, deferredCollapseFaces);
label nBadQualityFaces = checkPolyMeshQuality(points);
Info<< "Found " << nBadQualityFaces << " bad quality faces" << endl;
}
else
{
Info<< "A mesh could not be produced to satisfy the specified quality "
<< "criteria. The quality and the surface conformation controls "
<< "can be altered and the internalDelaunayVertices read in to try "
<< "again, or more cell size resolution and motion iterations can "
<< "be applied in areas where problems are occurring."
<< endl;
}
// Final dual face and owner neighbour construction
@ -148,7 +167,7 @@ void Foam::conformalVoronoiMesh::calcDualMesh
false
);
deferredCollapseFaceSet(owner, neighbour, deferredCollapseFaces);
// deferredCollapseFaceSet(owner, neighbour, deferredCollapseFaces);
removeUnusedCells(owner, neighbour);
@ -372,8 +391,7 @@ void Foam::conformalVoronoiMesh::mergeCloseDualVertices(const pointField& pts)
if (nPtsMerged > 0)
{
Info<< " Merged " << nPtsMerged << " points "
<< "closenessTolerance HARDCODED " << endl;
Info<< " Merged " << nPtsMerged << " points " << endl;
}
reindexDualVertices(dualPtIndexMap);