The nature of the world is known by people step by step with many powerful methods. Solving linear equation is the beginning of the process, meanwhile, the results do not well agree with the solution of the linear equation. People overcome the difference and find the nonlinear equation, which is used to describe many phenomena in various physics fields. Many different nonlinear equations appear gradually. Nonlinear Vakhnenko equation is a kind of nonlinear partial differential equation, which is proposed to describe long waves of small amplitudes broadcasting in nonlinear dispersive media. The equation and its variable coefficients play an important role in further study mathematics and knowing the physics of nature. However, obtaining the solutions of the nonlinear equation is a difficult story. With deeply studying the relation of nature, many powerful methods appear  - . The nonlinear Vakhnenko equation is investigated in many papers    ; some new solutions in the form of Jacobi elliptic function are given in this paper, which enrich the kinds of the solution for the equation.
Generally, the nonlinear Vakhnenko equation  is read as
where the subscripts denote the partial derivatives of x and t. The remaining structure of this paper is organized as follows: Section 2 is a brief introduction to the Jacobi elliptic function and its properties. In Section 3, by implementing the Jacobi elliptic function expansion method, some new traveling wave solutions for nonlinear Vakhnenko are reported. The conclusion is summarized in the last section.
2. The Jacobi Elliptic Function Expansion Method and Its Properties
The nonlinear partial differential equation with independent variables x and t is generally in the following form
The above equation is a function about , the subscripts denote the partial derivatives with x and t, respectively. The wave variable is applied to the Equation (2), which is changed into the following ordinary differential equation
where denotes the derivative with respect to the same variable .
Generally, the function u in terms of Jacobi elliptic function expansion method can be expressed with the first kind of Jacobi elliptic function as follow
where ai is constant parameter; m is determined by balancing the linear term of the highest order derivative with nonlinear term of Equation (3).
The three kinds of Jacobi elliptic functions have useful properties     , which are given as follow:
where m (0 < m < 1) is the modulus of the elliptic function. The modulus substantially affects the Jacobi elliptic solutions, which will asymptotically go into hyperbolic functions and trigonometric functions when the modulus and , respectively. The asymptotical functions are listed in the following
3. Jacobi Elliptic Function Expansion Method for Nonlinear Vakhnenko Equation
Consider the Nonlinear Vakhnenko equation in the form of Equation (1) with , we obtain
By using m to balance the highest order derivative term and the nonlinear term of Equation (9), we have m = 2 obtained from Equation (9). Then, Equation (4) reduced as
where is the first kind of Jacobi elliptic function.
Substituting Equation (10) into the equation Equation (9) and collecting all terms with the same power term together and equating each coefficient of to zero, a set of simultaneous algebraic equations are yielded as follows:
Solving the above algebraic equations, we get the results:
By using Equation (10), we obtained the solution of Equation (1) as the following:
The numerical simulation for Equations ((17), (18)) with , , and are shown in Figure 1 and Figure 2, respectively. The values of the parameters here and the following are determined randomly in the corresponding interval.
According to the asymptotically functions, the numerical simulation for Equation (18) with , , , and , are shown in Figure 3 and Figure 4, respectively. The Equation (18) is changed into hyperbolic function and trigonometric function when the modulus approaches 1 ( ) and 0 ( ), which are listed as follows
The derivative of the Equation (18) is given as
The numerical simulation with , , and is shown in Figure 5.
Figure 1. The numerical simulation for Equation (17) with , , and .
Figure 2. The numerical simulation for Equation (18) with , , and .
Figure 3. The numerical simulation for Equation (18) with , , and.
Figure 4. The numerical simulation for Equation (18) with, , and.
Figure 5. The numerical simulation for Equation (21) with, , and.
The nonlinear Vakhnenko equation in the form Equation (1) has been further studied by means of the Jacobi elliptic function expansion method. Several kinds of periodic evolution solutions (17)-(21) of Equation (1) are obtained in terms of Jacobi elliptic functions, hyperbolic solutions and trigonometric solutions. Those solutions can be used to discuss some interesting physical phenomena of amplitudes broadcasting in nonlinear dispersive media. The performance of Jacobi elliptic functions expansion method is powerful, which gives more solutions. From the simulation results, we can know the amplitudes of travelling wave solutions and period. Equation (21) is the derivative of the Equation (18), which shows the velocity of the wave solutions.
 Wang, Y. and Wang, B. (2005) High-Order Multi-Symplectic Schemes for the Nonlinear Klein-Gordon Equation. Applied Mathematics and Computation, 166, 608-632.
 Liu, S.K., Fu, Z.T. and Liu, S.D. (2001) A Simple Fast Method in Finding Particular Solutions of Some Nonlinear PDE. Applied Mathematics and Mechanics, 22, 326-331.
 Vu-Quoc, L. and Li, S. (1993) Invariant-Conserving Finite Difference Algorithms for the Nonlinear Klein-Gordon Equation. Computer Methods in Applied Mechanics and Engineering, 107, 341-391.
 Yusufolu, E. and Bekir, A. (2008) The Tanh and the Sine-Cosine Methods for Exact Solutions of the MBBM and the Vakhnenko Equations. Chaos Solitons & Fractals, 38, 1126-1133.
 Liu, G.-T. (2009) Lame Equation and Asymptotic Higher-Order Periodic Solutions to Nonlinear Evolution Equations. Applied Mathematics and Computation, 212, 312-317.