Traveling Wave Solution of the Olver–Rosenau Equation Solved by Dynamics System

Olver–Rosenau equations presented by Olver and Rosenau can be rewritten to the dynamic system by the wave transformation. The system is a Hamiltonian system with the ﬁrst integral, and its phase-space and equilibrium point analysis are given in diﬀerent parameter spaces in detail. On this basis, we can derive various solutions of the original equation relating these orbits in diﬀerent phase-space planes, and the theoretical basis of the numerical solution is provided for engineering application and production practice.


Introduction
In recent years, nonlinear partial differential equations (NLPDEs) are more and more extensively used to engineering application and production practice. Because of the difficulty of solving these equations, so many methods have been presented in the last decades.
e Camassa-Holm (CH) equation including the cubic term, that is, the Fokas-Olver-Rosenau-Qiao (FORQ) equation, has been studied for a long time and achieved so many results, such as the solution of Cauchy initial value problem [1], Holder continuous [2], the algebro-geometric solutions [3], the Cauchy problem of the generalized equation [4], and the nonuniqueness for the equation [5]. A lot of solving methods for the nonlinear partial differential equation are discussed to be applied in engineering and practice areas, such as the sine-Gordon expansion method [6] and the travelling wave method and its conservation laws [7], and so many examples are in this regard. In this paper, the Olver-Rosenau equation was discussed by Olver and Rosenau through a reshuffling procedure of the Hamiltonian operators in 1996, and they found that the equation was changed into the Hamiltonian system by the bitransformation structure of the mKdV equation [8]. Rosenau believed the nonanalytic solitary waves of the equation in 1997 [9]. Later, the equation was derived in 2013 [10]. e resulting Hamiltonian equations are considered by the dynamical system theory and a phasespace analysis of their singular points. ose results of the study proved that the equations can support double compacton solutions. ey found that the new Olver-Rosenau compactons are different from the well-known Rosenau-Hyman compacton and Cooper-Shepard-Sodano compacton. It was recently introduced by Li in [11], but Li did not give the solution because of the complexity of the integral. We will give more detailed discussion in this paper. In the second section, we discuss bifurcations and phase portraits of the system in all parameters. In the third section, smooth (or bright) solitary wave solutions of the system and their parametric representations are obtained in detail. Let Integrating once with respect to ζ and letting dϕ/dζ � ω, equation (3) can be rewritten as which has the first integral where k is an integral constant. Without loss of generality, we consider the case + instead of " ± ." Equation (5) is simplified as the dynamic system dϕ dζ � ω, with the Hamiltonian Obviously, system (6) is a planar dynamical system with three parameters depended on the parameter group (b, c, and k). All possible phase portraits of (6) in the (ϕ, ω) phase are discussed under the conditions of the different parameter group (b, c, and k).

Phase Portraits of the System
Imposing the transformation dζ � (ω 2 + ϕ 2 + 2c)dτ, when ω 2 + ϕ 2 + 2c ≠ 0, equation (6) leads to associated regular system dϕ dτ is system has the same first integral as equation (6). Apparently, the singular curve ω 2 + ϕ 2 + 2c � 0 is related to the singular solution of equation (8). Near the circumference, the variable τ is a fast variable, while the variable ζ is a slow variable in the sense of the geometric singular perturbation theory.
In order to find the equilibrium points of system (8), let e zero points of f(ϕ) are estimated by the monotonicity of the function ϕ based on zero points of the derivative function f ′ (ϕ). Let A � − 6(b + c), B � − 9k, and C � (2c + 2b) 2 ; then, the discriminant S � B 2 − 4AC of the 3 . It is easy to see that, for given b and c when 3 , we have S < 0. It follows that there exist three simple real roots ϕ 31 , ϕ 32 , and ϕ 33 of f(ϕ) satisfying ϕ 31 < ϕ 1 < ϕ 32 < ϕ 2 < ϕ 33 .
If k 2 � k 2 1 , there exist two simple real roots ϕ 21 and ϕ 22 On the singular circle line, there exist two equilibrium points By the theory of planar dynamical systems, for an equilibrium point of a planar integrable system, if J < 0, then the equilibrium point is a saddle point; if J > 0 and (traceM) 2 − 4J < 0( > 0), then it is a center point (a node point); and if J � 0 and the Poincare index of the equilibrium point is 0, then this equilibrium point is cusped (see [12]). Let where H comes from equation (5). (2) If A � B, there exists only one simple real root ϕ � 0, and then k � (2/3)(b + c).  Mathematical Problems in Engineering      (4) If S > 0, there exists one simple real root, and then and 3 ).
We obtain four class figures totally including 13 phase portraits under the conditions of different parameters. An orbit in a figure is related to Hamiltonian h of equation (7). An orbit derives a class solution of system (6), and at the same time, it is appropriative to the original partial differential equation. In the next section, we mainly consider how to solve system (6).

Smooth (Bright) Solitary Wave Solutions and Their Parametric Representations
As we know, system (6) for a fixed integral constant h has an explicit solution. However, the scope of value h is calculated by different orbits corresponding to different parameter conditions in the phase portraits from Figures 1 to 13.
To write this equation as integral form, As it is very difficult to solve the left integral directly, we believe that square root F � h + c 2 − bϕ 2 − kϕ plays a critical role in our discussion. If F is a perfect square of the function ϕ, we can obtain analytical solution because of the simplicity of the integral. If F is a complete square, then the algebraic expression in the first root in the above formula becomes rational, and the integration is relatively easy solution.

Analytical Solution. Supporting that
If we want to remove the second root sign of equation (11), we must make F to be complete square, and the necessary and sufficient condition for F to be completely square is also given: is not a rational function, the integral of equation (11) is very difficult to be calculated directly. If F cannot be expressed as a complete square, we can only carry out Taylor expansion of F; otherwise, equation (13) cannot calculate the integral, and the solution of the original equation cannot be obtained. While we need to approximate nonlinear function, our right choice is the Taylor display: where the remainder    Mathematical Problems in Engineering generally satisfying |A(ϕ + (k/2b)) 2 | < 1, Error |d| ≤ (1/8)(1 + Aθ(ϕ + (k/2b))) 1/2 .
that is, Let ; then the above equation of the integral becomes Approximation solution of system (6) is v � ϕ(x − c * t) is the approximation solution of the original equation.
If the order of Taylor expansion of F is higher, the accuracy of the solution will be higher.

Conclusion
We obtain some exact solutions and some approximation solutions of the Olver-Rosenau equation by the dynamic system method. e original equation considered by the first integral method and its phase-space analysis and equilibrium points are calculated under different parameter conditions. We can derive various solutions of the original equation relating these orbits in different phase-space planes. Nevertheless, as there are some troubles in the course of the calculation of these solutions because of the integral complexity, we need to find some numerical methods for the equation. e precision differs from the variable order of Taylor expansion.

Data Availability
is paper is mainly theoretical derivation and calculation, it does not involve data and its right to use.

Conflicts of Interest
e authors declare that they have no conflicts of interest.