diff --git a/src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C b/src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C
index 430d3883036972e6de14064c3df5fbab09ae04e0..9d053dc3de8f5b9b7c0d888284cf47e66b333adb 100644
--- a/src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C
+++ b/src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C
@@ -308,7 +308,8 @@ tmp<Foam::fvVectorMatrix> kinematicSingleLayer::solveMomentum
       + fvm::div(phi_, U_)
      ==
       - USp_
-      - fvm::SuSp(rhoSp_, U_)
+//      - fvm::SuSp(rhoSp_, U_)
+      - rhoSp_*U_
       + forces_.correct(U_)
     );