exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

BUGFIX: 552fd5fd requires calculation of deflated geometry

Browse source
This commit is contained in:
Henrik Rusche 2013-12-11 23:15:09 +01:00
parent c156c308f2
commit f3a9ce3f50
1 changed files with 5 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C
View file

@ -200,11 +200,13 @@ void Foam::processorPolyPatch::calcGeometry()
scalar nbrMagSf = mag(neighbFaceAreas_[facei]); scalar nbrMagSf = mag(neighbFaceAreas_[facei]);
scalar avSf = (magSf + nbrMagSf)/2.0; scalar avSf = (magSf + nbrMagSf)/2.0;
if (magSf < ROOTVSMALL && nbrMagSf < ROOTVSMALL) if (magSf < SMALL && nbrMagSf < SMALL)
{ {
// Undetermined normal. Use dummy normal to force separation // Undetermined normal. Use dummy normal to force separation
// check. (note use of sqrt(VSMALL) since that is how mag // check.
// scales) // (note use of sqrt(VSMALL) since that is how mag scales)
// HR, 11/12/2013: Found a face with area = 1e-21 before
// topo deflation. Hence must use SMALL here.
faceNormals[facei] = point(1, 0, 0); faceNormals[facei] = point(1, 0, 0);
nbrFaceNormals[facei] = faceNormals[facei]; nbrFaceNormals[facei] = faceNormals[facei];
} }