# Travelling Wave Analysis of a Diffusive COVID-19 Model

Authors: C. M. Wachira; G. O. Lawi; L. O. Omondi
Journal: Journal of Applied Mathematics
 (2022)

---

## Abstract

In this paper, a mathematical model based on a system of nonlinear parabolic partial differential equations is developed to investigate the effect of human mobility on the dynamics of coronavirus 2019 (COVID-19) disease. Positivity and boundedness of the model solutions are shown. The existence of the disease-free, the endemic equilibria, and the travelling wave solutions of the model are shown. From the numerical analysis, it is shown that human mobility plays a crucial role in the disease transmission. Therefore, interventions that affect diffusion (human mobility), such as lock-down, travel restrictions, and cessation of movement, may play a significant role in controlling and preventing the spread of COVID-19.


---

## Body

## 1. Introduction

COVID-19 disease, caused by SARS-COV-2 virus emerged from the city of Wuhan, Hubei Province, China in late December 2019 and quickly spread all over the world [1]. This global spread was greatly facilitated by movement of infectious individuals from Wuhan city (the epicenter) to other parts of the world where there were no COVID-19 infections. The disease was declared a global health pandemic by the World Health Organization (WHO) on 11th of March, 2020. The virus is primarily transmitted to humans through inhaling respiratory droplets produced by sneezing, coughing, and conversing with an infectious person (direct transmission) or indirectly through contact with contaminated surfaces (indirect transmission) [2]. COVID-19 incubation period, which is the time between virus exposure and symptom onset, is on average 5–6 days, but can be as long as 14 days [3]. Fever, a dry cough, and fatigue are the most common symptoms of COVID-19. Body aches and pains, a sore throat, diarrhoea, headache, loss of taste and smell, breathing difficulties, a rash on the skin, and discoloration of the fingers and toes are some of the other symptoms [1, 2].There were no antiviral prophylaxis or therapeutics available during the early phases of the COVID-19 outbreak, so all affected countries relied on nonpharmaceutical interventions such as isolation and quarantine of infected individuals, contact tracing, and population lockdown [4]. These interventions are aimed at containing and mitigating the epidemic, that was a threat to many health-care systems. These preventive measures appeared to be effective in lowering the number of fatalities and hospitalizations [4, 5]. COVID-19 vaccine candidates have been developed and are already in use [6]. Vaccines such as BNT162b2 (Pfizer), mRNA1273 (Moderna), ChAd0x1 nCoV-19 (Astra Zeneca), Sputnik, Sinopharm, and Johnson & Johnson vaccines are currently being administered worldwide.The dynamics of COVID-19 pandemic have been complex in both time and space. These dynamics are caused by a number of factors, including the spatial distribution of the population’s social and economic levels and the patterns of human mobility within a given country [7]. Human mobility may be due to travel, migration between localities, countries, or towns. When infected people travel to infection-free areas, there is a good risk that they will pass their infections on to locals. The transfer of infections from a region of high concentration to a region of low concentration causes a wave of transfer of infection.A variety of mathematical models have been used widely in various biological applications such as ecology, population dynamics, tumor growth (cancer), immunology, and epidemiology [8, 9]. A number of epidemics have been modeled using reaction-diffusion equations with the focus to demonstrate the existence of travelling wave solutions in work such as in references [10, 11]. Reaction-diffusion epidemic models are generally standard epidemic models that assume there is movement and interaction of individuals in a given population in space.Human mobility restrictions such as lock-down, travel restrictions, and cessation of movement are regarded as an effective strategy for controlling disease spread. However, it is still unclear whether mobility restriction is a proportional response to control the ongoing COVID-19 pandemic [12]. In this paper, the existence of travelling wave solutions as well as the effect of human mobility (diffusion) on the transmission dynamics of COVID-19 disease are investigated.


## 2. Model Formulation

In this study the transmission dynamics of COVID-19 is assumed to be determined by the spatial location,x, as well as time, t. Despite the fact that the underlying conditions mentioned in the previous section predisposes someone to a higher risk of contracting COVID-19 disease, because the focus is on human mobility, the study assumes the differential susceptibility. Individuals in this category are referred to as susceptible individuals, which are denoted by the symbol Sx,t. When susceptible individuals are exposed to the virus, they undergo an incubation period, which is the time between virus exposure and onset of symptoms. This is on average 5–6 days, but can be as long as 14 days [3]. The incubating individuals here are categorised as the exposed class, denoted by Ex,t. Following the end of the exposure period, the incubating individuals transit into the infection class denoted by Ix,t. The infectious potential of COVID-19 is significantly greater just before or within the first five days of symptom onset [13]. Therefore, the infected class Ix,t in this study will include both individuals with and without symptoms. Depending on the severity of symptoms and the intervention strategies in place, the infected individuals may die or recover from the disease. Thus, the class Rx,t denotes the number of individual who have recovered from the infection.It is assumed that susceptible humans are recruited at a constant rateΛ. The susceptible humans become infected through a force of infection βSI/1+ηI, where β is the transmission coefficient, and η is the half-maximal human saturation constant of the infected individuals. It is considered to be of Michaelis–Menten form to account for saturation of the human infection [14]. The saturated incidence given by βSI/1+ηI is reasonable due to the fact that as the infected individuals increase, they reach a saturation point. The number of infected individuals decreases due to psychological effects, behavioral changes, or preventive measures taken by the affected individuals [15]. The rate of progression from the exposed to infected classes is assumed to be ε. The rates of natural and disease-induced mortality are μ and δ, respectively. The recovery rate from the infection is taken as λ. The diffusivity of susceptible, exposed, infected, and recovered individuals are d1, d2, d3, and d4, respectively.The corresponding proposed mathematical equations arising from the preceding description can be described by the following system of partial differential equations:(1)∂S∂t−d1∇2S=Λ−βSI1+ηI−μS,(2)∂E∂t−d2∇2E=βSI1+ηI−μ+εE,(3)∂E∂t−d2∇2E=βSI1+ηI−μ+εE,(4)∂I∂t−d3∇2I=εE−μ+λ+δI,(5)∂R∂t−d4∇2R=λI−μR,where Sx,t,Ex,t,Ix,t,Rx,t∈Ψ×T⊂ℝ+4×ℝ and Ψ are bounded domains in ℝ+4. Human mobility is considered in one-dimensional space i.e., ∇2=∂2∂x2. The result can be extended in higher-dimensional space, such as 2-dimensional space.


## 3. Well-Posedness of the Model

To ensure that the solutions of the model reflect the biological reality, the fundamental question of positivity and boundedness of the solutions is addressed in this section. The model is analyzed with zero flux boundary conditions (no movement across the boundary of∂Ψ).
(6)∂S∂n=∂E∂n=∂I∂n=∂R∂n,where ∂/∂n, denotes the outward normal derivative on ∂Ψ. Let the initial conditions be given by
(7)Sx,t=Sx,0≥0,Ex,t=Ex,0≥0,Ix,t=Ix,0≥0,Rx,t=Rx,0≥0.Forx∈−∞,+∞, Ψ is a bounded domain in ℝ+4 with smooth boundary ∂Ψ, and t≥0.Proposition 1.
Suppose that the initial condition (7) holds, then the solutions of model (5) are nonnegative in 0,+∞ for all t≥0.Proof.
Model (5) can be written as an abstract Banach space X=C¯Ψ×CΨ¯ in the form
(8)u′=Aut+Fut,t>0u0=u0∈X,where u=S,E,I,RT, u0=Sx,0,Ex,0,Ix,0Rx,0T, Aut=d1S,d2E,d3I,d4RT, and
(9)Fut=Λ−βSI1+ηI−μSβSI1+ηI−μ+εEεE−μ+λ+δIλI−μR.
Clearly,F is locally Lipschitz in X, hence model (5) has local solutions on the interval 0,Tmax, where Tmax is the minimum existence time for the solutions of model (5) [16].
Model (5) can also be written in the following form:
(10)dSdt−d1∇2S=F1S,E,I,RdEdt−d2∇2E=F2S,E,I,RdIdt−d3∇2I=F3S,E,I,RdRdt−d4∇2R=F4S,E,I,R.
The functionsFiS,E,I,R and i=1,2,3,4,5 are continuously differentiable and they satisfy the following conditions: F10,E,I,R=Λ≥0,F2S,0,I,R=0≥0,F3S,E,0,R=0≥0,F4S,E,I,0=0≥0, for all S,E,I,R≥0. Since S,E,I,R≥0 with positive initial conditions, then solutions of model (5) are positive. This completes the proof.Proposition 2.
The solutions of model (5) are bounded in the region Ψ×T for all t≥0.Proof.
To check for the boundedness, adding all the equations in model (5) and setting D=maxd1,d2,d3,d4 to obtain
(11)∂Nx,t∂t≤Λ−μN+D∇2N.
The inequality (11) has a unique solution of the form
(12)∂Nx,t∂t≤Λ−∫−∞∞μN4Dπte−x2/4Dtdx,where the fundamental solution of inequality (12) is given by
(13)Kt,x=14Dπte−x2/4Dt,x∈ℝt>0,0,x∈ℝt<0.
For one-dimensional reaction-diffusion equation, equation (13) satisfies the following
(14)Kt,x=∫−∞+∞14Dπte−x2/4Dtdx=1.
Hence(15)dNdt≤Λ−μN.
SinceNt,x>0, solving this inequality (15) and taking the limit as t⟶∞. Hence by variation of constants formula, yields
(16)limsupt⟶∞Nt,x≤Λμ
HenceNt,x is bounded.From Proposition1 and Proposition 2, it is clear that solutions of model (5) are positive and bounded for t≥0. Therefore, model (5) is mathematically and epidemiologically meaningful and it is now sufficient to consider its solutions in Ψ×T.


## 4. Existence of Equilibrium Points

The disease-free equilibrium point is the state at which there is no infection spreading in the population.Proposition 3.
Model (5) has a disease-free equilibrium denoted by E0=Λ/μ,0,0,0 where S=Λ/μ, and E=I=R=0.The endemic equilibrium is the state at which the disease remains persistent in a given population.Proposition 4.
The endemic equilibriumE∗ of the (5) exists whenever I∗>0 for R0>1 and is given by E∗=S∗,E∗,I∗,R∗Proof.
. Settingdi=0 for i=1,2,3,4 in model (5), then from first equation of model (5) becomes
(17)S∗=Λβ1+ηI∗βI∗+μ1+ηI∗.
Substituting equation (17) into second equation of model (5) yields
(18)E∗=βΛI∗μ+εβI∗+μ1+ηI∗.
Now substituting equation (18) into the third equation of model (5) yields
(19)I∗=μR0−1β+μη.
It is clear from equation (19) that I∗>0 when R0>1. This implies that there exists an endemic equilibrium for model (5) whenever R0>1.


## 5. Existence of the Travelling Wave Solution

Traveling wave solutions have a specific shape and incorporate time and spatial variables through a moving coordinate. A solution of this type describes the spatial transition from one steady state to another. The traveling wave fronts that connect disease-free and endemic states are determined in this section. Model (5) travelling wave solutions are shown using the approach used in references [10, 11].Define a new variable of the formz=kx−vt, v>0, where v is the propagating wave speed of model (5) in one-dimensional space and k is the wave number. Let Sx,t=Sz,Ex,t=Ez,Ix,t=Iz, and Rx,t=Rz. Then model (5) can be transformed into the following set:
(20)−vkS′−d1k2S″=Λ−βSI1+ηI−μS−vkE′−d2k2E″=βSI1+ηI−μ+εE−vkI′−d3k2I″=εE−μ+λ+δI−vkR′−d4k2R″=λI−μR,where ′ denotes the derivative with respect to the new variable z. System (20) can be expressed as follows:
(21)d1k2S″+vkS′+Λ−βSI1+ηI−μS=0d2k2E″+vkE′+βSI1+ηI−μ+εE=0d3k2I″+vkI′+εE−μ+λ+δI=0d4k2R″+vkR′+λI−μR=0,which can further be expressed as
(22)S″+H1S′+F1S,E,I,R=0E″+H2E′+F2S,E,I,R=0I″+H3I′+F3S,E,I,R=0R″+H4R′+F4S,E,I,R=0,where
(23)F1=1d1k2Λ−βSI1+ηI−μS,F2=1d2k2βSI1+ηI−μ+εE,F3=1d3k2εE−μ+λ+δI,F4=1d4k2λI−μR,Hi=vdik,i=1,2,3,4.Setting the variablesx1=S′, x2=E′, x3=I′, and x4=R′, then model (22) is transformed into a system of first order differential equations.
(24)Y=x1,S,x2,E,x3,I,x4,R∈ℝ8,(25)dYdz=fY=−H1x1−F1x1−H2x2−F2x2−H3x3−F3x3−H4x4−F4x4,with boundary conditions
(26)limz⟶−∞x1,S,x2,E,x3,I,x4,R=E0,limz⟶+∞x1,S,x2,E,x3,I,x4,R=E∗,where E0 is the disease-free equilibrium point and E∗ is the endemic equilibrium point. Then a travelling wave solution is a trajectory that joins E0 and E∗.Computing the Jacobian matrix of equation (25) yields
(27)J=−H1−1d1k2βI1+ηI−μ000−1d1k2βS1+ηI200100000000−1d1k2βI1+ηI−H2μ+εd1k20−βSd2k21+ηI20000100000000εd3k2−H3μ+λ+δd3k2000000100000000λd4k2−H4μd4k200000010.Evaluating matrix (28) at the Disease-free equilibrium E0=Λ/μ,0,0,0,0 yields
(28)JE0=−vd1kμd1k2000βΛμd1k2001000000000−vd2kμ+εd2k20−βΛμd2k20000100000000εd3k−vd3μ+λ+δd3k2000000100000000λd4k2−vd4kμd4k200000010.The eigenvalues of the Jacobian matrix (28) evaluated at the disease-free equilibrium, ξ1,⋯,ξ8 are
(29)ξ1,2=−v±v2+4μd12d1kξ3,4=−v±v2+4μ+εd22d2k,ξ5,6=−v±v2+4μ+λ+δd32d3k,ξ7,8=−v±v2+4μd42d4k.Clearly,ξ1,⋯,ξ8 are real-valued eigenvalues, since all the parameters of the model are positive. Therefore, a traveling wave profile, that travels at a speed v>0, exists that connects the disease-free and endemic equilibrium states. This means that if infectious individuals diffuse into a susceptible population, a transition zone of infection will form.


## 6. Numerical Simulations and Discussion

Simulation analysis of model (5) is presented to graphically illustrate the behaviour of the solutions of the model with varying rates of diffusivity. For purposes of simulation the parameter values are stated in Table 1 and unless otherwise stated the initial populations are taken to be S0,x=3000, E0,x=2000, I0,x=1500, and R0,x=1000.Table 1
The descriptive summary of the model parameters.
ParameterDescriptionUnit/value unitsSourceΛRecruitment rate3.178×10−5 day−1[17]μNatural death rate3.91×10−5 day−1[17]δDisease mortality rate1.03×10−6 day−1[17]λHuman recovery rate0.125 day−1[5]βTransmission probability0.02 day−1EstimatedεTransition from E to I0.0877 day−1[18]ηHuman saturation constant0.05Estimatedd1Susceptible diffusivity constant(0-1.0) km day−1Variabled2Exposed diffusivity constant(0-1.0) km day−1Variabled3Infected diffusivity constant(0-1.0) km day−1Variabled4Recovered diffusivity constant(0-1.0) km day−1VariableFigures1(a) and 1(b) show the simulations of the susceptible humans for varying values of diffusivity rate d1 at 0.1 and 0.9, respectively. The number of the susceptible humans reduces faster when diffusivity rate is high and also reduces with an increase in space. When d1=0.9, it would take less than 10 days for a pool of susceptible humans to be exposed to the virus. When diffusivity rate is low, d1=0.1, it would take about 20 days for a pool of susceptible humans to become exposed once they come into contact with the virus.Figure 1
Plots (a) and (b) show the simulations of the susceptible humans corresponding to rate of diffusion, (a)d1=0.1 km day−1, and (b) d1=0.9 km day−1..
(a)(b)Figure2(a) and 2(b) show the spatiotemporal simulations of the infected humans for varying values of diffusivity rate d3. They show that as the rate of diffusion (human mobility), d3 increases the number of infected individuals increases. This means that, human mobility has an effect in the transmission dynamics of COVID-19. Therefore, use of control interventions such as lock-down, travel restrictions, and cessation of movement that reduce human mobility may play a significant role in reducing COVID-19 spread. Infectious individuals will spread the infection to noninfected areas until an equilibrium is reached. The numerical results show that the disease persists in the population due to the movement of infectious individuals. The formation of a wave profile indicates that the introduction of an infective individual into a susceptible population results in infection propagation.Figure 2
Plots (a) and (b) show the simulations of the Infected humans corresponding to the rate of diffusion, (a)d3=0.1 km day−1, and (b) d3=0.9 km day−1.
(a)(b)


## 7. Conclusion

The effect of human mobility on COVID-19 dynamics is investigated in this paper using a nonlinear parabolic partial differential equation. Well-posedness of the model and existence of steady states are shown. Existence of traveling wave solutions of the model are shown. The model analysis revealed that when infected individuals are introduced into a fully susceptible population, there is a moving transition “a propagating wave” of infections that connects the disease-free equilibrium to the endemic equilibrium. These waves propagate at speedv>0, which connects the two equilibria points. From the numerical analysis, it is shown that human mobility is critical for disease transmission. As a result, control interventions such as lock-down, travel restrictions, and cessation of movement that reduce human mobility may play a significant role in reducing COVID-19 spread. These interventions are aimed at creating barriers that will slow the spread of the virus in the affected areas.



---

*Source: 6052274-2022-10-07.xml*