Abstract

The role of integer and noninteger order partial differential equations (PDE) is essential in applied sciences and engineering. Exact solutions of these equations are sometimes difficult to find. Therefore, it takes time to develop some numerical techniques to find accurate numerical solutions of these types of differential equations. This work aims to present a novel approach termed as fractional iteration algorithm-I for finding the numerical solution of nonlinear noninteger order partial differential equations. The proposed approach is developed and tested on nonlinear fractional-order Fornberg–Whitham equation and employed without using any transformation, Adomian polynomials, small perturbation, discretization, or linearization. The fractional derivatives are taken in the Caputo sense. To assess the efficiency and precision of the suggested method, the tabulated numerical results are compared with the standard variational iteration method and the exact solution as well. In addition, numerical results for different cases of the fractional-order are presented graphically, which show the effectiveness of the proposed procedure and revealed that the proposed scheme is very effective, suitable for fractional PDEs, and may be viewed as a generalization of the existing methods for solving integer and noninteger order differential equations.

1. Introduction

The nonlinear PDEs have become a hot topic in the field of nonlinear science, which has been used to describe the problems in many fields, such as quantum mechanics, image processing, ecology and economic system, and epidemiology. PDEs are broadly emerging in different physical applications like dispersing and propagation of waves, magnetic resonance imaging, computational fluid dynamics, magnetohydrodynamic move through pipes, phenomena of supersonic and turbulence flow, acoustic transmission, and traffic. More details can be found in [1] and the references therein. PDEs are used in population models, medical imaging, proper distribution of oxygen to the healing tissues, electrical signalling of nerves, etc. [2]. The popularity of PDE has been confirmed in a very actual prediction of the number of COVID cases [3, 4]. Using PDE, it is possible to make a model of the shape of COVID-19 [5]. However, for some complex problems in these fields, the fractional PDE is more accurate than integer-order partial differential equation. So, generating numerical solutions of fractional PDEs has become extremely important.

The noninteger order Calculus was created immediately after the integer order Calculus, which means that it dates back to the 17th century. The notation was originated around 1675 by Leibniz to represent the th derivative of a function , assuming integer values of . In 1695, L′ Hospital stated the problem: “What if is ?” Later, Leibniz used fractional derivatives (FD) symbols in his research. Accordingly, the concept of an FD appeared almost at the same time as the integer order derivative. Many famous researchers in the 18th and 19th century, for example Euler, Lagrange, Laplace, Fourier, and many others, contributed to the development of the fractional calculus (FC). The first problem modeled in terms of FC was the known problem of the tautochrone curve. In 1823, Abel found a solution in terms of an integral equation, which is based on the Riemann-Liouville definition of the fractional integration. However, the popularity of fractional partial differential equations (FPDE) was slowed down by some incomplete and conflicting definitions and difficulties in establishing consistent regulations for inverse operations. Satisfactory definitions for the integration and differentiation with fractional power were developed in the middle of the 19th century. FC is still an active of research in engineering and sciences [610].

Due to the importance of numerical solution of fractional PDEs (FPDEs) in science and engineering, some powerful numerical techniques are developed in the literature, and many prominent researchers have made contributions in this area. These include finite element methods [11], mixed finite element (MFE) methods [12], (local) discontinuous Galerkin methods [13], finite volume methods [14], Jacobi collocation [15], variational iteration method [16], generalized Kudryashov method [17], finite Hankel transform procedure [18], modified Khater method [19], residual power series method [20], modified auxiliary equation method [21], local meshless method [22], RBF collocation method [23], operational matrix of B-spline functions [24], and some other advanced numerical methods as well. Zhang and Xu [25], Lin et al. [26] considered some spectral approximations for the time-fractional water wave model, time-fractional diffusion equation. Most nonlinear TFPDEs do not have accurate exact solutions. Therefore, direct and iterative approaches are applied. Most of the methods for nonlinear TFPDEs give infinite series solutions or mesh-based method. These methods are time-consuming due to repeated calculation in series solution and mesh creation in mesh-based methods.

Opposite to great popularity and applicability, the shortcoming is obvious in the available numerical techniques of reasonably high orders for finding approximate solutions of FPDEs. Such a need for a general method, which is usable in solving linear, nonlinear, homogenous, nonhomogeneous, and multivariable FPDEs problems without major changes, is the inspiration for the current research. Recently, many scholars have investigated the numerical solutions of the fractional PDEs, which has effectively promoted the advancement of the field of nonlinear PDEs. However, in general, numerical techniques have some limitations such as low accuracy, mesh generation, transformations, stability, convergence, and difficulty of implementing in complex geometries. In recent years, the variational iteration methods (VIM) have become popular in the field of numerical approximations. In this technique, the discretization of the domain and linearization of given differential equations is not required. We simply need to calculate the Lagrange multiplier of the given differential equation by restricting the nonlinear terms and in series form analytical solution of the given differential equations can be obtained. VIM can be implemented in the easiest way and is more flexible than other techniques available in the literature. It was applied for the first time to FPDEs by He in [27], and later on by Odibat and Momani in [28]. Inc [29] employed it for the numerical approximation of fractional Burgers equations. Yulita et al. [30] applied VIM for the analytical treatment of fractional Zakharov-Kuznetsov equations, while Safari et al. [31] utilized it for fractional KdV Burger’s Kuramoto equation. Das in [32] investigated fractional diffusion equations by using VIM and obtained the exact solution.

In this study, we present a novel approach, termed as fractional iteration algorithm-I, for solving noninteger order differential equations. The proposed approach is developed and tested on nonlinear fractional-order Fornberg–Whitham equation. This equation is defined assubject to the beneath initial condition

The constant in 1 lies in the interval , is the fluid velocity, is the spatial, and is the time coordinate.

The numerical solution of the Fornberg–Whitham type equations (1) is numerically challenging. To solve these equations numerically, several methods were proposed in the literature. Lu in [33] solved Fornberg–Whitham type equations by He’s variational iteration method. The approximate analytical solutions to this type of equations were presented by Abidi and Omrani in [34] employing homotopy analysis technique. Zhou and Tian [35] obtained a special type of travelling wave solution of the Fornberg–Whitham type equations by applying the bifurcation method. Yin et al. [36] classified all the travelling solutions obtained by an improved qualitative method of Fornberg–Whitham equation. Feng and Wu [37] reduced this equation to simple ODE and solved it by factorization technique successfully. Jiang and Bi [38] presented the bifurcation method and smooth travelling wave solutions are obtained. Reduced differential transform method has been utilized by Hesam [39] for this type of equations and in convergent power series, the results are obtained. Biazar and Eslami [40] used He’s HPM for the analytical solution of equation (1). Three different types of Fornberg–Whitham equation were investigated by Boutarfa et al. [41] who used the reproducing kernel Hilbert space method. A novel analytical approach for solving fractional Fornberg-Whitham equation was proposed in [42]. An iterative method for finding approximate analytical solutions to the fractional Fornberg-Whitham equation was proposed in [43]. An analytical solution of Fornberg-Whitham type equations in view of the fractional Caputo operator was considered in [44]. Our aim in this work is to apply the MVIA-I for these two types of Fornberg–Whitham type equations.

The remaining sections of the paper are organized on the basis of the following organization. In Section 2, some special and basic functions of fractional Calculus are defined. In Section 3, the fractional iteration algorithm-I is explained. Its convergence analysis is discussed in Section 4, and its implementation is illustrated for nonlinear Fornberg–Whitham type equations in Section 5. In Section 6, the results are discussed, and some utilizations of the proposed scheme are given. Conclusion is explained in the last Section 7.

2. Basic Functions and Properties of Fractional Calculus

This section aims to present some basic definitions and notions, which have significant importance in the FC. First of all, we discuss some basic special functions. Basic definitions of fractional derivatives were introduced thereafter.

2.1. Special Functions

In fractional/fractal calculus, the gamma, Beta, and Mittag–Leffler functions are the fundamental and key tools to understand the origin of its computational challenges.

2.1.1. Gamma Function

The gamma function is defined by the Euler integral of the second kind [45]:This integral is convergent in the right half of the complex plane .

2.1.2. Beta Function

The Beta function is defined by the Euler integral of the first kind [45]:

The relationship between Beta function and Gamma function can be established as

2.1.3. Mittag–Leffler Function

We describe the one-parameter and two-parameter Mittag–Leffler function, respectively, as [46]

2.2. Different Definitions of Fractional Derivatives

We present some main definitions of the fractional derivatives as follows:

Definition 1. The right-modified Riemann–Liouville derivative is defined as [47]

Definition 2. The Caputo fractional derivative is defined as [48]

Definition 3. The ABC derivative is introduced as [49]where is the normalization function.

Definition 4. The Ji-Huan He’s fractional derivative is defined by the rule [50]

Definition 5. The Ji-Huan He’s fractal derivative is defined as [50]

Definition 6. The Grünwald–Letnikov fractional derivative is presented as [45]

3. Fractional Iteration Algorithm-I

In this section, the main idea of the fractional iteration algorithm-I is illustrated by considering a nonlinear differential equation of the generic formwhere and denote linear and nonlinear operators, respectively, whereas is a nonhomogeneous term. For an appropriate given initial condition , series , which approximates the solution of equation (13), can be obtained aswhere is a restricted term, which gives with respect to the variation , and and are an auxiliary parameter and Lagrange multiplier, respectively, which can be optimally determined. The first one is used to accelerate the convergence to the exact solution [5155], while the second one is used to construct the correction function [56].

A suitable value of can be achieved applying on both sides of recurrent relation (14) with respect to , which leads to

For nonlinear problems, the nonlinear terms have to be considered as restricted variations for obtaining the value of Lagrange multiplier, and a correction functional can be easily constructed after determining the identified value of corresponding nonlinear terms. The beneath Lagrange multipliers can be obtained in the following way:

Also, the following general formula for the Lagrange multiplier in the cases is available:

After finding the value of , an iteration formula is constructed by using this value in the corrective function (14) as follows:

The iterative sequence can be obtained starting from a proper initial approximation and using the iterative formula (18). It is convenient to repeat iterations many times to arrive at the given accuracy for the advanced computer technique. An exact solution is obtained as the limiting value

It is worth mentioning that the proposed algorithm may be considered as a nice refinement in existing analytical and numerical methods, where the discretization, transformations, and linearization are not required, and the numerical solution of the given differential equations can be obtained in series form as

We employ this procedure for finding the analytical/numerical solution of nonlinear fractional-order Fornberg–Whitham equations. When , this procedure given in equation (20) becomes the standard variational iteration algorithm-I. Equation (20) has two obvious advantages; one is the limited step, which is needed for better accuracy, while the other is an auxiliary parameter , which ensures the convergence, and a more accurate solution can be gained after a higher iteration process.

4. Convergence Analysis

The convergence of the algorithm proposed in Section 3 for solving nonlinear fractional partial differential equations will be examined in this section. This algorithm can be performed in an effective and trustworthy way and can handle fractional differential equation (13) as well. When the fractional iteration algorithm is implemented for the numerical investigation of the fractional-order Fornberg–Whitham equation, the linear operator is defined as . First, in (21), we define the operator for solving problems of such type:where and and are defined by

In general, for , it can be written as

Accordingly,

The initial iteration can be chosen uninhibitedly, but it needs to fulfil the corresponding initial-boundary conditions. The determination of appropriate initial approximation will give productive and accurate results. The -order truncated series can be used to approximate the solution. The unknown parameter in ensures that the hypothesis is fulfilled by utilizing 2-norm error of the residual function. The error analysis and convergence criteria of VIA-I with an auxiliary parameter are revealed using the following theorems [57, 58].

Theorem 1. The operator defined in (21) maps a Hilbert space to . The solution given in (24) can be given in the following form of series:

It converges if , such that

Lemma 1. Let be a function from a Hilbert space to , and the operator required in (13) be defined as and the Lagrange multiplier be defined optimally by the variation theory, then

Theorem 2. Let the operator needed in 14 be defined as . If we have the series solution (24) defined bythen is an exact solution to the nonlinear partial differential equation (13).

Theorem 3. Let us suppose that the solution , given in (24), converges to the exact solution of the model equation (1). Also, assume that if the approximate solution is the truncated series then the maximum error norm can be assessed as

5. Implementation of the Algorithm and Examples

In this section, to clarify step by step solution procedure of the fractional iteration algorithm-I, the following time-fractional Fornberg–Whitham equation (1) is considered:

The initial condition isand the exact solution taken from [41] is equal to

Numerical solution of equation (13) for the provided initial condition can be achieved by means of

The Lagrange multiplier can be obtained using the variation theory. Multiplying both the sides of the equation (33) by , one obtains

For nonlinear problems, the nonlinear terms are considered to be restricted variations for obtaining the value of Lagrange multiplier; i.e., is a restricted term. This implies and gives the beneath value of the Lagrange multiplier . After using the value of in the recurrence relation (34), the beneath recurrent dependence is obtained:

Starting with a proper initial approximation and using the iterative formula (35), values in other iterations can be obtained. We stop the process at the third iteration. A residual function used before for variational iteration algorithm-II in [59] can be defined here for approximated solution to get optimal value of unknown parameter:And error of norm 2 of the above residual function with respect to the parameter for can be defined as

The value of is approximated using a numerical integration. The minimal value of the function (37) is the value of the auxiliary term . Here, for different values of , we get different values of auxiliary parameters using Maple software. For , , , , , the corresponding values of the auxiliary parameter are 1.21413823669280, 0.461366433098273, , and 1.36918865275675, respectively.

It is convenient to iterate many times to arrive at the given accuracy for the advanced computer technique, and the exact solution is obtained by the limit

The numerical solutions generated using the proposed technique with an optimal value of the auxiliary parameter are reported in Table 1. To show the applicability and efficiency of the proposed technique, the comparisons are carried out with the exact solutions as well as the standard variational iteration method from [60], in terms of absolute errors reported in Table 2 for different values of and .

A full agreement between the results of fractional variational iteration algorithm and exact solution can be observed, which confirms applicability and accuracy of the proposed algorithm. One can observe that the proposed algorithm is very rapid, effective, and accurate, and this is proved by comparing the solutions obtained through the proposed method with the results in [60] as well as with the exact solution. Results are also shown graphically by assigning different values of with under the given condition in each case, which can be seen in Figure 1.

The space and time surface graphs of approximate and exact solutions are observable in Figures 27.

6. Discussion

The nonlinear PDEs have become a hot topic in the field of nonlinear science and have been used in modeling miscellaneous problems in many fields of science and engineering, including plasma physics, fluid mechanics, quantum mechanics, fluid dynamics, image processing and economic system. However, for some complex problems in these fields, the fractional partial differential equation is more accurate than integer partial differential equation. So, finding the numerical solutions of PDEs, as well as of fractional PDEs, has emerged as an extremely important strategy.(1)Many engineering problems and other physical phenomena are modeled by nonlinear PDEs. On the other hand, the proposed iteration algorithms can be applied to various nonlinear PDEs. As a consequence, the proposed iterative scheme can be used in solving various problems.(2)The fractional-order operator plays an important role not only in mathematics but also in other fields such as mechanics, physics, biomedical engineering, and finance. There are many nonlinear problems in real life. Therefore, there is a need to study the fractional PDEs.(3)The numerical approaches are used for deeper understanding to predict the anomalies, which are not possible in the analytical methods because the analytical method can solve only two or three unknown variables, while numerical methods are applicable to much more variables very accurately.(4)Analytical methods, if available, are always the best. However, these are impossible to achieve in some cases. If the analytical approaches are unable to generate the exact solution, numerical methods give the possibility to obtain an approximate numerical solution.(5)The extreme benefit of the proposed algorithm is based on its straightforward applicability and conciseness. Furthermore, it can handle all types of nonlinear fractional PDEs.

7. Conclusion

In this work, a novel approach named as fractional iteration algorithm-I is introduced for various types of nonlinear PDEs as well as fractional PDEs and has been discussed in detail, including figures and tabulated numerical results. The proposed method is able to implement without the use of any transformation, linearization, discretization, or restrictive assumptions and thus is particularly perfect with the flexible and expanded nature of the physical problems. The coupling of Caputo fractional derivative and variational iteration algorithm-I to tackle both time and space derivatives accurately makes the solution process easiest, and the evaluation of fractional term becomes simpler, and FPDEs can be solved in a tremendous way. The proposed technique is able to solve all kinds of linear and nonlinear problems of physical nature arising in applied sciences and engineering.

Data Availability

Data will be provided on request to the first author.

Conflicts of Interest

The authors declare that there are no conflicts of interest associated with this publication.

Acknowledgments

Predrag Stanimirović acknowledges support from the Ministry of Education, Science and Technological Development, Republic of Serbia, under Grant no. 174013. The research was supported by the National Natural Science Foundation of China (Grant nos. 11971142, 11871202, 61673169, 11701176, 11626101, and 11601485).