On Backlund transformation of Riccati equation method and its application to nonlinear partial differential equations and differential-difference equations

Author(s): Reham Hassan1, Mustafa El-Agamy1, Mohamed Soror Abdel Latif1, Hamed Nour1
1Mathematics and Engineering Physics Department, Faculty of Engineering, Mansoura University, Mansoura, 35516, Egypt.
Copyright © Reham Hassan, Mustafa El-Agamy, Mohamed Soror Abdel Latif, Hamed Nour. 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

In this paper, we investigate the equivalence between the Backlund transformation of Riccati equation method and the extended tanh-function method. It is proved that the two methods are equivalent when applying them to partial differential equations and differential-difference equations. Two examples are introduced to justify our results.

Keywords: Extended tanh-function method, Backlund transformation of Riccati equation method, partial differential equations, differential-difference equations.

1. Introduction

Many physical, biological and chemical phenomena can be modeled using partial differential equations (PDEs) and differential-difference equations (DDEs). So, in the last decades, many researchers have been interested in obtaining exact solutions of PDEs and DDEs. Many methods were proposed for achieving this task. Some of these methods are: the tanh method [1], the extended tanh-function method (ETM) [2], the simplest equation method [3], the integral bifurcation method [4]Bin}, the extended mapping transformation method [5,6] and the Backlund transformation of Riccati equation method (BTREM) [7,8,9,10,11,12].Our objective in this paper is to investigate the equivalence between the BTREM and the ETM.

2.Description of the two methods

In the following two subsections we give a brief description of the two methods.

2.1.The extended tanh-function method [2]

Consider a given partial differential equation with some independent variables, say, x and t and dependent variable u:
H(u,ut,ux,uxx,uxt,)=0,
(1)
where u=u(x,t) is an unknown function, ux and ut are the derivatives of u with respect to x and t respectively. It is assumed that the Equation(1) has the following traveling wave solution:
u=u(z),z=kx+ct+z0,
(2)
where k,c and z0 are some constants. Substituting Equation(2) into Equation (1), we get the following reduced ordinary differential equation:
H(u,u,u,)=0,
(3)
where the primes denote the derivative with respect to z. The solution of the Equation (3) can be expressed as:
u=i=0naiϕi(z),
(4)
where ai,i=0,1,2,,n are some constants that will be computed later, n is a positive integer computed by the balance between the highest-order derivative term and the nonlinear terms in the Equation(3) and ϕ satisfies the following Riccati equation:
ϕ(z)=σ+ϕ2(z),
(5)
where σ is a constant. The Riccati Equation (5) has the following solutions:
  1. If σ<0
    ϕ(z)=σtanh(σz),
    (6)
    ϕ(z)=σcoth(σz),
    (7)
  2. If σ=0
    ϕ(z)=1z+ω,ω=const.
    (8)
  3. If σ>0
    ϕ(z)=σtan(σz),
    (9)
    ϕ(z)=σcot(σz).
    (10)
Substituting Equation (4)into Equation (3) and making use of Equation (5), then setting the coefficients of ϕi(z),i=0,1, to zero, we get a set of algebraic equations for ai,i=0,1,2,,n . Solving this obtained system will lead to the values of ai,i=0,1,2,,n.

2.2.Backlund transformation of Riccati equation method [10]

In this method the solution of the Equation (3) is given in the form:
u=i=0nbiΦi(z),
(11)
where Φ(z) is given by:
Φ(z)=σB+Dϕ(z)D+Bϕ(z),
(12)
bi,i=0,1,2,,n are some constants that will be computed later, n is a positive integer computed by the balance between the highest-order derivative term and nonlinear terms in the Equation (3), ϕ(z) are the known solutions of Ricatti Equation (5), B and D are arbitrary constants. Substituting Equation (11) into Equation (3), then setting the coefficients of ϕ(z) to zero, we get some algebraic equations for bi,i=0,1,2,,n. Solving this system of algebraic equations will lead to the values of bi,i=0,1,2,,n.

3.Equivalence of the two methods

Case 1 :when ϕ(z)=σtanh(σz). In this case, we have

Φ(z)=σBDσtanh(σz)DBσtanh(σz)=σBσDtanh(σz)1BσDtanh(σz).
(13)
By assuming that (BσD)=tanh(k1), k1 is a constant, we get k1=tanh1(BσD). Therefore,
Φ(z)=σtanh(σz)+tanh(k1)1+tanh(k1)tanh(σz)=σtanh(σz+k1).
(14)
It is clear that Φ(z) and ϕ(z) are only differed by the constant phase shift k1.

Case 2: when ϕ(z)=σcoth(σz). In this case, we get

Φ(z)=σBDσcoth(σz)DBσcoth(σz)=σ1DBσcoth(σz)DBσcoth(σz).
(15)
By setting (DBσ)=coth(k2), k2 is a constant, we obtain k2=coth1(DBσ). Therefore,
Φ(z)=σ1+coth(k2)coth(σz)coth(σz)+coth(k2)=σcoth(σz+k2).
(16)
It is clear that Φ(z) and ϕ(z) are only differed by the constant phase shift k2.

Case 3: when ϕ(z)=1z+ω, σ=0. In this case, we have

Φ(z)=D(1z+ω)D+B(1z+ω)=DDz+DωB=1z+ωBD.
(17)
It is clear that Φ(z) and ϕ(z) are only differed by the constant BD.

Case 4: when ϕ(z)=σtan(σz). In this case, we get

Φ(z)=σB+Dσtan(σz)D+Bσtan(σz)=σBσD+tan(σz)1+BσDtan(σz).
(18)
Assuming that (BσD)=tan(k3), k3 is a constant, we get k3=tan1(BσD). Therefore,
Φ(z)=σtan(σz)+tan(k3)1tan(k3)tan(σz)=σtan(σz+k3).
(19)
It is clear that Φ(z) and ϕ(z) are only differed by the constant phase shift k3.

Case 5: when ϕ(z)=σcot(σz). In this case, we get

Φ(z)=σBDσcot(σz)DBσcot(σz)=σ1+DBσcot(σz)DBσcot(σz).
(20)
By setting (DBσ)=cot(k4), k4 is a constant, we get k4=cot1(DBσ). Therefore,
Φ(z)=σcot(k4)cot(σz)1cot(σz)+cot(k4)=σcot(σz+k4).
(21)
It is clear that Φ(z) and ϕ(z) are only differed by the constant phase shift k4

.

4.The Drinfeld-Sokolov-Wilson equation

The Drinfeld-Sokolov-Wilson equation is given by [10]
ut+Pvvx=0,vt+ruvx+suxv+qvxxx=0,
(22)
where p,q,r and s are some nonzero constants. The authors in [10] have introduced the traveling wave transformation:
u(x,t)=U(z),v(x,t)=V(z),z=k(xct),
(23)
where k and c are constants. Substituting Equation (23) into Equation (22), we obtain the following ordinary differential equations:
kcU+pkVV=0,kcV+rkUV+skUV+qk3V=0.
(24)

After applying the BTREM in [10], the authors have obtained four solutions for the Equation (22). These four solutions are the same solutions obtained in [13] as will be shown in the following discussion.

The first solution is given by
(25)u1(x,t)=6cr+2s(c2qk2BDtanh(c2qk2(k(xct)))Dc2qk2Btanh(c2qk2(k(xct))))2=6cr+2s(c2qk2BD+tanh(c2qk2(k(xct)))1c2qk2BDtanh(c2qk2(k(xct))))2.
(25)
Assume that (c2qk2BD)=tanh(k5), k5 is a constant. Therefore,
(26)u1(x,t)=6cr+2s(tanh(k5)+tanh(c2qk2(k(xct)))1+tanh(k5)tanh(c2qk2(k(xct))))2=6cr+2s(tanh(c2qk2(k(xct)+k5)))2,
(26)
v1(x,t)=±12c2p(r+2s)(c2qk2BDtanh(c2qk2(k(xct)))Dc2qk2Btanh(c2qk2(k(xct))))=±12c2p(r+2s)(tanh(c2qk2(k(xct)+k5))),
(27)
which is the same solution given in [13]. They are only differed by the phase shift constant k5. The second solution is given by:
u2(x,t)=6cr+2s(c2qk2BDcoth(c2qk2(k(xct)))Dc2qk2Bcoth(c2qk2(k(xct))))2=6cr+2s(12qk2cDBcoth(c2qk2(k(xct)))2qk2cDB+coth(c2qk2(k(xct))))2.
(28)
Let (2qk2cDB)=coth(k6), k6 is a constant. Therefore,
u2(x,t)=6cr+2s(1+coth(k6)coth(c2qk2(k(xct)))coth(k6)+coth(c2qk2(k(xct))))2=6cr+2s(coth(c2qk2(k(xct)+k6)))2,
(29)
v2(x,t)=±12c2p(r+2s)(c2qk2BDcoth(c2qk2(k(xct)))Dc2qk2Bcoth(c2qk2(k(xct))))=±12c2p(r+2s)(coth(c2qk2(k(xct)+k6))),
(30)
which is the same solution given in [13]. They are only differed by the phase shift constant k6. The third solution is given by:
u3(x,t)=6cr+2s(c2qk2B+Dtan(c2qk2(k(xct)))D+c2qk2Btan(c2qk2(k(xct))))2=6cr+2s(c2qk2BD+tan(c2qk2(k(xct)))1+c2qk2BDtan(c2qk2(k(xct))))2.
(31)
Assume that (c2qk2BD)=tan(k7), k7 is a constant. Therefore,
u3(x,t)=6cr+2s(tan(k7)+tan(c2qk2(k(xct)))1tan(k7)tanh(c2qk2(k(xct))))2=6cr+2s(tan(c2qk2(k(xct)+k7)))2,
(32)
v3(x,t)=±12c2p(r+2s)(c2qk2B+Dtan(c2qk2(k(xct)))D+c2qk2Btan(c2qk2(k(xct))))=±12c2p(r+2s)(tan(c2qk2(k(xct)+k7))),
(33)
which is the same solution given in [13]. They are only differed by the phase shift constant k7. The forth solution is given by:
u4(x,t)=6cr+2s(c2qk2B+Dcot(c2qk2(k(xct)))D+c2qk2Bcot(c2qk2(k(xct))))2=6cr+2s(12qk2cDBcot(c2qk2(k(xct)))2qk2cDB+cot(c2qk2(k(xct))))2.
(34)
Let (2qk2cDB)=cot(k8), k8 is a constant. Therefore, u4(x,t)=6cr+2s(1+cot(k8)cot(c2qk2(k(xct)))cot(k8)+cot(c2qk2(k(xct))))2
=6cr+2s(cot(c2qk2(k(xct)+k8)))2,
(35)
v4(x,t)=±12c2p(r+2s)(c2qk2B+Dcot(c2qk2(k(xct)))D+c2qk2Bcot(c2qk2(k(xct))))=±12c2p(r+2s)(cot(c2qk2(k(xct)+k8))),
(36)
which is the same solution given in [13]. They are only differed by the phase shift constant k8.

5.The equivalence between the two methods when solving differential-difference equations

In this section, we also prove that the BTREM is equivalent to the ETM when applied to differential-difference equations. To achieve this task we choose the following example.

5.1.The discrete mKdV equation

The discrete mKdV equation is given by [14]:
un(t)t=(θun2)(un+1un1),
(37)
where θ is a constant. To get the traveling wave solutions for Equation (37), the following transformation was introduced [14]:
un(t)=u(ξn),ξn=dn+c1t+c0,
(38)
to transform Equation (37) into:
c1u(ξn)=(θun2(ξn))(un+1(ξn)un1(ξn))
(39)
where d, c1 and c0 are constants. After using the BTREM, the authors in [14] have obtained the following solutions: The first solution is given by:
u1(ξn)=a0+a1rb+artan(rξn)a+brtan(rξn)=a0+a1rrbar+tan(rξn)1+bratan(rξn),
(40)
where a, b, r, a0 and a1 are constants. Let (rbar)=tan(m1), m1 is a constant. Therefore,
u1(ξn)=a0+a1rtan(m1)+tan(rξn)1tan(m1)tan(rξn)=a0+a1rtan(rξn+m1),
(41)
which is a solution in the form of the tan function only with a constant phase shift m1. The second solution is given by:
u2(ξn)=a0+a1rbarcot(rξn)abrcot(rξn)=a0a1rabrcot(rξn)1cot(rξn)abr.
(42)
Let (abr)=cot(m2), m2 is a constant. Therefore,
u2(ξn)=a0a1rcot(m2)cot(rξn)1cot(rξn)+cot(m2)=a0a1rcot(rξn+m2),
(43)
which is a solution in the form of the cot function only with a constant phase shift m2. The third solution is given by:
u3(ξn)=a0+a1rbarcoth(rξn)abrcoth(rξn)=a0a1rabrcoth(rξn)+1coth(rξn)abr.
(44)
Let (abr)=coth(m3), m3 is a constant. Therefore,
u3(ξn)=a0a1rcoth(m3)coth(rξn)+1coth(rξn)+coth(m3)=a0a1rcoth(rξn+m3),
(45)
which is a solution in the form of the coth function only with a constant phase shift m3.

6.Conclusion

We have proved that the the BTREM is equivalent to the ETM. We demonstrated this fact using two examples from partial differential equations and differential-difference equations.

Author Contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Conflict of Interests

The authors declare no conflict of interest.

References

  1. Bekir, A. (2008). Applications of the extended tanh method for coupled nonlinear evolution equations. Communications in Nonlinear Science and Numerical Simulation, 13, 1748-1757. [Google Scholor]
  2. Fan, E. (2000). Extended tanh-function method and its applications to nonlinear equations. Physics Letters A, 277 , 212-218. [Google Scholor]
  3. Ebaid, A., & Abd Elazem, N. Y. (2011). On the exact solutions of a nano boundary layer problem using the simplest equation method. Physica Scripta, 84 , 065005. [Google Scholor]
  4. Bin, H., Weiguo, R., Can, C., & Shaolin, L. (2008). Exact travelling wave solutions of a generalized Camassa-Holm equation using the integral bifurcation method. Applied Mathematics and Computation, 206 , 141-149.[Google Scholor]
  5. Zhao, H., & Cheng, L. B. (2005). Extended mapping transformation method and its applications to nonlinear partial differential equation(s). Communications in Theoretical Physics, 44 , 473-478. [Google Scholor]
  6. Abdel Latif, M. S. (2011). Some exact solutions of KdV equation with variable coefficients. Communications in Nonlinear Science and Numerical Simulation, 16 , 1783-1786. [Google Scholor]
  7. Hon, Y. C., Zhang, Y., & Mei, J. (2010). Exact solutions for differential-difference equations by backlund transformation of riccati equation. Modern Physics Letters B, 24 ,(27), 2713-2724. [Google Scholor]
  8. Arnous, A. H., Mirzazadeh, M., & Eslami, M. (2014). The Backlund transformation method of Riccati equation applied to coupled Higgs field and Hamiltonian amplitude equations. Computational Methods for Differential Equations, 2 (4), 216-226.[Google Scholor]
  9. Arnous, A. H., Mirzazadeh, M., Moshokoa, S., Medhekar, S., Zhou, Q., Mahmood, M. F., Biswas, A. & Belic, M. (2015). Solitons in optical metamaterials with trial solution approach and Backlund transform of Riccati equation. Journal of Computational and Theoretical Nanoscience, 12 (12), 5940-5948. [Google Scholor]
  10. Arnous, A. H., Mirzazadeh, M., & Eslami, M. (2016). Exact solutions of the Drinfel’d-Sokolov-Wilson equation using Backlund transformation of Riccati equation and trial function approach. Pramana, 86 , 1153-1160. [Google Scholor]
  11. El-Borai, M. M., El-Owaidy, H. M., Ahmed, H. M., Arnous, A. H., & Mirzazadeh, M. (2017). Solitons and other solutions to the coupled nonlinear Schrodinger type equations. Nonlinear Engineering, 6 (2), 115-121. [Google Scholor]
  12. Zayed, E. M., Alurrfi, K. A., & Al Nowehy, A. G. (2017). Many exact solutions of the nonlinear kpp equation using the backlund transformation of the Riccati equation. International Journal of Optics and Photonic Engineering, 2 (1). [Google Scholor]
  13. Bibi, S., & Mohyud-Din, S. T. (2014). New traveling wave solutions of Drinefel’d- Sokolov- Wilson Equation using Tanh and Extended Tanh methods. Journal of the Egyptian Mathematical Society, 22 , 517-523. [Google Scholor]
  14. Zhang, Y., Hon, Y. C., & Mei, J. (2010). A systematic method for solving differential-difference equations. Communications in Nonlinear Science and Numerical Simulation, 15 , 2791-2797. [Google Scholor]