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 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()) {