diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/turbulenceModel/kEpsilon.H b/applications/solvers/multiphase/twoPhaseEulerFoam/turbulenceModel/kEpsilon.H
index 016d1f4488..e70cf92e2a 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/turbulenceModel/kEpsilon.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/turbulenceModel/kEpsilon.H
@@ -16,7 +16,12 @@ if(turbulence)
     (
         fvm::ddt(beta, epsilon)
       + fvm::div(phib, epsilon)
-      - fvm::laplacian(alphaEps*nuEffb, epsilon)
+      - fvm::laplacian
+        (
+            alphaEps*nuEffb,
+            epsilon,
+            "laplacian(alphaEps*nuEffb,epsilon)"
+        )
       ==
          C1*beta*G*epsilon/k
        - fvm::Sp(C2*beta*epsilon/k, epsilon)