From 7145ee94f9413e1b26c2ee2ba1727a9809a627e2 Mon Sep 17 00:00:00 2001 From: Harry Moffat Date: Tue, 17 Mar 2009 16:37:05 +0000 Subject: [PATCH] Switched the sample problem over to the vcs_equilibrate solver. This solver works on all of the temperatures, while the old solver failed after the 10th temperature or so. --- Cantera/python/examples/equilibrium/multiphase_plasma.py | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/Cantera/python/examples/equilibrium/multiphase_plasma.py b/Cantera/python/examples/equilibrium/multiphase_plasma.py index 61de96ceb..8e44fa426 100644 --- a/Cantera/python/examples/equilibrium/multiphase_plasma.py +++ b/Cantera/python/examples/equilibrium/multiphase_plasma.py @@ -26,14 +26,15 @@ f = open('equil_koh.csv','w') writeCSV(f,['T']+mix.speciesNames()) # loop over temperature -for n in range(500): - t = 350.0 + 10.0*n +for n in range(100): + t = 350.0 + 50.0*n print 'T = ',t mix.set(T= t, P = OneAtm, Moles="K:1.03, H2:2.12, O2:0.9") # set the mixture to a state of chemical equilibrium holding # temperature and pressure fixed - mix.equilibrate("TP",maxsteps=10000,loglevel=1) + # mix.equilibrate("TP",maxsteps=10000,loglevel=1) + mix.vcs_equilibrate("TP",printLvl=1,maxsteps=10000,loglevel=1) # write out the moles of each species writeCSV(f,[t]+ list(mix.speciesMoles()))