Abstract
In the current work, a combination between a new integral transform and the homotopy perturbation method is presented. This combination allows to obtain analytic and numerical solutions for linear and nonlinear systems of partial differential equations.
1. Introduction
We know that the HPM, proposed first by He [1], for solving differential [2, 3] and integral equations [4], linear and nonlinear, has been the subject of extensive analytical and numerical studies.
The HPM is applied to singular nonlinear differential equations [5], nonlinear wave equations [7], nonlinear oscillators [6], bifurcation of delay-differential equations [8], boundary value problems [9], initial value problems [10], and nonlinear coupled equations [5, 22]. Furthermore the HPM yields every rapid convergence of the solution series in most cases.
On the other hand, the integral transformations played an essential role in many fields of science [11, 12], especially, engineering mathematics [13], mathematical physics [14], optics [15], image processing [16] and, few others because they have been successfully used in solving many problems in those fields. Many of these transforms have been used and applied on theory and applications, such as Sumudu [17, 18], Laplace [19, 20], Fourier [12], Elzaki et al. [21] and new integral transform [23]. Among these the most widely used is Laplace transform. Here, new integral transform is proposed to avoid the complexity of previous transforms [12, 17, 18].
In general, the nonlinear partial differential equations (NPDEs) have modeled nonlinear complex phenomena in various scientific fields [24ā35]. The investigation of analytical, approximate, and exact solutions of NPDEs will help better understand the complex phynomena.
Our method, which is a coupling of the new integral transform and homotopy perturbation technniques, deforms continuously to a simple problem which is easily solved. Also this presentation has proposed a new method for solving NPDEs.
This article is organized as follows: In Section 2, we introduce some basic definitions and proberities for the new integral transform. In Section 3, we discuss the method used in this work. Some applications are given in Section 4 to show the accuracy and advantage of the proposed method. Finally, numerical results are discussed in Section 5.
2. Basic Definition of the New Integral Transform (NT)
In this section, we mention the following basic definitions and theorems of the new transform used in the present paper.
2.1. Definition of the New Transform
The transform of a function is defined by
Theorem 1. (Sufficient condition). If a function is piecewise continuous on and of exponential order so, then the transform of exists for .
Theorem 2. (Linear combination). If transforms and of in Equation (1) the functions are well defined and are constants, then
Theorem 3. ( Derivatives). If the functions are well defined, then
3. Analysis of the Method
To illustrate the modification algorithm of the NTHPM, we consider the following nonlinear partial differential equation with time derivatives of any order
where, is linear differential operator , represents the general nonlinear differential operator and is the source term, subject to the initial conditions
In view of the homotopy technique, we can constract the following homotopy
where the homotopy parameter always changes from zero to unity. The changing process of is called deformation. When , Equation (6) becomes
and when , Equation (6) turns out to the original Equation (4). Since is a function of only, Equation (6) can be rewritten to be in the following form
According to the homotopy technique, the basic assumption is that the solution of Equation (8) can be written as a power series in as
where are unknown functions to be determined. Now, taking in mind the initial conditions (2), the NT for Equation (8) gives
again taking the inverse of the NT for Equation (10), we obtain
Substituting from Equation (9) into Equation (11), yields
Equating the identical powers of , therefore, after doing some calculations for the NT and the inverse of NT we get the unknown functions After substituting into Equation (10) with , we get the solution of the problem (1)ā(2).
4. Applications on NTHPM
Our method will be illustrated through examples in one-dimension for linear and nonlinear coupled systems of partial differential equations.
Example 4.1. Consider the one-dimensional linear system
subjected to the initial conditions
Assume that the solutions of Equations (13) and (14) can be written as a power series as follows
substituting from Equation (17) into Equation (12) for and and using the initial conditions (11), yields
by a similar way, substituting Equation (18) into Equation (12), for and and using the initial conditions (12), we obtain
On putting the coeffiecients to the power of equal to zero in Equations (19), (20), we can obtain a series of linear equations, which are easy to solve by using Mathematica software to give
and so on. Proceeding as before the rest of components were obtained, and then the two functions and in the closed form are readily found to be
Example 4.2. We consider the homogenuous form of coupled Burgers equations
with the initial conditions
As illustrated in Example 4.1, substituting from Equations (17) and (18) into Equation (12) but in this case for and using the initial conditions (21) and (22) respectively, we get
Putting the coefficients of the power of equal to zero in Equations (37) and (38), we obtain
again substituting the functions into the Equations (17) and (18), we obtain directly
Furthermore, in closed form Equation (48) takes the form
Example 4.3. Next, we consider the nonlinear Drinfeld Sokolov system in the following form
with the initial conditions
By applying the same steps used in Examples (1) and (2), we can easily get
Substituting into Equation (17) and (18), we obtain
5. Numerical Results and Discussion
The numerical behavior of the error between the exact solution and the solution obtained by NTHPM is shown in Figures 1ā3. The numerical results are obtained by using fourth order perturbation only from the series formulas (13), (14) with . From these figures, we achieved a very good approximation for the solution of our systems at the small values of time , but at the large values of the time, the error can be reduced by adding new terms from the iteration formulas.
(a)
(b)
(a)
(b)
6. Conclusion
In this work, we have proposed a modification to the homotopy perturbation technique by combining it with a new integral transform. The aim of this approach is to obtain exact solutions of linear as well as nonlinear coupled systems. The efficiency and accuracy of the present scheme are validated through some examples. The results show that NITHPM is a powerful and good technique for obtaining exact solutions of many systems of linear as well as nonlinear differential equations. The computations associated in this work are performed by using the Mathematica software.
Conflicts of Interest
The authors declare that they have no conflicts of interest.