Batch Reactor with Multiple Reactions

The concentrations of five species are plotted as a function of dimensionless time, for three irreversible, elementary, liquid-phase reactions in an isothermal batch reactor. Use sliders to change the dimensionless rate constants for each reaction. The initial amounts of A and B in the reactor are 10 mol and 5 mol, respectively. Use buttons to plot just one concentration or all five together. The reactors are:
A + B → C
C → 2E 
2A → D

Download the CDF file to view the simulation offline using the free Wolfram CDF player.

Three liquid-phase reactions take place in an isothermal batch reactor:
A + B → C, with reaction rate r1 = k1CACB
→ 2E with reaction rate r2 = k2CC
2A → D with reaction rate r3 = k3C2A
where Ci is the concentration of component i, ri is the rate of reaction, and ki is the rate constant.

Mass balances on each component:
where t is time, NA are moles of A, and V is reactor volume. Initially 10 moles of A and 5 moles of B are fed into the reactor.