Analytical Solutions of Fractional-Order Heat and Wave Equations by the Natural Transform Decomposition Method.

In the present article, fractional-order heat and wave equations are solved by using the natural transform decomposition method. The series form solutions are obtained for fractional-order heat and wave equations, using the proposed method. Some numerical examples are presented to understand the procedure of natural transform decomposition method. The natural transform decomposition method procedure has shown that less volume of calculations and a high rate of convergence can be easily applied to other nonlinear problems. Therefore, the natural transform decomposition method is considered to be one of the best analytical techniques, in order to solve fractional-order linear and nonlinear Partial deferential equations, particularly fractional-order heat and wave equation.

1. Introduction

The idea of entropy and fractional calculus are attractive, further prevalent for investigating the dynamics of complex schemes. In recent years, fractional calculus (FC) has been progressively applied in various fields of science. Natural development identified with viscoelasticity, models of porous electrodes, thermal stresses, electromagnetism, propagation of energy in dissipative systems, relaxation vibrations and thermoelasticity are effectively portrayed by fractional differential equations (FDE’s) [1]. The knowledge of entropy was presented in the field of thermodynamics by Clausius (1862) and Boltzmann (1896) and was further applied by Shannon (1948) and Jaynes (1957) in information theory. Newly, more universal entropy measures have been suggested for applications in numerous varieties of complex systems outstanding for the relaxation of the additives axiom [2]. The concept of entropy for analyzing the dynamics of multi-particle systems with an integer and fractional order behavior. The entropy production rate for the fractional diffusion procedure was considered in [3,4]. In [5], it has been shown that the total spectral entropy can be used as a measure of the data comfortable in a fractional order model of anomalous diffusion. Entropies based on fractional calculus [6], integer and fractional dynamical systems can be solved by entropy analysis [7], nonlinear partial differential equations and third-order dispersive [8,9] in entropy and convexity. Bifurcation and recurrent analysis of memristive circuits [10], density analysis of multi-wing and multi-scroll chaotic systems [11], and numerical integration in a simulation of conservative chaotic systems [12]. Fractional derivative advection–diffusion in two-dimensional semi-conductor systems and the dynamics of a national soccer league [13]. The exact solution to differential equations (DEs) of fractional order with mixed partial derivatives [14] and space-fractional diffusion equation and Tsallis relative entropy [15].

Joseph Fourier first developed the heat equation in 1822. The heat equation is connected with the study of Brownian motion. This motion was solved by Robert Brown. The irregular movement of particles suspended in a liquid (gas or fluid) come about because of their impact with the quick moving atoms in the liquid. Heat is the dynamic energy of particles that are being exchanged. When emotional vitality is exchanged from increased surrendering to a cooler scheme, faster moving atoms in the environment crash into the dividers of the scheme that exchanges a portion of this energy to the atoms of the system and makes them move faster [16,17,18]. The wave equation is a significant second order linear partial differential equations (PDE’s) for the description of waves as they occur in traditional physics such as fluid dynamics, mechanical waves and the electromagnetic wave equation is an important PDE that arises in fields like acoustics, electromagnetics and fluid dynamics. They are light, sound, gravity and even matter (in relativistic quantum mechanics, the Klein–Gordon equation). The one and three-dimensional wave equation was discovered by Alembert and Euler. The solutions of heat and wave equations have attracted the attention of various authors in mathematics, such as the optimal homotopy asymptotic method (OHAM) [19], the modified Adomian decomposition method (MADM) [20], the variational iteration method [21], the differential transform method (DTM) [22], the homotopy perturbation method (HPM), [23], Bernstein polynomials with the operational matrix [24], Elzaki transform and the projected differential transform method for nonlinear wave equations [25], the variational iteration method with the help of the Yang–Laplace transform fractional heat equation [26], and the Aboodh decomposition method [27].

In the present work, we are applying the natural transform decomposition method (NTDM), to solve the following types of fractional partial differential equations (FPDE’s). with initial condition

Natural transform and Adomian decomposition methods are two powerful methods that have been used to develop the natural transform decomposition method. Many physical phenomena which are modeled by PDE and FPDEs are solved by using NTDM, such as the analytical solution of a couple of systems of nonlinear PDE’s is suggested in [28], the solution nonlinear ODE’s are successfully presented in [29], nonlinear PDEs [30], fractional unsteady flow of a polytropic gas model [31], fractional telegraph equations [32], fractional Fokker–Plank equation and Schrödinger equation [33]. The accuracy of the proposed method is compared with the solutions obtained by HPM and Modified homotopy perturbation method (MHPM). The comparisons has shown that the proposed has a higher rate of convergence than HPM and MHPM. The rest of the article is structured as: in Section 2, we recall several basic properties and definitions from natural transform and fractional calculus. In Section 3, we present the idea of the natural transform decomposition method. In Section 4, we explain many problems with maintaining the accuracy and efficiency of the proposed method, while the last section is devoted to conclusions.

The natural transform of where s and u are the transform variables.

The inverse natural transform of a function is defined by where s and u are the natural transform variables and p is a real constant and the integral is taken along

Natural Transform of nth Derivative

If

If where

The Riemann–Liouville fractional integral where

In this study, Caputo et al. suggested a revised fractional derivative operator in order to overcome inconsistency measured in the Riemann–Liouville derivative. The above mathematical statement described a Caputo fractional derivative operator of initial and boundary conditions for fractional as well as integer order derivatives [

The Caputo operator of order γ for a fractional derivative is given by the following mathematical expression for

3. Idea of the Fractional Natural Transform Decomposition Method

In this section, the natural transform decomposition method to find the general solution fractional-order equations: where the Caputo Operator , where L and N are linear and nonlinear functions, and q is the source function.

The initial condition is

Applying the natural transform to Equation (1), we have and using the differentiation property of natural transform, we get

The NTDM solution is represented by the following infinite series: and the nonlinear terms (if any) in the problem are defined by the infinite series of Adomian polynomials, substitution Equation (5) and Equation (6) in Equation (4), we get

Applying the linearity of the natural transform,

Generally, we can write

Applying the inverse natural transform, in Equation (9),

Consider the one-dimensional fractional heat equation [ with initial condition

Taking the natural transform of Equation (

Applying inverse natural transform,

Using the ADM procedure, we get for

The subsequent terms are

The NTDM solution for Example 1 is when

This result is calculated to get the exact solution in a closed form:

4. Results

Figure 1a,b shows the behavior of obtained solution by the proposed method for different values of , , , and , then Figure 1c,d are error plots for . Figure 2 combine error plots the range of x .

Figure 1

(a) Comparison with numerical solution of by NTDM of Example 1, for different values of and (b) for ; (c) and (d) are error plots of Example 1.

Figure 2

Error plot of Example 1 for different values of .

Consider the two-dimensional fractional heat equation [ with initial condition

Taking the natural transform of Equation (

Applying inverse natural transform,

Using the ADM procedure, we get for

The subsequent terms are

The NTDM solution for Example 2 is

This result is calculated to get the exact solution in a closed form:

Figure 3a shows the behavior of obtained solution by the proposed method for different values of , , , and ; Figure 3b error plot for the range of x and y .

Figure 3

(a) Comparison with numerical solution of by NTDM of Example 2, for different values of ; (b) error plot for .

Consider the three-dimensional fractional heat equation [ with initial condition

Taking the natural transform of Equation (

Applying inverse natural transform

Using the ADM procedure, we get for

The subsequent terms are

The NTDM solution for Example 3 is when

This result is calculated to get the exact solution in a closed form:

Figure 4a,b shows the behavior of obtained solution by the proposed method for different values of , , , and ; Figure 4c error plot for the range of x and y .

Figure 4

(a) Comparison with numerical solution of by NTDM of Example 3, for different values of and (b) for ; (c) error plot of Example 3 for .

Consider the one-dimensional fractional heat equation [ with initial condition

Taking natural transform of Equation (

Applying inverse natural transform,

Using the ADM procedure, we get for

The subsequent terms are

The NTDM solution for Example 4 is when

This result is calculated to get the exact solution in a closed form:

Figure 5a show the behavior of obtained solution by the proposed method for different values of , , , and ; Figure 5b error plot for the range of x .

Figure 5

(a) Comparison with numerical solution of by NTDM of Example 4, for different values of ; (b) error plot for .

Consider the two-dimensional fractional wave equation [ with initial condition

Taking natural transform of Equation (

Applying inverse natural transform

Using the ADM procedure, we get for

The subsequent terms are

The NTDM solution for Example 5 is

This result is calculated to get the exact solution in a closed form:

Figure 6a shows the behavior of obtained solution by the proposed method for different values of , , , and ; Figure 6b error plot for the range of x and y .

Figure 6

(a) Comparison with numerical solution of by NTDM of Example 5, for different values of ; (b) error plot for .

Consider the three-dimensional fractional wave equation [

with initial condition

Taking natural transform of Equation (

Applying inverse natural transform,

Using the ADM procedure, we get for

The subsequent terms are

The NTDM solution for Example 6 is

This result is calculated to get the exact solution in a closed form:

Figure 7a,b shows the behavior of obtained solution by the proposed method for different values of , , , and ; Figure 7c error plot for the range of x and y .

Figure 7

(a) Comparison with numerical solution of by NTDM of Example 6, for different values of and (b) for ; (c) error plot of Example 6 for .

Consider the two-dimensional fractional-order nonlinear wave equation [ with initial condition

Taking natural transform of Equation (

Applying inverse natural transform,

Using the ADM procedure, we get where A and B are nonlinear terms. The few nonlinear terms are as follows: and so on for

The subsequent terms are

The NTDM solution for Example 7 is

This result is calculated to get the exact solution in a closed form:

5. Conclusions

In this paper, the analytical solutions of fractional-order heat and wave equations are determined, using NTDM. The NTDM solutions are obtained at fractional and integer orders for all problems. The results revealed the highest agreement with the exact solutions for the problems. The NTDM solutions for some numerical examples have shown the validity of the proposed method. It is also investigated that the fractional order solutions are convergent to the exact solution for the problems as fractional order approaches to integer order. The implementation of NTDM to illustrative examples have also confirmed that the fractional order mathematical model can be the best representation of any experimental data as compared to integer order model. In the future, NTDM can be used to find the analytical solution of other nonlinear FPDEs, which are frequently used in science and engineering. NTDM solutions for fractional order problems will prove better understanding of the real world problems represented by FPDEs.