From ad9cb13d57f49f58ec68f803d038d875bf8e24ec Mon Sep 17 00:00:00 2001
From: Yeongdo Park <ignis@postech.ac.kr>
Date: Sat, 16 Feb 2019 10:05:10 -0500
Subject: [PATCH] diffusivityModel mobility method

---
 .../diffusivityModel/diffusivityModel.C        | 18 ++++++++++++++++++
 .../diffusivityModel/diffusivityModel.H        |  3 +++
 2 files changed, 21 insertions(+)

diff --git a/diffusivityModel/diffusivityModel/diffusivityModel.C b/diffusivityModel/diffusivityModel/diffusivityModel.C
index 9105991..0c1e6b6 100644
--- a/diffusivityModel/diffusivityModel/diffusivityModel.C
+++ b/diffusivityModel/diffusivityModel/diffusivityModel.C
@@ -811,6 +811,24 @@ Foam::tmp<Foam::scalarField> Foam::diffusivityModel::kpure(const scalar T, const
 }
 
 
+Foam::tmp<Foam::volScalarField> Foam::diffusivityModel::mu(const label i, const volScalarField& T) const
+{
+    const speciesTable &species_(thermo_.composition().species());
+
+    const word namei(species_[i]);
+    const label zi = thermo_.composition().z(i);
+
+    // - Elementary charge (default in [C])
+    const dimensionedScalar eCharge = constant::electromagnetic::e;
+
+    //- Boltzmann constant
+    const dimensionedScalar kB = constant::physicoChemical::k;
+
+
+    return (this->D(i)/T * (eCharge*zi/kB));
+}
+
+
 // * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * * //
 
 void Foam::diffusivityModel::operator=(const diffusivityModel& rhs)
diff --git a/diffusivityModel/diffusivityModel/diffusivityModel.H b/diffusivityModel/diffusivityModel/diffusivityModel.H
index 7eec9e8..7f0d654 100644
--- a/diffusivityModel/diffusivityModel/diffusivityModel.H
+++ b/diffusivityModel/diffusivityModel/diffusivityModel.H
@@ -225,6 +225,9 @@ public:
 
         // Edit
 
+        tmp<volScalarField> mu(const label i, const volScalarField& T) const;
+
+
         // Write