From bb473d2ae3604b15ce8ecd4903bfd31b2f0a38fa Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Fri, 15 Feb 2013 16:55:53 +0000
Subject: [PATCH] ENH: Updated film contact angle force calculation

---
 .../kinematic/force/contactAngleForce/contactAngleForce.C  | 7 +------
 1 file changed, 1 insertion(+), 6 deletions(-)

diff --git a/src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.C b/src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.C
index 9a80bc9e76..140f710a97 100644
--- a/src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.C
+++ b/src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.C
@@ -149,8 +149,6 @@ tmp<fvVectorMatrix> contactAngleForce::correct(volVectorField& U)
 
     volVectorField gradAlpha(fvc::grad(alpha));
 
-    scalarField nHits(owner_.regionMesh().nCells(), 0.0);
-
     forAll(nbr, faceI)
     {
         const label cellO = own[faceI];
@@ -173,7 +171,6 @@ tmp<fvVectorMatrix> contactAngleForce::correct(volVectorField& U)
                 gradAlpha[cellI]/(mag(gradAlpha[cellI]) + ROOTVSMALL);
             scalar theta = cos(degToRad(distribution_->sample()));
             force[cellI] += Ccf_*n*sigma[cellI]*(1.0 - theta)/invDx;
-            nHits[cellI]++;
         }
     }
 
@@ -200,15 +197,13 @@ tmp<fvVectorMatrix> contactAngleForce::correct(volVectorField& U)
                         scalar theta = cos(degToRad(distribution_->sample()));
                         force[cellO] +=
                             Ccf_*n*sigma[cellO]*(1.0 - theta)/invDx[faceI];
-                        nHits[cellO]++;
                     }
                 }
             }
         }
     }
 
-    force /= (max(nHits, scalar(1.0))*magSf);
-    tForce().correctBoundaryConditions();
+    force /= magSf;
 
     if (owner_.regionMesh().time().outputTime())
     {