Journal of Applied Mathematics

Volume 2012 (2012), Article ID 427479, 30 pages

doi:10.1155/2012/427479

Research Article

The Modified Rational Jacobi Elliptic Functions Method for Nonlinear Differential Difference Equations

Khaled A. Gepreel,1,2 Taher A. Nofal,2,3 and Ali A. Al-Thobaiti²

¹Mathematics Department, Faculty of Science, Zagazig University, Zagazig, Egypt

²Mathematics Department, Faculty of Science, Taif University, Saudi Arabia

³Mathematics Department, Faculty of Science, El-Minia University, El-Minia, Egypt

Received 21 March 2012; Accepted 30 May 2012

Academic Editor: Jianke Yang

Copyright © 2012 Khaled A. Gepreel et al. 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.

Abstract.

We modified the rational Jacobi elliptic functions method to construct some new exact solutions for nonlinear differential difference equations in mathematical physics via the lattice equation, the discrete nonlinear Schrodinger equation with a saturable nonlinearity, the discrete nonlinear Klein-Gordon equation, and the quintic discrete nonlinear Schrodinger equation. Some new types of the Jacobi elliptic solutions are obtained for some nonlinear differential difference equations in mathematical physics. The proposed method is more effective and powerful to obtain the exact solutions for nonlinear differential difference equations.

Introduction

In recent years, the study of difference equations has acquired a new significance, due in large part to their use in the formulation and analysis of discrete-time systems, the numerical integration of differential equation by finite-difference schemes and the study of deterministic chaos.

Since the work of Fermi et al. in the 1960s [1], DDEs have been the focus of many nonlinear studies. On the other hand, a considerable number of well-known analytic methods are successfully extended to nonlinear DDEs by researchers [2–17]. However, no method obeys the strength and the flexibility for finding all solutions to all types of nonlinear DDEs. Zhang et al. [18] and Aslan [19] used the (-expansion function method to some physically important nonlinear DDEs. Qiong and Bin [12] constructed the Jacobi elliptic solutions for nonlinear DDEs. Recently, Zhang [20] and Gepreel [21] and Gepreel and Shehata [22] have used the Jacobi elliptic function method for constructing new and more general Jacobi elliptic function solutions of the nonlinear differential difference equations. The main objective of this paper is to modify the rational Jacobi elliptic functions method to obtain the exact wave solutions for nonlinear DDEs. We use this method to calculate the exact wave solutions for some nonlinear DDEs in mathematical physics via the lattice equation, the discrete nonlinear Schrodinger equation with a saturable nonlinearity, the discrete nonlinear Klein-Gordon equation, and the quintic discrete nonlinear Schrodinger equation.

2. Description of the Modified Rational Jacobi Elliptic Functions Method

In this section, we would like to outline an algorithm for using the modified rational Jacobi elliptic functions method to solve nonlinear DDEs. For a given nonlinear DDEs,

where , and are integers, , denotes the set of all rth order derivatives of , with respect to .

The main steps of the algorithm for the modified rational Jacobi elliptic functions method to nonlinear DDEs are outlined as follows.

Step 1. We seek the traveling wave solutions of the following form:

where

where , and the phase are constants to be determined later. The transformations (2.2) are reduced (2.1) to the following ordinary differential difference equations:

where .

Step 2. We suppose the modified rational Jacobi elliptic functions solutions of (2.4) in the following form:

where , , , are constants to be determined later while satisfies a discrete Jacobi elliptic differential equation:

and , , are arbitrary constants.

Step 3. Since the general solution of the proposed equation (2.6) is difficult to obtain and so the iteration relations corresponding to the general exact solutions we discuss the solutions of the proposed discrete Jacobi elliptic differential equation (2.6) at some special cases to , , and , to cover all the Jacobi elliptic functions as follows.

Type 1. If , , . In this case, (2.6) has the solution , where is the Jacobi elliptic sine function, and is the modulus. In this case from using the properties of Jacobi elliptic functions (see [22]), the series expansion solutions (2.5) take the following form:

Further using the properties of Jacobi elliptic functions, the iterative relations can be written in the following form:

where

, is the jth component of shift vector .

Type 2. If , , and . In this case, (2.6) has the solution . From using the properties of Jacobi elliptic functions, the series expansion solutions (2.5) take the following form:

Type 3. If , , and . In this case, (2.6) has the solution . From using the properties of Jacobi elliptic functions, the series expansion solutions (2.5) take the following form:

From the properties of the Jacobi elliptic functions, we can deduce the iterative relation to the above kind of solutions from Types 2 and 3 as we show in Type 1.

Step 4. Determining the degree of (2.5) by balancing the nonlinear term(s) and the highest order derivatives of in (2.4). It should be noted that the leading terms will not affect the balance because we are interested in balancing the terms of .

Step 5. Substituting in each type form (2.1)–(2.4) and the given values of into (2.4). Cleaning the denominator and collecting all terms with the same degree of , , together, the left-hand side of (2.4) is converted into a polynomial in , , . Setting each coefficient of this polynomial to zero, we derive a set of algebraic equations for , , , , , and .

Step 6. Solving the over determined system of nonlinear algebraic equations by using Maple or Mathematica. We end up with explicit expressions for , , , , , and .

Step 7. Substituting , , , , , and into , in the corresponding type from (2.1)–(2.4), we can finally obtain exact solutions for (2.1).

3. Applications

In this section, we use the proposed method to construct the rational Jacobi elliptic wave solutions for the nonlinear DDEs via the lattice equation, the discrete nonlinear Schrodinger equation with a saturable nonlinearity, the discrete nonlinear Klein-Gordon equation, and the quintic discrete nonlinear Schrodinger equation, which are very important in mathematical physics and have been paid attention to by many researchers.

3.1. Example 1. The Lattice Equation

In this section, we study the lattice equation which takes the following form [23–26]:

where , , and are nonzero constants. This equation contains hybrid lattice equation, mKdV lattice equation, modified Volterra lattice equation, and Langmuir chain equation for some special values of , , . According to the above steps, to seek traveling wave solutions of (3.1), we construct the transformation

where , , and are constants. The transformation (3.2) permits us converting equation (3.1) into the following form:

where . Considering the homogeneous balance between the highest order derivative and the nonlinear term in (3.3), we get . Thus, the solution of (3.3) has the following form:

where , , and are constants to be determined later and satisfies a discrete Jacobi elliptic ordinary differential equation (2.6). When we discuss the solutions of (2.6), we get the following types.

Type 1. If , , and . In this case, the series expansion solution of (3.3) has the form:

With the help of Maple, we substitute (3.5) and (2.8) into (3.3). Cleaning the denominator and collecting all terms with the same degree of , , together, the left-hand side of (3.3) is converted into polynomial in , , . Setting each coefficient of this polynomial to be zero, we derive a set of algebraic equations for , , , , , , , and . Solving the set of algebraic equations by using Maple or Mathematica software package, we have the following.

Family 1.

where , , , , and are arbitrary constants.

From (3.5) and (3.6), the solution of (3.3) takes the following form:

where . Figure 1 illustrates the behavior of the exact solution (3.7).

(a)

(b)

Figure 1: (a) represents the Jacobi elliptic solution (3.7) when , , , , . (b) represents the Jacobi elliptic solution (3.7) when , , , , .

Family 2.

where , , , , and are arbitrary constants.

From (3.5) and (3.8), the solution of (3.3) has the following form:

where .

Type 2. If , , and . In this case, the solution of (3.3) has the form:

With the help of Maple, we substitute (3.10) into (3.3). Cleaning the denominator and collecting all terms with the same degree of , , together, the left-hand side of (3.3) is converted into polynomial in , , . Setting each coefficient of this polynomial to zero, we derive a set of algebraic equations in , , , , , , , and . Solving the set of algebraic equations by using Maple or Mathematica software package, we get the following.

Family 1.

where , , , , and are arbitrary constants.

In this case, the solution of (3.3) takes the following form:

where .

The Jacobi elliptic functions could be generated into the hyperbolic functions when tends to one in the other hand, they are generated into trigonometrical functions when tends to zero.

When , the trigonometrical solution (3.12) takes the following form:

Also if , the hyperbolic solution (3.12) takes the following form:

Family 2.

where , ,, , and are arbitrary constants.

In this case, the solution of (3.3) takes the following form:

where .

When , the trigonometrical solution (3.16) takes the following form:

When , the hyperbolic solution (3.16) takes the following form:

Type 3. If , , and . In this case, the series expansion solution of (3.3) has the form:

Consequently, using the Maple or Mathematica we get the following results.

Family 1.

where , ,, , and are arbitrary constants.

In this case, the solution of (3.3) takes the following form:

where .

Family 2.

where , ,, , and are arbitrary constants.

In this case, the solution of (3.3) takes the following form:

where .

3.2. Example 2. The Discrete Nonlinear Schrodinger Equation

The discrete nonlinear Schrodinger equation (DNSE) is one of the most fundamental nonlinear lattice models [8]. Its arise in nonlinear optics as a model of infinite wave guide arrays [27] and has been recently implemented to describe Bose-Einstein condensates in optical lattices. The class of DNSE model with saturable nonlinearity is also of particular interest in their own right, due to a feature first unveiled in [28]. In this section, we study the DNSE with a saturable nonlinearity [29, 30] form:

which describes optical pulse propagations in various doped fibers, is a complex valued wave function at sites while and are real parameters. We make the transformation:

where , , , and are arbitrary constants. The transformation (3.25) permits us converting equation (3.24) into the following nonlinear difference equation:

We assume that (3.26) has a solution of the form:

where , , and are constants to be determined later and satisfies a discrete Jacobi elliptic differential equation (2.6). When we discuss the solutions of (3.26), we have the following types.

Type 1. If , , and . In this case, the series expansion solution of (3.26) has the form:

With the help of Maple, we substitute (3.28) and (2.8) into (3.26), cleaning the denominator and collecting all terms with the same order of , , together, the left-hand side of (3.26) is converted into polynomial in , , . Setting each coefficient of this polynomial to zero, we derive a set of algebraic equations for , , , , , , and . Solving the set of algebraic equations by using Maple or Mathematica software package, we obtain the following.

Family 1.

In this case, the solution of (3.24) takes the following form:

where denotes and . Figure 2 illustrates the behavior exact solution (3.30).

(a)

(b)

(c)

(d)

Figure 2: (a) and (b) represent the Jacobi elliptic solution (3.30) when , , , , . (c) and (d) represent the Jacobi elliptic solution (3.30) when , , , , .

Family 2.

In this case, the solution of (3.24) takes the following form:

where denotes and .

Type 2. If , , and . In this case, the solution of (3.26) has the form:

Consequently, using Maple or Mathematica we get the following results.

Family 1.

In this case, the solution of (3.24) takes the following form:

where .

In the special case, when , the trigonometrical solution (3.35) takes the following form:

Also if , the hyperbolic solution (3.35) takes the following form:

Family 2.

In this case, the solution of (3.24) takes the following form:

where denotes and .

If , the trigonometrical solution (3.39) takes the following form:

Also if , the hyperbolic solution (3.39) takes the following form:

Type 3. if , , and . In this case, the series expansion solution of (3.26) has the form:

Consequently, using Maple or Mathematica we get the following results.

Family 1.

In this case, the solution of (3.24) takes the following form:

where .

Family 2.

In this case, the solution of (3.24) takes the following form:

where .

3.3. Example 3. The Discrete Nonlinear Klein-Gordon Equation

In this section, we consider the following discrete nonlinear Klein-Gordon equation [31]:

The nonconstant (in contrast to the standard models of harmonic coupling and linear dispersion [32]) function ensures the presence of nonlinear dispersion, which is critical for the existence of compactly supported solutions and can take values in the interval []. Kevrekidis and Konotop [31] have obtained some exact compaction solutions and claim that this DDE does not have the traveling compact solution.

If we set where and are arbitrary constants and take the traveling transformation:

where , , and are constants. The transformation (3.48) permits us converting equation (3.47) into the following form:

where . Considering the homogeneous balance between the highest-order derivative and the nonlinear term in (3.49), we get . Thus, the solutions of (3.49) have the following form:

where , , and are constants to be determined later and satisfies a discrete Jacobi elliptic ordinary differential equation (2.6), we have the following types.

Type 1. If , , and . In this case, the series expansion solution of (3.49) has the form:

With the help of Maple, we substitute (3.51) and (2.8) into (3.49). Cleaning the denominator and collecting all terms with the same degree of , , together, the left-hand side of (3.49) is converted into polynomial in , , . Setting each coefficient of this polynomial to zero, we derive a set of algebraic equations for , , , , , , , and . Solving the set of algebraic equations by using Maple or Mathematica software package, we have the following.

Family 1.

where , , and are arbitrary constants and

In this case, the solution of (3.47) takes the following form:

where .

Figure 3 illustrates the behavior of the exact solution (3.54).

(a)

(b)

Figure 3: (a) represents the Jacobi elliptic solution (3.54) when , , , , . (b) represents the Jacobi elliptic solution (3.54) when , , , , .

Family 2.

where , , and are arbitrary constants and

In this case, the solution of (3.47) takes the following form:

where .

Type 2. If , , and . In this case, the series expansion solution of (3.49) has the form:

Consequently, using Maple or Mathematica we get the following results.

Family 1.

where , , and are arbitrary constants.

In this case, the solution of (3.47) takes the following form:

where .

In the special case, when , the trigonometrical solution (3.60) takes the following form:

Also when , the hyperbolic solution (3.60) takes the following form:

Family 2.

where , , and are arbitrary constants, while

In this case, the solution of (3.47) takes the following form:

where .

When , the trigonometrical solution (3.65) takes the following form:

Also if , the hyperbolic solution (3.65) takes the following form:

Type 3. If , , and . In this case, the series expansion solution of (3.49) has the form:

Consequently, using Maple or Mathematica we get the following results.

Family 1.

where , , and are arbitrary constants and

In this case, the solution of (3.47) takes the following form:

where .

Family 2.

where , , and are arbitrary constants and

In this case, the solution of (3.47) takes the following form:

where .

3.4. Example 4. The Quintic Discrete Nonlinear Schrodinger Equation

We discuss the quintic discrete nonlinear Schrodinger (QDNLS) equation [33]:

which describes the propagation of discrete self-trapped beams in an array of weakly coupled nonlinear optical waveguides. Equation (3.75) was presented for the first time in [33], together with its localized solutions. We are looking for solutions of the form:

where is a real function with respect to the discrete variable . Substitution of equation (3.76) into (3.75) yields the corresponding symmetric equation as follows:

We take the traveling wave of transformation

where and are constants. The transformation (3.78) permits us converting equation (3.77) into the following form:

We suppose that the solutions of (3.79) have the following form:

where , , and are constant to be determined later and satisfies a discrete Jacobi elliptic ordinary differential equation (2.6), we have the following types.

Type 1. If , , and . In this case, the series expansion solution of (3.79) has the form:

With the help of Maple, we substitute (3.81) and (2.8) into (3.79). Cleaning the denominator and collecting all terms with the same degree of , , together, the left-hand side of (3.79) is converted into polynomial in , , . Setting each coefficient of this polynomial to be zero, we derive a set of algebraic equations for , , , , , , , , and . Solving the set of algebraic equations by using Maple or Mathematica software package, we have the following.

Family 1.

where , , and are arbitrary constants and

In this case, the solution of (3.75) takes the following form:

where .

Figure 4 illustrates the behavior of the exact solution (3.84).

(a)

(b)

(c)

(d)

Figure 4: (a) and (b) represent the Jacobi elliptic solution (3.84) when , , , , , , . (c) and (d) represent the Jacobi elliptic solution (3.84) when , , , , , , .

Family 2.

where , , and are arbitrary constants and

In this case, the solution of (3.75) takes the following form:

where .

Type 2. If , , and . In this case, the solution of (3.79) has the form:

Consequently, using the Maple or Mathematica we get the following results.

Family 1.

where , , and are arbitrary constants.

In this case, the solution of (3.75) takes the following form:

where .

In the special case, when , the trigonometrical solution (3.90) takes the following form:

When , the hyperbolic solution (3.90) takes the following form:

where

and .

Family 2.

where , , and are arbitrary constants and

In this case, the solution of (3.75) takes the following form:

where .

When , the trigonometrical solution (3.96) takes the following form:

where .

When , the hyperbolic solution (3.96) takes the following form:

where

and .

Type 3. If , , and . In this case, the series expansion solution of (3.79) has the form:

Consequently, using the Maple or Mathematica we get the following results:

Family 1.

where , , and are arbitrary constants.

In this case, the solution of (3.75) takes the following form:

where .

Family 2.

where , , and are arbitrary constants.

In this case, the solution of (3.75) takes the following form:

where .

Remark 3.1. Wang and Ma [34] constructed the Jacobi elliptic function solution to some nonlinear DDEs in terms of Jacobi elliptic functions or or . In Zhang [20], an algorithm is devised to derive exact traveling wave solutions of nonlinear DDEs by means of Jacobi elliptic functions. Gepreel and Shehata [22] put a direct method to construct the rational Jacobi elliptic solution for nonlinear DDEs. In this paper, we modified the direct method which was disscussed in [22]. When we are constructing the solutions of nonlinear DDEs we neglected the case because if , we get the same solution which was discussed by [22]. The Jacobi elliptic functions could be generated into the hyperbolic functions when tends to one in the other hand, they are generated into trigonometrical functions when tends to zero as shown in Table 1.

Table 1: The relation between the Jacobi elliptic functions and the trigonometrical functions, hyperbolic functions.

0

1

1
tanh 
sech 
sech 
coth 
csch 
csch 

4. Conclusion

In this paper, we modified the rational Jacobi elliptic functions method to calculate some new exact solutions for the nonlinear difference differential equations via the lattice equation, the discrete nonlinear Schrodinger equation with a saturable nonlinearity, the discrete nonlinear Klein-Gordon equation, and the quintic discrete nonlinear Schrodinger equation. As a result, many new and more rational Jacobi elliptic solutions are obtained, from which hyperbolic function solutions and trigonometric function solutions are derived when the moduli and .

References

E. Fermi, J. Pasta, and S. Ulam, Collected Papers of Enrico Fermi II, University of Chicago Press, Chicago, Ill, USA, 1965.

W. P. Su, J. R. Schrieffer, and A. J. Heeger, “Solitons in polyacetylene,” Physical Review Letters, vol. 42, no. 25, pp. 1698–1701, 1979.

A. S. Davydov, “The theory of contraction of proteins under their excitation,” Journal of Theoretical Biology, vol. 38, no. 3, pp. 559–569, 1973.

P. Marquié, J. M. Bilbault, and M. Remoissenet, “Observation of nonlinear localized modes in an electrical lattice,” Physical Review E, vol. 51, no. 6, pp. 6127–6133, 1995.

M. Toda, Theory of Nonlinear Lattices, Springer, Berlin, Germany, 1989.

M. Wadati, “Transformation theories for nonlinear discrete systems,” Progress of Theoretical Physics, vol. 59, pp. 36–63, 1976.

Y. Ohta and R. Hirota, “A discrete KdV equation and its Casorati determinant solution,” Journal of the Physical Society of Japan, vol. 60, no. 6, p. 2095, 1991.

M. J. Ablowitz and J. Ladik, “Nonlinear differential-difference equations,” Journal of Mathematical Physics, vol. 16, pp. 598–603, 1975.

X. B. Hu and W. X. Ma, “Application of Hirota's bilinear formalism to the Toeplitz lattice—some special soliton-like solutions,” Physics Letters A, vol. 293, no. 3-4, pp. 161–165, 2002.

D. Baldwin, U. Goktas, and W. Hereman, “Symbolic computation of hyperbolic tangent solutions for nonlinear differential-difference equations,” Computer Physics Communications, vol. 162, no. 3, pp. 203–217, 2004.

S. K. Liu, Z. T. Fu, Z. G. Wang, and S. D. Liu, “Periodic solutions for a class of nonlinear differential-difference equations,” Communications in Theoretical Physics, vol. 49, no. 5, pp. 1155–1158, 2008.

C. Qiong and L. Bin, “Applications of Jacobi elliptic function expansion method for nonlinear differential-difference equations,” Communications in Theoretical Physics, vol. 43, no. 3, pp. 385–388, 2005.

F. Xie, M. Jia, and H. Zhao, “Some solutions of discrete sine-Gordon equation,” Chaos, Solitons & Fractals, vol. 33, no. 5, pp. 1791–1795, 2007.

S. D. Zhu, “Exp-function method for the Hybrid-Lattice system,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 3, pp. 461–464, 2007.

I. Aslan, “A discrete generalization of the extended simplest equation method,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 8, pp. 1967–1973, 2010.

P. Yang, Y. Chen, and Z. B. Li, “ADM-Padé technique for the nonlinear lattice equations,” Applied Mathematics and Computation, vol. 210, no. 2, pp. 362–375, 2009.

S. D. Zhu, Y. M. Chu, and S. L. Qiu, “The homotopy perturbation method for discontinued problems arising in nanotechnology,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2398–2401, 2009.

S. Zhang, L. Dong, J. Ba, and Y. Sun, “The -expansion method for nonlinear differential-difference equations,” Physics Letters A, vol. 373, no. 10, pp. 905–910, 2009.

I. Aslan, “The Ablowitz-Ladik lattice system by means of the extended -expansion method,” Applied Mathematics and Computation, vol. 216, no. 9, pp. 2778–2782, 2010.

S. Zhang, “Discrete Jacobi elliptic function expansion method for nonlinear difference equation,” Physica Scripta, vol. 80, pp. 045002–045010, 2009.

A. Khaled Gepreel, “Rational Jacobi elliptic solutions for nonlinear difference differential equations,” Nonlinear Science Letters A, vol. 2, pp. 151–158, 2011.

K. A. Gepreel and A. R. Shehata, “Rational Jacobi elliptic solutions for nonlinear differential difference lattice equations,” Applied Mathematics Letters, vol. 25, no. 9, pp. 1173–1178, 2012.

G. Wu and T. Xia, “A new method for constructing soliton solutions to differential-difference equation with symbolic computation,” Chaos, Solitons & Fractals, vol. 39, no. 5, pp. 2245–2248, 2009.

F. Xie and J. Wang, “A new method for solving nonlinear differential-difference equation,” Chaos, Solitons & Fractals, vol. 27, no. 4, pp. 1067–1071, 2006.

C. Liu, “Exponential function rational expansion method for nonlinear differential-difference equations,” Chaos, Solitons & Fractals, vol. 40, no. 2, pp. 708–716, 2009.

Q. Wang and Y. Yu, “New rational formal solutions for -dimensional Toda equation and another Toda equation,” Chaos, Solitons & Fractals, vol. 29, no. 4, pp. 904–915, 2006.

A. Trombettoni and A. Smerzi, “Discrete solitons and breathers with dilute Bose-Einstein condensates,” Physical Review Letters, vol. 86, no. 11, pp. 2353–2356, 2001.

L. Hadžievski, A. Maluckov, M. Stepić, and D. Kip, “Power controlled soliton stability and steering in lattices with saturable nonlinearity,” Physical Review Letters, vol. 93, no. 3, Article ID 033901, 4 pages, 2004.

S. Gate and J. Herrmann, “Soliton propagation in materials with saturable nonlinearity,” Journal of the Optical Society of America B, vol. 9, pp. 2296–2302, 1991.

S. Gate and J. Herrmann, “Soliton propagation and soliton collision in double doped fibers with a non- Kerr- like nonlinear refractive- index change,” Optics Letters, vol. 17, pp. 484–486, 1992.

P. Kevrekidis and V. Konotop, “Bright compact breathers,” Physical Review E, vol. 65, no. 6, pp. 066614–066619, 2002.

P. Kevrekidis, V. Konotop, and A. Bishop, “Discrete compactons: some exact results,” Journal of Physics A, vol. 35, no. 45, pp. L641–L652, 2002.

K. Maruno, Y. Ohta, and N. Joshi, “Exact localized solutions of quintic discrete nonlinear Schrödinger equation,” Physics Letters A, vol. 311, no. 2-3, pp. 214–220, 2003.

Z. Wang and W.-X. Ma, “Discrete Jacobi sub-equation method for nonlinear differential-difference equations,” Mathematical Methods in the Applied Sciences, vol. 33, no. 12, pp. 1463–1472, 2010.