From 993b7af513fd7de219c0c1ad6eaf1f0941446c1e Mon Sep 17 00:00:00 2001
From: Kutalmis Bercin <kutalmis.bercin@esi-group.com>
Date: Tue, 6 Dec 2022 15:52:13 +0000
Subject: [PATCH] BUG: contactAngleForce: avoid suppresing the film entering
 the domain (fixes #2651)

---
 .../contactAngleForce/contactAngleForce.C         | 15 +++++++++++++++
 1 file changed, 15 insertions(+)

diff --git a/src/regionFaModels/liquidFilm/subModels/kinematic/force/contactAngleForces/contactAngleForce/contactAngleForce.C b/src/regionFaModels/liquidFilm/subModels/kinematic/force/contactAngleForces/contactAngleForce/contactAngleForce.C
index ce46bad3d8..8a8edd0f8a 100644
--- a/src/regionFaModels/liquidFilm/subModels/kinematic/force/contactAngleForces/contactAngleForce/contactAngleForce.C
+++ b/src/regionFaModels/liquidFilm/subModels/kinematic/force/contactAngleForces/contactAngleForce/contactAngleForce.C
@@ -167,6 +167,21 @@ tmp<faVectorMatrix> contactAngleForce::correct(areaVectorField& U)
         }
     }
 
+    for (const faPatchScalarField& sigmaBf : sigma.boundaryField())
+    {
+        const faPatch& p = sigmaBf.patch();
+
+        if (!p.coupled())
+        {
+            const labelUList& faces = p.edgeFaces();
+
+            forAll(sigmaBf, edgei)
+            {
+                const label face0 = faces[edgei];
+                force[face0] = Zero;
+            }
+        }
+    }
 
     if (film().regionMesh().time().writeTime())
     {