Synchronization of generalized fractional complex networks with partial subchannel losses

Changping Dai; Weiyuan Ma; Ling Guo; Changping Dai; Weiyuan Ma; Ling Guo

doi:10.3934/math.2024344

AIMS Mathematics

2024, Volume 9, Issue 3: 7063-7083. doi: 10.3934/math.2024344

Previous Article Next Article

Research article Special Issues

Synchronization of generalized fractional complex networks with partial subchannel losses

1.
School of Mathematics and Computer Science, Northwest Minzu University, Lanzhou 730000, China
2.
College of Electrical Engineering, Northwest Minzu University, Lanzhou 730000, China

Received: 24 November 2023 Revised: 26 January 2024 Accepted: 30 January 2024 Published: 19 February 2024
MSC : 26A33, 34K24, 34K37, 65R10, 68R10

This article focuses on the synchronization problem for two classes of complex networks with subchannel losses and generalized fractional derivatives. Initially, a new stability theorem for generalized fractional nonlinear system is formulated using the properties of generalized fractional calculus and the generalized Laplace transform. This result is also true for classical fractional cases. Subsequently, synchronization criteria for the generalized fractional complex networks are attained by the proposed stability theorem and the state layered method. Lastly, two numerical examples with some new kernel functions are given to validate the synchronization results.

Keywords:

Citation: Changping Dai, Weiyuan Ma, Ling Guo. Synchronization of generalized fractional complex networks with partial subchannel losses[J]. AIMS Mathematics, 2024, 9(3): 7063-7083. doi: 10.3934/math.2024344

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]	H. Jeong, B. Tombor, R. Albert, Z. N. Oltvai, The large-scale organization of metabolic networks, Nature, 407 (2000), 651–654. https://doi.org/10.1038/35036627 doi: 10.1038/35036627
[2]	R. Albert, H. Jeong, A. L. Barabsi, Diameter of the world-wide web, Nature, 401 (1999), 130–131. https://doi.org/10.1038/43601 doi: 10.1038/43601
[3]	D. Knoke, S. Yang, Social network analysis, London: SAGE Publications, 2020. https://doi.org/10.4135/9781506389332
[4]	D. Lohr, P. Venkov, J. Zlatanova, Transcriptional regulation in the yeast GAL gene family: A complex genetic network, Faseb. J., 9 (1995), 777–787. https://doi.org/10.1096/fasebj.9.9.7601342 doi: 10.1096/fasebj.9.9.7601342
[5]	S. Strogatz, Exploring complex network, Nature, 410 (2001), 268–276. https://doi.org/10.1038/35065725
[6]	W. Y. Ma, Z. M. Li, N. R. Ma, Synchronization of discrete fractional-order complex networks with and without unknown topology, Chaos, 32 (2022), 013112. https://doi.org/10.1063/5.0072207 doi: 10.1063/5.0072207
[7]	H. Zhang, T. Ma, G. B. Huang, Z. Wang, Robust global exponential synchronization of uncertain chaotic delayed neural networks via dual-stage impulsive control, IEEE Trans. Syst. Man. Cybern. B, 40 (2009), 831–844. https://doi.org/10.1109/TSMCB.2009.2030506 doi: 10.1109/TSMCB.2009.2030506
[8]	T. Yang, Y. Niu, J. Yu, Clock synchronization in wireless sensor networks based on bayesian estimation, IEEE Access, 8 (2020), 69683–69694. https://doi.org/10.1109/ACCESS.2020.2984785 doi: 10.1109/ACCESS.2020.2984785
[9]	H. Zhang, D. Liu, Y. Luo, D. Wang, Adaptive dynamic programming for control-algorithms and stability, London: Springer, 2013. https://doi.org/10.1007/978-1-4471-4757-2
[10]	M. Wu, N. Xiong, A. V. Vasilakos, V. C. Leung, RNN-K: A reinforced newton method for consensus-based distributed optimization and control over multiagent systems, IEEE Trans. Cybern., 52 (2022), 4012–4026. https://doi.org/10.1109/TCYB.2020.3011819 doi: 10.1109/TCYB.2020.3011819
[11]	I. Podlubny, Fractional differential equations, Cambridge: Academic Press, 1999.
[12]	V. V. Uchaikin, Fractional derivatives for physicists and engineers, Berlin: Springer, 2013. https://doi.org/10.1007/978-3-642-33911-0
[13]	W. Y. Ma, N. R. Ma, C. P. Dai, Y. Q. Chen, X. Wang, Fractional modeling and optimal control strategies for mutated COVID- $19$ pandemic, Math. Method. Appl. Sci., 2023, 1–25. https://doi.org/10.1002/mma.9313
[14]	L. Ma, B. W. Wu, On the fractional Lyapunov exponent for Hadamard-type fractional differential system, Chaos, 33 (2023), 013117. https://doi.org/10.1063/5.0131661 doi: 10.1063/5.0131661
[15]	H. J. Li, J. D. Cao, Event-triggered group consensus for one-sided Lipschitz multi-agent systems with input saturation, Commun. Nonlinear Sci., 121 (2023), 107234. https://doi.org/10.1016/j.cnsns.2023.107234 doi: 10.1016/j.cnsns.2023.107234
[16]	C. P. Li, W. H. Deng, Remarks on fractional derivatives, Appl. Math. Comput., 187 (2007), 777–784. https://doi.org/10.1016/j.amc.2006.08.163 doi: 10.1016/j.amc.2006.08.163
[17]	N. G. N'Gbo, J. Tang, On the bounds of Lyapunov exponents for fractional differential systems with an exponential kernel, Int. J. Bifurcat. Chaos, 32 (2022), 2250188. https://doi.org/10.1142/S0218127422501887 doi: 10.1142/S0218127422501887
[18]	J. Hadamard, Essai sur l'étude des fonctions donnes par leur développement de Taylor, J. Math. Pure. Appl., 8 (1892), 101–186.
[19]	F. Jarad, D. Baleanu, A. Abdeljawad, Caputo-type modification of the Hadamard fractional derivatives, Adv. Differ. Equ., 2012 (2012), 1–8. https://doi.org/10.1186/1687-1847-2012-142 doi: 10.1186/1687-1847-2012-142
[20]	T. J. Osler, The fractional derivatives of a composite function, SIAM J. Math. Anal., 1 (1970), 288–293. https://doi.org/10.1137/0501026 doi: 10.1137/0501026
[21]	R. Almeida, A Caputo fractional derivative of a function with respect to another function, Commun. Nonlinear Sci., 44 (2017), 460–481. https://doi.org/10.1007/s11868-021-00421-y doi: 10.1007/s11868-021-00421-y
[22]	W. Ma, C. Dai, X. Li, X. Bao, On the kinetics of $\psi$ -fractional differential equations, Fract. Calc. Appl. Anal., 2023. https://doi.org/10.1007/s13540-023-00210-y
[23]	G. Mahmoud, M. Ahmed, T. Abed-Elhameed, Active control technique of fractional-order chaotic complex systems, Eur. Phys. J. Plus, 131 (2016), 1–11. https://doi.org/10.1140/epjp/i2016-16200-x doi: 10.1140/epjp/i2016-16200-x
[24]	W. Zheng, Y. Q. Chen, X. Wang, M. Lin, A neural network-based design method of the fractional order PID controller for a class of motion control systems, Asian J. Control, 24 (2022), 3378–3393. https://doi.org/10.1002/asjc.2727 doi: 10.1002/asjc.2727
[25]	X. Yang, J. Cao, Finite-time stochastic synchronization of complex networks, Appl. Math. Model., 34 (2010), 3631–3641. https://doi.org/10.1016/j.apm.2010.03.012 doi: 10.1016/j.apm.2010.03.012
[26]	L. Duan, J. Li, Fixed-time synchronization of fuzzy neutral-type BAM memristive inertial neural networks with proportional delays, Inf. Sci., 576 (2021), 522–541. https://doi.org/10.1016/j.ins.2021.06.093 doi: 10.1016/j.ins.2021.06.093
[27]	W. Zhang, C. Li, X. He, H. Li, Finite-time synchronization of complex networks with non-identical nodes and impulsive disturbances, Mod. Phys. Lett. B, 32 (2018), 1850002. https://doi.org/10.1142/S0217984918500021 doi: 10.1142/S0217984918500021
[28]	Y. Wang, X. He, T. Li, Asymptotic and pinning synchronization of fractional-order nonidentical complex dynamical networks with uncertain parameters, Fractal Fract., 7 (2023), 571. https://doi.org/10.3390/fractalfract7080571 doi: 10.3390/fractalfract7080571
[29]	P. F. Xia, S. L. Zhou, G. B. Giannakis, Adaptive MIMO-OFDM based on partial channel state information, IEEE Trans. Signal. Process., 52 (2004), 202–213. https://doi.org/10.1109/TSP.2003.819986 doi: 10.1109/TSP.2003.819986
[30]	C. Huang, D. W. C. Ho, J. Lu, Partial-information-based distributed filtering in two-targets tracking sensor networks, IEEE Trans. Circ. Syst. I, 59 (2012), 820–832. https://doi.org/10.1109/TCSI.2011.2169912 doi: 10.1109/TCSI.2011.2169912
[31]	Q. Wu, H. Zhang, L. Xu, Q. Yan, Finite-time synchronization of general complex dynamical networks, Asian J. Control, 17 (2015), 1643–1653. https://doi.org/10.1002/asjc.985 doi: 10.1002/asjc.985
[32]	C. Zhou, L. Zemanová, G. Zamora-Lopez, Structure-function relationship in complex brain networks expressed by hierarchical synchronization, New J. Phys., 9 (2007), 178. https://doi.org/10.1088/1367-2630/9/6/178 doi: 10.1088/1367-2630/9/6/178
[33]	L. Li, X. Liu, W. Huang, Event-based bipartite multi-agent consensus with partial information transmission and communication delays under antagonistic interactions, Sci. China Inf. Sci., 63 (2020), 150204. https://doi.org/10.1007/s11432-019-2693-x doi: 10.1007/s11432-019-2693-x
[34]	Y. Li, J. Zhang, J. Lu, J. Lou, Finite-time synchronization of complex networks with partial communication channels failure, Inf. Sci., 634 (2023), 539–549. https://doi.org/10.1016/j.ins.2023.03.077 doi: 10.1016/j.ins.2023.03.077
[35]	Q. Fan, G. C. Wu, H. Fu, A note on function space and boundedness of the general fractional integral in continuous time random walk, J. Nonlinear Math. Phys., 29 (2022), 95–102. https://doi.org/10.1007/s44198-021-00021-w doi: 10.1007/s44198-021-00021-w
[36]	R. Almeida, A. B. Malinowska, T. Odzijewicz, On systems of fractional differential equations with the $\psi$ -Caputo derivative and their applications, Math. Methods Appl. Sci., 44 (2021), 8026–8041. https://doi.org/10.1002/mma.5678 doi: 10.1002/mma.5678
[37]	F. Jarad, T. Abdeljawad, Generalized fractional derivatives and Laplace transform, Discrete Cont. Dyn. S, 13 (2020), 709–722. https://doi.org/10.3934/dcdss.2020039 doi: 10.3934/dcdss.2020039
[38]	A. Ahmadova, N. Mahmudov, Asymptotic stability analysis of Riemann-Liouville fractional stochastic neutral differential equations, Miskolc Math. Notes, 22 (2021), 503–520. https://doi.org/10.18514/MMN.2021.3600 doi: 10.18514/MMN.2021.3600
[39]	B. K. Lenka, S. N. Bora, Lyapunov stability theorems for $\psi$ -Caputo derivative systems, Fract. Calc. Appl. Anal., 26 (2023), 220–236. https://doi.org/10.1007/s13540-022-00114-3 doi: 10.1007/s13540-022-00114-3
[40]	S. Liu, W. Jiang, X. Li, X. F. Zhou, Lyapunov stability analysis of fractional nonlinear systems, Appl. Math. Lett., 51 (2016), 13–19. https://doi.org/10.1016/j.aml.2015.06.018 doi: 10.1016/j.aml.2015.06.018
[41]	W. Yu, G. Chen, J. Lü, On pinning synchronization of complex dynamical networks, Automatica, 45 (2009), 429–435. https://doi.org/10.1016/j.automatica.2008.07.016 doi: 10.1016/j.automatica.2008.07.016

Reader Comments

Your name:*

Email:*
© 2024 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)