The Global dynamics of a SIR model considering competitions among multiple strains in patchy environments

1.   Introduction

Many problems in the real world, such as in the area of physics, engineering, natural sciences, economics and so on, can be represented by many nonlinear partial differential equations. Therefore, it is of great significance to determine the exact solutions of such equations to better understand their physical properties. Generally speaking, no direct expression formulas are available for these solutions, even for simple nonlinear partial differential equations. For the past few decades, many effective techniques have been developed by researchers for searching explicit solutions including the Lie symmetry analysis method ^[1,2,3,4], the Darboux transformation method ^[5], the inverse scattering transform method ^[6], the Jacobi elliptic function expansion method ^[7], the Kudryashov method ^[8], the $(G'/G)-$expansion method ^[9,10], the $(G'/G^{2})-$expansion method ^[11], the Painlevé analysis method ^[12] and so on.

In the middle of the 19th century, Norwegian mathematician Sophus Lie (1842–1899) pioneered the Lie symmetry analysis. Lie inaugurated the theory of Lie symmetries and applied it to differential equations. These research results have formed an important branch of mathematics in the 20th century, which is known as a Lie group and Lie algebra theory. By utilizing the invariance of solutions under a symmetry group for a partial differential equation, one can obtain solutions of a PDE via solving a reduced differential equation with lesser independent variables. Unfortunately, Lie symmetry analysis involves in a large number of complex calculations. With the rapid development of computing software, the Lie symmetry analysis method has become the most effective and powerful technique for obtaining closed form solutions to nonlinear partial differential equations. Many exact solutions with an important physical significance, such as the similarity solutions and travelling wave solutions, can be obtained by using the Lie symmetry analysis method ^{[13,14,15,16]}.

During the investigation of differential equations, conservation laws play an important role and have been receiving increased attention ^[17]. On one hand, they can be used to judge whether a partial differential equation is complete integrable. On the other hand, conservation laws can also be used to verify the validity of numerical solution methods. Moreover, exact solutions of partial differential equations can be constructed by using conservation laws ^[18,19,20]. Therefore, the most important thing is how to derive conservation laws for a given differential equation. With the development of research, some methods for constructing conservation laws of equations have emerged. For example, Noether's theorem provides an efficient approach to construct conservation laws for systems with the Lagrangian formulation ^[21,22]. Furthermore, the Ibragimov theorem ^[23,24,25] and the multiplier method ^[26] have more wider applications in deriving conservation laws, as they can be used for arbitrary differential equations without the limitation of the Lagrangian formulation.

The investigations of shallow water waves and solitary have gained a great significance to describe characteristics of nonlinear wave phenomena. Stable dynamics and excitations of single- and double-hump solitons in the Kerr nonlinear media have been achieved in Ref. ^[27]. Bidirectional solitons and interaction solutions for a new integrable fifth-order nonlinear equation have been investigated in Ref. ^[28]. Nonlinear dynamic behaviors of the generalized (3+1)-dimensional KP equation have been studied in Ref. ^[29]. Ten years ago, Wazwaz ^[30] derived a (2+1)-dimensional Korteweg–de Vries 4 (KdV4) equation by using the recursion operator of the KdV equation, which reads

As a KdV type equation, Eq (1.1) can be used to model the shallow-water waves, surface, and internal waves. In Ref. ^[30], Wazwaz derived multiple soliton solutions by invoking Hirota's bilinear method and other travelling wave solutions by using hyperbolic functions and trigonometric function methods. The Bäcklund transformations and soliton solutions have been investigated in Ref. ^[31]. In Ref. ^[32], the authors derived breather wave solutions by using the extend homoclinic test method and some travelling wave solutions by using the $(G'/G^{2})-$expansion method.

In this work, we aim to investigate the (2+1)-dimensional KdV4 equation (1.1) by invoking the Lie symmetry analysis method and consider its conservation laws. In Section 2, we perform the Lie symmetry analysis method to Eq (1.1) and present two symmetry reductions. In Section 3, we present solutions of symmetry reductions obtained in Section 2 by different methods. In Section 4, we present the conservation laws of the equation by invoking the multiplier method and Noether's theorem, respectively. A few concluding remarks will be given in the final section.

2.   Lie symmetry analysis

First, we employ the Lie classical method on Eq (1.1) to find its symmetries. In this end, we consider a parameter $(\epsilon)$ Lie symmetry group of infinitesimal transformations as:

in which $\varphi, \omega, \tau$ and $\eta$ are infinitesimal generators. The corresponding vector field on transformation (2.1) is given by

Noticeably, symmetries of Eq (1.1) can be derived according to the symmetry conditions

where $pr^{(4)}R$ is the fourth prolongation formula of Eq (2.2), which is given by

By utilizing the invariance condition (2.3), we have

As the original equation is $v_{xy} + v_{xxxt} + v_{xxxx} + 3(v_{x}^{2})_{x} + 4v_{x}v_{xt} + 2v_{xx}v_{t} = 0,$ when substituting differential prolongations of $\eta$ into (2.4) and separating on different derivatives of $v$, we set the coefficients of $v_{xy}$, $v_{xxxt}$, $v_{xxxx}$, $3(v_{x}^{2})_{x}$, $4v_{x}v_{xt}$, $2v_{xx}v_{t}$ to be equal and other derivatives of $v$ to be zero. Thereafter, we obtain the following homogeneous PDEs:

By calculating the above PDEs, the symmetries of the KdV4 equation (1.1) can be given as

According to the results of the infinitesimal transformation, the corresponding single parameter transformation groups are:

As a result, if $v = \rho(x, y, t)$ is a solution of the KdV4 equation (1.1), so are

In what follows, we implement symmetry reductions of the KdV4 equation (1.1) by employing its translation symmetries. We choose two combinations $k\sigma_{1} + \sigma_{2}$ and $2\sigma_{3} + \sigma_{4} - \sigma_{5}$ as examples to perform symmetry reductions.

(ⅰ) For symmetry $k\sigma_{1} + \sigma_{2}$, by taking $f(y) = 1, g(y) = 0$, then it gives three invariants

By substituting (2.5) into Eq (1.1), we have

By calculation, $\Gamma_{1} = (k-1)\frac{\partial}{\partial X} + k\frac{\partial}{\partial Y}$ and $\Gamma_{2} = \frac{\partial}{\partial X}$ are Lie point symmetries of Eq (2.6). Considering the symmetry $\Gamma = \Gamma_{1} + l\Gamma_{2}$, where $l$ being an arbitrary constant, it provides us with two invariants

By substituting (2.7) into Eq (2.6), we have

(ⅱ) For symmetry $-2\sigma_{4} + \sigma_{5}$, it admits three invariants

By substituting (2.9) into Eq (1.1), we have

By calculation, $\Gamma = -\frac{1}{2}\alpha\frac{\partial}{\partial \alpha} + \beta\frac{\partial}{\partial \beta} + (\frac{1}{16}\beta + \frac{1}{2}P)\frac{\partial}{\partial P}$ is a Lie point symmetrie of Eq (2.10), and here $\Gamma$ provides us with two invariants

By substituting (2.11) into Eq (2.10), we have

3.   Solutions of KdV4 equation

3.1.   Solutions of Eq (2.8) via direct integration method

Upon integrating Eq (2.8) twice, we get

where $m$ and $n$ are constants. Denote $F_{\chi}$ as $\Omega$, then Eq (3.1) can be written as

If the polynomial equation

has three roots $p_{1}, p_{2}, p_{3}$ and $p_{1} > p_{2} > p_{3}$, then Eq (3.2) can be rewritten as

whose solution is

in which $\mathrm{cn}$ denotes cosine elliptic function. As $\Omega = F_{\chi}$, integrating (3.4) once and going back to variables $x, y$ and $t$, we acquire a periodic solution of the KdV4 equation

where $\chi = x- \frac{l+k-1}{k}y+(l-1)t$, and $C$ is a constant of integration. Comparing Eqs (3.2) and (3.3), we can easily find that the following relationship between $p_{1}$, $p_{2}$, $p_{3}$ and $k$, $l$ needs to satisfy

Furthermore

denotes the incomplete elliptic integral. By selecting $k = 2$, $l = 3$, $p_{1} = 5$, $p_{2} = 3$, $p_{3} = 1$, $C = 0$ in (3.7), we present the profiles of Eq (3.5) in Figure 1 with $t = 0$.

By analysis, we can see the physical characteristics in Figure 1. On one hand, the features of the elliptic function are shown in Figure 1(a), (c); on the other hand, it exhibits certain periodicities in Figure 1(b).

If the roots satisfy $p_{2} = p_{3}$, then the solution (3.5) can be reduced to the following form

When $p_{1} = 2h^{2}$($h$ is a constant) and $p_{2} = p_{3} = 0$, then the solution (3.5) can be reduced by

which is similar to the solution obtained in ^[30]. By selecting $h = 2$, $k = 3$, $l = 4$, $C = 2$ in (3.5), we present the profiles of Eq (3.7) in Figure 2 with $t = 1$. By observing Figure 2, we can see that Eq (3.7) is a kink-type hyperbolic function solution.

3.2.   Solutions of Eq (2.8) via the $(G^{'}/G)$ expansion method

In what follows, we investigate the solutions of Eq (2.8) by invoking the $(G^{'}/G)$ expansion method ^[9]. By balancing the highest order derivative $F_{\chi\chi\chi\chi}$ and the nonlinear term of the highest order $F_{\chi}F_{\chi\chi}$, we assume the solution of Eq (2.8) in this form

where $G(\chi)$ satisfies the following Riccati equation

Substituting Eq (3.8) into Eq (2.8), using Eq (3.9) for reduction, collecting all the powers of $\left(\frac{G^{'}(\chi)}{G(\chi)}\right)^{-1}$ and $\left(\frac{G^{'}(\chi)}{G(\chi)}\right)$, then equating all the obtained coefficients to zero, eleven algebraic equations can be obtained:

Solving above equations with Maple, two kinds of solutions are obtained:

Case 1:

Case 2:

As a result, six kinds of solutions to the KdV4 equation (1.1) can be obtained accordingly.

For the first case, when $\lambda^{2} - 4\mu > 0$, we get one hyperbolic function solution

where $\Lambda_{1} = \frac{\sqrt{\lambda^{2} - 4\mu}}{2}$, $\chi = x- \frac{l+k-1}{k}y+(l-1)t$ with $k = \frac{l +1}{l \, \lambda^{2}-4 l \mu +1}$ and $C_{1}, C_{2}$ are arbitrary constants. By selecting $\lambda = 3$, $\mu = 2$, $l = 3$, $c_{1} = \frac{1}{3}$, $c_{2} = \frac{1}{4}$, $a_{0} = 1$, we present the profiles of Eq (3.10) in Figure 3 with $t = 1$.

Noticeably, Eq (3.10) is the combination of hyperbolic sine function and hyperbolic cosine function. When choosing special parameters, its expression may have discontinuities, and multiple peaks and valleys appeared in Figure 3.

When $\lambda^{2} - 4\mu < 0$, we acquire one trigonometric function solution

where $\Lambda_{2} = \frac{\sqrt{4\mu - \lambda^{2}}}{2}$, $\chi = x- \frac{l+k-1}{k}y+(l-1)t$ with $k = \frac{l +1}{l \, \lambda^{2}-4 l \mu +1}$ and $C_{1}, C_{2}$ are arbitrary constants. By selecting $\lambda = 2$, $\mu = 2$, $l = 2$, $c_{1} = -1$, $c_{2} = 3$, $a_{0} = 0$, we present the profiles of Eq (3.11) in Figure 4 with $t = 1$.

Noticeably, Eq (3.11) is the combination of sine function and cosine function. When choosing special parameters, its denominator is similar to its numerator, so features of cotangent function appeared in Figure 4.

When $\lambda^{2} - 4\mu = 0$, we achieve one rational function solution

where $C_{1}, C_{2}$ are arbitrary constants. Noticeably, Eq (3.12) is a plane wave solution.

For the second case, when $\lambda^{2} - 4\mu > 0$, we get the other hyperbolic function solution

where $\Lambda_{1} = \frac{\sqrt{\lambda^{2} - 4\mu}}{2}$, $\chi = x- \frac{l+k-1}{k}y+(l-1)t$ with $k = \frac{l +1}{l \, \lambda^{2}-4 l \mu +1}$ and $C_{1}, C_{2}$ are arbitrary constants. By selecting $\lambda = -3$, $\mu = 2$, $l = 3$, $c_{1} = \frac{1}{4}$, $c_{2} = \frac{1}{3}$, $a_{0} = 2$, we present the profiles of Eq (3.13) in Figure 5 with $t = 0$.

Noticeably, Eq (3.13) is the combination of hyperbolic sine function and hyperbolic cosine function and is seeming like the inverse of Eq (3.10). When choosing special parameters, its expression may also have discontinuities, and multiple peaks and valleys appeared in Figure 5.

When $\lambda^{2} - 4\mu < 0$, we acquire the other trigonometric function solution

where $\Lambda_{2} = \frac{\sqrt{4\mu - \lambda^{2}}}{2}$, $\chi = x- \frac{l+k-1}{k}y+(l-1)t$ with $k = \frac{l +1}{l \, \lambda^{2}-4 l \mu +1}$ and $C_{1}, C_{2}$ are arbitrary constants. By selecting $\lambda = 2$, $\mu = 2$, $l = 2$, $c_{1} = -1$, $c_{2} = 3$, $a_{0} = 0$, we present the profiles of Eq (3.14) in Figure 6 with $t = 1$.

Noticeably, Eq (3.14) is the combination of sine function and cosine function and is seeming like the inverse of Eq (3.11). When choosing special parameters, its denominator is similar to its numerator, so features of cotangent function appeared in Figure 6.

When $\lambda^{2} - 4\mu = 0$, we achieve the other rational function solution

where $C_{1}, C_{2}$ are arbitrary constants. Evidently the denominator of Eq (3.15) is a plane, so Eq (3.15) is a curved surface.

3.3.   Solutions of Eq (2.12) via the power series technique

First, we suppose the solution of (2.12) is given by

where $Q_{0}, Q_{1}, \cdots, Q_{k}, \cdots$ are constants to be determined. By substituting the expressions of $Q_{\gamma}, Q_{\gamma\gamma},$ $Q_{\gamma\gamma\gamma},$ $Q_{\gamma\gamma\gamma\gamma}$ into Eq (2.12), we have

By equating the coefficients of the powers of $\gamma$ in Eq (3.16) to be zero, we get

Therefore, the power series solution of Eq (2.12) is

Hence, the particular solution of Eq (1.1) is given by

where $Q(\frac{(x-t)^{2}t}{y})$ is determined by Eq (3.17).

4.   Conservation laws

In this section, we derive the conservation laws of the KdV4 equation (1.1) by means of the multiplier method and Noether's theorem, respectively.

4.1.   Conservation laws by multiplier method

In the beginning, we invoke the multiplier method to determine conserved quantities of the KdV4 equation (1.1). In order to obtain the first order multipliers $R$, namely

These multipliers are derived by

in which

denotes the Euler operator and $D_{x}, D_{y}, D_{t}$ are total derivative operators. Expanding (4.1) and collecting on derivatives of $v$ from second order to fourth order, then equating their coefficients to be zero, we get twenty-four PDEs:

By solving above equations with Maple, we obtain

where $C_{1}$, $C_{2}$, $C_{3}$, $C_{4}$ are arbitrary constants and $g, f$ are functions of $y$. The conserved quantities of Eq (1.1) can be derived by employing the following divergence expression

with $C^{x}, C^{y}$ are spatial fluxes and $C^{t}$ is the conserved density. Therefore, we obtain the following conserved vectors $(C^{x}, C^{y}, C^{t})$ with regard to the six multipliers:

Case 1. For $R_{1} = tv_{x} + 2tv_{t} + xv_{x} + 6yv_{x} + 4yv_{y} + v - x,$ we get the corresponding conserved vector corresponding to $R_{1}$ as follows

Case 2. For $R_{2} = v_{x},$ we get the corresponding conserved vector corresponding to $R_{2}$ as follows

Case 3. For $R_{3} = v_{y},$ we get the corresponding conserved vector corresponding to $R_{3}$ as follows

Case 4. For $R_{4} = v_{t},$ we get the corresponding conserved vector corresponding to $R_{4}$ as follows

Case 5. For $R_{5} = g'(y)t - 2v_{x}g(y),$ we get the corresponding conserved vector corresponding to $R_{5}$ as follows

Case 6. For $R_{6} = f(y),$ we get the corresponding conserved vector corresponding to $R_{6}$ as follows

4.2.   Conservation laws by Noether's Theorem

Next, we construct the conservation laws of Eq (1.1) by utilizing the classical Noether's theorem ^[20]. The second order Lagrangian of equation (1.1) is

Therefore, the Noether symmetries

of the KdV4 equation (1.1) are established by applying the Lagrangian equation (4.2) on the following equation

where $R^{[2]}$ is the second prolongation of $R$ and $B^{x}$, $B^{y}$, $B^{t}$ are gauge functions. Expanding the above equation (4.3) and seeking solutions for the resulting system of partial differential equations, the following Noether symmetries as well as their corresponding gauge functions can be obtained:

By using the above Noether symmetries and invoking the following formula ^[21]

We obtain the corresponding conserved vectors, which are given by

5.   Remarks

In this paper, we performed the Lie symmetry analysis method to a (2+1)-dimensional KdV4 equation and derived Lie symmetries of this equation, then used these symmetries to perform symmetry reductions. By using translation symmetries, two fourth-order ordinary differential equations were obtained. Different methods were adopted in order to derive the solutions of the obtained differential equations. For one fourth-order ordinary differential equation, we presented its solutions by using the direct integration method and the $(G'/G)$-expansion method, respectively. Furthermore, their corresponding solutions were shown with the appropriate graphical representations. For the other fourth-order ordinary differential equation, we derived its solutions by using the power series technique. Finally, two kinds of conserved vectors of this equation are presented by invoking the multiplier method and Noether's theorem, respectively. Six multipliers were obtained from the multiplier method, and thus six local conservation laws for the KdV4 equation (1.1) were given, while five local conservation laws were obtained by Noether's theorem. It is necessary to point out that the multiplier method has more wider applications in deriving conservation laws than Noether's theorem. Although we successfully presented the power series solution for the KdV4 equation (1.1), we had difficulties in discussing the corresponding convergence analysis. Its convergence analysis is worth studying in our future work. For the Lie symmetry analysis method, many attentions have been concentrated on (1+1)-dimensional and (2+1)-dimensional differential equations. For (3+1)-dimensional differential equations, there are not many relevant documents. We will investigate the (3+1)-dimensional differential equations by using the Lie point symmetry analysis method in the future.

Acknowledgments

This work was supported by Henan Youth Backbone Teacher Project (No. 2019GGJS204) and the Natural Science Foundation of Henan Province (No.232300420361).

Conflict of interest

The author declares there is no conflict of interest.