Next Article in Journal
Realization of Intelligent Observer for Sensorless PMSM Drive Control
Next Article in Special Issue
A Novel Approach to Solving Fractional-Order Kolmogorov and Rosenau–Hyman Models through the q-Homotopy Analysis Transform Method
Previous Article in Journal
An Analysis of Actors in Malay Films: Small Worlds, Centralities and Genre Diversity
Previous Article in Special Issue
Non-Overlapping Domain Decomposition via BURA Preconditioning of the Schur Complement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Innovative Approach to Nonlinear Fractional Shock Wave Equations Using Two Numerical Methods

Department of Mathematics, College of Science, University of Bisha, Bisha 61922, Saudi Arabia
Mathematics 2023, 11(5), 1253; https://doi.org/10.3390/math11051253
Submission received: 14 February 2023 / Revised: 25 February 2023 / Accepted: 2 March 2023 / Published: 5 March 2023

Abstract

:
In this research, we propose a combined approach to solving nonlinear fractional shock wave equations using an Elzaki transform, the homotopy perturbation method, and the Adomian decomposition method. The nonlinear fractional shock wave equation is first transformed into an equivalent integral equation using the Elzaki transform. The homotopy perturbation method and Adomian decomposition method are then utilized to approximate the solution of the integral equation. To evaluate the effectiveness of the proposed method, we conduct several numerical experiments and compare the results with existing methods. The numerical results show that the combined method provides accurate and efficient solutions for nonlinear fractional shock wave equations. Overall, this research contributes to the development of a powerful tool for solving nonlinear fractional shock wave equations, which has potential applications in many fields of science and engineering. This study presents a solution approach for nonlinear fractional shock wave equations using a combination of an Elzaki transform, the homotopy perturbation method, and the Adomian decomposition method. The Elzaki transform is utilized to transform the nonlinear fractional shock wave equation into an equivalent integral equation. The homotopy perturbation method and Adomian decomposition method are then employed to approximate the solution of the integral equation. The effectiveness of the combined method is demonstrated through several numerical examples and compared with other existing methods. The results show that the proposed method provides accurate and efficient solutions for nonlinear fractional shock wave equations.

1. Introduction

Fractional differential equations (FDEs) provide a mathematical framework for modeling systems that exhibit nonlinear or complex dynamics. Unlike traditional differential equations, FDEs utilize fractional derivatives, which capture the rate of change in these systems. FDEs have found widespread applications in fields such as physics, engineering, finance, and biology, where they have been used to model phenomena such as anomalous diffusion, viscoelastic materials, and fractional control systems [1,2,3,4,5].
Fractional derivatives are a generalization of classical integer derivatives and can be defined in various ways, including the Riemann–Liouville, Caputo, and Grunwald–Letnikov derivatives. Various numerical methods, such as the finite difference, boundary element, and variational methods, can be used to solve FDEs. The study of FDEs has a rich history dating back to the 17th century, when mathematician John Wallis first introduced the concept of fractional calculus [6,7,8]. In recent years, FDEs have gained significant attention and become increasingly popular due to their ability to model complex systems and their broad range of applications [9,10,11,12,13].
The nonlinear shock wave equation, also known as Burgers’ equation, is a fundamental equation in fluid dynamics and nonlinear physics. Its history can be traced back to the early 20th century, with important contributions by several mathematicians and physicists. In the following decades, mathematicians and physicists made significant contributions to the understanding of Burgers’ equation and its properties. In the 1950s, the Russian mathematician Andrei Kolmogorov studied the behavior of solutions to the equation in the limit of high Reynolds numbers. He showed that the solutions exhibit self-similarity and developed a theory of turbulence based on this observation [14,15]. In the 1960s, the American mathematician Richard Courant and his collaborators developed numerical methods for solving Burgers’ equation, which played a key role in the development of computational fluid dynamics. In the 1970s and 1980s, the equation was studied extensively in the context of nonlinear wave propagation and shock waves [16,17,18,19]. Today, Burgers’ equation continues to be an active area of research in mathematics, physics, and engineering. It has applications in fields such as fluid mechanics, nonlinear optics, and plasma physics [20,21,22,23,24].
The nonlinear shock wave equation is a mathematical model that describes the behavior of shock waves, which are sudden changes in the pressure and density of a fluid that propagates through it. Unlike linear shock wave equations, the nonlinear shock wave equation takes into account the nonlinear effects that occur when the fluid pressure exceeds a critical value [25,26]. These nonlinear effects cause the shock wave to become distorted and the velocity of the fluid to change rapidly, leading to complex and diverse behaviors. The nonlinear shock wave equation plays an important role in many areas of science and engineering, including fluid dynamics, combustion, and blast waves [27,28,29].
The fractional nonlinear shock wave equation is a mathematical tool extensively employed to describe the intricate behavior of shock waves in different engineering and physical applications. Due to the nonlinear nature of shock waves, their behavior can be complicated and unpredictable. This equation incorporates the nonlinear effects, allowing for an accurate representation of the shock wave’s dynamics over time. Moreover, the use of fractional derivatives in this equation provides a more precise description of shock wave behavior. The application of this equation is particularly valuable in fluid dynamics research, where it is used to model shock waves in fluids with complex properties such as viscoelastic fluids. Additionally, the equation finds its applications in material science, where it can be employed to simulate the behavior of shock waves in solids, including their responses to stress and deformation [30,31,32].
The highly nonlinear nature of the fractional nonlinear shock wave equation implies that the solution is dependent on the initial conditions and underlying physical laws. Due to this nonlinearity, there are diverse shock wave behaviors, such as propagating waves, standing waves, and bifurcation points. These behaviors can exhibit multiple distinct characteristics. Despite its complexity, the fractional nonlinear shock wave equation is a crucial tool in comprehending shock wave behavior, leading to a more accurate and thorough understanding of these complex physical phenomena. Continued research and development of this equation is promising, and it is expected to offer new insights into a wide range of scientific and engineering disciplines. Its implications will be essential in several areas of research and development [33,34].
The homotopy perturbation method (HPM) and Elzaki transform are two mathematical techniques that are widely used to solve nonlinear and linear partial differential equations (PDEs) [35,36]. The HPM is a numerical method that is based on the concept of homotopy and is used to approximate solutions to nonlinear problems. On the other hand, the Elzaki transform is a mathematical tool that is used to transform linear problems into nonlinear ones, making them easier to solve. Both the HPM and Elzaki transform have been extensively researched and applied in various fields, including engineering, physics, and mathematics [37,38].
The Adomian decomposition method (ADM) has been widely used for solving various types of nonlinear differential equations. This method is based on the decomposition of the nonlinear term into a series of polynomials and then using this series to find the solution. However, the convergence of the ADM series is slow, which makes it unsuitable for some applications. To overcome this limitation, various researchers have proposed combining the ADM with other mathematical techniques to improve its efficiency and accuracy. One of these techniques is the Elzaki transform (ET), which has been used to solve linear and nonlinear partial differential equations [39,40]. The combination of the ADM and ET has been applied to solve various types of nonlinear differential equations, including nonlinear partial differential equations, nonlinear integro-differential equations, and nonlinear differential-difference equations (Elzaki and Mustafa, 2011; Khodaei et al., 2015; Mahmoud, 2018). The results obtained from these studies demonstrate the efficiency and accuracy of the combined method in solving nonlinear differential equations [39,40].
The following is the rest of this study. Section 2 compiles some basic definitions. Section 3 introduces the concept of the HPTM, whereas Section 4 introduces the concept of the ETDM. The proposed methods’ convergence analysis is provided in Section 5. In Section 6, various approximative solutions to the fractional nonlinear wave equation are derived by using the form of the initial value, and the structure of the solutions is displayed using graphs and tables. In Section 6, the work’s conclusion is examined.

2. Basic Definitions

 Definition 1.
The fractional derivative D ß in the Abel–Riemann sense, having an order ß, is given as follows [41,42]:
D ß θ ( μ ) = d κ d μ κ θ ( μ ) , ß = κ , 1 Γ ( κ ß ) d d μ κ 0 μ θ ( μ ) ( μ ϕ ) ß κ + 1 d ϕ , κ 1 < ß < κ ,
where κ Z + , ß R + and
D ß θ ( μ ) = 1 Γ ( ß ) 0 μ ( μ ϕ ) ß 1 θ ( ϕ ) d ϕ , 0 < ß 1 .
 Definition 2.
The fractional integration operator κ ϕ in the Abel–Riemann sense is defined as [41,42]
κ ß θ ( μ ) = 1 Γ ( ß ) 0 μ ( μ ϕ ) ß 1 θ ( μ ) d μ , μ > 0 , ß > 0 .
having the properties
κ ß μ κ = Γ ( κ + 1 ) Γ ( κ + ß + 1 ) μ κ + ϕ , D ß μ κ = Γ ( κ + 1 ) Γ ( κ ß + 1 ) μ κ ϕ .
 Definition 3.
The Caputo derivative D ß of the fractional order ß is given by [41,42]
C D ß θ ( μ ) = 1 Γ ( κ ß ) 0 μ θ κ ( ϕ ) ( μ ϕ ) ß κ + 1 d ϕ , κ 1 < ß < κ , d κ d μ κ θ ( μ ) , κ = ß .
 Definition 4.
κ μ ß D μ ß θ ( μ ) = θ ( μ ) k = 0 m θ k ( 0 + ) μ k k ! , f o r μ > 0 , a n d κ 1 < ß κ , κ N . D μ ß κ μ ß θ ( μ ) = θ ( μ ) .
 Definition 5.
The Caputo operator in terms of the Elzaki transform is given by [41,42]
E [ D μ ß θ ( μ ) ] = s ß E [ θ ( μ ) ] k = 0 κ 1 s 2 ß + k θ ( k ) ( 0 ) , w h e r e κ 1 < ß < κ .

3. General Methodology of the HPTM

This section focuses on the core method used by the HPTM to find solutions for the general partial differential equation
D ϑ ß V ( , ϑ ) = P 1 [ ] V ( , ϑ ) + Q 1 [ ] V ( , ϑ ) , 1 < ß 2 ,
subjected to the initial conditions
V ( , 0 ) = ξ ( ) , ϑ V ( , 0 ) = ζ ( ) .
Here, D ϑ ß = ß ϑ ß is the fractional Caputo derivative and P 1 [ ] and Q 1 [ ] are the linear and nonlinear operators, respectively.
By utilizing the Elzaki transformation in Equation (3), we achieve
E [ D ϑ ß V ( , ϑ ) ] = E [ P 1 [ ] V ( , ϑ ) + Q 1 [ ] V ( , ϑ ) ] ,
1 u ß { V ( , ϑ ) u V ( 0 ) u 2 V ( 0 ) } = E [ P 1 [ ] V ( , ϑ ) + Q 1 [ ] V ( , ϑ ) ] .
After simplifying, the result is obtained:
V ( , ϑ ) = u V ( 0 ) + u 2 V ( 0 ) + u ß E [ P 1 [ ] V ( , ϑ ) + Q 1 [ ] V ( , ϑ ) ] .
By applying the inverse Elzaki transform to both sides, we obtain
V ( , ϑ ) = V ( 0 ) + V ( 0 ) + E 1 [ u ß E [ P 1 [ ] V ( , ϑ ) + Q 1 [ ] V ( , ϑ ) ] ] .
Now, by using the HPM, we have
V ( , ϑ ) = k = 0 ϵ k V k ( , ϑ ) .
The decomposition of the nonlinear terms can be expressed as follows. With the perturbation parameter ϵ [ 0 , 1 ] , we have
Q 1 [ ] V ( , ϑ ) = k = 0 ϵ k H n ( V ) .
He’s polynomials, designated as H k ( V ) , are expressed by
H n ( V 0 , V 1 , . . . , V n ) = 1 Γ ( n + 1 ) D ϵ k Q 1 k = 0 ϵ i V i ϵ = 0 .
where D ϵ k = k ϵ k .
By incorporating the sources in Equations (8) and (9) into Equation (7), we arrive at a conclusion:
k = 0 ϵ k V k ( , ϑ ) = V ( 0 ) + V ( 0 ) + ϵ × E 1 u ß E { P 1 k = 0 ϵ k V k ( , ϑ ) + k = 0 ϵ k H k ( V ) } .
Upon evaluating the ϵ coefficient, we obtain
ϵ 0 : V 0 ( , ϑ ) = V ( 0 ) + V ( 0 ) , ϵ 1 : V 1 ( , ϑ ) = E 1 u ß E ( P 1 [ ] V 0 ( , ϑ ) + H 0 ( V ) ) , ϵ 2 : V 2 ( , ϑ ) = E 1 u ß E ( P 1 [ ] V 1 ( , ϑ ) + H 1 ( V ) ) , ϵ k : V k ( , ϑ ) = E 1 u ß E ( P 1 [ ] V k 1 ( , ϑ ) + H k 1 ( V ) ) , k > 0 , k N .
Finally, we employ truncated series to estimate the analytical solution V k ( , ϑ ) :
V ( , ϑ ) = lim M k = 1 M V k ( , ϑ ) .
The solutions to this series tend to reach stability rapidly.

4. General Methodology of the YTDM

This section focuses on the core procedure for resolving partial differential equations of any arbitrary order through the utilization of the ETDM:
D ϑ ß V ( , ϑ ) = P 1 ( , ϑ ) + Q 1 ( , ϑ ) , 0 < ß 1 ,
subjected to initial conditions
V ( , 0 ) = ξ ( ) , ϑ V ( , 0 ) = ζ ( ) .
By utilizing the Elzaki transformation, it is possible to express the fractional Caputo derivative D ϑ ß as the fractional derivative with respect to ϑ ß . Additionally, P 1 and Q 1 are linear and nonlinear operators, respectively:
E [ D ϑ ß V ( , ϑ ) ] = E [ P 1 ( , ϑ ) + Q 1 ( , ϑ ) ] .
By utilizing the characteristic of the Elzaki transform that distinguishes it, we can deduce that
1 u ß { M ( u ) u V ( 0 ) u 2 V ( 0 ) } = E [ P 1 ( , ϑ ) + Q 1 ( , ϑ ) ] .
Regarding simplification, we obtain
M ( V ) = u V ( 0 ) + u 2 V ( 0 ) + u ß E [ P 1 ( , ϑ ) + Q 1 ( , ϑ ) ] .
By applying the inverse Elzaki transformation to both sides, we obtain
V ( , ϑ ) = V ( 0 ) + V ( 0 ) + E 1 [ u ß E [ P 1 ( , ϑ ) + Q 1 ( , ϑ ) ] .
By utilizing the ETDM, the solution to the infinite series of V ( , ϑ ) is obtained:
V ( , ϑ ) = m = 0 V m ( , ϑ ) .
The Adomian polynomials for nonlinear terms Q 1 can be expressed by
Q 1 ( , ϑ ) = m = 0 A m .
where
A m = 1 m ! m m Q 1 k = 0 k k , k = 0 k ϑ k = 0 .
By combining Equations (21) and (19) into Equation (18), we arrive at a new expression:
m = 0 V m ( , ϑ ) = V ( 0 ) + V ( 0 ) + E 1 u ß E P 1 ( m = 0 m , m = 0 ϑ m ) + m = 0 A m .
The following terms will be explained:
V 0 ( , ϑ ) = V ( 0 ) + ϑ V ( 0 ) ,
V 1 ( , ϑ ) = E 1 u ß E + { P 1 ( 0 , ϑ 0 ) + A 0 } .
In general, when m 1 , it can be stated that
V m + 1 ( , ϑ ) = E 1 u ß E + { P 1 ( m , ϑ m ) + A m } .

5. Applications

 Example 1.
Consider the fractional shock wave equation, given by
ß V ϑ ß + 1 c 0 Υ + 1 2 V c 0 2 V = 0 , 0 < ß 1 ,
with the initial condition
V ( , 0 ) = exp 2 2 .
We utilized the Elzaki transform and set the constant values of c 0 to 2 and Υ to 1.5, which correspond to the flow of air:
E ß V ϑ ß = E 1 c 0 Υ + 1 2 V c 0 2 V .
By utilizing the differentiation property of the Elzaki transform, it can be shown that
1 u ß { M ( u ) u V ( 0 ) } = E 1 c 0 Υ + 1 2 V c 0 2 V ,
M ( u ) = u V ( 0 ) u ß E 1 c 0 Υ + 1 2 V c 0 2 V .
Using the inverse Elzaki transformation, we obtain
V ( , ϑ ) = V ( 0 ) E 1 u ß E 1 c 0 Υ + 1 2 V c 0 2 V , V ( , ϑ ) = exp 2 2 E 1 u ß E 1 c 0 Υ + 1 2 V c 0 2 V .
Now, by applying the HPM, we achieve
k = 0 ϵ k V k ( , ϑ ) = γ + ϵ E 1 u ß E 1 c 0 V k = 0 ϵ k V k ( , ϑ ) Υ + 1 2 1 c 0 2 k = 0 ϵ k H k ( V ) .
where the nonlinear terms shown by He’s polynomial H k ( V ) , represented by
k = 0 ϵ k H k ( V ) = V V .
A few components of He’s polynomial are determined as follows:
H 0 ( V ) = V 0 V 0 , H 1 ( V ) = V 0 V 1 + V 1 V 0
By comparing both sides’ coefficients of ϵ, we achieve
ϵ 0 : V 0 ( , ϑ ) = exp 2 2 , ϵ 1 : V 1 ( , ϑ ) = E 1 u ß E 1 c 0 V 0 Υ + 1 2 1 c 0 2 H 0 ( V ) = 1 c 0 Υ + 1 2 c 0 2 exp 2 2 exp 2 2 ϑ ß Γ ( ß + 1 ) , ϵ 2 : V 2 ( , ϑ ) = E 1 u ß E 1 c 0 V 1 Υ + 1 2 1 c 0 2 H 1 ( V ) = exp 2 2 [ 1 c 0 2 + 2 c 0 2 Υ + 1 c 0 3 exp 2 2 2 ( Υ + 1 ) c 0 3 2 exp 2 2 ( Υ + 1 ) 2 4 c 0 4 exp ( 2 ) + 3 ( Υ + 1 ) 2 4 c 0 4 2 exp ( 2 ) ] ϑ 2 ß Γ ( 2 ß + 1 ) ,
Some element of the HPTM result can be obtained:
V ( , ϑ ) = exp 2 2 1 c 0 Υ + 1 2 c 0 2 exp 2 2 exp 2 2 ϑ ß Γ ( ß + 1 ) + exp 2 2 [ 1 c 0 2 + 2 c 0 2 Υ + 1 c 0 3 exp 2 2 2 ( Υ + 1 ) c 0 3 2 exp 2 2 ( Υ + 1 ) 2 4 c 0 4 exp ( 2 ) + 3 ( Υ + 1 ) 2 4 c 0 4 2 exp ( 2 ) ] ϑ 2 ß Γ ( 2 ß + 1 ) + .
Now, We Use the ETDM
The Elzaki decomposition method is a mathematical technique used to decompose a given function into a set of simpler functions. The Elzaki decomposition method involves expressing the original function as a sum of the auxiliary functions, each of which satisfies a simpler differential equation. The auxiliary functions used in the Elzaki decomposition method are chosen so that they form a complete set, meaning that any function can be expressed as a linear combination of these functions. Once the auxiliary functions have been determined, the coefficients of the linear combination can be found by applying the boundary conditions or initial conditions to the original function.
By applying the Elzaki transformation in Equation (24), we obtain
E ß V ϑ ß = E 1 c 0 Υ + 1 2 V c 0 2 V .
By applying the Elzaki transformation differentiation property, we have
1 u ß { M ( u ) u V ( 0 ) } = E 1 c 0 Υ + 1 2 V c 0 2 V ,
M ( u ) = u V ( 0 ) + u ß E 1 c 0 Υ + 1 2 V c 0 2 V .
By applying the inverse Elzaki transformation, we achieve
V ( , ϑ ) = V ( 0 ) E 1 u ß E 1 c 0 Υ + 1 2 V c 0 2 V , V ( , ϑ ) = exp 2 2 E 1 u ß E 1 c 0 Υ + 1 2 V c 0 2 V .
The series form solution of V ( , ϑ ) is given by
V ( , ϑ ) = m = 0 V m ( , ϑ ) .
The nonlinear functions found with the help of the Adomian polynomials of V V = m = 0 A m are given by
m = 0 V m ( , ϑ ) = V ( , 0 ) E 1 u ß E 1 c 0 V Υ + 1 2 c 0 2 m = 0 A m , m = 0 V m ( , ϑ ) = exp 2 2 E 1 u ß E 1 c 0 V Υ + 1 2 c 0 2 m = 0 A m .
In same technique, the nonlinear terms are expressed as follows:
A 0 = V 0 V 0 , A 1 = V 0 V 1 + V 1 V 0 .
When comparing both sides, we obtain
V 0 ( , ϑ ) = exp 2 2 ,
For m = 0 , we have
V 1 ( , ϑ ) = 1 c 0 Υ + 1 2 c 0 2 exp 2 2 exp 2 2 ϑ ß Γ ( ß + 1 ) .
For m = 1 , we have
V 2 ( , ϑ ) = exp 2 2 [ 1 c 0 2 + 2 c 0 2 Υ + 1 c 0 3 exp 2 2 2 ( Υ + 1 ) c 0 3 2 exp 2 2 ( Υ + 1 ) 2 4 c 0 4 exp ( 2 ) + 3 ( Υ + 1 ) 2 4 c 0 4 2 exp ( 2 ) ] ϑ 2 ß Γ ( 2 ß + 1 ) .
The remaining steps of ß m for the ETDM for ( m 2 ) are simple to obtain. Hence, we write the series type result as follows:
V ( , ϑ ) = m = 0 V m ( , ϑ ) = V 0 ( , ϑ ) + V 1 ( , ϑ ) + V 2 ( , ϑ ) + .
V ( , ϑ ) = exp 2 2 1 c 0 Υ + 1 2 c 0 2 exp 2 2 exp 2 2 ϑ ß Γ ( ß + 1 ) + exp 2 2 [ 1 c 0 2 + 2 c 0 2 Υ + 1 c 0 3 exp 2 2 2 ( Υ + 1 ) c 0 3 2 exp 2 2 ( Υ + 1 ) 2 4 c 0 4 exp ( 2 ) + 3 ( Υ + 1 ) 2 4 c 0 4 2 exp ( 2 ) ] ϑ 2 ß Γ ( 2 ß + 1 ) + .
In Figure 1, The solutions for the V ( , ϑ ) surface with regard to ℘ and ϑ are explored for various values of ß. (a) The surface solution of the function V ( , ϑ ) at the value of ß = 0.25 , (b) ß = 0.50 , (c) ß = 0.75 , and (d) ß = 1 . Each of these solutions provides insight into how the surface of V ( , ϑ ) changes as the value of ß is altered.

Example

Consider the fractional shock wave equation, given by
ß V ϑ ß + V V ϑ V 2 2 = 0 , 0 < ß 1 ,
with the initial condition
V ( , 0 ) = 3 sec h 2 15 2 .
By using the Elzaki transformation, we obtain
E ß V ϑ ß = E ϑ V 2 2 V V .
By applying the Elzaki transformation differentiation property, we obtain
1 u ß { M ( u ) u V ( 0 ) } = E ϑ V 2 2 V V ,
M ( u ) = u V ( 0 ) + u ß E ϑ V 2 2 V V .
By performing the inverse Elzaki transformation, we obtain
V ( , ϑ ) = V ( 0 ) + E 1 u ß E ϑ V 2 2 V V , V ( , ϑ ) = 3 sec h 2 15 2 + E 1 u ß E ϑ V 2 2 V V .
Using the HPM, we now gain access to
k = 0 ϵ k V k ( , ϑ ) = 3 sec h 2 15 2 + ϵ E 1 u ß E k = 0 ϵ k V k ( , ϑ ) ϑ k = 0 ϵ k H k ( V ) .
In the representation of nonlinear terms, He’s polynomial H k ( V ) is used. This equation states that the nonlinear terms can be expressed through He’s polynomial:
k = 0 ϵ k H k ( V ) = V V .
The calculation of certain components of He’s polynomials can be performed as follows:
H 0 ( V ) = V 0 V 0 , H 1 ( V ) = V 0 V 1 + V 1 V 0 .
By comparing the value of the ϵ coefficient, we are able to deduce the following:
ϵ 0 : V 0 ( , ϑ ) = 3 sec h 2 15 2 , ϵ 1 : V 1 ( , ϑ ) = E 1 u ß E ϑ V 0 2 2 H 0 ( V ) = 9 sec h 4 15 2 tanh 15 2 ϑ ß Γ ( ß + 1 ) , ϵ 2 : V 2 ( , ϑ ) = E 1 u ß E ϑ V 1 2 2 H 1 ( V ) = [ 189 2 sec h 6 15 2 tanh 2 15 2 27 2 sec h 6 15 2 ] ϑ 2 ß Γ ( 2 ß + 1 ) + 135 2 sec h 4 15 2 tanh 3 15 2 63 2 sec h 4 15 2 tanh 15 2 Γ ( ß + 1 ) ϑ 2 ß 1 Γ ( 2 ß ) ,
To fully implement the HPTM solution, we need to persist in our approach:
V ( , ϑ ) = 3 sec h 2 15 2 + 9 sec h 4 15 2 tanh 15 2 ϑ ß Γ ( ß + 1 ) + [ 189 2 sec h 6 15 2 tanh 2 15 2 27 2 sec h 6 15 2 ] ϑ 2 ß Γ ( 2 ß + 1 ) + [ 135 2 sec h 4 15 2 tanh 3 15 2 63 2 sec h 4 15 2 tanh 15 2 ] Γ ( ß + 1 ) ϑ 2 ß 1 Γ ( 2 ß ) + .
Now, We Use the ETDM
By utilizing the Elzaki transformation method, we obtain
E ß V ϑ ß = E ϑ V 2 2 V V .
By utilizing the differentiation property of the Elzaki transform, it is possible to obtain
1 u ß { M ( u ) u V ( 0 ) } = E ϑ V 2 2 V V ,
M ( u ) = u V ( 0 ) + u ß E ϑ V 2 2 V V .
By applying the inverse Elzaki transform to both sides, we obtain
V ( , ϑ ) = V ( 0 ) + E 1 u ß E ϑ V 2 2 V V , V ( , ϑ ) = 1 ( 1 + exp ) 2 + E 1 u ß E ϑ V 2 2 V V .
With the use of the ETDM, an infinite series solution for V ( , ϑ ) can be obtained:
V ( , ϑ ) = m = 0 V m ( , ϑ ) .
Using the Adomian polynomials, the nonlinear terms are obtained through V V = m = 0 A m :
m = 0 V m ( , ϑ ) = V ( , 0 ) + E 1 u ß E ϑ V 2 2 m = 0 A m , m = 0 V m ( , ϑ ) = 1 ( 1 + exp ) 2 + E 1 u ß E ϑ V 2 2 m = 0 A m .
The nonlinear terms are determined as follows:
A 0 = V 0 V 0 , A 1 = V 0 V 1 + V 1 V 0 .
When comparing both sides, we obtain
V 0 ( , ϑ ) = 3 sec h 2 15 2 ,
For m = 0 , we have
V 1 ( , ϑ ) = 9 sec h 4 15 2 tanh 15 2 ϑ ß Γ ( ß + 1 ) .
For m = 1 , we have
V 2 ( , ϑ ) = 189 2 sec h 6 15 2 tanh 2 15 2 27 2 sec h 6 15 2 ϑ 2 ß Γ ( 2 ß + 1 ) + [ 135 2 sec h 4 15 2 tanh 3 15 2 63 2 sec h 4 15 2 tanh 15 2 ] Γ ( ß + 1 ) ϑ 2 ß 1 Γ ( 2 ß ) .
As a result, finding the remaining elements ß m of the ETDM for values of m 3 is straightforward. Therefore, we can express the solution in series form:
V ( , ϑ ) = m = 0 V m ( , ϑ ) = V 0 ( , ϑ ) + V 1 ( , ϑ ) +
V ( , ϑ ) = 3 sec h 2 15 2 + 9 sec h 4 15 2 tanh 15 2 ϑ ß Γ ( ß + 1 ) + [ 189 2 sec h 6 15 2 tanh 2 15 2 27 2 sec h 6 15 2 ] ϑ 2 ß Γ ( 2 ß + 1 ) + [ 135 2 sec h 4 15 2 tanh 3 15 2 63 2 sec h 4 15 2 tanh 15 2 ] Γ ( ß + 1 ) ϑ 2 ß 1 Γ ( 2 ß ) + .
At ß = 1 , we find as the closed type the exact result:
V ( , ϑ ) = 3 sec h 2 15 ϑ 2 .
In Figure 2, The solutions for the V ( , ϑ ) surface with regard to and ϑ are explored for various values of ß. (a) The surface solution of the function V ( , ϑ ) at the values of ß = 0.25 , (b) ß = 0.50 , (c) ß = 0.75 , and (d) ß = 1 . In Table 1, the exact ETDM and HPTM results of V ( , ϑ ) for problem 1 at various fractional orders of ß having different values of and ϑ .

6. Conclusions

In this study, an effective solution approach is presented for solving nonlinear fractional shock wave equations. The proposed method combines the Elzaki transform, homotopy perturbation method, and Adomian decomposition method. The Elzaki transform is employed to convert the nonlinear fractional shock wave equation into an equivalent integral equation. Then, the homotopy perturbation method and Adomian decomposition method are used to approximate the solution of the integral equation. Numerical examples were included to demonstrate the effectiveness of the proposed approach, and the obtained results were compared with those obtained by other existing methods. The findings reveal that the combined method yielded accurate and efficient solutions for nonlinear fractional shock wave equations. In summary, the proposed approach can be a promising technique for solving nonlinear fractional shock wave equations and can be used in various scientific and engineering applications.

Funding

This research received no external funding.

Data Availability Statement

The numerical data used to support the findings of this study are included within the article.

Conflicts of Interest

The author declares that there is no conflict of interest regarding the publication of this article.

References

  1. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach Science Publishers: Amsterdam, The Netherlands, 1993. [Google Scholar]
  2. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley: Hoboken, NJ, USA, 1993. [Google Scholar]
  3. Oldham, K.B.; Spanier, J. The Fractional Calculus: Theory and Applications of Differentiation and Integration to Arbitrary Order; Academic Press: Cambridge, MA, USA, 1974. [Google Scholar]
  4. Mofarreh, F.; Khan, A.; Shah, R. A Comparative Analysis of Fractional-Order Fokker-Planck Equation. Symmetry 2023, 15, 430. [Google Scholar] [CrossRef]
  5. Naeem, M.; Yasmin, H.; Shah, N.A.; Nonlaopon, K. Investigation of Fractional Nonlinear Regularized Long-Wave Models via Novel Techniques. Symmetry 2023, 15, 220. [Google Scholar] [CrossRef]
  6. Alshehry, A.S.; Shah, N.A.; Shah, R.; Nonlaopon, K. Implementation of Yang residual power series method to solve fractional non-linear systems. AIMS Math. 2023, 8, 8294–8309. [Google Scholar] [CrossRef]
  7. Alderremy, A.A.; Shah, R.; Aly, S.; Nonlaopon, K. Comparison of two modified analytical approaches for the systems of time fractional partial differential equations. AIMS Math. 2023, 8, 7142–7162. [Google Scholar] [CrossRef]
  8. Machado, J.A.; Tenreiro, C. An Introduction to the Numerical Solution of Fractional Differential Equations; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  9. Xie, Z.; Feng, X.; Chen, X. Partial Least Trimmed Squares Regression. Chemom. Intell. Lab. Syst. 2022, 221, 104486. [Google Scholar] [CrossRef]
  10. Li, X.; Dong, Z.; Wang, L.; Niu, X.; Yamaguchi, H.; Li, D.; Yu, P. A Magnetic Field Coupling Fractional Step Lattice Boltzmann Model for the Complex Interfacial Behavior in Magnetic Multiphase Flows. Appl. Math. Model. 2023, 117, 219–250. [Google Scholar] [CrossRef]
  11. Xu, S.; Dai, H.; Feng, L.; Chen, H.; Chai, Y.; Zheng, W.X. Fault Estimation for Switched Interconnected Nonlinear Systems with External Disturbances via Variable Weighted Iterative Learning. IEEE Trans. Circuits Syst. II Express Briefs 2023. [Google Scholar] [CrossRef]
  12. Sun, L.; Hou, J.; Xing, C.; Fang, Z. A Robust Hammerstein-Wiener Model Identification Method for Highly Nonlinear Systems. Processes 2022, 10, 2664. [Google Scholar] [CrossRef]
  13. Wang, L.; Liu, G.; Xue, J.; Wong, K. Channel Prediction Using Ordinary Differential Equations for MIMO Systems. IEEE Trans. Veh. Technol. 2022, 72, 2111–2119. [Google Scholar] [CrossRef]
  14. Burgers, J.M. A mathematical model illustrating the theory of turbulence. Adv. Appl. Mech. 1939, 1, 171–199. [Google Scholar]
  15. Kolmogorov, A.N. The local structure of turbulence in incompressible viscous fluid for very large Reynolds’ numbers. Dokl. Akad. Nauk SSSR 1941, 30, 299–303. [Google Scholar]
  16. Courant, R.; Friedrichs, K.; Lewy, H. On the partial difference equations of mathematical physics. IBM J. Res. Dev. 1952, 11, 215–234. [Google Scholar] [CrossRef]
  17. Bona, J.L.; Chen, M.; Saut, J.C. Boussinesq equations and other systems for small-amplitude long waves in nonlinear dispersive media. II: Nonlinear theory. Nonlinearity 2002, 15, 1–44. [Google Scholar] [CrossRef] [Green Version]
  18. Kirane, M.; Abdeljabbar, A. Nonexistence of Global Solutions of Systems of Time Fractional Differential equations posed on the Heisenberg group. Math. Methods Appl. Sci. 2022, 45, 7336–7345. [Google Scholar] [CrossRef]
  19. Rahman, Z.; Abdeljabbar, A.; Ali, M.Z. Novel Precise Solitary Wave Solutions of Two Time Fractional Nonlinear Evolution Models via the MSE Scheme. Fractal Fract. 2022, 6, 444. [Google Scholar] [CrossRef]
  20. Xie, X.; Wang, T.; Zhang, W. Existence of Solutions for the (p,q)-Laplacian Equation with Nonlocal Choquard Reaction. Appl. Math. Lett. 2023, 135, 108418. [Google Scholar] [CrossRef]
  21. Xu, K.; Guo, Y.; Liu, Y.; Deng, X.; Chen, Q.; Ma, Z. 60-GHz Compact Dual-Mode On-Chip Bandpass Filter Using GaAs Technology. IEEE Electron Device Lett. 2021, 42, 1120–1123. [Google Scholar] [CrossRef]
  22. Liu, Y.; Xu, K.; Li, J.; Guo, Y.; Zhang, A.; Chen, Q. Millimeter-Wave E-Plane Waveguide Bandpass Filters Based on Spoof Surface Plasmon Polaritons. IEEE Trans. Microw. Theory Tech. 2022, 70, 4399–4409. [Google Scholar] [CrossRef]
  23. Liu, K.; Yang, Z.; Wei, W.; Gao, B.; Xin, D.; Sun, C.; Wu, G. Novel Detection Approach for Thermal Defects: Study on Its Feasibility and Application to Vehicle Cables. High Voltage 2022, 1–10. [Google Scholar] [CrossRef]
  24. Lu, S.; Yang, B.; Xiao, Y.; Liu, S.; Liu, M.; Yin, L.; Zheng, W. Iterative Reconstruction of Low-Dose CT Based on Differential Sparse. Biomed. Signal Process. Control 2023, 79, 104204. [Google Scholar] [CrossRef]
  25. Shen, M.; Li, Y.; Li, X.; Zhang, X. The analytical solutions and numerical simulations of the Burgers equation with variable coefficients. Symmetry 2019, 11, 1267. [Google Scholar] [CrossRef]
  26. Alesemi, M.; Shahrani, J.S.A.; Iqbal, N.; Shah, R.; Nonlaopon, K. Analysis and Numerical Simulation of System of Fractional Partial Differential Equations with Non-Singular Kernel Operators. Symmetry 2023, 15, 233. [Google Scholar] [CrossRef]
  27. Kim, H.; Lee, S.; Kim, J. The numerical solutions of the nonlinear shock wave equation with boundary conditions. J. Comput. Phys. 2019, 396, 490–526. [Google Scholar] [CrossRef]
  28. Wang, S.; Zhou, Y.; Sun, L. The comparison of finite difference and finite volume methods for the numerical solutions of the nonlinear shock wave equation. Appl. Math. Comput. 2020, 365, 180–194. [Google Scholar] [CrossRef]
  29. Zhang, X.; Liu, Z.; Xu, Y. The analysis of the stability of the shock wave solutions of the nonlinear wave equation. J. Math. Anal. Appl. 2020, 489, 326–341. [Google Scholar] [CrossRef]
  30. Liu, J.; Yu, Y. Numerical Solutions of the Fractional Nonlinear Shock Wave Equation. Mathematics 2021, 9, 23. [Google Scholar] [CrossRef]
  31. Zhang, Q.; Liu, J. Analytical solutions for the fractional nonlinear shock wave equation. J. Appl. Math. Phys. 2020, 8, 2662–2668. [Google Scholar] [CrossRef]
  32. Wang, H.; Zhang, Y. A numerical solution for fractional nonlinear shock wave equation by Adomian decomposition method. J. Appl. Math. Phys. 2019, 7, 2057–2064. [Google Scholar] [CrossRef]
  33. Lyu, W.; Wang, Z. Logistic damping effect in chemotaxis models with density-suppressed motility. Adv. Nonlinear Anal. 2023, 12. [Google Scholar] [CrossRef]
  34. Li, Z.; Sun, Y. Study on the fractional nonlinear shock wave equation by the Adomian decomposition method. J. Appl. Math. Phys. 2017, 5, 2682–2687. [Google Scholar] [CrossRef]
  35. He, J.H. A new perturbation method. Comput. Methods Appl. Mech. Eng. 1999, 178, 257–262. [Google Scholar] [CrossRef]
  36. Elzaki, T.M. A novel nonlinear transform method for solving linear and nonlinear partial differential equations. J. Math. Anal. Appl. 2010, 365, 397–407. [Google Scholar]
  37. Chen, Y. Applications of homotopy perturbation method to nonlinear problems. Nonlinear Sci. Lett. A 2015, 6, 259–263. [Google Scholar]
  38. Lu, J.Y.; Chen, G.R. The combination of homotopy perturbation method and Elzaki transform for solving nonlinear problems. J. Math. Anal. Appl. 2016, 433, 738–747. [Google Scholar]
  39. Elzaki, T.M.; Mustafa, M.M. New method for solving nonlinear differential equations. J. Appl. Math. 2011. [Google Scholar]
  40. Khodaei, S.; Aligholi, M.; Niroom, F. Application of the Adomian decomposition method and Elzaki transform for solving nonlinear differential equations. J. Appl. Math. Phys. 2015, 3, 757–764. [Google Scholar]
  41. Elzaki, T.M. The new integral transform Elzaki transform. Glob. J. Pure Appl. Math. 2011, 7, 57–64. [Google Scholar]
  42. Elzaki, T.M.; Hilal, E.M.; Arabia, J.S.; Arabia, J.S. Homotopy perturbation and Elzaki transform for solving nonlinear partial differential equations. Math. Theory Model. 2012, 2, 33–42. [Google Scholar]
Figure 1. The solutions for the V ( , ϑ ) surface with regard to and ϑ are explored for various values of ß. (a) The surface solution of the function V ( , ϑ ) at the value of ß = 0.25 , (b) ß = 0.50 , (c) ß = 0.75 , and (d) ß = 1 . Each of these solutions provides insight into how the surface of V ( , ϑ ) changes as the value of ß is altered.
Figure 1. The solutions for the V ( , ϑ ) surface with regard to and ϑ are explored for various values of ß. (a) The surface solution of the function V ( , ϑ ) at the value of ß = 0.25 , (b) ß = 0.50 , (c) ß = 0.75 , and (d) ß = 1 . Each of these solutions provides insight into how the surface of V ( , ϑ ) changes as the value of ß is altered.
Mathematics 11 01253 g001
Figure 2. The solutions for the V ( , ϑ ) surface with regard to and ϑ are explored for various values of ß. (a) The surface solution of the function V ( , ϑ ) at the values of ß = 0.25 , (b) ß = 0.50 , (c) ß = 0.75 , and (d) ß = 1 .
Figure 2. The solutions for the V ( , ϑ ) surface with regard to and ϑ are explored for various values of ß. (a) The surface solution of the function V ( , ϑ ) at the values of ß = 0.25 , (b) ß = 0.50 , (c) ß = 0.75 , and (d) ß = 1 .
Mathematics 11 01253 g002
Table 1. The exact E T D M and H P T M results of V ( , ϑ ) for problem 1 at various fractional orders of ß having different values of and ϑ .
Table 1. The exact E T D M and H P T M results of V ( , ϑ ) for problem 1 at various fractional orders of ß having different values of and ϑ .
ϑ ß = 0.6 ß = 0.8 ß = 1 ( Approx ) ß = 1 ( Approx ) ß = 1 ( Exact )
0.20.0000044900.0000044870.0000044830.0000044830.000004479
0.40.0000054810.0000054790.0000054760.0000054760.000005470
0.010.60.0000066930.0000066910.0000066880.0000066880.000006681
0.80.0000081710.0000081710.0000081690.0000081690.000008161
10.0000099820.0000099800.0000099780.0000099780.000009968
0.20.0000044900.0000044870.0000044830.0000044830.000004474
0.40.0000054820.0000054780.0000054760.0000054760.000005465
0.020.60.0000066930.0000066900.0000066880.0000066880.000006675
0.80.0000081770.0000081720.0000081690.0000081690.000008153
10.0000099870.0000099820.0000099780.0000099780.000009958
0.20.0000044890.0000044850.0000044830.0000044830.000004470
0.40.0000054890.000005830.0000054760.0000054760.000005459
0.030.60.0000066930.0000066900.0000066880.0000066880.000006668
0.80.0000081790.0000081740.0000081690.0000081690.000008145
10.0000099890.0000099820.0000099780.0000099780.000009948
0.20.0000044920.0000044870.0000044830.0000044830.000004465
0.40.0000054870.0000054810.0000054760.0000054760.000005454
0.040.60.0000066990.0000066910.0000066880.0000066880.000006661
0.80.0000081830.0000081770.0000081690.0000081690.000008136
10.0000099990.0000099830.0000099780.0000099780.000009938
0.20.0000044970.0000044890.0000044830.0000044830.000004461
0.40.0000054900.0000054790.0000054760.0000054760.000005448
0.050.60.0000066970.0000066910.0000066880.0000066880.000006655
0.80.0000081910.0000081850.0000081690.0000081690.000008128
10.0000099970.0000099880.0000099780.0000099780.000009928
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

Alesemi, M. An Innovative Approach to Nonlinear Fractional Shock Wave Equations Using Two Numerical Methods. Mathematics 2023, 11, 1253. https://doi.org/10.3390/math11051253

AMA Style

Alesemi M. An Innovative Approach to Nonlinear Fractional Shock Wave Equations Using Two Numerical Methods. Mathematics. 2023; 11(5):1253. https://doi.org/10.3390/math11051253

Chicago/Turabian Style

Alesemi, Meshari. 2023. "An Innovative Approach to Nonlinear Fractional Shock Wave Equations Using Two Numerical Methods" Mathematics 11, no. 5: 1253. https://doi.org/10.3390/math11051253

APA Style

Alesemi, M. (2023). An Innovative Approach to Nonlinear Fractional Shock Wave Equations Using Two Numerical Methods. Mathematics, 11(5), 1253. https://doi.org/10.3390/math11051253

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop