diff --git a/applications/solvers/lagrangian/reactingParcelFilmFoam/YEqn.H b/applications/solvers/lagrangian/reactingParcelFilmFoam/YEqn.H
index 5eed16aef407f694fd67b644eabba11a14c8a005..db79376614e29b4289c41f54e76c7bdf5c25e9b6 100644
--- a/applications/solvers/lagrangian/reactingParcelFilmFoam/YEqn.H
+++ b/applications/solvers/lagrangian/reactingParcelFilmFoam/YEqn.H
@@ -23,7 +23,7 @@ tmp<fv::convectionScheme<scalar> > mvConvection
             (
                 fvm::ddt(rho, Yi)
               + mvConvection->fvmDiv(phi, Yi)
-              - fvm::laplacian(turbulence->muEff(), Yi)
+              - fvm::laplacian(turbulence->alphaEff(), Yi)
               ==
                 parcels.SYi(i, Yi)
               + surfaceFilm.Srho(i)
diff --git a/applications/solvers/lagrangian/reactingParcelFilmFoam/reactingParcelFilmFoam.C b/applications/solvers/lagrangian/reactingParcelFilmFoam/reactingParcelFilmFoam.C
index fdedebb78550ef5a0be066be06f7be306464668b..9ec100f3fd1d32c5ca46ca13aacd44665b395800 100644
--- a/applications/solvers/lagrangian/reactingParcelFilmFoam/reactingParcelFilmFoam.C
+++ b/applications/solvers/lagrangian/reactingParcelFilmFoam/reactingParcelFilmFoam.C
@@ -67,7 +67,7 @@ int main(int argc, char *argv[])
     while (runTime.run())
     {
         #include "readTimeControls.H"
-        #include "readPISOControls.H"
+        #include "readPIMPLEControls.H"
         #include "compressibleCourantNo.H"
         #include "setMultiRegionDeltaT.H"
         #include "setDeltaT.H"
diff --git a/applications/solvers/lagrangian/reactingParcelFilmFoam/setMultiRegionDeltaT.H b/applications/solvers/lagrangian/reactingParcelFilmFoam/setMultiRegionDeltaT.H
index ed3b584675f94ae862e5e1db06ade682fda0654a..037a09f93d46ec92b2f277039260b0f2dde93cf4 100644
--- a/applications/solvers/lagrangian/reactingParcelFilmFoam/setMultiRegionDeltaT.H
+++ b/applications/solvers/lagrangian/reactingParcelFilmFoam/setMultiRegionDeltaT.H
@@ -43,14 +43,11 @@ if (adjustTimeStep)
 
     const scalar dt0 = runTime.deltaTValue();
 
-    const scalar TFactorFluidLim =
-        min(min(TFactorFluid, 1.0 + 0.1*TFactorFluid), 1.2);
-
     runTime.setDeltaT
     (
         min
         (
-            dt0*min(TFactorFluidLim, TFactorFilm),
+            dt0*min(min(TFactorFluid, TFactorFilm), 1.2),
             maxDeltaT
         )
     );