Dynamical behavior of a hepatitis B epidemic model and its NSFD scheme

DOI: https://doi.org/10.1007/s12190-024-02103-6

2024-05-14

Journal of Applied Mathematics and Computing

Abstract:Hepatitis is inflammation of the liver, and one of its types, hepatitis B, is a contagious infection. Using mathematical models, the nature of the spread of the Hepatitis B virus can be predicted. In the present paper, a hepatitis B epidemic model with a Beddington–DeAngelis type incidence rate and a constant vaccination rate is considered. Some dynamical properties of this model, such as non-negativity, boundedness character, the basic reproduction number , stability nature, and the bifurcation phenomenon, are investigated. By the Bendixson theorem, it is demonstrated that the disease-free equilibrium is globally asymptotically stable. It is shown that a transcritical bifurcation phenomenon occurs when . It is concluded that the endemic equilibrium is globally asymptotically stable when , by utilizing Dulac's criteria. Also, a discrete system of difference equations is obtained by constructing a non-standard finite difference (NSFD) scheme for the continuous model. It is shown that the solutions of this discrete system are dynamically consistent for all finite step sizes. The theoretical results obtained are also supported and visualized by numerical simulations. These simulations also demonstrate that the NSFD scheme produces much more efficient results than the Euler or RK4 schemes, as shown in the theoretical results obtained.

mathematics, applied

What problem does this paper attempt to address?

### Problems Addressed by the Paper

This paper primarily explores the dynamic behavior of the Hepatitis B infectious disease model and its non-standard finite difference (NSFD) scheme. Specifically, the paper aims to address the following issues:

1. **Model Establishment**:
- Propose a Hepatitis B infectious disease model that includes a Beddington-DeAngelis type contact rate and a constant vaccination rate.
- Study the non-negativity and boundedness of the model.

2. **Basic Reproduction Number (R₀)**:
- Determine the basic reproduction number R₀ of the model and study its impact on disease transmission.
- Prove the global asymptotic stability of the disease-free equilibrium (DFE) when R₀ < 1; prove the local asymptotic stability of the endemic equilibrium (EE) when R₀ > 1.

3. **Dynamic Characteristics Analysis**:
- Analyze the dynamic characteristics of the model, including non-negativity, boundedness, stability, and bifurcation phenomena.
- Prove the global asymptotic stability of the disease-free equilibrium using Bendixson's theorem and the global asymptotic stability of the endemic equilibrium using Dulac's criterion.

4. **Non-Standard Finite Difference (NSFD) Scheme**:
- Construct a non-standard finite difference (NSFD) scheme for the continuous model and verify the dynamic consistency of the discrete model.
- Compare the performance of the NSFD scheme with the Euler method and the Runge-Kutta method in numerical simulations, demonstrating the effectiveness of the NSFD scheme.

5. **Numerical Simulation**:
- Verify the theoretical results through numerical simulations and demonstrate the dynamic consistency of the NSFD scheme under different time steps.

In summary, this paper aims to study the transmission dynamics of Hepatitis B through mathematical modeling and numerical simulations, and proposes an effective non-standard finite difference scheme to improve the accuracy and stability of the model.