diff --git a/src/fvMotionSolver/fvMotionSolvers/displacement/surfaceAlignedSBRStress/surfaceAlignedSBRStressFvMotionSolver.C b/src/fvMotionSolver/fvMotionSolvers/displacement/surfaceAlignedSBRStress/surfaceAlignedSBRStressFvMotionSolver.C
index 098f4fb7a4eb773cb719d8d897f326e94a60b5ff..a888fd225df333e02629a7dd6ee592f1001c0a95 100644
--- a/src/fvMotionSolver/fvMotionSolvers/displacement/surfaceAlignedSBRStress/surfaceAlignedSBRStressFvMotionSolver.C
+++ b/src/fvMotionSolver/fvMotionSolvers/displacement/surfaceAlignedSBRStress/surfaceAlignedSBRStressFvMotionSolver.C
@@ -384,11 +384,9 @@ void Foam::surfaceAlignedSBRStressFvMotionSolver::solve()
         sigmaD_ = magNewSigmaD;
     }
 
-    const surfaceScalarField Df
-    (
-        dimensionedScalar("viscosity", dimViscosity, 1.0)
-       *diffusivity().operator()()
-    );
+    const dimensionedScalar oneViscosity("viscosity", dimViscosity, 1.0);
+
+    const surfaceScalarField Df(oneViscosity*diffusivity().operator()());
 
     pointDisplacement_.boundaryFieldRef().updateCoeffs();
 
@@ -417,7 +415,7 @@ void Foam::surfaceAlignedSBRStressFvMotionSolver::solve()
                 )
             )
           ==
-            fvc::div(sigmaD_)
+            oneViscosity*fvc::div(sigmaD_)
           + fvOptions(cellDisp)
         );