Cantera/ResidenceTimeReactor

from Cantera import *
from Cantera.Reactor import *

class ResidenceTimeReactor(Reactor):
    def __init__(self, g, tau):

        Reactor.__init__(self,g)

        self.upstream = Reservoir(g)
        self.downstream = Reservoir(g)

        self.inlet = Valve(self.upstream, self, Kv=100)
        self.outlet = MassFlowController(self, self.downstream)


if __name__=="__main__":

    g = importPhase('h2o2.cti')
    g.set(T=1200,P=OneAtm,X="H2:2.0, O2:1.0, H:0.000001")

    r = ResidenceTimeReactor(g,0.01)
    n = ReactorNet([r])

    t = 0
    dt = 1.0e-4
    for step in range(100):
        t = t + dt
        n.advance(t)

    print r