Stability and stabilization for models of chemostats with multiple limiting substrates.

We study chemostat models in which multiple species compete for two or more limiting nutrients. First, we consider the case where the nutrient flow and species removal rates and input nutrient concentrations are all given as positive constants. In that case, we use Brouwer degree theory to give conditions guaranteeing that the models admit globally asymptotically stable componentwise positive equilibrium points, from all componentwise positive initial states. Then we use the results to develop stabilization theory for a class of controlled chemostats with two or more limiting nutrients. For cases where the dilution rate and input nutrient concentrations can be selected as controls, we prove that many different componentwise positive equilibria can be made globally asymptotically stable. This extends the existing control results for chemostats with one limiting nutrient. We demonstrate our methods in simulations.

On the stability of periodic solutions in the perturbed chemostat.

We study the chemostat model for one species competing for one nutrient using a Lyapunov-type analysis. We design the dilution rate function so that all solutions of the chemostat converge to a prescribed periodic solution. In terms of chemostat biology, this means that no matter what positive initial levels for the species concentration and nutrient are selected, the long-term species concentration and substrate levels closely approximate a prescribed oscillatory behavior. This is significant because it reproduces the realistic ecological situation where the species and substrate concentrations oscillate. We show that the stability is maintained when the model is augmented by additional species that are being driven to extinction. We also give an input-to-state stability result for the chemostat-tracking equations for cases where there are small perturbations acting on the dilution rate and initial concentration. This means that the long-term species concentration and substrate behavior enjoys a highly desirable robustness property, since it continues to approximate the prescribed oscillation up to a small error when there are small unexpected changes in the dilution rate function.