MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2015/414808 414808 Research Article A New Exact Solution of Burgers’ Equation with Linearized Solution Kuo Chun-Ku 1, 2 Lee Sen-Yung 1 Rodriguez-Lopez Rosana 1 Department of Mechanical Engineering, National Cheng Kung University, Tainan 701 Taiwan ncku.edu.tw 2 Department of Mechanical Engineering, Air Force Institute of Technology, Kangshan 820 Taiwan afats.khc.edu.tw 2015 3082015 2015 24 03 2015 13 08 2015 3082015 2015 Copyright © 2015 Chun-Ku Kuo and Sen-Yung Lee. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper considers a general Burgers’ equation with the nonlinear term coefficient being an arbitrary constant. Two identical solutions of the general Burgers’ equation are separately derived by a direct integration method and the simplest equation method with the Bernoulli equation being the simplest equation. The proposed exact solutions overcome the long existing problem of discontinuity and can be successfully reduced to linearity, while the nonlinear term coefficient approaches zero. In addition, a general Cole-Hopf transform is introduced. Finally, the proposed derived solution is compared with the perturbation solution and other existing exact solutions. A new phenomenon, which we named “kink sliding,” is observed.

1. Introduction

Burgers’ equation was first introduced by Bateman  in 1915 and later analyzed by Burgers  in 1948. The equation is used as a model in many fields such as acoustics , continuous stochastic processes , dispersive water , shock waves , heat conduction , and turbulence . Burgers’ equation can also be considered as a simplified form of the Navier-Stokes equation  due to the form of the nonlinear convection term and the occurrence of a viscosity term.

Burgers’ equation is one of the very few nonlinear partial differential equations that can be solved exactly. When the absolute value of the nonlinear term coefficient of Burgers’ equation is 2 or 1, the exact solutions can be derived by the modified extended tanh-function method , the Exp-function method , the tanh-coth method , and the Cole-Hopf transform method . In this paper, the nonlinear term coefficient of Burgers’ equation is considered an arbitrary constant. The exact solutions of thisgeneral Burgers’ equation are derived by the four methods mentioned above . In addition, they are also separately derived by the simplest equation method with the Riccati equation being the simplest equation and a newly developed general Cole-Hopf transform method [11, 12]. It is shown that when the coefficient of the nonlinear term is the same as the regular Burgers’ equation, the derived exact solutions are the same as those in the existing literature. However, all these exact solutions do not satisfy the continuity condition and will not be reducible to linear solutions when the coefficient of the nonlinear term in the Burgers’ equation approaches zero.

In this paper, new exact solutions of the general Burgers’ equation are separately derived by a direct integration method and the simplest equation method with the Bernoulli equation being the simplest equation, respectively. The two exact solutions are shown to be the same. In addition, the newly derived solution can be successfully reduced to linearity, while the nonlinear term coefficient approaches zero. Finally, the proposed derived solution is compared with the perturbation solution and other existing exact solutions. Several numerical results are presented and illustrated.

2. General Burgers’ Equation

Consider the general Burgers’ equation:(1)ut+αuux-νuxx=0,where the nonlinear term coefficient α is an arbitrary constant and ν is the coefficient of the kinematics viscosity of fluids and is equal to 1/R. Further, R is the Reynolds number, and when it is large the equation describes shock wave behavior .

2.1. Existing Exact Solutions

Following the solution methods and procedures of the modified extended tanh-function method , the Exp-function method , and the tanh-coth method , the following sets of solutions have been obtained.

One set, when c>0, is a front wave (kink) solution,(2)ux,t=cα-cαtanhc2νx-ct,while the other is a traveling wave solution,(3)ux,t=cα-cαcothc2νx-ct.For c<0, the two traveling wave solutions are(4)ux,t=cα+cαtanc2νx-ct,(5)ux,t=cα-cαcotc2νx-ct.Here, c is the wave speed. When α=1, solution (2) is the same as that obtained by Ebaid ; however, when α=-2, ν=1, solutions (2) to (5) are the same as those obtained by Wazwaz . Many other solutions can also be found in Wazwaz’s book .

2.2. The Simplest Equation Method with Riccati Equation being the Simplest Equation

The simplest equation method with Riccati equation being the simplest equation is a method commonly used to develop the exact solution of some nonlinear equations. It can also be used to develop exact solutions of the general Burgers’ equation.

Herein, a new variable is introduced, ξ=x-ct, so that U(ξ)=u(x,t). General Burgers’ equation (1) is thus transformed into (6)-cUξ+αUUξ-νUξξ=0.

To solve the problem by the simplest equation method [17, 18], the solution can be assumed to be in the form of(7)Uξ=i=0LaiYξi,where ai are constants to be determined; Y is the exact solution of a chosen nonlinear ordinary differential equation, called the simplest equation; and L is a constant to be determined and the power of the specified solution function finite series, Y.

Now, consider the following Riccati equation as the simplest equation:(8)Yξ=μ1-Y2,where μ is wave number and μ>0. The exact solution for the chosen Riccati equation (8) is(9)Y=tanhμξ.

To determine the parameters, substitute (7) and (8) back into (6) and balance the linear terms of the derivative’s highest order, which yields Uξξ with the highest order nonlinear term UUξ in (6). Setting coefficients of like power of Y to be zeros yields the following equations:(10)L+22-1=L+L+2-1,Y0:-c+a0α=0,Y1:  a1α+2νμ=0,Y2:-c-a0α=0,Y3:-a1α-2νμ=0.

The equations give L=1 and a0=c/α, a1=-2νμ/α. Consequently, the exact solution of the Burgers’ equation is obtained:(11)u=cα-2νμαtanhμx-ct.

If one sets μ=c/2ν, then c>0. It can accordingly be found that solution (11) is exactly the same as that given in (2).

2.3. General Cole-Hopf Transformation Solution

When α=-2 and α=-1, regular Cole-Hopf transformations, u=ν/xlnF, and u=2ν/xlnF were introduced to study the associated Burgers’ equations [11, 12, 15], respectively. In this paper, a general Cole-Hopf transformation solution is developed to derive an exact solution of the general Burgers’ equation.

The regular Cole-Hopf transformation in the following form is first modified:(12)u=-βνxlnF,where β is a constant to be specified. Substituting (12) into (1) yields(13)-βνxFtF+βν2FxxFxF2αβ-3+FxF32-αβ+FxxxFF2=0.

Obviously, if αβ=2, (13) can be reduced to (14)-βνxFtF+βν2-FxxFxF2+FxxxFF2=βνxFt-νFxxF=0.Consequently, the following linear equation is generated:(15)Ft-νFxx=0.Moreover, it is a heat conduction equation with a solution in the form of [11, 12](16)F=1+ek1x+νk12t,where k1 is an arbitrary constant.

Therefore, the general Cole-Hopf transformation for the Burgers’ equation is (17)u=-2ναxlnF,and an exact solution of the general Burgers’ equation can be derived as(18)u=-2ναk1ek1x+νk12t1+ek1x+νk12t.

From the solutions above, it can be observed that the tanh-coth solutions, (2)–(5), the Cole-Hopf solution, (18), the solution derived from the simplest equation method with the Riccati equation being the simplest equation, and (11) are all proportional to 1/α. In particular, when ν=1, the Cole-Hopf solution, namely (18), with k1=-1 is the same as one of the tanh-coth solutions, namely (2), with c=1.

It is well known that, for solution continuity, if the nonlinear coefficient α in a nonlinear differential equation (1) approaches zero, the equation is reduced to a linear one and the nonlinear solution should be reducible to a linear solution. However, from the existing exact solutions above, it can be observed that when the nonlinear coefficient α is reduced to zero, the solutions become infinite. Therefore, none will satisfy the solution continuity condition. In the following, a new exact solution of the general Burgers’ equation with a linearized solution is developed.

3. New Exact Solution 3.1. Direct Integration Method

Integrating (6) once with respect to ξ and setting the integration constant to be zero yields(19)νUξ-α2U2+cU=0.Equation (19) can then be expressed as(20)Uξ=Up+qU,where p=-c/ν, q=α/2ν.

Integrating (20) once with respect to ξ and using the result given in the integration table , it becomes(21)-1plogp+qUU=ξ+ξ0.Therefore(22)U=pepξ+ξ01-qepξ+ξ0,where ξ0 is an integral constant. As a result, a new exact solution of the general Burgers’ equation can be derived as (23)ux,t=-2c2νec/νx-ct+ξ0-α.

When the nonlinear coefficient α in nonlinear differential equation (1) approaches zero, the exact solution, namely (23), is reduced to (24)ux,t=-cνe-c/νx-ct+ξ0.This solution is consistent with the exact solution of the following linearized equation from Burgers’ equation:(25)ut-νuxx=0.

3.2. The Simplest Equation Method with the Bernoulli Equation being the Simplest Equation

Following the simplest equation method as discussed in Section 2.2, the same new exact solution can also be obtained by employing the simplest equation method with the following Bernoulli equation being the simplest equation :(26)Yξ=aY+bY2.Equation (26) has the exact solution in the form of (27)Y=aeaξ+ξ01-beaξ+ξ0.

After substituting (7) and (26) back into (6) and balancing the linear terms of the derivative’s highest order, Uξξ, with the highest order nonlinear term UUξ in (6), L=1, and the following equations are obtained(28)Y1:-ca1a+αa0a1a-νa1a2=0,Y2:-ca1b+αa0a1+αa12a-3νa1ab=0,Y3:  αa12b-2νa1b2=0.

These equations give a0=0, a1=1, and a=-c/ν, b=α/2ν. Consequently, the solution is derived as(29)Uξ=a0+a1Y=-c/νe-c/νξ+ξ01-α/2νe-c/νξ+ξ0.In terms of x and t variables, the new exact solution of the general Burgers’ equation can be expressed as (30)ux,t=-c/νe-c/νx-ct+ξ01-α/2νe-c/νx-ct+ξ0.After multiplying the numerator and denominator of (30) by 2νec/ν(x-ct+ξ0), the same solution as that given in (23) is obtained.

3.3. The Perturbation Solution

The perturbation method is commonly used to study nonlinear equations and obtain approximate analytic solutions. The transformed general Burgers’ equation, (6), can be set into the form (31)-cUξ+εUUξ-νUξξ=0,where ε=α is a perturbation term. The perturbation solution is assumed to be in the following form:(32)U=U0+εU1+ε2U2+.

Substituting (32) into (31) and setting like power coefficients of ε to be zeros, the following equations are obtained:(33)ε0:  cU0ξ+νU0ξξ=0,ε1:  cU1ξ+νU1ξξ=U0U0ξ,ε2:  cU2ξ+νU2ξξ=U0U1ξ+U1U0ξ.The approximated solution is(34)U=-C1cνe-c/νξ+ξ0+ε-cνe-c/νξ+ξ0+12ce-2c/νξ+ξ0+ε2-cνe-c/νξ+ξ0+14c2e-3c/νξ+ξ0+,where C1 is a constant. When the coefficient of the perturbation term ε is zero, (31) is reduced to a common linear differential equation. Solution (34) is reduced to being linear and is the same as the solution reduced from (24).

4. Numerical Results and Discussions

To illustrate the analysis, several numerical results are presented and illustrated.

Figure 1 illustrates the difference between the new solution and the perturbation solution, with ν=1 and c=1, at ξ=0, as the nonlinear term coefficient α(ε) varies from −1 to 1. For comparison, the constant ξ0 in (23) and (34) is set to zero. It can be found that when α=0, both solutions lead to linear solutions and are the same; however, when the absolute value of α is increased to one, the perturbation solution will have significant errors. Figure 2 shows the difference between the new exact solution and the perturbation solution with the absolute value of α being one, ν=1 and c=1, at various ξ. It can be observed that the difference between the two solutions becomes negligible when ξ is large (ξ>4).

Influence of the nonlinear coefficient α(ε) on the difference between the new exact solution and the perturbation solution [c=1, ν=1, and ξ=0; –⋅–: the perturbation solution; —: the new exact solution].

Difference between the new exact solution and the perturbation solution with various ξ [ν=1, c=1, and (a) α=1; (b) α=-1; –⋅–: the perturbation solution; —: the new exact solution].

A kink wave is defined as a traveling wave which (1) rises or descends from one asymptotic state to another and (2) approaches a constant at infinity . Moreover, when x approaches positive and negative infinities, the corresponding tanh-coth solutions from (2) are (35)Utanh-coth=cα-cαtanhc2ν=0,c>0;Utanh-coth-=cα-cαtanhc2ν-=2cα,c>0.The corresponding Cole-Hopf solutions from (18), at x=±, are(36)UCole-Hopf=-2ναk1ek1+νk12t1+ek1+νk12t=-2νk1α,k1>0;UCole-Hopf-=-2ναk1ek1-+νk12t1+ek1-+νk12t=0,k1>0;UCole-Hopf=-2ναk1ek1+νk12t1+ek1+νk12t=0,k1<0;UCole-Hopf-=-2ναk1ek1-+νk12t1+ek1-+νk12t=-2νk1α,k1<0.The corresponding new exact solutions from (23), at x=±, are(37)Unew=-2c2νec/ν-α=0,c>0;Unew-=-2c2νec/ν--α=2cα,c>0;Unew=-2c2νec/ν-α=2cα,c<0;Unew-=-2c2νec/ν--α=0,c<0.Hence, they satisfy the second condition of a kink wave. After numerical calculation using the software MATLAB, it can be found that the solutions of (1) with α<0, (2), and (18) all satisfy the first condition of a kink wave. The kinds of waves of the tanh-coth solutions, (2)–(5), the Cole-Hopf solution, (18), and the new exact solution, (23), are summarized in Table 1.

Waves types of the tanh-coth solutions, the Cole-Hopf solution, and the new exact solution of the general Burgers’ equation.

Parameters Exact solutions
α c The tanh-coth solutions The Cole-Hopf solution The new exact solution
Equation (2) Equation (3) Equation (4) Equation (5)
>0 >0 Kink Traveling wave NA NA Kink Traveling wave
<0 >0 Kink Traveling wave NA NA Kink Kink
>0 <0 NA NA Traveling wave Traveling wave Kink Traveling wave
<0 <0 NA NA Traveling wave Traveling wave Kink Kink

In the following, the three different kinds of solutions, (2), (18), and (23), are examined and compared. As the nonlinear term coefficient of α varies from −1 to 1, the difference between the new exact solution and the Cole-Hopf solution, (18), with ν=1, c=1, k1=-1, x=1, and t=1, is shown in Figure 3. It should be mentioned that under this case, the tanh-coth solution, (2), is the same as the Cole-Hopf solution. It can be found that the new exact solution is a continuous function of α while the other solution has discontinuity at α=0. This phenomenon can be easily observed from the solution forms, namely, (2) and (18). When the nonlinear coefficient α in the tanh-coth solution and the Cole-Hopf solution is reduced to zeros, the solutions approach infinity.

Comparison among the new exact solution, the tanh-coth solution, and the Cole-Hopf solution as the nonlinear coefficient α varies from −1 to 1 [t=1, x=1, c=1, ν=1, and k1=-1; —: the new exact solution; – –: the Cole-Hopf solution and the tanh-coth solution].

Figure 4 shows the influence of the kinematics viscosity coefficient ν on the three solutions as x=1, c=1, t=1, k1=-1, and α=-2 or -1. As can be seen, the coefficient ν has no influence on the tanh-coth solution. As the coefficient ν is increased, the new exact solution increases and the Cole-Hopf solution decreases, respectively. Under the unique case with c=1, ν=1, k1=-1, and α=-2, the three different solutions are the same.

Influence of the kinematics viscosity coefficient ν on the three solutions [x=1, c=1, t=1, k1=-1, (a) α=-2, and (b) α=-1; —: the new exact solution; ⋯: the tanh-coth solution; – –: the Cole-Hopf solution].

Figure 5 shows the new exact solution, the tanh-coth solution, and the Cole-Hopf solution with various x and different negative nonlinear coefficients α. It can be found that these three solutions are all kink solutions and are consistent with the summarized results in Table 1. As mentioned in Section 2.3, with constants x, c, t, and k1, the tanh-coth solution and the Cole-Hopf solution are both proportional to 1/α and become infinite as α approaches zero. The steep slope in the solutions is kept at the same location; however, this is not the case for the new exact solution. Under the same physical conditions, the new exact solution will not become infinite as α approaches zero. Moreover, the location of the steep slope in the new exact solution will slide and change as α is changed. As shown in Figure 5, the location of the steep slope in the new exact solution slides to the left as α is changed from −1 to −0.1. We have named this peculiarity the kink-sliding phenomenon.

The new exact solution, the tanh-coth solution, and the Cole-Hopf solution with various x and different nonlinear coefficients α [c=1.2, t=1, ν=1, k1=-1, (a) α=-1, and (b) α=-0.1; —: the new exact solution; ⋯: the tanh-coth solution; – –: the Cole-Hopf solution].

5. Conclusions

In this paper, a new exact solution of the general Burgers’ equation was derived by a direct integration method and the simplest equation method with the Bernoulli equation as the simplest equation. It was shown that the proposed derived exact solution can be successfully reduced to linearity, while the nonlinear term coefficient approaches zero. In addition, a general Cole-Hopf transformation solution was introduced. The differences among the newly derived exact solution, the other existing exact solutions, and the perturbation solution were elaborated. It was shown that when α=0, both the new solution and perturbation solution lead to linearity and are the same; however, when the absolute value of α is increased to one, the perturbation solution has significant errors. With constants x, c, t, and k1, the tanh-coth solutions and the Cole-Hopf solution are proportional to 1/α and become infinite as α approaches zero. Although the steep slope in the solutions is maintained at the same location, under the same physical conditions, the new exact solution will not become infinite as α approaches zero. Moreover, the location of the steep slope in the new exact solution slides, which was termed the kink-sliding phenomenon. Extensions of the proposed methods to study different nonlinear differential equations without linearized solutions are expected in future works.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

It is gratefully acknowledged that this work was supported by the Ministry of Science and Technology of Taiwan, under Grants NSC 102-2221-E-006-050 and MOST 103-2221-E-006-048.

Bateman H. Some recent researches on the motion of fluids Monthly Weather Review 1915 43 4 163 170 10.1175/1520-0493(1915)4360;163:srrotm62;2.0.co;2 Burgers J. M. A mathematical model illustrating the theory of turbulence Advances in Applied Mechanics 1948 1 171 199 Sugimoto N. Burgers equation with a fractional derivative; hereditary effects on nonlinear acoustic waves Journal of Fluid Mechanics 1991 225 631 653 10.1017/s0022112091002203 2-s2.0-0026135818 Weinan E. Khanin K. Mazel A. Sinai Y. Invariant measures for Burgers equation with stochastic forcing Annals of Mathematics 2000 151 3 877 960 10.2307/121126 ZBL0972.35196 2-s2.0-0034347507 Whitham G. B. Linear and Nonlinear Waves 2011 42 John Wiley & Sons Pure and Applied Mathematics 10.1002/9781118032954 MR1699025 Bluman G. W. Cole J. D. The general similarity solution of the heat equation Journal of Mathematics and Mechanics 1969 18 1025 1042 MR0293257 Bec J. Khanin K. Burgers turbulence Physics Reports 2007 447 1-2 1 66 10.1016/j.physrep.2007.04.002 2-s2.0-34250782224 Bell J. B. Colella P. Glaz H. M. A second-order projection method for the incompressible navier-stokes equations Journal of Computational Physics 1989 85 2 257 283 10.1016/0021-9991(89)90151-4 2-s2.0-0001608079 Soliman A. A. The modified extended tanh-function method for solving Burgers-type equations Physica A: Statistical Mechanics and Its Applications 2006 361 2 394 404 10.1016/j.physa.2005.07.008 2-s2.0-31044444121 Ebaid A. Exact solitary wave solutions for some nonlinear evolution equations via Exp-function method Physics Letters A 2007 365 3 213 219 10.1016/j.physleta.2007.01.009 2-s2.0-34147152832 Wazwaz A.-M. Multiple-front solutions for the Burgers equation and the coupled Burgers equations Applied Mathematics and Computation 2007 190 2 1198 1206 10.1016/j.amc.2007.02.003 2-s2.0-34250643151 Wazwaz A. M. Partial Differential Equations and Solitary Waves Theory 2010 Springer Hopf E. The partial differential equation ut+uux=uxx Communications on Pure and Applied Mathematics 1950 3 201 230 Cole J. D. On a quasi-linear parabolic equation occurring in aerodynamics Quarterly of Applied Mathematics 1951 9 225 236 MR0042889 Ohwada T. Cole-Hopf transformation as numerical tool for the burgers equation Applied and Computational Mathematics 2009 8 1 107 113 ZBL1187.65116 2-s2.0-70249089103 Vaganan B. M. Cole-Hopf transformations for higher dimensional Burgers equations with variable coefficients Studies in Applied Mathematics 2012 129 3 300 308 10.1111/j.1467-9590.2012.00551.x 2-s2.0-84867321022 Kudryashov N. A. Seven common errors in finding exact solutions of nonlinear differential equations Communications in Nonlinear Science and Numerical Simulation 2009 14 9-10 3507 3529 10.1016/j.cnsns.2009.01.023 2-s2.0-63449123795 Kudryashov N. A. One method for finding exact solutions of nonlinear differential equations Communications in Nonlinear Science and Numerical Simulation 2012 17 6 2248 2253 10.1016/j.cnsns.2011.10.016 ZBL1250.35055 2-s2.0-84855223308 Spiegel M. R. Seymour L. John L. Schaum's Outlines Mathematical Handbook of Formulas and Tables 1968 Schaum Publications Vitanov N. K. Application of simplest equations of Bernoulli and Riccati kind for obtaining exact traveling-wave solutions for a class of PDEs with polynomial nonlinearity Communications in Nonlinear Science and Numerical Simulation 2010 15 8 2050 2060 10.1016/j.cnsns.2009.08.011 MR2592618 2-s2.0-76049089677 Vitanov N. K. Modified method of simplest equation: powerful tool for obtaining exact and approximate traveling-wave solutions of nonlinear PDEs Communications in Nonlinear Science and Numerical Simulation 2011 16 3 1176 1185 10.1016/j.cnsns.2010.06.011 2-s2.0-77957345115 Vitanov N. K. On modified method of simplest equation for obtaining exact and approximate solutions of nonlinear PDEs: the role of the simplest equation Communications in Nonlinear Science and Numerical Simulation 2011 16 11 4215 4231 10.1016/j.cnsns.2011.03.035 2-s2.0-79957892557