**Journal of Applied Mathematics and Physics**

Vol.05 No.09(2017), Article ID:79359,18 pages

10.4236/jamp.2017.59155

The Qualitative and Quantitative Methods of Kovalevskys Case

Fawzy Mohamed Fahmy El-Sabaa, Alshimaa Abdelbasit Mohamed^{*}, Salma Khalel Zakria^{ }

Department of Mathematics, Faculty of Education, Ain Shams University, Cairo, Egypt

Copyright © 2017 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: July 1, 2017; Accepted: September 24, 2017; Published: September 27, 2017

ABSTRACT

The phase portrait of the functions obtained by Hamilton-Jacobi equations is substantiated, and the classification of singular points is found, and the bifurcation diagram for the problem is studied. The numerical calculation by using Poincaré surface section is used to get the invariant tori for our problem.

**Keywords:**

Hamilton-Jacobi’s Equations, Phase Portrait, Poincaré Surface-Section, Qualitative and Quantitative Methods

1. Introduction

The problem of the motion of a rigid body with the fixed point in the field of constant gravity is one of the oldest in mechanics [1] . The system of the motion contains three degrees of freedom (Euler’s angles). The system possesses three first integrals, and according to Jacobi’s last multiplier, the system can be integrated if there are other first integrals. Kovalevskaya [2] proved that the only one-valued first integral can be found in three cases: first, in Euler’s case, where the fixed point is located at the center of mass and rotation occurs freely without the influence of torque; second, in the Lagrange’s case when the rigid body is a symmetric top with the center of mass on the axis of symmetry; thirdly, in Kovalevskaya’s case in which the ellipsoid of inertia about a fixed is symmetric, and the moments of inertia at the supporting point are equal, and each one of them is twice as much as the third one and the center of the mass in the plane of equal moments of inertia. Kovalevskaya [3] showed that other integral cannot be single valued over the whole t-complex plane and Hussen [4] found that, when ellipsoid of inertia is not symmetric, a new algebraic integral cannot exist, for arbitrary initial conditions, except in the three cases mentioned above.

In the case of Kovalevskaya, the problem reduced to quadrature and the integral was obtained through a Remain θ-function of two variables [5] . These functions however are not single-valued, and on the other hand, the functions do not have branch point and hence they are root functions [6] . Consequently, the qualitative and quantitative investigations give us more understanding of the motion of the problem.

2. The Equations of Motion

The Euler-Poisson equations for the Kovalevskaya case where two moments of inertia are equal and the third is half of one of them and the center of mass in the equatorial of the ellipsoid of inertia has the form

$\begin{array}{l}2\frac{\text{d}p}{\text{d}t}=qr,\\ 2\frac{\text{d}q}{\text{d}t}=-pr-{\gamma}_{3},\\ \frac{\text{d}r}{\text{d}t}={\gamma}_{2},\\ \frac{\text{d}{\gamma}_{1}}{\text{d}t}=r{\gamma}_{2}-q{\gamma}_{3},\\ \frac{\text{d}{\gamma}_{2}}{\text{d}t}=p{\gamma}_{3}-r{\gamma}_{1},\\ \frac{\text{d}{\gamma}_{1}}{\text{d}t}=q{\gamma}_{1}-p{\gamma}_{2},\end{array}$ (1)

where $p\mathrm{,}q\mathrm{,}r\mathrm{,}{\gamma}_{1}\mathrm{,}{\gamma}_{2}\mathrm{,}{\gamma}_{3}$ are the Euler-Poisson variables.

The first integrals are taken the following form and for more details [7] .

$2\left({p}^{2}+{q}^{2}\right)+{r}^{2}+2a{\gamma}_{1}=3{l}_{1},$

$2\left(p{\gamma}_{1}+q{\gamma}_{2}\right)+r{\gamma}_{3}=l,\text{'}$

${\gamma}_{1}^{2}+{\gamma}_{2}^{2}+{\gamma}_{3}^{2}=1,$

${\left({p}^{2}-{q}^{2}-a{\gamma}_{1}\right)}^{2}+2{\left(pq-a{\gamma}_{2}\right)}^{2}={k}^{2},$ (2)

and the system (1) can be reduced to quadrature

$\frac{\text{d}{s}_{1}}{\sqrt{\Phi \left({s}_{1}\right)}}+\frac{\text{d}{s}_{2}}{\sqrt{\Phi \left({s}_{2}\right)}}=0,$

$\frac{{s}_{1}\text{d}{s}_{1}}{\sqrt{\Phi \left({s}_{1}\right)}}+\frac{{s}_{2}\text{d}{s}_{2}}{\sqrt{\Phi \left({s}_{2}\right)}}=\frac{\text{d}t}{2},$ (3)

where the Kovalevskaya polynomial $\Phi \left(s\right)$ is defined as

$\Phi \left(s\right)=\left[s\left[{\left(s-6{l}_{1}\right)}^{2}+1-{k}^{2}\right]-2{l}^{2}\right]\left[\left(s-3{l}_{1}-k\right)\left(s-3{l}_{1}+k\right)\right].$ (4)

Kolossoff [8] introduced the potential

$V=\frac{{x}^{2}+{y}^{2}-kx+1}{\sqrt{{x}^{2}+{y}^{2}}},$ (5)

and with the new time $\tau $ through the relation

$i\text{d}\tau =\frac{2\left(r{x}_{1}+{\gamma}_{3}\right)\left(r{x}_{2}+{\gamma}_{3}\right)}{{\left({x}_{1}-{x}_{2}\right)}^{2}}\text{d}t,$ (6)

and then system (1) can be converted to the plan motion of massless point defined as

${x}^{\u2033}=\frac{\partial V}{\partial x},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{y}^{\u2033}=\frac{\partial V}{\partial y},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(\text{'}=\frac{\text{d}}{\text{d}\tau}\right)$ (7)

and possess the energy integral

$\frac{1}{2}\left({\stackrel{\dot{}}{x}}^{2}+{\stackrel{\dot{}}{y}}^{2}\right)+V=h,$ (8)

under the condition that the constants of energy of both systems are equal

$3{l}_{1}=h,$ (9)

Introducing the elliptic coordinates $\lambda \mathrm{,}\mu $ such that

$\begin{array}{l}x=\frac{\lambda \mu}{k}+k,\\ y=\frac{1}{k}\sqrt{\left[\left({\lambda}^{2}-{k}^{2}\right)\left({k}^{2}-{\mu}^{2}\right)\right]},\end{array}$ (10)

with Jacobian matrix $M$ defined as

$M=\left(\begin{array}{cc}\frac{\partial x}{\partial \lambda}& \frac{\partial x}{\partial \mu}\\ \frac{\partial y}{\partial \lambda}& \frac{\partial y}{\partial \mu}\end{array}\right),$ (11)

one can give the momenta ${p}_{x}\mathrm{,}{p}_{y}$ as function of ${p}_{\lambda}\mathrm{,}{p}_{\mu}$ , such that

$\left(\begin{array}{c}{p}_{x}\\ {p}_{y}\end{array}\right)={M}^{-1}\left(\begin{array}{c}{p}_{\lambda}\\ {p}_{\mu}\end{array}\right).$ (12)

Then the canonical variables $\left({p}_{\lambda}\mathrm{,}{p}_{\mu}\mathrm{,}\lambda \mathrm{,}\mu \right)$ are given by

$\begin{array}{l}{p}_{x}=\frac{\left({\lambda}^{2}-{k}^{2}\right)\mu {p}_{\lambda}-\left({\mu}^{2}-{k}^{2}\right)\lambda {p}_{\mu}}{k\left({\lambda}^{2}-{\mu}^{2}\right)},\\ {p}_{y}=\frac{\sqrt{\left[\left({\lambda}^{2}-{k}^{2}\right)\left({k}^{2}-{\mu}^{2}\right)\right]}\left(\lambda {p}_{\lambda}-\mu {p}_{\mu}\right)}{k\left({\lambda}^{2}-{\mu}^{2}\right)}.\end{array}$ (13)

The system is a Liouville system, where the constants of integration are founded from

$\begin{array}{l}\left({\lambda}^{2}-{k}^{2}\right){p}_{\lambda}^{2}+2{\lambda}^{3}+2\left(1-{k}^{2}\right)\lambda -2{\lambda}^{2}h\\ =\left({\mu}^{2}-{k}^{2}\right){p}_{\mu}^{2}+2{\mu}^{3}+2\left(1-{k}^{2}\right)\mu -2{\mu}^{2}h=C,\end{array}$ (14)

and the momenta ${p}_{\lambda}$ and ${p}_{\mu}$ are given by

${p}_{\lambda}=\sqrt{\frac{\Phi \left(\lambda \right)}{{\lambda}^{2}-{k}^{2}}},$

${p}_{\mu}=\sqrt{\frac{\Phi \left(\mu \right)}{{\mu}^{2}-{k}^{2}}}.$ (15)

By using Hamilton-Jacob method, we obtain

$\begin{array}{l}\frac{\text{d}\lambda}{\sqrt{\varphi \left(\lambda \right)\left({\lambda}^{2}-{k}^{2}\right)}}+\frac{\text{d}\mu}{\sqrt{\varphi \left(\mu \right)\left({\mu}^{2}-{k}^{2}\right)}}=0,\\ \frac{\lambda \text{d}\lambda}{\sqrt{\varphi \left(\lambda \right)\left({\lambda}^{2}-{k}^{2}\right)}}+\frac{\mu \text{d}\mu}{\sqrt{\varphi \left(\mu \right)\left({\mu}^{2}-{k}^{2}\right)}}=\frac{\text{d}\tau}{\sqrt{2}\left(\lambda +\mu \right)},\end{array}$ (16)

and with time changing

$\frac{\text{d}\tau}{\sqrt{2}\left(\lambda +\mu \right)}=\text{d}t\mathrm{,}$ (17)

the Equation (16) are the same as the Kovalevskaya equations, where the polynomial

$p\left(u\right)=\left({u}^{2}-{k}^{2}\right)\varphi \left(u\right)\mathrm{,}$ (18)

corresponds to the Kovalevskaya polynomial $\Phi \left(s\right)$ .

Accordingly the problem of the integration of Kovalevskaya top is equivalent to the problem of the plane motion of the massless point, and then we study the problem in Kolossoff’s variable $\left({p}_{\lambda}\mathrm{,}{p}_{\mu}\mathrm{,}\lambda \mathrm{,}\mu \right)$ .

3. Phase Portrait of the Separated Functions

Consider the function

$f={q}^{3}+q\left(1-{k}^{2}\right)+h\left({q}^{2}+1-{k}^{2}\right)-{p}^{2}\left({k}^{2}-{q}^{2}\right).$ (19)

We construct the lines of constant $f$ on the plane $\left(p\mathrm{,}q\right)$ , which is called the phase portrait of $f$ . The phase portrait helps us to find the topological interpretation of the trajectory as follows: if the roots of the function are distinct for given initial values of $\left(p\mathrm{,}q\right)$ , then $P\left({p}_{1}\mathrm{,}{p}_{2}\right)$ and $Q\left({q}_{1}\mathrm{,}{q}_{2}\right)$ will change in a periodic manner, but if the function has multiple roots, then we have an infinity of motion, which gives an asymptotic solution of the canonical equations. Thus the study of lines of constant $f$ provides a complete picture of the bifurcations of our problem.

We first study the singular points and it is distinguished by using Gaussian curvature of $f$ . These points can be found from the equations

$\begin{array}{l}\frac{\partial f}{\partial p}=-p\left({k}^{2}-{q}^{2}\right)=0,\\ \frac{\partial f}{\partial q}=3{q}^{2}+2hq+1-{k}^{2}+2{p}^{2}q=0,\end{array}$ (20)

and hence we have the following: where $p=0$ ,we get

$3{q}^{2}+2hq+1-{k}^{2}=0,$ (21)

and when $q=\pm k$ , we get the two equations

${p}^{2}=-h-k-\frac{1}{2k},$

${p}^{2}=-h+k+\frac{1}{2k},$ (22)

The discriminant of (21) is

${f}_{1}={h}^{2}+3{k}^{2}-3$ (23)

and from (22) the two functions are ${f}_{2}\mathrm{,}{f}_{3}$

$\begin{array}{l}{f}_{2}=-2hk-2{k}^{2}-1,\\ {f}_{3}=-2hk+2{k}^{2}+1.\end{array}$ (24)

It is clear that the curve ${f}_{1}=0$ is tangent to the branches of the curves ${f}_{2}$ and ${f}_{3}$ at the points $k=\pm \frac{1}{2}$ . This can be found from the consideration that

the equation ${f}_{1}={f}_{2}$ gives the two roots $\left({k}_{1}\mathrm{,}{h}_{1}\right)$ , $\left({k}_{2}\mathrm{,}{h}_{2}\right)$ such that $sign\text{\hspace{0.17em}}{k}_{1}\cdot sign\text{\hspace{0.17em}}{k}_{2}<0$ .

From ${f}_{2}=0$ and ${f}_{1}=0$ , we get

$h=-k-\frac{1}{2k}$ (25)

$4{k}^{2}+\frac{1}{4{k}^{2}}-2=0$ (26)

Now, we study the points in the domain ${D}_{i}$ where $i=1,2,\cdots ,16$ , Figure 1.

Figure 1. The regions ${D}_{i}$ of real motions on the $\left(k\mathrm{,}h\right)$ plane.

1) The Domain ${D}_{1}:{f}_{1}<0,\text{\hspace{0.17em}}{f}_{2}<0,\text{\hspace{0.17em}}{f}_{3}>0$

Equations ((21) and (22)) are not solved when $p=0$ , $q=k$ and then there are no singular points of $f$ on the line $p=0$ and $q=k$ , but when $q=-k$ , there are two singular points with $p$ coordinates,

$p=\pm {\left(\frac{2{k}^{2}-2hk+1}{2k}\right)}^{\frac{1}{2}}.$ (27)

To get these types of points, put

$\begin{array}{l}q=-k+y,\\ p=\pm {\left(\frac{2{k}^{2}-2hk+1}{2k}\right)}^{\frac{1}{2}}+x,\end{array}$ (28)

in the function $f$ , neglecting terms of degree greater than 2, then we have

$f=\frac{1-4{k}^{2}}{2k}{y}^{2}\mp 4k{\left(\frac{2{k}^{2}-2hk+1}{2k}\right)}^{\frac{1}{2}}xy+{A}_{0},$ (29)

where ${A}_{0}$ contains the terms of zeros and first degree of $x\mathrm{,}y$ .

The singular points are hyperbolic points, where

${\left[\begin{array}{cc}\frac{{\partial}^{2}f}{\partial {x}^{2}}& \frac{{\partial}^{2}f}{\partial x\partial y}\\ \frac{{\partial}^{2}f}{\partial x\partial y}& \frac{{\partial}^{2}f}{\partial {y}^{2}}\end{array}\right]}_{x=y=0}<0$ (30)

and when $p=0$ , we have the two singular points: $\left(\frac{-h\pm \sqrt{{h}^{2}+3\left({k}^{2}-1\right)}}{3}\mathrm{,0}\right)$ .

In the same manner we can get the type of points in the domains ${D}_{i}\text{\hspace{0.17em}}\left(i=1,\cdots ,16\right)$ , they are shown in Figures 2-20, and summarize that in Table 1.

4. The Integration of the Problem

To complete the picture of our problem, we solve the kolossoff system numerically by using a surface section introducing by Poincaré [6] .

A surface of section is a way of presenting a trajectory in n-dimension phase space in an (n − 1) dimension space. By picking one phase element constant and plotting the value of the other element each time the selected element has the desired value, an intersection surface is obtained. Accordingally the solution of system (7) represents itself the trajectory in the phase space $$ , and on the energy surface h = constant, the trajectory may be treated in the three dimensional invariant hypersurface $$ . The successive points with $$ in at which an orbit crosses the plane $x=0$ in phase space represent a locus of two-dimensional area-preserving mapping on the $$ plane, which is called Poincaré mapping.

Figure 2. The four-hyperbolic points in domain D_{1}.

Figure 3. The Three-hyperbolic points and one-elliptic point in domain D_{2}.

Figure 4. The two-hyperbolic points in domain D_{3}.

Figure 5. The three-hyperbolic points and one elliptic point in domain D_{4}.

Figure 6. The two-elliptic points and two-hyperbolic points in domain D_{5}.

Figure 7. The three-hyperbolic points and one-elliptic point in domain D_{6}.

Figure 8. The two-hyperbolic points and two-elliptic points in domain D_{7}.

Figure 9. The two-hyperbolic points and one-elliptic point in domain D_{8}.

Figure 10. The one-elliptic point in domain D_{8}.

Figure 11. The two-hyperbolic points and two-elliptic points in domain D_{9}.

Figure 12. The one-hyperbolic point and one-elliptic point in domain D_{10}.

Figure 13. The two-elliptic points and two-hyperbolic points in domain D_{11}.

Figure 14. The two-hyperbolic points in domain D_{12}.

Figure 15. The two-elliptic points in domain D_{12}.

Figure 16. The two-hyperbolic points in domain D_{13}.

Figure 17. The one-hyperbolic point and one-elliptic point in domain D_{14}.

Figure 18. The two-hyperbolic points in domain D_{15}.

Figure 19. The one-hyperbolic point in domain D_{16}.

Figure 20. The one-hyperbolic point in domain D_{16}.

If we put $$ in (8), then we have the equipotential lines

$V\left(x,y\right)=h.$ (31)

Figure 21 shows the equipotential lines for different values of h and k. Putting $\stackrel{\dot{}}{x}=0$ in integral energy equation, we have

$\begin{array}{l}\stackrel{\dot{}}{y}=\sqrt{2\left({h}^{}-{V}^{}{\mathrm{(0,y)}}^{}\right)}\end{array}$ (32)

where the curve of zero velocity is

$h-V\left(0,y\right)=0.$ (33)

Figure 22 shows the curve of zero velocity, its clear that the curves is independent of Kovalevsky constant k.

Now the initial values of $\tau ={\tau}_{0}$ is $\left({\mathrm{x\; =}}_{}\stackrel{\dot{}}{x}\mathrm{=0,}{y}_{0}{\stackrel{\dot{}}{y}}_{0}\right)$ and then we get the point ${p}_{1}\left({y}_{1}\mathrm{,}{\stackrel{\dot{}}{y}}_{1}\right)$ at $\tau ={\tau}_{1}$ . The point ${p}_{0}$ is the mapped into ${p}_{1}$ , or ${T}_{{p}_{0}}\to {p}_{1}$ , so if we plot the successive intersections of the motion with the surface of section, they will in general occur anywhere within a bounded area of the plane and must lie on a unique curve.

Hénon-Heiles H-H [9] are used the Poincar ${e}^{\prime}$ section-surface to prove that they exist, the third integral of a non-linear motion of star around a galactic center with the motion restricted to plane, for some values of the constant of energy, the motion is ordered, and with increasing the energy, the ordered motion becomes stochastic. The H-H model [10] becomes a criterion for the existence of the integral of any conservative system. By using these criteria in our problem, we find that, the image points of the initial points, lie on invariant curve. Therefore, these orbits lie on invariant tori and the motion is ordered. Figures 23-28 show the invariant curves for different values of k and h and they are symmetric with respect to y and $$ .

Table 1. The type of points in the domains ${D}_{i}$ .

Figure 21. The equipotential line with a different value of h.

Figure 22. The curve of zero velocity, where $-2\ge h\ge 2$ .

Figure 23. The invariant curves for $h=k=0$ .

Figure 24. The invariant curves for $h<0$ and $k=0$ .

Figure 25. The invariant curves for $0<h\ll 1$ and $0<k\ll 1$ and $h=3k$ .

Figure 26. The invariant curves for $0<h\ll 1$ and $0<k\ll 1$ and $k=5000h$ .

Figure 27. The invariant curves for $0<h\ll 1$ and $0<k\ll 1$ and $h=50k$ .

Figure 28. The invariant curves for $h\ll 0,k<0$ for small value of h and large value of k such that k = 100 times of h, the curves seem to be a closed curves.

5. Conclusions

We conclude some results related to the behavior of the motion:

1) The elliptic points in the figures are stable in the Lyapunov sense, because a small disturbance will result in a closed trajectory that surrounds it and along which the state of the system remains close to these points.

2) The hyperbolic points are unstable because any small disturbance will result in a trajectory on which the state of the system deviates more and more from these points as t goes to infinity.

Cite this paper

El-Sabaa, F.M.F., Mohamed, A.A. and Zakria, S.K. (2017) The Qualitative and Quantitative Methods of Kovalevskys Case. Journal of Applied Mathematics and Physics, 5, 1837-1854. https://doi.org/10.4236/jamp.2017.59155

References

- 1. Leimanis, E. (1965) The General Problem of the Motion of Coupled Rigid Bodies about a Fixed Point. Springer-Verlag, Berlin-Heidelberg, New York. https://doi.org/10.1007/978-3-642-88412-2
- 2. Kovalevskaya, S. (1889) Sur le probléme de la rotation d’un corps solide autour d’un point fixe. Acta Mathematica, 12, 77-232.
- 3. Golubev, V.V. (1953) Lectures on Integration of the Equations of Motion of Rigid Body about a Fixed Point. GITTL, Moscow.
- 4. Kozlov, V.V. (1983) Integrability and Non-Integrability in Hamiltonian Mechanics. Russian Mathematical Surveys, 38, 1-76.
- 5. El-Sabaa, F. (1983) Solution of the Equations of the Problem of Motion of a Heavy Rigid Body about a Fixed Point in the Kovaleveskaya Case Using θ-Function, Celestial Mechanics, 29, 249-253. https://doi.org/10.1007/BF01229138
- 6. Poincaré, H. (1957) Les méthodes nouvelles de la mécanique céleste.
- 7. Fomenko, A.T. (1988) Integrability and Nonintegrability in Geometry and Mechanics. Kluwer Academic Publishers, Dordrecht. https://doi.org/10.1007/978-94-009-3069-8
- 8. Kozlov, V.V. (1976) Splitting of the Separatrices in the Perturbed Euler-Poinsot Problem. Moskovskii Universitet, Vestnik, Seriia I: Matematika, Mekhanika, 31, 99-104.
- 9. Henon, M. and Heiles, C. (1964) The Applicability of the Third Integral of Motion: Some Numerical Experiments. The Astronomical Journal, 69, 73.
- 10. Tabor, M. (1988) Chaos and Integrability in Nonlinear Dynamics. A Wiley-Interscience Publication, New York.