Global stability of a continuous bioreactor model under persistent variation of the dilution rate

In this work, the global stability of a continuous bioreactor model is studied, with the concentrations of biomass and substrate as state variables, a general non-monotonic function of substrate concentration for the specific growth rate, and constant inlet substrate concentration. Also, the dilution rate is time varying but bounded, thus leading to state convergence to a compact set instead of an equilibrium point. Based on the Lyapunov function theory with dead-zone modification, the convergence of the substrate and biomass concentrations is studied. The main contributions with respect to closely related studies are: i) The convergence regions of the substrate and biomass concentrations are determined as function of the variation region of the dilution rate (D) and the global convergence to these compact sets is proved, considering monotonic and non-monotonic growth functions separately; ii) several improvements are proposed in the stability analysis, including the definition of a new dead zone Lyapunov function and the properties of its gradient. These improvements allow proving convergence of substrate and biomass concentrations to their compact sets, while tackling the interwoven and nonlinear nature of the dynamics of biomass and substrate concentrations, the non-monotonic nature of the specific growth rate, and the time-varying nature of the dilution rate. The proposed modifications are a basis for further global stability analysis of bioreactor models exhibiting convergence to a compact set instead of an equilibrium point. Finally, the theoretical results are illustrated through numerical simulation, showing the convergence of the states under varying dilution rate.