In this paper, new homotopy perturbation method (NHPM) (Biazar et al., 2007) is applied to obtain solutions of system of Burgers equations. In this method, the solution is considered as an infinite series expansion where it converges rapidly to the exact solution. In order to show the ability and reliability of the method some examples are provided. The results reveal that the method is very effective and simple. The modified method accelerates the rapid convergence of the series solution and reduces the size of work.
System of Burger equation; New homotopy perturbation method
The investigation of the numerical and exact solutions for nonlinear partial differential equations (NLPDEs) plays an important role in the study of nonlinear physical phenomena. Nonlinear wave phenomena appear in various scientific and engineering fields, such as fluid mechanics, plasma physics, optical fibers, biology, solid-state physics, chemical kinematics, chemical physics and geochemistry. The dispersion, dissipation, diffusion, reaction and convection of nonlinear wave phenomena are very important in nonlinear wave equations. In the past several decades, there have been significant improvements in the study of exact solutions. In the recent years, many authors mainly had paid attention to obtain solutions of NLPDEs by using various methods.
The purpose of this paper was to apply new homotopy perturbation method (NHPM) [10] to coupled Burgers equations.
General form of system of Burgers equations can be considered as the following forms
|
(1) |
with initial conditions
|
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \alpha \mbox{,}\beta \mbox{,}\eta }
are constants, and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f\mbox{,}g} are inhomogeneous terms.
The Burgers model of turbulence is a very important fluid dynamic model and the study of this model and the theory of shock waves have been considered by many authors, both to obtain a conceptual understanding of a class of physical flows and for testing various numerical methods. The distinctive feature of Eq. (1) is that it is the simplest mathematical formulation of the competition between nonlinear advection and viscous diffusion. Homotopy perturbation method (HPM) was proposed by He [1] and [2] which is, in fact, a coupling of the traditional perturbation method and Homotopy in topology. This method has been applied to solve linear and nonlinear equations of heat transfer [3], sixth-order boundary value problems [4], quadratic Riccati differential equation of fractional order [5], conservative truly nonlinear oscillators [6], and linear and nonlinear weakly singular VIE of second kind [7], and many others [8], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21], [22], [23], [24], [25], [26], [27], [28] and [29].
This paper is arranged as follows. In Section 2, the new modification of HPM for solving systems of Burgers equations is presented. The efficiency of this new method is verified by the numerical results for two examples in Section 3. Conclusion will appear in Section 4.
For solving system (1) by NHPM we construct the following homotopies:
|
(2) |
or
|
(3) |
Applying the inverse operator, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle L^{-1}={\int }_{t_0}^t(.)dt}
to both sides of Eq. (3), we obtain
|
(4) |
where
|
Let’s present the solution of the system (4) as the following
|
(5) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U_i\mbox{,}V_i\mbox{,}\quad i=1\ldots n} , are functions which should be determined.
Suppose that the initial approximation of Eq. (1) is in the following form
|
(6) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle a_{i\mbox{,}j}(x_1\mbox{,}x_2\mbox{,}\ldots \mbox{,}x_{n-1})\mbox{,}i=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 1\ldots n\mbox{,}j=0\ldots n , are unknown coefficients and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle P_0(t)\mbox{,}P_1(t)\mbox{,}P_2(t)\mbox{,}\ldots }
are specific functions.
Substituting (5) and (6) into (4) and equating the coefficients of p with the same powers lead to
|
(7) |
Now if we solve these equations in such a way that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U_{i\mbox{,}1}(x_1x_2\mbox{,}\ldots \mbox{,}x_{n-1}\mbox{,}t)=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 0 , then Eq. (7) results in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U_{i\mbox{,}1}(x_1x_2\mbox{,}\ldots \mbox{,}x_n\mbox{,}t)=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): U_{i\mbox{,}2}(x_1x_2\mbox{,}\ldots \mbox{,}x_n\mbox{,}t)= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \cdots =0 .
Therefore the exact solution may be obtained as the following.
|
(8) |
It is worthwhile to mention that if Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle g(x_1\mbox{,}x_2\mbox{,}\ldots \mbox{,}x_{n-1}\mbox{,}t)} , and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u_{i\mbox{,}0}(x_1\mbox{,}x_2\mbox{,}\ldots \mbox{,}x_{n-1}\mbox{,}t)} , are analytic around Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle t=t_0} , then their Taylor series can be defined as
|
(9) |
can be used in Eq. (7), where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle a_j(x)\mbox{,}\quad j=0\mbox{,}\ldots \mbox{,}n\mbox{,}b_j(x)\mbox{,}\quad j=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 0\mbox{,}\ldots \mbox{,}n , are unknown coefficients which must be computed, and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle a_j^{{_\ast}}(x)\mbox{,}\quad j=0\mbox{,}\ldots \mbox{,}n.b_j^{{_\ast}}(x)\mbox{,}\quad j=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 0\mbox{,}\ldots \mbox{,}n . are known ones.
The two most important steps in application of the new Homotopy perturbation method are to construct a suitable Homotopy equation and to choose a suitable initial guess.
To show the capability of the method, NHPM has been applied to some examples in the next section.
To demonstrate the effectiveness of the method two examples of system of Burgers equations are presented.
Consider the following system of Burgers equations
|
(10) |
with initial conditions
|
The exact solutions are
|
To solve Eq. (10), by the NHPM, we construct the following homotopies
|
(11) |
Applying the inverse operator, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle L^{-1}={\int }_0^t(.)dt}
to both sides of these equations, we obtain
|
(12) |
Suppose the solutions of system (12) are, as assumed in (5), substituting Eq. (5) into Eq. (12), collecting the same powers of p, and equating each coefficient of p to zero, results in
|
Assuming Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \begin{array}{l} u_0(x\mbox{,}t)={\sum }_{n=0}^{\infty }a_n(x)P_n(t)\mbox{,}\quad P_k(t)=t^k\mbox{,}\quad U(x\mbox{,}0)=u(x\mbox{,}0)\mbox{,}\\ v_0(x\mbox{,}t)={\sum }_{n=0}^{\infty }b_n(x)P_n(t)\mbox{,}\quad P_k(t)=t^k\mbox{,}\quad V(x\mbox{,}0)=v(x\mbox{,}0)\mbox{.} \end{array}}
And solving the above equations for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U_1(x\mbox{,}y\mbox{,}t)\mbox{,}V_1(x\mbox{,}y\mbox{,}t)}
leads to the result
|
|
By vanishing Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U_1(x\mbox{,}t)\mbox{,}V_1(x\mbox{,}t)}
coefficients Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle a_n(x)\mbox{,}b_n(x)\mbox{,}(n=1\mbox{,}2\mbox{,}3\mbox{,}\ldots )} are determined as
|
Therefore we gain the solution of Eq. (14) as
|
|
which is an exact solution (see Fig. 1).
|
Figure 1. Numerical result of Example 1. |
Consider the following system of coupled Burgers equations
|
(13) |
with initial conditions
|
The exact solution of Eq. (13) is
|
To solve Eq. (13), by the NHPM, we construct the following homotopy
|
(14) |
Applying the inverse operator, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle L^{-1}={\int }_0^t(.)dt}
to the both sides of the above equations, we obtain
|
(15) |
Suppose the solutions of system (15) have the form (5), substituting Eq. (5) into Eq. (15), collecting the same powers of p, and equating each coefficient of p to zero, results in
|
By assuming Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \begin{array}{ll} u_0(x\mbox{,}t)={\sum }_{n=0}^{\infty }a_n(x)t^k\mbox{,} & U(x\mbox{,}0)=u(x\mbox{,}0)\mbox{,}\\ v_0(x\mbox{,}t)={\sum }_{n=0}^{\infty }b_n(x)t^k\mbox{,} & V(x\mbox{,}0)=v(x\mbox{,}0)\mbox{,} \end{array}} and solving equations Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U_1(x\mbox{,}t)\mbox{,}V_1(x\mbox{,}t)}
leads to the following results
|
|
Vanishing Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U_1(x\mbox{,}t)\mbox{,}V_1(x\mbox{,}t)} , lets the coefficients Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle a_n(x)\mbox{,}b_n(x)\mbox{,}(n=1\mbox{,}2\mbox{,}3\mbox{,}\ldots )}
to take the following values
|
This implies that
|
In this example we have also derived exact solutions (see Fig. 2).
|
Figure 2. Numerical result of Example 2. |
In this paper, the modified HPM is applied to approximate solutions of Burgers equations. The present method reduces the computational difficulties of the other methods and all the calculations can be made with simple manipulations. The solutions introduced in this study can be used to obtain the closed form of the solutions if they are required. Thus, it can be concluded that the MHPM is very powerful and efficient in finding analytical as well as numerical solutions for wide classes of partial differential equations. The computations associated with the examples were performed using Maple 15.
The authors would like to thank the referees for their useful comments which lead to some improvements of the current paper.
Published on 11/04/17
Licence: Other
Are you one of the authors of this document?