# American Institute of Mathematical Sciences

March  2020, 13(3): 957-974. doi: 10.3934/dcdss.2020056

## A fractional order HBV model with hospitalization

 1 Department of Mathematics, University of Peshawar, Khyber Pakhtunkhwa, 25120, Pakistan 2 Department of Mathematics, City University of Science and Information Technology, Peshawar, Khyber Pakhtunkhwa, 25000, Pakistan 3 Department of Mathematics, Abbottabad University of Science and Technology, Abbottabad, Khyber Pakhtunkhwa, Pakistan

* Corresponding author: altafdir@gmail.com

Received  April 2018 Revised  June 2018 Published  March 2019

Hepatitis B is a viral infection that can cause both acute and chronic disease and mainly attacks the liver. The present paper describes the dynamics of HBV with hospitalization. Due to the fatal nature of this disease, it is necessary to formulate a new mathematical model in order to reduce the burden of HBV. Therefore, we formulate a new HBV model with fractional order derivative. The fractional order model is formulated in Caputo sense. Two equilibria for the model exist: the disease-free and the endemic equilibriums. It is shown, that the disease-free equilibrium is both locally and globally asymptotically stable if $\mathcal{R}_0<1$ for any $\alpha\in(0,1)$. The sensitivity analysis of the model parameters are calculated and their results are depicted. The numerical results for the stability of the endemic equilibrium are presented. The complex dynamics of the disease can be best described by using the fractional derivative and this is illustrated through many graphical results.

Citation: Saif Ullah, Muhammad Altaf Khan, Muhammad Farooq, Taza Gul, Fawad Hussain. A fractional order HBV model with hospitalization. Discrete & Continuous Dynamical Systems - S, 2020, 13 (3) : 957-974. doi: 10.3934/dcdss.2020056
##### References:

show all references

##### References:
The effect of $\delta$ and $h_1$ on $\mathcal{R}_0$
Contour plot of $\delta$ and $h_1$
The effect of $\delta$ and $h_2$ on $\mathcal{R}_0$
Contour plot of $\delta$ and $h_2$
The effect of $\mu$ and $h_2$ on $\mathcal{R}_0$
Contour plot of $\mu$ and $h_2$
The effect of $\mu$ and $h_1$ on $\mathcal{R}_0$
Contour plot of $\mu$ and $h_1$
The plot shows the susceptible individuals when $\mathcal{R}_0 = 1.0248>1$ for different values of $\alpha$.
The plot shows the exposed individuals when $\mathcal{R}_0 = 1.0248>1$ for different values of $\alpha$.
The plot shows the acute individuals when $\mathcal{R}_0 = 1.0248>1$ for different values of $\alpha$.
The plot shows the carrier individuals when $\mathcal{R}_0 = 1.0248>1$ for different values of $\alpha$.
The plot shows the hospitalized individuals when $\mathcal{R}_0 = 1.0248>1$ for different values of $\alpha$.
The plot shows the recovered individuals when $\mathcal{R}_0 = 1.0248>1$ for different values of $\alpha$.
The plot shows the total number of infected individuals when $\mathcal{R}_0 = 1.0248>1$ for different values of $\alpha$
The plot shows the total number of infected individuals when $\mathcal{R}_0 = 1.0248>1$ for different values of $\alpha$ and $h_1$
The plot shows the total number of infected individuals when $\mathcal{R}_0 = 1.0248>1$ for different values of $\alpha$ and $h_2$.
Values of parameters used for numerical simulations
 parameters description of parameter Values $b$ Birth rate 0.4 $d$ Natural death rate 0.01 $h_1$ The acute individuals to be hospitalized 0.01 $h_2$ Flow rate from carrier class to the hospitalized class 0.01 $\beta$ The transmission coefficient 0.0002 $\delta$ Rate of flow from exposed to carrier 0.01 $d_A$ mortality rate due to acute infection 0.001 $d_C$ carrier individuals death rate 0.002 $\gamma$ the rate by which acute individuals move to carries class 0.01 $\xi$ The rate of recovery 0.02 $\psi$ Un-immunized children born to carrier mothers 0.2 $\mu$ Carriers infectiousness related to acute infection 0.2
 parameters description of parameter Values $b$ Birth rate 0.4 $d$ Natural death rate 0.01 $h_1$ The acute individuals to be hospitalized 0.01 $h_2$ Flow rate from carrier class to the hospitalized class 0.01 $\beta$ The transmission coefficient 0.0002 $\delta$ Rate of flow from exposed to carrier 0.01 $d_A$ mortality rate due to acute infection 0.001 $d_C$ carrier individuals death rate 0.002 $\gamma$ the rate by which acute individuals move to carries class 0.01 $\xi$ The rate of recovery 0.02 $\psi$ Un-immunized children born to carrier mothers 0.2 $\mu$ Carriers infectiousness related to acute infection 0.2
Sensitivity indices of $\mathcal{R}_0$ with respect to the model parameters
 Parameter Sensitivity index $\beta$ +0.7165 $\delta$ +0.5000 $\gamma$ +0.3680 $h_1$ -0.0455 $h_2$ -0.3739 $\mu$ +0.0161 $\psi$ +0.8064 $d_A$ -0.0454 $d_C$ -0.0748
 Parameter Sensitivity index $\beta$ +0.7165 $\delta$ +0.5000 $\gamma$ +0.3680 $h_1$ -0.0455 $h_2$ -0.3739 $\mu$ +0.0161 $\psi$ +0.8064 $d_A$ -0.0454 $d_C$ -0.0748
 [1] Reza Chaharpashlou, Abdon Atangana, Reza Saadati. On the fuzzy stability results for fractional stochastic Volterra integral equation. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020432 [2] Christopher S. Goodrich, Benjamin Lyons, Mihaela T. Velcsov. Analytical and numerical monotonicity results for discrete fractional sequential differences with negative lower bound. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2020269 [3] A. M. Elaiw, N. H. AlShamrani, A. Abdel-Aty, H. Dutta. Stability analysis of a general HIV dynamics model with multi-stages of infected cells and two routes of infection. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020441 [4] Omid Nikan, Seyedeh Mahboubeh Molavi-Arabshai, Hossein Jafari. Numerical simulation of the nonlinear fractional regularized long-wave model arising in ion acoustic plasma waves. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020466 [5] Laurence Cherfils, Stefania Gatti, Alain Miranville, Rémy Guillevin. Analysis of a model for tumor growth and lactate exchanges in a glioma. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020457 [6] Yining Cao, Chuck Jia, Roger Temam, Joseph Tribbia. Mathematical analysis of a cloud resolving model including the ice microphysics. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 131-167. doi: 10.3934/dcds.2020219 [7] Hai-Feng Huo, Shi-Ke Hu, Hong Xiang. Traveling wave solution for a diffusion SEIR epidemic model with self-protection and treatment. Electronic Research Archive, , () : -. doi: 10.3934/era.2020118 [8] Chao Xing, Jiaojiao Pan, Hong Luo. Stability and dynamic transition of a toxin-producing phytoplankton-zooplankton model with additional food. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2020275 [9] Leilei Wei, Yinnian He. A fully discrete local discontinuous Galerkin method with the generalized numerical flux to solve the tempered fractional reaction-diffusion equation. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020319 [10] Claudianor O. Alves, Rodrigo C. M. Nemer, Sergio H. Monari Soares. The use of the Morse theory to estimate the number of nontrivial solutions of a nonlinear Schrödinger equation with a magnetic field. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2020276 [11] Min Chen, Olivier Goubet, Shenghao Li. Mathematical analysis of bump to bucket problem. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5567-5580. doi: 10.3934/cpaa.2020251 [12] Marco Ghimenti, Anna Maria Micheletti. Compactness results for linearly perturbed Yamabe problem on manifolds with boundary. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020453 [13] Scipio Cuccagna, Masaya Maeda. A survey on asymptotic stability of ground states of nonlinear Schrödinger equations II. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020450 [14] Qianqian Han, Xiao-Song Yang. Qualitative analysis of a generalized Nosé-Hoover oscillator. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020346 [15] Vieri Benci, Sunra Mosconi, Marco Squassina. Preface: Applications of mathematical analysis to problems in theoretical physics. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020446 [16] Zexuan Liu, Zhiyuan Sun, Jerry Zhijian Yang. A numerical study of superconvergence of the discontinuous Galerkin method by patch reconstruction. Electronic Research Archive, 2020, 28 (4) : 1487-1501. doi: 10.3934/era.2020078 [17] Vieri Benci, Marco Cococcioni. The algorithmic numbers in non-archimedean numerical computing environments. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020449 [18] Zhilei Liang, Jiangyu Shuai. Existence of strong solution for the Cauchy problem of fully compressible Navier-Stokes equations in two dimensions. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020348 [19] Thabet Abdeljawad, Mohammad Esmael Samei. Applying quantum calculus for the existence of solution of $q$-integro-differential equations with three criteria. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020440 [20] Xavier Carvajal, Liliana Esquivel, Raphael Santos. On local well-posedness and ill-posedness results for a coupled system of mkdv type equations. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020382

2019 Impact Factor: 1.233