New Numerical Approach of Solving Highly Nonlinear Fractional Partial Differential Equations via Fractional Novel Analytical Method

: In this work, the fractional novel analytic method (FNAM) is successfully implemented on some well-known, strongly nonlinear fractional partial differential equations (NFPDEs), and the results show the approach’s efﬁciency. The main purpose is to show the method’s strength on FPDEs by minimizing the calculation effort. The novel numerical approach has shown to be the simplest technique for obtaining the numerical solution to any form of the fractional partial differential equation (FPDE).


Introduction
In the 18th century, the work of Euler, D'Alembert, Laplace, and Lagrange heavily relied on the use of partial differential equations (PDEs) as a key tool for examining physical models and describing the mechanics of continua. Similarly, PDEs were a crucial tool for scientists throughout the middle of the nineteenth century, particularly in the work of Riemann. They are now essential for characterizing dynamical systems in the physical and mathematics sciences [1]. PDEs may be used to simulate many circumstances and phenomena in Physics, Fluid Mechanics, Heat and Wave Analysis, Quantum Mechanics, Engineering, Chemistry, and Electrodynamics [2,3]. In general, the dynamics of any physical system can be addressed by either linear or nonlinear PDEs [4,5].
Since there are not really any strategies that work for all equation types and each equation frequently needs to be addressed as a separate issue, it is challenging to solve nonlinear EDPs in general. With analytic solutions to PDEs, we may test and comprehend the numerical methods as well as the physical processes under consideration [6]. A robust technique is illustrated, and the Finite Approximation methodology is detailed in depth in [7]. The study in [8] discusses the Differential Transform approach to the physical application of electrical circuits. In [9], the Fishers Reaction Diffusion Equation is solved by utilizing the Least Squares Finite Element Approximation approach, whereas the authors of [10] applied the Haar Wavelet Method to the Fishers Reaction-Diffusion Equation.
The research in [11] contains a detailed description of how to solve Eigen Value Problems and Boundary Value Problems using Fourier Spectral Methods, Galerkin Methods, Chebyshev technique, Pseudospectral Methods, the Tau Method, Domain Decomposition Methods, and methods for unbounded intervals. The Homotopy Perturbation approach is used on both homogeneous and inhomogeneous PDE systems [12].
For semiconductor devices, the Poisson-Boltzmann Equation is discussed by Homotopy Analysis Method in [13]. In [14], the applicability of Adomain Decomposition Method (ADM) is demonstrated by solving non-integer PDEs in infinite domains.
PDEs of Fractional Order have also piqued the interest of mathematicians and physicists [15]. For example, in [16], a numerical approach known as PIA was used for several highly nonlinear PDEs, and quick convergence of these equations to their precise solutions was demonstrated. The interpolating element-free Galerkin method was applied to 2D generalized Benjamin-Bona-Mahony-Burgers (BBMB) and regularized Long-Wave Equations on non-rectangular domains in [17]. The Meshless approach using radial basis functions was used in [18] to determine a numerical solution to the nonlinear, high-dimensional, generalized BBMB problem. The test problems for different geometries of the two-dimensional and three-dimensional cases of the generalized Benjamin-Bona-Mahony-Burgers Equation (BBMBE) have been solved in this study.
The VIM and HPM are used to solve the Fornberg-Whitham Equation (FWE) in [19]. On Burgers Equation (BE), [20] implements a multi-symplectic box approach with two kinds of box schemes. Three approaches, namely ADM, VIM, and HPM, are used for the Fitzhugh-Nagumo problem in [21], and numerical results are produced. The material presented in [16][17][18][19][20][21] contains thorough information about the PDEs described and piques the reader's curiosity. In [22,23], the authors have offered a full introduction to different types of PDEs and their solutions; also, several numerical techniques with different examples are offered in these two books. Several transformations are employed in the aforementioned numerical procedures.
In this article, some highly nonlinear fractional partial differential equations (NF-PDEs) will be solved by the Fractional Novel Analytic Method (FNAM), introduced by [24][25][26]. In [26], this method was introduced as Novel Analytical Method. This approach is simple to use and saves time over other approaches. This paper is organized as follows: Section 3 discusses the mathematical formulation of this strategy, Section 4 contains the convergence of this method, Section 5 illustrates numerical examples, and the conclusion is presented in Section 6.

Preliminaries of Fractional Calculus
In this section, we introduce the fundamental concepts and features of fractional calculus [27][28][29].
and it is said to be in the space C ℵ υ , ℵ ∈ N. Definition 2. The Riemann-Liouville (RL) Fractional Integral (FI) operator of order α, of a function G ∈ C υ , υ > −1, is defined as: Numerous scholars have recently researched various RL-FI inequalities; for more details, see [27][28][29]. The features of operator J α that we need in this study are: Definition 3. The Fractional Derivative (FD) of G(τ) in the Caputo sense [30] is defined as: In Caputo FD, an ordinary derivative is estimated, followed by a FI, to attain the desired order of FD. The RL-FI operator is a linear operation, defined as: where {C } h are constants. The FDs are considered in the Caputo sense in this study.

Highly Nonlinear Partial Differential Equations of Fractional Order via Fractional Novel Analytical Method
In this section, we will go through the fundamental ideas behind implementing an FNAM for the NFPDE. Let us consider the following general Fractional Order PDE: with initial condition Taking the Fractional Integral (FI) for both sides of Equation (1) from 0 to ξ, we get , · · · . Then, again taking the FI from 0 to ξ, on both sides of Equation (3). We obtain, Thus, Substituting Equation (5) with Equation (4), we obtain where such that the highest derivative of is k. The endorsement of Equation (6) is to extend FTS for about ξ = 0. It means that, thus we simply obtain our wanted numerical solution.

Novel Analytical Method Convergence Analysis
Let us examine the following PDE where a nonlinear operator is G. Using the approach outlined, the achieved solution is identical to the following sequence Theorem 1. Let G be an operator from H → H (Hilbert space) and be the exact solution of Equation (7). The estimated solution Proof of Theorem 1. We seek to demonstrate that {S n } ∞ n=0 is a converged Cauchy Sequence, Now for n, m ∈ N, n ≥ m, we obtain Hence lim is converged to when 0 ≤ σ i < 1, i = 0, 1, 2, 3, · · · for further details, see [22,23].

Numerical Application for Fractional Novel Analytical Method for Highly Nonlinear Partial Differential Equations
Four problems are addressed in this section utilizing the suggested Fractional Novel Analytical Method.

Example 1. The Gardner Equation,
which is used to describe internal solitary waves in shallow water, was generated by combining the KdV and modified KdV Equations [32]. Consider the Nonlinear Time Fractional Gardner Equation [33], where 0 < α ≤ 1. By carefully following the steps involved in the Fractional Novel Analytical approach, we obtain, When α = 1, the exact solution of this problem is where β, λ > 0. Comparisons between exact and numerical solutions are plotted in Figure 1

Conclusions
The numerical approach provided here has been utilized to solve highly nonlinear fractional partial differential equations. The purposed method is very efficient at finding and comparing the results of nonlinear FPDEs with the exact solution. Based on the numerical solutions obtained, we can affirm that there is an excellent agreement with the existing solutions and demonstrate that this technique can be utilized to tackle the proposed problems effectively. By adding appropriate terms to the solution's series, the error can be reduced. The solution obtained at each fractional order is found to be converging to the integer order. We have compared the provided numerical solutions with exact solutions to demonstrate the validity and vast potential of the suggested numerical approach. The proposed technique is more reliable, efficient, and easy to use; it is also less computationally intense to find the solution to the nonlinear partial differential equation of the non-integer order. This method could be applied to a system of fractional differential equations, fractional oscillators, fractional integro-differential equations, etc.