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 comment to collapse exis determined by eigenvector and

Browse Source 
reordering.
This commit is contained in:
graham
2010-01-08 11:42:20 +00:00
parent 67d2158c2f
commit 7800ceac99
1 changed files with 6 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -821,14 +821,18 @@ bool Foam::conformalVoronoiMesh::collapseFaceToEdge
{ {
vector eVals = eigenValues(J); vector eVals = eigenValues(J);
// The maximum eigenvalue (z()) must be the direction of the
// normal, as it has the greatest value. The minimum eigenvalue
// is the dominant collapse axis for high aspect ratio faces.
collapseAxis = eigenVector(J, eVals.x());
// The inertia calculation describes the mass distribution as a // The inertia calculation describes the mass distribution as a
// function of distance squared to the axis, so the square root of // function of distance squared to the axis, so the square root of
// the ratio of face-plane moments gives a good indication of the // the ratio of face-plane moments gives a good indication of the
// aspect ratio. // aspect ratio.
aspectRatio = sqrt(eVals.y()/max(eVals.x(), SMALL)); aspectRatio = sqrt(eVals.y()/max(eVals.x(), SMALL));
collapseAxis = eigenVector(J, eVals.x());
} }
if (magSqr(collapseAxis) < VSMALL) if (magSqr(collapseAxis) < VSMALL)