lagrangian solver updates

applications/solvers/combustion/coalChemistryFoam/hEqn.H

@@ -2,7 +2,7 @@
     fvScalarMatrix hEqn
     (
         fvm::ddt(rho, h)
-      + fvm::div(phi, h)
+      + mvConvection->fvmDiv(phi, h)
       - fvm::laplacian(turbulence->alphaEff(), h)
      ==
         DpDt

applications/solvers/combustion/dieselEngineFoam/YEqn.H

@@ -14,7 +14,7 @@ tmp<fv::convectionScheme<scalar> > mvConvection
     label inertIndex = -1;
     volScalarField Yt = 0.0*Y[0];
 
-    for(label i=0; i<Y.size(); i++)
+    for (label i=0; i<Y.size(); i++)
     {
         if (Y[i].name() != inertSpecie)
         {
@@ -39,8 +39,7 @@ tmp<fv::convectionScheme<scalar> > mvConvection
             inertIndex = i;
         }
     }
-    
+
     Y[inertIndex] = scalar(1) - Yt;
     Y[inertIndex].max(0.0);
-
 }

applications/solvers/lagrangian/reactingParcelFoam/YEqn.H

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