diff --git a/src/twoPhaseModels/interfaceProperties/contactAngleModels/dynamic/dynamicContactAngle.C b/src/twoPhaseModels/interfaceProperties/contactAngleModels/dynamic/dynamicContactAngle.C
index 68280e3296..5d0346042f 100644
--- a/src/twoPhaseModels/interfaceProperties/contactAngleModels/dynamic/dynamicContactAngle.C
+++ b/src/twoPhaseModels/interfaceProperties/contactAngleModels/dynamic/dynamicContactAngle.C
@@ -81,8 +81,8 @@ Foam::contactAngleModels::dynamic::cosTheta
     return cos
     (
         theta0_
-      + thetaRec_ - theta0_*max(uCoeff, scalar(0))
-      - thetaAdv_ - theta0_*min(uCoeff, scalar(0))
+      + (thetaRec_ - theta0_)*max(uCoeff, scalar(0))
+      - (thetaAdv_ - theta0_)*min(uCoeff, scalar(0))
     );
 }
 
diff --git a/src/twoPhaseModels/interfaceProperties/contactAngleModels/gravitational/gravitationalContactAngle.C b/src/twoPhaseModels/interfaceProperties/contactAngleModels/gravitational/gravitationalContactAngle.C
index 9ab5ea66ff..a0fe6e12e9 100644
--- a/src/twoPhaseModels/interfaceProperties/contactAngleModels/gravitational/gravitationalContactAngle.C
+++ b/src/twoPhaseModels/interfaceProperties/contactAngleModels/gravitational/gravitationalContactAngle.C
@@ -72,9 +72,9 @@ Foam::contactAngleModels::gravitational::cosTheta
 
     return cos
     (
-        degToRad(theta0_)
-      + degToRad(thetaRec_ - theta0_)*max(uCoeff, scalar(0))
-      - degToRad(thetaAdv_ - theta0_)*min(uCoeff, scalar(0))
+        theta0_
+      + (thetaRec_ - theta0_)*max(uCoeff, scalar(0))
+      - (thetaAdv_ - theta0_)*min(uCoeff, scalar(0))
     );
 }