A simple model of nutrient recycling and dormancy in a chemostat: Mathematical analysis and a second-order nonstandard finite difference method

Abstract

A chemostat is an apparatus that sustains a homogeneous environment through continuous inflow and outflow. Presented is a chemostat model that characterizes the dynamics of dormancy-capable microorganisms. This model of coupled systems of nonlinear ordinary differential equations (ODEs) can apply to various types of organisms, such as different species of bacteria, archaea, algae, fungi, viruses, and protozoa. However, these species reside in different environments and rely on different sets of nutrients. Thus, the model adapts to each species’ limiting nutrient through nutrient recycling. This paper includes a complete stability analysis, supporting phase plane portraits, and accompanying bifurcation diagrams. The paper also proposes an advanced second-order, positivity-preserving, and elementary stable nonstandard finite difference method for solving the mathematical model. Series of numerical simulations are presented that support the theoretical results and explore different biological scenarios. The stability analysis reveals that (1) if the overall dilution, death, and conversion are less than the overall growth, both the dormant and active populations persist when introduced to a chemostat; and (2) if the overall dilution, death, and conversion are more than the overall growth, the microorganism population in its entirety will die off when introduced to a chemostat. Furthermore, the model study suggests that neither dormancy nor nutrient recycling provides substantial survival advantages in a basic chemostat when no threat to the active microorganism is present.