Consider a cascade of two continuous stirredtank reactors with recycle undergoing exothermic firstorder chemical reactions. The steady state is described by four nonlinear equations for the reaction conversions , , dimensionless temperature , , and parameter (Damköhler number). The dimensionless equations that describe the steady state are given by [2] (see also [1]) , , , . The derivation of these equations (which represent the species and energy balances for each tank) is given in [2]. A steady solution for is , , , and (see snapshot 4). The computational method used to generate is based on arc length continuation. The Demonstration shows the dependence of on the parameter α. The solution path has six turning points (indicated in green in snapshot 3) and for there are seven steady state solutions! (These are indicated in blue in snapshot 1.) Again, not all these steady state solutions are stable to small perturbations. The stability of the steady states is determined by the sign of the real part of the eigenvalues of the Jacobian matrix. The results are in quantitative agreement with the calculations in the references. Though not shown, one can readily plot the dependence of the remaining variables , , and on .
[1] M. Kubíček, H. Hofmann, V. Hlaváček, and J. Sinkule, "Multiplicity and Stability in a Sequence of Two Nonadiabatic Nonisothermal CSTR," Chemical Engineering Science, 35(4), 1980 pp. 987–996. doi: 10.1016/00092509(80)850883. [2] M. Kubíček and M. Marek, Computational Methods in Bifurcation Theory and Dissipative Structures, New York: Springer–Verlag, 1984. [3] B. G. Higgins and H. Binous, "A Simple Method for Tracking Turning Points in Parameter Space," Journal of Chemical Engineering of Japan, 43(12), 2010 pp. 1035–1042. doi: 10.1252/jcej.10we122.
