[plasma] add safe guard to prevent negative Te in updateElectronTempe…

…ratureFromEnergyDist
Cantera · Nov 19, 2024 · b1ac6e9 · b1ac6e9
1 parent 72998bb
commit b1ac6e9
Showing 1 changed file with 11 additions and 0 deletions.
diff --git a/src/thermo/PlasmaPhase.cpp b/src/thermo/PlasmaPhase.cpp
@@ -180,6 +180,17 @@ void PlasmaPhase::updateElectronTemperatureFromEnergyDist()
     Eigen::ArrayXd eps52 = m_electronEnergyLevels.pow(5./2.);
     double epsilon_m = 2.0 / 5.0 * numericalQuadrature(m_quadratureMethod,
                                                        m_electronEnergyDist, eps52);
+    if (epsilon_m < 0.0 && m_quadratureMethod == "simpson") {
+        // try trapezoidal method
+        epsilon_m = 2.0 / 5.0 * numericalQuadrature(
+            "trapezoidal", m_electronEnergyDist, eps52);
+    }
+
+    if (epsilon_m < 0.0) {
+        throw CanteraError("PlasmaPhase::updateElectronTemperatureFromEnergyDist",
+            "The electron energy distribution produces negative electron temperature.");
+    }
+
     m_electronTemp = 2.0 / 3.0 * epsilon_m * ElectronCharge / Boltzmann;
 }