Modern Applications of Bioconvection with Fractional Derivatives
View this Special IssueResearch Article  Open Access
Kamsing Nonlaopon, Muhammad Naeem, Ahmed M. Zidan, Rasool Shah, Ahmed Alsanad, Abdu Gumaei, "Numerical Investigation of the TimeFractional Whitham–Broer–Kaup Equation Involving without Singular Kernel Operators", Complexity, vol. 2021, Article ID 7979365, 21 pages, 2021. https://doi.org/10.1155/2021/7979365
Numerical Investigation of the TimeFractional Whitham–Broer–Kaup Equation Involving without Singular Kernel Operators
Abstract
This paper aims to implement an analytical method, known as the Laplace homotopy perturbation transform technique, for the result of fractionalorder Whitham–Broer–Kaup equations. The technique is a mixture of the Laplace transformation and homotopy perturbation technique. Fractional derivatives with MittagLeffler and exponential laws in sense of Caputo are considered. Moreover, this paper aims to show the Whitham–Broer–Kaup equations with both derivatives to see their difference in a realworld problem. The efficiency of both operators is confirmed by the outcome of the actual results of the Whitham–Broer–Kaup equations. Some problems have been presented to compare the solutions achieved with both fractionalorder derivatives.
1. Introduction
In engineering and applied sciences and technology, fractional partial differential equations (FPDEs) containing nonlinearities define many phenomena, ranging from gravitation to dynamics. The nonlinear FPDEs are significant tools analyzed to model nonlinear dynamical behaviour in many areas such as plasma physics, mathematical biology, fluid dynamics, and solidstate physics. The widely held dynamical schemes can be denoted by an appropriate set of FPDEs. Moreover, it is well identified that FPDEs solve mathematical models, such as Poincare conjecture and Calabi conjecture models [1–12].
It has been determined that the nonlinear development of shallowwater waves in fluid mechanics is described by a coupled system of Whitham–Broer–Kaup (WBK) equations [13]. Whitham [14], Broer [15], and Kaup [16] proposed the coupled scheme of the aforementioned equations. The aforementioned equations define shallowwater wave propagation with various spreading relations, as shown in [17]. The governing equations for the respective phenomena in classic order are provided bywhere and show the height and horizontal velocity that diverges from the liquid’s equilibrium position, respectively, and are constant, signified in different diffusion powers. Over the last few decades, there has been a lot of research into solutions to such nonlinear PDEs. So, many researchers have created a variety of mathematical methods to investigate the analytical results of nonlinear PDEs. The HPM was used by Biazar and Khah [18] to solve the coupled schemes of the Burger and Brusselator problems. Amjad et al. [19] applied the solution of standard order coupled with fractionalorder Whitham–Broer–Kaup equation by Laplace decomposition technique. Noor et al. [20] used the homotopy perturbation technique to investigate the results of much classical order of PDEs. Whitham–Broer–Kaup equations are used by other scholars who implemented several numerical techniques, such as residual power series technique [4], reduced differential transformation technique [21], Adomian decomposition technique [7], homotopy perturbation technique [22, 23], Lie symmetry analysis [24, 25], expfunction technique [26], expansion technique [27], and homotopy analysis technique [28].
Fractional calculus (FC) is a new mathematical approach for describing models of nonlocal behaviour. Fractional derivatives have mathematically described many other physical problems in recent years; these representations have yielded excellent outcomes in the simulation of realworld issues. Some basic definitions of fractional operators were given by Riesz, Coimbra, Hadamard, Riemann–Liouville, Grunwald–Letnikov, Weyl, Caputo, Fabrizio, and Atangana Baleanu, among others [29–31]. To investigate the solutions of nonlinear FPDEs, some wellknown techniques for finding actual results have been develop, for example, the homotopy perturbation transform technique [32, 33], the invariant subspace technique [34], the Hermite colocation technique [35], the homotopy analysis transform method [36], the optimal homotopy asymptotic technique [37], the homotopy analysis Sumudu transform technique [38], the Adomian decomposition technique [39], the Pade approximation and homotopyPade method [39], and the Sumudu transform series expansion technique [40]. The Laplace homotopy perturbation transform technique is a mixture of the homotopy perturbation technique introduced by Liao [41] and of the Laplace transformation [42].
The rest of the paper is arranged as follows. Section 2 discusses the basic definitions from fractional calculus. Section 3 is introduced to the fundamental methodology of HPTM. Sections 4 and 5 are the implementations of the techniques for the CF and AN operators. The conclusion of the work is written in Section 6.
2. Fractional Calculus
This section provides some fundamental concepts of fractional calculus.
Definition 1 (see [42]). The Liouville–Caputo operator is given bywhere is the derivative of integer th order of and . If , then we defined the Laplace transformation for the Caputo fractional derivative as follows:
Definition 2 (see [42]). The Caputo–Fabrizio operator (CF) is defined bywhere is a normalization form and . The exponential law is used as the nonsingular kernel in this fractional operator.
If , then we define the Caputo–Fabrizio of Laplace transformation for the fractional derivative is given as
Definition 3 (see [42]). The fractional generalized MittagLeffler law with the sense of Atangana–Baleanu operator is defined as follows:where is a normalization function with . The MittagLeffler law is used as a nonsingular and nonlocal kernel in this fractional operator.
If , then we express the Laplace transformation for the Atangana–Baleanu operator fractional derivative as
3. Implementation of the LHPTM for the Solution of Fractional Partial Differential Equation
The LHPTM is general methodology and can be written as follows.
The main procedure of this technique is defined as follows: Step 1: let us consider the following equation: under the initial condition where is a known function, is the order of the derivative, represents a linear differential operator, and is the general nonlinear differential operator. Step 2: using both sides Laplace transformation operator of (8), we obtain Laplace transformation is applied to Caputo–Fabrizio (5) and Atangana–Baleanu (7) operators. Step 3: on both sides, using the inverse Laplace transformation of equation (10), we obtain Step 4: applying the homotopy producer, the result of the above equations in a series form is defined byand the nonlinear terms can be expressed aswhere is an embedding parameter and are He’s polynomials that can be provided by
Finally, the LHPTM is achieved by coupling the decomposition technique which is defined bywhere .
The terms, comparing with the same powers of , produce results of many orders. The initial estimated of the approximation is , which is actually the Taylor series for the exact result of order .
Using the aforementioned technique, we solved the fractionalorder Whitham–Broer–Kaup equations in the Atangana–Baleanu and Caputo–Fabrizio senses using the LHPM.
4. Implementation of Caputo–Fabrizio Operator
Example 1. Let us consider the coupled system of fractionalorder WBKEs in the CFC sense:under the initial conditions,Using the Laplace transformation to equation (16), we obtainSimplify the above equation and use the initial conditions (17), and we obtainThe inverse Laplace transformation is implemented to (19), and we obtainThe LHPTM is used in (20), and we obtainThe nonlinear can be found with the help of He’s polynomial and can be defined asComparing the coefficient of , we haveWe can calculate few terms of (16) which can be written asThe exact solution of (16) isFigure 1 shows the actual and approximate solutions of at , and Figure 2 shows the actual and approximate solutions of at . Figures 3 and 4 show that the first graph has a different fractional order with respect to and the second graph has a different fractional order with respect to of Example 1. Tables 1 and 2 show the different fractionalorder of and . Tables 3 and 4 show the comparisons with different methods.
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)




Example 2. Let us consider the coupled system of fractionalorder WBKEs in the CFC sense:under the initial conditions,Using the Laplace transformation to equation (26), we obtainSimplify the above equation and use the initial conditions (27), and we obtainThe inverse Laplace transformation is implemented to (29), and we obtainThe LHPTM is used in (30), and we obtainThe nonlinear can be found with the help of He’s polynomial and can be defined asComparing the coefficient of , we haveWe can calculate few terms of (26) which can be written asThe exact solution of (26) is
5. Implementation of Atangana–Baleanu Operator
Example 3. Let us consider the coupled system of fractionalorder WBKEs in the ABC sense:under the initial conditions,Using the Laplace transformation to (36), we obtainSimplify the above equation and use the initial conditions (37), and we obtainThe inverse Laplace transformation implement to (39), we obtainThe LHPTM is used in (40), and we obtainThe nonlinear can be found with the help of He’s polynomial and can be defined asComparing the coefficient of , we haveWe can calculate few terms of (36) which can be written as