Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Surjective Identifications of Convex Noetherian Separations in Topological (C, R) Space
Next Article in Special Issue
The Existence, Uniqueness, and Stability Analysis of the Discrete Fractional Three-Point Boundary Value Problem for the Elastic Beam Equation
Previous Article in Journal
Exploring the High-Temperature Electrical Performance of Ca3−xLaxCo4O9 Thermoelectric Ceramics for Moderate and Low Substitution Levels
Previous Article in Special Issue
General Fractional Integrals and Derivatives of Arbitrary Order
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Comparative Study for Solving Fractional-Order Fornberg–Whitham Equation via ρ-Laplace Transform

1
Department of Mathematics and Computer Science, Faculty of Science and Technology, Rajamangala University of Technology Thanyaburi (RMUTT), Pathumthani 12110, Thailand
2
Department of Mathematics, College of Science, King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia
3
Department of Mathematics, Faculty of Science, Al-Azhar University, Assuit 71511, Egypt
4
School of Mathematics and Information Science, Henan Polytechnic University, Jiaozuo 454000, China
5
Department of Mathematics, Abdul Wali khan University, Mardan 23200, Pakistan
6
Department of Computer Engineering, Biruni University, Istanbul 34096, Turkey
7
Department of Mathematics, Science Faculty, Firat University, Elazig 23119, Turkey
8
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40402, Taiwan
*
Authors to whom correspondence should be addressed.
Symmetry 2021, 13(5), 784; https://doi.org/10.3390/sym13050784
Submission received: 29 March 2021 / Revised: 14 April 2021 / Accepted: 16 April 2021 / Published: 1 May 2021
(This article belongs to the Special Issue Applied Mathematics and Fractional Calculus)

Abstract

:
In this article, we also introduced two well-known computational techniques for solving the time-fractional Fornberg–Whitham equations. The methods suggested are the modified form of the variational iteration and Adomian decomposition techniques by ρ -Laplace. Furthermore, an illustrative scheme is introduced to verify the accuracy of the available methods. The graphical representation of the exact and derived results is presented to show the suggested approaches reliability. The comparative solution analysis via graphs also represented the higher reliability and accuracy of the current techniques.

1. Introduction

With engineering and science development, non-linear evolution models have been analyzed as the problems to define physical phenomena in plasma waves, fluid mechanics, chemical physics, solid-state physics, etc. For the last few years, therefore, a lot of interest has been paid to the result (both numerical and analytical) of these significant models [1,2,3,4]. Different methods are available in the literature for the approximate and exact results of these models. In current years, fractional calculus (FC) applied in many phenomena in applied sciences, fluid mechanics, physics and other biology can be described as very effective using mathematical tools of FC. The fractional derivatives have occurred in many applied sciences equations such as reaction and diffusion processes, system identification, velocity signal analysis, relaxation of damping behaviour fabrics and creeping of polymer composites [5,6,7,8].
The investigation of non-linear wave models and their application is significant in different areas of engineering. Travelling wave notions are between the most attractive results for non-linear fractional-order partial differential equations (NLFPDEs). NLFPDEs are usually identified as mechanical processes and complex physical. Therefore, it is important to get exact results for non-linear time-fractional partial differential equations [9,10,11,12]. Overall, travelling wave results are between the exciting forms of products for NFPDEs. On the other hand, other NLFPDEs, such as the Camassa–Holm or the Kortewegde–Vries equations, have been well-known to have some moving wave solutions. These are non-linear multi-directional dispersive waves in shallow water design problems [13,14,15,16].
The FWE study is of crucial significance in different areas of mathematical physics. The FWEs [15,16] is defined as
D δ μ D φ φ μ + D φ μ = μ D φ φ φ μ μ D φ μ + 3 D φ μ D φ φ μ .
The quantities performance of wave deformation, a non-linear dispersive wave model, is shown in the investigation. The FWE is presented as a mathematical model for limiting wave heights and wave breaks, allowing peakon results as a numerical model. In 1978, Fornberg and Whitham achieved a measured outcome of the form μ ( φ , η ) = C e φ 2 4 3 , where C is constant. The investigation of FWEs has been carried out by several analytical and numerical techniques, such as Adomian decomposition transform method [17], variational iteration technique [18], Lie Symmetry [19], new iterative method [20], differential transformation method [21], homotopy analysis transformation technique [22] and homotopy-perturbation technique [23].
Recently, Abdeljawad and Fahd [24] introduced the Laplace transformation of the fractional-order Caputo derivatives. We suggested a new iterative technique with ρ -Laplace transformation to investigate fractional-order ordinary and partial differential equations with fractional-order Caputo derivative. We apply this novel method for solving many fractional-order differential equations such as linear and non-linear diffusion equation, fractional-order Zakharov–Kuznetsov equation and Fokker–Planck equations. We analyzed the impact of δ and ρ in the process. The Variational iteration method (VIM) was first introduced by He [25,26] and was effectively implemented to the autonomous ordinary differential equation in [27], to non-linear polycrystalline solids [28], and other areas. Similarly, this technique is modified with ρ -Laplace transformation, so the modified method is called the ρ -Laplace variational iteration method. Many types of differential equations and partial differential equations have solved VITM. For example, this technique is analyzed for solving the time-fractional differential equation (FDEs) in [27]. In [28], this technique is applied to solve non-linear oscillator models. Compared to Adomian’s decomposition process, VITM solves the problem without the need to compute Adomian’s polynomials. This scheme provides a quick result to the equation, whereas the [29] mesh point techniques provide an analytical solution. This method can also be used to get a close approximation of the exact result. G. Adomian, an American mathematician, developed the Adomian decomposition technique. It focuses on finding series-like results and decomposing the non-linear operator into a sequence, with the terms presently computed using Adomian polynomials [30]. This method is modified with ρ -Laplace transform, so the modified approach is the ρ -Laplace decomposition method. This technique is used for the non-homogeneous FDEs [31,32,33,34,35,36].
This paper has implemented the ρ -Laplace variational iteration method and ρ -Laplace decomposition method to solve the time-fractional Fornberg–Whitham equations with the Caputo fractional derivative operator. The ρ -LDM and ρ -LVIM achieve the approximate results in the form of series results.

2. Basic Definitions

In this section, the fractional generalized derivative, the fractional generalized integral, the Mittag-Lefller function the ρ -Laplace transform have been discussed.
Definition 1.
The generalized fractional-order integral δ of a continuous function f : [ 0 , + ] R is expressed as [24]
( I δ , ρ f ) ( ζ ) = 1 Γ ( δ ) 0 ζ ζ ρ s ρ ρ δ 1 f ( s ) d s s 1 ρ ,
the gamma function denote by Γ, ρ > 0 , ζ > 0 and 0 < δ < 1 .
Definition 2.
The generalized fractional-order derivative of δ of a continuous function f : [ 0 , + ] R is given as [24].
( D δ , ρ f ) ( ζ ) = ( I 1 δ , ρ f ) ( ζ ) = 1 Γ ( 1 δ ) d d ζ 0 ζ ζ ρ s ρ ρ δ f ( s ) d s s 1 ρ .
where define the gamma function Γ, ρ > 0 , ζ > 0 and 0 < δ < 1 .
Definition 3.
The Caputo fractional-order derivative δ of a continuous function f : [ 0 , + ] R is expressed as [24]
( D δ , ρ f ) ( ζ ) = 1 Γ ( 1 δ ) 0 ζ ζ ρ s ρ ρ δ β n f ( s ) d s s 1 ρ .
where n = 1 , ρ > 0 , ζ > 0 , β = ζ 1 ρ d d ζ and 0 < δ < 1 .
Definition 4.
The ρ-Laplace transformation of a continuous function f : [ 0 , + ] R is given as [24]
L ρ { f ( ζ ) } = 0 e s ζ ρ ρ f ( ζ ) d ζ ζ 1 ρ .
The Caputo generalized fractional-order ρ-Laplace transform derivative of a continuous function f is defined by [24].
L ρ { D δ , ρ f ( ζ ) } = s δ L ρ { f ( ζ ) } k = 0 n 1 s δ k 1 ( I δ , ρ β n f ) ( 0 ) n = 1 .

3. The General Methodology of ρ -LDM

The ρ -LDM is a combination of the Laplace decomposition method and the ρ -Laplace transformation. In this section, we solve the ρ -LDM solution of fractional partial differential equation. The main steps of this method are described as follows:
D δ , ρ ω ( φ , ) + L ¯ ( φ , ) + N ( φ , ) H ( φ , ) = 0 , 0 < δ 1 ,
where L ¯ and N are linear and nonlinear functions, H is the sources function.
The initial condition is
ω ( φ , 0 ) = f ( φ ) ,
Apply ρ -Laplace transform to Equation (2),
L ρ [ D δ , ρ ω ( φ , ) ] + L ρ [ L ¯ ( φ , ) + N ( φ , ) H ( φ , ) ] = 0 .
Applying the ρ -Laplace transformation differentiation property, we get
L ρ [ ω ( φ , ) ] = 1 s ω ( φ , 0 ) + 1 s δ L ρ [ H ( φ , ) ] 1 s δ L ρ { L ¯ ( φ , ) + N ( φ , ) } ] .
ρ -LDM solution of infinite series ω ( φ , ) ,
ω ( φ , ) = j = 0 ω m ( φ , ) .
The N is the nonlinear term defined as
N ( φ , ) = j = 0 A m .
So the with the help of Adomian polynomial we can define the nonlinear terms
A m = 1 j ! m λ m N k = 0 λ k ω k λ = 0 .
Putting Equations (6) and (7) into (5), we get
L ρ j = 0 ω m ( φ , ) = 1 s ω ( φ , 0 ) + 1 s δ S { H ( φ , ) } 1 s δ L ρ L ¯ ( j = 0 ω m ) + j = 0 A m .
Using the inverse ρ -Laplace transform with Equation (9),
j = 0 ω m ( φ , ) = L ρ 1 1 s ω ( φ , 0 ) + 1 s δ L ρ { H ( φ , ) } 1 s δ L ρ L ¯ j = 0 ω m + j = 0 A m .
we define the next terms,
ω 0 ( φ , ) = L ρ 1 1 s ω ( φ , 0 ) + 1 s δ L ρ { H ( φ , ) } ,
ω 1 ( φ , ) = L ρ 1 1 s δ L ρ { L ¯ 1 ( ω 0 ) + A 0 } .
For m 1 , is expressed as
ω j + 1 ( φ , ) = L ρ 1 1 s δ L ρ { L ¯ ( ω m ) + A m } .

4. Convergence Analysis

Theorem 1.
[37] (Uniqueness theorem) Equation has a unique solution whenever 0 < ε < 1 where ε = ( h 1 + h 2 + h 3 ) δ + 1 ( δ 1 ) ! .
Theorem 2.
[37] (Convergence Theorem) The series solution (11) and (12) of the problem (3) using ρ -LTADM and ρ -LTVIM converges if 0 < ε < 1 .
Proof. 
Let S be the mth partial sum, i.e., S = j = 0 m ω ( φ , ) . We shall prove that S is a Cauchy sequence in Banach space E. By using a new formulation of Adomian polynomials we get [37]
R ( S ) = A ^ + j = 0 m 1 A j ^ ( S ) = A ^ + n = 0 m 1 A n ^
| | S S m 1 | | = max I | S S m 1 | = max I j = n + 1 m ω j ^ ( φ , ) , j = 0 , 1 , 2 max I L ρ 1 { 1 s δ L ρ { [ j = n + 1 m k [ ω j 1 ( φ , ) ] ] } + L ρ 1 { 1 s δ L ρ { [ j = n + 1 m M [ ω j 1 ( φ , ) ] ] } + L ρ 1 { 1 s δ L ρ { [ j = n + 1 m [ A j 1 ( φ , ) ] ] } , max I L ρ 1 { 1 s δ L ρ { [ j = n + 1 m k [ ω j ( φ , ) ] ] } + L ρ 1 { 1 s δ L ρ { [ j = n + 1 m M [ ω j ( φ , ) ] ] } + L ρ 1 { 1 s δ L ρ { [ j = n + 1 m [ A j ( φ , ) ] ] } , max I L ρ 1 { 1 s δ L ρ { [ j = n + 1 m k [ S m 1 1 S m 2 1 ] ] } + L ρ 1 { 1 s δ L ρ { [ j = n + 1 m M [ S m 1 1 S m 2 1 ] ] } + L ρ 1 { 1 s δ L ρ { [ j = n + 1 m [ S m 1 1 S m 2 1 ] ] } , max I L ρ 1 { 1 s δ L ρ { [ k [ S m 1 1 S m 2 1 ] ] } + L ρ 1 { 1 s δ L ρ { [ M [ S m 1 1 S m 2 1 ] ] } + L ρ 1 { 1 s δ L ρ { [ [ S m 1 1 S m 2 1 ] ] } , k 1 max I L ρ 1 { 1 s δ L ρ { [ S m 1 1 S m 2 1 ] } , +   k 2 max I L ρ 1 { 1 s δ L ρ { [ S m 1 1 S m 2 1 ] } , +   k 3 max I L ρ 1 1 s δ L ρ { [ S m 1 1 S m 2 1 ] } , = ( k 1 + k 2 + k 3 ) δ 1 ( δ 1 ) ! S m 1 1 S m 2 1 .
Letting m 1 = m 2 + 1 , we get
S m 2 + 1 S m 2   ε S m 2 S m 2 1   ε 2 S m 2 1 S m 2 2   ε m 2 S 1 S 0 ,
where ε = ( k 1 + k 2 + k 3 ) δ 1 ( δ 1 ) ! similarly, we have from the triangle inequality we get
S m 1 1 S m 2 1     S m 1 + 1 S m 2   +   S m 1 + 2 S m 2 + 1   + +   S m 1 S m 1 1 ,   ε m 2 + ε m 2 + 1 + + ε m 1 1     S 1 + S 0 , ε m 2 ( 1 ε m 1 m 2 ε ) ω 1 .
Since 0 < ε < 1 we have 1 ε m 1 m 2 < 1
S m 1 + S m 2   ε m 2 1 ε max I ω .
However | ω | < so, as m 2 then | | S m 1 S m 2 | | 0 , hence S m 1 is a Cauchy sequence, the series m 1 = 0 ω m 1 converges and the proof is complete.  □
Theorem 3.
[37] (Error estimate) The maximum absolute error of the series solution can be given the following formula
max I ω ( φ , ) = 1 ω ( φ , ) ε m 2 1 ε max I ω 1 .

5. The General Methodology of ρ -Laplace Variational Iteration Method

In this section we show the general methodology of the ρ -Laplace variational iteration method solution for fractional partial differential equations.
D δ , ρ ω ( φ , ) + L ¯ ( φ , ) + N ( φ , ) H ( φ , ) = 0 , 0 < δ 1 ,
with the initial condition
ω ( φ , 0 ) = f ( φ ) ,
The using ρ -Laplace transformation to Equation (12),
L ρ [ D δ , ρ ω ( φ , ) ] + L ρ [ L ¯ ( φ , ) + N ( φ , ) H ( φ , ) ] = 0 .
Applying the differentiation property of ρ -Laplace transform, we get
s δ L ρ [ ω ( φ , ) ] s δ 1 ω ( φ , 0 ) = L ρ L ¯ ( φ , ) + N ( φ , ) H ( φ , ) .
The Lagrange multiplier is used in the iterative method
L ρ [ ω j + 1 ( φ , ) ] = L ρ [ ω j ( φ , ) ] + λ ( s ) s δ L ρ [ ω j ( φ , ) ] s δ 1 ω j ( φ , 0 ) L ρ { L ¯ ( φ , ) + N ( φ , ) } L ρ [ H ( φ , ) ] .
The Lagrange multiplier is
λ ( s ) = 1 s δ ,
using inverse ρ -Laplace transform L 1 , Equation (16), we get
ω j + 1 ( φ , ) = ω j ( φ , ) L ρ 1 1 s δ L ρ { L ¯ ( φ , ) + N ( φ , ) } L ρ [ H ( φ , ) ] ,
the initial value can be defined as
ω 0 ( φ , ) = L ρ 1 1 s δ s δ 1 ω ( φ , 0 ) .

6. Implementation of Techniques

We now proceed to derive an approximate solution to the time-fractional nonlinear FW equations using suggested techniques with generalized Caputo fractional derivative.

6.1. Problem

Consider the time-fractional nonlinear FWE is given as
D δ , ρ ω D φ φ ω + D φ ω = ω D φ φ φ ω ω D φ ω + 3 D φ ω D φ φ ω , 0 < δ 1 ,
the initial condition is
ω ( φ , 0 ) = e φ 2 .
Taking ρ -Laplace transform of (20),
s δ L ρ ω ( φ , ) s δ 1 ω ( φ , 0 ) = L ρ D φ φ ω D φ ω + ω D φ φ φ ω ω D φ ω + 3 D φ ω D φ φ ω .
Applying inverse ρ -Laplace transform
ω ( φ , ) = L ρ 1 ω ( φ , 0 ) s 1 s δ L ρ D φ φ ω D φ ω + ω D φ φ φ ω ω D φ ω + 3 D φ ω D φ φ ω .
Using ADM procedure, we get
ω 0 ( φ , ) = L ρ 1 ω ( φ , 0 ) s = L ρ 1 e φ 2 s ,
ω 0 ( φ , t ) = e φ 2 ,
= 0 ω + 1 ( φ , ) = L ρ 1 1 s δ L ρ = 0 ( D φ φ ω ) = 0 ( D φ ω ) + = 0 A = 0 B + 3 = 0 C , = 0 , 1 , 2 ,
A 0 ( ω D φ φ φ ω ) = ω 0 D φ φ φ ω 0 , A 1 ( ω D φ φ φ ω ) = ω 0 D φ φ φ ω 1 + ω 1 D φ φ φ ω 0 , A 2 ( ω D φ φ φ ω ) = ω 1 D φ φ φ ω 2 + ω 1 D φ φ φ ω 1 + ω 2 D φ φ φ ω 0 ,
B 0 ( ω D φ ω ) = ω 0 D φ ω 0 , B 1 ( ω D φ ω ) = ω 0 D φ ω 1 + ω 1 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 ( φ , ) = L ρ 1 1 s δ L ρ D φ φ ω 0 D φ ω 0 + A 0 B 0 + 3 C 0 , ω 1 ( φ , t ) = 1 2 L ρ 1 e φ 2 s δ + 1 = 1 2 e φ 2 ρ ρ δ Γ ( δ + 1 ) .
for = 2
ω 2 ( φ , ) = L ρ 1 1 s δ L ρ D φ φ ω 1 D φ ω 1 + A 1 B 1 + 3 C 1 , ω 2 ( φ , ) = 1 8 e φ 2 ρ ρ 2 δ 1 Γ ( 2 δ ) + 1 4 e φ 2 ρ ρ 2 δ Γ ( 2 δ + 1 ) ,
for = 3
ω 3 ( φ , ) = L ρ 1 1 s δ L ρ 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 ρ ρ 3 δ Γ ( 3 δ + 1 ) ,
The ρ -LDM result of Example 1 is
ω ( φ , ) = ω 0 ( φ , ) + ω 1 ( φ , ) + ω 2 ( φ , ) + ω 3 ( φ , ) + ω 4 ( φ , ) + ,
ω ( φ , ) = e φ 2 1 2 e φ 2 ρ ρ δ Γ ( δ + 1 ) 1 8 e φ 2 ρ ρ 2 δ 1 Γ ( 2 δ ) + 1 4 e φ 2 ρ ρ 2 δ Γ ( 2 δ + 1 ) 1 32 e φ 2 ρ ρ 3 δ 2 Γ ( 3 δ 1 ) + 1 8 e φ 2 3 δ 1 Γ ( 3 δ ) 1 8 e φ 2 ρ ρ 3 δ Γ ( 3 δ + 1 ) .
The simplify we can write Equation (26), we get
ω ( φ , ) = e φ 2 1 ρ ρ δ 2 Γ ( δ + 1 ) 1 8 ρ ρ 2 δ 1 Γ ( 2 δ ) + 1 4 ρ ρ 2 δ Γ ( 2 δ + 1 ) 1 32 ρ ρ 3 δ 2 Γ ( 3 δ 1 ) + 1 8 ρ ρ 3 δ 1 Γ ( 3 δ ) 1 8 ρ ρ 3 δ Γ ( 3 δ + 1 ) + .
The analytical result by ρ -LVIM.
The iteration method apply for Equation (20), we get
ω + 1 ( φ , ) = ω ( φ , ) L ρ 1 1 s δ L ρ s δ D ω D φ φ ω + D φ ω ω D φ φ φ ω + ω D φ ω 3 D φ ω D φ φ ω ,
where
ω 0 ( φ , ) = e φ 2 .
For = 0 , 1 , 2 ,
ω 1 ( φ , ) = ω 0 ( φ , ) L ρ 1 1 s δ L ρ s δ D ω 0 D φ φ ω 0 + D φ ω 0 ω 0 D φ φ φ ω 0 + ω 0 D φ ω 0 3 D φ ω 0 D φ φ ω 0 , ω 1 ( φ , ) = 1 2 e φ 2 ρ ρ δ Γ ( δ + 1 ) ,
ω 2 ( φ , ) = ω 1 ( φ , ) L ρ 1 1 s δ L ρ s δ D ω 1 D φ φ ω 1 + D φ ω 1 ω 1 D φ φ φ ω 1 + ω 1 D φ ω 1 3 D φ ω 1 D φ φ ω 1 , ω 2 ( φ , ) = 1 8 e φ 2 ρ ρ 2 δ 1 Γ ( 2 δ ) + 1 4 e φ 2 ρ ρ 2 δ Γ ( 2 δ + 1 ) ,
ω 3 ( φ , ) = ω 2 ( φ , ) L ρ 1 1 s δ L ρ s δ D ω 2 D φ φ ω 2 + D φ ω 2 ω 2 D φ φ φ ω 2 + ω 2 D φ ω 2 3 D φ ω 2 D φ φ ω 2 , ω 3 ( φ , ) = 1 32 e φ 2 ρ ρ 3 δ 2 Γ ( 3 δ 1 ) + 1 8 e φ 2 ρ ρ 3 δ 1 Γ ( 3 δ ) 1 8 e φ 2 ρ ρ 3 δ Γ ( 3 δ + 1 ) ,
ω ( φ , ) = m = 0 ω m ( φ ) = e φ 2 1 2 e φ 2 ρ ρ δ Γ ( δ + 1 ) 1 8 e φ 2 ρ ρ 2 δ 1 Γ ( 2 δ ) + 1 4 e φ 2 ρ ρ 2 δ Γ ( 2 δ + 1 ) 1 32 e φ 2 ρ ρ 3 δ 2 Γ ( 3 δ 1 ) + 1 8 e φ 2 3 δ 1 Γ ( 3 δ ) 1 8 e φ 2 ρ ρ 3 δ Γ ( 3 δ + 1 ) .
The exact result of Equation (20) at δ = 1 ,
ω ( φ , ) = e φ 2 2 3 .
Figure 1 shows the ρ -LDM and ρ -LVIM solution of the fractional Fornberg–Whitham defined by generalized fractional-order Caputo derivative in the space coordinate and time 0 < 0.5 , ρ = 1 and δ = 1 . Figure 2, the 3D graph shows approximate and exact solutions graph at δ = 1 and ρ = 0.9 ; the figure shows that different fractional-order at δ . Similarly, in Figure 3, the 2D graph of exact and approximate solutions plot at δ = 1 and ρ = 0.9 the figure shows that different fractional-order at δ .

6.2. Problem

Consider the time-fractional non-linear FWE given as
D δ , ρ ω D φ φ ω + D φ ω = ω D φ φ φ ω ω D φ ω + 3 D φ ω D φ φ ω , > 0 , 0 < δ 1 ,
with the initial condition
ω ( φ , 0 ) = cosh 2 φ 4 .
Taking ρ -Laplace transform of (35),
s δ L ρ ω ( φ , ) s δ 1 ω ( φ , 0 ) = L ρ D φ φ ω D φ ω + ω D φ φ φ ω ω D φ ω + 3 D φ ω D φ φ ω .
Applying inverse ρ -Laplace transform
ω ( φ , ) = L ρ 1 ω ( φ , 0 ) s 1 s δ L ρ D φ φ ω D φ ω + ω D φ φ φ ω ω D φ ω + 3 D φ ω D φ φ ω .
Using ADM procedure, we get
ω 0 ( φ , ) = L ρ 1 ω ( φ , 0 ) s = L ρ 1 cosh 2 φ 4 s ,
ω 0 ( φ , ) = cosh 2 φ 4 ,
= 0 ω + 1 ( φ , ) = L ρ 1 1 s δ L ρ = 0 ( D φ φ ω ) = 0 ( D φ ω ) + = 0 A = 0 B + 3 = 0 C , = 0 , 1 , 2 ,
for = 0
ω 1 ( φ , ) = L ρ 1 1 s δ L ρ D φ φ ω 0 D φ ω 0 + A 0 B 0 + 3 C 0 , ω 1 ( φ , ) = 11 32 L ρ 1 sinh x 2 s δ + 1 = 11 32 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) .
for = 1
ω 2 ( φ , ) = L ρ 1 1 s δ L ρ D φ φ ω 1 D φ ω 1 + A 1 B 1 + 3 C 1 , ω 2 ( φ , ) = 11 28 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) + 121 1024 cosh φ 4 ρ ρ 2 δ Γ ( 2 δ + 1 ) ,
for = 2
ω 3 ( φ , ) = L ρ 1 1 s δ L ρ D φ φ ω 2 D φ ω 2 + A 2 B 2 + 3 C 2 , ω 3 ( φ , ) = 11 512 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) + 121 2048 cosh φ 4 ρ ρ 2 δ Γ ( 2 δ + 1 ) 1331 49152 sinh φ 4 ρ ρ 3 δ Γ ( 3 δ + 1 ) ,
The ρ -LDM result for problem 2 is
ω ( φ , ) = ω 0 ( φ , ) + ω 1 ( φ , ) + ω 2 ( φ , ) + ω 3 ( φ , ) + ω 4 ( φ , ) + ,
ω ( φ , ) = cosh 2 φ 4 11 32 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) 11 28 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) + 121 1024 cosh φ 4 ρ ρ 2 δ Γ ( 2 δ + 1 ) 11 512 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) + 121 2048 cosh φ 4 ρ ρ 2 δ Γ ( 2 δ + 1 ) 1331 49152 sinh φ 4 ρ ρ 3 δ Γ ( 3 δ + 1 ) .
The analytical solution by ρ -LVIM.
The iteration method is apply by Equation (35), we get
ω + 1 ( φ , ) = ω ( φ , ) L ρ 1 1 s δ L ρ s δ D ω D φ φ ω + D φ ω ω D φ φ φ ω + ω D φ ω 3 D φ ω D φ φ ω ,
where
ω 0 ( φ , t ) = cosh 2 φ 4 .
For = 0 , 1 , 2 ,
ω 1 ( φ , ) = ω 0 ( φ , ) L ρ 1 1 s δ L ρ s δ D ω 0 D φ φ ω 0 + D φ ω 0 ω 0 D φ φ φ ω 0 + ω 0 D φ ω 0 3 D φ ω 0 D φ φ ω 0 , ω 1 ( φ , ) = cosh 2 φ 4 11 32 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) ,
ω 2 ( φ , ) = ω 1 ( φ , ) L ρ 1 1 s δ L ρ s δ D ω 1 D φ φ ω 1 + D φ ω 1 ω 1 D φ φ φ ω 1 + ω 1 D φ ω 1 3 D φ ω 1 D φ φ ω 1 , ω 2 ( φ , ) = cosh 2 φ 4 11 32 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) 11 28 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) + 121 1024 cosh φ 4 ρ ρ 2 δ Γ ( 2 δ + 1 ) ,
ω 3 ( φ , ) = ω 2 ( φ , ) L ρ 1 1 s δ L ρ s δ D ω 2 D φ φ ω 2 + D φ ω 2 ω 2 D φ φ φ ω 2 + ω 2 D φ ω 2 3 D φ ω 2 D φ φ ω 2 , ω 3 ( φ , ) = cosh 2 φ 4 11 32 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) 11 28 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) + 121 1024 cosh φ 4 ρ ρ 2 δ Γ ( 2 δ + 1 ) , 11 512 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) + 121 2048 cosh φ 4 ρ ρ 2 δ Γ ( 2 δ + 1 ) 1331 49152 sinh φ 4 ρ ρ 3 δ Γ ( 3 δ + 1 ) ,
ω ( φ , ) = m = 0 ω m ( φ ) = cosh 2 φ 4 11 32 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) 11 28 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) + 121 1024 cosh φ 4 ρ ρ 2 δ Γ ( 2 δ + 1 ) , 11 512 sinh φ 4 ρ ρ δ Γ ( δ + 1 ) + 121 2048 cosh φ 4 ρ ρ 2 δ Γ ( 2 δ + 1 ) 1331 49152 sinh φ 4 ρ ρ 3 δ Γ ( 3 δ + 1 ) .
The exact result of Equation (35) at δ = 1 ,
ω ( φ , ) = cosh 2 φ 4 11 24 .
Figure 4 shows the ρ -LDM and ρ -LVIM solution of the fractional Fornberg–Whitham defined by generalized Caputo fractional-order derivative in the space coordinate and time 0 < 0.5 , ρ = 1 and δ = 1 . Figure 5, the 3D graph shows exact and approximate solutions plot at δ = 1 and ρ = 0.9 ; the figure shows that different fractional-order at δ . Similarly, in Figure 6, the 2D graph of exact and approximate solutions plot at δ = 1 and ρ = 0.9 the figure shows that different fractional-order at δ .

7. Conclusions

In this article, different semi-analytical techniques are implemented to solve time-fractional Fornberg–Whitham equation. The approximate solution of the equations is evaluated to confirm the validity and reliability of the proposed methods. Graphs of the solutions are plotted to display the closed relation between the obtained and exact results. In addition, the suggested techniques provide easily computable components for the series-form tests. It is investigated that the results achieved in the series form have a higher convergence rate towards the exact results. The proposed methods have a small number of calculations to achieve the approximate solution. In conclusion, it is found that the proposed technique is a sophisticated method for solving other NLFPDEs. In the future, the analytical result of non-linear fractional-order boundary values problems achieved using this technique is in the form of uniform convergence series.

Author Contributions

Conceptualization, R.S. and M.I.; methodology, R.S.; software, P.S. and A.M.Z.; validation, R.S.; formal analysis, M.I. and P.S.; data curation, R.S.; writing—original draft preparation R.S.; writing—review and editing, A.M.Z.; supervision, M.I.; project administration, A.M.Z.; funding acquisition, S.-W.Y. All authors have read and agreed to the published version of the manuscript.

Funding

National Natural Science Foundation of China (No. 71601072), Key Scientific Research Project of Higher Education Institutions in Henan Province of China (No. 20B110006) and the Fundamental Research Funds for the Universities of Henan Province (No. NSFRF210314).

Data Availability Statement

Not applicable.

Acknowledgments

One of the Co-authors (A. M. Zidan) extend his appreciation to the Deanship of Scientific Research at King Khalid University, Abha 61413, Saudi Arabia for funding this work through research groups program under grant number R.G.P.1/30/42.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qureshi, S.; Yusuf, A.; Shaikh, A.A.; Inc, M.; Baleanu, D. Fractional modeling of blood ethanol concentration system with real data application. Chaos Interdiscip. J. Nonlinear Sci. 2019, 29, 013143. [Google Scholar] [CrossRef]
  2. Yusuf, A.; Qureshi, S.; Inc, M.; Aliyu, A.I.; Baleanu, D.; Shaikh, A.A. Two-strain epidemic model involving fractional derivative with Mittag-Leffler kernel. Chaos Interdiscip. J. Nonlinear Sci. 2018, 28, 123121. [Google Scholar] [CrossRef] [PubMed]
  3. Kilic, B.; Inc, M. The first integral method for the time fractional Kaup-Boussinesq system with time dependent coefficient. Appl. Math. Comput. 2015, 254, 70–74. [Google Scholar] [CrossRef]
  4. Inc, M.; Akgul, A.; Kilicman, A. Numerical solutions of the second-order one-dimensional telegraph equation based on reproducing kernel Hilbert space method. Abstr. Appl. Anal. 2013, 2013, 768963. [Google Scholar] [CrossRef] [Green Version]
  5. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000; Volume 35, pp. 87–130. [Google Scholar]
  6. Cesarano, C.; Pierpaolo, N.; Paolo, E.R. Pseudo-Lucas Functions of Fractional Degree and Applications. Axioms 2021, 10, 51. [Google Scholar] [CrossRef]
  7. Sabatier, J.A.T.M.J.; Agrawal, O.P.; Machado, J.T. Advances in Fractional Calculus; Springer: Dordrecht, The Netherlands, 2007; Volume 4. [Google Scholar]
  8. Al-luhaibi, M.S. An analytical treatment to fractional Fornberg-Whitham equation. Math. Sci. 2017, 11, 1–6. [Google Scholar] [CrossRef] [Green Version]
  9. Agarwal, P.; Ahs, S.; Akbare, M.; Nawaz, R.; Cesarano, C. A Reliable Algorithm for solution of Higher Dimensional Nonlinear (1+1) and (2+1) Dimensional Volterra-Fredholm Integral Equations. Dolomites Res. Notes Approx. 2021, 14, 18–25. [Google Scholar]
  10. Zakarya, M.; Altanji, M.; AlNemer, G.; Abd El-Hamid, H.A.; Cesarano, C.; Rezk, H.M. Fractional Reverse Coposn’s Inequalities via Conformable Calculus on Time Scales. Symmetry 2021, 13, 542. [Google Scholar] [CrossRef]
  11. Ahmad, I.; Ahmad, H.; Thounthong, P.; Chu, Y.-M.; Cesarano, C. Solution of Multi-Term Time-Fractional PDE Models Arising in Mathematical Biology and Physics by Local Meshless Method. Symmetry 2020, 12, 1195. [Google Scholar] [CrossRef]
  12. Bazighifan, O.; Cesarano, C. A Philos-type oscillation criteria for fourth-order neutral differential equations. Symmetry 2020, 12, 379. [Google Scholar] [CrossRef] [Green Version]
  13. Camacho, J.C.; Rosa, M.; Garias, M.L.; Bruzon, M.S. Classical symmetries, travelling wave solutions and conservation laws of a generalized Fornberg-Whitham equation. J. Comput. Appl. Math. 2017, 318, 149–155. [Google Scholar] [CrossRef]
  14. Bruzon, M.S.; Marquez, A.P.; Garrido, T.M.; Recio, E.; de la Rosa, R. Conservation laws for a generalized seventh order KdV equation. J. Comput. Appl. Math. 2019, 354, 682–688. [Google Scholar] [CrossRef]
  15. Whitham, G.B. Variational methods and applications to water waves. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1967, 299, 6–25. [Google Scholar]
  16. Fornberg, B.; Whitham, G.B. A numerical and theoretical study of certain nonlinear wave phenomena. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Sci. 1978, 289, 373–404. [Google Scholar]
  17. Kumar, D.; Singh, J.; Baleanu, D. A new analysis of the Fornberg-Whitham equation pertaining to a fractional derivative with Mittag-Leffler-type kernel. Eur. Phys. J. Plus 2018, 133, 1–10. [Google Scholar] [CrossRef]
  18. Lu, J. An analytical approach to the Fornberg-Whitham type equations by using the variational iteration method. Comput. Math. Appl. 2011, 61, 2010–2013. [Google Scholar] [CrossRef] [Green Version]
  19. Hashemi, M.S.; Haji-Badali, A.; Vafadar, P. Group invariant solutions and conservation laws of the Fornberg-Whitham equation. Z. Naturforschung A 2014, 69, 489–496. [Google Scholar] [CrossRef]
  20. Ramadan, M.A.; Al-luhaibi, M.S. New iterative method for solving the fornberg-whitham equation and comparison with homotopy perturbation transform method. J. Adv. Math. Comput. Sci. 2014, 4, 1213–1227. [Google Scholar] [CrossRef]
  21. Merdan, M.; Gokdogan, A.; Yildirim, A.; Mohyud-Din, S.T. Numerical simulation of fractional Fornberg-Whitham equation by differential transformation method. Abstr. Appl. Anal. 2012, 2012, 965367. [Google Scholar] [CrossRef]
  22. Wang, K.; Liu, S. Application of new iterative transform method and modified fractional homotopy analysis transform method for fractional Fornberg-Whitham equation. J. Nonlinear Sci. Appl. 2016, 9, 2419–2433. [Google Scholar] [CrossRef] [Green Version]
  23. 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]
  24. Jarad, F.; Abdeljawad, T. A modifi ed Laplace transform for certain generalized fractional operators. Results Nonlinear Anal. 2018, 1, 88–98. [Google Scholar]
  25. He, J.H. Approximate solution of nonlinear differential equations with convolution product nonlinearities. Comput. Methods Appl. Mech. Eng. 1998, 167, 69–73. [Google Scholar] [CrossRef]
  26. He, J.H. Variational iteration method for autonomous ordinary differential systems. Appl. Math. Comput. 2000, 114, 115–123. [Google Scholar] [CrossRef]
  27. Wu, G.C.; Baleanu, D. Variational iteration method for fractional calculus-a universal approach by Laplace transform. Adv. Differ. Equ. 2013, 2013, 18. [Google Scholar] [CrossRef] [Green Version]
  28. Anjum, N.; He, J.H. Laplace transform: Making the variational iteration method easier. Appl. Math. Lett. 2019, 92, 134–138. [Google Scholar] [CrossRef]
  29. Dehghan, M. Finite difference procedures for solving a problem arising in modeling and design of certain optoelectronic devices. Math. Comput. Simul. 2006, 71, 16–30. [Google Scholar] [CrossRef]
  30. Adomian, G. Solving Frontier Problems of Physics: The Decomposition Method; Kluwer Academic Publishers: Boston, MA, USA, 1994. [Google Scholar]
  31. Khalouta, A.; Kadem, A. A New Method to Solve Fractional Differential Equations: Inverse Fractional Laplace Transform Method. Appl. Appl. Math. 2019, 14, 926–941. [Google Scholar]
  32. Bokhari, A.; Baleanu, D.; Belgacem, R. Application of Laplace transform to Atangana-Baleanu derivatives. J. Math. Comput. Sci. 2019, 20, 101–107. [Google Scholar] [CrossRef] [Green Version]
  33. Belgacem, R.; Baleanu, D.; Bokhari, A. Laplace Transform and Applications to Caputo-Fractional Differential Equations. Int. J. Anal. Appl. 2019, 17, 917–927. [Google Scholar]
  34. Machado, J.; Baleanu, D.; Chen, W.; Sabatier, J. New trends in fractional dynamics. J. Vib. Control 2014, 20, 963. [Google Scholar] [CrossRef] [Green Version]
  35. Baleanu, D.; Guvenc, Z.; Machado, J. New Trends in Nanotechnology and Fractional Calculus Applications; Springer: Dordrecht, The Netherlands, 2010. [Google Scholar]
  36. Maitama, S.; Zhao, W. New integral transform: Laplace transform a generalization of Sumudu and Laplace transform for solving differential equations. arXiv 2019, arXiv:1904.11370. [Google Scholar]
  37. El-Kalla, I.L. Convergence of the Adomian method applied to a class of nonlinear integral equations. Appl. Math. Lett. 2008, 21, 372–376. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The graph of Exact and analytical solutions of δ = 1 and ρ = 1 of problem 1.
Figure 1. The graph of Exact and analytical solutions of δ = 1 and ρ = 1 of problem 1.
Symmetry 13 00784 g001
Figure 2. The first 3D graph of Exact and analytical solutions graph at δ = 1 and ρ = 0.9 and second plot of the approximate different fractional-order of δ = 1 of problem 1.
Figure 2. The first 3D graph of Exact and analytical solutions graph at δ = 1 and ρ = 0.9 and second plot of the approximate different fractional-order of δ = 1 of problem 1.
Symmetry 13 00784 g002
Figure 3. The first 2D graph of Exact and analytical solutions graph at δ = 1 and ρ = 0.9 and second plot of the approximate different fractional-order of δ = 1 of problem 1.
Figure 3. The first 2D graph of Exact and analytical solutions graph at δ = 1 and ρ = 0.9 and second plot of the approximate different fractional-order of δ = 1 of problem 1.
Symmetry 13 00784 g003
Figure 4. The graph of Exact and approximate solutions of δ = 1 and ρ = 1 of Example 2.
Figure 4. The graph of Exact and approximate solutions of δ = 1 and ρ = 1 of Example 2.
Symmetry 13 00784 g004
Figure 5. The first 3D graph of Exact and approximate solutions plot at δ = 1 and ρ = 0.9 and second plot of the approximate different fractional-order of δ = 1 of Example 2.
Figure 5. The first 3D graph of Exact and approximate solutions plot at δ = 1 and ρ = 0.9 and second plot of the approximate different fractional-order of δ = 1 of Example 2.
Symmetry 13 00784 g005
Figure 6. The first 2D graph of Exact and approximate solutions plot at δ = 1 and ρ = 0.9 and second plot of the approximate different fractional-order of δ = 1 of Example 2.
Figure 6. The first 2D graph of Exact and approximate solutions plot at δ = 1 and ρ = 0.9 and second plot of the approximate different fractional-order of δ = 1 of Example 2.
Symmetry 13 00784 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sunthrayuth, P.; Zidan, A.M.; Yao, S.-W.; Shah, R.; Inc, M. The Comparative Study for Solving Fractional-Order Fornberg–Whitham Equation via ρ-Laplace Transform. Symmetry 2021, 13, 784. https://doi.org/10.3390/sym13050784

AMA Style

Sunthrayuth P, Zidan AM, Yao S-W, Shah R, Inc M. The Comparative Study for Solving Fractional-Order Fornberg–Whitham Equation via ρ-Laplace Transform. Symmetry. 2021; 13(5):784. https://doi.org/10.3390/sym13050784

Chicago/Turabian Style

Sunthrayuth, Pongsakorn, Ahmed M. Zidan, Shao-Wen Yao, Rasool Shah, and Mustafa Inc. 2021. "The Comparative Study for Solving Fractional-Order Fornberg–Whitham Equation via ρ-Laplace Transform" Symmetry 13, no. 5: 784. https://doi.org/10.3390/sym13050784

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