PairCollision: stabilize normalization of the force direction for low force

Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1948
2015-12-12 09:00:42 +00:00 · 2015-12-12 09:00:42 +00:00 · 35439753ef
commit 35439753ef
parent 71f13ec9dd
1 changed files with 1 additions and 1 deletions
--- a/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/PairCollision.C
+++ b/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/PairCollision.C
@ -246,7 +246,7 @@ void Foam::PairCollision<CloudType>::wallInteraction()

                    vector pW = nearPt - pos;

-                    scalar normalAlignment = normal & pW/mag(pW);
+                    scalar normalAlignment = normal & pW/(mag(pW) + ROOTSMALL);

                    // Find the patchIndex and wallData for WallSiteData object
                    label patchI = patchID[realFaceI - mesh.nInternalFaces()];