From 8c70f4af8829096fd5dc639f91151fddab2ec2b8 Mon Sep 17 00:00:00 2001
From: ignis <ignis@postech.ac.kr>
Date: Thu, 13 Oct 2016 17:13:03 +0900
Subject: [PATCH] plasmaReactingFoam with Sc != 1 diffusivity

---
 applications/solvers/combustion/plasmaReactingFoam/YEqn.H  | 7 ++++++-
 .../solvers/combustion/plasmaReactingFoam/createFields.H   | 1 +
 2 files changed, 7 insertions(+), 1 deletion(-)

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-");