#### Abstract

In this paper, we investigate the numerical solution of the Fornberg-Whitham equations with the help of two powerful techniques: the modified decomposition technique and the modified variational iteration technique involving fractional-order derivatives with Mittag-Leffler kernel. To confirm and illustrate the accuracy of the proposed approach, we evaluated in terms of fractional order the projected models. Furthermore, the physical attitude of the results obtained has been acquired for the fractional-order different value graphs. The results demonstrated that the future method is easy to implement, highly methodical, and very effective in analyzing the behavior of complicated fractional-order linear and nonlinear differential equations existing in the related areas of applied science.

#### 1. Introduction

The analysis of the Fornberg-Whitham equation (FWE) is a significant mathematical equation of mathematical physics. The Fornberg-Whitham equation is defined as [1, 2]

This model was invented to evaluate the nonlinear breaking dispersive ocean waves. The Fornberg-Whitham equation is shown to yield peakon solutions as a physical equation for waves of restricting height and the occurrences of wave breaking. Fractional calculus is now widely used and accepted, owing to its well-known uses in a variety of fields of seemingly disparate sectors of science and engineering [3, 4]. Many scholars, including Gupta and Singh [5] and Alderremy et al. [6], have examined the fractional of the Fornberg-Whitham equation relevant to the fractional Caputo derivative, Sunthrayuth et al. [7], Singh et al. [8], etc. Because of the singular kernel of the fractional Caputo derivative, its implementations are limited. Caputo and Fabrizio [9] created derivatives of any (real or complex) order with a nonsingular kernel. Caputo and Fabrizio’s derivative has been used to a variety of real-world situations, including fractional nonhomogeneous heat models [10], El Nino-Southern fractional oscillations models [11], and arbitrary-order system of smoking models [12]. Atangana and Baleanu [13] devised a novel fractional-order derivative called the Atangana-Baleanu (AB) fractional derivative, which has the kernel of a Mittag-Leffler-type function. Kumar et al. [14] investigated the regularised long-wave equation with a Mittag-Leffler-type kernel incorporating the fractional operator. A Mittag-Leffler-type kernel is used in the chemical kinetics system connected with a fractional derivative which was investigated by Singh et al. [15]. Baleanu et al. [16] recently proposed optimal fractional models with nonsingular Mittag-Leffler kernels. As we all know, the Mittag-Leffler function is more beneficial in expressing physical difficulties than the power function or the exponential function; as a result, the AB fractional derivative is well suited to unravelling material heterogeneities and structures or media with different scales.

George Adomian was introduced to and established a technique for “solving integro-differential, differential equations, delay differential, and partial differential equations” [17, 18]. The result is discovered as an infinite sequence that quickly converges to precise solution. This method has been shown to be effective in solving both linear and nonlinear models. The method for solving a nonlinear operator problem is to use decomposition equations in a series of functions. Each expression’s sequence is derived from a polynomial derived from the expansion of an approximate solution into power sequences. The Adomian decomposition method technique is really simple in theory, but the difficulty arises when it comes to determining polynomials and illustrating the convergence of a series of functions [19]. Lesnic [20] analyzed the convergent of the Adomian decomposition method when using heat and wave models for both backward and forward time evolution. Gaber and El-Sayed used the Adomian method of solving fractal-order partial differential equations on a finite domain in [21]. Ghoreishi et al. [22] investigated the Adomian decomposition method’s ability to investigate nonlinear wave problems with changing coefficients, demonstrating that the Adomian decomposition method can solve these equations without the need for dissertation, linearization, transformation, or perturbation.

The variational iteration approach [23, 24] was published in the late 1990s to solve a seepage flow with fractional derivatives and a nonlinear oscillator, and it has since been widely utilized as a primary analytical tool for solving a variety of nonlinear problems. It has fully grown into a fully fledged mathematical approach as a result of considerable research by a number of authors, including He [25, 26], Ganji and Sadighi [27], Ozis and Yildirim [28], and Noor and Mohyud-Din [29]. On November 24, 2018, we searched Clarivate’s Web of Science for “variational iteration approach” and got 3761 hits. The technique’s identification of the Lagrange multiplier necessitates the understanding of variational theory [30], and the technique’s sophisticated identification process may limit its implementation to real-world issues.

#### 2. Preliminary Concepts

*Definition 1. *The Caputo fractional derivative is given as [31]

*Definition 2. *The Laplace transformation connected with fractional Caputo derivative is expressed by [31]

*Definition 3. *In Caputo sense, the Atangana-Baleanu derivative is defined as [31]
where is a normalization function such that , and represent the Mittag-Leffler function.

*Definition 4. *The Atangana-Baleanu derivative in the Riemann-Liouville sense is defined as [31]

*Definition 5. *The Laplace transform connected with the Atangana-Baleanu operator is defined as [31]

*Definition 6. *Consider , and is a function of ; then, the fractional-order integral operator of is given as [31]

#### 3. The Methodology of Variational Iteration Method

This section introduces the solution of fractional partial differential equations with the help of the variational iteration method. The initial condition is where is the fractional derivative Caputo order , and are linear and nonlinear terms, respectively, and is the source function.

The Laplace transformation is applied to equation (8); we get

The Lagrange multiplier iterative method is

A Lagrange multiplier is as

Applying inverse Laplace transform , equation (11) can be written as

#### 4. The Conceptualization of MDM

In this section, we discuss the solution of fractional partial differential equations with the help of the modified decomposition method. The initial condition is where is the fractional derivative of Caputo order , and are linear and nonlinear terms, respectively, and is the source term.

Using Laplace transformation to equation (14), we get

Taking the Laplace transform of differentiation property, we have

MDM result of infinite series ,

nonlinear function is defined as

The nonlinear terms can be analyzed with the aid of Adomian polynomials. So the Adomian polynomial formula is expressed as

Then, put equations (18) and (19) into (17), which gives

Applying the inverse Laplace transformation to equation (21), we get

Define the terms as follows:

In general, is defined as

#### 5. Application of Techniques

*Example 7. *Consider the time-fractional nonlinear Fornberg-Whitham equation
with the initial condition

Taking Laplace transformation of (25),

Using inverse Laplace transformation

Applying Adomian procedure, we have for , for , for ,

The modified decomposition method solution of example (1) is

The simplification of equation (33)

Apply the variational method to obtain series form solution. The iteration formulas for equation (25), we get

where

For ,

The exact solution of equation (25) at ,

In Figure 1, the analytical results of MDM/MVITM example 1 graphs show close contact with each other at and 0.8. It is investigated that analytical results are in close relation with the actual results of example 1. In Figure 2, the results of example 1 at different fractional-order of the derivative are plotted at and 0.4. Figure 3 shows the different fractional of two and three dimensional. The graphical representation has shown the convergence phenomena of fractional-order results towards the result at integer-order of example 1.

*Example 8. *Consider the time-fractional nonlinear Fornberg-Whitham equation is given as
The initial condition is

Applying Laplace transformation of (39), we get

Using inverse Laplace transformation,

Using ADM procedure, we get for , for , for ,

The MDM solution of example (8) is

Apply the variational method to find the analytical solution.

The iteration formulas for equation (39), we get where for ,

The exact solution of equation (39) at ,

In Figure 4, the analytical results of MDM/MVITM example 2 graphs show that close contact with each other at and 0.8. It is investigated that analytical results are in close relation with the actual results of example 2. In Figure 5, the results of example 2 at different fractional-order of the derivative are plotted at and 0.4. Figure 6 shows the different fractional of two and three dimensional. The graphical representation has shown the convergence phenomena of fractional-order results towards the result at integer-order of example 2.

#### 6. Conclusion

In this paper, we have been successfully applied two modified methods to investigate the approximate solutions of fractional Fornberg-Whitham equations. Agreement between numerical results obtained by the modified decomposition method and modified variational iteration method involving fractional-order derivatives with Mittag-Leffler kernel with exact result appears very appreciable by means of illustrative results in figures. The proposed techniques are easy to implement, effective, and suitable for achieving the results of nonlinear fractional Fornberg-Whitham equations. Moreover, both the modified decomposition method and variational iteration method provide the convergent series results with easily calculated components without applying any linearization, perturbation, or limiting assumptions. Finally, we can conclude the suggested methods are more accurate and highly methodical and which can be applied to investigate nonlinear models that arise in applied sciences.

#### Data Availability

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

#### Conflicts of Interest

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

#### Acknowledgments

This research has been funded by Scientific Research Deanship at University of Ha’il, Saudi Arabia, through project number RG-21 005.