Extended Jacobian Elliptic Function Expansion Method and Its Applications in Biology ()
1. Introduction
The nonlinear partial differential equations of mathematical physics are major subjects in physical science [1] . Exact solutions for these equations play an important role in many phenomena in physics such as fluid mechanics, hydrodynamics, Optics, Plasma physics and so on. Recently many new approaches for finding these solutions have been proposed, for example, tanh-sech method [2] -[4] , extended tanh-method [5] -[7] , [8] - [11] , homogeneous balance method [12] , F-expansion method [13] -[15] , exp-function method [16] [17] , trigono-
metric function series method [18] , -expansion method [19] -[22] , Jacobi elliptic function method [23] -[26]
and so on.
The objective of this article is to apply the extended Jacobian elliptic function expansion method for finding the exact traveling wave solution of Dynamical system in a new Double-Chain Model of DNA and a diffusive predator-prey system which play an important role in biology and mathematical physics.
The rest of this paper is organized as follows: In Section 2, we give the description of the extended Jacobi elliptic function expansion method In Section 3, we use this method to find the exact solutions of the nonlinear evolution equations pointed out above. In Section 4, conclusions are given.
2. Description of Method
Consider the following nonlinear evolution equation
(1)
where F is polynomial in and its partial derivatives in which the highest order derivatives and nonlinear terms are involved. In the following, we give the main steps of this method [23] -[26]
Step 1. Using the transformation
(2)
where k and c are the wave number and wave speed, to reduce Equation (1) to the following ODE:
(3)
where P is a polynomial in and its total derivatives, while.
Step 2. Making good use of ten Jacobian elliptic functions, we assume that (3) has the solutions in these forms:
(4)
With
(5)
where, , , are the Jacobian elliptic sine function, The jacobian elliptic cosine function and the Jacobian elliptic function of the third kind and other Jacobian functions which is denoted by Glaisher’s symbols and are generated by these three kinds of functions, namely
(6)
That have the relations
(7)
With the modulus m. In addition we know that
(8)
The derivatives of other Jacobian elliptic functions are obtained by using Equation (8). To balance the highest order linear term with nonlinear term we define the degree of u as which gives rise to the degrees of other expressions as
(9)
According the rules, we can balance the highest order linear term and nonlinear term in Equation (3) so that n in Equation (4) can be determined.
In addition we see that when, , , and degenerate as, , , re- spectively, while when therefore Equation (5) degenerate as the following forms
(10)
(11)
(12)
(13)
Therefore the extended Jacobian elliptic function expansion method is more general than sine-cosine method, the tan-function method and Jacobian elliptic function expansion method.
3. Application
3.1. Example 1: Dynamical System in a New Double-Chain Model of DNA
An attractive nonlinear model for the nonlinear science in the deoxyribonucleic acid (DNA). The dynamics of DNA molecules is one of the most fascinating problems of modern biophysics because it is at the basis of life. The DNA structure has been studied during last decades. The investigation of DNA dynamics has successfully predicted the appearance of important nonlinear structures. It has been shown that the non linearity is respon- sible for forming localized waves. These localized waves are interesting because they have the capability to transport energy without dissipation [27] -[35] . In Ref. [34] [35] , it is given that a new double-chain model of DNA consists of two long elastic homogeneous strands which represent two poly nucleotide chains of the DNA molecule, connected with each other by an elastic membrane representing the hydrogen bonds between the base pair of the two chains. Under some appropriate approximation, the new double-chain model of DNA can be described by the following two general nonlinear dynamical system:
(14)
(15)
where
(16)
where, , Y and F denote respectively the mass density, the area of transverse cross-section, the Young’s modulus and tension density of each strand; is the rigidity of the elastic membrane; h is the distance be- tween the two strands, and is the height of the membrane in the equilibrium positive. In Equations (14) and (15), u is the difference of the longitudinal displacements of the bottom and top strands, while v is the difference of the transverse displacements of the bottom and top strands.
we first introduce the transformation
(17)
where a and b are constants, to reduce Equations (14) and (15) to the following system of equations:
(18)
and
(19)
Comparing Equations (18) and (19) and using (17) we deduce that and. Now Equations (18)
and (19) can be written as
(20)
where
(21)
The wave transformation, , reduce Equation (20) to the following ODE:
(22)
where. Balancing and yields,. Consequently, we have the for- mal solution:
(23)
where, and are constant such that or. From (23), it is easy to see that
(24)
(25)
Substituting Equations (23) and (25) into Equation (22) and equating all coefficients of, , , , , , to zero, we obtain
(26)
(27)
(28)
(29)
(30)
(31)
(32)
Solving the above system with the aid of Maple or Mathematica, we have the following solitary wave solution:
Case 1.
(33)
Case 2.
(34)
Case 3.
(35)
Sothat solution of Equation (22) has the form
Case 1.
(36)
Case 2.
(37)
Case 3.
(38)
3.2. Example 2. A Diffusive Predator-Prey System
Consider a system of two coupled nonlinear partial differential equations describing the spatio-temporal dy- namics of a predator-prey system [36] ,
(39)
where, , m and are positive parameters. The solutions of predator-prey system have been studied in various aspects [36] -[38] . The dynamics of the diffusive predator-prey system have assumed the following
relations between the parameters, namely and. Under there assumptions, Equation (39)
can be rewritten in the form:
(40)
We use the wave transformation to reduce Equation (40) to the following non- linear system of ordinary differential equations:
(41)
where c is a nonzero constant.
In order to solve Equation (41), let us consider the following transformation
(42)
Substituting the transformation (42) into Equation (41), we get
(43)
Balancing with in Equation (43) yields,. Consequently, we get the same for- mal solution (23). Substituting (23)-(25) into (43), setting the coefficients of (, , , , , ,) to zero, we obtain the following under determined system of algebraic equations for (, ,).
(44)
(45)
(46)
(47)
(48)
(49)
(50)
solving Equations (44)-(50) using the maple or mathematica program to get solitary wave solution of equations we get
So we get
(51)
(52)
when hyperbolic solution
(53)
4. Conclusion
We establish exact solutions for the dynamics of DNA molecules which is one of the most fascinating problems of modern biophysics because it is at the basis of life. The DNA structure has been studied during last decades. The investigation of DNA dynamics has successfully predicted the appearance of important nonlinear structures and a system of two coupled nonlinear partial differential equations describing the spatio-temporal dynamics of a predator-prey system where the prey per capita growth rate is subject to the All effect. The extended Jacobian elliptic function expansion method has been successfully used to find the exact traveling wave solutions of some nonlinear evolution equations. As an application, the traveling wave solutions for Dynamical system in a new Double-Chain Model of DNA and a diffusive predator-prey system, which have been constructed using the extended Jacobian elliptic function expansion method. Let us compare between our results obtained in the present article with the well-known results obtained by other authors using different methods as follows: Our results of the system of shallow water wave equations and a diffusive predator-prey system, are new and different from those obtained in [34] -[38] . It can be concluded that this method is reliable and proposes a variety of exact solutions NPDEs. The performance of this method is effective and can be applied to many other nonlinear evolution equations.