Exp-function method for some nonlinear PDE’s and a nonlinear ODE’s

Kourosh Parand; Jamal Amani Rad

doi:10.1016/j.jksus.2010.08.004

View/Download PDF

Buy Reprints

PDF

Translate this page into:

ORIGINAL ARTICLE

24 (

); 1-10

doi:

10.1016/j.jksus.2010.08.004

Exp-function method for some nonlinear PDE’s and a nonlinear ODE’s

Kourosh Parand^*, Jamal Amani Rad

Department of Computer Sciences, Shahid Beheshti University, G.C., Tehran, Iran

*Corresponding author. Tel.: +98 21 22431653; fax: +98 21 22431650 k_parand@sbu.ac.ir (Kourosh Parand),

Received: 2010-6-25, Accepted: 2010-8-2, Published: 2010-01-01

Disclaimer:
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.

Available online 11 August 2010

Abstract

In this paper, we apply the Exp-function method to find some exact solutions for two nonlinear partial differential equations (NPDE) and a nonlinear ordinary differential equation (NODE), namely, Cahn-Hilliard equation, Allen-Cahn equation and Steady-State equation, respectively. It has been shown that the Exp-function method, with the help of symbolic computation, provides a very effective and powerful mathematical tool for solving NPDE’s and NODE’s. Mainly we try to present an application of Exp-function method taking to consideration rectifying a commonly occurring errors during some of recent works. The results of the other methods clearly indicate the reliability and efficiency of the used method.

Keywords

The Cahn-Hilliard equation

Allen-Cahn equation

Steady-State equation

Exp-function method

NPDE

NODE

Show Related Articles from PubMed

1 Introduction

The study of exact solutions of nonlinear partial differential equations (NPDE) plays an important role in mathematical physics, engineering and the other sciences. In the past several decades, various methods for obtaining solutions of NPDEs and ODE’s have been presented, such as, tanh-function method (Wazwaz, 2005, 2006a,b ), Adomian decomposition method (Hashim et al., 2006; Tatari et al., 2007), Homotopy perturbation method (Rashidi et al., 2009; Biazar et al., 2009; Berberler and Yildirim, 2009 ), variational iteration method (Shakeri and Dehghan, 2008; Soliman and Abdou, 2007; Yusufoglu and Bekir, 2007 ), spectral method (Parand and Taghavi, 2009; Parand et al., 2009, 2010 ), sine–cosine method (Tascan and Bekir, 2009; Wazwaz, 2007), radial basis method (Tatari and Dehghan, 2010; Dehghan and Shokri, 2009) and so on. Recently, He and Wu (2006) proposed a novel method, so called Exp-function method, which is easy, succinct and powerful to implement to nonlinear partial differential equations arising in mathematical physics. The Exp-function method has been successfully applied to many kinds of NPDEs, such as, KdV equation with variable coefficients (Zhang, 2007), Maccari’s system (Zhang, 2007), Boussinesq equations (Abdou et al., 2007), Burger’s equations (Ebaid, 2007; Biazar and Ayati, 2009; Ebaid, 2009 ), Double Sine-Gordon equation (Domairry et al., 2010; He and Abdou, 2007), Fisher equation (Ozis and Koroglu, 2008), Jaulent–Miodek equations (He and Zhang, 2008) and the other important nonlinear partial differential equations (Koroglu and zis, 2009; Shin et al., 2009; Zhang, 2008 ). Recently, in some of papers applying the Exp-function method (He and Wu, 2006) have been occurred with common errors. Seven common errors are formulated and classified by (Kudryashov, 2009). In this paper we try to apply this method taking to rectifying these common errors to look exact solutions of three nonlinear differential equations, namely, Cahn-Hilliard equation, Allen-Cahn equation and Steady-State equation given by $\begin{matrix} u_{t} = γ u_{x} + 6 u (u_{x})^{2} + (3 u^{2} - 1) u_{xx} - u_{xxxx}, \\ u_{t} = u_{xx} - u^{3} + u, \\ α u^{″} (x) = β u (x) (u (x) - m) (u (x) + m), \end{matrix}$ respectively, which $α, β, m$ and $γ$ are the constants.

The Cahn–Hilliard equation was proposed to describe phase separation phenomena in binary systems (Cahn et al., 1958). This equation is related with a number of interesting physical phenomena like the spinodal decomposition, phase separation and phase ordering dynamics. It is also very crucial in material sciences (Chan, 1961; Choo et al., 2004; Gurtin, 1996 ). On the other hand, this equation is very hard and difficult to solve. The Cahn-Hilliard equation has been extensively studied by Wang and Shi (1993), Jabbari and Peppas (1995), Puri and Binder (1991) for the study of interfaces. Global existence and uniqueness of the solution have been shown by Elliott and Zheng (1986). Jingxue (1992) has shown the existence of continuous solution for the problem with degenerate mobility. Recently, Dehghan and Mirzaei (2009) applied a numerical method based on the boundary integral equation and dual reciprocity methods for one-dimensional Cahn-Hilliard equation. Ugurlu and Kaya (2008) solved Cahn-Hilliard equation by tanh-function method. Furihata (2001) applied finite difference for Cahn-Hilliard equation. Many articles have investigated this equation mathematically and numerically this equation (Mello et al., 2005; Kim, 2007; Wells et al., 2006 ). Also, Allen-Cahn equation arise in many scientific applications such as mathematical biology, quantum mechanics and plasma physics. It is well known that wave phenomena of plasma media and fluid dynamics are modelled by kink shaped and tanh solution or bell shaped sech solutions (Wazwaz, 2007; Tascana and Bekir, 2009).

The rest of the paper is organized as follows: Section 2 describes Exp-function method for finding exact solutions to the NPDEs. The applications of the proposed analytical scheme presented in Section 3. The conclusions are discussed in the Section 4. exp-function calculations are provided in the end.

2 Basic idea of Exp-function method

We consider a general nonlinear PDE in the following form

(1)

N (u, u_{x}, u_{t}, u_{xx}, u_{tt}, u_{xt}, \dots) = 0,

where N is a polynomial function with respect to the indicated variables or some functions which can be reduced to a polynomial function by using some transformation. We introduce a complex variation as

(2)

u (x, t) = U (η), η = kx + ω t,

where k and

ω

are constants. We can rewrite Eq. (1) in the following nonlinear ordinary differential equations

N (U, {kU}^{'}, ω U^{'}, k^{2} U^{″}, \dots) = 0,

where the prime denotes the derivation with respect to

η

. According to the Exp-function method (He and Wu, 2006), we assume that the solution can be expressed in the form

(3)

U (η) = \frac{\sum_{i = - d}^{c} a_{i} \exp (i η)}{\sum_{j = - q}^{p} b_{j} \exp (j η)},

where

c, d, p

and q are positive integers which can be freely chosen,

a_{i}

and

b_{j}

are unknown constants to be determined. To determine the values of c and p, we balance the highest order linear term with the highest order nonlinear term in Eq. (3). Similarly to determine the values of d and q. So by means of the exp-function method, we obtain the generalized solitary solution and periodic solution for nonlinear evolution equations arising in mathematical physics.

3 Applications of the Exp-function method

Example 1

Let us consider the Cahn-Hilliard equation (Ugurlu and Kaya, 2008; Dahmani and Benbachir, 2009) in the form

(4)

u_{t} = γ u_{x} + 6 u (u_{x})^{2} + (3 u^{2} - 1) u_{xx} - u_{xxxx},

that, by using the complex variation

(5)

u (x, t) = U (η), η = kx + ω t,

and integrating with respect to

η

, Eq. (4) can be converted to the ODE (for

γ = 1

)

(6)

(ω - k) U + k^{4} U^{‴} - 3 k^{2} U^{2} U^{'} + k^{2} U^{'} = 0,

where the prime denotes the derivative with respect to

η

and also where the integration constant is chosen as zero. In other words, we are solved this problem for the case when integration constant is zero. In view of the Exp-function method, we assume that the solution of Eq. (6) can be expressed in the form

U (η) = \frac{a_{c} \exp (c η) + \dots + a_{- d} \exp (- d η)}{b_{p} \exp (p η) + \dots + b_{- q} \exp (- q η)},

where

c, d, p

and q are positive integers which are unknown to be determined later. In order to determine the values of c and p, we balance the linear term of the highest order with the highest order nonlinear terms in Eq. (6) i.e.

U^{‴}

and

U^{2} U^{'}

. By simple calculation, we have

(7)

U^{‴} = \frac{c_{1} \exp [(c + 7 p) η] + \dots}{c_{2} \exp [8 p η] + \dots},

and

(8)

U^{2} U^{'} = \frac{c_{3} \exp [(3 c + p) η] + \dots}{c_{4} \exp [4 p η] + \dots} = \frac{c_{3} \exp [(3 c + 5 p) η] + \dots}{c_{4} \exp [8 p η] + \dots},

where

c_{i}

are determined coefficients only for simplicity. By balancing the highest order of Exp-function in Eqs. (8) and (7), we have

c + 7 p = 3 c + 5 p,

which leads to the result

p = c .

Similarly, to determine the values of d and q, we balance the linear term of the lowest order in Eq. (6)

(9)

U^{‴} = \frac{\dots + d_{1} \exp [- (7 q + d) η]}{\dots + d_{2} \exp [- 8 q η]},

and

(10)

U^{2} U^{'} = \frac{\dots + d_{3} \exp [- (q + 3 d) η]}{\dots + d_{4} \exp [- 4 q η]} = \frac{\dots + d_{3} \exp [- (5 q + 3 d) η]}{\dots + d_{4} \exp [- 8 q η]},

where

d_{i}

are determined coefficients only for simplicity, we have

- (7 q + d) = - (5 q + 3 d),

which leads to results

q = d .

For simplicity, we set

p = c = 1

and

q = d = 1

, so Eq. (3) reduces to

(11)

U (η) = \frac{a_{1} \exp (η) + a_{0} + a_{- 1} \exp (- η)}{\exp (η) + b_{0} + b_{- 1} \exp (- η)} .

Substituting Eq. (11) into Eq. (6), equating to zero the coefficients of all powers of

\exp (n η)

yields a set of algebraic equations for

a_{0}, b_{0}, a_{- 1}, a_{1}, b_{1}, k

and

ω

(see Appendix A). By solving the system of algebraic equations with a professional mathematical software, we obtain

Case 1.

(12)

\{\begin{matrix} a_{- 1} = b_{- 1}, & b_{- 1} = b_{- 1}, \\ a_{0} = 0, & b_{0} = 0, \\ a_{1} = - 1, & k = ω = \frac{\sqrt{2}}{2} . \end{matrix}

where

b_{- 1}

is free parameter which can be determined by initial or boundary conditions. Substituting these results into (11), we obtain

(13)

u (x, t) = \frac{- \exp [\frac{\sqrt{2}}{2} (x + t)] + b_{- 1} \exp [- \frac{\sqrt{2}}{2} (x + t)]}{\exp [\frac{\sqrt{2}}{2} (x + t)] + b_{- 1} \exp [- \frac{\sqrt{2}}{2} (x + t)]} .

These results cover some of special solutions Eq. (4) regarding to initial value conditions. For example, if initial conditions

u (x, t)

becomes

u (x, 0) = 0

, then by substituting

t = 0

in the Eq. (13) we could obtain

b_{- 1}

and exact solution of Eq. (4). If Eq. (4) were in the following form:

(14)

\{\begin{matrix} u_{t} = u_{x} + 6 u (u_{x})^{2} + (3 u^{2} - 1) u_{xx} - u_{xxxx}, \\ u (x, 0) = - \tanh [\frac{\sqrt{2}}{2} x] . \end{matrix}

From Eqs. (14) and (13), we obtain

(15)

b_{- 1} = 1 .

Thus, from substituting Eq. (15) into Eq. (13), we obtain

u (x, t) = \frac{- \exp [\frac{\sqrt{2}}{2} (x + t)] + \exp [- \frac{\sqrt{2}}{2} (x + t)]}{\exp [\frac{\sqrt{2}}{2} (x + t)] + \exp [- \frac{\sqrt{2}}{2} (x + t)]} = - \tanh [\frac{\sqrt{2}}{2} (x + t)],

which is the solution obtained by tanh method in Ugurlu and Kaya, 2008. Also, if Eq. (4) were in the following form:

(16)

\{\begin{matrix} u_{t} = u_{x} + 6 u (u_{x})^{2} + (3 u^{2} - 1) u_{xx} - u_{xxxx}, \\ u (x, 0) = - \coth [\frac{\sqrt{2}}{2} x] . \end{matrix}

From Eqs. (16) and (13), we derive

(17)

b_{- 1} = - 1 .

Thus, from substituting Eq. (17) into Eq. (13), we obtain

u (x, t) = \frac{- \exp [\frac{\sqrt{2}}{2} (x + t)] - \exp [- \frac{\sqrt{2}}{2} (x + t)]}{\exp [\frac{\sqrt{2}}{2} (x + t)] - \exp [- \frac{\sqrt{2}}{2} (x + t)]} = - \coth [\frac{\sqrt{2}}{2} (x + t)],

which is the solution obtained by tanh method in Ugurlu and Kaya (2008).

Case 2.

(18)

\{\begin{matrix} a_{- 1} = - b_{- 1}, & b_{- 1} = b_{- 1}, \\ a_{0} = 0, & b_{0} = 0, \\ a_{1} = 1, & k = ω = \frac{\sqrt{2}}{2} . \end{matrix}

Substituting these result into Eq. (11), we obtain the following solution:

(19)

u (x, t) = \frac{\exp [\frac{\sqrt{2}}{2} (x + t)] - b_{- 1} \exp [- \frac{\sqrt{2}}{2} (x + t)]}{\exp [\frac{\sqrt{2}}{2} (x + t)] + b_{- 1} \exp [- \frac{\sqrt{2}}{2} (x + t)]} .

If Eq. (4) were in the following form:

(20)

\{\begin{matrix} u_{t} = u_{x} + 6 u (u_{x})^{2} + (3 u^{2} - 1) u_{xx} - u_{xxxx}, \\ u (x, 0) = \coth [\frac{\sqrt{2}}{2} x] . \end{matrix}

From Eq. (20) and Eq. (19), we obtain

(21)

b_{- 1} = 1 .

Thus, from substituting Eq. (21) into Eq. (19), we obtain

u (x, t) = \frac{\exp [\frac{\sqrt{2}}{2} (x + t)] - \exp [- \frac{\sqrt{2}}{2} (x + t)]}{\exp [\frac{\sqrt{2}}{2} (x + t)] + \exp [- \frac{\sqrt{2}}{2} (x + t)]} = \coth [\frac{\sqrt{2}}{2} (x + t)],

which is the same as Ugurlu’s solution (Ugurlu and Kaya, 2008). Also, if Eq. (4) were in the following form:

(22)

\{\begin{matrix} u_{t} = u_{x} + 6 u (u_{x})^{2} + (3 u^{2} - 1) u_{xx} - u_{xxxx}, \\ u (x, 0) = \tanh [\frac{\sqrt{2}}{2} x] . \end{matrix}

From Eq. (22) and Eq. (19), we obtain

(23)

b_{- 1} = - 1 .

Thus, from substituting Eq. (23) into Eq. (19), we obtain

(24)

u (x, t) = \frac{\exp [\frac{\sqrt{2}}{2} (x + t)] + \exp [- \frac{\sqrt{2}}{2} (x + t)]}{\exp [\frac{\sqrt{2}}{2} (x + t)] - \exp [- \frac{\sqrt{2}}{2} (x + t)]} = \tanh [\frac{\sqrt{2}}{2} (x + t)],

which is the exact solutions given by Ugurlu and Kaya (2008), Dahmani and Benbachir (2009). Fig. 1 depict the solution Eq. (24).

Case 3.

(25)

\{\begin{matrix} a_{- 1} = 0, & b_{- 1} = \frac{1}{8} a_{0}^{2}, \\ a_{0} = a_{0}, & b_{0} = 0, \\ a_{1} = 0, & k = ω = i . \end{matrix}

In this case, k and

ω

are imaginary numbers. Substituting these results into (11), we obtain

(26)

u (x, t) = \frac{a_{0}}{\exp [i (x + t)] + \frac{1}{8} a_{0}^{2} \exp [- i (x + t)]} .

If Eq. (4) were in the following form:

(27)

\{\begin{matrix} u_{t} = u_{x} + 6 u (u_{x})^{2} + (3 u^{2} - 1) u_{xx} - u_{xxxx}, \\ u (x, 0) = \frac{\sqrt{2}}{2} \tan [\frac{1}{2} x] + \frac{\sqrt{2}}{2} \cot [l \frac{1}{2} x] . \end{matrix}

From Eqs. (27) and (26), we obtain

(28)

a_{0} = 2 \sqrt{2} i .

Thus, from substituting Eq. (28) into Eq. (26), we have

(29)

u (x, t) = \frac{2 \sqrt{2} i}{\exp [i (x + t)] - \exp [- i (x + t)]},

we know

(30)

\exp [i (x + t)] - \exp [- i (x + t)] = 2 i \sin (x + t),

and

(31)

\frac{1}{\sin (x + t)} = \frac{1}{2} \tan [\frac{1}{2} (x + t)] + \frac{1}{2} \cot [\frac{1}{2} (x + t)] .

Substituting Eqs. (30) and (31) into Eq. (29), we obtain

(32)

u (x, t) = \frac{\sqrt{2}}{2} \tan [\frac{1}{2} (x + t)] + \frac{\sqrt{2}}{2} \cot [\frac{1}{2} (x + t)],

which is the solution obtained by tanh method in Ugurlu and Kaya (2008). Fig. 2 depict the solution Eq. (32). Also, if Eq. (4) were in the following form:

(33)

\{\begin{matrix} u_{t} = u_{x} + 6 u (u_{x})^{2} + (3 u^{2} - 1) u_{xx} - u_{xxxx}, \\ u (x, 0) = - \frac{\sqrt{2}}{2} \tan [\frac{1}{2} x] - \frac{\sqrt{2}}{2} \cot [\frac{1}{2} x] . \end{matrix}

From Eqs. (33) and (26), we obtain

(34)

a_{0} = - 2 \sqrt{2} i .

Thus, from substituting Eq. (34) into Eq. (26), we have

(35)

u (x, t) = \frac{- 2 \sqrt{2} i}{\exp [i (x + t)] - \exp [- i (x + t)]} .

Substituting Eqs. (30) and (31) into Eq. (35), we obtain

u (x, t) = - \frac{\sqrt{2}}{2} \tan [\frac{1}{2} (x + t)] - \frac{\sqrt{2}}{2} \cot [\frac{1}{2} (x + t)],

which is the solution obtained by tanh method in Ugurlu and Kaya (2008).

Case 4.

(36)

\{\begin{matrix} a_{- 1} = - \frac{1}{8} \frac{b_{0}^{2} (a_{1}^{2} - 1)}{a_{1}}, & b_{- 1} = - \frac{1}{8} \frac{b_{0}^{2} (a_{1}^{2} - 1)}{a_{1}^{2}}, \\ a_{0} = - \frac{(2 a_{1}^{2} - 1) b_{0}}{a_{1}}, & b_{0} = b_{0}, \\ a_{1} = a_{1}, & k = ω = \sqrt{3 a_{1}^{2} - 1}, \end{matrix}

where

b_{0}

and

a_{1}

are free parameters. Inserting these result into Eq. (11), we obtain

\begin{matrix} u (x, t) = \frac{a_{1} \exp (η) - \frac{(2 a_{1}^{2} - 1) b_{0}}{a_{1}} - \frac{1}{8} \frac{b_{0}^{2} (a_{1}^{2} - 1)}{a_{1}} \exp (- η)}{\exp (η) + b_{0} - \frac{1}{8} \frac{b_{0}^{2} (a_{1}^{2} - 1)}{a_{1}^{2}} \exp (- η)} \\ = a_{1} - \frac{3 a_{1} b_{0} - \frac{b_{0}}{a_{1}}}{\exp (η) + b_{0} - \frac{1}{8} \frac{b_{0}^{2} (a_{1}^{2} - 1)}{a_{1}^{2}} \exp (- η)}, \end{matrix}

with

η = \sqrt{3 a_{1}^{2} - 1} (x + t)

Case 5.

(37)

\{\begin{matrix} a_{- 1} = \frac{1}{4} (a_{0}^{2} - b_{0}^{2}), & b_{- 1} = - \frac{1}{4} (a_{0}^{2} - b_{0}^{2}), \\ a_{0} = a_{0}, & b_{0} = b_{0}, \\ a_{1} = 1, & k = ω = \sqrt{2}, \end{matrix}

where

b_{0}

and

a_{0}

are free parameters. Inserting these result into Eq. (11), we obtain exact solution

u (x, t) = \frac{\exp (η) + a_{0} + \frac{1}{4} (a_{0}^{2} - b_{0}^{2}) \exp (- η)}{\exp (η) + b_{0} - \frac{1}{4} (a_{0}^{2} - b_{0}^{2}) \exp (- η)},

with

η = \sqrt{2} (x + t)

Example 2

Now we consider the Allen-Cahn equation (Wazwaz, 2007; Tascana and Bekir, 2009)

(38)

u_{t} = u_{xx} - u^{3} + u .

Using the transformation (2), Eq. (38) becomes

(39)

k^{2} U^{″} - ω U^{'} - U^{3} + U = 0,

where prime denotes the differential with respect to

η

. The highest nonlinear term

U^{3}

is now given by

(40)

U^{3} = \frac{c_{1} \exp [(3 c + p) η] + \dots}{c_{2} \exp [4 p η] + \dots},

and the highest linear term

U^{″}

is given by

(41)

U^{″} = \frac{c_{3} \exp [(c + 3 p) η] + \dots}{c_{4} \exp [4 p η] + \dots} .

Balancing the highest order of exp-function in Eqs. (41) and (40), we have

3 c + p = c + 3 p

, so

p = c

. As since illustrated in the previous example, we can also obtain that

q = d

. Here, we only consider the simplest case

p = c = 1

and

q = d = 1

, so Eq. (3) reduces to

(42)

U (η) = \frac{a_{1} \exp (η) + a_{0} + a_{- 1} \exp (- η)}{\exp (η) + b_{0} + b_{- 1} \exp (- η)} .

Substituting Eq. (42) into Eq. (39), equating to zero the coefficients of all powers of

\exp (n η)

yields a set of algebraic equations for

a_{0}, b_{0}

a_{- 1}, a_{1}, b_{1}, k

and

ω

(see Appendix B). By solving the system of algebraic equations with a professional mathematical software, we obtain

Case 1.

(43)

\{\begin{matrix} a_{1} = 1, & b_{- 1} = 0, \\ a_{0} = 0, & b_{0} = b_{0}, \\ a_{- 1} = 0, & k = \pm \frac{\sqrt{2}}{2}, & ω = \frac{3}{2} . \end{matrix}

Substituting these result into Eq. (42), we obtain

(44)

u (x, t) = \frac{1}{1 + b_{0} \exp [\frac{\sqrt{2}}{2} (\mp x - \frac{3 \sqrt{2}}{2} t)]},

where

b_{0}

is free parameter. Replacing

b_{0} = \pm 1

in Eq. (44) gives the solutions obtained in Wazwaz (2007), Tascana and Bekir (2009). Fig. 3 depict the solution Eq. (44) with

b_{0} = 1

The solutions of Eq. (44) for the parameter value b 0 = 1 .

Case 2.

(45)

\{\begin{matrix} a_{1} = - 1, & b_{- 1} = 0, \\ a_{0} = 0, & b_{0} = b_{0}, \\ a_{- 1} = 0, & k = \pm \frac{\sqrt{2}}{2}, & ω = \frac{3}{2} . \end{matrix}

Substituting these result into Eq. (42), we obtain

(46)

u (x, t) = - \frac{1}{1 + b_{0} \exp [\frac{\sqrt{2}}{2} (\mp x - \frac{3 \sqrt{2}}{2} t)]},

where

b_{0}

is free parameter. Replacing

b_{0} = \pm 1

in Eq. (46) gives the solutions obtained in Wazwaz (2007), Tascana and Bekir (2009).

Example 3

Finally, we consider the solution of steady-state equation that is presented by a ordinary differential equation in the following form (Elliott and French, 1987):

(47)

α u^{″} (x) = β u (x) (u (x) - m) (u (x) + m),

where

α, β, m

are the constants. By making the transformation

v (x) = m^{- 1} (u (ɛ x) + m),

where

ɛ = {(\frac{α}{β m^{2}})}^{1 / 2}

, Eq. (47) becomes

(48)

v^{″} (x) = v (x) (v (x) - 1) (v (x) - 2) .

According to the Exp-function method (He and Wu, 2006), we assume that the solution of Eq. (48) can be expressed in the following from:

(49)

v (x) = \frac{a_{1} \exp [kx + ω] + a_{0} + a_{- 1} \exp [- kx - ω]}{\exp [kx + ω] + b_{0} + b_{- 1} \exp [- kx - ω]},

where

a_{1}, a_{0}, a_{- 1}, b_{0}, b_{- 1}, k

and

ω

are constants which are unknown to be further determined. Substituting Eq. (49) into Eq. (48) and equating the coefficients of all powers of

\exp [i (kx + ω)] (i = 0, \pm 1, \pm 2, \dots)

to zero yields a set of algebraic equations for

a_{1}, a_{0}, a_{- 1}, b_{0}, b_{- 1}, k

and

ω

(see Appendix C). By solving the system of algebraic equations with a professional mathematical software, we obtain

Case 1.

(50)

\{\begin{matrix} a_{- 1} = 0, & b_{- 1} = b_{- 1}, \\ a_{0} = 0, & b_{0} = 0, \\ a_{1} = 2, & k = \frac{\sqrt{2}}{2}, & ω = ω . \end{matrix}

Substituting these results into (49), we obtain

(51)

v (x) = \frac{2 \exp [\frac{\sqrt{2}}{2} x + ω]}{\exp [\frac{\sqrt{2}}{2} x + ω] + b_{- 1} \exp [- \frac{\sqrt{2}}{2} x - ω]} .

By the some manipulation in Eq. (51), we obtain

v (x) = \frac{\exp [\frac{\sqrt{2}}{2} x + ω] + b_{- 1} \exp [- \frac{\sqrt{2}}{2} x - ω] + \exp [\frac{\sqrt{2}}{2} x + ω] - b_{- 1} \exp [- \frac{\sqrt{2}}{2} x - ω]}{\exp [\frac{\sqrt{2}}{2} x + ω] + b_{- 1} \exp [- \frac{\sqrt{2}}{2} x - ω]},

or equivalently

(52)

v (x) = 1 + \frac{\exp [\frac{\sqrt{2}}{2} x + ω] - b_{- 1} \exp [- \frac{\sqrt{2}}{2} x - ω]}{\exp [\frac{\sqrt{2}}{2} x + ω] + b_{- 1} \exp [- \frac{\sqrt{2}}{2} x - ω]} .

If we choose

b_{- 1} = 1, ω = - \frac{\sqrt{2}}{2} ɛ

in our solution Eq. (52) gives

v (x) = 1 + \tanh [\frac{\sqrt{2}}{2} (x - ɛ)] .

(53)

u (x) = m \tanh [\frac{\sqrt{2}}{2} (m \sqrt{\frac{β}{α}} x - ɛ)],

which this the solution gives the solution obtained in Elliott and French (1987). Fig. 4 depict the solution of Eq. (47), when

m = α = β = 1

The solutions of Eq. (53) for the parameter value m = α = β = 1 .

Case 2.

(54)

\{\begin{matrix} a_{- 1} = 0, & b_{- 1} = 0, \\ a_{0} = 2 b_{0}, & b_{0} = b_{0}, \\ a_{1} = 0, & k = \sqrt{2}, & ω = ω . \end{matrix}

Substituting these results into (49), we obtain

(55)

v (x) = \frac{2 b_{0}}{\exp [\sqrt{2} x + ω] + b_{0}} .

By the some manipulation in Eq. (55), we obtain

v (x) = \frac{b_{0} + \exp [\sqrt{2} x + ω] + b_{0} - \exp [\sqrt{2} x + ω]}{\exp [\sqrt{2} x + ω] + b_{0}} .

Or equivalently

v (x) = 1 - \frac{\exp [\sqrt{2} x + ω] - b_{0}}{\exp [\sqrt{2} x + ω] + b_{0}} .

Or equivalently

(56)

v (x) = 1 - \frac{\exp [\frac{\sqrt{2}}{2} x + \frac{ω}{2}] - b_{0} \exp [- \frac{\sqrt{2}}{2} x - \frac{ω}{2}]}{\exp [\frac{\sqrt{2}}{2} x + \frac{ω}{2}] + b_{0} \exp [- \frac{\sqrt{2}}{2} x - \frac{ω}{2}]} .

If we choose

b_{0} = 1, ω = - \sqrt{2} ɛ

in our solution Eq. (56) gives

v (x) = 1 - \tanh [\frac{\sqrt{2}}{2} (x - ɛ)] .

u (x) = - m \tanh [\frac{\sqrt{2}}{2} (m \sqrt{\frac{β}{α}} x - ɛ)],

which this the solution gives the solution obtained in Elliott and French (1987).

4 Conclusion

In this paper, the Exp-function method has been tested by applying it successfully to the Cahn–Hilliard equations, Allen–Cahn equation and Steady-State equation. In the Exp-function method, the free parameters may imply some physical meaningful results for the problem considered. The free parameters, of course, might be related to initial conditions as well. The performance of the Exp-function method is reliable and effective. Many methods (such as Adomian decomposition method (ADM), Homotopy analysis method (HAM) or variational iteration method (VIM)) can only obtain a special equation with special boundary (initial) conditions; some obtained solutions are physically meaningless. This paper obtains solutions with free parameters that can be determined via boundary (initial) conditions. In applications of Exp-function method in past decade common errors in finding exact solutions of nonlinear problems have been omitted (53). In this paper we present an application of this method with tackling these common errors. The solving procedure reveals that the exp-function method is a straightforward, succinct and promising tool for solving nonlinear partial differential equation and nonlinear ordinary differential equation.

References

Abdou M.A., Soliman A.A., El-Basyony S.T., . New application of exp-function method for improved Boussinesq equation. Physics Letters, Section A: General, Atomic and Solid State Physics. 2007;369:469-475.
[Google Scholar]
Berberler M., Yildirim A., . He’s homotopy perturbation method for solving the shock wave equation. Applicable Analysis. 2009;88:997-1004.
[Google Scholar]
Biazar J., Ayati Z., . Extension of the exp-function method for systems of two-dimensional Burger’s equations. Computers and Mathematics with Applications. 2009;58:2103-2106.
[Google Scholar]
Biazar J., Badpeimaa F., Azimi F., . Application of the homotopy perturbation method to Zakharov–Kuznetsov equations. Computers and Mathematics with Applications. 2009;58:2391-2394.
[Google Scholar]
Cahn J.W., Hilliard J.E., . Free energy of a nonuniform system. I: Interfacial free energy. The Journal of Chemical Physics. 1958;28:258-267.
[Google Scholar]
Chan J., . On spinodal decomposition. Acta Metallurgica. 1961;9:795-801.
[Google Scholar]
Choo S., Chung S.K., Lee Y., . A conservative difference scheme for the viscous Cahn–Hilliard with a nonconstant gradient energy coefficient. Applied Numerical Mathematics. 2004;51:207-219.
[Google Scholar]
Dahmani Z., Benbachir M., . Solutions of the Cahn–Hilliard equation with time-and space-fractional derivatives. International Journal of Nonlinear Science. 2009;8:19-26.
[Google Scholar]
Dehghan M., Mirzaei D., . A numerical method based on the boundary integral equation and dual reciprocity methods for one-dimensional Cahn–Hilliard equation. Engineering Analysis with Boundary Elements. 2009;33:522-528.
[Google Scholar]
Dehghan M., Shokri A., . Numerical solution of the nonlinear Klein–Gordon equation using radial basis functions. Journal of Computational and Applied Mathematics. 2009;230:400-410.
[Google Scholar]
Domairry G., Davodi A.G., Davodi A.G., . Solutions for the double sine-Gordon equation by exp-function, tanh, and extended tanh methods. Numerical Methods for Partial Differential Equations. 2010;26:384-398.
[Google Scholar]
Ebaid A., . Exact solitary wave solutions for some nonlinear evolution equations via exp-function method. Physics Letters, Section A: General, Atomic and Solid State Physics. 2007;365:213-219.
[Google Scholar]
Ebaid A., . Generalization of He’s exp-function method and new exact solutions for Burgers equation. Zeitschrift fur Naturforschung Section A: Journal of Physical Sciences. 2009;64:604-608.
[Google Scholar]
Elliott C.M., French D.A., . Numerical studies of the Cahn–Hilliard equation for phase separation. IMA Journal of Applied Mathematics. 1987;38:97-128.
[Google Scholar]
Elliott C., Zheng S., . On the Cahn–Hilliard equation. Archive for Rational Mechanics and Analysis. 1986;96:339-357.
[Google Scholar]
Furihata D., . A stable and conservative finite difference scheme for the Cahn–Hilliard equation. Numerische Mathematik. 2001;87:675-699.
[Google Scholar]
Gurtin M., . Generalised Ginzburg–Landau and Cahn–Hilliard equations based on a microforce balance. Physica D. 1996;92:178-192.
[Google Scholar]
Hashim I., Noorani M.S.M., Hadidi M.R.S., . Solving the generalized Burgers–Huxley equation using the Adomian decomposition method. Mathematical and Computer Modelling. 2006;43:1404-1411.
[Google Scholar]
He J.H., Abdou M.A., . New periodic solutions for nonlinear evolution equations using exp-function method. Chaos, Solitons and Fractals. 2007;34:1421-1429.
[Google Scholar]
He J.H., Wu X.H., . Exp-function method for nonlinear wave equations. Chaos, Solitons and Fractals. 2006;30:700-708.
[Google Scholar]
He J.H., Zhang L., . Generalized solitary solution and compacton-like solution of the Jaulent–Miodek equations using the exp-function method. Physics Letters, Section A: General, Atomic and Solid State Physics. 2008;372:1044-1047.
[Google Scholar]
Jabbari E., Peppas N.A., . A model for interdiffusion at interfaces of polymers with dissimilar physical properties. Polymer. 1995;36:575-586.
[Google Scholar]
Jingxue Y., . On the existence of nonnegative continuous solutions of the Cahn–Hilliard equation. Journal of Differential Equations. 1992;97:310-327.
[Google Scholar]
Kim J., . A numerical method for the Cahn–Hilliard equation with a variable mobility. Communications in Nonlinear Science and Numerical Simulation. 2007;12:1560-1571.
[Google Scholar]
Koroglu C., zis T., . A novel traveling wave solution for Ostrovsky equation using exp-function method. Computers and Mathematics with Applications. 2009;58:2142-2146.
[Google Scholar]
Kudryashov M.A., . Seven common errors in finding exact solutions of nonlinear differential equations. Communications in Nonlinear Science and Numerical Simulation. 2009;14:3507-3529.
[Google Scholar]
Mello E.D., Filho O.T.D.S., . Numerical study of the Cahn–Hilliard equation in one, two and three dimensions. Physica A: Statistical Mechanics and its Applications. 2005;347:429-443.
[Google Scholar]
Ozis T., Koroglu C., . A novel approach for solving the Fisher equation using exp-function method. Physics Letters, Section A: General, Atomic and Solid State Physics. 2008;372:3836-3840.
[Google Scholar]
Parand K., Taghavi A., . Rational scaled generalized Laguerre function collocation method for solving the Blasius equation. Journal of Computational and Applied Mathematics. 2009;233:980-989.
[Google Scholar]
Parand K., Dehghan M., Pirkhedri A., . Sinc-collocation method for solving the Blasius equation. Physics Letters, Section A: General, Atomic and Solid State Physics. 2009;373:4060-4065.
[Google Scholar]
Parand, K., Dehghan, M., Rezaei, A.R., Ghaderi, S.M., An approximation algorithm for the solution of the nonlinear Lane–Emden type equations arising in astrophysics using hermite functions collocation method, Computer Physics Communications. doi:10.1016/j.cpc.2010.02.018.
Puri S., Binder K., . Phenomenological theory for the formation of interfaces via the interdiffusion of layers. Physical Review B. 1991;44:9735-9738.
[Google Scholar]
Rashidi M., Ganji D., Dinarvand S., . Explicit analytical solutions of the generalized Burger and Burger–Fisher equations by homotopy perturbation method. Numerical Methods for Partial Differential Equations. 2009;25:409-417.
[Google Scholar]
Shakeri F., Dehghan M., . Numerical solution of the Klein–Gordon equation via he’s variational iteration method. Nonlinear Dynamics. 2008;51:89-97.
[Google Scholar]
Shin B., Darvishi M., Barati A., . Some exact and new solutions of the Nizhnik–Novikov–Vesselov equation using the exp-function method. Computers and Mathematics with Applications. 2009;58:2147-2151.
[Google Scholar]
Soliman A., Abdou M., . Numerical solutions of nonlinear evolution equations using variational iteration method. Journal of Computational and Applied Mathematics. 2007;207:111-120.
[Google Scholar]
Tascana F., Bekir A., . Travelling wave solutions of the Cahn–Allen equation by using first integral method. Applied Mathematics and Computation. 2009;207:279-282.
[Google Scholar]
Tascan F., Bekir A., . Analytic solutions of the (2 + 1)-dimensional nonlinear evolution equations using the sine–cosine method. Applied Mathematics and Computation. 2009;215:3134-3139.
[Google Scholar]
Tatari M., Dehghan M.M., . A method for solving partial differential equations via radial basis functions: application to the heat equation. Engineering Analysis with Boundary Elements. 2010;34:206-212.
[Google Scholar]
Tatari M., Dehghan M., Razzaghi M., . Application of the Adomian decomposition method for the Fokker–Planck equation. Mathematical and Computer Modelling. 2007;45:639-650.
[Google Scholar]
Ugurlu Y., Kaya D., . Solutions of the Cahn–Hilliard equation. Computers and Mathematics with Applications. 2008;56:3038-3045.
[Google Scholar]
Wang S.Q., Shi Q., . Interdiffusion in binary polymer mixtures. Macromolecules. 1993;26:1091-1096.
[Google Scholar]
Wazwaz A., . The tanh method: exact solutions of the sine-Gordon and the sinh-Gordon equations. Applied Mathematics and Computation. 2005;167:1196-1210.
[Google Scholar]
Wazwaz A., . The variable separated ode and the tanh methods for solving the combined and the double combined sinh–cosh-Gordon equations. Applied Mathematics and Computation. 2006;177:745-754.
[Google Scholar]
Wazwaz A., . Exact solutions for the generalized sine-Gordon and the generalized sinh-Gordon equations. Chaos, Solitons and Fractals. 2006;28:127-135.
[Google Scholar]
Wazwaz A.M., . New solitary wave solutions to the modified Kawahara equation. Physics Letters, Section A: General, Atomic and Solid State Physics. 2007;360:588-592.
[Google Scholar]
Wazwaz A.M., . The tanh–coth method for solitons and kink solutions for nonlinear parabolic equations. Applied Mathematics and Computation. 2007;188:1467-1475.
[Google Scholar]
Wells G.N., Kuhl E., Garikipati K., . A discontinuous Galerkin method for the Cahn–Hilliard equation. Journal of Computational Physics. 2006;218:860-877.
[Google Scholar]
Yusufoglu E., Bekir A., . The variational iteration method for solitary patterns solutions of gbbm equation. Physics Letters A. 2007;367:461-464.
[Google Scholar]
Zhang S., . Application of exp-function method to a KdV equation with variable coefficients. Physics Letters, Section A: General, Atomic and Solid State Physics. 2007;365:448-453.
[Google Scholar]
Zhang S., . Exp-function method for solving Maccari’s system. Physics Letters, Section A: General, Atomic and Solid State Physics. 2007;371:65-71.
[Google Scholar]
Zhang S., . Application of exp-function method to high-dimensional nonlinear evolution equation. Chaos, Solitons and Fractals. 2008;38:270-276.
[Google Scholar]

Appendix A

$- ω a_{- 1} b_{0} b_{- 1}^{3} - k^{4} a_{- 1} b_{0} b_{- 1}^{3} - 3 k^{2} a_{- 1}^{2} a_{0} b_{- 1}^{2} + 3 k^{2} a_{- 1}^{3} b_{0} b_{- 1} - k^{2} a_{- 1} b_{0} b_{- 1}^{3} + {ka}_{- 1} b_{0} b_{- 1}^{3} + ω a_{0} b_{- 1}^{4} + k^{4} a_{0} b_{- 1}^{4} + k^{2} a_{0} b_{- 1}^{4} - {ka}_{0} b_{- 1}^{4} = 0,$ $- 3 {ka}_{0} b_{- 1}^{3} b_{0} - 3 ω a_{- 1} b_{0}^{2} b_{- 1}^{2} + 3 ω a_{0} b_{- 1}^{3} b_{0} + 11 k^{4} a_{- 1} b_{0}^{2} b_{- 1}^{2} - 11 k^{4} a_{0} b_{- 1}^{3} b_{0} - 12 k^{2} a_{- 1} a_{0}^{2} b_{- 1}^{2} - k^{2} a_{- 1} b_{0}^{2} b_{- 1}^{2} + k^{2} a_{0} b_{- 1}^{3} b_{0} + 3 {ka}_{- 1} b_{0}^{2} b_{- 1}^{2} - 16 k^{4} a_{- 1} b_{- 1}^{3} + 12 k^{2} a_{- 1}^{3} b_{- 1} + 4 k^{2} a_{1} b_{- 1}^{4} - 4 k^{2} a_{- 1} b_{- 1}^{3} - 2 {ka}_{1} b_{- 1}^{4} - 9 k^{2} a_{- 1}^{3} b_{0}^{2} - 2 ω a_{- 1} b_{- 1}^{3} + 16 k^{4} a_{1} b_{- 1}^{4} + 2 ω a_{1} b_{- 1}^{4} + 2 {ka}_{- 1} b_{- 1}^{3} - 12 k^{2} a_{- 1}^{2} a_{1} b_{- 1}^{2} + 21 k^{2} a_{- 1}^{2} b_{0} a_{0} b_{- 1} = 0,$ $2 ω a_{0} b_{- 1}^{3} - 2 {ka}_{0} b_{- 1}^{3} - 76 k^{4} a_{0} b_{- 1}^{3} - 4 k^{2} a_{0} b_{- 1}^{3} + 77 k^{4} a_{- 1} b_{0} b_{- 1}^{2} + 11 k^{2} a_{1} b_{0} b_{- 1}^{3} - 7 k^{2} a_{- 1} b_{0} b_{- 1}^{2} - 7 {ka}_{1} b_{0} b_{- 1}^{3} + 9 {ka}_{- 1} b_{0} b_{- 1}^{2} - 9 ω a_{- 1} b_{0} b_{- 1}^{2} + 7 ω a_{1} b_{0} b_{- 1}^{3} - 9 k^{2} a_{0}^{3} b_{- 1}^{2} - 33 k^{2} a_{- 1}^{3} b_{0} + 21 k^{2} a_{- 1} b_{0} a_{0}^{2} b_{- 1} - 54 k^{2} a_{- 1} a_{1} b_{- 1}^{2} a_{0} + 21 k^{2} a_{- 1}^{2} a_{1} b_{- 1} b_{0} - k^{4} a_{1} b_{0} b_{- 1}^{3} - 3 {ka}_{0} b_{- 1}^{2} b_{0}^{2} - 3 ω a_{- 1} b_{0}^{3} b_{- 1} + 3 ω a_{0} b_{- 1}^{2} b_{0}^{2} - 11 k^{4} a_{- 1} b_{0}^{3} b_{- 1} + 11 k^{4} a_{0} b_{- 1}^{2} b_{0}^{2} - 12 k^{2} a_{0} a_{- 1}^{2} b_{0}^{2} + 66 k^{2} a_{- 1}^{2} a_{0} b_{- 1} + k^{2} a_{- 1} b_{0}^{3} b_{- 1} - k^{2} a_{0} b_{- 1}^{2} b_{0}^{2} + 3 {ka}_{- 1} b_{0}^{3} b_{- 1} = 0,$ $- 36 k^{2} a_{- 1}^{3} + 18 k^{2} a_{0} a_{1} b_{- 1} a_{- 1} b_{0} - 51 k^{2} a_{- 1}^{2} b_{0} a_{0} + 84 k^{2} a_{- 1}^{2} a_{1} b_{- 1} - 3 k^{2} a_{1} a_{- 1}^{2} b_{0}^{2} + 11 k^{2} a_{1} b_{0}^{2} b_{- 1}^{2} - {ka}_{0} b_{- 1} b_{0}^{3} - 58 k^{4} a_{- 1} b_{0}^{2} b_{- 1} + 47 k^{4} a_{0} b_{- 1}^{2} b_{0} - 6 {ka}_{1} b_{- 1}^{3} - 4 k^{2} a_{- 1} b_{- 1}^{2} + 6 {ka}_{- 1} b_{- 1}^{2} + 4 k^{2} a_{1} b_{- 1}^{3} - 176 k^{4} a_{1} b_{- 1}^{3} - 6 ω a_{- 1} b_{- 1}^{2} + 176 k^{4} a_{- 1} b_{- 1}^{2} + 6 ω a_{1} b_{- 1}^{3} - ω a_{- 1} b_{0}^{4} + k^{4} a_{- 1} b_{0}^{4} + 3 k^{2} a_{0}^{3} b_{- 1} b_{0} + 2 k^{2} a_{- 1} b_{0}^{2} b_{- 1} - 13 k^{2} a_{0} b_{- 1}^{2} b_{0} - 9 {ka}_{1} b_{0}^{2} b_{- 1}^{2} + 12 {ka}_{- 1} b_{0}^{2} b_{- 1} - 3 {ka}_{0} b_{- 1}^{2} b_{0} + 9 ω a_{1} b_{0}^{2} b_{- 1}^{2} - 12 ω a_{- 1} b_{0}^{2} b_{- 1} + 3 ω a_{0} b_{- 1}^{2} b_{0} + 11 k^{4} a_{1} b_{0}^{2} b_{- 1}^{2} + ω a_{0} b_{- 1} b_{0}^{3} - k^{4} a_{0} b_{- 1} b_{0}^{3} - 48 k^{2} a_{1} b_{- 1}^{2} a_{0}^{2} + 84 k^{2} a_{- 1} a_{0}^{2} b_{- 1} - 48 k^{2} a_{- 1} a_{1}^{2} b_{- 1}^{2} - 3 k^{2} a_{0}^{2} a_{- 1} b_{0}^{2} - k^{2} a_{0} b_{- 1} b_{0}^{3} + k^{2} a_{- 1} b_{0}^{4} + {ka}_{- 1} b_{0}^{4} = 0,$ $15 {ka}_{- 1} b_{0} b_{- 1} + 5 ω a_{1} b_{0}^{3} b_{- 1} + 15 ω a_{1} b_{0} b_{- 1}^{2} - 15 {ka}_{1} b_{0} b_{- 1}^{2} - 10 k^{4} a_{0} b_{- 1} b_{0}^{2} + 5 k^{2} a_{1} b_{0}^{3} b_{- 1} + 5 k^{2} a_{1} b_{0} b_{- 1}^{2} + 5 k^{2} a_{- 1} b_{0} b_{- 1} - 10 k^{2} a_{0} b_{- 1} b_{0}^{2} - 5 {ka}_{1} b_{0}^{3} b_{- 1} - 115 k^{4} a_{- 1} b_{0} b_{- 1} - 75 k^{2} a_{0} a_{1}^{2} b_{- 1}^{2} - 5 ω a_{- 1} b_{0}^{3} + 5 k^{4} a_{- 1} b_{0}^{3} + 230 k^{4} a_{0} b_{- 1}^{2} + 30 k^{2} a_{0}^{3} b_{- 1} - 75 k^{2} a_{0} a_{- 1}^{2} - 15 k^{2} a_{1} a_{- 1}^{2} b_{0} - 15 k^{2} a_{- 1} a_{0}^{2} b_{0} + 180 k^{2} a_{0} a_{1} b_{- 1} a_{- 1} - 15 k^{2} a_{1}^{2} b_{- 1} a_{- 1} b_{0} - 15 k^{2} a_{1} b_{- 1} a_{0}^{2} b_{0} + 5 k^{2} a_{- 1} b_{0}^{3} + 5 {ka}_{- 1} b_{0}^{3} - 10 k^{2} a_{0} b_{- 1}^{2} - 15 ω a_{- 1} b_{0} b_{- 1} + 5 k^{4} a_{1} b_{0}^{3} b_{- 1} - 115 k^{4} a_{1} b_{0} b_{- 1}^{2} = 0,$ $47 k^{4} a_{0} b_{- 1} b_{0} - 3 k^{2} a_{1}^{2} b_{0}^{2} a_{- 1} + 2 k^{2} a_{1} b_{0}^{2} b_{- 1} - 13 k^{2} a_{0} b_{- 1} b_{0} - 12 {ka}_{1} b_{0}^{2} b_{- 1} + 3 {ka}_{0} b_{- 1} b_{0} + 12 ω a_{1} b_{0}^{2} b_{- 1} - 3 ω a_{0} b_{- 1} b_{0} - 58 k^{4} a_{1} b_{0}^{2} b_{- 1} + 84 k^{2} a_{1}^{2} b_{- 1} a_{- 1} + 84 k^{2} a_{1} b_{- 1} a_{0}^{2} + 9 {ka}_{- 1} b_{0}^{2} + ω a_{1} b_{0}^{4} - 176 k^{4} a_{- 1} b_{- 1} - 9 ω a_{- 1} b_{0}^{2} - 6 {ka}_{1} b_{- 1}^{2} + 11 k^{2} a_{- 1} b_{0}^{2} - 4 k^{2} a_{1} b_{- 1}^{2} - 51 k^{2} a_{1}^{2} b_{0} a_{0} b_{- 1} + 18 k^{2} a_{0} a_{1} b_{0} a_{- 1} - 3 k^{2} a_{0}^{2} a_{1} b_{0}^{2} + 4 k^{2} a_{- 1} b_{- 1} + 6 {ka}_{- 1} b_{- 1} + 176 k^{4} a_{1} b_{- 1}^{2} + 11 k^{4} a_{- 1} b_{0}^{2} - ω a_{0} b_{0}^{3} - k^{4} a_{0} b_{0}^{3} - 6 ω a_{- 1} b_{- 1} + 6 ω a_{1} b_{- 1}^{2} + k^{4} a_{1} b_{0}^{4} - 36 k^{2} a_{1}^{3} b_{- 1}^{2} + 3 k^{2} a_{0}^{3} b_{0} - 48 k^{2} a_{1} a_{- 1}^{2} - 48 k^{2} a_{- 1} a_{0}^{2} - k^{2} a_{0} b_{0}^{3} - {ka}_{1} b_{0}^{4} + {ka}_{0} b_{0}^{3} + k^{2} a_{1} b_{0}^{4} = 0,$ $- 76 k^{4} a_{0} b_{- 1} + 3 {ka}_{0} b_{0}^{2} - 3 ω a_{0} b_{0}^{2} - 7 ω a_{- 1} b_{0} + 11 k^{2} a_{- 1} b_{0} - k^{4} a_{- 1} b_{0} - 11 k^{4} a_{1} b_{0}^{3} + 3 ω a_{1} b_{0}^{3} + 7 {ka}_{- 1} b_{0} + 11 k^{4} a_{0} b_{0}^{2} + k^{2} a_{1} b_{0}^{3} - 9 {ka}_{1} b_{0} b_{- 1} - 7 k^{2} a_{1} b_{0} b_{- 1} + 9 ω a_{1} b_{0} b_{- 1} + 77 k^{4} a_{1} b_{0} b_{- 1} - 33 k^{2} a_{1}^{3} b_{0} b_{- 1} + 21 k^{2} a_{1}^{2} b_{0} a_{- 1} - 54 k^{2} a_{1} a_{- 1} a_{0} + 21 k^{2} a_{1} b_{0} a_{0}^{2} - 12 k^{2} a_{0} a_{1}^{2} b_{0}^{2} + 66 k^{2} a_{1}^{2} b_{- 1} a_{0} - 2 ω a_{0} b_{- 1} - 3 {ka}_{1} b_{0}^{3} + 2 {ka}_{0} b_{- 1} - k^{2} a_{0} b_{0}^{2} - 4 k^{2} a_{0} b_{- 1} - 9 k^{2} a_{0}^{3} = 0,$ $21 k^{2} a_{1}^{2} b_{0} a_{0} - 12 k^{2} a_{1}^{2} a_{- 1} + 12 k^{2} a_{1}^{3} b_{- 1} - 4 k^{2} a_{1} b_{- 1} - 2 {ka}_{1} b_{- 1} + 3 ω a_{1} b_{0}^{2} - 3 ω a_{0} b_{0} + 11 k^{4} a_{1} b_{0}^{2} - 11 k^{4} a_{0} b_{0} - 12 k^{2} a_{0}^{2} a_{1} - 9 k^{2} a_{1}^{3} b_{0}^{2} - k^{2} a_{1} b_{0}^{2} + k^{2} a_{0} b_{0} - 3 {ka}_{1} b_{0}^{2} + 3 {ka}_{0} b_{0} - 16 k^{4} a_{1} b_{- 1} + 2 ω a_{1} b_{- 1} - 2 ω a_{- 1} + 2 {ka}_{- 1} + 16 k^{4} a_{- 1} + 4 k^{2} a_{- 1} = 0,$ $k^{2} a_{0} - ω a_{0} + k^{4} a_{0} + {ka}_{0} + ω a_{1} b_{0} - k^{4} a_{1} b_{0} + 3 k^{2} a_{1}^{3} b_{0} - 3 k^{2} a_{0} a_{1}^{2} - k^{2} a_{1} b_{0} - {ka}_{1} b_{0} = 0 .$

Appendix B

$a_{- 1} b_{- 1}^{2} - a_{- 1}^{3} = 0, 2 a_{- 1} b_{0} b_{- 1} + k^{2} a_{0} b_{- 1}^{2} - ω a_{0} b_{- 1}^{2} - 3 a_{0} a_{- 1}^{2} + a_{0} b_{- 1}^{2} - k^{2} a_{- 1} b_{0} b_{- 1} + ω a_{- 1} b_{0} b_{- 1} = 0, 2 a_{0} b_{- 1} b_{0} + k^{2} a_{- 1} b_{0}^{2} + ω a_{- 1} b_{0}^{2} + 4 k^{2} a_{1} b_{- 1}^{2} - 4 k^{2} a_{- 1} b_{- 1} - 2 ω a_{1} b_{- 1}^{2} + 2 ω a_{- 1} b_{- 1} - 3 a_{0}^{2} a_{- 1} + a_{- 1} b_{0}^{2} - 3 a_{1} a_{- 1}^{2} + 2 a_{- 1} b_{- 1} + a_{1} b_{- 1}^{2} - k^{2} a_{0} b_{- 1} b_{0} - ω a_{0} b_{- 1} b_{0} = 0, 2 a_{1} b_{0} b_{- 1} - 6 k^{2} a_{0} b_{- 1} + 3 k^{2} a_{- 1} b_{0} + 3 ω a_{- 1} b_{0} - 6 a_{1} a_{0} a_{- 1} + 3 k^{2} a_{1} b_{0} b_{- 1} - 3 ω a_{1} b_{0} b_{- 1} + a_{0} b_{0}^{2} + 2 a_{0} b_{- 1} + 2 a_{- 1} b_{0} - a_{0}^{3} = 0, k^{2} a_{1} b_{0}^{2} - k^{2} a_{0} b_{0} - ω a_{1} b_{0}^{2} + ω a_{0} b_{0} - 2 ω a_{1} b_{- 1} - 4 k^{2} a_{1} b_{- 1} + a_{1} b_{0}^{2} + 2 a_{1} b_{- 1} + 2 ω a_{- 1} - 3 a_{1} a_{0}^{2} + 2 a_{0} b_{0} - 3 a_{1}^{2} a_{- 1} + a_{- 1} + 4 k^{2} a_{- 1} = 0, - k^{2} a_{1} b_{0} - ω a_{1} b_{0} + ω a_{0} + k^{2} a_{0} - 3 a_{1}^{2} a_{0} + 2 a_{1} b_{0} + a_{0} = 0, - a_{1}^{3} + a_{1} = 0 .$

Appendix C

$- a_{- 1}^{3} - 2 a_{- 1} b_{- 1}^{2} + 3 a_{- 1}^{2} b_{- 1} = 0, - k^{2} a_{- 1} b_{0} b_{- 1} + k^{2} a_{0} b_{- 1}^{2} + 6 a_{0} a_{- 1} b_{- 1} - 4 a_{- 1} b_{0} b_{- 1} + 3 a_{- 1}^{2} b_{0} - 3 a_{0} a_{- 1}^{2} - 2 a_{0} b_{- 1}^{2} = 0, - 4 a_{0} b_{- 1} b_{0} + 6 a_{0} a_{- 1} b_{0} + 3 a_{- 1}^{2} + 4 k^{2} a_{1} b_{- 1}^{2} - 4 k^{2} a_{- 1} b_{- 1} + 6 a_{1} a_{- 1} b_{- 1} + k^{2} a_{- 1} b_{0}^{2} - 2 a_{- 1} b_{0}^{2} - 3 a_{0}^{2} a_{- 1} + 3 a_{0}^{2} b_{- 1} - k^{2} a_{0} b_{- 1} b_{0} - 3 a_{1} a_{- 1}^{2} - 2 a_{1} b_{- 1}^{2} - 4 a_{- 1} b_{- 1} = 0, - a_{0}^{3} - 4 a_{1} b_{0} b_{- 1} + 3 k^{2} a_{- 1} b_{0} - 6 k^{2} a_{0} b_{- 1} - 6 a_{1} a_{0} a_{- 1} + 6 a_{1} a_{0} b_{- 1} + 6 a_{1} a_{- 1} b_{0} + 6 a_{0} a_{- 1} - 4 a_{- 1} b_{0} + 3 a_{0}^{2} b_{0} - 2 a_{0} b_{0}^{2} - 4 a_{0} b_{- 1} + 3 k^{2} a_{1} b_{0} b_{- 1} = 0, 3 a_{0}^{2} - k^{2} a_{0} b_{0} + 6 a_{1} a_{0} b_{0} - 2 a_{- 1} - 4 k^{2} a_{1} b_{- 1} + k^{2} a_{1} b_{0}^{2} - 2 a_{1} b_{0}^{2} - 4 a_{0} b_{0} - 3 a_{1} a_{0}^{2} + 3 a_{1}^{2} b_{- 1} + 6 a_{1} a_{- 1} - 4 a_{1} b_{- 1} + 4 a_{- 1} k^{2} - 3 a_{1}^{2} a_{- 1} = 0, - 2 a_{0} - k^{2} a_{1} b_{0} - 4 a_{1} b_{0} + k^{2} a_{0} - 3 a_{1}^{2} a_{0} + 3 a_{1}^{2} b_{0} + 6 a_{1} a_{0} = 0, 3 a_{1}^{2} - 2 a_{1} - a_{1}^{3} = 0 .$

Introduction

Basic idea of Exp-function method

Applications of the Exp-function method

Conclusion

Fulltext Views
14

PDF downloads
7

View/Download PDF
Download Citations

BibTeX
RIS

Show Sections

[1] Abdou M.A., Soliman A.A., El-Basyony S.T., . New application of exp-function method for improved Boussinesq equation. Physics Letters, Section A: General, Atomic and Solid State Physics. 2007;369:469-475.
[Google Scholar]

[2] Berberler M., Yildirim A., . He’s homotopy perturbation method for solving the shock wave equation. Applicable Analysis. 2009;88:997-1004.
[Google Scholar]

[3] Biazar J., Ayati Z., . Extension of the exp-function method for systems of two-dimensional Burger’s equations. Computers and Mathematics with Applications. 2009;58:2103-2106.
[Google Scholar]

[4] Biazar J., Badpeimaa F., Azimi F., . Application of the homotopy perturbation method to Zakharov–Kuznetsov equations. Computers and Mathematics with Applications. 2009;58:2391-2394.
[Google Scholar]

[5] Cahn J.W., Hilliard J.E., . Free energy of a nonuniform system. I: Interfacial free energy. The Journal of Chemical Physics. 1958;28:258-267.
[Google Scholar]

[6] Chan J., . On spinodal decomposition. Acta Metallurgica. 1961;9:795-801.
[Google Scholar]

[7] Choo S., Chung S.K., Lee Y., . A conservative difference scheme for the viscous Cahn–Hilliard with a nonconstant gradient energy coefficient. Applied Numerical Mathematics. 2004;51:207-219.
[Google Scholar]

[8] Dahmani Z., Benbachir M., . Solutions of the Cahn–Hilliard equation with time-and space-fractional derivatives. International Journal of Nonlinear Science. 2009;8:19-26.
[Google Scholar]

[9] Dehghan M., Mirzaei D., . A numerical method based on the boundary integral equation and dual reciprocity methods for one-dimensional Cahn–Hilliard equation. Engineering Analysis with Boundary Elements. 2009;33:522-528.
[Google Scholar]

[10] Dehghan M., Shokri A., . Numerical solution of the nonlinear Klein–Gordon equation using radial basis functions. Journal of Computational and Applied Mathematics. 2009;230:400-410.
[Google Scholar]

[11] Domairry G., Davodi A.G., Davodi A.G., . Solutions for the double sine-Gordon equation by exp-function, tanh, and extended tanh methods. Numerical Methods for Partial Differential Equations. 2010;26:384-398.
[Google Scholar]

[12] Ebaid A., . Exact solitary wave solutions for some nonlinear evolution equations via exp-function method. Physics Letters, Section A: General, Atomic and Solid State Physics. 2007;365:213-219.
[Google Scholar]

[13] Ebaid A., . Generalization of He’s exp-function method and new exact solutions for Burgers equation. Zeitschrift fur Naturforschung Section A: Journal of Physical Sciences. 2009;64:604-608.
[Google Scholar]

[14] Elliott C.M., French D.A., . Numerical studies of the Cahn–Hilliard equation for phase separation. IMA Journal of Applied Mathematics. 1987;38:97-128.
[Google Scholar]

[15] Elliott C., Zheng S., . On the Cahn–Hilliard equation. Archive for Rational Mechanics and Analysis. 1986;96:339-357.
[Google Scholar]

[16] Furihata D., . A stable and conservative finite difference scheme for the Cahn–Hilliard equation. Numerische Mathematik. 2001;87:675-699.
[Google Scholar]

[17] Gurtin M., . Generalised Ginzburg–Landau and Cahn–Hilliard equations based on a microforce balance. Physica D. 1996;92:178-192.
[Google Scholar]

[18] Hashim I., Noorani M.S.M., Hadidi M.R.S., . Solving the generalized Burgers–Huxley equation using the Adomian decomposition method. Mathematical and Computer Modelling. 2006;43:1404-1411.
[Google Scholar]

[19] He J.H., Abdou M.A., . New periodic solutions for nonlinear evolution equations using exp-function method. Chaos, Solitons and Fractals. 2007;34:1421-1429.
[Google Scholar]

[20] He J.H., Wu X.H., . Exp-function method for nonlinear wave equations. Chaos, Solitons and Fractals. 2006;30:700-708.
[Google Scholar]

[21] He J.H., Zhang L., . Generalized solitary solution and compacton-like solution of the Jaulent–Miodek equations using the exp-function method. Physics Letters, Section A: General, Atomic and Solid State Physics. 2008;372:1044-1047.
[Google Scholar]

[22] Jabbari E., Peppas N.A., . A model for interdiffusion at interfaces of polymers with dissimilar physical properties. Polymer. 1995;36:575-586.
[Google Scholar]

[23] Jingxue Y., . On the existence of nonnegative continuous solutions of the Cahn–Hilliard equation. Journal of Differential Equations. 1992;97:310-327.
[Google Scholar]

[24] Kim J., . A numerical method for the Cahn–Hilliard equation with a variable mobility. Communications in Nonlinear Science and Numerical Simulation. 2007;12:1560-1571.
[Google Scholar]

[25] Koroglu C., zis T., . A novel traveling wave solution for Ostrovsky equation using exp-function method. Computers and Mathematics with Applications. 2009;58:2142-2146.
[Google Scholar]

[26] Kudryashov M.A., . Seven common errors in finding exact solutions of nonlinear differential equations. Communications in Nonlinear Science and Numerical Simulation. 2009;14:3507-3529.
[Google Scholar]

[27] Mello E.D., Filho O.T.D.S., . Numerical study of the Cahn–Hilliard equation in one, two and three dimensions. Physica A: Statistical Mechanics and its Applications. 2005;347:429-443.
[Google Scholar]

[28] Ozis T., Koroglu C., . A novel approach for solving the Fisher equation using exp-function method. Physics Letters, Section A: General, Atomic and Solid State Physics. 2008;372:3836-3840.
[Google Scholar]

[29] Parand K., Taghavi A., . Rational scaled generalized Laguerre function collocation method for solving the Blasius equation. Journal of Computational and Applied Mathematics. 2009;233:980-989.
[Google Scholar]

[30] Parand K., Dehghan M., Pirkhedri A., . Sinc-collocation method for solving the Blasius equation. Physics Letters, Section A: General, Atomic and Solid State Physics. 2009;373:4060-4065.
[Google Scholar]

[31] Parand, K., Dehghan, M., Rezaei, A.R., Ghaderi, S.M., An approximation algorithm for the solution of the nonlinear Lane–Emden type equations arising in astrophysics using hermite functions collocation method, Computer Physics Communications. doi:10.1016/j.cpc.2010.02.018.

[32] Puri S., Binder K., . Phenomenological theory for the formation of interfaces via the interdiffusion of layers. Physical Review B. 1991;44:9735-9738.
[Google Scholar]

[33] Rashidi M., Ganji D., Dinarvand S., . Explicit analytical solutions of the generalized Burger and Burger–Fisher equations by homotopy perturbation method. Numerical Methods for Partial Differential Equations. 2009;25:409-417.
[Google Scholar]

[34] Shakeri F., Dehghan M., . Numerical solution of the Klein–Gordon equation via he’s variational iteration method. Nonlinear Dynamics. 2008;51:89-97.
[Google Scholar]

[35] Shin B., Darvishi M., Barati A., . Some exact and new solutions of the Nizhnik–Novikov–Vesselov equation using the exp-function method. Computers and Mathematics with Applications. 2009;58:2147-2151.
[Google Scholar]

[36] Soliman A., Abdou M., . Numerical solutions of nonlinear evolution equations using variational iteration method. Journal of Computational and Applied Mathematics. 2007;207:111-120.
[Google Scholar]

[37] Tascana F., Bekir A., . Travelling wave solutions of the Cahn–Allen equation by using first integral method. Applied Mathematics and Computation. 2009;207:279-282.
[Google Scholar]

[38] Tascan F., Bekir A., . Analytic solutions of the (2 + 1)-dimensional nonlinear evolution equations using the sine–cosine method. Applied Mathematics and Computation. 2009;215:3134-3139.
[Google Scholar]

[39] Tatari M., Dehghan M.M., . A method for solving partial differential equations via radial basis functions: application to the heat equation. Engineering Analysis with Boundary Elements. 2010;34:206-212.
[Google Scholar]

[40] Tatari M., Dehghan M., Razzaghi M., . Application of the Adomian decomposition method for the Fokker–Planck equation. Mathematical and Computer Modelling. 2007;45:639-650.
[Google Scholar]

[41] Ugurlu Y., Kaya D., . Solutions of the Cahn–Hilliard equation. Computers and Mathematics with Applications. 2008;56:3038-3045.
[Google Scholar]

[42] Wang S.Q., Shi Q., . Interdiffusion in binary polymer mixtures. Macromolecules. 1993;26:1091-1096.
[Google Scholar]

[43] Wazwaz A., . The tanh method: exact solutions of the sine-Gordon and the sinh-Gordon equations. Applied Mathematics and Computation. 2005;167:1196-1210.
[Google Scholar]

[44] Wazwaz A., . The variable separated ode and the tanh methods for solving the combined and the double combined sinh–cosh-Gordon equations. Applied Mathematics and Computation. 2006;177:745-754.
[Google Scholar]

[45] Wazwaz A., . Exact solutions for the generalized sine-Gordon and the generalized sinh-Gordon equations. Chaos, Solitons and Fractals. 2006;28:127-135.
[Google Scholar]

[46] Wazwaz A.M., . New solitary wave solutions to the modified Kawahara equation. Physics Letters, Section A: General, Atomic and Solid State Physics. 2007;360:588-592.
[Google Scholar]

[47] Wazwaz A.M., . The tanh–coth method for solitons and kink solutions for nonlinear parabolic equations. Applied Mathematics and Computation. 2007;188:1467-1475.
[Google Scholar]

[48] Wells G.N., Kuhl E., Garikipati K., . A discontinuous Galerkin method for the Cahn–Hilliard equation. Journal of Computational Physics. 2006;218:860-877.
[Google Scholar]

[49] Yusufoglu E., Bekir A., . The variational iteration method for solitary patterns solutions of gbbm equation. Physics Letters A. 2007;367:461-464.
[Google Scholar]

[50] Zhang S., . Application of exp-function method to a KdV equation with variable coefficients. Physics Letters, Section A: General, Atomic and Solid State Physics. 2007;365:448-453.
[Google Scholar]

[51] Zhang S., . Exp-function method for solving Maccari’s system. Physics Letters, Section A: General, Atomic and Solid State Physics. 2007;371:65-71.
[Google Scholar]

[52] Zhang S., . Application of exp-function method to high-dimensional nonlinear evolution equation. Chaos, Solitons and Fractals. 2008;38:270-276.
[Google Scholar]