Kinetics for multiphase reaction

[arthur gilles]

Dear all,

I try to simulate a plug flow reactor for methane pyrolysis. The reaction is a gas that is converted into H2 + C (solid).

I cannot find how to put a gas/solid mixture in a reactor?

Is it possible ?

My attempt :

gas2 = ct.Solution('gri30.xml');
carbon = ct.Solution('graphite.cti');
gas2.TPX=T_0, ct.one_atm, 'CH4:1'
mix_phases = [(gas2, 1.0), (carbon, 0.0)]
mix = ct.Mixture(mix_phases)
mix.T = T_0
mix.P = ct.one_atm
mass_flow_rate2 = u_0 * gas2.density * area
dz = length / n_steps
r_vol = area * dz


r2 = ct.Reactor(gas2, energy='off')==> here how to introduce carbon?
r2.volume = r_vol

Thanks