ENH: Updated reactingParcelFilmFoam solver

applications/solvers/lagrangian/reactingParcelFilmFoam/YEqn.H

@@ -26,6 +26,7 @@ tmp<fv::convectionScheme<scalar> > mvConvection
               - fvm::laplacian(turbulence->muEff(), Yi)
               ==
                 parcels.Srho(i)
+              + surfaceFilm.Srho(i)
               + kappa*chemistry.RR(i)().dimensionedInternalField(),
                 mesh.solver("Yi")
             );

applications/solvers/lagrangian/reactingParcelFilmFoam/hsEqn.H

@@ -6,9 +6,10 @@
       - fvm::laplacian(turbulence->alphaEff(), hs)
      ==
         DpDt
-     +  parcels.Sh()
-     +  radiation->Shs(thermo)
-     +  chemistrySh
+      + parcels.Sh()
+      + surfaceFilm.Sh()
+      + radiation->Shs(thermo)
+      + chemistrySh
     );
 
     hsEqn.relax();

applications/solvers/lagrangian/reactingParcelFilmFoam/pEqn.H

@@ -24,6 +24,7 @@ if (transonic)
           - fvm::laplacian(rho*rUA, p)
          ==
             parcels.Srho()
+          + surfaceFilm.Srho()
         );
 
         pEqn.solve();
@@ -52,6 +53,7 @@ else
           - fvm::laplacian(rho*rUA, p)
          ==
             parcels.Srho()
+          + surfaceFilm.Srho()
         );
 
         pEqn.solve();

applications/solvers/lagrangian/reactingParcelFilmFoam/rhoEqn.H

@@ -36,6 +36,7 @@ Description
       + fvc::div(phi)
       ==
         parcels.Srho()
+      + surfaceFilm.Srho()
     );
 }