Next Article in Journal
Oracle-Preserving Latent Flows
Previous Article in Journal
Modulated Lacunary Statistical and Strong-Cesàro Convergences
Previous Article in Special Issue
An Improved Symmetric Numerical Approach for Systems of Second-Order Two-Point BVPs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Investigation of Fractional-Order Fornberg–Whitham Equations in the Framework of Aboodh Transformation

1
Department of Basic Sciences, Preparatory Year Deanship, King Faisal University, Al Ahsa 31982, Saudi Arabia
2
Department of Mathematics, Al-Zaytoonah University of Jordan, Amman 11733, Jordan
3
Department of Mathematics, Abdul Wali Khan University Mardan, Mardan 23200, Pakistan
4
Department of Physics, College of Science, Princess Nourah bint Abdulrahman University, Riyadh 11671, Saudi Arabia
5
Department of Physics, Faculty of Science, Port Said University, Port Said 42521, Egypt
6
Research Center for Physics (RCP), Department of Physics, Faculty of Science and Arts, Al-Mikhwah, Al-Baha University, Al-Baha 65431, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(7), 1353; https://doi.org/10.3390/sym15071353
Submission received: 20 May 2023 / Revised: 20 June 2023 / Accepted: 27 June 2023 / Published: 3 July 2023
(This article belongs to the Special Issue Numerical Analysis, Approximation Theory, Differential Equations)

Abstract

:
In this investigation, the fractional Fornberg–Whitham equation (FFWE) is solved and analyzed via the variational iteration method (VIM) and Adomian decomposition method (ADM) with the help of the Aboodh transformation (AT). The FFWE is an important model for describing several nonlinear wave propagations in various fields of science and plasma physics. The AT provides a powerful tool for transforming fractional-order differential equations (DEs) into integer-order ones, making them more amenable to analytical solutions. Accordingly, the main objective of this investigation is to demonstrate the effectiveness and accuracy of ADM and VIM in deriving some approximations for the FFWE. Furthermore, we highlight the advantages and potential applications of these methods in solving other fractional-order nonlinear problems in several scientific fields, especially in plasma physics and some engineering problems.

1. Introduction

The Fornberg–Whitham equation (FWE) holds significant physical meaning in the field of fluid dynamics and wave propagation. It is a partial differential equation (PDE) that combines dispersion and nonlinearity to describe the evolution of unidirectional wavepackets in dispersive media. Dispersion occurs when waves of different frequencies travel at different speeds, and this equation incorporates this effect, allowing for the examination of dispersive wavepacket evolution. Nonlinearity, which refers to the interaction between waves resulting in a wave that is not a simple superposition, is also included in the equation, enabling the study of wave interactions and complex wave phenomena [1,2]. This equation finds important applications in various physical systems. For example, in water waves, it can be used to analyze surface wave behavior, including wave stability, breaking, and the formation of rogue waves [3,4,5,6,7]. In optics, the FWE is applied to investigate the propagation of pulses and understand phenomena such as dispersion, self-focusing, and soliton formation. Overall, the Fornberg–Whitham equation provides valuable insights into the physical behavior of waves in dispersive media with nonlinear interactions, making it a crucial tool in the study of wave dynamics [8,9,10].
The fractional FWE (FFWE) is a nonlinear PDE (NLPDE) that describes the propagation of nonlinear structures in different mediums with fractional derivatives. This equation is a generalization of the classical FWE, which was widely used in the study of water waves, solitons, and many other nonlinear phenomena. The FFWE has attracted a lot of attention in recent years due to its wide range of applications in various fields such as physics, mathematics, engineering, and finance. The main reason for this is the fact that the FFWE can describe the behavior of a variety of complex systems more accurately than traditional models based on integer derivatives [11,12,13]. The relation between symmetry and the Fractional-Order Fornberg-Whitham equations lies in the preservation of certain symmetry properties in the fractional-order formulation. Symmetry is a fundamental concept in mathematics and physics, representing invariance under transformations. In the context of fractional calculus, the Fornberg-Whitham equations are extended to include fractional derivatives, which introduce additional degrees of freedom and complexity. The fractional-order Fornberg-Whitham equations can exhibit symmetry properties such as time-reversibility, translation invariance, or scaling invariance. These symmetries play a crucial role in understanding the behavior and properties of the equations, enabling the identification of specific solutions, conservation laws, and physical interpretations. By studying the symmetries of the fractional-order Fornberg-Whitham equations, researchers can gain insights into the underlying dynamics and develop efficient numerical methods or analytical techniques for their analysis.
The study of fractional differential equations (FFEs) has gained a lot of attention in recent years due to its numerous applications in various scientific fields, and physical and engineering problems [14,15,16,17,18,19]. In particular, the FFWE has been the subject of intense research in the last decade. A lot of progress has been made in the theoretical and numerical analysis of this equation, as well as in its applications. One of the earliest works on the FFWE was done by Deng and Li in 2010, where they studied the existence and uniqueness of solutions to the FFWE with initial and boundary conditions. Later, Liu and Anh investigated the well-posedness of the FFWE with different types of boundary conditions [20]. In 2013, Zhang and Deng introduced a finite difference scheme for solving the FFWE and proved its convergence and stability [21].
Several researchers have also studied the numerical simulation of the FFWE. For example, Liu et al. proposed a numerical scheme based on the finite volume method for solving the FFWE and demonstrated its accuracy and efficiency [22,23]. In 2018, Li et al. developed a spectral method for solving the FFWE and showed that it has higher accuracy and a higher convergence rate than other numerical methods [24]. The FFWE has also been applied in various fields, such as finance and image processing. For example, Hu et al. used the FFWE to model the dynamics of stock prices and showed that it can provide more accurate predictions than traditional models [25]. In 2015, Wang et al. used the FFWE to remove noise from digital images and demonstrated that it can achieve better results than other image denoising techniques [26]. In conclusion, the FFWE is a fascinating and important equation that has attracted a lot of attention from researchers in various fields. Theoretical and numerical studies of the FFWE have made significant progress in recent years, and its applications are growing rapidly [27,28,29,30,31].
The variational iteration transform method (VITM) is a powerful mathematical tool for solving nonlinear differential equations (NLDEs). The VITM was first introduced by Professor J.H. He in 2006, and since then it has been applied in many different fields of science and engineering. The VITM is based on the concept of a variational iteration, which involves constructing an auxiliary linear operator and using it to iteratively approach the solution of some NLDEs [32,33,34]. The VITM has several advantages over other methods for solving nonlinear differential equations. It is simple to use, computationally efficient, and can be applied to a wide range of NLDE. Additionally, the VITM can be used to obtain analytical solutions, which can provide valuable insights into the behavior of the system being studied. The VITM has been used to solve a variety of problems in science and engineering, including fluid dynamics, quantum mechanics, and finance [35,36,37]. It has also been used to model and analyze complex systems such as chaotic systems, fractional differential equations (FDEs), and systems with delay. Overall, the VITM is a powerful and versatile tool for solving NLDEs, and its wide range of applications makes it an important tool for researchers and engineers in many different fields [38,39,40].
The Adomian decomposition transform method (ADTM) and its family is a powerful mathematical tool used to solve NLDEs. This method was developed by George Adomian in the early 1980s and has since gained popularity due to its effectiveness in finding analytical solutions to complex problems. ADTM involves breaking down a NLDE into a series of linear equations, which are then solved sequentially. The solution is then obtained by summing the series of solutions obtained [41,42,43]. This method has been widely applied in various fields of science and engineering, including physics, chemistry, biology, and finance. It has proven to be an efficient alternative to traditional numerical methods as it provides closed-form solutions that are easy to interpret and can offer insights into the underlying dynamics of the problem being studied. In this article, we will discuss the principles and applications of the ADTM and explore some examples to illustrate its effectiveness in solving NLDEs [44,45,46,47]. Additionally, we will make a comparison between the approximations that will be obtained using the suggested methods and the exact solutions for the integer case in order to verify the accuracy of the used methods and their approximations. We will apply the suggested methods for analyzing two different types of Fornberg–Whitham equations (FWEs).
The rest of the paper is divided into the following sections: Section 2 briefly describes the Fundamental definitions. The general discussion of the proposed methods is presented in Section 3. Numerical implementations are presented in Section 4, and Section 5 gives the conclusion.

2. Fundamental Definitions

Definition 1. 
The Aboodh transformation (AT) of a function θ ( ß ) reads [48,49]
C = θ : | θ ( ß ) | < B e p j | ß | , if ß ( 1 ) i × [ 0 , ) , j = 1 , 2 ; B , p 1 , p 2 > 0 ,
can be defined as
A [ θ ( ß ) ] = M ( ψ ) ,
which leads to
A [ θ ( ß ) ] = 1 ψ 0 θ ( ß ) e ψ ß d ß = M ( ψ ) , p 1 ψ p 2 ,
Definition 2. 
The inverse of the AT of the function θ ( ß ) reads [48,49]
θ ( ß ) = A 1 [ M ( ψ ) ] .
Definition 3. 
Let θ ( ß ) E ; then, the Laplace transformation (LT) of the function θ ( ß ) reads [48,49]
L [ θ ( ß ) ] = 0 θ ( ß ) e s ß d ß = θ ( s ) .
Theorem 1. 
If θ ( ß ) C with the AT A [ θ ( ß ) ] and LT L [ θ ( ß ) ] is given as [48,49]
M ( ψ ) = 1 ψ θ ( ψ ) .
Definition 4. 
The Mittag–Leffler function frequently arises in the solution of fractional calculus and can be expressed as a special term [48,49]
E ( Z ) = ρ = 0 Z ρ Γ ( ρ + 1 ) , , Z C , R e ( ) 0 .
In the general form, we have
E , γ ξ = ρ = 0 Z ρ ( ξ ) ρ Γ ( γ + ρ ) ρ ! , , γ , Z C , R ( ) 0 , Re ( γ ) 0 .
Definition 5. 
The fractional AB derivative is a concept that pertains to the function θ H 1 ( 0 , 1 ) , where 0 < < 1 . The definition of the fractional AB derivative is as follows [48,49]:
0 A B C D ß θ ( ß ) = N ( ) 1 0 ß θ ( x ) E ( ß x ) 1 d x .
Definition 6. 
Let θ be an element in the Sobolev space H 1 ( 0 , 1 ) , and let 0 < ω < 1 ; then, the fractional AB derivatives can be defined using the Riemann-Liouville approach [48,49]
0 A B R D ß θ ( ß ) = N ( ) 1 d d ß 0 ß θ ( x ) E ( ß x ) 1 d x .
The normalization function N ( ) satisfies the requirement of being positive and has the values of N ( 0 ) = 1 and N ( 1 ) = 1 .
Theorem 2. 
The fractional AB operator of the LT in the presence of Caputo is given by [48,49]
L 0 A B C D ß θ ( ß ) = N ( ) 1 × s F ( s ) s 1 f ( 0 ) s + 1 .
Additionally, the LT of the fractional AB derivative when utilizing the Riemann-Liouville method reads
L 0 A B R D ß θ ( ß ) = N ( ) 1 × s F ( s ) s + 1 .
Theorem 3. 
If Ω , C , with Re ( ) > 0 , then the AT of E Ω ß is given by [48,49]
M E Ω ß = 1 ψ 2 1 Ω ψ 1 ,
where Ω ψ < 1 .
Theorem 4. 
The AT of ß γ 1 E , ξ ξ Ω ß can be represented by the following: let ℘ and γ be complex numbers such that their real part is positive, namely, R e ( ) > 0 and R e ( γ ) > 0 [48,49]
ß γ 1 E , ξ ξ Ω ß = 1 ψ γ + 1 1 Ω ψ ξ , Ω ψ < 1 .
Theorem 5. 
The AT of a fractional AB operator in the presence of Caputo can be defined as follows: if M ( ψ ) represents the AT of θ ( ß ) C and the LT of θ ( ß ) C is θ ( s ) [48,49]
M 0 A B C D ß θ ( ß ) = N ( ) M ( ψ ) ψ 2 θ ( 0 ) 1 + ψ .
Theorem 6. 
The AT of a fractional AB operator in the context of Riemann-Liouville is defined as follows: let M ( ψ ) represent the AT of θ ( ß ) , which is an element of C . Additionally, let θ ( s ) be the LT of θ ( ß ) C [48,49]
M 0 A B R D ß θ ( ß ) = N ( ) M ( ψ ) 1 + ψ .

3. The General Application of ADTM

Here, we will explore how the ADTM can be utilized for analyzing the following FPDEs
A B C D ß θ ( ξ , ß ) + G ¯ ( ξ , ß ) + N ( ξ , ß ) P ( ξ , ß ) = 0 , m 1 < m ,
the initial condition is
θ ( ξ , 0 ) = g ( ξ ) .
Let A B C D ß denote the Caputo fractional derivative of order with respect to ß, and let G ¯ and N represent the linear and nonlinear terms, respectively. Additionally, P is the source term.
Using AT to Equation (1), we obtain
A [ A B C D ß θ ( ξ , ß ) ] + A [ G ¯ ( ξ , ß ) + N ( ξ , ß ) P ( ξ , ß ) ] = 0 .
Applying AT to the property of differentiation, we obtain
A [ θ ( ξ , ß ) ] = 1 ϖ 2 θ ( ξ , 0 ) + 1 + ϖ N ( ) A [ P ( ξ , ß ) ] 1 + ϖ N ( ) A { G ¯ ( ξ , ß ) + N ( ξ , ß ) } ] .
The MDM result of infinite series θ ( ξ , ß ) ,
θ ( ξ , ß ) = ϕ = 0 θ ϕ ( ξ , ß ) ,
where N denotes the nonlinear function, which is defined by
N ( ξ , ß ) = ϕ = 0 A ϕ .
The nonlinear terms can be analyzed using Adomian polynomials. Thus, the expression for the Adomian polynomial formula is:
A ϕ = 1 j ! ϕ λ ϕ N ϕ = 0 λ ϕ θ ϕ λ = 0 .
Then, putting Equations (5) and (6) into (4) gives
A ϕ = 0 θ ϕ ( ξ , ß ) = 1 ϖ 2 θ ( ξ , 0 ) + 1 + ϖ N ( ) A { P ( ξ , ß ) } 1 + ϖ N ( ) A G ¯ ( ϕ = 0 θ ϕ ) + ϕ = 0 A ϕ .
By performing the inverse of the AT on Equation (8), we have
ϕ = 0 θ ϕ ( ξ , ß ) = A 1 1 ϖ 2 θ ( ξ , 0 ) + 1 + ϖ N ( ) A { P ( ξ , ß ) } 1 + ϖ N ( ) A G ¯ ϕ = 0 θ ϕ + ϕ = 0 A ϕ .
Introducing the following terms:
θ 0 ( ξ , ß ) = A 1 1 ϖ 2 θ ( ξ , 0 ) + 1 + ϖ N ( ) A { P ( ξ , ß ) } ,
θ 1 ( ξ , ß ) = A 1 1 + ϖ N ( ) A { G ¯ 1 ( θ 0 ) + A 0 } .
In general, for ϕ 1 , we obtain
θ ϕ + 1 ( ξ , ß ) = A 1 1 + ϖ N ( ) A { G ¯ ( θ ϕ ) + A ϕ } .

4. Numerical Results

Example 1. 
Consider the following nonlinear FFWE [50]:
D ß θ D ξ ξ ß θ + D ξ θ = θ D ξ ξ ξ θ θ D ξ θ + 3 D ξ θ D ξ ξ θ , 0 < 1 ,
with the initial condition
θ ( ξ , 0 ) = e ξ 2 .
Taking the AT to Equation (11), we obtain
N ( ) 1 + ϖ A [ θ ( ξ , ß ) ] 1 ϖ 2 θ ( ξ , 0 ) = A D ξ ξ ß θ D ξ θ + θ D ξ ξ ξ θ θ D ξ θ + 3 D ξ θ D ξ ξ θ .
Using the inverse of AT
θ ( ξ , ß ) = A 1 θ ( ξ , 0 ) ϖ 2 1 + ϖ N ( ) A D ξ ξ ß θ D ξ θ + θ D ξ ξ ξ θ θ D ξ θ + 3 D ξ θ D ξ ξ θ .
Applying Adomain procedure, we obtain
θ 0 ( ξ , ß ) = A 1 θ ( ξ , 0 ) ϖ 2 = A 1 e ξ 2 ϖ 2 ,
θ 0 ( ξ , ß ) = e ξ 2 ,
ϕ = 0 θ ϕ + 1 ( ξ , ß ) = A 1 1 + ϖ N ( ) A ϕ = 0 ( D ξ ξ ß θ ) ϕ ϕ = 0 ( D ξ θ ) ϕ + ϕ = 0 A ϕ ϕ = 0 B ϕ + 3 ϕ = 0 C ϕ , ϕ = 0 , 1 , 2 ,
A 0 ( θ D ξ ξ ξ θ ) = θ 0 D ξ ξ ξ θ 0 , B 0 ( θ D ξ θ ) = θ 0 D ξ θ 0 , A 1 ( θ D ξ ξ ξ θ ) = θ 0 D ξ ξ ξ θ 1 + θ 1 D ξ ξ ξ θ 0 , B 1 ( θ D ξ θ ) = θ 0 D ξ θ 1 + θ 1 D ξ θ 0 , A 2 ( θ D ξ ξ ξ θ ) = θ 1 D ξ ξ ξ θ 2 + θ 1 D ξ ξ ξ θ 1 + θ 2 D ξ ξ ξ θ 0 , B 2 ( θ D ξ θ ) = θ 1 D ξ θ 2 + θ 1 D ξ θ 1 + θ 2 D ξ θ 0 ,
C 0 ( D ξ θ D ξ ξ θ ) = D ξ θ 0 D ξ ξ θ 0 , C 1 ( D ξ θ D ξ ξ θ ) = D ξ θ 0 D ξ ξ θ 1 + D ξ θ 1 D ξ ξ θ 0 , C 2 ( D ξ θ D ξ ξ θ ) = D ξ θ 1 D ξ ξ θ 2 + D ξ θ 1 D ξ ξ θ 1 + D ξ θ 2 D ξ ξ θ 0 ,
for ϕ = 1
θ 1 ( ξ , ß ) = A 1 1 + ϖ N ( ) A D ξ ξ ß θ 0 D ξ θ 0 + A 0 B 0 + 3 C 0 = 1 2 e ξ 2 ( 1 ) + ß Γ ( + 1 ) .
for ϕ = 2
θ 2 ( ξ , ß ) = A 1 1 + ϖ N ( ) A D ξ ξ ß θ 1 D ξ θ 1 + A 1 B 1 + 3 C 1 , θ 2 ( ξ , ß ) = 1 8 e ξ 2 ß 2 1 Γ ( 2 ) + 1 4 e ξ 2 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) ,
for ϕ = 3
θ 3 ( ξ , ß ) = A 1 1 + ϖ N ( ) A D ξ ξ ß θ 2 D ξ θ 2 + A 2 B 2 + 3 C 2 , θ 3 ( ξ , ß ) = 1 32 e ξ 2 ß 3 2 Γ ( 3 1 ) + 1 8 e ξ 2 ß 3 1 Γ ( 3 ) 1 8 e ξ 2 ( 1 ) 3 + ( 1 ) ( 1 + + 2 2 ) ß Γ ( + 1 ) + 3 2 ( 1 ) ß 2 Γ ( 2 + 1 ) + 3 Γ ( 2 + 1 ) ß 3 Γ ( 3 + 1 ) .
The solution of example (1) using the MDM reads
θ ( ξ , ß ) = θ 0 ( ξ , ß ) + θ 1 ( ξ , ß ) + θ 2 ( ξ , ß ) + θ 3 ( ξ , ß ) + θ 4 ( ξ , ß ) + ,
θ ( ξ , ß ) = e ξ 2 1 2 e ξ 2 ( 1 ) + ß Γ ( + 1 ) 1 8 e ξ 2 ß 2 1 Γ ( 2 ) + 1 4 e ξ 2 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) 1 32 e ξ 2 ß 3 2 Γ ( 3 1 ) + 1 8 e ξ 2 3 1 Γ ( 3 ) 1 8 e ξ 2 { ( 1 ) 3 + ( 1 ) ( 1 + + 2 2 ) ß Γ ( + 1 ) + 3 2 ( 1 ) ß 2 Γ ( 2 + 1 ) + 3 Γ ( 2 + 1 ) ß 3 Γ ( 3 + 1 ) } .
The simplification form to Equation (17) reads
θ ( ξ , ß ) = e ξ 2 1 1 2 ( 1 ) + ß Γ ( + 1 ) 1 8 ß 2 1 Γ ( 2 ) + 1 4 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) 1 32 ß 3 2 Γ ( 3 1 ) + 1 8 ß 3 1 Γ ( 3 ) 1 8 { ( 1 ) 3 + ( 1 ) ( 1 + + 2 2 ) ß Γ ( + 1 ) + 3 2 ( 1 ) ß 2 Γ ( 2 + 1 ) + 3 Γ ( 2 + 1 ) ß 3 Γ ( 3 + 1 ) } + .
To obtain a solution in series form, the variational method can be applied.
By deriving the iteration formulas for Equation (11), we obtain:
θ ϕ + 1 ( ξ , ß ) = θ j ( ξ , ß ) A 1 1 + ϖ N ( ) A D ξ ξ ß θ ϕ + D ξ θ ϕ θ ϕ D ξ ξ ξ θ ϕ + θ ϕ D ξ θ ϕ 3 D ξ θ ϕ D ξ ξ θ ϕ ,
with
θ 0 ( ξ , ß ) = e ξ 2 .
For ϕ = 0 , 1 , 2 ,
θ 1 ( ξ , ß ) = θ 0 ( ξ , ß ) A 1 1 + ϖ N ( ) A D ξ ξ ß θ 0 + D ξ θ 0 θ 0 D ξ ξ ξ θ 0 + θ 0 D ξ θ 0 3 D ξ θ 0 D ξ ξ θ 0 , θ 1 ( ξ , ß ) = e ξ 2 1 2 e ξ 2 ( 1 ) + ß Γ ( + 1 ) ,
θ 2 ( ξ , ß ) = θ 1 ( ξ , ß ) A 1 1 + ϖ N ( ) A D ξ ξ ß θ 1 + D ξ θ 1 θ 1 D ξ ξ ξ θ 1 + θ 1 D ξ θ 1 3 D ξ θ 1 D ξ ξ θ 1 , θ 2 ( ξ , ß ) = e ξ 2 1 2 e ξ 2 ( 1 ) + ß Γ ( + 1 ) 1 8 e ξ 2 ß 2 1 Γ ( 2 ) + 1 4 e ξ 2 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) ,
θ 3 ( ξ , ß ) = θ 2 ( ξ , ß ) A 1 1 + ϖ N ( ) A D ξ ξ ß θ 2 + D ξ θ 2 θ 2 D ξ ξ ξ θ 2 + θ 2 D ξ θ 2 3 D ξ θ 2 D ξ ξ θ 2 , θ 3 ( ξ , ß ) = e ξ 2 1 2 e ξ 2 ( 1 ) + ß Γ ( + 1 ) 1 8 e ξ 2 ß 2 1 Γ ( 2 ) + 1 4 e ξ 2 ( 1 ) 2 + γ 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) 1 32 e ξ 2 ß 3 2 Γ ( 3 1 ) + 1 8 e ξ 2 ß 3 1 Γ ( 3 ) 1 8 e ξ 2 { ( 1 ) 3 + ( 1 ) ( 1 + + 2 2 ) ß Γ ( + 1 ) + 3 2 ( 1 ) ß 2 Γ ( 2 + 1 ) + 3 Γ ( 2 + 1 ) ß 3 Γ ( 3 + 1 ) } ,
θ ( ξ , ß ) = e ξ 2 1 2 e ξ 2 ( 1 ) + ß Γ ( + 1 ) 1 8 e ξ 2 ß 2 1 Γ ( 2 ) + 1 4 e ξ 2 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) 1 32 e ξ 2 ß 3 2 Γ ( 3 1 ) + 1 8 e ξ 2 3 1 Γ ( 3 ) 1 8 e ξ 2 { ( 1 ) 3 + ( 1 ) ( 1 + + 2 2 ) ß Γ ( + 1 ) + 3 2 ( 1 ) ß 2 Γ ( 2 + 1 ) + 3 Γ ( 2 + 1 ) ß 3 Γ ( 3 + 1 ) } .
The exact solution of Equation (11) at = 1 , reads
θ ( ξ , ß ) = e ξ 2 2 ß 3 .
In Figure 1, we make a comparison between both the approximate solution (18) or (24) at = 1 using ADTM/VITM and the exact solution (25) to Equation (11). Moreover, the absolute error for this case is estimated at t = 0.1 , as demonstrated in Table 1. Almost perfect agreement between the two solutions is noted, which confirms the high accuracy and high efficiency of the used approximate methods (ADTM/VITM). On the other hand, the obtained approximations (18) or (24) are simulated numerically at different values of the fractional-order , as shown in Figure 2. From the later figure, we can observe the effect of changing on the profile of the solutions. This graphical representation allows us to visually observe any changes in the shape, magnitude, or other properties of the solution as varies. These graphical discussions provide a visual representation of the solutions obtained using ADTM/VITM for Example 1 at different values of 1 . They help us to understand the impact of changing the fractional order on the behavior of the solution, allowing for a more comprehensive analysis of the problem.
Example 2. 
Consider the following nonlinear FFWE [51]:
D ß θ D ξ ξ ß θ + D ξ θ = θ D ξ ξ ξ θ θ D ξ θ + 3 D ξ θ D ξ ξ θ , ß > 0 , 0 < 1 ,
with the initial condition
θ ( ξ , 0 ) = cosh 2 ξ 4 .
Applying the AT to (26), we obtain
N ( ) 1 + ϖ A [ θ ( ξ , ß ) ] 1 ϖ 2 θ ( ξ , 0 ) = A D ξ ξ ß θ D ξ θ + θ D ξ ξ ξ θ θ D ξ θ + 3 D ξ θ D ξ ξ θ .
Using the inverse of AT
θ ( ξ , ß ) = A 1 θ ( ξ , 0 ) ϖ 2 1 + ϖ N ( ) A D ξ ξ ß θ D ξ θ + θ D ξ ξ ξ θ θ D ξ θ + 3 D ξ θ D ξ ξ θ .
Using the ADM procedure, we find
θ 0 ( ξ , ß ) = A 1 θ ( ξ , 0 ) ϖ 2 = A 1 exp cosh 2 ξ 4 ϖ 2 ,
θ 0 ( ξ , ß ) = cosh 2 ξ 4 ,
ϕ = 0 θ ϕ + 1 ( ξ , ß ) = A 1 1 + ϖ N ( ) A ϕ = 0 ( D ξ ξ ß θ ) ϕ ϕ = 0 ( D ξ θ ) ϕ + ϕ = 0 A ϕ ϕ = 0 B ϕ + 3 ϕ = 0 C ϕ , ϕ = 0 , 1 , 2 ,
for ϕ = 0
θ 1 ( ξ , ß ) = A 1 1 + ϖ N ( ) A D ξ ξ ß θ 0 D ξ θ 0 + A 0 B 0 + 3 C 0 = 11 32 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) .
for ϕ = 1
θ 2 ( ξ , ß ) = A 1 1 + ϖ N ( ) A D ξ ξ ß θ 1 D ξ θ 1 + A 1 B 1 + 3 C 1 , θ 2 ( ξ , ß ) = 11 28 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) + 121 1024 cosh ξ 4 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) ,
for ϕ = 2
θ 3 ( ξ , ß ) = A 1 1 + ϖ N ( ) A D ξ ξ ß θ 2 D ξ θ 2 + A 2 B 2 + 3 C 2 , θ 3 ( ξ , ß ) = 11 512 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) + 121 2048 cosh ξ 4 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) 1331 49,152 sinh ξ 4 ( 1 ) 3 + ( 1 ) ( 1 + + 2 2 ) ß Γ ( + 1 ) + 3 2 ( 1 ) ß 2 Γ ( 2 + 1 ) + 3 Γ ( 2 + 1 ) ß 3 Γ ( 3 + 1 ) ,
The MDM solution of example (2) reads
θ ( ξ , ß ) = θ 0 ( ξ , ß ) + θ 1 ( ξ , ß ) + θ 2 ( ξ , ß ) + θ 3 ( ξ , ß ) + θ 4 ( ξ , ß ) + ,
θ ( ξ , ß ) = cosh 2 ξ 4 11 32 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) 11 28 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) + 121 1024 cosh ξ 4 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) 11 512 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) + 121 2048 cosh ξ 4 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) 1331 49,152 sinh ξ 4 { ( 1 ) 3 + ( 1 ) ( 1 + + 2 2 ) ß Γ ( + 1 ) + 3 2 ( 1 ) ß 2 Γ ( 2 + 1 ) + 3 Γ ( 2 + 1 ) ß 3 Γ ( 3 + 1 ) } .
To obtain an analytical solution, the variational method can be employed.
We can obtain the iteration formulas for Equation (26), as follows:
θ ϕ + 1 ( ξ , ß ) = θ j ( ξ , ß ) A 1 1 + ϖ N ( ) A D ξ ξ ß θ ϕ + D ξ θ ϕ θ ϕ D ξ ξ ξ θ ϕ + θ ϕ D ξ θ ϕ 3 D ξ θ ϕ D ξ ξ θ ϕ ,
with
θ 0 ( ξ , ß ) = cosh 2 ξ 4 .
For ϕ = 0 , 1 , 2 ,
θ 1 ( ξ , ß ) = θ 0 ( ξ , ß ) A 1 1 + ϖ N ( ) A D ξ ξ ß θ 0 + D ξ θ 0 θ 0 D ξ ξ ξ θ 0 + θ 0 D ξ θ 0 3 D ξ θ 0 D ξ ξ θ 0 , θ 1 ( ξ , ß ) = cosh 2 ξ 4 11 32 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) ,
θ 2 ( ξ , ß ) = θ 1 ( ξ , ß ) A 1 1 + ϖ N ( ) A D ξ ξ ß θ 1 + D ξ θ 1 θ 1 D ξ ξ ξ θ 1 + θ 1 D ξ θ 1 3 D ξ θ 1 D ξ ξ θ 1 , θ 2 ( ξ , ß ) = cosh 2 ξ 4 11 32 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) 11 28 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) + 121 1024 cosh ξ 4 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) ,
θ 3 ( ξ , ß ) = θ 2 ( ξ , ß ) A 1 1 + ϖ N ( ) A D ξ ξ ß θ 2 + D ξ θ 2 θ 2 D ξ ξ ξ θ 2 + θ 2 D ξ θ 2 3 D ξ θ 2 D ξ ξ θ 2 , θ 3 ( ξ , ß ) = cosh 2 ξ 4 11 32 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) 11 28 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) + 121 1024 cosh ξ 4 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) 11 512 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) + 121 2048 cosh ξ 4 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) 1331 49,152 sinh ξ 4 { ( 1 ) 3 + ( 1 ) ( 1 + + 2 2 ) ß Γ ( + 1 ) + 3 2 ( 1 ) ß 2 Γ ( 2 + 1 ) + 3 Γ ( 2 + 1 ) ß 3 Γ ( 3 + 1 ) } ,
θ ( ξ , ß ) = cosh 2 ξ 4 11 32 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) 11 28 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) + 121 1024 cosh ξ 4 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) 11 512 sinh ξ 4 ( 1 ) + ß Γ ( + 1 ) + 121 2048 cosh ξ 4 ( 1 ) 2 + 2 ß 2 Γ ( 2 + 1 ) + 2 ( 1 ) ß Γ ( + 1 ) 1331 49,152 sinh ξ 4 { ( 1 ) 3 + ( 1 ) ( 1 + + 2 2 ) ß Γ ( + 1 ) + 3 2 ( 1 ) ß 2 Γ ( 2 + 1 ) + 3 Γ ( 2 + 1 ) ß 3 Γ ( 3 + 1 ) } + .
We put = 1 ; then, the series solution is given as
θ ( ξ , ß ) = cosh 2 ξ 4 11 32 sinh ξ 4 ß 11 28 sinh ξ 4 ß + 121 1024 cosh ξ 4 2 2 11 512 sinh ξ 4 + 121 2048 cosh ξ 4 2 2 1331 49,152 sinh ξ 4 3 3 + .
The exact result of Equation (26) at = 1 reads
θ ( ξ , ß ) = cosh 2 ξ 4 11 ß 24 .
In Figure 3, a comparison between both the approximate solution (32) or (38) at = 1 using ADTM/VITM and the exact solution (40) to Equation (26) is presented. Further, the absolute error for this case is estimated at t = 0.1 , as illustrated in Table 2. It is noted from the comparison results that the two solutions are compatible to a large degree, which enhances the accuracy and efficiency of the used numerical techniques (ADTM/VITM). Additionally, the obtained approximations (32) or (38) are displayed in Figure 4 at different values of the fractional-order . This graphical representation allows us to visually observe any changes in the shape, magnitude, or other properties of the solution as varies. These graphical discussions provide a visual representation of the solutions obtained using ADTM/VITM for Example 2 at different values of . These figures demonstrate how the wavepacket evolves and interacts within the dispersive medium according to the Fornberg–Whitham equation. Comparing this figure with the previous one at = 1 , noticeable changes in the wavepacket’s behavior can be observed. The alterations might include variations in amplitude, shape, and propagation speed. These changes indicate the influence of the fractional order on the dynamics of the wavepacket. By examining both Figure 3 and Figure 4, it becomes evident that altering the fractional order affects the characteristics of the solutions to Equation (26). The plots provide graphical insights into how changes in the fractional order parameter influence the wavepacket’s behavior, allowing for a better understanding of the physical significance of the equation.

5. Conclusions

In summing up, this study has presented a comprehensive comparison between the Adomian decomposition method (ADM) and the variational iteration method (VIM) for solving the fractional Fornberg–Whitham equation (FFWE) in the sense of the Aboodh transformation. Both methods have demonstrated their effectiveness and accuracy in solving this complex nonlinear fractional partial differential equation (FPDE). The ADM and VIM were employed to obtain approximate solutions of the FFWE, taking into account their inherent advantages and limitations. The ADM, a semi-analytical method, was shown to provide a systematic approach to constructing the solution in the form of a convergent series, while the VIM, an iterative method, provided a more straightforward and simpler approach for obtaining approximate solutions. A thorough analysis of the numerical results revealed that both methods have the potential to deliver accurate and reliable solutions to the FFWE. However, the choice of the most appropriate method may depend on the specific problem at hand, the computational resources, and the desired level of accuracy. Future research could focus on extending the application of these methods to other types of FPDEs, as well as exploring the combination of these methods with other numerical or analytical techniques to enhance their efficiency and accuracy. Furthermore, the implementation of these methods in high-performance computing environments could be investigated to tackle more complex problems and improve computational efficiency.

Author Contributions

Methodology, S.N.; Validation, M.A.H.; Formal analysis, R.S. and A.W.A.; Investigation, S.A.E.-T.; Resources, S.N. and R.S.; Data curation, S.A.E.-T.; Funding acquisition, S.N. All authors have read and agreed to the published version of the manuscript.

Funding

The authors express their gratitude to Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R378), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. This work was supported by the Deanship of Scientific Research, the Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (Grant No. 3739).

Data Availability Statement

The numerical data used to support the findings of this study are included within the article. Mathematica codes for drawing the figures are available, which can be requested from El-Tantawy.

Conflicts of Interest

The authors declare that there are no conflict of interest regarding the publication of this article.

References

  1. Johnson, R.S. Fornberg-Whitham equation. In Encyclopedia of Mathematics and Its Applications; Cambridge University Press: Cambridge, UK, 1997; Volume 60, pp. 35–37. [Google Scholar]
  2. Choi, W.; Camassa, R. Fully nonlinear internal waves in a two-fluid system. J. Fluid Mech. 2007, 581, 369–380. [Google Scholar] [CrossRef]
  3. He, H.M.; Peng, J.G.; Li, H.Y. Iterative approximation of fixed point problems and variational inequality problems on Hadamard manifolds. UPB Bull. Ser. A 2022, 84, 25–36. [Google Scholar]
  4. Xie, X.; Huang, L.; Marson, S.M.; Wei, G. Emergency response process for sudden rainstorm and flooding: Scenario deduction and Bayesian network analysis using evidence theory and knowledge meta-theory. Nat. Hazards 2023, 117, 3307–3329. [Google Scholar] [CrossRef]
  5. Jin, H.Y.; Wang, Z.A. Global stabilization of the full attraction-repulsion Keller-Segel system. Discret. Contin. Dyn. Syst. Ser. A 2020, 40, 3509–3527. [Google Scholar] [CrossRef] [Green Version]
  6. Guo, C.; Hu, J. Fixed-Time Stabilization of High-Order Uncertain Nonlinear Systems: Output Feedback Control Design and Settling Time Analysis. J. Syst. Sci. Complex. 2023. [Google Scholar] [CrossRef]
  7. Lyu, W.; Wang, Z. Global classical solutions for a class of reaction-diffusion system with density-suppressed motility. Electron. Res. Arch. 2022, 30, 995–1015. [Google Scholar] [CrossRef]
  8. Shah, N.A.; Hamed, Y.S.; Abualnaja, K.M.; Chung, J.D.; Khan, A. A comparative analysis of fractional-order kaup-kupershmidt equation within different operators. Symmetry 2022, 14, 986. [Google Scholar] [CrossRef]
  9. Ostrovsky, L.A.; Pelinovsky, E.N.; Shrira, V.I. Rogue waves in nonlinear dispersive media: Physical mechanisms, models, and applications. Phys. Rep. 2008, 443, 1–53. [Google Scholar]
  10. Stolen, R.H.; Gordon, J.P. Self-phase-modulation and small-scale filaments in nonlinear fibers. Opt. Lett. 1982, 7, 28–33. [Google Scholar]
  11. Fornberg, B.; Whitham, G.B. A numerical and theoretical study of certain nonlinear wave phenomena. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 1978, 289, 373–404. [Google Scholar]
  12. Fornberg, B.; Whitham, G.B. A numerical and theoretical study of certain nonlinear wave phenomena. II. Nonlinear geometrical optics. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 1979, 292, 385–409. [Google Scholar]
  13. Fornberg, B. Numerical solution of the Fornberg-Whitham equation. J. Comput. Phys. 1980, 36, 362–381. [Google Scholar]
  14. Zayed, E.M.; Rahman, H.M.A. On using the modified variational iteration method for solving the nonlinear coupled equations in the mathematical physics. Ric. Mat. 2010, 59, 137–159. [Google Scholar] [CrossRef]
  15. Zayed, E.M.E.; Nofal, T.A.; Gepreel, K.A. The travelling wave solutions for non-linear initial-value problems using the homotopy perturbation method. Int. J. Control. 2009, 88, 617–634. [Google Scholar] [CrossRef]
  16. Zhang, K.; Alshehry, A.S.; Aljahdaly, N.H.; Shah, N.A.; Ali, M.R. Efficient computational approaches for fractional-order Degasperis-Procesi and Camassa-Holm equations. Results Phys. 2023, 50, 106549. [Google Scholar] [CrossRef]
  17. Abu Hammad, M. Conformable Fractional Martingales and Some Convergence Theorems. Mathematics 2021, 10, 6. [Google Scholar] [CrossRef]
  18. Dahmani, Z.; Anber, A.; Gouari, Y.; Kaid, M.; Jebril, I. Extension of a Method for Solving Nonlinear Evolution Equations Via Conformable Fractional Approach. In Proceedings of the 2021 International Conference on Information Technology (ICIT 2021), Amman, Jordan, 14–15 July 2021; pp. 38–42. [Google Scholar]
  19. Batiha, I.M.; Oudetallah, J.; Ouannas, A.; Al-Nana, A.A.; Jebril, I.H. Tuning the fractional-order pid-controller for blood glucose level of diabetic patients. Int. J. Adv. Soft Comput. Its Appl. 2021, 13, 1–10. [Google Scholar]
  20. Deng, W.; Li, C. Existence and uniqueness of solutions for the fractional Fornberg-Whitham equation with initial and boundary conditions. Appl. Math. Lett. 2010, 23, 937–942. [Google Scholar] [CrossRef] [Green Version]
  21. Liu, F.; Anh, V. Well-posedness of the fractional Fornberg-Whitham equation with different types of boundary conditions. Comput. Math. Appl. 2011, 62, 1295–1303. [Google Scholar]
  22. Zhang, H.; Deng, W. A finite difference scheme for the fractional Fornberg-Whitham equation. J. Comput. Appl. Math. 2013, 239, 12–23. [Google Scholar]
  23. Liu, F.; Li, X.; Zhao, X. A finite volume method for the fractional Fornberg-Whitham equation. J. Comput. Phys. 2015, 295, 336–353. [Google Scholar]
  24. Li, C.; Deng, W.; Zhu, M. A spectral method for the fractional Fornberg-Whitham equation. Numer. Algorithms 2018, 79, 377–392. [Google Scholar]
  25. Hu, X.; Li, C.; Deng, W. Fractional Fornberg-Whitham equation for the dynamics of stock prices. J. Appl. Math. Comput. 2016, 50, 601–612. [Google Scholar]
  26. Wang, Y.; Zhang, C.; Song, W. Image denoising using the fractional Fornberg-Whitham equation. J. Comput. Appl. Math. 2015, 279, 152–161. [Google Scholar]
  27. Zhang, J.; Xie, J.; Shi, W.; Huo, Y.; Ren, Z.; He, D. Resonance and bifurcation of fractional quintic Mathieu-Duffing system. Chaos Interdiscip. J. Nonlinear Sci. 2023, 33, 23131. [Google Scholar] [CrossRef] [PubMed]
  28. Qi, M.; Cui, S.; Chang, X.; Xu, Y.; Meng, H.; Wang, Y.; Arif, M. Multi-region Nonuniform Brightness Correction Algorithm Based on L-Channel Gamma Transform. Secur. Commun. Netw. 2022, 2022, 2675950. [Google Scholar] [CrossRef]
  29. Zhu, H.; Xue, M.; Wang, Y.; Yuan, G.; Li, X. Fast Visual Tracking with Siamese Oriented Region Proposal Network. IEEE Signal Process. Lett. 2022, 29, 1437. [Google Scholar] [CrossRef]
  30. Guo, F.; Zhou, W.; Lu, Q.; Zhang, C. Path extension similarity link prediction method based on matrix algebra in directed networks. Comput. Commun. 2022, 187, 83–92. [Google Scholar] [CrossRef]
  31. Song, J.; Mingotti, A.; Zhang, J.; Peretto, L.; Wen, H. Accurate Damping Factor and Frequency Estimation for Damped Real-Valued Sinusoidal Signals. IEEE Trans. Instrum. Meas. 2022, 71, 6503504. [Google Scholar] [CrossRef]
  32. He, J.H. variational iteration method-a kind of nonlinear analytical technique: Some examples. Int. J. Non-Linear Mech. 2007, 34, 699–708. [Google Scholar] [CrossRef]
  33. He, J.H. variational iteration method for autonomous ordinary differential systems. Appl. Math. Comput. 2010, 217, 869–877. [Google Scholar] [CrossRef]
  34. Khader, M.M.; Hashim, I. Numerical methods for solving fractional differential equations: A comparative study. J. Comput. Appl. Math. 2016, 305, 195–210. [Google Scholar]
  35. Gao, G.H.; Li, X.Z.; He, J.H. Chaos in the fractional order Chen system and its control. Chaos Solitons Fractals 2004, 22, 549–554. [Google Scholar] [CrossRef]
  36. Hu, Y.; Sun, Z. Variational iteration transform method for solving the coupled Burgers’ equations with time-fractional derivatives. Appl. Math. Comput. 2017, 303, 132–141. [Google Scholar]
  37. Shah, N.A.; Alyousef, H.A.; El-Tantawy, S.A.; Chung, J.D. Analytical investigation of fractional-order Korteweg-De-Vries-type equations under Atangana-Baleanu-Caputo operator: Modeling nonlinear waves in a plasma and fluid. Symmetry 2022, 14, 739. [Google Scholar] [CrossRef]
  38. Xu, L.; Cao, X. The variational iteration transform method for solving the time-space fractional Fisher equation. Appl. Math. Comput. 2017, 305, 188–194. [Google Scholar]
  39. Wang, J.; Tian, J.; Zhang, X.; Yang, B.; Liu, S.; Yin, L.; Zheng, W. Control of Time Delay Force Feedback Teleoperation System with Finite Time Convergence. Front. Neurorobot. 2022, 16, 877069. [Google Scholar] [CrossRef]
  40. Jafari, H.; Seifi, S. Analytical solution of a nonlinear differential equation using the Variational Iteration Transform Method. J. Math. Anal. Appl. 2017, 446, 1261–1275. [Google Scholar]
  41. Adomian, G. A review of the decomposition method and some recent results for nonlinear equations. Math. Comput. Model. 1988, 13, 17–43. [Google Scholar] [CrossRef]
  42. Wazwaz, A.M. A First Course in Integral Equations; World Scientific: Singapore, 2002. [Google Scholar]
  43. Momani, S.; Odibat, Z. Analytical solution of a time-fractional Navier-Stokes equation by Adomian decomposition method. Appl. Math. Comput. 2007, 177, 488–494. [Google Scholar] [CrossRef]
  44. Abbasbandy, S.; Shirzadi, A. Application of the Adomian decomposition method for solving a system of nonlinear fractional differential equations. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 210–219. [Google Scholar]
  45. Eftekhari, G.; Alhuthali, M.S. Solving fractional partial differential equations using the Adomian decomposition method. J. Comput. Appl. Math. 2018, 339, 318–328. [Google Scholar]
  46. Cakir, M.; Arslan, D. The Adomian Decomposition Method and the Differential Transform Method for Numerical Solution of Multi-Pantograph Delay Differential Equations. Appl. Math. 2015, 6, 1332. [Google Scholar] [CrossRef] [Green Version]
  47. Bhrawy, A.H.; Alofi, A.S. Solving nonlinear differential equations by the modified Adomian decomposition method with application to wave equation. Results Phys. 2021, 26, 104708. [Google Scholar]
  48. Benattia, M.E.; Belghaba, K. Application of the Aboodh transform for solving fractional delay differential equations. Univers. J. Math. Appl. 2020, 3, 93–101. [Google Scholar] [CrossRef]
  49. Awuya, M.A.; Subasi, D. Aboodh transform iterative method for solving fractional partial differential equation with Mittag-Leffler Kernel. Symmetry 2021, 13, 2055. [Google Scholar] [CrossRef]
  50. Gupta, P.K.; Singh, M. Homotopy perturbation method for fractional Fornberg-Whitham equation. Comput. Math. Appl. 2011, 61, 250–254. [Google Scholar] [CrossRef]
  51. Abidi, F.; Omrani, K. Numerical solutions for the nonlinear Fornberg-Whitham equation by He’s methods. Int. J. Mod. Phys. B 2011, 25, 4721–4732. [Google Scholar] [CrossRef]
Figure 1. In this figure, both the approximate solution (18) or (24) at γ = = 1 using ADTM/VITM and the exact solution (25) are considered for Example 1.
Figure 1. In this figure, both the approximate solution (18) or (24) at γ = = 1 using ADTM/VITM and the exact solution (25) are considered for Example 1.
Symmetry 15 01353 g001
Figure 2. In this figure, both the approximate solution (18) or (24) using ADTM/VITM for Example 1 at different values of is considered: (a) = 0.2 , (b) = 0.5 (c) = 0.8 , and (d) two-dimensional plot for the cases (ac) at t = 0.1 .
Figure 2. In this figure, both the approximate solution (18) or (24) using ADTM/VITM for Example 1 at different values of is considered: (a) = 0.2 , (b) = 0.5 (c) = 0.8 , and (d) two-dimensional plot for the cases (ac) at t = 0.1 .
Symmetry 15 01353 g002
Figure 3. In this figure, both the approximate solution (32) or (38) at = 1 using ADTM/VITM and the exact solution (40) are considered for Example 2.
Figure 3. In this figure, both the approximate solution (32) or (38) at = 1 using ADTM/VITM and the exact solution (40) are considered for Example 2.
Symmetry 15 01353 g003
Figure 4. In this figure, both the approximate solution (32) or (38) using ADTM/VITM for Example 2 at different values of is considered: (a) = 0.2 , (b) = 0.5 (c) = 0.8 , and (d) two-dimensional plot for the cases (ac) at t = 0.1 .
Figure 4. In this figure, both the approximate solution (32) or (38) using ADTM/VITM for Example 2 at different values of is considered: (a) = 0.2 , (b) = 0.5 (c) = 0.8 , and (d) two-dimensional plot for the cases (ac) at t = 0.1 .
Symmetry 15 01353 g004
Table 1. Numerical values for both approximate and exact solutions for Example 1 and the absolute error.
Table 1. Numerical values for both approximate and exact solutions for Example 1 and the absolute error.
xADTM/VITMExact SolutionAbsolute Error
−1.00.5674140.5678390.000425389
−0.90.5965060.5969530.000447199
−0.80.6270890.6275590.000470128
−0.70.6592410.6597350.000494232
−0.60.6930410.693560.000519572
−0.50.7285740.729120.000546211
−0.40.7659280.7665030.000574215
−0.30.8051980.8058020.000603656
−0.20.8464820.8471160.000634606
−0.10.8898820.8905490.000667143
0.00.9355070.9362080.000701348
0.10.9834710.9842090.000737307
0.21.03391.034670.00077511
0.31.08691.087720.00081485
0.41.142631.143490.000856629
0.51.201211.202120.000900549
Table 2. Numerical values of both approximate and exact solutions for Example 2 and the absolute error.
Table 2. Numerical values of both approximate and exact solutions for Example 2 and the absolute error.
xADTM/VITMExact SolutionAbsolute Error
0.01.00211.000890.00121594
0.11.000430.9996160.000818105
0.21.000020.9995970.000420427
0.31.000851.000830.0000183526
0.41.002941.003330.000392684
0.51.006281.00710.000817276
0.61.010891.012150.00126006
0.71.016781.01850.00172572
0.81.023961.026180.00221902
0.91.032451.035190.0027448
1.01.042271.045570.003308
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

Noor, S.; Hammad, M.A.; Shah, R.; Alrowaily, A.W.; El-Tantawy, S.A. Numerical Investigation of Fractional-Order Fornberg–Whitham Equations in the Framework of Aboodh Transformation. Symmetry 2023, 15, 1353. https://doi.org/10.3390/sym15071353

AMA Style

Noor S, Hammad MA, Shah R, Alrowaily AW, El-Tantawy SA. Numerical Investigation of Fractional-Order Fornberg–Whitham Equations in the Framework of Aboodh Transformation. Symmetry. 2023; 15(7):1353. https://doi.org/10.3390/sym15071353

Chicago/Turabian Style

Noor, Saima, Ma’mon Abu Hammad, Rasool Shah, Albandari W. Alrowaily, and Samir A. El-Tantawy. 2023. "Numerical Investigation of Fractional-Order Fornberg–Whitham Equations in the Framework of Aboodh Transformation" Symmetry 15, no. 7: 1353. https://doi.org/10.3390/sym15071353

APA Style

Noor, S., Hammad, M. A., Shah, R., Alrowaily, A. W., & El-Tantawy, S. A. (2023). Numerical Investigation of Fractional-Order Fornberg–Whitham Equations in the Framework of Aboodh Transformation. Symmetry, 15(7), 1353. https://doi.org/10.3390/sym15071353

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