Phase-Plane Analysis of a Biochemical Reactor

A biochemical reactor can show complex dynamic behavior, including multiple steady states. This Demonstration plots the phase-plane diagram and computes the stability of the steady states. You can choose either the Monod model or the substrate inhibition model for the growth rate expression. Parameters of the specific growth rate expression are taken from [1]. You can select numerical values for the yield, dilution rate, and substrate concentration in the feed stream.
The Monod model has two steady states, which are indicated in the phase-plane diagram by the green and blue dots. As you can see in Snapshot 1, the green and the blue steady states are a saddle point and a stable node, respectively. The eigenvalues of the Jacobian matrix confirm this observation (see Snapshot 3, where and are the two eigenvalues for the steady states). Each steady state () is defined by its two coordinates ( and ). If both eigenvalues are negative, then the steady state is a stable node. If the eigenvalues have opposite signs, then the steady state is a saddle point and is unstable. Finally, if both eigenvalues are positive, the steady state is an unstable node. The green steady state corresponds to washout (i.e., no biomass is left in the bioreactor).
When the substrate inhibition model is used, up to three steady states are obtained. These steady states are indicated by the red, green, and blue dots. The stability analysis (i.e., determined by the signs of the Jacobian matrix eigenvalues) indicates that the steady states given by the blue and red dots are stable, while the one represented by the green dot is a saddle point (see Snapshot 2). Again, a simple eigenvalue computation for the Jacobian matrix confirms this observation (see Snapshot 4, where and are the two eigenvalues). Here, the red steady state corresponds to the washout.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


In 1942, Monod proposed the following form of the specific growth coefficient:
The specific growth coefficient with the substrate inhibition is given by:
The biochemical reactor is governed by the coupled equations:
where is the biomass concentration, is the substrate concentration, is the dilution rate, is the yield, is the feed substrate concentration, and is the specific growth coefficient.
Here, , , are taken from [1].
[1] B. Wayne Bequette, Process Dynamics, Modeling, Analysis, and Simulation, Upper Saddle River, NJ: Prentice Hall, 1998.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2014 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+