Abstract

In this paper, we are presenting our work where the noninteger order partial differential equation is studied analytically and numerically using the noninteger power series technique, proposed to solve a noninteger differential equation. We are familiar with a coupled system of the nonlinear partial differential equation (NLPDE). Noninteger derivatives are considered in the Caputo operator. The fractional-order power series technique for finding the nonlinear fractional-order partial differential equation is found to be relatively simple in implementation with an application of the direct power series method. We obtained the solution of nonlinear dispersive equations which are used in electromagnetic and optics signal transformation. The proposed approach of using the noninteger power series technique appears to have a good chance of lowering the computational cost of solving such problems significantly. How to paradigm an initial representation plays an important role in the subsequent process, and a few examples are provided to clarify the initial solution collection.

1. Introduction

Mathematical equations containing two or more independent variables are called partial differential equations (PDEs). The highest derivative in the PDEs is called the order of the equation. The specific solution of the PDEs is a function that solves the equation, which can be verified by substituting the solution into the equation and obtaining an identity. The solution containing all the particular solutions of the equation in question is called the general solution. The term used to express a specific solution for second-order and higher-order nonlinear partial differential equations is called an exact solution. Problems involving the function of many variables, such as heat or sound propagation, static electricity, electrodynamics, and fluid flow, are solved by PDE [1, 2].

Generally, the solution with unknown coefficients is assumed to be a power series, and then the solution is substituted into the differential equation to find the recurrence relationship of the coefficients. New applications in fluid dynamics, viscoelasticity, mathematical biology, electrochemistry, and physics have led to the latest development in fractional DE. Using data from experiments, fractional DEs have recently proved to be useful tools for simulating many physical phenomena and fractional PDEs have been suggested for filtering flows in porous media. The examples that have been studied and solved include the space-time fractional diffusion wave equation, fractional advection scattering equation, fractional telegraph equation, fractional KdV equation, and linear nonhomogeneous fractional partial differential equations [2,3].

Nowadays, the uses of nonlinear partial (or ordinary) differential equations in the practice of multidisciplinary research are considered to be an important technique for solving real-life problems that have social and public-health implications. Researchers in different scientific areas like biology, physiology, cell biology, chemistry, chemical physics, engineering, and systematic natural sciences are including applied mathematical methods in their practices. These applications are not only useful for technological advances in computer science to study cancers, hereditary diseases, tumor formation, kinematic processes, and DNA communication processes but also in mathematical controls to supply a successful cancer therapy. Indeed, there are many biomedical engineering applications for studying biological problems involving highly nonlinear and sophisticated collaborative processes, where the most promising theory available to conduct such studies is complex dynamical systems theory. The description of this theory almost entirely depends on nonlinear partial differential equations (NLPDE). Power series solution (PSS) technique has been restricted to unravel the linear differential equations, i.e., ordinary differential equations (ODEs) [1, 2] and partial differential equations (PDEs) [4, 5].

Linear partial differential equations have traditionally been overcome using the variable separation method because it creates an ODE system that is easier to decipher with PSSM. Examples of them are the spherical harmonics used and the Legendre polynomials in the Bessel equation in cylindrical coordinates or the Laplace equation in spherical coordinates [4, 5]. It is understood that in NLPDE, as we remember, variable separation technology cannot be used. The way to unlock NLPDE is to use approximate analysis techniques to get the answer, that is, semianalytical or nondigital, direct or indirect. Directly, there are methods such as loose operator formalism [6] or inverse dispersion transformation [7]. In a direct way, for instance, the PSS can be used in an asymptotic calculation of the Hirota technique linking a bilinear operator method [8] and the homotopy analysis method (HAM) [9, 10] and the Adomian decomposition technique [11]. The latest approach requires a sequence of expansions with a nonsmall perturbation parameter estimate to change the convergence. This approach is distinct from the standard perturbation theory. Methods, also simpler, to estimate the representation in NLPDE are the Taylor polynomial approximation (TPA) technique [12, 13] and therefore for the power series solution method (PSSM). In each technique, a semianalytical representation is attained to implement the PSS technique. However, no use has been made of the PSSM to decipher the nonlinear ODE [1417] or NLPDE [1820]. Because NLPDEs do not have exact solutions in general, it is hard to define the most efficient method to be used; indeed, the only way to establish the efficiencies of the relevant methods involves employing experimental methods.

Further information on how to solve NLPDEs by using a point contact, Bäcklund, hodograph, Legendre, or Euler transformations or by applying Lie algebras and groups can be found in reference [21]. There is an immense work on noninteger calculus, and this has been grown exponentially in recent years, with noninteger ODEs becoming commonly used as mathematical prototypes [22, 23]. Recent studies consist of the solution of fractal media [24], viscoelastic materials [23], economics and finance [22, 25], compartment models [26], in porous media [27], epidemiology [28], nerve cell signaling [29], and anomalous diffusion system [30]. From a modeling perspective, the interpretation of the representation as an algebraic countenance involving functions of structure strictures is commonly pursued. Important progress has been made concerning linear noninteger ODEs by expanding proven techniques of elucidation for linear nonfractional ODEs. The Laplace transform technique can be extended to noninteger linear ODEs, with constant coefficients (see [3133]).

The disadvantage of this method depends on the complexity of the Laplace transform of such functions, regardless of whether the equations are homogeneous or not, or by inverting the Laplace transform based on a calibrated computable function. In the noninteger linear ODE, with variable coefficients, the Laplace transform method has an additional limitation, that is, the Laplace transform of the product of the function is only available in special cases. In nonfractional regular linear ODEs, with different coefficients, the series expansion method is considered to help clarify the problem. Examples of power series can be seen in [3448]. The description of this theory is also related to the use of the residual power series method [49] to solve space-time fractional PDEs. The comparison of Aboodh transformation and differential transformation method (DTM) numerically in result solution obtained is compared the solution by DTM is rapid convergent [50]. Several papers have been dedicated to the application of series expansion approaches to linear noninteger ODEs which basically is the subject of this article.

2. Fundamental Concepts

Definition 1. The Caputo noninteger derivative operator of order with respect to “t” is defined in the subsequent arrangement.where .
In a comparable manner of nonfractional-order differentiation, the Caputo noninteger derivative operator is a linear operation.where and are constants. We have for the Caputo’s derivative, if is constant.We practice the maximum function to signify the least integer larger than or equivalent to , and . Evoke that for , the Caputo differential operator corresponds with the normal differential operator of nonfractional order.

Definition 2. For the variable and coefficients , if , the fractional power series (FPS) about the point is defined aswhere .

Theorem 1. Let the radius of convergence for the function with fractional power series (FPS) representation , , be greater than zero (i.e., ).

Then, for and , the following expression holds true:

Proof. From the linearity concept of Caputo derivative and the idea of power series derivative, we haveFrom the power rule of the Caputo derivative, we haveFrom equations (3) and (5), we arrive atwhere .

2.1. Application

In order to validate the high degree of efficiency and precision of the projected FPS approach for unraveling fractional-order systems, numerical forms and instances are pragmatic. The reader can discover a sketch and applications for this technique in [42]. Computations were accomplished by using MATLAB. Examples are taken from [32].

Example 1. Consider the succeeding noninteger order structure:subject to the initial conditionsThe closed-form result of this arrangement when isSolution:So using Theorem 1, equations (12) and (13) becomeSimilarly,Apply (12)–(17) on the first and second equation of system (9), respectively, system (9) becomesComparing the different powers of to both sides of system (18), for ,But , , , and , so (19) becomesFor ,But , , , , so (21) becomesFor ,But , , , ,, , so (26) becomesorFor ,In a similar way, we can compare the values for , ,
So,If , then equations (27) and (28) give

Example 2. Consider the following fractional-order system:subject to the initial conditionsThe closed-form result of this system when isSolution:
LetApply (12)–(17) on the first and second equations of the system (30), respectively. So, system (30) becomesComparing the different powers of to both sides of the system (35), for ,After substituting the values , , and , we obtainFor ,After substituting the values , , , , , and , we obtainFor ,After substituting the values , , , , , , , and , we obtainFor ,After inserting the values , , , , , , , , and , we obtainIn a similar way, we can compare the values for , ,
So, by inserting the values of , , , , , , , , , and in (33 and 34), we haveIf , then equations (44) and (45) give

Example 3. Consider the following fractional-order system:subject to the initial conditionsThe exact solution of this system when isSolution:
LetApply (12)–(17) on the first and second equations of the system (47), respectively, system (47) becomesComparing the different powers of to both sides of the system (52), for ,After substituting the values , , and , we obtainFor ,After substituting the values , , , , , and , we obtainFor ,After substituting the values , , , , , , , and , we obtainFor ,After inserting the values , , , , , , , , and , we obtainIn a similar way, we can compare the values for , , …...
So, by inserting the values of , , , , , , , , , and in (51) and (52), we haveIf , then equations (61) and (62) give

3. Results and Discussion

Tables 1and 2 show the values of and for Example 1 with different values of , and . Afterward, the graphical simulations are drawn with the help of MATLAB. The graphs are shown in 2D and 3D, and these can be seen in Figures 1 and 2. Tables 3 and 4 elaborate the values of and for Example 2 and graphical imitations can be seen in Figures 3 and 4 for different values of , and . Besides these, Tables 5 and 6 represent the values of , , and and simulations can be viewed in Figures 5 and 6.

Fractional physical equations can be solved successfully by using analytical and approximate solutions method called fractional power series method (FPSM). The Caputo operator is presented as fractional derivatives. As compared to other methods to solve nonlinear equations, the power series method is employed efficiently to obtain the solution. Results thus obtained show that FPSM can be applied to solve the system of nonlinear partial differential equations (PDEs) with accuracy and effectiveness.

4. Conclusion

Numerical simulation results prove the validity and reliability of the power series method. Therefore, it is more convenient to describe that PSM is very powerful in constructing fractional power series solutions for fractional PDEs in scientific problems of any order and time of space. As per the work done, we had offered the fractional-order power series method for finding the nonlinear fractional-order partial differential equation with comparable easy implementation. In the tables, reasonable values are presented for fractional and integer cases. 2D and 3D graphical presentations are given for fractional-order for , , , and 1. The proposed method is quite promising in reducing the computational cost of solving such problems to a great extent.

Data Availability

All the data are available in the article and cited wherever required.

Conflicts of Interest

The authors declare that they have no conflicts of interest.