diff --git a/src/dynamicMesh/motionSmoother/polyMeshGeometry/polyMeshGeometry.C b/src/dynamicMesh/motionSmoother/polyMeshGeometry/polyMeshGeometry.C index f42f2a4ef4130059d8aad52c65824d7e0e8b3a79..f777081d03097dfc72b6307ba8d7306441a1fae7 100644 --- a/src/dynamicMesh/motionSmoother/polyMeshGeometry/polyMeshGeometry.C +++ b/src/dynamicMesh/motionSmoother/polyMeshGeometry/polyMeshGeometry.C @@ -824,7 +824,7 @@ bool Foam::polyMeshGeometry::checkFaceSkewness // (i.e. treat as if mirror cell on other side) vector faceNormal = faceAreas[faceI]; - faceNormal /= mag(faceNormal) + VSMALL; + faceNormal /= mag(faceNormal) + ROOTVSMALL; vector dOwn = faceCentres[faceI] - cellCentres[own[faceI]]; @@ -834,7 +834,7 @@ bool Foam::polyMeshGeometry::checkFaceSkewness scalar skewness = mag(faceCentres[faceI] - faceIntersection) - /(2*mag(dWall) + VSMALL); + /(2*mag(dWall) + ROOTVSMALL); // Check if the skewness vector is greater than the PN vector. // This does not cause trouble but is a good indication of a poor