Homework 7.2: Dynamics of a chemostat (55 pts)


In this problem, we will investigate the dynamics of a chemostat. Here, a solution of bacteria in medium is held in a well-mixed, well-aerated vessel with a total solution volume \(V\). There is an inflow of sterile medium that has necessary nutrients for the bacteria to grow. To maintain the constant volume, there is an equal outflow of the bacterial medium. The outflowing medium contains depleted nutrients, living and dead bacteria, and waste.

Let \(Q\) be the volumetric inflow rate of fresh medium (and also the outflow rate). Let \(c\) be the concentration of living bacterial cells in the vessel and let \(n\) be the concentration of nutrient in the vessel. We denote by \(n_0\) the concentration of nutrients in the incoming fresh medium.

a) Explain why the following equation is a reasonable model for bacterial concentration.

\begin{align} \frac{\mathrm{d}c}{\mathrm{d}t} = \beta\,\frac{n/k}{1 + n/k}\,c - \frac{Q}{V}\,c. \end{align}

b) Taking the above equation as given, and defining by \(\nu\) the number of nutrient molecules needed to make one bacterium, derive an expression for \(\mathrm{d}n/\mathrm{d}t\).

c) Nondimensionalize your equations. You should end up with two dimensionless parameters. One is a dimensionless version of the inlet nutrient concentration and another is a dimensionless growth rate.

d) Analyze the circuit. Include at least the following in your analysis:

  • The number of fixed points for various parameter values.

  • Phase portraits with nullclines and fixed points.

Based on your analysis, propose an inlet flow rate that will ensure a stable bacterial population.