[1D] Prevent two-point control from finding solutions with reverse flow

Fixes #1787
Cantera · Nov 16, 2024 · a3755b5 · a3755b5
1 parent 70521b7
commit a3755b5
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/src/oneD/Flow1D.cpp b/src/oneD/Flow1D.cpp
@@ -1258,6 +1258,9 @@ void Flow1D::enableTwoPointControl(bool twoPointControl)
 {
     if (isStrained()) {
         m_twoPointControl = twoPointControl;
+        // Prevent finding spurious solutions with negative velocity (outflow) at either
+        // inlet.
+        setBounds(c_offset_V, -1e-5, 1e20);
     } else {
         throw CanteraError("Flow1D::enableTwoPointControl",
             "Invalid operation: two-point control can only be used"