Study of Fractional Order Dynamical System of Viral Infection Disease under Piecewise Derivative

Kamal Shah,Hafsa Naz,Thabet Abdeljawadand Bahaaeldin Abdalla

1Department of Mathematics and Sciences,Prince Sultan University,P.O.Box 66833,Riyadh,11586,Saudi Arabia

2Department of Mathematics,University of Malakand,Chakdara Dir(L),Khyber Pakhtunkhwa,18000,Pakistan

3Department of Medical Research,China Medical University,Taichung,40402,Taiwan

ABSTRACT This research aims to understand the fractional order dynamics of the deadly Nipah virus(NiV)disease.We focus on using piecewise derivatives in the context of classical and singular kernels of power operators in the Caputo sense to investigate the crossover behavior of the considered dynamical system. We establish some qualitative results about the existence and uniqueness of the solution to the proposed problem.By utilizing the Newtonian polynomials interpolation technique, we recall a powerful algorithm to interpret the numerical findings for the aforesaid model.Here,we remark that the said viral infection is caused by an RNA type virus which can transmit from animals and also from an infected person to person. Fruits bats which are also known as flying foxes are one of the sources of transmission of NiV disease.Here in this work,we investigate its transmission mechanism through some new concepts of fractional calculus for further analysis and prediction.We present the approximate results for different compartments using different fractional orders.By using the piecewise derivative concept,we detect the crossover or multi-steps behavior in the transmission dynamics of the mentioned disease.Therefore,the considered form of the derivative is used to deal with problems exhibiting crossover behaviors.

KEYWORDS NiV disease;fractional calculus;piecewise derivative;qualitative results;newton polynomial;RNA virus

1 Introduction

From ancient times,viral diseases have remained a great threat to human as well as animal life on this globe. The said infections have caused millions of death in the past.Sometimes these infections have given birth to various terrible outbreaks in which millions of people lost their lives.For instance,viral infections like swine flu,MERS,bird flu,Nipah,and Henipa have triggered global public health disasters in the last many decades [1]. These infections have reappeared after some period of time and spread dramatically in most parts of the world. During their transmission, millions of people have died(see[2,3]).NiV disease is a highly contagious viral disease and can transmit among animals and humans. NiV infection causes encephalitis (brain swelling), which can vary from severe chronic sickness to even death.The said virus was identified for the first time in the town of Malaysia Nipah in 1998.Therefore,this virus has been named NiV.It was initially identified in a patient in the said town of Malaysia. From September 1998 to May 1999, it was initially identified in a massive epidemic of 276 cases in Peninsular Malaysia and Singapore[4].Migratory fruit bats have been reported as natural reservoirs of this disease.Bats that are infected transmit viruses through their wastes and mucus,such as spitting, urinating, and reproductive cells. Hence, they are asymptomatic carriers of the disease[5,6].The NIV is extremely infectious in pigs and transmitted through coughing.Direct contact with infected pigs has been found to be the primary mechanism of transfer in individuals when it was initially detected in a massive epidemic in Malaysia in 1999 (see [7]). More than eighty percent of people during the 1998 to 1999 outbreaks were belonging to the pigs’farmers’community or were exposed to pigs(see[8]).

The awareness of NiV transmission has made tremendous progress over the last two decades,especially in ways to investigate and study human Nipah microbial infections. Also, WHO declared it as one of the most dangerous and highly contagious diseases in 2015. Diagnostic equipment has been invented in the context of the Malaysian pandemic which is helping us to detect cases much more easily.Two NiV epidemics arose in India and Bangladesh in 2001(see detail[9]).After outbreaks in the Philippines in 2014 as well as in Kerala India in 2018,the geographic range of human infections of the NiV virus has continued to expand.NiV is a disease that can be passed from person to person,but no cases have been reported in Malaysia or Singapore in this context.Incidents of humans to humans and animals to humans have been reported in Bangladesh,the Philippines,and India.More research on the dynamics of NiV transmission, however, is required. Clear symptoms of NiV illness include fever,memory problems,migraines,dysentery,asthma,chest tightness,extreme weakness,and tremors[10]. Approximately 10% of the total number of Nipah-infected people transmit the virus to others.The most efficient technique to prevent the infection from spreading is quarantine. The majority of viral infections behave in an unanticipated manner. However, there are a number of elements that contribute to a viral disease’s becoming a pandemic[11].In the case of the said disease,the primary reservoir,genetic mutations,infectious transmitters,speed of transmission,and overflow from animal sources are some of the crucial components that should be addressed with it and comprehended in order to control the next pandemic outbreak [12]. Recently some updated work has been published on NiV disease. For example, authors [13] have analyzed the immunology of the virus. Similarly on new analysis and information about the said disease, authors have published good results in [14,15]recently.

Here, we demonstrate that mathematical models are essential tools for understanding how to effectively investigate the dynamics of transmissions and treatment of these infections. We refer to the importance of the mathematical model of the book for the reader in [16]. The literature on NiV contains a wide diversity of research. A few researchers have worked on the medical trial of NIV disease (see [17]). The authors have explained the characteristic symptoms of the disease encephalitis patients among Malaysian pigs and its psychological effect on human life (we refer to[18]). Further, some researchers worked on controlling of the disease in the community (see [19]).Some authors have worked on clinical trials.For instance,authors[20]have investigated the clinical features of NiV encephalitis among pig farmers in Malaysia.In the same line,authors[21]have studied differences in epidemiologic and clinical features of Nipah virus encephalitis between the Malaysian and Bangladesh outbreaks.Additionally,authors[22]have analyzed the clinical presentation of NiV disease in Bangladesh. Biswas established a mathematical model consisting of ordinary differential equations for the transmission of the NiV disease.The said model has been solved numerically by the said author.Further,some authors investigated the nature of the infection using an optimal control method[23].

So far,ordinary derivatives have been used as powerful tools to study the dynamical behavior of real-world problems.Since most of the dynamical problems suffer from abrupt or sudden changes in their state of rest or motions and hence show multi-behavior which is known as crossover behavior.It has been investigated that the said behavior cannot be explained by the concept of traditional calculus properly.Also,ordinary fractional order derivatives involving exponential,Mittag-Leffler,and power law kernel do not explain these terms significantly. It is because many real-world models exhibit crossover phenomena that are difficult to describe using traditional fractional order derivatives.For instance,earthquakes,pendulum motion,and contemporary economic fluctuations in less developed nations,experience rapid shifts in their state of dynamics.Instead of utilizing continuous derivatives,piecewise derivatives have been proved an ideal way to show this crossover phenomenon. In this context,the mentioned derivative’s key characteristics have recently been examined in[24].The authors have defined the classical and global piecewise derivatives,as well as several examples.Here it should be kept in mind that traditional fractional order derivative has been increasingly used in mathematical models of various infectious diseases like COVID-19. A mathematical model of COVID-19 using fractional order derivative has been investigated in[25].Moreover,some mathematical models under the fractional order derivatives have been investigated in [26]. In the same line, the author has also studied a model of COVID-19 using the concept of fractional calculus. Authors [27] have studied attractors of chaotic dynamical systems with fractal-fractional operators.In the same fashion,authors[28] have established a fractional order model for the investigation of COVID-19. Furthermore,authors[29]have investigated a fractional order susceptible,infected,and recovered epidemic model of childhood disease numerically. Authors [30] have developed and analyzed a mathematical model for the dynamics of COVID-19 with re-infection. In the same line, authors [31] have considered the fractional order predator-prey model with the harvesting rate.Keeping the importance of fractional calculus, authors [32] have investable transmission dynamics of Nipah virus disease under Caputo fractional derivative.

2 Formulation of NiV Model

Inspired by the importance of fractional calculus,we considered the integer order model of NiV which has been studied under some new concept of fractional order derivative. Authors [33] have presented a novel NiV disease model under traditional order derivative as presented in Fig.1 and described as

where various compartments, parameters and variables involved in the model are explained in the Table 1. Here, the total population of the community at timetis N, which we divide into five compartments including susceptible,exposed,infected,treated,and recovered classes respectively.Let the susceptible class S population be recruited at the rateΠ, and diminished by the class who gets an infection and hence follows contact with an infected class at a rateDue to the natural death rate,the susceptible population further decays by the amountμS.In the same fashion,the class E is generated by the susceptible class when it comes in contact with infected people at the rateThe exposed population is further reduced by the infection rateΦand a natural death rateμ. Due to infection rateΦin the exposed class and individuals who get an infection in the treating class at a rateη,the population of the infected class is generated asΦE+ηT.The population of the infected class is reduced by the natural death rateμ, treatment rateΥ, and natural recovery rates in infected and treatment groups at ratesε1andρdeath rate of infection respectively.Furthermore,the treated class is generated from the infected class with rateΥand reduced by ratesa,η,θρ,ε2,μ,respectively.The last recovered class is generated at the ratea,ε1, andε2respectively from the treated class, and the population is reduced by natural death rateμ.Hence briefly,we describe the transmission mechanism of the proposed model (1) in the flow chart 1. Authors have investigated NiV disease in [34,35] by using the concept provided in[36].

Figure 1:Flow chart of the model(1)

Table 1: The nomenclature and their description of NiV model

Table 1 (continued)Nomenclature Descriptions θ Death rate due to infection of T N Number of total population

Motivated by the importance of piecewise derivatives, we extend model (1) under the piecewise derivative with fractional order 0＜α≤1 as

denotes the piecewise classical Caputo derivative in[0,t2].Ins pired by the above discussion,we examine a SEITR-type mathematical model for NiV disease using the piece-wise classical and Caputo fractional order derivatives. Here, we state that the mentioned derivative has been recently utilized for studying various problems,which we refer to as[37,38].Eventually,the model’s formulation is presented.Second,the suggested model includes the fundamental reproduction number as well as disease-free and endemic equilibrium.We develop an attempt at existence theory for the model using the fixed point technique.Before studying dynamical models,it is necessary to examine the existence theory.The existence of real-world problems in the form of mathematical formulation is investigated properly by fixed point theory and other tools of nonlinear functional analysis.For some applications of fixed point theory, we refer to [39,40]. The said fixed point results have used for some problems in [41]. We also use the numerical procedure based on Adam Bashforth method to investigate the numerical analysis of the proposed model at different fractional orders. The mentioned numerical method is a powerful tool. The convergence and stability analysis of the Adams-Bashforth method have been discussed in detail in[42].

3 Basic Results

Here,some basic results are recollected from fractional calculus.

Definition 3.1.[41].The fractional derivative of g of order 0＜α≤1 in Caputo sense is given as

such that right hand side exists.

Definition 3.2.[41].Fractional integration of g w.r.t t of orderαis

such that right hand side exists.

Definition 3.3.[24]. Let g be differentiable function, then the traditional and arbitrary order piecewise derivative is recalled as

withα∈(0,1].

Definition 3.4.Let g be differentiable,then the classical and fractional order piecewise integration[24]is given as

Lemma 3.1.[24]The solution of

is given by

4 Qualitative Analysis

Some existence results are derived here. It should be kept in mind that prior to dealing with a dynamical system,existence theory is important to investigate for the prosed problem.As the model is used to estimate the population of humans,it is necessary that all of its parameters and variables for all time t be nonnegative.

Theorem 4.1.[33]. Let {S(0),E(0),I(0),T(0),R(0)≥0} ∈Then, the region of solution is invariant and all solutions are falling in the region representing by

Proof.By adding all equations of model(2),where N is the total papulation,we have

Take Laplace transform of and use N(0)=N0,as

yields after simplification

After evaluating,we have

whereEdenoted Mittag-Leffler function of one parameter.Further,att→∞in(3),we get

Also,it is straightforward to show that each compartment is non-negative.

Remark 1.Disease free equilibrium refers to a stable condition in which there is no infection or disease.The said point has been computed in[33]as

and the fundamental threshold number is given by

The said number predicts the transmission dynamics of disease in the community. If R0＜1, it indicates a decrease in disease in the community while R0＞1 shows a raise in infection in the society.The mentioned number can be computed by using the next generation matrix method which has been discussed in detail in[36].

Some 3D profiles of R0against different parameters are given in Fig.2.

Figure 2:3D configuration of R0 against different parameters of the model with N=100 millions

Lemma 4.1.Inview of Lemma 3.1,the solution of

is given by

where

Proof.Following the fashion of Lemma 3.1,we can easily obtain the solution above.LetEbe the space of all piecewise differentiable functions from[0,t2]to R,then obviously it is a complete normed space and hence the Banach space is expressed byE=C[0,t2],where 0＜t2＜∞with the norm

We use the growth and Lipschtiz conditions on the non-linear operator F to achieve the desired objectives.

(D1)Let LF＞0 be the Lipschitz constant and g,∈E,then the following holds

(D2)If CF＞0 and MF＞0,then the following holds

Theorem 4.2.Let piece-wise continuous function F be on the sub-intervals 0＜t ≤t1and t1＜t ≤t2of the interval[0,T]and also if the hypothesis(D2)holds,then the system(2)will have at least one solution.

Proof.Using Schauder fixed point theorem, let denote a closed and bounded subsetBofEasB={g ∈E:‖g‖≤R1,2,R1,2＞0}.

LetT:B→Busing(4),then

For g ∈B,it follows that

Thus,it proved that‖T(g)‖≤R1,2,which further implies thatTis bounded.Also,T(B)⊂B.Hence,Tmaps bounded set.For continuity,we use tm＜tn∈[0,t1]and take a classical derivative case first

By(7),we obtain tm→tnwhich yields

|T(g)(tn)-T(g)(tm)|→0, as tm→tn.

SinceTis uniformly continuous and bounded on the interval[0,t1]interval,it implies

‖T(g)(tn)-T(g)(tm)‖→0, as tm→tn.

It follows thatTis equi-continuous in this case.

Next,we take the other interval tm＜tn∈(t1,t2]as

We see from(8),that

|T(g)(tn)-T(g)(tm)|→0, as tm→tn.

Tis bounded on the interval(t1,t2]and also uniformly continuous.Thus,

‖T(g)(tn)-T(g)(tm)‖→0, as tm→tn.

Hence,Tis equi-continuous on the interval(t1,t2].

Therefore,Tis equi-continuous mapping over[0,t1]∪(t1,t2].On using the Arzelá-Ascoli theorem[39],which states if the operator is relatively compact,then it will be completely continuous.Hence,Tis completely continuous. Moreover, by using Schauder’s fixed point theorem [40], the suggested model(2)has at least one solution.

Theorem 4.3.By hypothesis(D1),if,then the proposed model under piecewise derivative has a unique approximate solution.

Proof.Consider the mappingT:E→E,and if g,∈E,then in case of[0,t1],one has

Using(9),one has

Fort∈(t1,t2]in the Caputo sense of derivative,we get from(10)

From Eq.(11),we get

Therefore,Tis a contraction operator in both cases.As a conclusion,the suggested model has a unique approximate solution under piecewise derivative.

5 Numerical Scheme

For piecewise derivable problem(2),we develop a numerical technique for the two sub-intervals of [0,t2] in order to perform numerical results. For the piecewise problem, the numerical algorithm will be the same as the integer order numerical algorithm defined in[24].We remark that the Adams-Bashforth method preserves the basic type of numerical stability associated with the usual one-step numerical methods, including Euler, backward Euler, trapezoidal, etc. Further, it is a(l+ 1)-step explicit method, whose truncation error is of sizeO(hl+2). For detailed convergence and stability of the Adams-Bashforth method,we refer to the reference[42]for further detail.

The solution of(2)under the piecewise derivative can be expressed as

We will first develop the strategy for the first equation in System(12)and then apply it to the rest of the equations.At t=tn+1,then we have

Eq.(13)is expressed in the Newton interpolation formula given in the reference as

For the rest of three equations,we can write the Newton interpolation scheme as shown below:

and

6 Numerical Simulations and Discussion

We illustrate the numerical simulation in Figures utilizing the acquire scheme of Newton Polynomial of classical and global piece-wise derivative notion.We divide the interval into two sub-intervals and analyze the first for integer order derivatives,whereas the second is tested for different fractional orders in the sense of Caputo,using the data in a table.Here,we simulate our model corresponding to the numerical data given in Table 2.Here,we present approximate solutions under piecewise fractional order differentiation of the proposed model using various fractional orders.The concerned graphs for approximate solutions when 0＜α≤0.6 are given in Figs.3-7.Those graphs where fractional order is 0.6＜α≤1 are presented in Figs.8-12,respectively.By simulating the results with clear crossover(multi-steps)behaviors in various compartments and by splitting the whole domain[0,30]into subintervals as 0＜t ≤5 and 5＜t ≤30, we see during treatment the corresponding population of susceptible, exposed, infected, treated, and recovered in the corresponding point appear to start a multi-step behavior.The decay process in smaller orders is faster while growth or raise occurs rapidly in larger fractional orders.The nature of the model under our study indicated a clear date when the infection transmitted from some infected sources to the healthy ones and also there is a termination if the infected ones got rid of disease either by treatment or death.Hence,the date of the evolution can be well known,more precisely as there is no peculiarity in the dynamic of the transmission,but a crossover behavior can be clearly observed in this case.The process is non-Markovian and has a crossover(multistep)behavior.These remarks can be clearly understood from the Figs.3-12.Initially,the population dynamics of each compartment are progressing normally, but after some sudden or abrupt changes occur,the dynamics of various compartments start multi-step behavior.Here,it should be kept in mind that such scenarios cannot be achieved by using usual fractional order or integer order derivatives.

Table 2: Numerical values with description and sources

Figure 3:Graphical presentation of approximate solutions of susceptible compartment against different values of fractional orders such that 0 ＜α ≤0.6

Figure 4:Graphical presentation of approximate solutions of exposed compartment against different values of fractional orders such that 0 ＜α ≤0.6

Figure 5:Graphical presentation of approximate solutions of infected compartment against different values of fractional orders such that 0 ＜α ≤0.6

Figure 6: Graphical presentation of approximate solutions of treated compartment against different values of fractional orders such that 0 ＜α ≤0.6

Figure 7:Graphical presentation of approximate solutions of recovered compartment against different values of fractional orders such that 0 ＜α ≤0.6

Figure 8:Graphical presentation of approximate solutions of susceptible compartment against different values of fractional orders such that 0.6 ＜α ≤1

Figure 9:Graphical presentation of approximate solutions of exposed compartment against different values of fractional orders such that 0.6 ＜α ≤1

Figure 10:Graphical presentation of approximate solutions of infected compartment against different values of fractional orders such that 0.6 ＜α ≤1

Figure 11:Graphical presentation of approximate solutions of treated compartment against different values of fractional orders such that 0.6 ＜α ≤1

Figure 12:Graphical presentation of approximate solutions of recovered compartment against different values of fractional orders such that 0.6 ＜α ≤1

7 Conclusion

In this article,we have developed a qualitative and computational analysis of a NiV disease.Based on the notion of the newly explored piecewise derivative of fractional order,we have investigated the model to study the crossover behavior in the dynamics of the disease.First of all,the invariance of the model and the feasible region has been established in the sense of the Caputo fractional derivative via using Laplace transform.Some remarks about thresholds number and their 3D profile corresponding to different parameters have been given. We have developed necessary and sufficient conditions for the existence and uniqueness of approximate solutions via fixed point results of Banach and Schauder.Utilizing the Newtonian polynomials of numerical interpolation,we have established algorithms for simulation purposes of the model. We have constructed the required scheme for numerical findings and their graphical presentation. We have simulated the model on different sets of fractional order by dividing the domain of time into two sub-intervals. Graphical presentations have indicated the crossover(multi-step)behavior in the dynamical study of the aforesaid model.In the future,this model can be investigated by using anABCtype derivative which has a nonlocal and nonsingular kernel.It would be expected that more dynamic properties may appear for better understanding.

Acknowledgement:The authors K. Shah, T. Abdeljawad, and B. Abdalla thank Prince Sultan University for support through the TAS Research Lab.

Funding Statement:The research was financially supported by Prince Sultan University.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.