In [61]:

    

"""
Constant-pressure, adiabatic kinetics simulation.
"""

import sys
import numpy as np

import cantera as ct

gas = ct.Solution('cti_one.cti')
gas.TPX = 300, 101325, 'PLIGC:1.7, PLIGH:0.26, PLIGO:0.3'
r = ct.IdealGasReactor(gas)

sim = ct.ReactorNet([r])
time = 0.0
states = ct.SolutionArray(gas, extra=['t'])

print('%10s %10s %10s %14s' % ('t [s]','T [K]','P [Pa]','u [J/kg]'))
for n in range(100):
    time += 1.e-5
    sim.advance(time)
    states.append(r.thermo.state, t=time*1e3)
    print('%10.3e %10.3f %10.3f %14.6e' % (sim.time, r.T,
                                           r.thermo.P, r.thermo.u))


import matplotlib.pyplot as plt
plt.clf()
plt.subplot(2, 2, 1)
plt.plot(states.t, states.T)
plt.xlabel('Time (ms)')
plt.ylabel('Temperature (K)')
plt.subplot(2, 2, 2)
plt.plot(states.t, states.X[:,gas.species_index('PLIGC')])
plt.xlabel('Time (ms)')
plt.ylabel('PLIGC Mole Fraction')
plt.subplot(2, 2, 3)
plt.plot(states.t, states.X[:,gas.species_index('PLIGH')])
plt.xlabel('Time (ms)')
plt.ylabel('PLIGH Mole Fraction')
plt.subplot(2, 2, 4)
plt.plot(states.t, states.X[:,gas.species_index('PLIGO')])
plt.xlabel('Time (ms)')
plt.ylabel('PLIGO Mole Fraction')
plt.tight_layout()
plt.show()


    

    




     t [s]      T [K]     P [Pa]       u [J/kg]
 1.000e-05     33.584 101325.002  -7.509018e+03
 2.000e-05     32.890 101325.002  -7.509018e+03
 3.000e-05     32.570 101325.002  -7.509018e+03
 4.000e-05     32.355 101325.002  -7.509018e+03
 5.000e-05     32.189 101325.002  -7.509018e+03
 6.000e-05     32.051 101325.002  -7.509018e+03
 7.000e-05     31.932 101325.002  -7.509018e+03
 8.000e-05     31.826 101325.002  -7.509018e+03
 9.000e-05     31.730 101325.002  -7.509018e+03
 1.000e-04     31.643 101325.002  -7.509018e+03
 1.100e-04     31.563 101325.002  -7.509018e+03
 1.200e-04     31.488 101325.002  -7.509018e+03
 1.300e-04     31.419 101325.002  -7.509018e+03
 1.400e-04     31.355 101325.002  -7.509018e+03
 1.500e-04     31.294 101325.002  -7.509018e+03
 1.600e-04     31.237 101325.002  -7.509018e+03
 1.700e-04     31.183 101325.002  -7.509018e+03
 1.800e-04     31.132 101325.002  -7.509018e+03
 1.900e-04     31.084 101325.002  -7.509018e+03
 2.000e-04     31.038 101325.002  -7.509018e+03
 2.100e-04     30.995 101325.002  -7.509018e+03
 2.200e-04     30.953 101325.002  -7.509018e+03
 2.300e-04     30.914 101325.002  -7.509018e+03
 2.400e-04     30.876 101325.002  -7.509018e+03
 2.500e-04     30.840 101325.002  -7.509018e+03
 2.600e-04     30.805 101325.002  -7.509018e+03
 2.700e-04     30.772 101325.002  -7.509018e+03
 2.800e-04     30.740 101325.002  -7.509018e+03
 2.900e-04     30.710 101325.002  -7.509018e+03
 3.000e-04     30.680 101325.002  -7.509018e+03
 3.100e-04     30.652 101325.002  -7.509018e+03
 3.200e-04     30.625 101325.002  -7.509018e+03
 3.300e-04     30.599 101325.002  -7.509018e+03
 3.400e-04     30.574 101325.002  -7.509018e+03
 3.500e-04     30.549 101325.002  -7.509018e+03
 3.600e-04     30.526 101325.002  -7.509018e+03
 3.700e-04     30.503 101325.002  -7.509018e+03
 3.800e-04     30.481 101325.002  -7.509018e+03
 3.900e-04     30.460 101325.002  -7.509018e+03
 4.000e-04     30.439 101325.002  -7.509018e+03
 4.100e-04     30.420 101325.002  -7.509018e+03
 4.200e-04     30.400 101325.002  -7.509018e+03
 4.300e-04     30.382 101325.002  -7.509018e+03
 4.400e-04     30.363 101325.002  -7.509018e+03
 4.500e-04     30.346 101325.002  -7.509018e+03
 4.600e-04     30.329 101325.002  -7.509018e+03
 4.700e-04     30.312 101325.002  -7.509018e+03
 4.800e-04     30.296 101325.002  -7.509018e+03
 4.900e-04     30.280 101325.002  -7.509018e+03
 5.000e-04     30.265 101325.002  -7.509018e+03
 5.100e-04     30.250 101325.002  -7.509018e+03
 5.200e-04     30.236 101325.002  -7.509018e+03
 5.300e-04     30.222 101325.002  -7.509018e+03
 5.400e-04     30.208 101325.002  -7.509018e+03
 5.500e-04     30.195 101325.002  -7.509018e+03
 5.600e-04     30.182 101325.002  -7.509018e+03
 5.700e-04     30.169 101325.002  -7.509018e+03
 5.800e-04     30.157 101325.002  -7.509018e+03
 5.900e-04     30.145 101325.002  -7.509018e+03
 6.000e-04     30.133 101325.002  -7.509018e+03
 6.100e-04     30.122 101325.002  -7.509018e+03
 6.200e-04     30.111 101325.002  -7.509018e+03
 6.300e-04     30.100 101325.002  -7.509018e+03
 6.400e-04     30.089 101325.002  -7.509018e+03
 6.500e-04     30.078 101325.002  -7.509018e+03
 6.600e-04     30.068 101325.002  -7.509018e+03
 6.700e-04     30.058 101325.002  -7.509018e+03
 6.800e-04     30.049 101325.002  -7.509018e+03
 6.900e-04     30.039 101325.002  -7.509018e+03
 7.000e-04     30.030 101325.002  -7.509018e+03
 7.100e-04     30.020 101325.002  -7.509018e+03
 7.200e-04     30.012 101325.002  -7.509018e+03
 7.300e-04     30.003 101325.002  -7.509018e+03
 7.400e-04     29.994 101325.002  -7.509018e+03
 7.500e-04     29.986 101325.002  -7.509018e+03
 7.600e-04     29.977 101325.002  -7.509018e+03
 7.700e-04     29.969 101325.002  -7.509018e+03
 7.800e-04     29.961 101325.002  -7.509018e+03
 7.900e-04     29.954 101325.002  -7.509018e+03
 8.000e-04     29.946 101325.002  -7.509018e+03
 8.100e-04     29.939 101325.002  -7.509018e+03
 8.200e-04     29.931 101325.002  -7.509018e+03
 8.300e-04     29.924 101325.002  -7.509018e+03
 8.400e-04     29.917 101325.002  -7.509018e+03
 8.500e-04     29.910 101325.002  -7.509018e+03
 8.600e-04     29.903 101325.002  -7.509018e+03
 8.700e-04     29.897 101325.002  -7.509018e+03
 8.800e-04     29.890 101325.002  -7.509018e+03
 8.900e-04     29.884 101325.002  -7.509018e+03
 9.000e-04     29.877 101325.002  -7.509018e+03
 9.100e-04     29.871 101325.002  -7.509018e+03
 9.200e-04     29.865 101325.002  -7.509018e+03
 9.300e-04     29.859 101325.002  -7.509018e+03
 9.400e-04     29.853 101325.002  -7.509018e+03
 9.500e-04     29.847 101325.002  -7.509018e+03
 9.600e-04     29.842 101325.002  -7.509018e+03
 9.700e-04     29.836 101325.002  -7.509018e+03
 9.800e-04     29.830 101325.002  -7.509018e+03
 9.900e-04     29.825 101325.002  -7.509018e+03
 1.000e-03     29.820 101325.002  -7.509018e+03






    

In [35]:

    




    

In [38]:

    




    

In [ ]: