Cascade of Two Continuous Stirred-Tank Reactors with Recycle

Initializing live version
Download to Desktop

Requires a Wolfram Notebook System

Interact on desktop, mobile and cloud with the free Wolfram Player or other Wolfram Language products.

Consider a cascade of two continuous stirred-tank reactors with recycle undergoing exothermic first-order 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 .


Contributed by: Housam Binous, Ahmed Bellagi, and Brian G. Higgins (January 2014)
Open content licensed under CC BY-NC-SA




[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/0009-2509(80)85088-3.

[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.

Feedback (field required)
Email (field required) Name
Occupation Organization
Note: Your message & contact information may be shared with the author of any specific Demonstration for which you give feedback.