diff --git a/applications/solvers/combustion/plasmaReactingFoam/YEqn.H b/applications/solvers/combustion/plasmaReactingFoam/YEqn.H
index 69fb948a..70dc185d 100644
--- a/applications/solvers/combustion/plasmaReactingFoam/YEqn.H
+++ b/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)
diff --git a/applications/solvers/combustion/plasmaReactingFoam/createFields.H b/applications/solvers/combustion/plasmaReactingFoam/createFields.H
index b3f3d315..074b0cdd 100644
--- a/applications/solvers/combustion/plasmaReactingFoam/createFields.H
+++ b/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-");