Changed the initial condition slightly. This was to avoid a

problem with convergence that occurred in the chemical
equilibrium solver.
This commit is contained in:
Harry Moffat 2006-05-03 18:12:58 +00:00
parent 90a83569d0
commit bd2c97399e

View file

@ -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;