time += 5.0e-6
sim.advance(time)
tim[n] = time
data[n,0] = r.temperature()
data[n,1] = r.moleFraction('OH')
data[n,2] = r.moleFraction('H')
data[n,3] = r.moleFraction('CH4')
# sensitivity of OH to reaction 2
data[n,4] = sim.sensitivity('OH',2)
# sensitivity of OH to reaction 3
data[n,5] = sim.sensitivity('OH',3)
print '%10.3e %10.3f %10.3f %14.6e' % (sim.time(), r.temperature(),
r.pressure(), r.intEnergy_mass())
#sim.sensitivity("OH",0))
# plot the results if matplotlib is installed.
# see http://matplotlib.sourceforge.net to get it
args = sys.argv
if 1==1:
try:
from matplotlib.pylab import *
clf
subplot(2,2,1)
plot(tim,data[:,0])
xlabel('Time (s)');
ylabel('Temperature (K)');
subplot(2,2,2)
sim = ReactorNet([r])
time = 0.0
tim = zeros(100, 'd')
data = zeros([100, 5], 'd')
for n in range(100):
time += 1.e-5
sim.advance(time)
tim[n] = time
data[n, 0] = r.temperature()
data[n, 1] = r.moleFraction('OH')
data[n, 2] = r.moleFraction('H')
data[n, 3] = r.moleFraction('H2')
print '%10.3e %10.3f %10.3f %14.6e' % (sim.time(), r.temperature(),
r.pressure(), r.intEnergy_mass())
# plot the results if matplotlib is installed.
# see http://matplotlib.sourceforge.net to get it
args = sys.argv
if len(args) > 1 and args[1] == '-plot':
try:
from matplotlib.pylab import *
clf
subplot(2, 2, 1)
plot(tim, data[:, 0])
xlabel('Time (s)')
ylabel('Temperature (K)')
subplot(2, 2, 2)
plot(tim, data[:, 1])
xlabel('Time (s)')
time += 5.0e-6
sim.advance(time)
tim[n] = time
data[n,0] = r.temperature()
data[n,1] = r.moleFraction('OH')
data[n,2] = r.moleFraction('H')
data[n,3] = r.moleFraction('CH4')
# sensitivity of OH to reaction 2
data[n,4] = sim.sensitivity('OH',2)
# sensitivity of OH to reaction 3
data[n,5] = sim.sensitivity('OH',3)
print '%10.3e %10.3f %10.3f %14.6e %10.3f %10.3f' % (sim.time(), r.temperature(),
r.pressure(), r.intEnergy_mass(), data[n,4], data[n,5])
# plot the results if matplotlib is installed.
# see http://matplotlib.sourceforge.net to get it
args = sys.argv
if len(args) > 1 and args[1] == '-plot':
try:
from matplotlib.pylab import *
clf
subplot(2,2,1)
plot(tim,data[:,0])
xlabel('Time (s)');
ylabel('Temperature (K)');
subplot(2,2,2)
plot(tim,data[:,1])
c = (4.0-3.0*R)/(4.0*(1.0-R))
e = c/phi
d = 3.7*e
comp = 'CH4:0.5, H2:%(b)f, O2:%(e)f, N2:%(d)f'% vars()
print '#'+comp
gas = GRI30()
gas.set(T = T0, P = OneAtm, X = comp)
r = Reactor(gas)
env = Reservoir(Air())
w = Wall(r,env)
w.set(K = 0) # set expansion parameter. dV/dt = KA(P_1 - P_2)
w.set(A = 1.0)
sim = ReactorNet([r])
time = 0.0
#Told = r.temperature()
#print '%10.3e %10.3f %10.3f %14.6e' % (sim.time(), r.temperature(), r.pressure(), r.intEnergy_mass())
for n in range(36):
time = (n+1)*0.0005
sim.advance(time)
print '%10.3e %10.3f %10.3f %14.6e' % (sim.time(), r.temperature(), r.pressure(), r.intEnergy_mass())
#while time <= tFinal:
# time = sim.step(tFinal)
# print '%10.3e %10.3f %10.3f %14.6e' % (sim.time(), r.temperature(),
# r.pressure(), r.intEnergy_mass())#print MC_Gemisch(3,1)
time += 5.0e-6
sim.advance(time)
tim[n] = time
data[n, 0] = r.temperature()
data[n, 1] = r.moleFraction('OH')
data[n, 2] = r.moleFraction('H')
data[n, 3] = r.moleFraction('CH4')
# sensitivity of OH to reaction 2
data[n, 4] = sim.sensitivity('OH', 2)
# sensitivity of OH to reaction 3
data[n, 5] = sim.sensitivity('OH', 3)
print '%10.3e %10.3f %10.3f %14.6e %10.3f %10.3f' % (
sim.time(), r.temperature(), r.pressure(), r.intEnergy_mass(),
data[n, 4], data[n, 5])
# plot the results if matplotlib is installed.
# see http://matplotlib.sourceforge.net to get it
args = sys.argv
if len(args) > 1 and args[1] == '-plot':
try:
from matplotlib.pylab import *
clf
subplot(2, 2, 1)
plot(tim, data[:, 0])
xlabel('Time (s)')
ylabel('Temperature (K)')
subplot(2, 2, 2)
plot(tim, data[:, 1])
