diff --git a/applications/solvers/multiphase/interPhaseChangeFoam/alphaEqnSubCycle.H b/applications/solvers/multiphase/interPhaseChangeFoam/alphaEqnSubCycle.H
index 4acc3c137cc5c35ad3b1aab9ab4855c99a7c48ff..f4643f1df431259d0cd963384c20908d380f544e 100644
--- a/applications/solvers/multiphase/interPhaseChangeFoam/alphaEqnSubCycle.H
+++ b/applications/solvers/multiphase/interPhaseChangeFoam/alphaEqnSubCycle.H
@@ -9,7 +9,12 @@
         phic += (interface.cAlpha()*icAlpha)*fvc::interpolate(mag(U));
     }
 
-    volScalarField divU(fvc::div(phi));
+    volScalarField divU
+    (
+        mesh.moving()
+      ? fvc::div(phi + mesh.phi())
+      : fvc::div(phi)
+    );
 
     if (nAlphaSubCycles > 1)
     {