Abstract

In this study, we geometrically analyze the relation between a point vortex system and deviation curvatures on the Jacobi field. First, eigenvalues of deviation curvatures are calculated from relative distances of point vortices in a three point vortex system. Afterward, based on the assumption of self-similarity, time evolutions of eigenvalues of deviation curvatures are shown. The self-similar motions of three point vortices are classified into two types, expansion and collapse, when the relative distances vary monotonously. Then, we find that the eigenvalues of self-similarity are proportional to the inverse fourth power of relative distances. The eigenvalues of the deviation curvatures monotonically convergent to zero for expansion, whereas they monotonically diverge for collapse, which indicates that the strengths of interactions between point vortices related to the time evolution of spatial geometric structure in terms of the deviation curvatures. In particular, for collapse, the collision point becomes a geometric singularity because the eigenvalues of the deviation curvature diverge. These results show that the self-similar motions of point vortices are classified by eigenvalues of the deviation curvature. Further, nonself-similar expansion is numerically analyzed. In this case, the eigenvalues of the deviation curvature are nonmonotonous but converge to zero, suggesting that the motion of the nonself-similar three point vortex system is also classified by eigenvalues of the deviation curvature.

1. Introduction

Geometric theories are often applied to analyze dynamical system behavior. The relationship between a dynamical system and geometry is such that the Euler-Lagrange equation corresponds to the geodesic equation as a second-order differential equation. An invariant theory of a system of second-order differential equations has been established by Kosambi [1], Cartan and Kosambi [2], and Chern [3] (KCC theory). From the KCC theory, there are five geometric invariants. The first invariant is related to the existence of external force. The second invariant is a deviation curvature; eigenvalues of a deviation curvature define the Jacobi stability of whole trajectories. The third invariant is obtained from the differential of a deviation curvature. The fourth invariant is called the Riemann curvature and is associated with deviation curvature. The fifth invariant is called the Douglas tensor.

In previous studies, the theory of the Jacobi field has been applied to various dynamical systems, and concrete phenomena in biology, oscillating system, geophysics, astrophysics, and chaos models have been analyzed using the KCC invariants [428]. In biology, symbiogenesis modeled by the Volterra-Hamilton system has been studied in the Jacobi field [6]. From this analysis, a characteristic of the model has been clarified by the invariants. For the oscillating system, the study of van der Pol has compared the Jacobi and linear stabilities [12, 14], and the study of overhead cranes has investigated the relation between the Jacobi stability and Hopf bifurcation [15]. In geophysics, the Lorenz and Rikitake systems have been studied. The Jacobi stability of the Lorenz system has been analyzed around equilibrium points [17]. Moreover, the nonperiodic behavior of the Lorenz and Rikitake systems has been discussed on the basis of the third invatiant [21, 22]. In astrophysics, the KCC theory has been applied to braneworld models [23]. For chaos models, the Rabinovich-Fabrikant model has been examined from the viewpoint of the KCC theory [27].

In the theory of dynamical systems, it is known that the Hamilton formalism was used to model a behavior of system. From the viewpoint of the Jacobi field theory, geometric structures of the Hamilton system have only been studied in terms of a general form [29]. The relationship between geometric objects and the behavior of specific phenomena remains unclear. As a concrete Hamilton system, the motion of point vortices in fluid dynamics has been studied widely [3036]. Particularly, a two point vortex system can be classified into two types, translation and rigid rotation. In a three point vortex system, chaotic behavior does not develop because the system is integrable [32]. Moreover, four kinds of motion, translation, rigid rotation, expansion, and collapse, are found in three point vortices that satisfy a condition of self-similarity. Here, the self-similarity means that the similarity of initial positions of point vortices is retained whenever. Translation and rigid rotation are classified by circulation conditions, whereas expansion and collapse are determined by circulation conditions and initial positions [30, 33, 35]. In such a point vortex system, equations of motion are expressed by a system of first-order differential equations, which means that the theory of the Jacobi field can be applied to the point vortex system by differentiating the equations of motion. Therefore, geometric objects, such as deviation curvatures can be obtained from the second-order differential equations of the point vortex system. Then, we consider the Jacobi field of the three point vortex system and discuss a classification of motion based on the deviation curvature.

The remainder of this article is structured as follows. In Section 2, we review the KCC theory. In Section 3, we explain the motions of point vortices and their conditions of self-similarity. In Section 4, we calculate eigenvalues of deviation curvatures on self-similarity and discuss a classification of the motions of point vortices by the eigenvalues of the deviation curvature. Finally, Section 5 concludes this study.

2. KCC Theory and Jacobi Stability

In this section, we review the KCC theory and Jacobi stability based on [4, 12, 14, 37, 38]. Notably, the Einstein convention is used throughout this section. Let be a smooth -dimensional manifold, and be a tangent bundle. A projection is a mapping from to . Local coordinates on are denoted by , where . Local coordinates on are defined by , where is the time parameter, and . We consider the following system of second-order differential equations:

where is a smooth function. The first invariant is given by

where is the coefficient of nonlinear connection:

A trajectory of Equation (1) is perturbed to a trajectory , where is a vector field, and is an infinitesimal parameter. From Equation (1), a variation equation is obtained as follows :

Then, the KCC-covariant derivative is defined by

Using (5), Equation (4) is rewritten as follows:

where is the second invariant and is called a deviation curvature tensor

Here, is a coefficient of the Berwald connection given by

The third, fourth, and fifth invariants are defined, respectively, as follows:

The deviation curvature tensor gives the stability of whole trajectories via the following definition [12, 14]: The trajectories of the system (1) are Jacobi stable if and only if the real parts of the eigenvalues of are strictly negative everywhere, and Jacobi unstable otherwise.

3. Review of Point Vortex and Its Self-Similarity

In this section, we review a motion of point vortices that satisfies a self-similarity condition based on previous studies [3034]. The self-similar motion means that the similarity of initial positions of point vortices is retained whenever, namely, the proportion of relative distances are constant. Moreover, we mention the linear stability of a motion of point vortices using relative distances at the end of this section. First, coordinates of three point vortices are denoted by , , and . Their respective relative distances are defined as follows:

Generally, equations of motions of point vortices are obtained from complex velocity potential. From [32], the equations of motion of three point vortices are given by

Using Equations (11), (12), and (13), the motions of the three point vortices are shown on a two-dimensional plane. Further, to describe a three point vortex system as Hamilton’s canonical form, the Hamiltonian of three point vortices is defined as follows:

where , , and are circulations of first, second, and third point vortices, respectively. Thus, from (11), (12), (13), and (14), equations of motions of point vortices are expressed as the following Hamilton’s canonical form:

where , and each left hand sides of Equation (15) are not summed. For all motions of point vortices, the Hamiltonian is constant. Then, Equation (15) is rewritten in terms of the relative distances , , and as follows:

where is given by

is an area of the triangle formed by the three point vortices. Moreover, Equations (16), (17), and (18) can be, respectively, rewritten as follows:

Here, all initial values of , , and should be positive because , , and are relative distances. Subsequently, we use Equations (20), (21), and (22) to analyze geometric quantities of point vortices.

Next, we review the self-similarity of the three point vortices based on previous studies [30, 33, 35]. When the initial positions of the three point vortices are an equilateral triangle, all relative distances are constant since all right hand sides of Equations (20), (21), and (22) become zero; thus, the Hamiltonian (14) is constant even without circulation conditions. Then, the motions of point vortices are translation or rigid rotation. These two kinds of motions are classified by the circulation conditions. From [30], for translation, the circulations satisfy

Meanwhile, from [30], the motion of rigid rotation satisfies

When the initial positions of the three point vortices are not an equilateral triangle, all self-similar motions of the point vortex system satisfy which is because the Hamiltonian has to be constant for a motion of point vortices. Moreover, the initial position of the third point vortex of all self-similar motions satisfies

whereas the first and second point vortices lie on the -axis . Here, is an angle between the -axis and the line segment from the center of the circle to the position of the third point vortex (Figure 1), where the point of the center of circle is given by

Finally, we mention the linear stability of a motion of point vortices that satisfies the conditions (26) and (27) using relative distances and its Jacobian matrix. The equilibrium point is which all of left hand sides of Equations (20), (21), and (22) equal zero. Thus, for the equilibria, the relative distances satisfy the following:

where is relative distances and . From Equations (20), (21), and (22), the Jacobian matrix at the equilibrium point is given by

The eigenvalues of (30) is

where , , and are the first, second, and third eigenvalues of the matrix (30). From (31), since , the linear stability around the equilibrium point of the three point vortex system is not asymptotically stable but is Lyapunov stable when .

In this article, we handle expansion and collapse as the self-similar motion of point vortices.

4. Relationship between Three Point Vortices and a Deviation Curvature

In this section, a motion of three point vortices is discussed from the viewpoint of differential geometry. The KCC theory is applied to the point vortex system, and the deviation curvature is obtained. Then, the expansion and collapse of point vortices are classified by the deviation curvature for both self- and nonself-similar motions.

4.1. Eigenvalues of Deviation Curvature of a Three Point Vortex System

First, we calculate eigenvalues of deviation curvatures in a three point vortex system that exhibits general motions including self- and nonself-similarity. For applying the KCC theory, Equations (20), (21), and (22) are differentiated by time parameter as follows:

Here, the initial values of Equations (32), (33), and (34) have to be positive, and are determined by Equations (20), (21), and (22). From Equations (32), (33), and (34), the function , nonlinear connection and Berwald connection can be obtained (Appendix A). Then, a deviation curvature of the three point vortex system becomes a matrix , and the characteristic equation is written as follows:

where is an eigenvalue of . Let , , , and be the first, second, third, and th eigenvalues of a deviation curvature. From the characteristic equation (35) and Cardano’s formula of a cubic equation [39], the th eigenvalue of the deviation curvature of a general motion system including self- and nonself-similarity is given by where

The condition in which the eigenvalues (36) become double or triple repeated roots is given by discriminant :

When , the eigenvalues (36) become double or triple repeated roots. For and , the eigenvalues (36) become triple repeated roots.

4.2. Classification of Self-Similar Motions of a Three Point Vortex System by Deviation Curvature

In this section, we consider a motion of self-similarity of point vortices based on [30, 33, 35]:

where , , and are constants that are determined by the initial positions of point vortices. From [33], the particular form of the function is given by

where is the collapse time if , and if , the function increases monotonically.

For simplicity, using (25), let be as follows:

Then, , , and can be expressed as follows: where the details of (46) are described in Appendix B.

Using (41) and (46), , , and can be rewritten as follows:

where , , and are determined by , , , , and . Therefore, from (36) and (47), eigenvalues of a deviation curvature of self-similar motions can be rewritten as follows:

where , , and are determined by , , , , and ; thus, the term is constant that means independent of the time . From (48), eigenvalues of self-similarity are proportional to the inverse fourth power of relative distances. Accordingly, since of (42) is monotonous, for expansion, the eigenvalues of the deviation curvature monotonically converge to zero, whereas for collapse, the eigenvalues of the deviation curvature monotonically diverge. Thus, when the strengths of interactions between vortices decrease, the relative distances are increase, that is expansion, and eigenvalues of the deviation curvature monotonically converge to zero from (48). When the strengths of interactions between vortices increase, relative distances decrease, that is collapse, and eigenvalues of the deviation curvature monotonically diverge from (48). These indicate that the strengths of interactions between vortices are related to the time evolution of spatial geometric structure in terms of the deviation curvature. In particular, for collapse, the collision point becomes a geometric singularity because the eigenvalues of the deviation curvature are divergence. These results mean that motions of self-similarity are classified by the eigenvalues of the deviation curvature.

Next, we consider specific examples of expansion and collapse based on the discriminant . When the third point vortex is on the circle of Figure 1, using , and , the discriminant (40) is zero ( =0), and of (39) is given by

According to [30, 33, 35], when , , , and , the motions of point vortices are not expansion and collapse. This means for expansion and collapse. Therefore, for expansion and collapse, the eigenvalues (36) always have double repeated roots. This result suggests that eigenvalues of a deviation curvature are affected by the assumption of self-similarity, which implies a decrease in degrees of freedom. The results for expansion and collapse can be confirmed numerically. Initial conditions and circulations for the expansion are given by where expression (53) satisfies conditions (26) and (27), and . Moreover, initial position and circulation conditions for the collapse are given by where (58) satisfies conditions (26) and (27), and . Then, variations in relative distances and time evolutions of the deviation curvature are shown in Figures 2 and 3. The relative distances , , and are represented by black, blue, and red, respectively. For the eigenvalues of the deviation curvature, double repeated roots are represented by a solid line, and the other roots are represented by a dashed line. For expansion in Figure 2, the eigenvalues of the deviation curvature monotonically converge to zero, whereas for collapse in Figure 3, the eigenvalues of the deviation curvature monotonically diverge. The other KCC invariants are shown in Appendix C.

4.3. Numerical Example of a Nonself-Similar Motion

In this section, the variations in relative distances and time evolutions of eigenvalues of the deviation curvature are analyzed numerically for a nonself-similar expansion. The initial positions and circulations are given by

The initial position (61) does not satisfy the self-similar condition shown in [30, 33, 35]. Figure 4 shows the variation of relative distances and time evolutions of eigenvalues of the deviation curvature. The relative distances , , and are represented by black, blue, and red, respectively. In this case, the eigenvalues of the deviation curvature are nonmonotonous but converge to zero. This result suggests that the motions of a nonself-similarity system of three point vortices can be classified by eigenvalues of the deviation curvature as well as motions of a self-similarity system.

5. Conclusions

The KCC theory has been applied to a three point vortex system, and eigenvalues of a deviation curvature have been calculated using relative distances. For self-similar motions, the eigenvalues of deviation curvature are proportional to the inverse fourth power of relative distances. We consider the classification of self-similar motion of point vortices by eigenvalues of the deviation curvature. The details of the classification are as follows. For expansion, the eigenvalues are double repeated roots, and the eigenvalues of the deviation curvature monotonically converge to zero. For collapse, the eigenvalues are double repeated roots, and the eigenvalues of the deviation curvature monotonically diverge. These results mean that self-similar motions are classified by eigenvalues of the deviation curvature. Further, in the case of nonself-similar expansion, eigenvalues of the deviation curvature are not monotonous but convergent to zero. This numerical result suggests that motions of a nonself-similar system of three point vortices can be classified by eigenvalues of a deviation curvature as well as motions of a self-similarity system.

Appendix

A. Details of Nonlinear Connection and Berwald Connection

In this section, we show the details of the function , nonlinear connection , and Berwald connection mentioned in Section 4. From Equation (1), is rewritten as :

From expressions (A.1), (A.2), and (A.3), coefficients of nonlinear connection are given by where

Then, from expressions (A.4)–(A.12), , , and are not included in , therefore .

B. Details of Trace and Determinant of Deviation Curvature

In this section, we describe in detail , , and on the assumptions of self-similarity and (43) mentioned in Section 4, , , and are given by where

C. Details of Other KCC invariants

In this section, we describe in detail the first, third, fourth, and fifth KCC invariants on the assumptions of self-similarity, and . The components of the first invariant are given by

Moreover, all components of the third, fourth, and fifth invariants are zero:

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.