Time-dependent non-homogeneous stochastic epidemic model of SIR type

Mireia Besalú; Giulia Binotto; Mireia Besalú; Giulia Binotto

doi:10.3934/math.20231181

AIMS Mathematics

2023, Volume 8, Issue 10: 23218-23246. doi: 10.3934/math.20231181

Previous Article Next Article

Research article

Time-dependent non-homogeneous stochastic epidemic model of SIR type

Mireia Besalú ¹,
Giulia Binotto ^{2
,
,}

1.
Departament de Genètica, Microbiologia i Estadística, Universitat de Barcelona, Avinguda Diagonal 643, 08028-Barcelona, Spain
2.
Department de Matemàtiques, Universitat Autònoma de Barcelona, 08193-Bellaterra, Spain

Received: 05 May 2023 Revised: 20 June 2023 Accepted: 24 June 2023 Published: 20 July 2023
MSC : 60H10, 60J10, 92D30

To better describe the spread of a disease, we extend a discrete time stochastic SIR-type epidemic model of Tuckwell and Williams. We assume the dependence on time of the number of daily encounters and include a parameter to represent a possible quarantine of the infectious individuals. We provide an analytic description of this Markovian model and investigate its dynamics. Both a diffusion approximation and the basic reproduction number are derived. Through several simulations, we show how the evolution of a disease is affected by the distribution of the number of daily encounters and its dependence on time. Finally, we show how the appropriate choice of this parameter allows a suitable application of our model to two real diseases.

Keywords:

Citation: Mireia Besalú, Giulia Binotto. Time-dependent non-homogeneous stochastic epidemic model of SIR type[J]. AIMS Mathematics, 2023, 8(10): 23218-23246. doi: 10.3934/math.20231181

Related Papers:

[1]	Yuanfu Shao . Dynamics and optimal harvesting of a stochastic predator-prey system with regime switching, S-type distributed time delays and Lévy jumps. AIMS Mathematics, 2022, 7(3): 4068-4093. doi: 10.3934/math.2022225
[2]	Xiaodong Wang, Kai Wang, Zhidong Teng . Global dynamics and density function in a class of stochastic SVI epidemic models with Lévy jumps and nonlinear incidence. AIMS Mathematics, 2023, 8(2): 2829-2855. doi: 10.3934/math.2023148
[3]	Shuo Ma, Jiangman Li, Qiang Li, Ruonan Liu . Adaptive exponential synchronization of impulsive coupled neutral stochastic neural networks with Lévy noise and probabilistic delays under non-Lipschitz conditions. AIMS Mathematics, 2024, 9(9): 24912-24933. doi: 10.3934/math.20241214
[4]	Zhengwen Yin, Yuanshun Tan . Threshold dynamics of stochastic SIRSW infectious disease model with multiparameter perturbation. AIMS Mathematics, 2024, 9(12): 33467-33492. doi: 10.3934/math.20241597
[5]	Tian Xu, Ailong Wu . Stabilization of nonlinear hybrid stochastic time-delay neural networks with Lévy noise using discrete-time feedback control. AIMS Mathematics, 2024, 9(10): 27080-27101. doi: 10.3934/math.20241317
[6]	Yassine Sabbar, Aeshah A. Raezah . Modeling mosquito-borne disease dynamics via stochastic differential equations and generalized tempered stable distribution. AIMS Mathematics, 2024, 9(8): 22454-22485. doi: 10.3934/math.20241092
[7]	Chuanfu Chai, Yuanfu Shao, Yaping Wang . Analysis of a Holling-type IV stochastic prey-predator system with anti-predatory behavior and Lévy noise. AIMS Mathematics, 2023, 8(9): 21033-21054. doi: 10.3934/math.20231071
[8]	Yassine Sabbar, Aeshah A. Raezah, Mohammed Moumni . Enhancing epidemic modeling: exploring heavy-tailed dynamics with the generalized tempered stable distribution. AIMS Mathematics, 2024, 9(10): 29496-29528. doi: 10.3934/math.20241429
[9]	Hong Qiu, Yanzhang Huo, Tianhui Ma . Dynamical analysis of a stochastic hybrid predator-prey model with Beddington-DeAngelis functional response and Lévy jumps. AIMS Mathematics, 2022, 7(8): 14492-14512. doi: 10.3934/math.2022799
[10]	Yassine Sabbar, Aeshah A. Raezah . Threshold analysis of an algae-zooplankton model incorporating general interaction rates and nonlinear independent stochastic components. AIMS Mathematics, 2024, 9(7): 18211-18235. doi: 10.3934/math.2024889

Abstract

1. Introduction

In the last few decades, an outstanding advancement has been witnessed in nonlinear sciences and engineering fields. Many scientists showed keen interest in finding the exact and numerical solutions for the nonlinear PDEs. Numerous techniques were devised in this regard which includes the GERFM method ^[1], modified variational iteration algorithm-II ^[3], enhanced ${(}\frac{{G}^{^{{\prime }}}}{G}$ )-expansion method ^[4], the direct algebraic method ^[5], the extended trial equation method ^[6], the generalized fractional integral conditions method ^[7], the modified simple equation method ^[8], the Monch's theorem method ^[9], the extended modified mapping method ^[10], the reductive perturbation method ^[11], the new probability transformation method ^[12] and the differential transformation method ^[13]. In future work for more related extensions or generalizations of the results these references may be very helpful ^{[2,32,33,34,35,36]}.

The time fractional derivatives in the fractional reaction-diffusion model describes the process relating to the physical phenomena, physically known as the historical dependence. The space fractional derivative explains the path dependence and global correlation properties of physical processes, that is, the global dependence. The reaction diffusion equation has a dynamic role in dissipative dynamical systems as studied by various biologists ^[38], scientists and engineers ^[14]. The nonlinear form of this model has found a number of applications in numerous branches of biology, physics and chemistry ^[14,15]. This model has also been useful for other areas of science and effectively generalized by employing the theory of fractional calculus, for instance see ^[15,16]. Diffusion-wave equations involving Caputo's derivative ^[17,18], Riemann-Liouville derivatives ^[19] have been discussed by various researchers. Anomalous dispersion equations can be explained by fractional derivative ^[20]. An extensive variety of exact methods which have been applied for exact solutions of the fractional nonlinear reaction diffusion equation, for example see ^{[14,15,16,17,18,19,20]} and references there in.

To interpret numerous physical phenomena in some special fields of science and engineering, nonlinear evolution equations are extensively used as models especially in solid-state physics and plasma physics. Finding the exact solutions of NLEEs is a key role in the study of these physical phenomena ^[39]. A lot of research work has been carried out during the past decades for evaluating the exact and numerical solutions of many nonlinear evolution equations. Among them are homotopy analysis method ^[21], modified exp-function method ^[22], ${(}\frac{{G}^{^{{\prime } }}}{G}$ )-expansion method ^[23], exp-function method ^[24], homotopy perturbation method ^[25], Jacobi elliptic function method ^[26], sub equation function method ^[27], kudryashov method ^[28], and so on.We can be expressed the exact solutions of FPDE via exp $(-\varphi \left(\eta \right))$ .

$\begin{equation} \left( \varphi ^{\prime }\left( \eta \right) \right) = \mathrm{exp}\left( \mathrm{-}\varphi \left( \eta \right) \right) +\mu \ \mathrm{exp}\mathrm{(} \varphi \left( \eta \right) \mathrm{)}+\lambda \end{equation}$

(1.1)

The article is arranged as follows: Section 1 represents the introduction of the article. In section 2, we have explained the Caputo's fractional derivative. In section 3, we have interpreted the exp $(-\varphi \left(\eta \right))$ -expansion method. In section 4, we use this method to explore the reaction-diffusion model. In section 5 and 6, graphical representation and physical interpretation are explained. In section 7 and 8, we have interpreted the results, discussions and conclusion.

1.1. Caputo's fractional derivative

Property 1: ^[29] A function $f(x, t),$ where $x > 0$ is considered as $C_{\alpha }$ . Here $\alpha \in \mathfrak{R},$ if $\exists$ a ${\mathbb{R}}$ and $\left(p > \alpha \right),$ such that

$\begin{equation} f\left( x\right) = x^{p}f_{1}\left( x\right) \end{equation}$

(1.2)

$f_{1}\left( x\right) \in C[0,\infty ).$ Where

$f_{1}\left( x\right) \in C[0,\infty )$

Property2: ^[29] A function $f\left(x, t\right),$ where $x > 0$ is considered to be in space $C_{\alpha }^{m}$ . Here $m\in { \mathbb{N}}\cup \left\{ 0\right\},$ if $f^{(m)}\in C_{\alpha }.$

$\begin{equation} I_{t}^{\mu }f\left( x,t\right) = \frac{\mathrm{1}}{\Gamma \left( \mu \right) } \int_{0}^{t}{{\left( t-\mathrm{T}\right) }^{\mu -1}}f\left( x,\mathrm{T} \right) \mathrm{d}\mathrm{T},\ t > 0. \end{equation}$

(1.3)

Property3: ^[29] Suppose $f\in C_{\alpha }$ and $\alpha \geq -1$ , then the Riemann Liouville integral $\mu,$ where $\mu > 0$ is given by

Property4: ^[29] A fractional Caputo derivative of $f$ with respect to $t$ , where ${f\in C}_{-1}^{m}, \ m\in {\mathbb{N}}\cup \left\{ 0\right\},$ is given as

$\begin{equation} D_{t}^{\mu }f\left( x,t\right) = \frac{{\partial }^{m}}{\partial t^{m}} f\left( x,t\right) ,\ \mu = m \end{equation}$

(1.4)

$\begin{equation} = I_{t}^{m-\mu }\frac{{\partial }^{m}}{\partial t^{m}}f\left( x,t\right) ,m-1\leq \mu < m \end{equation}$

(1.5)

Note that

$\begin{equation} I_{t}^{\mu }D_{t}^{\mu }f\left( x,t\right) = f\left( x,t\right) -\sum\limits_{k = 0}^{m-1}{\frac{{\partial }^{k}\ f}{\partial t^{k}}\left( x,0\right) }\frac{t^{k}}{k!},\ \ m-1 < \mu \leq m,\ m\in {\mathbb{N}} \end{equation}$

(1.6)

$\begin{equation} I_{t}^{\mu }t^{\nu } = \frac{\Gamma \left( \nu +1\right) }{\Gamma \left( \mu +\nu +1\right) }t^{\mu +\nu }. \end{equation}$

(1.7)

1.2. Interpretation of the method

Consider the fractional partial differential equation,

$\begin{equation} \varphi \left( u,D_{t}^{\alpha }u,u_{x},u_{xx},D_{t}^{2\alpha }u,D_{t}^{\alpha }u_{x},\dots \right) = 0,\;t > 0,x\in R,0\leq \alpha \leq 1, \end{equation}$

(1.8)

where $D_{t}^{\alpha }u, \ D_{x}^{\alpha }u, \ D_{xx}^{\alpha }u$ are derivatives, $u\left(\eta \right) = u\left(x, t\right)$ . For solving Eq 1.8, we follow:

Step 1: Using a transformation, we get,

$\begin{equation} \eta = x\pm V\frac{t^{\alpha }}{\Gamma \left( 1+\alpha \right) },\; u = u\left( \eta \right) , \end{equation}$

(1.9)

where constant $V$ is a nonzero. By substituting Eq 1.9 in Eq 1.8 yields ODE:

$\begin{equation} \varphi \left( u,\pm Vu^{\mathrm{^{\prime }}},ku^{\mathrm{^{\prime }} },V^{2}u^{\mathrm{^{\prime \prime }}},k^{2}u^{\mathrm{^{\prime \prime }} },\dots \right) = 0. \end{equation}$

(1.10)

Step 2: Assuming the traveling wave solution

$\begin{equation} u\left( \eta \right) = \sum\limits_{n = 0}^{M}{a_{n}{\left( {\mathrm{e}}^{-\varphi \left( \eta \right) }\right) }^{n}}, \end{equation}$

(1.11)

where $\varphi \left(\eta \right)$ satisfies the following equation:

$\begin{equation} {\varphi }^{^{{\prime }}}\left( \eta \right) = {\mathrm{e}}^{-\varphi \left( \eta \right) }+\lambda +\mu {\mathrm{e}}^{\varphi \left( \eta \right) } \end{equation}$

(1.12)

$\begin{equation*} {\varphi }^{^{\prime \prime }}\left( \eta \right) = -\lambda {\mathrm{e}} ^{-\varphi \left( \eta \right) }-{\mathrm{e}}^{-2\varphi \left( \eta \right) }+\mu \lambda {\mathrm{e}}^{\varphi \left( \eta \right) }+{\mu }^{2}{\mathrm{ e}}^{2\varphi \left( \eta \right) }, \end{equation*}$

where prime indicates derivative w.r.t. $\eta.$ The solutions of Eq 1.12 are written in the form of different cases.

Class1: when $\mu \neq 0$ and ${\lambda }^{2}-4\mu > 0,$ we have

$\begin{equation} \varphi \left( \eta \right) \mathrm{ = ln}\left\{ \frac{\mathrm{1}}{\mathrm{2} \mu }\left( \mathrm{-}\lambda \mathrm{-}\lambda \mathrm{tanh}\left( \frac{ \lambda }{\mathrm{2}}\mathrm{(}c_{1}+\eta )\right) \right) \right\} . \end{equation}$

(1.13)

where $\lambda = \sqrt{\left(-4\mu +{\lambda }^{2}\right) }$

Class 2: when, $\mu \neq 0$ and ${\lambda }^{2}-4\mu < 0,$ we have

$\begin{equation} \varphi \left( \eta \right) \mathrm{ = ln}\left\{ \frac{\mathrm{1}}{\mathrm{2} \mu }\left( \mathrm{-}\lambda \mathrm{+}\lambda \mathrm{tan}\left( \frac{ \lambda }{\mathrm{2}}\mathrm{(}c_{1}+\eta )\right) \right) \right\} . \end{equation}$

(1.14)

where, $\lambda = \sqrt{\left(-4\mu +{\lambda }^{2}\right) }$

Class 3: when, $\lambda \neq 0, \mu = 0$ and ${\lambda }^{2}-4\mu > 0,$ we have

$\begin{equation*} \varphi \left( \eta \right) \mathrm{ = -ln}\left\{ \frac{\lambda }{\mathrm{-} \mathrm{1+}{\exp \left( \lambda \left( k+\eta \right) \right) \ }}\right\} \end{equation*}$

$\varphi \left(\eta \right) \mathrm{ = -ln}\left\{ \frac{\lambda }{\mathrm{-} \mathrm{1+}{\exp \left(\lambda \left(k+\eta \right) \right) \ }}\right\}.$

Class 4: when, $\lambda, \mu \neq 0$ and ${\lambda }^{2}-4\mu = 0,$ we have

$\begin{equation*} \varphi \left( \eta \right) \mathrm{ = ln}\left\{ 2\frac{\left( 2+\lambda \left( k+\eta \right) \right) }{\left( {\lambda }^{2}\left( \eta +k\right) \right) }\right\}. \end{equation*}$

Class 5: when, $\lambda, \mu = 0$ and ${\lambda }^{2}-4\mu = 0,$ we have

$\begin{equation} \varphi \left( \eta \right) \mathrm{ = }\left[ \mathrm{ln}\left\{ \eta +k\right\} \right] , \end{equation}$

(1.15)

Step 3: Using the homogeneous balancing principal, in (10), we attain M. In view of Eq (11), Eq (10) and Eq (12), we obtain a system of equations with these parameters, a_n, $\lambda, \mu$ . We substitute the values in Eq (11) and Eq (12) and obtained the results of Eq (8).

2. Solution procedure

Suppose the reaction-diffusion equation is,

$\begin{equation} D_{t}^{2\alpha }u+\delta u_{xx}+\beta u+\gamma u^{3} = 0,\;0 < \alpha \leq 1, \end{equation}$

(2.1)

where $\delta$ , $\beta$ and $\gamma$ are parameters without zero, setting, $\delta = a$ , $\beta = b$ and $\gamma = c$ and changing Eq 2.1 into an ODE.

$\begin{equation} V^{2}u^{^{\prime \prime }}+au^{^{\prime \prime }}+bu+cu^{3} = 0, \end{equation}$

(2.2)

where prime represents the derivative w. r. t. $\eta$ .With the help of balancing principal, $u^{\prime \prime }$ and $u^{3}$ , we attain, $M = 1$ .

Rewriting the solution of Eq 2.2 we get,

$\begin{equation} u\left( \eta \right) = \left[ a_{0}+a_{1}\left( \mathrm{exp}\left( -\varphi \left( \eta \right) \right) \right) \right] , \end{equation}$

(2.3)

where $a_{0}$ , $a_{1}\neq 0$ are constants, while $\lambda, \mu$ are some constants.

Substituting $u, u^{^{\prime \prime }}$ and ${u}^{3}$ in Eq 2.2, we get the solution sets as.

Solution set 1

$\begin{equation*} \lambda = 0,\ C = C,V = -\frac{1}{2}\frac{\sqrt{2}\sqrt{-\mu \left( 2a\mu +b\right) }}{\mu },a_{0} = 0,\ a_{1} = \frac{\sqrt{cb\mu }}{c} \end{equation*}$

Substituting in Eq 2.3, we get,

$\begin{equation} u\left( \eta \right) = a_{1}\left( \mathrm{exp}\left( -\varphi \left( \eta \right) \right) \right) , \end{equation}$

(2.4)

Substituting all the five classes in Eq 2.4, we get the solutions.

Class 1: When, ${\lambda }^{2}-4\mu > 0$ and $\mu \neq 0$ ,

$\begin{equation*} v_{1} = -\frac{1}{2}\frac{\sqrt{cb\mu }\ tanh\left( \frac{1}{2}\left( \frac{1}{ 2}\frac{\sqrt{2}\sqrt{-\mu \left( 2a\mu +b\right) }t^{\alpha }}{\mu \Gamma \mathrm{(}\alpha +1)}+x\right) \sqrt{-4\mu }\right) \sqrt{-4\mu }}{c\mu } \end{equation*}$

Class 2: When, ${\lambda }^{2}-4\mu < 0$ and $\mu \neq 0,$

$\begin{equation} v_{2} = \frac{1}{2}\frac{\sqrt{cb\mu }\ tan\left( \frac{1}{2}\left( \frac{1}{2} \frac{\sqrt{2}\sqrt{-\mu \left( 2a\mu +b\right) }t^{\alpha }}{\mu \Gamma \mathrm{(}\alpha +1)}+x\right) \sqrt{4}\sqrt{\mu }\right) \sqrt{4}}{c\sqrt{ \mu }} \end{equation}$

(2.5)

Solution set 2:

$\begin{equation*} \left\{ \lambda = 0,C = C,\ V = \frac{1}{2}\frac{\sqrt{2}\sqrt{-\mu (2a\mu +b)}}{ \mu },a_{0} = 0,\ a_{1} = -\frac{\sqrt{cb\mu }}{c}\right\} \end{equation*}$

Substituting in Eq 2.3, we get,

$\begin{equation} u\left( \eta \right) = a_{1}\left( \mathrm{exp}\left( -\varphi \left( \eta \right) \right) \right) , \end{equation}$

(2.6)

Substituting equations 1.13 and 1.14 in Eq 2.6, we get the solutions.

Class 1: When, $\mu \neq 0$ and ${\lambda }^{2}-4\mu > 0,$

$\begin{equation} v_{3} = \frac{1}{2}\frac{\sqrt{cb\mu }\ tanh\left( \frac{1}{2}\left( -\frac{1}{ 2}\frac{\sqrt{2}\sqrt{-\mu \left( 2a\mu +b\right) }t^{\alpha }}{\mu \Gamma \mathrm{(}\alpha +1)}+x\right) \sqrt{-4\mu }\right) \sqrt{-4\mu }}{c\mu } \end{equation}$

(2.7)

Class 2: When, $\mu \neq 0$ and ${\lambda }^{2}-4\mu < 0,$

$\begin{equation} v_{4} = -\frac{1}{2}\frac{\sqrt{cb\mu }\ tan\left( \frac{1}{2}\left( -\frac{1}{ 2}\frac{\sqrt{2}\sqrt{-\mu \left( 2a\mu +b\right) }t^{\alpha }}{\mu \Gamma \mathrm{(}\alpha +1)}+x\right) \sqrt{4}\sqrt{\mu }\right) \sqrt{4}}{c\sqrt{ \mu }} \end{equation}$

(2.8)

Solution Set 3

$\begin{equation*} \left\{ \begin{array}{c} \lambda = \lambda ,\ C = C,\ V = \frac{\sqrt{-\left( {\lambda }^{2}-4\mu \right) \left( a{\lambda }^{2}-4a\mu -2b\right) }}{{\lambda }^{2}-4\mu },\ a_{0} = \frac{b\lambda }{\sqrt{c\left( {\lambda }^{2}-4\mu \right) b}},\ \\ a_{1} = -\frac{2\sqrt{-c\left( {\lambda }^{2}-4\mu \right) b}\ \mu }{c({ \lambda }^{2}-4\mu )} \end{array} \right\} \end{equation*}$

Substituting in Eq 2.3, we get,

$\begin{equation} u\left( \eta \right) = a_{1}\left( \mathrm{exp}\left( -\varphi \left( \eta \right) \right) \right) +a_{0}, \end{equation}$

(2.9)

Substituting equations 1.13 to 1.15 in Eq 2.9, we get the solutions.

Class 1: When, $\mu \neq 0$ and $-4\mu +{\lambda }^{2} > 0$ ,

$\begin{eqnarray} v_{5} & = &-\frac{b\lambda }{\sqrt{c\left( {\lambda }^{2}-4\mu \right) b}}-\frac{1}{c({\lambda }^{2}-4\mu )} \\ &&\left( \sqrt{c\left( {\lambda }^{2}-4\mu \right) b}\left( \begin{array}{c} -\lambda -tanh \\ \left( \frac{1}{2}\left( -\frac{\sqrt{\left( {\lambda }^{2}-4\mu \right) \left( a{\lambda }^{2}-4a\mu -2b\right) }t^{\alpha }}{{(\lambda }^{2}-4\mu )\Gamma \mathrm{(}\alpha +1)}+x\right) \sqrt{{\lambda }^{2}-4\mu }\right) \\ \sqrt{{\lambda }^{2}-4\mu } \end{array} \right) \right) \end{eqnarray}$

(2.10)

Class 2: When, $\mu \neq 0$ and $-4\mu +{\lambda }^{2} < 0$ ,

$\begin{eqnarray} v_{6} & = &-\frac{b\lambda }{\sqrt{c\left( {\lambda }^{2}-4\mu \right) b}} \\ &&-\left( \begin{array}{c} \sqrt{c\left( {\lambda }^{2}-4\mu \right) b} \\ \left( \begin{array}{c} -\lambda +tan \\ \left( \frac{1}{2}\left( -\frac{\sqrt{\left( {\lambda }^{2}-4\mu \right) \left( a{\lambda }^{2}-4a\mu -2b\right) }t^{\alpha }}{{(\lambda }^{2}-4\mu )\Gamma \mathrm{(}\alpha +1)}+x\right) \sqrt{-{\lambda }^{2}+4\mu }\right) \\ \sqrt{-{\lambda }^{2}+4\mu } \end{array} \right) \end{array} \right) \end{eqnarray}$

(2.11)

Class 3: When ${\lambda }\neq 0$ , $-4\mu +{\lambda } ^{2} > 0$ and $\mu = 0$ ,

$\begin{eqnarray} v_{7} & = &-\frac{b\lambda }{\sqrt{c\left( {\lambda }^{2}-4\mu \right) b}} \\ &&-\frac{2\sqrt{c\left( {\lambda }^{2}-4\mu \right) b}\mu \left( e^{\lambda }\left( -\frac{\sqrt{\left( {\lambda }^{2}-4\mu \right) \left( a{\lambda } ^{2}-4a\mu -2b\right) }t^{\alpha }}{({\lambda }^{2}-4\mu )\Gamma \mathrm{(} \alpha +1)}+x\right) -1\right) }{c\left( {\lambda }^{2}-4\mu \right) \lambda } \end{eqnarray}$

(2.12)

Class 4: When, ${\lambda }\neq 0$ , $\mu = 0$ and $-4\mu +{ \lambda }^{2} = 0,$

$\begin{eqnarray} v_{8} & = &-\frac{b\lambda }{\sqrt{c\left( {\lambda }^{2}-4\mu \right) b}} \\ &&+\frac{2\sqrt{c\left( {\lambda }^{2}-4\mu \right) b}\mu \left( 2\lambda \left( -\frac{\sqrt{\left( {\lambda }^{2}-4\mu \right) \left( a{\lambda } ^{2}-4a\mu -2b\right) }t^{\alpha }}{({\lambda }^{2}-4\mu )\Gamma \mathrm{(} \alpha +1)}+x\right) +2\right) }{c\left( {\lambda }^{2}-4\mu \right) { \lambda }^{2}\left( -\frac{\sqrt{\left( {\lambda }^{2}-4\mu \right) \left( a{ \lambda }^{2}-4a\mu -2b\right) }t^{\alpha }}{({\lambda }^{2}-4\mu )\Gamma \mathrm{(}\alpha +1)}+x\right) } \end{eqnarray}$

(2.13)

Case 5: When, ${\lambda } = 0$ , $\mu = 0$ and $-4\mu +{\lambda }^{2} = 0,$

$\begin{eqnarray} v_{9} & = &\frac{b\lambda }{\sqrt{-c\left( {\lambda }^{2}-4\mu \right) b}}- \frac{2}{c({\lambda }^{2}-4\mu )} \\ &&\left( \sqrt{-c\left( {\lambda }^{2}-4\mu \right) b}\ \mu \left( -\frac{ \sqrt{-\left( {\lambda }^{2}-4\mu \right) \left( a{\lambda }^{2}-4a\mu -2b\right) }t^{\alpha }}{({\lambda }^{2}-4\mu )\Gamma \mathrm{(}\alpha +1)} +x\right) \right) \end{eqnarray}$

(2.14)

3. Graphical demonstration

Physical interpretation

With some free parameters the proposed technique provides solitary wave solutions. By setting the specific parameters we have explained the miscellaneous wave solutions. In this study, we would explain the physical interpretation of the solutions for reaction-diffusion equation taking solution $v_{\mathrm{1}}$ for $\mu = 20 \; \lambda = 1 \; b = 11 \; a = 10 \; a_{1} = 12 \; c = -11 \; \alpha = 1,$ shows the solitary wave solution in . shows Soliton wave solution with paremeters $\mu = 910 \; \lambda = 1 \; b = 11 \; a = 10 \; a_{1} = 12 \; c = -11 \; \alpha = .5.$ Figures 3, and interprets the singular kink solution of $v_{ \mathrm{3}}, \; v_{ \mathrm{4}}$ , $v_{ \mathrm{7}}$ for $\mu = .010 \; \lambda = 1 \; b = .11 \; a = 10 \; a_{1} =$ 12 $c = -11 \; \alpha = .1, \; \mu = .0010 \; \lambda = -991 \; b = 11 \; a = 10 \; a_{1} = 102 \; c = -1 \; \alpha$ = .1, $\mu = 20 \; \lambda = 1 \; b = 11 \; a = 10 \; a_{1} = 12 \; c = -11 \; \alpha = .25.$ Finally kink wave results have been obtained from $\mathrm{v}_{\mathrm{5}}, \; \mathrm{v}_{\mathrm{6}},$ $\mathrm{v}_{\mathrm{8}}$ by setting the parameters, $\mu = 3.0 \; \lambda = 1 \; b = 11 \; a = 10 \; a_{1} =$ 12 $c = -11 \; \alpha = 0.001, \; \mu = 2.0 \; \lambda = 1 \; b = 11 \; a = 10 \; a_{1} = 12 \; c =$ -11 $\alpha = 0.001, \; \mu = 20 \; \lambda = 1 \; b = 11 \; a = 10 \; a_{1} = 12 \; c = -11 \; \alpha = 0.01,$ which is presented in Figures 5, 6 and 8. The solutions gained in this article have been checked by putting them back into the original equation and found correct. From the above obtained results we have many potential applications in fluid mechanics, quantum field theory, plasma physics and nonlinear optics.

Figure 1. Solitary wave solusion

$\nu_1(\eta)$ .

DownLoad: Full-Size Img PowerPoint

Figure 2. Soliton wave solusion

$\nu_2(\eta)$ .

DownLoad: Full-Size Img PowerPoint

Figure 3. Singular kink wave solusion

$\nu_3(\eta)$ .

DownLoad: Full-Size Img PowerPoint

Figure 4. Singular kink wave solusion

$\nu_4(\eta)$ .

DownLoad: Full-Size Img PowerPoint

Figure 5. kink wave solusion

$\nu_5(\eta)$ .

DownLoad: Full-Size Img PowerPoint

Figure 6. kink wave solusion

$\nu_6(\eta)$ .

DownLoad: Full-Size Img PowerPoint

Figure 7. Singular kink wave solusion

$\nu_7(\eta)$ .

DownLoad: Full-Size Img PowerPoint

Figure 8. kink wave solusion

$\nu_8(\eta)$ .

DownLoad: Full-Size Img PowerPoint

When

$\begin{equation*} \mu = 20,\ \lambda = 1,\ b = 11,\ a = 10,a_{1} = 12,c = -11,\alpha = 1 \end{equation*}$

When

$\begin{equation*} \mu = 910,\ \lambda = 1,\ b = 11,\ a = 10,a_{1} = 12,c = -11,\alpha = 1 \end{equation*}$

When

$\begin{equation*} \mu = .010,\ \lambda = 1,\ b = .11,\ a = 10,a_{1} = 12,c = -11,\alpha = .1 \end{equation*}$

When

$\begin{equation*} \mu = .0010,\ \lambda = -991,\ b = 11,\ a = 10,a_{1} = 102,c = -11,\alpha = .1 \end{equation*}$

When

$\begin{equation*} \mu = 3.0,\ \lambda = 1,\ b = 11,\ a = 10,a_{1} = 12,c = -11,\alpha = .001 \end{equation*}$

When

$\begin{equation*} \mu = 2.0,\ \lambda = 1,\ b = 11,\ a = 10,a_{1} = 12,c = -11,\alpha = .001 \end{equation*}$

When

$\begin{equation*} \mu = 20,\ \lambda = 1,\ b = 11,\ a = 10,a_{1} = 12,c = -11,\alpha = .25 \end{equation*}$

When

$\begin{equation*} \mu = 20,\ \lambda = 1,\ b = 11,\ a = 10,a_{1} = 12,c = -11,\alpha = .01 \end{equation*}$

4. Results and discussions

If we set $b = \beta, \; c = \gamma, \; \mu = r$ and $\eta = \xi$ in the obtaining solution $v_{2}$ and $v_{3}$ in this article is equal to $u_{5}$ and $u_{1}$ for case 5 respectively found in ^[31] see Table 1.

Table 1. Comparing the results of ^[31], with our results.

Attained Results	^[31] results
$\left(i\right)$ If we set $b=\beta, \;c=\gamma, \;\mu =r$ and $\eta =\xi$ then our solution $v_{3}$ becomes $v_{3}=\sqrt{\frac{-\beta }{\gamma }}\tanh \;\left(\sqrt{-r}\xi \right)$	$\left(i\right)$ The solution $u_{1}$ is as $u_{1}=\sqrt{\frac{-\beta }{\gamma }}\tanh \left(\sqrt{-r}\xi \right)$
$\left(ii\right)$ If we set $b=\beta, \;c=\gamma, \;\mu =r$ and $\eta=\xi$ then our solution $v_{2}$ becomes $v_{2}=\sqrt{\frac{\beta }{\gamma }}\tan \;\left(\sqrt{r}\xi \right)$	$\left(ii\right)$ The solution $u_{5}$ is as $u_{5}=\sqrt{\frac{\beta }{\gamma }}\tan \;\left(\sqrt{r}\xi \right)$

| Show Table

DownLoad: CSV

If we set $b = \beta, \; c = \gamma$ and $\eta = \xi,$ in the obtaining solution $v_{2}$ and $v_{4}$ in this article is equal to $u$ and $u$ for $c_{1}\neq 0, \; c_{2} = 0, \; \lambda = 0$ and $\mu > 0$ in ^[23] see Table 2.

Table 2. Comparing the results of ^[23], with our results.

Attained Results	^[23] results
$\left(i\right)$ If we set $b=\beta, \;c=\gamma$ and $\eta =\xi$ then our solution $v_{2}$ becomes $v_{2}=\sqrt{\frac{\beta }{\gamma }}\tan \;\left(\sqrt{\mu }\xi \right)$	$\left(i\right)$ The solution $u$ is as $u\left(\xi \right) =\pm \sqrt{\frac{\beta }{\gamma }}\tan \;\left(\sqrt{\mu }\xi \right)$
$\left(ii\right)$ If we set $b=\beta, \;c=\gamma$ and $\eta =\xi$ then our solution $v_{3}$ becomes $v_{3}=-\sqrt{\frac{\beta }{\gamma }}\tan \;\left(\sqrt{\mu }\xi \right)$	$\left(ii\right)$ The solution $u$ is as $u\left(\xi \right) =\pm \sqrt{\frac{\beta }{\gamma }}\tan \;\left(\sqrt{\mu }\xi \right)$

| Show Table

DownLoad: CSV

If we set $b = \beta, \; c = \gamma$ in the obtaining solution $v_{2}$ and $v_{4}$ are equal to $v_{9}$ and $v_{9}$ for $\lambda = 0$ and $\mu$ is positive in $v_{9}$ found in ^[32] see Table 3.

Table 3. Comparing the results of ^[32], with our results.

Attained Results	^[32] results
$\left(i\right)$ If we set $b=\beta, \;c=\gamma$ then our solution $v_{2}$ becomes $v_{2}=\sqrt{\frac{\beta }{\gamma }}\tan \;\left(\sqrt{\mu }\eta \right)$	$\left(i\right)$ The solution $u$ is as $u\left(\xi \right) =\pm \sqrt{\frac{\beta }{\gamma }}\tan \;\left(\sqrt{\mu }\eta \right)$
$\left(ii\right)$ If we set $b=\beta, \;c=\gamma$ then our solution $v_{4}$ becomes $v_{3}=-\sqrt{\frac{\beta }{\gamma }}\tan \;\left(\sqrt{\mu }\eta \right)$	$\left(ii\right)$ The solution $u$ is as $u\left(\xi \right) =\pm \sqrt{\frac{\beta }{\gamma }}\tan \;\left(\sqrt{\mu }\eta \right)$

| Show Table

DownLoad: CSV

If we set $b = -\beta, \; c = \gamma$ and $\eta = \xi$ in the obtaining solution $v_{2}$ and $v_{4}$ in this article are equal to $u_{2}$ and $u_{2}$ for $k > 0, \beta > 0$ and our $v_{1}$ and $v_{3}$ are equal to $u_{4}$ and $u_{4}$ for $k < 0, \beta < 0$ respectively founded in ^[30] see Table 4.

Table 4. Comparing the results of ^[30], with our results.

Attained Results	^[30] results
$\left(i\right)$ If we set $b=-\beta, \;c=\gamma, \;\mu =k$ and $\eta=\xi$ then our solution $v_{2}$ becomes $v_{2}=\sqrt{\frac{-\beta }{\gamma }}\tan \;\left(\sqrt{k}\xi \right)$	$\left(i\right)$ The solution $u_{2}$ is as $u_{1}=\pm \sqrt{\frac{-\beta }{\gamma }}\tan \;\left(\sqrt{k}\xi \right)$
$\left(ii\right)$ If we set $b=-\beta, \;c=\gamma, \;\mu =k$ and $\eta=\xi$ then our solution $v_{4}$ becomes $v_{2}=-\sqrt{\frac{\beta }{\gamma }}\tan \;\left(\sqrt{k}\xi \right)$	$\left(ii\right)$ The solution $u_{2}$ is as $u_{2}=\pm \sqrt{\frac{\beta }{\gamma }}\tan \;\left(\sqrt{k}\xi \right)$
$\left(iii\right)$ If we set $b=-\beta, \;c=\gamma, \;\mu =-k$ and $\eta =\xi$ then our solution $v_{1}$ becomes $v_{1}=-\sqrt{\frac{-\beta }{\gamma }}\tanh \;\left(\sqrt{-k}\xi \right)$	$\left(iii\right)$ The solution $u_{4}$ is as $u_{4}=\pm \sqrt{\frac{-\beta }{\gamma }}\tanh \;\left(\sqrt{-k}\xi \right)$
$\left(iv\right)$ If we set $b=\beta, \;c=\gamma, \;\mu =-k$ and $\eta=\xi$ then our solution $v_{4}$ becomes $v_{4}=-\sqrt{\frac{\beta }{\gamma }}\tanh \;\left(\sqrt{-k}\xi \right)$	$\left(iv\right)$ The solution $u_{4}$ is as $u_{4}=\pm \sqrt{\frac{\beta }{\gamma }}\tanh \;\left(\sqrt{-k}\xi \right)$

| Show Table

DownLoad: CSV

5. Conclusions

In the current paper, we explore that the proposed method is effective and capable to find exact solutions of reaction-diffusion equation. The obtained solutions indicate that the suggested method is direct, constructive and simple. The proposed technique can be implemented to the other NLPDEs of fractional order to establish new reliable solutions. The exact solutions are different and new along with different values of parameters. The reduction in the magnitude of computational part and the consistency of the technique give a broader applicability to the technique. The reaction diffusion equation has a dynamic role in dissipative dynamical systems as studied by various biologists, scientists and engineers. This model has found a number of applications in biology, physics (neutron diffusi0n theory), ecology and chemistry. It has also been claimed that reaction-diffusion processes have crucial basis for procedures associated to morphogenesis in biology and may even be connected to skin pigmentation and animal coats. Other applications of this model contain spread of epidemics, ecological invasions, wound healing and tumors growth. Another aim for the consideration in reaction-diffusion systems is that although they are nonlinear partial differential equations, there are often visions for an analytical treatment. My main contribution is programming and comparisons.

Conflict of interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

References

[1]	D. Bernoulli, Essai d'une nouvelle analyse de la mortalité causée par la petite vérole, et des avantages de l'Inoculation pour la prévenir, Hist. Acad. R. Sci. Paris, 1760, 1–45.
[2]	N. T. Bailey, The mathematical theory of infectious diseases and its applications, London: Charles Griffin & Company Limited, 1975.
[3]	R. M. Anderson, R. M. May, Infectious diseases of humans. Dynamics and control, Oxford: Oxford University Press, 1991. https://doi.org/10.1002/hep.1840150131
[4]	O. Diekmann, H. Heesterbeek, T. Britton, Mathematical tools for understanding infectious disease dynamics, Princeton: Princeton University Press, 2013.
[5]	Y. Cai, Y. Kang, W. Wang, A stochastic SIRS epidemic model with nonlinear incidence rate, Appl. Math. Comput., 305 (2017), 221–240. https://doi.org/10.1016/j.amc.2017.02.003 doi: 10.1016/j.amc.2017.02.003
[6]	X. Bardina, M. Ferrante, C. Rovira, A stochastic epidemic model of COVID-19 disease, AIMS Math., 5 (2020), 7661–7677. https://doi.org/10.3934/math.2020490 doi: 10.3934/math.2020490
[7]	F. Flandoli, E. La Fauci, M. Riva, Individual-based Markov model of virus diffusion: Comparison with COVID-19 incubation period, serial interval and regional time series, Math. Mod. Meth. Appl. Sci., 31 (2021), 907–939. https://doi.org/10.1142/S0218202521500226 doi: 10.1142/S0218202521500226
[8]	A. McKendrick, Application of mathematics to medical problems, Proc. Edinburgh Math. Soc., 14 (1926), 98–130. https://doi.org/10.1017/S0013091500034428 doi: 10.1017/S0013091500034428
[9]	W. Kermack, A. McKendrick, A contribution to the mathematical theory of epidemics, Proc. Roy. Soc. Lond., 115 (1927), 700–721. https://doi.org/10.1098/rspa.1927.0118 doi: 10.1098/rspa.1927.0118
[10]	H. Abbey, An examination of the Reed-Frost theory of epidemics, Hum. Biol., 24 (1952), 201–233.
[11]	L. J. S. Allen, An introduction to stochastic epidemic models, In: Mathematical Epidemiology, Berlin, Heidelberg: Springer-Verlag, 2008, 81–130. https://doi.org/10.1007/978-3-540-78911-6
[12]	H. C. Tuckwell, R. J. Williams, Some properties of a simple stochastic epidemic model of SIR type, Math. Biosci., 208 (2007), 76–97. https://doi.org/10.1016/j.mbs.2006.09.018 doi: 10.1016/j.mbs.2006.09.018
[13]	M. Ferrante, E. Ferraris, C. Rovira, On a stochastic epidemic SEIHR model and its diffusion approximation, Test, 25 (2016), 482–502. https://doi.org/10.1007/s11749-015-0465-z doi: 10.1007/s11749-015-0465-z
[14]	A. Gray, D. Greenhalgh, L. Hu, X. Mao, J. Pan, A stochastic differential equation SIS epidemic model, SIAM J. Appl. Math., 71 (2011), 876–902. https://doi.org/10.1137/10081856X doi: 10.1137/10081856X
[15]	T. Tsutsui, N. Minamib, M. Koiwai, T. Hamaokaa, I. Yamanea, K. Shimura, A stochastic-modeling evaluation of the foot-and-mouth-disease survey conducted after the outbreak in Miyazaki, Japan in 2000, Prev. Vet. Med., 61 (2003), 45–58. https://doi.org/10.1016/s0167-5877(03)00160-0 doi: 10.1016/s0167-5877(03)00160-0
[16]	S. Z. Huang, A new SEIR epidemic model with applications to the theory of eradication and control of diseases, and to the calculation of R0, Math. Biosci., 215 (2008), 84–104. https://doi.org/10.1016/j.mbs.2008.06.005 doi: 10.1016/j.mbs.2008.06.005
[17]	S. E. A. Mohammed, Stochastic differential systems with memory: Theory, examples and applications, In: Stochastic Analysis and related topics Ⅵ, Boston: Birkhaüser, 1998. https://doi.org/10.1007/978-1-4612-2022-0_1
[18]	H. Wahlstrom, L. Englund, T. Carpenter, U. Emanuelson, A. Engvall, I. Vagsholm, A Reed-Frost model of the spread of tuberculosis within seven Swedish extensive farmed fallow deer herds, Prev. Vet. Med., 35 (1998), 181–193. https://doi.org/10.1016/s0167-5877(98)00061-0 doi: 10.1016/s0167-5877(98)00061-0
[19]	J. Ng, E. J. Orav, A generalized chain-binomial model with application to HIV infection, Math. Biosci., 101 (1990), 99–119. https://doi.org/10.1016/0025-5564(90)90104-7 doi: 10.1016/0025-5564(90)90104-7
[20]	O. C. Ibe, Markov processes for stochastic modeling, Elsevier, 2013.

Reader Comments

Your name:*

Email:*
© 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)