From 94a4439bd46216b911ae5d23353dfd7ec1d9d083 Mon Sep 17 00:00:00 2001 From: Ray Speth Date: Mon, 15 May 2017 18:12:26 -0400 Subject: [PATCH] [1D] Impose upper temperature bound based on thermo data This fixes some platform-specific test failures when compiling in debug mode, where estimated temperatures during the Newton iterations were too high, resulting in non-finite reaction rates. --- src/oneD/StFlow.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/oneD/StFlow.cpp b/src/oneD/StFlow.cpp index 7ba697197..d667822c7 100644 --- a/src/oneD/StFlow.cpp +++ b/src/oneD/StFlow.cpp @@ -66,7 +66,7 @@ StFlow::StFlow(IdealGasPhase* ph, size_t nsp, size_t points) : //-------------- default solution bounds -------------------- setBounds(0, -1e20, 1e20); // no bounds on u setBounds(1, -1e20, 1e20); // V - setBounds(2, 200.0, 1e9); // temperature bounds + setBounds(2, 200.0, 2*m_thermo->maxTemp()); // temperature bounds setBounds(3, -1e20, 1e20); // lambda should be negative // mass fraction bounds