From bd2c97399ec4d6342fc9679c93d654d9122d2874 Mon Sep 17 00:00:00 2001 From: Harry Moffat Date: Wed, 3 May 2006 18:12:58 +0000 Subject: [PATCH] Changed the initial condition slightly. This was to avoid a problem with convergence that occurred in the chemical equilibrium solver. --- examples/cxx/equil_example1.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/examples/cxx/equil_example1.cpp b/examples/cxx/equil_example1.cpp index 928bf87fa..8145c2252 100755 --- a/examples/cxx/equil_example1.cpp +++ b/examples/cxx/equil_example1.cpp @@ -85,7 +85,7 @@ int equil_example1(int job) { for (int i = 0; i < ntemps; i++) { temp = tlow + dt*i; if (temp > gas.maxTemp()) break; - gas.setState_TPX(temp, pres, "SIH4:0.01, H2:1.0"); + gas.setState_TPX(temp, pres, "SIH4:0.01, H2:0.99"); equilibrate(gas,"TP"); output(0,i) = temp;