diff --git a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
index b4259a00be9d6af8939a83d6f3790cc5249649af..e30881ac436fccde3f8f6f14c61c9352fcf55e36 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
+++ b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
@@ -123,8 +123,6 @@ int main(int argc, char *argv[])
 
                 // Make the fluxes relative to the mesh motion
                 fvc::makeRelative(phi, U);
-
-                mesh.topoChanging(false);
             }
         }
 
diff --git a/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C b/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
index bf7f7ec1b4f40a75f3f32775f18c63970f187288..366369797614ed08b930373d570a5dd3a2827a2d 100644
--- a/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
+++ b/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
@@ -135,7 +135,7 @@ int main(int argc, char *argv[])
                 if ((mesh.changing() && correctPhi) || mesh.topoChanging())
                 {
                     // Calculate absolute flux from the mapped surface velocity
-                    // SAF: temporary fix until mapped Uf is assessed
+                    // Note: temporary fix until mapped Uf is assessed
                     Uf = fvc::interpolate(U);
 
                     phi = mesh.Sf() & Uf;
@@ -146,8 +146,6 @@ int main(int argc, char *argv[])
                     fvc::makeRelative(phi, U);
 
                     mixture.correct();
-
-                    mesh.topoChanging(false);
                 }
 
                 if (mesh.changing() && checkMeshCourantNo)