From 617e030a5ec36d2036ee07e476e449f57d5f59e2 Mon Sep 17 00:00:00 2001
From: andy <a.heather@opencfd.co.uk>
Date: Fri, 17 Sep 2010 13:20:53 +0100
Subject: [PATCH] ENH: Updated reactingParcelFilmFoam solver

---
 .../solvers/lagrangian/reactingParcelFilmFoam/YEqn.H       | 1 +
 .../solvers/lagrangian/reactingParcelFilmFoam/hsEqn.H      | 7 ++++---
 .../solvers/lagrangian/reactingParcelFilmFoam/pEqn.H       | 2 ++
 .../solvers/lagrangian/reactingParcelFilmFoam/rhoEqn.H     | 1 +
 4 files changed, 8 insertions(+), 3 deletions(-)

diff --git a/applications/solvers/lagrangian/reactingParcelFilmFoam/YEqn.H b/applications/solvers/lagrangian/reactingParcelFilmFoam/YEqn.H
index c687f2035b..0a2c36b7b5 100644
--- a/applications/solvers/lagrangian/reactingParcelFilmFoam/YEqn.H
+++ b/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")
             );
diff --git a/applications/solvers/lagrangian/reactingParcelFilmFoam/hsEqn.H b/applications/solvers/lagrangian/reactingParcelFilmFoam/hsEqn.H
index 156d5505ab..3c76f1384c 100644
--- a/applications/solvers/lagrangian/reactingParcelFilmFoam/hsEqn.H
+++ b/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();
diff --git a/applications/solvers/lagrangian/reactingParcelFilmFoam/pEqn.H b/applications/solvers/lagrangian/reactingParcelFilmFoam/pEqn.H
index d27afcce8a..93764f87de 100644
--- a/applications/solvers/lagrangian/reactingParcelFilmFoam/pEqn.H
+++ b/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();
diff --git a/applications/solvers/lagrangian/reactingParcelFilmFoam/rhoEqn.H b/applications/solvers/lagrangian/reactingParcelFilmFoam/rhoEqn.H
index b9b2d6a13f..2c5d41f4fe 100644
--- a/applications/solvers/lagrangian/reactingParcelFilmFoam/rhoEqn.H
+++ b/applications/solvers/lagrangian/reactingParcelFilmFoam/rhoEqn.H
@@ -36,6 +36,7 @@ Description
       + fvc::div(phi)
       ==
         parcels.Srho()
+      + surfaceFilm.Srho()
     );
 }