diff --git a/applications/solvers/multiphase/compressibleInterFoam/TEqn.H b/applications/solvers/multiphase/compressibleInterFoam/TEqn.H
index 13e1feb..b5012b2 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/TEqn.H
+++ b/applications/solvers/multiphase/compressibleInterFoam/TEqn.H
@@ -5,7 +5,7 @@
       + fvm::div(rhoPhi, T)
       - fvm::laplacian(mixture.alphaEff(turbulence->mut()), T)
       + (
-            fvc::div(fvc::absolute(phi, U), p)
+            fvc::div(phi, p)
           + fvc::ddt(rho, K) + fvc::div(rhoPhi, K)
         )
        *(
diff --git a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/TEqn.H b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/TEqn.H
new file mode 100644
index 0000000..18476c1
--- /dev/null
+++ b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/TEqn.H
@@ -0,0 +1,23 @@
+{
+    fvScalarMatrix TEqn
+    (
+        fvm::ddt(rho, T)
+      + fvm::div(rhoPhi, T)
+      - fvm::laplacian(mixture.alphaEff(turbulence->mut()), T)
+      + (
+            fvc::div(Uf & mesh.Sf(), p, "div(phi,p)")
+          + fvc::ddt(rho, K) + fvc::div(rhoPhi, K)
+        )
+       *(
+           alpha1/mixture.thermo1().Cv()
+         + alpha2/mixture.thermo2().Cv()
+        )
+    );
+
+    TEqn.relax();
+    TEqn.solve();
+
+    mixture.correct();
+
+    Info<< "min(T) " << min(T).value() << endl;
+}
diff --git a/tutorials/multiphase/compressibleInterDyMFoam/ras/sloshingTank2D/system/fvSchemes b/tutorials/multiphase/compressibleInterDyMFoam/ras/sloshingTank2D/system/fvSchemes
index 197d5f4..72d9af0 100644
--- a/tutorials/multiphase/compressibleInterDyMFoam/ras/sloshingTank2D/system/fvSchemes
+++ b/tutorials/multiphase/compressibleInterDyMFoam/ras/sloshingTank2D/system/fvSchemes
@@ -35,7 +35,7 @@ divSchemes
     div(phi,thermo:rho.air) Gauss linear;
     div(rhoPhi,T)  Gauss vanLeer;
     div(rhoPhi,K)  Gauss linear;
-    div((phi+meshPhi),p)  Gauss linear;
+    div(phi,p)     Gauss linear;
 
     div((muEff*dev2(T(grad(U))))) Gauss linear;
 }