,

,

,

,

,

,

where represents dimensionless time, is the dimensionless concentration of reactant , and and stand for the dimensionless temperatures of the reactor and cooling fluid, respectively. The following parameters are also dimensionless: , , and for the Damköhler number, the activation energy, and the heat transfer coefficient; , , and stand for the adiabatic temperature rise, the ratio of heat capacity of the reactor flow rate to the coolant flow rate, and the ratio of the heat capacity of the coolant volume to the reactor volume; and is the inlet coolant temperature. The system of equations is solved with

,

and initial values

.

In this high-dimensional nonlinear system, the entanglement of input and output between the two reactors leads to chaotic behavior. As , the Damköhler number of the first reactor, is increased, the trajectories change from oscillations leading to a steady state to periodic, then chaotic oscillations, and again to oscillations leading to a steady state.