muf -> muEff for consistency

applications/solvers/multiphase/compressibleInterDyMFoam/UEqn.H

-    surfaceScalarField muf
+    surfaceScalarField muEff
     (
-        "muf",
+        "muEff",
         twoPhaseProperties.muf()
       + fvc::interpolate(rho*turbulence->nut())
     );
@@ -9,8 +9,8 @@
     (
         fvm::ddt(rho, U)
       + fvm::div(rhoPhi, U)
-      - fvm::laplacian(muf, U)
-      - (fvc::grad(U) & fvc::grad(muf))
+      - fvm::laplacian(muEff, U)
+      - (fvc::grad(U) & fvc::grad(muEff))
     //- fvc::div(muf*(mesh.Sf() & fvc::interpolate(fvc::grad(U)().T())))
     );
 

applications/solvers/multiphase/interPhaseChangeFoam/UEqn.H

-    surfaceScalarField muf
+    surfaceScalarField muEff
     (
-        "muf",
+        "muEff",
         twoPhaseProperties->muf()
       + fvc::interpolate(rho*turbulence->nut())
     );
@@ -10,9 +10,9 @@
         fvm::ddt(rho, U)
       + fvm::div(rhoPhi, U)
       - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
-      - fvm::laplacian(muf, U)
-      - (fvc::grad(U) & fvc::grad(muf))
-    //- fvc::div(muf*(fvc::interpolate(dev2(fvc::grad(U))) & mesh.Sf()))
+      - fvm::laplacian(muEff, U)
+      - (fvc::grad(U) & fvc::grad(muEff))
+    //- fvc::div(muEff*(fvc::interpolate(dev2(fvc::grad(U))) & mesh.Sf()))
     );
 
     UEqn.relax();