diff --git a/applications/solvers/multiphase/interDyMFoam/UEqn.H b/applications/solvers/multiphase/interDyMFoam/UEqn.H
index 7e4a29468e6dac4c28d416bfed9e2607a8abf5fc..4c14afe1a64410d8f621e10d51899e18a6433d8c 100644
--- a/applications/solvers/multiphase/interDyMFoam/UEqn.H
+++ b/applications/solvers/multiphase/interDyMFoam/UEqn.H
@@ -1,29 +1,34 @@
-    surfaceScalarField muf = twoPhaseProperties.muf();
+surfaceScalarField muEff
+(
+    "muEff",
+    twoPhaseProperties.muf()
+  + fvc::interpolate(rho*turbulence->nut())
+);
 
-    fvVectorMatrix UEqn
-    (
-        fvm::ddt(rho, U)
-      + fvm::div(rhoPhi, U)
-      - fvm::laplacian(muf, U)
-      - (fvc::grad(U) & fvc::grad(muf))
-    //- fvc::div(muf*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf()))
-    );
+fvVectorMatrix UEqn
+(
+    fvm::ddt(rho, U)
+  + fvm::div(rhoPhi, U)
+  - fvm::laplacian(muEff, U)
+  - (fvc::grad(U) & fvc::grad(muEff))
+//- fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf()))
+);
 
-    UEqn.relax();
+UEqn.relax();
 
-    if (momentumPredictor)
-    {
-        solve
+if (momentumPredictor)
+{
+    solve
+    (
+        UEqn
+      ==
+        fvc::reconstruct
         (
-            UEqn
-         ==
-            fvc::reconstruct
             (
-                (
-                    fvc::interpolate(interface.sigmaK())*fvc::snGrad(gamma)
-                  - ghf*fvc::snGrad(rho)
-                  - fvc::snGrad(pd)
-                ) * mesh.magSf()
-            )
-        );
-    }
+                fvc::interpolate(interface.sigmaK())*fvc::snGrad(gamma)
+              - ghf*fvc::snGrad(rho)
+              - fvc::snGrad(pd)
+            )*mesh.magSf()
+        )
+    );
+}