From c48de45842fc9db88e0b2b849bb3ba86df105050 Mon Sep 17 00:00:00 2001
From: ignis <ignis@postech.ac.kr>
Date: Thu, 3 Dec 2015 00:24:02 +0900
Subject: [PATCH] plasmaReactingFoam PhiEqn is added to solve Poisson's
 equation

---
 .../solvers/combustion/plasmaReactingFoam/PhiEqn.H        | 8 ++++++++
 .../combustion/plasmaReactingFoam/plasmaReactingFoam.C    | 1 +
 2 files changed, 9 insertions(+)
 create mode 100644 applications/solvers/combustion/plasmaReactingFoam/PhiEqn.H

diff --git a/applications/solvers/combustion/plasmaReactingFoam/PhiEqn.H b/applications/solvers/combustion/plasmaReactingFoam/PhiEqn.H
new file mode 100644
index 00000000..71e4e1df
--- /dev/null
+++ b/applications/solvers/combustion/plasmaReactingFoam/PhiEqn.H
@@ -0,0 +1,8 @@
+{
+    volScalarField& qc = thermo.qc();
+
+    solve
+    (
+        fvm::laplacian(Phi) + fvc::Sp(rho, qc)/epsilon0
+    );
+}
diff --git a/applications/solvers/combustion/plasmaReactingFoam/plasmaReactingFoam.C b/applications/solvers/combustion/plasmaReactingFoam/plasmaReactingFoam.C
index 831df5c9..5efbb114 100644
--- a/applications/solvers/combustion/plasmaReactingFoam/plasmaReactingFoam.C
+++ b/applications/solvers/combustion/plasmaReactingFoam/plasmaReactingFoam.C
@@ -66,6 +66,7 @@ int main(int argc, char *argv[])
         runTime++;
         Info<< "Time = " << runTime.timeName() << nl << endl;
 
+        #include "PhiEqn.H"
         #include "rhoEqn.H"
 
         while (pimple.loop())