ENH: Updated coalChemistryFoam solver

applications/solvers/lagrangian/coalChemistryFoam/UEqn.H

@@ -5,8 +5,8 @@
       + turbulence->divDevRhoReff(U)
      ==
         rho.dimensionedInternalField()*g
-      + coalParcels.SU()
-      + limestoneParcels.SU()
+      + coalParcels.SU(U)
+      + limestoneParcels.SU(U)
     );
 
     UEqn.relax();

applications/solvers/lagrangian/coalChemistryFoam/YEqn.H

@@ -25,7 +25,7 @@ tmp<fv::convectionScheme<scalar> > mvConvection
               + mvConvection->fvmDiv(phi, Yi)
               - fvm::laplacian(turbulence->muEff(), Yi)
               ==
-                coalParcels.Srho(i)
+                coalParcels.SYi(i, Yi)
               + kappa*chemistry.RR(i)().dimensionedInternalField()
             );
 

applications/solvers/lagrangian/coalChemistryFoam/hsEqn.H

@@ -6,8 +6,8 @@
       - fvm::laplacian(turbulence->alphaEff(), hs)
      ==
         DpDt
-     +  coalParcels.Sh()
-     +  limestoneParcels.Sh()
+     +  coalParcels.Sh(hs)
+     +  limestoneParcels.Sh(hs)
      +  enthalpySource.Su()
      +  radiation->Shs(thermo)
      +  chemistrySh

applications/solvers/lagrangian/coalChemistryFoam/rhoEqn.H

@@ -35,7 +35,7 @@ Description
         fvm::ddt(rho)
       + fvc::div(phi)
       ==
-        coalParcels.Srho()
+        coalParcels.Srho(rho)
     );
 }