Best Operating Conditions for Biogas Production in Some Simple Anaerobic Digestion Models

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.

Operating diagrams for a three-tiered microbial food web in the chemostat

In this paper, we consider a three-tiered food web model in a chemostat, including chlorophenol, phenol, and hydrogen substrates and their degraders. The model takes into account the three substrate inflowing concentrations, as well as maintenance, that is, decay terms of the species. The operating diagrams give the asymptotic behavior of the model with respect to the four operating parameters, which are the dilution rate and the three inflowing concentrations of the substrates. These diagrams were obtained only numerically in the existing literature. Using the mathematical analysis of this model obtained in our previous studies, we construct the operating diagrams, by plotting the curves that separate their various regions. Hence, the regions of the operating diagrams are constructed analytically and there is no requirement for time-consuming algorithms to generate the plots, as in the numerical method. Moreover, our method reveals behaviors that have not been detected in the previous numerical studies. The growth functions are of Monod form with the inclusion of a product inhibition term. However, our method applies for a large class of growth functions. We construct operating diagrams with and without maintenance showing the role of maintenance on the stability of the system.

Special Issue: Mathematical Modeling and Control of Bioprocesses

This Special Issue (SI) of Processes on Mathematical Modeling and Control of Bioprocesses (MMCB) contains papers focusing, on the one hand, on mathematical modeling of biological processes at different scales ranging from microscopic to macroscopic levels and, on the other hand, on model-based estimation, optimization and control of these processes [...]