diff --git a/applications/solvers/combustion/plasmaReactingFoam/PhiEqn.H b/applications/solvers/combustion/plasmaReactingFoam/PhiEqn.H
index b73c80e9..8b78a8b2 100644
--- a/applications/solvers/combustion/plasmaReactingFoam/PhiEqn.H
+++ b/applications/solvers/combustion/plasmaReactingFoam/PhiEqn.H
@@ -3,7 +3,7 @@
     (
         fvm::laplacian(Phi)
         + fvc::Sp(rho, qc)/epsilon0
-        + fvc::Sp(eCharge, ne)/epsilon0
+        - fvc::Sp(eCharge, ne)/epsilon0
     );
 
     E = -fvc::grad(Phi);
diff --git a/applications/solvers/combustion/plasmaReactingFoam/UEqn.H b/applications/solvers/combustion/plasmaReactingFoam/UEqn.H
index 36eb9e6a..aa678de6 100644
--- a/applications/solvers/combustion/plasmaReactingFoam/UEqn.H
+++ b/applications/solvers/combustion/plasmaReactingFoam/UEqn.H
@@ -6,7 +6,6 @@
      ==
         rho*g
       + rho*qc*E
-      + eCharge*ne*E
       + fvOptions(rho, U)
     );