BUG: Identify OPEN edges in featureEdgeMesh - cvMesh conformation issues

otherwise.
This commit is contained in:
graham
2011-02-03 11:37:17 +00:00
parent 7bcf9dda01
commit a05a4ddc4c
4 changed files with 14 additions and 8 deletions

View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2004-2010 OpenCFD Ltd.
\\ / A nd | Copyright (C) 2004-2011 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -120,10 +120,16 @@ Foam::featureEdgeMesh::edgeStatus Foam::featureEdgeMesh::classifyEdge
const vector n0(norms[edNorms[0]]);
const vector n1(norms[edNorms[1]]);
if ((n0 & n1) > cosNormalAngleTol_)
scalar n0n1 = n0 & n1;
if (n0n1 > cosNormalAngleTol_)
{
return FLAT;
}
else if (n0n1 < -cosNormalAngleTol_)
{
return OPEN;
}
else if ((fC0tofC1 & n0) > 0.0)
{
return INTERNAL;

View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2004-2010 OpenCFD Ltd.
\\ / A nd | Copyright (C) 2004-2011 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License

View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2009-2010 OpenCFD Ltd.
\\ / A nd | Copyright (C) 2009-2011 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License

View File

@ -185,8 +185,6 @@ void Foam::conformalVoronoiMesh::buildSurfaceConformation
<< " Number of edge hits " << featureEdgeHits.size()
<< endl;
timeCheck("After initial conformation");
insertSurfacePointPairs
(
surfaceHits,
@ -201,6 +199,8 @@ void Foam::conformalVoronoiMesh::buildSurfaceConformation
"edgeConformationLocations_initial.obj"
);
timeCheck("After initial conformation");
initialTotalHits = surfaceHits.size() + featureEdgeHits.size();
}
@ -321,8 +321,6 @@ void Foam::conformalVoronoiMesh::buildSurfaceConformation
<< " Number of edge hits " << featureEdgeHits.size()
<< endl;
timeCheck("Conformation iteration " + name(iterationNo));
totalHits = surfaceHits.size() + featureEdgeHits.size();
if (totalHits > 0)
@ -345,6 +343,8 @@ void Foam::conformalVoronoiMesh::buildSurfaceConformation
);
}
timeCheck("Conformation iteration " + name(iterationNo));
iterationNo++;
if (iterationNo == maxIterations)