plasmaReactingFoam with Sc != 1 diffusivity

applications/solvers/combustion/plasmaReactingFoam/YEqn.H

@@ -15,6 +15,8 @@ tmp<fv::convectionScheme<scalar> > mvConvection
     label inertIndex = -1;
     volScalarField Yt(0.0*Y[0]);
 
+    composition.calculateDiffusivities(p, T);
+
     forAll(rho, celli)
     {
         const scalar Eni = En[celli] * EnFac;
@@ -33,6 +35,8 @@ tmp<fv::convectionScheme<scalar> > mvConvection
 
     forAll(Y, i)
     {
+        const volScalarField& Di = D[i];
+
         if (Y[i].name() == electronSpecie)
         {
             volScalarField& Yi = Y[i];
@@ -94,7 +98,8 @@ tmp<fv::convectionScheme<scalar> > mvConvection
             (
                 fvm::ddt(rho, Yi)
               + mvConvection->fvmDiv(phi, Yi)
-              - fvm::laplacian(turbulence->muEff(), Yi)
+              // - fvm::laplacian(turbulence->muEff(), Yi)
+              - fvm::laplacian(rho*Di, Yi)
              ==
                 reaction->R(Yi)
               + fvOptions(rho, Yi)

applications/solvers/combustion/plasmaReactingFoam/createFields.H

@@ -118,6 +118,7 @@ const dimensionedScalar R = constant::physicoChemical::R;
 const scalar pi = constant::mathematical::pi;
 
 PtrList<volScalarField>& Y = composition.Y();
+const PtrList<volScalarField>& D = composition.D();
 
 word inertSpecie(thermo.lookup("inertSpecie"));
 word electronSpecie("E-");