Next Article in Journal
New Version of Fractional Pachpatte-Type Integral Inequalities via Coordinated ℏ-Convexity via Left and Right Order Relation
Next Article in Special Issue
A Two-Temperature Fractional DPL Thermoelasticity Model with an Exponential Rabotnov Kernel for a Flexible Cylinder with Changeable Properties
Previous Article in Journal
Topological Properties of Polymeric Networks Modelled by Generalized Sierpiński Graphs
Previous Article in Special Issue
Dynamic Compressive Mechanical Property Characteristics and Fractal Dimension Applications of Coal-Bearing Mudstone at Real-Time Temperatures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of the Time Fractional Higher-Dimensional Nonlinear Modified Equation of Wave Propagation

1
School of Mathematics and Statistics, Changshu Institute of Technology, Changshu 215500, China
2
Qin Institute of Mathematics, Shanghai Hanjing Centre for Science and Technology, Yexie Town, Songjiang District, Shanghai 201609, China
3
School of Mathematics and Statistics, SuZhou University, Suzhou 234000, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(3), 124; https://doi.org/10.3390/fractalfract8030124
Submission received: 31 January 2024 / Revised: 14 February 2024 / Accepted: 17 February 2024 / Published: 20 February 2024

Abstract

:
In this article, we analyzed the time fractional higher-dimensional nonlinear modified model of wave propagation, namely the (3 + 1)-dimensional Benjamin–Bona–Mahony-type equation. The fractional sense was defined by the classical Riemann–Liouville fractional derivative. We derived firstly the existence of symmetry of the time fractional higher-dimensional equation. Next, we constructed the one-dimensional optimal system to the time fractional higher-dimensional nonlinear modified model of wave propagation. Subsequently, it was reduced into the lower-dimensional fractional differential equation. Meanwhile, on the basis of the reduced equation, we obtained its similarity solution. Through a series of analyses of the time fractional high-dimensional model and the results of the above obtained, we can gain a further understanding of its essence.

1. Introduction

It is known that differential equations originated from the real natural life, such as the RLC circuit model [1], population model [2], realistic ecological model [3], KdV model [4], diffusion type model [5], etc. By analyzing these models, we can better understand, comprehend, and utilize them. Meanwhile, since the higher-dimensional nonlinear models contain more spatial dimensions x ,   y and z, they contain more information for the natural world. As a result, the nonlinear higher-dimensional models is one of the important topics in the field of mathematical physics. On the basis of the above idea, A.M. Wazwaz [6] derived five categories of Benjamin–Bona–Mahony-type equations of wave propagation and fluids, and obtained their exact solutions, such as the soliton, kink, and periodic solutions. At the same time, there is another type of fractional derivative that appears at the same time as integer order. We known that the fractional derivative is used to describe the complex nonlinear mechanism. There exist many different fractional derivative definitions, such as the Riemann–Liouville fractional derivative [7,8], Caputo fractional derivative [7,8], general fractional derivative [9], He’s fractional derivative [10], etc. They have been widely used in diffusion processes in solids, biological media, and fractal media [11]. Here, we considered one of the five categories of the following form:
D t α u + u x + u 2 u y u x z t = 0 , 0 < α 1 ,
where:
D t α u ( x , y , z , t ) = 1 Γ ( n α ) n t n 0 t u ( x , y , z , τ ) ( t τ ) α + 1 n d τ , n 1 < α < n , n N , n u ( x , y , z , t ) t n , α = n N ,
is the Riemann–Liouville fractional derivative. A.R. Seadawy et al. [12] obtained the periodic wave and kink wave of Equation (1). Liu et al. [13] studied the dynamical behavior of the traveling systems and deep neural network. S. Muhammad et al. [14] derived the diversity of soliton solutions of Equation (1). B. Muhammad et al. [15,16] found many solitary waves solutions by using the analytical method.
Due to the more detailed characterization of real life by fractional derivatives and the fact that high-dimensional nonlinear models contain more information, the main motivation of this paper is to explore the mechanism of water wave propagation by studying the high-dimensional model. The research results will fill the gap in the water wave theory in this regard. The aims of this article are to apply the symmetry scheme to investigate the time fractional higher-dimensional nonlinear modified equation of wave propagation (1). As a direct result, a series of new results were obtained. The fractional symmetry analysis scheme is an effective tool to deal with fractional differential equations, specifically the higher-dimensional nonlinear models [17,18,19,20,21,22,23,24]. For example, Adeyemo et al. [17] studied the time fractional (3 + 1)-dimensional generalized Zakharov-Kuznetsov equation type I, Liu et al. [18,19,20,21] researched the higher-dimensional KdV-type equation, dissipation Burgers equation, and diffusion equation, Sahoo et al. [22] discussed the (3 + 1)-dimensional time-fractional mKdV-ZK equation, and Zhuo et al. [23] analyzed the generalized (4 + 1)-dimensional time-fractional Fokas equation, Zhu et al. [24] considered the time-fractional (2 + 1)-dimensional Hirota-Satsuma-Ito equations. Therefore, we hope to find more novel results through this effective tool, laying the foundation for our deeper understanding of this model.
The compositions of this article as follows: the symmetries of Equation (1) through the symmetry scheme are derived in Section 2. In Section 3 and Section 4, we discuss, respectively, the one-dimensional optimal system and one-parameter Lie transformation group of Equation (1). With the help of the Erdélyi–Kober fractional operators, Equation (1) is reduced into the lower fractional differential equation in Section 5. Then, the similarity solution of the reduced equation is found. The conclusions and discussions of the entire text and future prospects are presented in Section 6.

2. Symmetry Scheme of Equation (1)

In the first place, the symmetries of the time fractional higher-dimensional nonlinear modified equation of wave propagation (1) were obtained. As a direct conclusion, we have Theorem 1.
Theorem 1. 
The time fractional higher-dimensional nonlinear modified equation of wave propagation (1) has the infinitesimal generators of the forms:
X 1 = 2 x x + 2 α y y 2 α z z + 2 α α u u + 2 α t t , X 2 = x , X 3 = y , X 4 = z .
Proof. 
It is assumed that Equation (1) is invariant under a one-parameter Lie transformation group as follows:
x * x + ε ξ x ( x , y , z , t , u ) + O ( ε 2 ) , y * y + ε ξ y ( x , y , z , t , u ) + O ( ε 2 ) , z * z + ε ξ z ( x , y , z , t , u ) + O ( ε 2 ) , t * t + ε ξ t ( x , y , z , t , u ) + O ( ε 2 ) , u * u + ε η u ( x , y , z , t , u ) + O ( ε 2 ) , D t α u * D t α u + ε η u α ( x , y , z , t , u ) + O ( ε 2 ) , u * x * u x + ε η u x ( x , y , z , t , u ) + O ( ε 2 ) ,
where ε 1 is an infinitesimal parameter.
In detail:
η u x = D x ( η u ξ x u x ξ y u y ξ z u z ξ t u t ) + ξ x u x x + ξ y u x y + ξ z u x z + ξ t u x t , η u y = D y ( η u ξ x u x ξ y u y ξ z u z ξ t u t ) + ξ x u x y + ξ y u y y + ξ z u z y + ξ t u y t , η u z = D z ( η u ξ x u x ξ y u y ξ z u z ξ t u t ) + ξ x u x z + ξ y u y z + ξ z u z z + ξ t u z t , η u t = D t ( η u ξ x u x ξ y u y ξ z u z ξ t u t ) + ξ x u x t + ξ y u y t + ξ z u z t + ξ t u t t ,
where D x , D y , D z , and D t are total derivatives [25,26].
At the same time, it is supposed that Equation (1) has the following general infinitesimal generator:
X = ξ x x + ξ y y + ξ z z + ξ t t + η u u .
Utilizing the infinitesimal invariance criterion [11,12,13,14,15,16,17,18,19,20], one has:
P r ( α , t ) X ( Ξ ) | Ξ = 0 = 0 ,
where P r ( α , t ) X is the prolongation of X and:
Ξ = D t α u + u x + u 2 u y u x z t .
For Equation (1), it can be written as:
p r ( α , t ) X = X + D t α u + η u x u x + η u y u y + η u x z t u x z t .
Applying the invariance condition (4) to Equation (1), we can obtain the symmetry determining equation as follows:
η u ( α , t ) + η u x + 2 η u u y + u 2 η u y η u x z t = 0 ,
where the α -extended infinitesimal [17,18,19,20,21,22,23,24] is:
η u ( α , t ) = α η u t α + ( η u , u α D t ( ξ t ) ) α u t α u α η u , u t α + m = 1 [ α m α η u , u t α α m + 1 D t m + 1 ( ξ t ) ] D t α m ( u ) m = 1 α m D t m ( ξ x ) D t α m ( u x ) m = 1 α m D t m ( ξ y ) D t α m ( u y ) m = 1 α m D t m ( ξ z ) D t α m ( u z ) + Δ ,
with:
Δ = m = 2 n = 2 m j = 2 n r = 0 j 1 α m m n j r 1 j ! t m α Γ ( m + 1 α ) ( u ) r n t n ( u j r ) m n + j η u t m n u j .
Inserting the Equations (7) and (2) into Equation (6) and collecting the coefficients of different powers of u to zero, we have:
ξ x = 2 c 2 α x + c 3 , ξ y = 2 c 2 y + c 5 , ξ z = 2 c 2 z + c 4 , ξ t = 2 c 2 t + c 1 , η u = c 2 α u + 2 c 2 u ,
where c i ( i = 1 , 2 , 3 , 4 , 5 ) are arbitrary parameters.
In addition, the constraint condition ξ t ( x , y , z , t , u ) | t = 0 = 0 to hold the Riemann–Liouville fractional derivative structure. We get c 1 = 0 .
Hence, the infinitesimal generator X of Equation (1) has:
X = ( 2 c 2 α x + c 3 ) x + ( 2 c 2 y + c 5 ) y + ( 2 c 2 z + c 4 ) z + 2 c 2 t t + ( c 2 α u + 2 c 2 u ) u .
It can be spanned by the following four vector fields:
X 1 = 2 x x + 2 α y y 2 α z z + 2 α α u u + 2 α t t , X 2 = x , X 3 = y , X 4 = z .
This concludes the proof. □
The above four vector fields (10) can construct a closed Lie algebra. Their commutation relationships can be found in Table 1.

3. Optimal System of Equation (1)

In order to reduce the original equation as much as possible and construct the exact solutions, we need to find the one-dimensional optimal system of Equation (1). With the help of the Lie series [25,26], the action of the adjoint operator can given by:
A d ( exp ( ε X i ) ) · X j = X j ε [ X i , X j ] + ε 2 2 [ X i , [ X i , X j ] ] ,
where [ X i , X j ] is Lie bracket and ε is a constant.
Specifically:
A d ( exp ( ε · X 1 ) ) · X 2 = X 2 ε · [ X 1 , X 2 ] + ε 2 2 [ X 1 , [ X 1 , X 2 ] ] = e 2 α ε X 2 ,
A d ( exp ( ε · X 1 ) ) · X 3 = X 3 ε · [ X 1 , X 3 ] + ε 2 2 [ X 1 , [ X 1 , X 3 ] ] = 2 X 3 ,
A d ( exp ( ε · X 2 ) ) · X 1 = X 1 ε · [ X 2 , X 1 ] + ε 2 2 [ X 2 , [ X 2 , X 1 ] ] = X 1 e 2 α ε X 2 ,
A d ( exp ( ε · X 2 ) ) · X 1 = X 1 ε · [ X 2 , X 1 ] + ε 2 2 [ X 2 , [ X 2 , X 1 ] ] = X 1 2 ε X 3 .
The adjoint actions of the infinitesimal generators (10) are seen as Table 2.
Hence, we have the following result of the one-dimensional optimal system of Equation (1).
Theorem 2. 
An optimal system of one-dimensional Lie subalgebras operators (10) to Equation (1) can be spanned by:
V 1 = X 1 , V 2 = k 2 X 2 + k 3 X 3 + k 4 X 4 ,
where k i ( i = 2 , 3 , 4 ) are not all zero.
The simple proof process of this Theorem 2 can be found in references [25,26]. Here, we omit it.
Remark 1. 
If making linear operator L = a 1 X 1 + a 2 X 2 + a 3 X 3 + a 4 X 4 , a j 0 ( j = 1 , 2 , 3 , 4 ) , then the Killing form is 4 α 2 a 1 2 + 8 a 1 2 . It can be seen that the value α = 1 is conform to Equation (1) of integer order for the Killing form 12 a 1 2 . That is to say, fractional differential equations have a wider research domain.

4. One-Parameter Lie Transformation Group of Equation (1)

In this subsection, we derived the one-parameter Lie transformation group [25,26] of Equation (1). It keep the solution set unchanged or transform the solution into a solution to differential equations. Hence, we may obtain new solutions from the known ones of this considered model.
In addition to constructing the one-parameter Lie transformation group, we solved the following initial problems [25,26]:
d ( x * ( ε ) ) d ε = ξ x ( x * ( ε ) , y * ( ε ) , z * ( ε ) , t * ( ε ) , u * ( ε ) ) , x * ( 0 ) = x , d ( y * ( ε ) ) d ε = ξ y ( x * ( ε ) , y * ( ε ) , z * ( ε ) , t * ( ε ) , u * ( ε ) ) , y * ( 0 ) = y , d ( z * ( ε ) ) d ε = ξ z ( x * ( ε ) , y * ( ε ) , z * ( ε ) , t * ( ε ) , u * ( ε ) ) , z * ( 0 ) = z , d ( t * ( ε ) ) d ε = ξ t ( x * ( ε ) , y * ( ε ) , z * ( ε ) , t * ( ε ) , u * ( ε ) ) , t * ( 0 ) = t , d ( u * ( ε ) ) d ε = η u ( x * ( ε ) , y * ( ε ) , z * ( ε ) , t * ( ε ) , u * ( ε ) ) , u * ( 0 ) = u ,
where ε is a small parameter.
Solving the above system (13) with vector fields (10), we can obtain the one-parameter Lie transformation groups g i ( ε ) , ( i = 1 , 2 , 3 , 4 ) of the forms:
g 1 ( ε ) : ( x , y , z , t , u ) ( e 2 ε x , e 2 α ε y , e 2 α ε z , e 2 α ε t , e 2 α α ε u ) , g 2 ( ε ) : ( x , y , z , t , u ) ( x + ε , y , z , t , u ) , g 3 ( ε ) : ( x , y , z , t , u ) ( x , y + ε , z , t , u ) , g 4 ( ε ) : ( x , y , z , t , u ) ( x , y , z + ε , t , u ) .
As a result, we obtained Theorem 3.
Theorem 3. 
If u = f ( x , y , z , t ) is an exact solution of the time fractional higher-dimensional nonlinear modified equation of wave propagation (1), then:
u 1 = e α 2 α ε f 1 ( e 2 ε x , e 2 α ε y , e 2 α ε z , e 2 α ε t ) , u 2 = f 2 ( x ε , y , z , t ) , u 3 = f 3 ( x , y ε , z , t ) , u 4 = f 4 ( x , y , z ε , t ) ,
are also solutions to Equation (1).
Remark 2. 
It is can be seen that the Lie transformation group g 2 , g 3 , and g 4 are invariant on spaces x, y, and z. The well known scaling symmetry action is g 1 .

5. Similarity Reduction and Similarity Solution of Equation (1)

In this section, we reduce Equation (1) into the lower dimensional fractional differential equation with the Erdélyi–Kober fractional operators [27]. In addition, we obtained the similarity solution of the reduced equation.

5.1. Similarity Reduction of Equation (1)

Case 1.  V 1 = X 1 = 2 x x + 2 α y y 2 α z z + 2 α α u u + 2 α t t
From V 1 , Equation (1) was reduced into the following form.
Theorem 4. 
The similarity variables ξ 1 = x t α , ξ 2 = y t 1 , ξ 3 = z t and group invariant solution u = t 1 1 2 α f ( ξ 1 , ξ 2 , ξ 3 ) were used to reduce Equation (1) into the following result:
( P 1 α , 1 , 1 2 3 2 α , α f ) ( ξ 1 , ξ 2 , ξ 3 ) + f ξ 1 + f f ξ 2 2 ( 1 3 2 α ) t f ξ 1 ξ 2 = 0 ,
where the ( P 1 α , 1 , 1 2 3 2 α , α f ) ( ξ 1 , ξ 2 , ξ 3 ) represents the three parameters Erdélyi–Kober fractional differential operator [27]
( P χ 1 , χ 2 , χ 3 τ , α F ) ( ξ 1 , ξ 2 , ξ 3 ) = j = 0 n 1 ( τ + j 1 χ 1 ξ 1 ξ 1 1 χ 2 ξ 2 ξ 2 1 χ 3 ξ 3 ξ 3 ) ( K χ 1 , χ 2 , χ 3 τ + α , n α F ) ( ξ 1 , ξ 2 , ξ 3 ) , χ i R ,
where:
( K χ 1 , χ 2 , χ 3 τ , α F ) ( ξ 1 , ξ 2 , ξ 3 ) = 1 Γ ( α ) 1 ( ϑ 1 ) α 1 ϑ ( τ + α ) F ( ξ 1 ϑ 1 χ 1 , ξ 2 ϑ 1 χ 2 , ξ 3 ϑ 1 χ 3 ) d ϑ , α > 0 , F ( ξ 1 , ξ 2 , ξ 3 ) , α = 0 .
is the three parameters Erdélyi–Kober fractional integral operator [27].
Proof. 
For V 1 , it has the characteristic equation of the form:
d x 2 x = d z 2 α z = d y 2 α y = d t 2 α t = d u 2 α α u .
Solving characteristic Equation (16), we obtained, respectively, the similarity variables and group invariant solution:
ξ 1 = x t α , ξ 2 = y t 1 , ξ 3 = z t ,
and
u = t 1 1 2 α f ( ξ 1 , ξ 2 , ξ 3 ) .
Next, for the Riemann–Liouville fractional derivative definition, we have:
D t α u = n t n [ 1 Γ ( n α ) 0 t ( t s ) n α 1 s 1 1 2 α f ( x s α , y s 1 , z s ) d s ] , n 1 < α < n , ( n = 1 , 2 , 3 , ) .
If considering v = t s , then expression (19) becomes:
D t α u = n t n [ t n + 1 3 2 α Γ ( n α ) 1 ( v 1 ) n α 1 v ( n + 2 3 2 α ) f ( ξ 1 v α , ξ 2 v , ξ 3 v 1 ) d v ] = n t n [ t n + 1 3 2 α ( K 1 α , 1 , 1 2 1 2 α , n α f ) ( ξ 1 , ξ 2 , ξ 3 ) ] .
Meanwhile, we note the following fact:
t t ϕ ( ξ 1 , ξ 2 , ξ 3 ) = α ξ 1 ϕ ξ 1 ξ 2 ϕ ξ 2 + ξ 3 ϕ ξ 3 .
Hence, Equation (20) can be further reduced into the following form:
n t n [ t n + 1 3 2 α ( K 1 α , 1 , 1 2 1 2 α , n α f ) ( ξ 1 , ξ 2 , ξ 3 ) ] = n 1 t n 1 [ t ( t n 3 2 α ( K 1 α , 1 , 1 2 1 2 α , n α f ) ( ξ 1 , ξ 2 , ξ 3 ) ) ] = n 1 t n 1 [ t n 3 2 α ( n + 1 3 2 α α ξ 1 ξ 1 ξ 2 ξ 2 + ξ 3 ξ 3 ) · ( K 1 α , 1 , 1 2 1 2 α , n α f ) ( ξ 1 , ξ 2 , ξ 3 ) ] = ( n 1 ) t i m e s = t 1 3 2 α j = 0 n 1 ( 2 3 2 α + j α ξ 1 ξ 1 ξ 2 ξ 2 + ξ 3 ξ 3 ) · ( K 1 α , 1 , 1 2 1 2 α , n α f ) ( ξ 1 , ξ 2 , ξ 3 ) = t 1 3 2 α ( P 1 α , 1 , 1 2 3 2 α , α f ) ( ξ 1 , ξ 2 , ξ 3 ) .
In summary, Equation (1) was reduced into the following form:
( P 1 α , 1 , 1 2 3 2 α , α f ) ( ξ 1 , ξ 2 , ξ 3 ) + f ξ 1 + f f ξ 2 2 ( 1 3 2 α ) t f ξ 1 ξ 2 = 0 .
Proving the theorem. □
Case 2.  V 2 = X 2 = x
Similar to Case 1, Equation (1) easily was reduced into the following expression:
D t α u + u 2 u y = 0 .
For Equation (24), the infinitesimal generators can be obtained through the fractional Lie group scheme:
X 1 ^ = 2 t t + u ( α 1 ) u + 2 y y , X 2 ^ = y .
Similarity, the infinitesimal generators (25) were used to reduce Equation (24) as Theorem 4 processes. Here, we ignore it.
Case 3.  V 3 = X 3 = y
Similar to Case 1, Equation (1) can be reduced into the following form:
D t α u + u x u x z t = 0 .
If Theorem 1 processes are performed, then the infinitesimal generators are as follows:
X 1 ¯ = z , X 2 ¯ = u u , X 3 ¯ = x , X 4 ¯ = α x x z z + t t .
Similarly, the infinitesimal generators (27) were used to reduce Equation (26). Here, we ignore it.
Case 4.  V 4 = X 4 = z
Similar to Case 1, Equation (1) can be reduced into the following form
D t α u + u x + u 2 u y = 0 .
If perform Theorem 1 processes, then we have
X 1 ¯ = 2 t t + u ( 1 α ) u + 2 α x x + 2 y y , X 2 ¯ = x , X 3 ¯ = y .
Similarity, the infinitesimal generators (29) were used to reduce Equation (27) as Theorem 4 processes. Here we ignore it.
Case 5.  V 5 = x + y
In a similar way, Equation (1) was reduced into the following form:
D t α u = 0 .
It is easy to find an exact solution as follows:
u = c Γ ( α ) t α 1 ,
where c is a free constant.
Taking the values c = ± 1 , α = 0.3 , α = 0.6 , and α = 0.9 into rational solution (31), Figure 1 was plotted.
Case 6.  V 6 = x + z
From V 6 , Equation (1) was reduced into Equation (24).
Case 7.  V 7 = y + z
In a similar way, Equation (1) was reduced into the following form:
D t α u + u x = 0 .
This is known as the one-dimensional linear Burger equation.
Case 8.  V 7 = x + y + z
From V 8 , Equation (1) was reduced into Equation (30).

5.2. Similarity Solution of Equation (24)

In this subsection, we proved firstly the reduced Equation (24) existence the similarity solution. Then, we presented specific expressions of similar solution.
It is supposed that Equation (24) is invariant under the following scaling transformation:
t = λ t ˜ , y = λ p y ˜ , u = λ q u ˜ ( y ˜ , t ˜ ) ,
where p and q are arbitrary parameters to be determined later.
Inserting transformations (33) into Equation (24), we have:
λ q α α u ˜ t α ˜ + λ 3 q p u 2 ˜ u ˜ y ˜ = 0 ,
which implies that:
p α = 2 q
has similarity solution.
We can use the following to find traveling wave similarity solution for Equation (24).
( t t ˜ ) = ( y y ˜ ) 1 p , ( u u ˜ ) = ( y y ˜ ) q p
Let the similarity transformation be as follows:
u = y q p U ( ξ ) , ξ = c t y 1 ,
where p, q, and c are arbitrary parameters to be determined later.
According to the definition of the Riemann–Liouville, one has:
α u t α = 1 Γ ( 1 α ) t 0 t y q p ( c t y 1 ) ( t x ) α d x .
Let x = c t y 1 , then Equation (37) becomes:
α u t α = 1 Γ ( 1 α ) c y d d ξ 1 ξ y q p ( x ) ( ξ x ) α ( y c ) α y c d x = c α y q α q p d α U ( ξ ) d ξ .
Furthermore:
u y = y q p p [ q p U ( ξ ) ( ξ + 1 ) U ( ξ ) ] .
Substituting (39) and (38) into Equation (24), we have:
c α y q α q p d α U ( ξ ) d ξ + y 3 q p p U 2 ( ξ ) T ( ξ ) = 0 ,
which yields:
q p = 1 α 2 ,
where:
T ( ξ ) = q p U ( ξ ) ( ξ + 1 ) U ( ξ ) .
Thereby, Equation (40) becomes the following ordinary differential equation of fractional order:
c α d α U ( ξ ) d ξ + U 2 ( ξ ) T ( ξ ) = 0 .
Now, we consider that Equation (42) has the special solutions of the following form:
U ( ξ ) = k ξ ρ , ξ 0 , 0 , ξ < 0 .
Plugging (43) into Equation (42), one has:
c α k Γ ( 1 + ρ ) Γ ( 1 + ρ α ) ξ ρ α + k 3 ξ 3 ρ 1 [ q p ξ ( ξ + 1 ) ρ ] = 0 ,
which derives:
ρ = 1 α 2 .
Inserting the comparison expressions (45) and (41) into Equation (44), we have:
k = [ c α Γ ( 1 2 1 2 α ) Γ ( 3 2 3 2 α ) ] 1 2 .
Conversely:
c = [ k 2 Γ ( 3 2 3 2 α ) Γ ( 1 2 1 2 α ) ] 1 α .
Therefore, we obtained the traveling wave similarity solution of Equation (24) of the form:
u ( y , t ) = y q p k ξ ρ = k { [ k 2 Γ ( 3 2 3 2 α ) Γ ( 1 2 1 2 α ) ] 1 α t y } 1 α 2 .

6. Conclusions and Discussion

In this paper, we studied the time fractional higher-dimensional nonlinear modified equation of wave propagation. As a result of obtaining novelty, the symmetries, one-dimensional optimal system, and one-parameter Lie transformation group of this considered model were derived. Next, the original Equation (1) carried out the similarity reduction of the lower differential equation. On the basis of the reduced equation, we constructed a similarity solution. Although the similarity solutions have been obtained for the reduced model, it is still difficult to seek other types of exact solutions, such as the soliton solutions, rational solutions, and Kink solutions. Therefore, how to obtain the exact solutions of the reduced equation is one of our future key tasks. Finally, these beautiful results can help us discover more evolutionary mechanisms of this considered model.

Author Contributions

Methodology, J.-G.L. formal analysis, J.-G.L. and Y.-Y.F.; writing original draft preparation, J.-G.L.; writing review and editing, J.-G.L. and Y.-Y.F.; funding acquisition, J.-G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (No. 11971475), the Natural Science Foundation of Jiangsu Province (No. BK20230708), and the Natural Science Foundation for the Universities in Jiangsu Province (No. 23KJB110003).

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors are grateful to the anonymous referees for their useful comments and suggestions.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Atangana, A.; Nieto, J.J. Numerical solution for the model of RLC circuit via the fractional derivative without singular kernel. Adv. Mech. Eng. 2015, 7, 1687814015613758. [Google Scholar] [CrossRef]
  2. Schaub, M.; Abadi, F. Integrated population models: A novel analysis framework for deeper insights into population dynamics. J. Ornithol. 2011, 152, 227–237. [Google Scholar] [CrossRef]
  3. Letellier, C.; Aziz-Alaoui, M.A. Analysis of the dynamics of a realistic ecological model. Chaos Solitons Fractals 2002, 13, 95–107. [Google Scholar] [CrossRef]
  4. Lax, P.D. Periodic solutions of the KdV equation. Commun. Pure Appl. Math. 1975, 28, 141–188. [Google Scholar] [CrossRef]
  5. Caraballo Garrido, T.; Ngoc, T.B.; Thach, T.N.; Tuan, N.H. On initial value and terminal value problems for subdiffusive stochastic Rayleigh-Stokes equation. Discret. Contin. Dyn. Syst. Ser. B 2020, 26, 4299–4323. [Google Scholar] [CrossRef]
  6. Wazwaz, A.-M. Exact soliton and kink solutions for new (3 + 1)-dimensional nonlinear modified equations of wave propagation. Open Eng. 2017, 7, 169–174. [Google Scholar] [CrossRef]
  7. Samko, S.; Kilbas, A.A.; Marichev, O. Fractional Integral and Derivatives: Theory and Applications; Gordon and Breach Science: Yverdon, Switzerland, 1993. [Google Scholar]
  8. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  9. Yang, X.J. General Fractional Derivatives: Theory, Methods and Applications; CRC Press: New York, NY, USA, 2019. [Google Scholar]
  10. Wang, K.-J. A fractal modification of the unsteady Korteweg–de Vries model and its generalized fractal variational principle and diverse exact solutions. Fractals 2022, 30, 2250192. [Google Scholar] [CrossRef]
  11. Henry, B.I.; Langlands, T.A.; Straka, P. An introduction to fractional diffusion. Complex Phys. Biophys. Econophys. Syst. 2010, 9, 37–89. [Google Scholar]
  12. Seadawy, A.R.; Ali, K.K.; Nuruddeen, R.I. Avariety of soliton solutions for the fractional Wazwaz-Benjamin-Bona-Mahony equations. Results Phys. 2019, 12, 2234–2241. [Google Scholar] [CrossRef]
  13. Liu, Q.; Zhou, Y.; Li, K.; Zhang, S. Application of the dynamical system method and the deep learning method to solve the new (3 + 1)-dimensional fractional modified Benjamin-Bona-Mahony equation. Nonlinear Dyn. 2022, 110, 3737–3750. [Google Scholar] [CrossRef]
  14. Shakeel, M.; Turki, N.B.; Shah, N.A.; Tag, S.M. Diversity of Soliton Solutions to the (3 + 1)-Dimensional Wazwaz-Benjamin-Bona-Mahony Equations Arising in Mathematical Physics. Results Phys. 2023, 51, 106624. [Google Scholar] [CrossRef]
  15. Bilal, M.; Younas, U.; Baskonus, H.M.; Younis, M. Investigation of shallow water waves and solitary waves to the conformable 3D-WBBM model by an analytical method. Phys. Lett. A 2021, 403, 127388. [Google Scholar] [CrossRef]
  16. Rezazadeh, H.; Inc, M.; Baleanu, D. New solitary wave solutions for variants of (3 + 1)-dimensional Wazwaz-Benjamin-Bona-Mahony equations. Front. Phys. 2020, 8, 332. [Google Scholar] [CrossRef]
  17. Adeyemo, O.D.; Zhang, L.; Khalique, C.M. Optimal solutions of Lie subalgebra, dynamical system, travelling wave solutions and conserved currents of (3 + 1)-dimensional generalized Zakharov-Kuznetsov equation type I. Eur. Phys. J. Plus 2022, 137, 954. [Google Scholar] [CrossRef]
  18. Liu, J.-G.; Zhang, Y.-F.; Wang, J.-J. Investigation of the time fractional generalized (2 + 1)-dimensional Zakharov-Kuznetsov equation with single-power law nonlinearity. Fractals 2023, 31, 2350033. [Google Scholar] [CrossRef]
  19. Liu, J.-G.; Yang, X.-J.; Feng, Y.-Y.; Cui, P.; Geng, L.-L. On integrability of the higher-dimensional time fractional KdV-type equation. J. Geom. Phys. 2021, 160, 104000. [Google Scholar] [CrossRef]
  20. Liu, J.-G.; Yang, X.-J.; Geng, L.-L.; Yu, X.-J. On fractional symmetry group scheme to the higher dimensional space and time fractional dissipative Burgers equation. Int. J. Geom. Methods Mod. Phys. 2022, 19, 2250173. [Google Scholar] [CrossRef]
  21. Liu, J.-G.; Geng, F.-Z.; Li, X. Supplement a high-dimensional time fractional diffusion equation. Alex. Eng. J. 2023, 15, 459–464. [Google Scholar] [CrossRef]
  22. Sahoo, S.; Ray, S.S. Analysis of Lie symmetries with conservation laws for the (3+1) dimensional time-fractional mKdV-ZK equation in ion-acoustic waves. Nonlinear Dyn. 2017, 90, 1105–1113. [Google Scholar] [CrossRef]
  23. Jiang, Z.; Zhang, Z.G.; Li, J.J.; Yang, H.W. Analysis of Lie symmetries with conservation laws and solutions of generalized (4 + 1)-dimensional time-fractional Fokas equation. Fractal Fract. 2022, 6, 108. [Google Scholar] [CrossRef]
  24. Zhu, H.-M.; Zhang, Z.-Y.; Zheng, J. The time-fractional (2 + 1)-dimensional Hirota-Satsuma-Ito equations: Lie symmetries, power series solutions and conservation laws. Commun. Nonlinear Sci. Numer. Simul. 2022, 115, 106724. [Google Scholar] [CrossRef]
  25. Bluman, G.W.; Anco, S. Symmetry and Integration Methods for Differential Equations; Springer: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  26. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer: Berlin/Heidelberg, Germany, 1986. [Google Scholar]
  27. Kiryakova, V.S. Generalized Fractional Calculus and Applications; Longman Scientific & Technical, Longman Group: Harlow, UK, 1994. [Google Scholar]
Figure 1. Plots of evolution of rational solution (31) with the above parameter values.
Figure 1. Plots of evolution of rational solution (31) with the above parameter values.
Fractalfract 08 00124 g001
Table 1. The commutation table of Lie algebra.
Table 1. The commutation table of Lie algebra.
[ X i , X j ] X 1 X 2 X 3 X 4
X 1 0 2 α X 2 2 X 3 2 X 4
X 2 2 α X 2 000
X 3 2 X 3 000
X 4 2 X 4 000
Table 2. The adjoint representation of infinitesimal generators (10).
Table 2. The adjoint representation of infinitesimal generators (10).
Ad ( exp ( ε X i ) ) · X j X 1 X 2 X 3 X 4
X 1 X 1 e 2 α ε X 2 e 2 ε X 3 e 2 ε X 4
X 2 X 1 e 2 α ε X 2 X 2 X 3 X 4
X 3 X 1 2 ε X 3 X 2 X 3 X 4
X 4 X 1 + 2 ε X 4 X 2 X 3 X 4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, J.-G.; Feng, Y.-Y. Investigation of the Time Fractional Higher-Dimensional Nonlinear Modified Equation of Wave Propagation. Fractal Fract. 2024, 8, 124. https://doi.org/10.3390/fractalfract8030124

AMA Style

Liu J-G, Feng Y-Y. Investigation of the Time Fractional Higher-Dimensional Nonlinear Modified Equation of Wave Propagation. Fractal and Fractional. 2024; 8(3):124. https://doi.org/10.3390/fractalfract8030124

Chicago/Turabian Style

Liu, Jian-Gen, and Yi-Ying Feng. 2024. "Investigation of the Time Fractional Higher-Dimensional Nonlinear Modified Equation of Wave Propagation" Fractal and Fractional 8, no. 3: 124. https://doi.org/10.3390/fractalfract8030124

APA Style

Liu, J. -G., & Feng, Y. -Y. (2024). Investigation of the Time Fractional Higher-Dimensional Nonlinear Modified Equation of Wave Propagation. Fractal and Fractional, 8(3), 124. https://doi.org/10.3390/fractalfract8030124

Article Metrics

Back to TopTop