WOLFRAM|DEMONSTRATIONS PROJECT

Multiple Steady States in Continuous Culture with Substrate Inhibition

​
inhibition constant
0.25
dilution rate
0.1
phase plane analysis
nutrient and biomass versus time
specific growth rate
The dynamic behavior of a chemostat can be described with the following system of ordinary differential equations:​
dX
dt
=(μ-D)X
​​
dS
dt
=-
μX
Y
+D(
S
in
-S)
,where
D
is the dilution rate in
-1
s
,
X
is the biomass concentration in g/l,
S
is the substrate concentration in g/l,
S
in
is the feed substrate concentration,
Y=-
dX
dS
is the yield, and
μ=
μ
max
S
K
S
+S+
2
S
K
I
is the specific growth rate (where
K
S
is the half-saturation constant,
K
I
is the inhibition constant, and
μ
max
is the maximum value of the specific growth rate, with
μ
max
=1
here).
​This Demonstration lets you plot the substrate (magenta curve) and biomass (blue curve) concentrations versus time. A phase plane analysis shows that there can be up to three steady states (SS) depending on the values of the dilution rate and inhibition coefficient. With the help of the specific growth rate curve, you can easily see when there are only three or one SS (in the latter case we have
μ<D
for all
S
). One SS always corresponds to the wash out case (
X=0
), while the other two, if present, correspond to the case where
D=μ
. Indeed, as can be seen in the snapshot, the red dots are stable SS and the green dot is unstable SS. You can drag the locator to change the initial conditions. The cyan curve is a parametric plot of substrate versus biomass.