Mathematical analysis of a two-tiered microbial food-web model for the anaerobic digestion process

In this research paper, we presented a four-dimensional mathematical system modeling the anaerobic mineralization of phenol in a two-step microbial food-web. The inflowing concentrations of the hydrogen and the phenol are considered in our model. We considered the case of general class of nonlinear growth kinetics, instead of Monod kinetics. Due to some conservative relations, the proposed model was reduced to a two-dimensional system. The stability of the steady states was carried out. Based on the species growth rates and the three main operating parameters of the model, represented by the dilution rate and input concentrations of the phenol and the hydrogen, we showed that the system can have up to four steady states. We gave the necessary and sufficient conditions ensuring the existence and the stability for each feasible equilibrium state. We showed that in specific cases, the positive steady state exists and is stable. We gave numerical simulations validating the obtained results.

Exploring HIV Dynamics and an Optimal Control Strategy

In this paper, we propose a six-dimensional nonlinear system of differential equations for the human immunodeficiency virus (HIV) including the B-cell functions with a general nonlinear incidence rate. The compartment of infected cells was subdivided into three classes representing the latently infected cells, the short-lived productively infected cells, and the long-lived productively infected cells. The basic reproduction number was established, and the local and global stability of the equilibria of the model were studied. A sensitivity analysis with respect to the model parameters was undertaken. Based on this study, an optimal strategy is proposed to decrease the number of infected cells. Finally, some numerical simulations are presented to illustrate the theoretical findings.

A mathematical investigation of an "SVEIR" epidemic model for the measles transmission

A generalized "SVEIR" epidemic model with general nonlinear incidence rate has been proposed as a candidate model for measles virus dynamics. The basic reproduction number $ \mathcal{R} $, an important epidemiologic index, was calculated using the next generation matrix method. The existence and uniqueness of the steady states, namely, disease-free equilibrium ($ \mathcal{E}_0 $) and endemic equilibrium ($ \mathcal{E}_1 $) was studied. Therefore, the local and global stability analysis are carried out. It is proved that $ \mathcal{E}_0 $ is locally asymptotically stable once $ \mathcal{R} $ is less than. However, if $ \mathcal{R} > 1 $ then $ \mathcal{E}_0 $ is unstable. We proved also that $ \mathcal{E}_1 $ is locally asymptotically stable once $ \mathcal{R} > 1 $. The global stability of both equilibrium $ \mathcal{E}_0 $ and $ \mathcal{E}_1 $ is discussed where we proved that $ \mathcal{E}_0 $ is globally asymptotically stable once $ \mathcal{R}\leq 1 $, and $ \mathcal{E}_1 $ is globally asymptotically stable once $ \mathcal{R} > 1 $. The sensitivity analysis of the basic reproduction number $ \mathcal{R} $ with respect to the model parameters is carried out. In a second step, a vaccination strategy related to this model will be considered to optimise the infected and exposed individuals. We formulated a nonlinear optimal control problem and the existence, uniqueness and the characterisation of the optimal solution was discussed. An algorithm inspired from the Gauss-Seidel method was used to resolve the optimal control problem. Some numerical tests was given confirming the obtained theoretical results.

Mathematical analysis and optimal control for Chikungunya virus with two routes of infection with nonlinear incidence rate

Chikungunya fever, caused by Chikungunya virus (CHIKV) and transmitted to humans by infected Aedes mosquitoes, has posed a global threat in several countries. In this paper, we investigated a modified within-host Chikungunya virus (CHIKV) infection model with antibodies where two routes of infection are considered. In a first step, the basic reproduction number [Formula: see text] was calculated and the local and global stability analysis of the steady states is carried out using the local linearization and the Lyapunov method. It is proven that the CHIKV-free steady-state [Formula: see text] is globally asymptotically stable when [Formula: see text], and the infected steady-state [Formula: see text] is globally asymptotically stable when [Formula: see text]. In a second step, we applied an optimal strategy in order to optimize the infected compartment and to maximize the uninfected one. For this, we formulated a nonlinear optimal control problem. Existence of the optimal solution was discussed and characterized using some adjoint variables. Thus, an algorithm based on competitive Gauss–Seidel-like implicit difference method was applied in order to resolve the optimality system. The theoretical results are confirmed by some numerical simulations.

Mutual inhibition in presence of a virus in continuous culture

In this paper, we consider two species competing for a limiting substrate such that each species impedes the growth of the other one (Mutual inhibition) in presence of a virus inhibiting one bacterial species. A system of ordinary differential equations is proposed as a mathematical model for this competition. A detailed local qualitative analysis of the system is carried out. We proved that for a general nonlinear growth rates, the Competitive Exclusion Principle still valid, that at least one species goes extinct. For some cases where we have two locally stable equilibrium points, initial species concentrations are important in determining which is the winning species. Obtained results were confirmed by some numerical simulations using Matlab software.

Interspecific density-dependent model of predator–prey relationship in the chemostat

The objective of this study is to analyze a model of competition for one resource in the chemostat with general interspecific density-dependent growth rates, taking into account the predator–prey relationship. This relationship is characterized by the fact that the prey species promotes the growth of the predator species which in turn inhibits the growth of the first species. The model is a three-dimensional system of ordinary differential equations. With the same dilution rates, the model can be reduced to a planar system where the two models have the same local and even global behavior. The existence and stability conditions of all steady states of the reduced model in the plane are determined according to the operating parameters. Using the nullcline method, we present a geometric characterization of the existence and stability of all equilibria showing the multiplicity of coexistence steady states. The bifurcation diagrams illustrate that the steady states can appear or disappear only through saddle-node or transcritical bifurcations. Moreover, the operating diagrams describe the asymptotic behavior of this system by varying the control parameters and show the effect of the inhibition of predation on the emergence of the bistability region and the reduction until the disappearance of the coexistence region by increasing this inhibition parameter.

Modelling and optimal control for Chikungunya disease

A generalized model of intra-host CHIKV infection with two routes of infection has been proposed. In a first step, the basic reproduction number [Formula: see text] was calculated using the next-generation matrix method and the local and global stability analyses of the steady states are carried out using the Lyapunov method. It is proven that the CHIKV-free steady state [Formula: see text] is globally asymptotically stable when [Formula: see text] and the infected steady state [Formula: see text] is globally asymptotically stable when [Formula: see text]. In a second step, we applied an optimal strategy via the antibodies' flow rate in order to optimize infected compartment and to maximize the uninfected one. For this, we formulated a nonlinear optimal control problem. Existence of the optimal solution was discussed and characterized using an adjoint variables. Thus, an algorithm based on competitive Gauss-Seidel-like implicit difference method was applied in order to resolve the optimality system. The theoretical results are confirmed by some numerical simulations.