From f3a9ce3f50e26c54a276bc26cb15732bb8341968 Mon Sep 17 00:00:00 2001 From: Henrik Rusche Date: Wed, 11 Dec 2013 23:15:09 +0100 Subject: [PATCH] BUGFIX: 552fd5fd requires calculation of deflated geometry --- .../polyPatches/constraint/processor/processorPolyPatch.C | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) diff --git a/src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C b/src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C index d83ce96e4..ffc9b39cb 100644 --- a/src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C +++ b/src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C @@ -200,11 +200,13 @@ void Foam::processorPolyPatch::calcGeometry() scalar nbrMagSf = mag(neighbFaceAreas_[facei]); scalar avSf = (magSf + nbrMagSf)/2.0; - if (magSf < ROOTVSMALL && nbrMagSf < ROOTVSMALL) + if (magSf < SMALL && nbrMagSf < SMALL) { // Undetermined normal. Use dummy normal to force separation - // check. (note use of sqrt(VSMALL) since that is how mag - // scales) + // check. + // (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); nbrFaceNormals[facei] = faceNormals[facei]; }