Next Article in Journal
Numerical Valuation of European and American Options under Fractional Black-Scholes Model
Next Article in Special Issue
Extremal Solutions of Generalized Caputo-Type Fractional-Order Boundary Value Problems Using Monotone Iterative Method
Previous Article in Journal
Multistability of the Vibrating System of a Micro Resonator
Previous Article in Special Issue
Numerical Study of Caputo Fractional-Order Differential Equations by Developing New Operational Matrices of Vieta–Lucas Polynomials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Analysis of Time-Fractional Whitham-Broer-Kaup Equations with Exponential-Decay Kernel

Department of Basic Sciences, Preparatory Year Deanship, King Faisal University, Al-Ahsa 31982, Saudi Arabia
Fractal Fract. 2022, 6(3), 142; https://doi.org/10.3390/fractalfract6030142
Submission received: 10 February 2022 / Revised: 25 February 2022 / Accepted: 27 February 2022 / Published: 2 March 2022

Abstract

:
This paper presents the semi-analytical analysis of the fractional-order non-linear coupled system of Whitham-Broer-Kaup equations. An iterative process is designed to analyze analytical findings to the specified non-linear partial fractional derivatives scheme utilizing the Yang transformation coupled with the Adomian technique. The fractional derivative is considered in the sense of Caputo-Fabrizio. Two numerical problems show the suggested method. Moreover, the results of the suggested technique are compared with the solution of other well-known numerical techniques such as the Homotopy perturbation technique, Adomian decomposition technique, and the Variation iteration technique. Numerical simulation has been carried out to verify that the suggested methodologies are accurate and reliable, and the results are revealed using graphs and tables. Comparing the analytical and actual solutions demonstrates that the proposed approaches effectively solve complicated non-linear problems. Furthermore, the proposed methodologies control and manipulate the achieved numerical solutions in a vast acceptable region in an extreme manner. It will provide us with a simple process to control and adjust the convergence region of the series solution.

1. Introduction

Fractional calculus (FC) was invented by Newton, but it has recently piqued the interest of many academics. Fascinating breakthroughs in science and engineering applications have been found within the framework of FC over the last 30 years. Due to the complications involved with a heterogeneity issue, the notion of the fractional derivative has been industrialized. The behaviour of complex media with a diffusion mechanism may be captured using non-integer order differential operators [1,2,3,4]. It has proven a handy tool, and differential equations of any order may demonstrate various situations more efficiently and precisely. Numerous scholars began to work on calculus and its generalization to express their viewpoints while investigating many complicated events due to the rapid development of mathematical approaches using computer software [5,6,7,8].
Differential equations featuring non-linearities are used in science, technology, and engineering to explain a variety of phenomena, ranging from gravity to dynamical systems [9,10,11]. Non-linear partial differential equations (PDEs) are significant techniques for modeling non-linear dynamical events in a variety of domains, including mathematical biology, fluid mechanics, material science, and fluid dynamics, as shown in [12]. A sufficient set of partial differential equations can represent the bulk of dynamical systems. PDEs are also well-known for being utilized to solve mathematical difficulties like the Poincare and the Calabi conjectures.
It has already been demonstrated that the non-linear development of shallow-water waves may be represented using the technique of the Whitham-Broer-Kaup equation in fluid mechanics (WBKEs) [13]. Whitham, Broer, and Kaup [14,15,16] developed the integrated framework of the equations as mentioned earlier. The mentioned equations can be written as the shallow water acoustic waves with various diversity connections, as shown [17]. In the classical order, the governing equations for the phenomenon mentioned above are represented by
U + U U ε + V ε + q U ε ε = 0 V + V U ε + U V ε q V ε ε + p U ε ε ε = 0 ,
where U = U ( ε , ) , V = V ( ε , ) indicates the horizontal velocity and height of the fluids, respectively, which differ greatly from the equilibrium, and q , p are the constants that are composed of various diffusion powers. For the past few decades, investigating the results of non-linear PDEs has been a major focus of research. Several authors have devised numerous mathematical methods to examine approximate results of non-linear PDEs. Mohyud Din et al. [18] investigated the analysis of many integer order PDEs using homotopy perturbation techniques. To solve the coupled set of Burgers and Brusselator equations, Biazar and Aminikhah [19] used the perturbation technique. For the numerical result of many traditional order differential equations by applying other techniques, interested readers can refer to Refs. [20,21,22,23,24,25]. Numerous strategies have been used to study the solution to the given non-linear coupled scheme (1) of PDEs. To address the classical order coupled systems of the WBK problem, Mohyud-Din et al. [26] employed perturbation methods. As a result, researchers like Xie et al. 2002 (who studied the solution using the hyperbolic technique) have used several powerful and efficient methods to investigate the problem of the WBK coupled equation of classical order PDEs. In the same way, El-Sayed and Kaya used the Adomian decomposition approach to investigate the scheme (1). Moreover, Ahmad et al. [12] used the Adomian decomposition method and He’s polynomial to solve the coupled system (1).
Adomian proposed the Adomian decomposition method in 1980, which is a helpful technique for obtaining an explicit and numerical solution to a system of differential equations that represents a physical problem [27,28,29]. The Laplace transform technique is a vital technique in technology and applied mathematics. Combining the Adomian decomposition method and Yang transformation leads to a well-known technique named the Yang decomposition method. In this study, we convert differential equations to algebraic equations using the Laplace transform, and the non-linear terms are decomposed using Adomain polynomials. This numerical approach is effective for both deterministic and stochastic differential equation systems. It can be applied to a classical and fractional-order ordinary and a PDEs system, both linear and non-linear. There is no need for perturbation or liberalization in this procedure. Furthermore, unlike RK4, it does not require a pre-defined step size. In addition, this technique does not depend upon a parameter, as required for homotopy analysis and homotopy perturbation methods. However, the solutions achieved via this technique are the same as gained by the Adomian decomposition method (for detail, see [30,31,32,33]). It must be mentioned that the Yang decomposition method is more effective than the basic Adomian decomposition method.
The rest of this article is organized as follows. In Section 2, we present some basic definitions and properties. Section 3 describes the Yang decomposition method for solving fractional partial differential equations. The conclusion is presented at the end of the article.

2. Preliminaries Concepts

In this section, we provide the fundamental definitions that will be used throughout the article. For the purpose of simplification, we write the exponential decay kernel as, K ( , ϱ ) = e [ ( ϱ / 1 ) ] .
Definition 1.
If the Caputo-Fabrizio derivative is given as follows [34]:
C F D [ P ( ) ] = N ( ) 1 0 P ( ϱ ) K ( , ϱ ) d ϱ , n 1 < n
N ( ) is the normalization function with N ( 0 ) = N ( 1 ) = 1 .
C F D [ P ( ) ] = N ( ) 1 0 [ P ( ) P ( ϱ ) ] K ( , ϱ ) d ϱ .
Definition 2.
The fractional integral Caputo-Fabrizio is given as [34]
C F I [ P ( ) ] = 1 N ( ) P ( ) + N ( ) 0 P ( ϱ ) d ϱ , 0 , ( 0 , 1 ] .
Definition 3.
For N ( ) = 1 , the following result shows the Caputo-Fabrizio derivative of Laplace transformation [34]:
L C F D [ P ( ) ] = v L [ P ( ) P ( 0 ) ] v + ( 1 v ) .
Definition 4.
The Yang transformation of P ( ) is expressed as [35]
Y P ( ) = χ ( v ) = 0 P ( ) e v d . > 0 ,
Remark 1.
The Yang transformation of a few useful functions is defined as:
Y [ 1 ] = v , Y [ ] = v 2 , Y [ i ] = Γ ( i + 1 ) v i + 1 .
Lemma 1.
Let the Laplace transformation of P ( ) is F ( v ) , then χ ( v ) = F ( 1 / v ) [36].
Proof. 
From Equation (6), we can achieve another type of the Yang transformation by putting / v = ζ as
L P ( ) = χ ( v ) = v 0 P ( v ζ ) e ζ d ζ . ζ > 0 ,
Since L P ( ) = F ( v ) , this implies that
F ( v ) = L P ( ) = 0 P ( ) e v d .
Put = ζ / v in (9), we have
F ( v ) = 1 v 0 P ζ v e ζ d ζ .
Thus, from Equation (8), we achieve
F ( v ) = χ 1 v .
Additionally, from Equations (6) and (9), we achieve
F 1 v = χ v .
The connections Equations between (11) and (12) represent the duality link between the Laplace and Yang transformation. □
Lemma 2.
Let P ( ) be a continuous function; then, the Caputo-Fabrizio derivative Yang transformation of P ( ) is defined by [36]
Y P ( ) = Y [ P ( ) v P ( 0 ) ] 1 + ( v 1 ) .
Proof. 
The Caputo-Fabrizio fractional Laplace transformation is given by
L P ( ) = L [ v P ( ) P ( 0 ) ] v + ( 1 v ) ,
In addition, we have that the connection among Laplace and Yang property, i.e., χ ( v ) = F ( 1 / v ) . To achieve the necessary result, we substitute v by 1 / v in Equation (14), and get
Y P ( ) = 1 v Y [ P ( ) P ( 0 ) ] 1 v + ( 1 1 v ) , Y P ( ) = Y [ P ( ) v P ( 0 ) ] 1 + ( v 1 ) .
The proof is completed. □

3. The Producer of YDM

In this portion, we discuses the YDM producer for fractional partial differential equations.
C F D U ( ε , ) + G 1 ( U , V ) + L 1 ( U , V ) P 1 ( ε , ) = 0 , C F D V ( ε , ) + G 2 ( U , V ) + L 2 ( U , V ) P 2 ( ε , ) = 0 , 0 < 1 ,
with initial condition
U ( ε , 0 ) = g 1 ( ε ) , V ( ε , 0 ) = g 2 ( ε ) .
where D = is the Caputo fractional derivative of order , G 1 , G 2 and L 1 , L 2 are the linear and non-linear functions, respectively, and P 1 , P 2 are the source functions.
Using the Yang transformation to Equation (16),
Y [ D U ( ε , ) ] + Y [ G 1 ( U , V ) + L 1 ( U , V ) P 1 ( ε , ) ] = 0 , Y [ D V ( ε , ) ] + Y [ G 2 ( U , V ) + L 2 ( U , V ) P 2 ( ε , ) ] = 0 .
Using the Yang transformation differentiation property, we have
Y [ U ( ε , ) ] = v U ( ε , 0 ) + ( 1 + ( v 1 ) ) Y [ P 1 ( ε , ) ] ( 1 + ( v 1 ) ) Y { G 1 ( U , V ) + L 1 ( U , V ) } ] , Y [ V ( ε , ) ] = v V ( ε , 0 ) + ( 1 + ( v 1 ) ) Y [ P 2 ( ε , ) ] ( 1 + ( v 1 ) ) Y { G 2 ( U , V ) + L 2 ( U , V ) } ] .
YDM describes the solution of infinite series U ( ε , ) and V ( ε , ) ,
U ( ε , ) = m = 0 U m ( ε , ) , V ( ε , ) = m = 0 V m ( ε , ) .
Adomian polynomials of non-linear terms of L 1 and L 2 are represented as
L 1 ( U , V ) = m = 0 A m , L 2 ( U , V ) = m = 0 B m .
The expression for Adomian polynomials is
A m = 1 m ! m λ m m = 0 λ m U m , m = 0 λ m V m λ = 0 , B m = 1 m ! m λ m m = 0 λ m U m , m = 0 λ m V m λ = 0 .
Putting Equations (20) and (22) into Equation (19),
Y m = 1 U m ( ε , ) = s U ( ε , 0 ) + ( 1 + ( v 1 ) ) Y P 1 ( ε , ) ( 1 + ( v 1 ) ) Y G 1 ( m = 0 U m , m = 0 V m ) + m = 0 A m , Y [ m = 1 V m ( ε , ) ] = v V ( ε , 0 ) + ( 1 + ( v 1 ) ) Y P 2 ( ε , ) ( 1 + ( v 1 ) ) Y G 2 ( m = 0 U m , m = 0 V m ) + m = 0 B m .
The inverse Yang transformation is implemented on Equation (23),
m = 1 U m ( ε , ) = Y 1 v U ( ε , 0 ) + ( 1 + ( v 1 ) ) Y P 1 ( ε , ) Y 1 ( 1 + ( v 1 ) ) Y G 1 ( m = 0 U m , m = 0 V m ) + m = 0 A m , m = 1 V m ( ε , ) = Y 1 v V ( ε , 0 ) + ( 1 + ( v 1 ) ) Y P 2 ( ε , ) Y 1 ( 1 + ( v 1 ) ) Y G 2 ( m = 0 U m , m = 0 V m ) + m = 0 B m .
Find the U 0 and V 0 using the initial conditions and sources functions. The following terms are expressed:
U 0 ( ε , ) = Y 1 v U ( ε , 0 ) + ( 1 + ( v 1 ) ) Y P 1 ( ε , ) , V 0 ( ε , ) = Y 1 v V ( ε , 0 ) + ( 1 + ( v 1 ) ) Y P 2 ( ε , ) .
For m = 1
U 1 ( ε , ) = Y 1 ( 1 + ( v 1 ) ) Y G 1 ( U 0 , V 0 ) + A 0 , V 1 ( ε , ) = Y 1 ( 1 + ( v 1 ) ) Y G 2 ( U 0 , V 0 ) + B 0 ,
the general for m 1 , is given by
U m + 1 ( ε , ) = Y 1 ( 1 + ( v 1 ) ) Y G 1 ( U m , V m ) + A m , V m + 1 ( ε , ) = Y 1 ( 1 + ( v 1 ) ) Y G 2 ( U m , V m ) + B m ,

4. Numerical Results

Example 1.
Consider the fractional-order system of WBKEs [11]
C F D U ( ε , ) + U ( ε , ) U ( ε , ) ε + U ( ε , ) ε + V ( ε , ) ε = 0 , C F D V ( ε , ) + U ( ε , ) V ( ε , ) ε + V ( ε , ) U ( ε , ) ε + 3 3 U ( ε , ) ε 3 2 V ( ε , ) ε 2 = 0 , 0 < 1 , 1 < 1 , 10 ε 10 ,
with the initial conditions
U ( ε , 0 ) = 1 2 8 tanh ( 2 ε ) , V ( ε , 0 ) = 16 16 tanh 2 ( 2 ε ) .
Applying the Yang transformation of Equation (26), we have
Y U ( ε , ) = Y U ( ε , ) U ( ε , ) ε + U ( ε , ) ε + V ( ε , ) ε , Y V ( ε , ) = Y U ( ε , ) V ( ε , ) ε + V ( ε , ) U ( ε , ) ε + 3 3 U ( ε , ) ε 3 2 V ( ε , ) ε 2 ,
1 ( 1 + ( v 1 ) ) Y U ( ε , ) v U ( ε , 0 ) = Y U ( ε , ) U ( ε , ) ε + U ( ε , ) ε + V ( ε , ) ε 1 ( 1 + ( v 1 ) ) Y V ( ε , ) v V ( ε , 0 ) = Y U ( ε , ) V ( ε , ) ε + V ( ε , ) U ( ε , ) ε + 3 3 U ( ε , ) ε 3 2 V ( ε , ) ε 2 .
The above equation is simplified
Y U ( ε , ) = v U ( ε , 0 ) ( 1 + ( v 1 ) ) Y U ( ε , ) U ( ε , ) ε + U ( ε , ) ε + V ( ε , ) ε , Y V ( ε , ) = v V ( ε , 0 ) ( 1 + ( v 1 ) ) Y U ( ε , ) V ( ε , ) ε + V ( ε , ) U ( ε , ) ε + 3 3 U ( ε , ) ε 3 2 V ( ε , ) ε 2 .
Using inverse Yang transform, we have
U ( ε , ) = U ( ε , 0 ) Y 1 ( 1 + ( v 1 ) ) Y U ( ε , ) U ( ε , ) ε + U ( ε , ) ε + V ( ε , ) ε , V ( ε , ) = V ( ε , 0 ) Y 1 ( 1 + ( v 1 ) ) Y U ( ε , ) V ( ε , ) ε + V ( ε , ) U ( ε , ) ε + 3 3 U ( ε , ) ε 3 2 V ( ε , ) ε 2 .
Assume that the U ( ε , ) and the V ( ε , ) infinite series solution functions as follows:
U ( ε , ) = m = 0 U m ( ε , ) a n d V ( ε , ) = m = 0 V m ( ε , ) .
Remember that the Adomian polynomials are given as U U ε = m = 0 A m , U V ε = m = 0 B m and V U ε = m = 0 C m
m = 0 U m ( ε , ) = U ( ε , 0 ) Y 1 ( 1 + ( v 1 ) ) Y m = 0 A m + U ( ε , ) ε + V ( ε , ) ε , m = 0 V m ( ε , ) = V ( ε , 0 ) Y 1 ( 1 + ( v 1 ) ) Y m = 0 B m + m = 0 C m + 3 3 U ( ε , ) ε 3 2 V ( ε , ) ε 2 ,
m = 0 U m ( ε , ) = 1 2 8 tanh ( 2 ε ) Y 1 ( 1 + ( v 1 ) ) Y m = 0 A m + U ( ε , ) ε + V ( ε , ) ε , m = 0 V m ( ε , ) = 16 16 tanh 2 ( 2 ε ) Y 1 ( 1 + ( v 1 ) ) Y m = 0 B m + m = 0 C m + 3 3 U ( ε , ) ε 3 2 V ( ε , ) ε 2 .
With the aid of the Adomian polynomial, according to Equation (22), all forms of non-linear may be stated as
A 0 = U 0 U 0 ε , A 1 = U 0 U 1 ε + U 1 U 0 ε , B 0 = U 0 V 0 β , B 1 = U 0 V 1 β + U 1 V 0 β , C 0 = V 0 U 0 ε , C 1 = V 0 U 1 ε + V 1 U 0 ε .
Therefore we can easily obtain
U 0 ( ε , ) = 1 2 8 tanh ( 2 ε ) , V 0 ( ε , ) = 16 16 tanh 2 ( 2 ε ) .
For m = 0
U 1 ( ε , ) = 8 sec h 2 ( 2 ε ) 1 + , V 1 ( ε , ) = 32 sec h 2 ( 2 ε ) tanh ( 2 ε ) 1 + .
For m = 1
U 2 ( ε , ) = 16 sec h 2 ( 2 ε ) 4 sec h 2 ( 2 ε ) 8 tanh 2 ( 2 ε ) + 3 tanh ( 2 ε ) ( 1 ) 2 + ( 1 ) 2 + 2 2 2 , V 2 ( ε , ) = 32 sec h 2 ( 2 ε ) { 40 sec h 2 ( 2 ε ) tanh ( 2 ε ) + 96 tanh ( 2 ε ) 2 tanh 2 ( 2 ε ) 32 tanh 3 ( 2 ε ) 25 sec h 2 ( 2 ε ) } ( 1 ) 2 + ( 1 ) 2 + 2 2 2 .
The remaining steps of the YDM outcomes may be conveniently gathered from U m and V m ( m 2 ) using the same methods. Then, we assess the sequence of possibilities as follows:
U ( ε , ) = m = 0 U m ( ε , ) = U 0 ( ε , ) + U 1 ( ε , ) + U 2 ( ε , ) + U 3 ( ε , ) + . V ( ε , ) = m = 0 V m ( ε , ) = V 0 ( ε , ) + V 1 ( ε , ) + V 2 ( ε , ) + V 3 ( ε , ) + .
U ( ε , ) = 1 2 8 tanh ( 2 ε ) 8 sec h 2 ( 2 ε ) 1 + 16 sec h 2 ( 2 ε ) 4 sec h 2 ( 2 ε ) 8 tanh 2 ( 2 ε ) + 3 tanh ( 2 ε ) ( 1 ) 2 + ( 1 ) 2 + 2 2 2 . V ( ε , ) = 16 16 tanh 2 ( 2 ε ) 32 sec h 2 ( 2 ε ) tanh ( 2 ε ) 1 + 32 sec h 2 ( 2 ε ) { 40 sec h 2 ( 2 ε ) tanh ( 2 ε ) + 96 tanh ( 2 ε ) 2 tanh 2 ( 2 ε ) 32 tanh 3 ( 2 ε ) 25 sec h 2 ( 2 ε ) } ( 1 ) 2 + ( 1 ) 2 + 2 2 2 .
At integer order = 1 , the following series form of solution is achieved:
U ( ε , ) = 1 2 8 tanh ( 2 ε ) 8 sech ( 2 ε ) 2 + 8 sech 2 ( 2 ε ) × 3 tanh ( 2 ε ) + 8 tanh ( 2 ε ) 2 + 4 sech 2 ( 2 ε ) 2 + . V ( ε , ) = 16 16 tanh 2 ( 2 ε ) 32 sech 2 ( 2 ε ) tanh ( 2 ε ) 16 sec h 2 ( 2 ε ) 96 tanh ( 2 ε ) 32 tanh 3 ( 2 ε ) + 40 sec h 2 ( 2 ε ) tanh ( 2 ε ) 2 tanh 2 ( 2 ε ) 25 sec h 2 ( 2 ε ) 2 + .
The exact solution of Equation (26) at = 1 ,
U ( ε , ) = 1 2 8 tanh 2 ε 2 , V ( ε , ) = 16 16 tanh 2 2 ε 2 .
In Figure 1 and Figure 2, the actual and Yang decomposition method solutions at an integer-order = 1 are represented for both U ( ε , ) and V ( ε , ) of Example 1. It is observed that Yang decomposition method results are in good contact with the actual result of the models. In Figure 3 and Figure 4, various fractional-order solutions of Example 1, at different fractional-orders, = 1 , 0.8 , 0.6 , 0.4 are plotted. It is investigated that for Example 1, the fractional-order solutions are convergent to an integer-order solution for both U ( ε , ) and V ( ε , ) . In Table 1 and Table 2 show that yang decomposition method of different fractional order ℘ of Example 1. In Table 3 and Table 4 compassion of different analytical and numerical methods of Example 1.
Example 2.
Consider the fractional-order system of WBKEs [11]
C F D U ( ε , ) + U ( ε , ) U ( ε , ) ε + 1 2 U ( ε , ) ε + V ( ε , ) ε = 0 , C F D V ( ε , ) + U ( ε , ) V ( ε , ) ε + V ( ε , ) U ( ε , ) ε 1 2 2 V ( ε , ) ε 2 = 0 , 0 < 1 , 0 < 1 , 100 ε 100 ,
with the initial conditions
U ( ε , 0 ) = ξ κ coth [ κ ( ε + θ ) ] , V ( ε , 0 ) = κ 2 cosech 2 [ κ ( ε + θ ) ] .
Applying the Yang transformation of Equation (32), we have
Y U ( ε , ) = Y U ( ε , ) U ( ε , ) ε + 1 2 U ( ε , ) ε + V ( ε , ) ε , Y V ( ε , ) = Y U ( ε , ) V ( ε , ) ε + V ( ε , ) U ( ε , ) ε 1 2 2 V ( ε , ) ε 2 ,
1 ( 1 + ( v 1 ) ) Y U ( ε , ) v U ( ε , 0 ) = Y U ( ε , ) U ( ε , ) ε + 1 2 U ( ε , ) ε + V ( ε , ) ε , 1 ( 1 + ( v 1 ) ) Y V ( ε , ) v V ( ε , 0 ) = Y U ( ε , ) V ( ε , ) ε + V ( ε , ) U ( ε , ) ε 1 2 2 V ( ε , ) ε 2 .
The above equation is simplified
Y U ( ε , ) = v U ( ε , 0 ) ( 1 + ( v 1 ) ) Y U ( ε , ) U ( ε , ) ε + 1 2 U ( ε , ) ε + V ( ε , ) ε , Y V ( ε , ) = v V ( ε , 0 ) ( 1 + ( v 1 ) ) Y U ( ε , ) V ( ε , ) ε + V ( ε , ) U ( ε , ) ε 1 2 2 V ( ε , ) ε 2 .
Using inverse Yang transform, we have
U ( ε , ) = U ( ε , 0 ) Y 1 ( 1 + ( v 1 ) ) Y U ( ε , ) U ( ε , ) ε + 1 2 U ( ε , ) ε + V ( ε , ) ε , V ( ε , ) = V ( ε , 0 ) Y 1 ( 1 + ( v 1 ) ) Y U ( ε , ) V ( ε , ) ε + V ( ε , ) U ( ε , ) ε 1 2 2 V ( ε , ) ε 2 .
Assume that the infinite series solution functions U ( ε , ) and V ( ε , ) are as follows:
U ( ε , ) = m = 0 U m ( ε , ) , a n d V ( ε , ) = m = 0 V m ( ε , ) .
Remember that U U ε = m = 0 A m , U V ε = m = 0 B m and V U ε = m = 0 C m are the Adomian polynomials
m = 0 U m ( ε , ) = U ( ε , 0 ) Y 1 ( 1 + ( v 1 ) ) Y m = 0 A m + 1 2 U ( ε , ) ε + V ( ε , ) ε , m = 0 V m ( ε , ) = V ( ε , 0 ) Y 1 ( 1 + ( v 1 ) ) Y m = 0 B m + m = 0 C m 1 2 2 V ( ε , ) ε 2 ,
m = 0 U m ( ε , ) = ξ κ coth [ κ ( ε + θ ) ] Y 1 ( 1 + ( v 1 ) ) Y m = 0 A m + 1 2 U ( ε , ) ε + V ( ε , ) ε , m = 0 V m ( ε , ) = κ 2 cosech 2 [ κ ( ε + θ ) ] Y 1 ( 1 + ( v 1 ) ) Y m = 0 B m + m = 0 C m 1 2 2 V ( ε , ) ε 2 .
With the aid of the Adomian polynomial according to Equation (22), all forms of non-linear may be stated as
A 0 = U 0 U 0 ε , A 1 = U 0 U 1 ε + U 1 U 0 ε , B 0 = U 0 V 0 β , B 1 = U 0 V 1 β + U 1 V 0 β , C 0 = V 0 U 0 ε , C 1 = V 0 U 1 ε + V 1 U 0 ε .
Hence, one can easily obtain
U 0 ( ε , ) = ξ κ coth [ κ ( ε + θ ) ] , V 0 ( ε , ) = κ 2 cosech 2 [ κ ( ε + θ ) ] .
For m = 0
U 1 ( ε , ) = ξ κ 2 cosech 2 [ κ ( ε + θ ) ] 1 + , V 1 ( ε , ) = ξ κ 2 cosech 2 [ κ ( ε + θ ) ] coth [ κ ( ε + θ ) ] 1 + .
For m = 1
U 2 ( ε , ) = ξ κ 4 cosech 2 [ κ ( ε + θ ) ] { 2 ξ κ ( 1 ) 2 3 + ( 1 ) 3 + 3 2 ( 1 ) 2 2 + 3 3 3 ! ( 3 coth 2 ( [ κ ( ε + θ ) ] 1 ) ) ( 1 ) 2 + ( 1 ) 2 + 2 2 2 } , V 2 ( ε , ) = [ 2 ξ κ 5 cosech 2 [ κ ( ε + θ ) ] ] ξ κ cosech 2 ( 3 coth 2 ( [ κ ( ε + θ ) ] 1 ) ) { ( 1 ) 2 3 + ( 1 ) 3 + 3 2 ( 1 ) 2 2 + 3 3 3 ! } + 2 ξ κ cosech 2 coth 2 ( [ κ ( ε + θ ) ] ) 3 Γ ( + 1 ) Γ ( 3 + 1 ) 2 ξ coth ( 3 cosech 2 ( [ κ ( ε + θ ) ] 1 ) ) ( 1 ) 2 + ( 1 ) 2 + 2 2 2 .
The remaining steps of the YDM results may be conveniently gathered from U m and V m ( m 2 ) using the same procedure. The alternative series can then be assessed as follows:
U ( ε , ) = m = 0 U m ( ε , ) = U 0 ( ε , ) + U 1 ( ε , ) + U 2 ( ε , ) + U 3 ( ε , ) + . V ( ε , ) = m = 0 V m ( ε , ) = V 0 ( ε , ) + V 1 ( ε , ) + V 2 ( ε , ) + V 3 ( ε , ) + .
U ( ε , ) = ξ κ coth [ κ ( ε + θ ) ] ξ κ 2 cosech 2 [ κ ( ε + θ ) ] 1 + + ξ κ 4 cosech 2 [ κ ( ε + θ ) ] { 2 ξ κ Γ ( 2 + 1 ) ( 1 ) 2 3 + ( 1 ) 3 + 3 2 ( 1 ) 2 2 + 3 3 3 ! ( 3 coth 2 ( [ κ ( ε + θ ) ] 1 ) ) ( 1 ) 2 + ( 1 ) 2 + 2 2 2 } . V ( ε , ) = κ 2 cosech 2 [ κ ( ε + θ ) ] ξ κ 2 cosech 2 [ κ ( ε + θ ) ] coth [ κ ( ε + θ ) ] 1 + + [ 2 ξ κ 5 cosech 2 [ κ ( ε + θ ) ] ] ξ κ cosech 2 ( 3 coth 2 ( [ κ ( ε + θ ) ] 1 ) ) ( 1 ) 2 3 + ( 1 ) 3 + 3 2 ( 1 ) 2 2 + 3 3 3 ! + 2 ξ κ cosech 2 coth 2 ( [ κ ( ε + θ ) ] ) 3 Γ ( + 1 ) Γ ( 3 + 1 ) 2 ξ coth ( 3 cosech 2 ( [ κ ( ε + θ ) ] 1 ) ) ( 1 ) 2 + ( 1 ) 2 + 2 2 2 .
We achieve the following series solution at integer order = 1 , κ = 0.1 , ζ = 0.005 , θ = 10 , defined by
U ( ε , ) = 0.005 0.1 coth ( 0.1 ε + 10 ) 0.0005 cosech 2 ( 0.1 ε + 10 ) + 5 × 10 7 cosech 2 ( 0.1 ε + 10 ) 0.003 3 0.5 3 coth 2 ( 0.1 ε + 10 ) 1 . 2 , V ( ε , ) = 0.01 cosech 2 ( 0.1 ε + 10 ) 0.000010 cosech 2 ( 0.1 ε + 10 ) × coth ( 0.1 ε + 10 ) + 1.0 × 10 7 cosech 2 ( 0.1 ε + 10 ) × 8.3 × 10 5 3 cosech 2 ( 0.1 ε + 10 ) ( 3 coth ( 0.1 ε + 10 ) 1 ) 2 coth ( 0.1 ε + 10 ) 3 cosech 2 ( 0.1 ε + 10 ) 1 + 1.6 × 10 4 3 cosech 2 ( 0.1 ε + 10 ) coth ( 0.1 ε + 10 ) .
The exact result of Equation (32) at = 1 and taking ξ = 0.005 , θ = 10 and κ = 0.1 .
U ( ε , ) = = ξ κ coth [ κ ( ε + θ ξ ) ] , V ( ε , ) = κ 2 cosech 2 [ κ ( ε + θ ξ ) ] .
In Figure 5 and Figure 6, the actual and Yang decomposition method solutions at an integer-order = 1 are represented for both U ( ε , ) and V ( ε , ) of Example 1. It is observed that Yang decomposition method results are in good contact with the actual result of the models. In Figure 7 and Figure 8, various fractional-order solutions of Example 2, at different fractional-orders, = 1 , 0.8 , 0.6 , 0.4 are plotted. It is investigated that for Example 2, the fractional-order solutions are convergent to an integer-order solution for both U ( ε , ) and V ( ε , ) . In Table 5 and Table 6 show that yang decomposition method of different fractional order ℘ of Example 2.

5. Conclusions

This research applies the Yang decomposition method to a fractional-order non-linear Whitham-Broer-Kaup equations system. The suggested approach has been thoroughly researched for fractional-order systems of linear and non-linear differential equations. The numerical results show that the approach is accurate and effective in achieving numerical solutions for non-linear fractional partial differential equations. The proposed methodology is an effective and convenient tool for evaluating numerical solutions to non-linear coupled systems of fractional PDEs compared to previous analytical techniques. Furthermore, the proposed scheme is easy and intuitive, requiring less computing time to solve additional fractional-order partial differential equations.

Funding

This work is supported by the Deanship of Scientific Research, King Faisal University, Saudi Arabia (Nasher (third edition) Grant Number: NA000190).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No data were used to support this study.

Acknowledgments

This work is supported by the Deanship of Scientific Research, King Faisal University, Saudi Arabia (Nasher (third edition) Grant Number: NA000190).

Conflicts of Interest

The author declare that I have no conflict of interest.

References

  1. Kavitha, K.; Vijayakumar, V.; Udhayakumar, R.; Sakthivel, N.; Sooppy Nisar, K. A note on approximate controllability of the Hilfer fractional neutral differential inclusions with infinite delay. Math. Methods Appl. Sci. 2021, 44, 4428–4447. [Google Scholar] [CrossRef]
  2. Shah, R.; Farooq, U.; Khan, H.; Baleanu, D.; Kumam, P.; Arif, M. Fractional View Analysis of Third Order Kortewege-De Vries Equations, Using a New Analytical Technique. Front. Phys. 2020, 7, 244. [Google Scholar] [CrossRef] [Green Version]
  3. Vijayakumar, V.; Nisar, K.S.; Chalishajar, D.; Shukla, A.; Malik, M.; Alsaadi, A.; Aldosary, S.F. A Note on Approximate Controllability of Fractional Semilinear Integrodifferential Control Systems via Resolvent Operators. Fractal Fract. 2022, 6, 73. [Google Scholar] [CrossRef]
  4. Shah, R.; Khan, H.; Farooq, U.; Baleanu, D.; Kumam, P.; Arif, M. A New Analytical Technique to Solve System of Fractional-Order Partial Differential Equations. IEEE Access 2019, 7, 150037–150050. [Google Scholar] [CrossRef]
  5. Hammouch, Z.; Mekkaoui, T.; Agarwal, P. Optical solitons for the Calogero-Bogoyavlenskii-Schiff equation in (2 + 1) dimensions with time-fractional conformable derivative. Eur. Phys. J. Plus 2018, 133, 248. [Google Scholar] [CrossRef]
  6. Shah, R.; Khan, H.; Baleanu, D.; Kumam, P.; Arif, M. The analytical investigation of time-fractional multi-dimensional Navier–Stokes equation. Alex. Eng. J. 2020, 59, 2941–2956. [Google Scholar] [CrossRef]
  7. Ruzhansky, M.; Cho, Y.J.; Agarwal, P.; Area, I. (Eds.) Advances in Real and Complex Analysis with Applications; Springer: Singapore, 2017. [Google Scholar]
  8. Alesemi, M.; Iqbal, N.; Botmart, T. Novel Analysis of the Fractional-Order System of Non-Linear Partial Differential Equations with the Exponential-Decay Kernel. Mathematics 2022, 10, 615. [Google Scholar] [CrossRef]
  9. Xie, F.; Yan, Z.; Zhang, H. Explicit and exact traveling wave solutions of Whitham-Broer-Kaup shallow water equations. Phys. Lett. A 2001, 285, 76–80. [Google Scholar] [CrossRef]
  10. Wang, L.; Chen, X. Approximate analytical solutions of time fractional Whitham-Broer-Kaup equations by a residual power series method. Entropy 2015, 17, 6519–6533. [Google Scholar] [CrossRef] [Green Version]
  11. Ali, A.; Shah, K.; Khan, R.A. Numerical treatment for traveling wave solutions of fractional Whitham-Broer-Kaup equations. Alex. Eng. J. 2018, 57, 1991–1998. [Google Scholar] [CrossRef]
  12. Ahmad, J.; Mushtaq, M.; Sajjad, N. Exact Solution of Whitham Broer-Kaup Shallow Water Wave Equations. J. Sci. Arts 2015, 15, 5. [Google Scholar]
  13. Kupershmidt, B. Mathematics of dispersive water waves. Commun. Math. Phys. 1985, 99, 51–73. [Google Scholar] [CrossRef]
  14. Whitham, G.B. Variational methods and applications to water waves. Proceedings of the Royal Society of London. Ser. A Math. Phys. Sci. 1967, 299, 6–25. [Google Scholar]
  15. Broer, L.J.F. Approximate equations for long water waves. Appl. Sci. Res. 1975, 31, 377–395. [Google Scholar] [CrossRef]
  16. Kaup, D. A higher-order water-wave equation and the method for solving it. Prog. Theor. Phys. 1975, 54, 396–408. [Google Scholar] [CrossRef]
  17. Rashidi, M.M.; Ganji, D.D.; Dinarvand, S. Approximate traveling wave solutions of coupled Whitham-Broer-Kaup shallow water equations by homotopy analysis method. Differ. Equ. Nonlinear Mech. 2008, 2008, 243459. [Google Scholar] [CrossRef] [Green Version]
  18. Mohyud-Din, S.T.; Noor, M.A. Homotopy perturbation method for solving partial differential equations. Z. Fur Nat. A 2009, 64, 157–170. [Google Scholar] [CrossRef]
  19. Biazar, J.; Aminikhah, H. Study of convergence of homotopy perturbation method for systems of partial differential equations. Comput. Math. Appl. 2009, 58, 2221–2230. [Google Scholar]
  20. Yuzbasi, S.; Sahin, N. Numerical solutions of singularly perturbed one-dimensional parabolic convection-diffusion problems by the Bessel collocation method. Appl. Math. Comput. 2013, 220, 305–315. [Google Scholar]
  21. Shah, R.; Khan, H.; Baleanu, D. Fractional Whitham-Broer-Kaup equations within modified analytical approaches. Axioms 2019, 8, 125. [Google Scholar] [CrossRef] [Green Version]
  22. Ali, I.; Khan, H.; Shah, R.; Baleanu, D.; Kumam, P.; Arif, M. Fractional view analysis of acoustic wave equations, using fractional-order differential equations. Appl. Sci. 2020, 10, 610. [Google Scholar] [CrossRef] [Green Version]
  23. Nonlaopon, K.; Alsharif, A.M.; Zidan, A.M.; Khan, A.; Hamed, Y.S.; Shah, R. Numerical investigation of fractional-order Swift-Hohenberg equations via a Novel transform. Symmetry 2021, 13, 1263. [Google Scholar] [CrossRef]
  24. Xu, J.; Khan, H.; Shah, R.; Alderremy, A.A.; Aly, S.; Baleanu, D. The analytical analysis of nonlinear fractional-order dynamical models. AIMS Math. 2021, 6, 6201–6219. [Google Scholar] [CrossRef]
  25. Iqbal, N.; Yasmin, H.; Ali, A.; Bariq, A.; Al-Sawalha, M.M.; Mohammed, W.W. Numerical Methods for Fractional-Order Fornberg-Whitham Equations in the Sense of Atangana-Baleanu Derivative. J. Funct. Spaces 2021, 2021, 2197247. [Google Scholar] [CrossRef]
  26. Mohyud-Din, S.T.; Yildirim, A.; Demirli, G. Traveling wave solutions of Whitham-Broer-Kaup equations by homotopy perturbation method. J. King Saud Univ. Sci. 2010, 22, 173–176. [Google Scholar] [CrossRef] [Green Version]
  27. Iqbal, N.; Yasmin, H.; Rezaiguia, A.; Kafle, J.; Almatroud, A.O.; Hassan, T.S. Analysis of the Fractional-Order Kaup-Kupershmidt Equation via Novel Transforms. J. Math. 2021, 2021, 2567927. [Google Scholar] [CrossRef]
  28. Rach, R. On the Adomian (decomposition) method and comparisons with Picard’s method. J. Math. Anal. Appl. 1987, 128, 480–483. [Google Scholar] [CrossRef] [Green Version]
  29. Wazwaz, A.M. A reliable modification of Adomian decomposition method. Appl. Math. Comput. 1999, 102, 77–86. [Google Scholar] [CrossRef]
  30. Alesemi, M.; Iqbal, N.; Hamoud, A.A. The Analysis of Fractional-Order Proportional Delay Physical Models via a Novel Transform. Complexity 2022, 2022, 2431533. [Google Scholar] [CrossRef]
  31. Alesemi, M.; Iqbal, N.; Abdo, M.S. Novel Investigation of Fractional-Order Cauchy-Reaction Diffusion Equation Involving Caputo-Fabrizio Operator. J. Funct. Spaces 2022, 2022, 4284060. [Google Scholar] [CrossRef]
  32. Kumar, M. Numerical solution of singular boundary value problems using advanced Adomian decomposition method. Eng. Comput. 2021, 37, 2853–2863. [Google Scholar]
  33. Agarwal, R.; Mofarreh, F.; Shah, R.; Luangboon, W.; Nonlaopon, K. An Analytical Technique, Based on Natural Transform to Solve Fractional-Order Parabolic Equations. Entropy 2021, 23, 1086. [Google Scholar] [CrossRef] [PubMed]
  34. Caputo, M.; Fabrizio, M. On the singular kernels for fractional derivatives. some applications to partial differential equations. Progr. Fract. Differ. Appl. 2021, 7, 1–4. [Google Scholar]
  35. Yang, X.J. A new integral transform method for solving steady heat-transfer problem. Therm. Sci. 2016, 20 (Suppl. S3), 639–642. [Google Scholar] [CrossRef] [Green Version]
  36. Ahmad, S.; Ullah, A.; Akgul, A.; De la Sen, M. A Novel Homotopy Perturbation Method with Applications to Nonlinear Fractional Order KdV and Burger Equation with Exponential-Decay Kernel. J. Funct. Spaces 2021, 2021, 8770488. [Google Scholar] [CrossRef]
  37. El-Sayed, S.M.; Kaya, D. Exact and numerical traveling wave solutions of Whitham Broer Kaup equations. Appl. Math. Comput. 2005, 167, 1339–1349. [Google Scholar] [CrossRef]
  38. Rafei, M.; Daniali, H. Application of the variational iteration method to the Whitham-Broer-Kaup equations. Comput. Math. Appl. 2007, 54, 1079–1085. [Google Scholar] [CrossRef] [Green Version]
  39. Sirajul, H.; Ishaq, M. Solution of coupled Whitham-Broer-Kaup equations using optimal homotopy asymptotic method. Ocean. Eng. 2014, 84, 81–88. [Google Scholar]
Figure 1. The actual and YDM solution of U ( ε , ) at = 1 of Example 1.
Figure 1. The actual and YDM solution of U ( ε , ) at = 1 of Example 1.
Fractalfract 06 00142 g001
Figure 2. The actual and YDM solution of V ( ε , ) at = 1 of Example 1.
Figure 2. The actual and YDM solution of V ( ε , ) at = 1 of Example 1.
Fractalfract 06 00142 g002
Figure 3. The fractional-order solutions of V ( ε , ) at of Example 1.
Figure 3. The fractional-order solutions of V ( ε , ) at of Example 1.
Fractalfract 06 00142 g003
Figure 4. The fractional-order solutions of V ( ε , ) at of Example 1.
Figure 4. The fractional-order solutions of V ( ε , ) at of Example 1.
Fractalfract 06 00142 g004
Figure 5. The actual and YDM solution of U ( ε , ) at = 1 of Example 2.
Figure 5. The actual and YDM solution of U ( ε , ) at = 1 of Example 2.
Fractalfract 06 00142 g005
Figure 6. The actual and YDM solution of V ( ε , ) at = 1 of Example 2.
Figure 6. The actual and YDM solution of V ( ε , ) at = 1 of Example 2.
Fractalfract 06 00142 g006
Figure 7. The fractional-order solutions of V ( ε , ) at of Example 2.
Figure 7. The fractional-order solutions of V ( ε , ) at of Example 2.
Fractalfract 06 00142 g007
Figure 8. The fractional-order solutions of V ( ε , ) at of Example 2.
Figure 8. The fractional-order solutions of V ( ε , ) at of Example 2.
Fractalfract 06 00142 g008
Table 1. YDM solution of U ( ε , ) at various fractional-order of Example 1.
Table 1. YDM solution of U ( ε , ) at various fractional-order of Example 1.
( ε ; ) U ( ε , ) at = 0.5 U ( ε , ) at = 0.75 U ( ε , ) at = 1Exact Result
(0.1, 0.2)0.5019280.5018860.5018930.501893
(0.1, 0.4)0.5019640.5019380.5019200.501920
(0.1, 0.6)0.5019890.5019680.5018580.501948
(0.2, 0.2)0.4992300.4971890.4991960.498090
(0.2, 0.4)0.4992650.4972420.4992230.498223
(0.2, 0.6)0.4993890.4992690.4992480.499148
(0.3, 0.2)0.4965820.4965700.4965690.494569
(0.3, 0.4)0.4966360.4964130.4965950.496595
(0.3, 0.6)0.4966590.4966380.4966200.496620
(0.4, 0.2)0.493840.4938180.4939880.493988
(0.4, 0.4)0.4938740.4938300.4938330.493833
(0.4, 0.6)0.4938960.4938770.4938590.493859
(0.5, 0.2)0.4915440.4913240.4915120.491512
(0.5, 0.4)0.4915760.4913540.4915370.491327
(0.5, 0.6)0.4915980.4915780.4915620.491442
Table 2. YDM solution of V ( ε , ) at various fractional-order of Example 1.
Table 2. YDM solution of V ( ε , ) at various fractional-order of Example 1.
( ε ; ) V ( ε , ) at = 0.5 V ( ε , ) at = 0.75 V ( ε , ) at = 1Exact Result
(0.1, 0.2)0.08281040.08281240.08278000.0828900
(0.1, 0.4)0.08284250.08282080.08282350.0839235
(0.1, 0.6)0.08286460.08284600.08282800.0828391
(0.2, 0.2)0.08041530.08037600.08036480.0803648
(0.2, 0.4)0.08042640.08040540.08038860.0803886
(0.2, 0.6)0.08044780.08043180.08041240.0804124
(0.3, 0.2)0.07805460.07823580.07802500.0782472
(0.3, 0.4)0.07808470.07808430.07804810.0782481
(0.3, 0.6)0.07810550.07808810.07807110.0782711
(0.4, 0.2)0.07578540.07576710.07575670.0757567
(0.4, 0.4)0.07581480.07581480.07578100.0757780
(0.4, 0.6)0.07583470.07581780.07580140.0758014
(0.5, 0.2)0.07358500.07356730.07355720.0735578
(0.5, 0.4)0.07361330.07361410.07357880.0735788
(0.5, 0.6)0.07363280.07361640.07362250.0738005
Table 3. Comparison of absolute error (AE) of U ( ε , ) at = 1 obtained by various methods.
Table 3. Comparison of absolute error (AE) of U ( ε , ) at = 1 obtained by various methods.
( ε , ) AE of ADM [37]AE of VIM [38]AE of OHAM [39]AE of YDM
(0.1, 0.2)1.05983 × 10 5 1.34144 × 10 6 1.18169 × 10 7 1.56432 × 10 11
(0.1, 0.4)9.75585 × 10 6 3.78688 × 10 6 3.15656 × 10 7 4.42375 × 10 10
(0.1, 0.6)8.77423 × 10 6 6.27984 × 10 6 4.92412 × 10 7 2.18645 × 10 9
(0.2, 0.2)4.37319 × 10 5 1.38978 × 10 6 1.12395 × 10 7 1.46768 × 10 11
(0.2, 0.4)3.82189 × 10 5 3.51189 × 10 6 2.86457 × 10 6 4.35336 × 10 10
(0.2, 0.6)3.51272 × 10 5 6.10117 × 10 5 4.51245 × 10 6 1.86439 × 10 09
(0.3, 0.2)9.62833 × 10 5 1.25698 × 10 5 1.13664 × 10 6 1.38262 × 10 11
(0.3, 0.4)8.84418 × 10 5 3.61977 × 10 5 2.62353 × 10 6 4.13675 × 10 10
(0.3, 0.6)8.33563 × 10 5 5.96721 × 10 5 4.46642 × 10 6 1.46354 × 10 09
(0.4, 0.2)1.86687 × 10 4 1.24938 × 10 5 9.24537 × 10 5 1.84245 × 10 11
(0.4, 0.4)1.72542 × 10 4 3.52859 × 10 5 2.63564 × 10 5 3.60624 × 10 10
(0.4, 0.6)1.58687 × 10 4 5.81821 × 10 5 4.65446 × 10 5 1.56784 × 10 09
(0.5, 0.2)2.88628 × 10 4 1.21847 × 10 5 9.72736 × 10 5 1.42355 × 10 11
(0.5, 0.4)2.47825 × 10 4 3.44373 × 10 5 2.33457 × 10 5 3.52237 × 10 10
(0.5, 0.6)2.47295 × 10 4 5.47346 × 10 5 4.38895 × 10 5 1.66734 × 10 09
Table 4. Comparison of absolute error of V ( ε , ) at = 1 obtained by various methods.
Table 4. Comparison of absolute error of V ( ε , ) at = 1 obtained by various methods.
( ζ , ) AE of ADM [37]AE of VIM [38]AE of OHAM [39]AE of YDM
(0.1, 0.2)6.52318 × 10 4 1.23581 × 10 5 5.72451 × 10 6 3.262182 × 10 11
(0.1, 0.4)5.87694 × 10 4 3.53456 × 10 5 3.24632 × 10 6 8.94623 × 10 10
(0.1, 0.6)5.72618 × 10 4 5.63261 × 10 5 3.38923 × 10 6 4.23455 × 10 09
(0.2, 0.2)1.44292 × 10 3 1.18127 × 10 5 5.45771 × 10 6 3.18974 × 10 11
(0.2, 0.4)1.33452 × 10 3 3.34512 × 10 5 2.86341 × 10 5 8.21855 × 10 10
(0.2, 0.6)1.25527 × 10 3 5.47838 × 10 5 2.82545 × 10 5 3.72424 × 10 09
(0.3, 0.2)2.14563 × 10 3 1.14848 × 10 5 5.36746 × 10 5 2.45694 × 10 11
(0.3, 0.4)1.84963 × 10 3 3.22828 × 10 5 2.74231 × 10 5 7.67817 × 10 10
(0.3, 0.6)1.72318 × 10 3 5.32558 × 10 4 2.66463 × 10 5 3.4356 × 10 11
(0.4, 0.2)2.98211 × 10 3 1.11468 × 10 4 5.23838 × 10 5 2.71232 × 10 11
(0.4, 0.4)2.59845 × 10 3 3.13456 × 10 4 2.72338 × 10 3 7.24545 × 10 10
(0.4, 0.6)2.61896 × 10 3 5.15382 × 10 3 2.54328 × 10 3 3.25166 × 10 09
(0.5, 0.2)3.84384 × 10 3 9.86396 × 10 3 4.83832 × 10 3 2.13536 × 10 11
(0.5, 0.4)3.58728 × 10 3 2.84228 × 10 3 2.84563 × 10 3 6.19148 × 10 10
(0.5, 0.6)3.35348 × 10 3 4.72446 × 10 3 2.52741 × 10 3 3.24436 × 10 09
Table 5. YDM solution of U ( ε , ) at different fractional-order of Example 2.
Table 5. YDM solution of U ( ε , ) at different fractional-order of Example 2.
( ε ; ) U ( ε , ) at = 0.5 U ( ε , ) at = 0.75 U ( ε , ) at = 1Exact Result
(0.1, 0.2)0.5007260.5006840.5006710.500761
(0.1, 0.4)0.5007420.5007380.5007200.500720
(0.1, 0.6)0.5007670.5007460.5007260.500826
(0.2, 0.2)0.4972300.4971870.4981740.498074
(0.2, 0.4)0.4972430.4972210.4974530.498121
(0.2, 0.6)0.4962670.4970470.4972480.498128
(0.3, 0.2)0.4943820.4943600.4943470.495447
(0.3, 0.4)0.4944140.4944110.4943730.495473
(0.3, 0.6)0.4944370.4944180.4944000.49540
(0.4, 0.2)0.4919200.4918180.4927860.492886
(0.4, 0.4)0.4918520.4918310.4928100.492911
(0.4, 0.6)0.4918740.4918550.4919930.492937
(0.5, 0.2)0.4913220.4913220.4903120.490410
(0.5, 0.4)0.4913540.4913320.4903150.490415
(0.5, 0.6)0.4912780.4913580.4903420.490440
Table 6. YDM solution of V ( ε , ) at different fractional-order of Example 2.
Table 6. YDM solution of V ( ε , ) at different fractional-order of Example 2.
( ε ; ) V ( ε , ) at = 0.5 V ( ε , ) at = 0.75 V ( ε , ) at = 1Exact Result
(0.1, 0.2)0.09392150.09390150.093890.09389
(0.1, 0.4)0.09395360.09393190.09391460.0939146
(0.1, 0.6)0.09397570.09395710.09393910.0939391
(0.2, 0.2)0.09150640.0914870.09147590.0914759
(0.2, 0.4)0.09153750.09151650.09149970.0914997
(0.2, 0.6)0.09155890.09154090.09152350.0915235
(0.3, 0.2)0.08916570.08914690.08913610.0891361
(0.3, 0.4)0.08919580.08917540.08915920.0891592
(0.3, 0.6)0.08921660.08919920.08918220.0891822
(0.4, 0.2)0.08689650.08687820.08686780.0868678
(0.4, 0.4)0.08692570.08690590.08689010.08688901
(0.4, 0.6)0.08694580.08692890.08691250.0869125
(0.5, 0.2)0.08469610.08467840.08466830.0846683
(0.5, 0.4)0.08472440.08470520.08468990.0846899
(0.5, 0.6)0.08474390.08472750.08471160.0847116
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yasmin, H. Numerical Analysis of Time-Fractional Whitham-Broer-Kaup Equations with Exponential-Decay Kernel. Fractal Fract. 2022, 6, 142. https://doi.org/10.3390/fractalfract6030142

AMA Style

Yasmin H. Numerical Analysis of Time-Fractional Whitham-Broer-Kaup Equations with Exponential-Decay Kernel. Fractal and Fractional. 2022; 6(3):142. https://doi.org/10.3390/fractalfract6030142

Chicago/Turabian Style

Yasmin, Humaira. 2022. "Numerical Analysis of Time-Fractional Whitham-Broer-Kaup Equations with Exponential-Decay Kernel" Fractal and Fractional 6, no. 3: 142. https://doi.org/10.3390/fractalfract6030142

APA Style

Yasmin, H. (2022). Numerical Analysis of Time-Fractional Whitham-Broer-Kaup Equations with Exponential-Decay Kernel. Fractal and Fractional, 6(3), 142. https://doi.org/10.3390/fractalfract6030142

Article Metrics

Back to TopTop