ENH: Updated porousExplicitSource...Foam solver

applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/YEqn.H

@@ -11,6 +11,7 @@ tmp<fv::convectionScheme<scalar> > mvConvection
 );
 
 
+if (solveSpecies)
 {
     label inertIndex = -1;
     volScalarField Yt(0.0*Y[0]);

applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/hsEqn.H

@@ -16,17 +16,16 @@
         )
     );
 
-    if (dpdt)
+    if (pressureWork)
     {
-        pWork() += fvc::ddt(p);
-    }
-    if (eWork)
-    {
-        pWork() = -p*fvc::div(phi/fvc::interpolate(rho));
-    }
-    if (hWork)
-    {
-        pWork() += fvc::div(phi/fvc::interpolate(rho)*fvc::interpolate(p));
+        surfaceScalarField phiU("phiU", phi/fvc::interpolate(rho));
+
+        pWork() += fvc::div(phiU*fvc::interpolate(p)) - p*fvc::div(phiU);
+
+        if (pressureWorkTimeDerivative)
+        {
+            pWork() += fvc::ddt(p);
+        }
     }
 
     {

applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/readAdditionalSolutionControls.H

 dictionary additional = mesh.solutionDict().subDict("additional");
 
-bool dpdt = additional.lookupOrDefault("dpdt", true);
-bool eWork = additional.lookupOrDefault("eWork", true);
-bool hWork = additional.lookupOrDefault("hWork", true);
+// pressure work term for enthalpy equation
+bool pressureWork = additional.lookupOrDefault("pressureWork", true);
+bool pressureWorkTimeDerivative =
+    additional.lookupOrDefault("pressureWorkTimeDerivative", true);
+
+// flag to activate solve transport for each specie (Y vector)
+bool solveSpecies = additional.lookupOrDefault("solveSpecies", true);