Next Article in Journal
Implicit-Explicit Methods for a Convection-Diffusion-Reaction Model of the Propagation of Forest Fires
Next Article in Special Issue
Mathematical Modeling of Vortex Interaction Using a Three-Layer Quasigeostrophic Model. Part 1: Point-Vortex Approach
Previous Article in Journal
The Square-Zero Basis of Matrix Lie Algebras
Previous Article in Special Issue
Sub-Mesoscale Frontal Instabilities in the Omani Coastal Current
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Effects of Circulation around a Circle on the Stability of a Thomson Vortex N-gon

by
Leonid Kurakin
1,2,3,*,† and
Irina Ostrovskaya
3,†
1
Water Problems Institute, RAS, 119333 Moscow, Russia
2
Southern Mathematical Institute, Vladikavkaz Scientific Center of RAS, 362027 Vladikavkaz, Russia
3
Institute for Mathematics, Mechanics and Computer Sciences, Southern Federal University, Milchakova st., 8a, 344090 Rostov on Don, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2020, 8(6), 1033; https://doi.org/10.3390/math8061033
Submission received: 19 May 2020 / Revised: 18 June 2020 / Accepted: 22 June 2020 / Published: 24 June 2020
(This article belongs to the Special Issue Vortex Dynamics: Theory and Application to Geophysical Flows)

Abstract

:
The stability problem of the stationary rotation of N identical point vortices is considered. The vortices are located on a circle of radius R 0 at the vertices of a regular N-gon outside a circle of radius R. The circulation Γ around the circle is arbitrary. The problem has three parameters N, q, Γ , where q = R 2 / R 0 2 . This old problem of vortex dynamics is posed by Havelock (1931) and is a generalization of the Kelvin problem (1878) on the stability of a regular vortex polygon (Thomson N-gon) on the plane. In the case of Γ = 0 , the problem has already been solved: in the linear setting by Havelock, and in the nonlinear setting in the series of our papers. The contribution of this work to the solution of the problem consists in the analysis of the case of non-zero circulation Γ 0 . The linearization matrix and the quadratic part of the Hamiltonian are studied for all possible parameter values. Conditions for orbital stability and instability in the nonlinear setting are found. The parameter areas are specified where linear stability occurs and nonlinear analysis is required. The nonlinear stability theory of equilibria of Hamiltonian systems in resonant cases is applied. Two resonances that lead to instability in the nonlinear setting are found and investigated, although stability occurs in the linear approximation. All the results obtained are consistent with those known for Γ = 0 . This research is a necessary step in solving similar problems for the case of a moving circular cylinder, a model of vortices inside an annulus, and others.

1. Introduction

The classical hydrodynamics of an ideal fluid uses a model of point vortices (or parallel vortex threads) that goes back to H. Helmholtz. Let’s list the reviews that contain various tasks that use this model [1,2,3,4,5,6,7]. G. Kirchhoff obtained the Hamiltonian form of the motion equations of point vortices on the plane and indicated possible first integrals. These integrals were sufficient for the complete integrability of the problem in the case of two and three vortices.
The point vortices model was used by Kelvin (1878) in the construction of his theory of the atom. Although this theory was rejected the mathematical model survived. Now it took on new urgency in connection with the theory of vortices in superfluid [8] and the electronic columns [9]. Kelvin posed the stability problem for the stationary rotation of the system of N equal vortices located in the vertices of regular N-gon (Thomson N-gon). A long history of research on this problem is described in [10], where a review of the experimental and theoretical work is given. The stability analysis in linear approximation was done by J. Thomson (1883) and T. Havelock (1931). Nonlinear stability analysis was performed much later [10].
From the general papers on the movement of point vortices in a fluid bounded by absolutely smooth walls, it should be noted the Lin’s paper [11], in which, in particular, it is shown that the motion equations of vortices inside and outside the circle are Hamiltonian. The model of vortices outside the circle takes into account that the circulation Γ around the circle can be arbitrary one [11,12]. The case of Γ = 0 is of particular interest because when the radius of the circle tends to zero, we come to the Kirchhoff’s model on the plane.
Havelock investigated the stability problem of a regular vortex polygon inside and outside the circle in the case of Γ = 0 in a linear setting [13]. He obtained conditions under which the corresponding linearized system has exponentially growing solutions. Thus, after this paper, the question of stability in the case of Γ = 0 remained open when all the eigenvalues of the linearization matrix are purely imaginary or null. Nonlinear analysis, in this case, is performed for the problem inside the circle in the papers [14,15] and for the problem outside the circle in the series of the papers [16,17,18,19]. The general theory of nonlinear stability analysis of solutions of Hamiltonian systems by methods of Kolmogorov–Arnold–Moser theory (KAM-theory) was used [20,21,22,23]. In particular, the quadratic part of the Hamiltonian was studied. Resonant cases requiring analysis of the third and fourth-order terms of the Taylor series expansion of the Hamiltonian were considered separately. It turned out that two of them lead to instability, although in the linear setting there is stability.
The purpose of this paper is to study the stability of a regular vortex N-gon in the case of an arbitrary circulation Γ 0 . The linearized system and the quadratic part of the Hamiltonian are investigated. This allowed us to divide the parameter space of the problem into three areas: the areas of orbital stability and instability in the nonlinear setting, and the area where linear stability takes place and nonlinear analysis is required.
In this paper, the stability theory of Hamiltonian systems in resonant cases has been used [23,24]. We applied it to study the influence of circulation Γ on stability in two resonant cases, the instability of which was found at Γ = 0 in [17,18]. These are critical case of the double-zero resonance at N = 3 and the case of resonance 1 : 2 at N = 5 .
The model of point vortices outside the circular domain considered by us is a limiting case of more complex models, for example, the movement of discrete vortices in an annular, in a deformation flow [25], the model of moving cylinder interacting with the configuration of N point vortices [26]. Therefore, our research is a necessary step in the analysis of similar problems in all these models.

2. Problem Formulation

The motion of the system N identical point vortices of intensity ϰ outside the circle of radius R is described by the system (see, for example [12])
z ¯ ˙ k = ϰ 2 π i j k j = 1 N 1 z k z j ϰ 2 π i j = 1 N 1 z k z ^ j + 1 2 π i 1 z k Γ + N ϰ , k = 1 , , N .
Here z k = x k + i y k , k = 1 , , N are complex variables, x k , y k are Cartesian coordinates of kth vortex, z ^ k = R 2 z ¯ k is the reflection of the kth vortex from the boundary of the circle, parameter Γ is circulation around circle.
Let us assume that ϰ = 1 .
The system (1) can be written as
z ¯ ˙ k = 2 i H ( z , z ¯ ) z k , z ˙ k = 2 i H ( z , z ¯ ) z ¯ k , k = 1 , , N ,
where z = ( z 1 , , z N ) , H = H ( z , z ¯ ) is a Hamiltonian:
H = 1 4 π 1 j < k N ln z j z k 2 + 1 8 π j , k = 1 N ln R 2 z j z ¯ k 2 1 4 π Γ + N k = 1 N ln z k 2 .
In the case, Γ = 0 and R 0 the system of Equation (1) transforms to the Kirchhoff equations of motion of point vortices on the plane.
The system (2) and (3) has an exact solution
z k = e i Ω N t u k , u k = R 0 e 2 π i ( k 1 ) / N , k = 1 , , N ,
Ω N = 1 4 π R 0 2 3 N 1 2 N 1 q N + 2 Γ , q = R 2 R 0 2 < 1 .
Thus, a system of N identical vortices located on a circle of radius R 0 at the vertices of a regular N-gon rotates with a constant angular velocity Ω N ( q , Γ ) (see Figure 1).
The angular velocity Ω N ( q , Γ ) decreases monotonously over q ( 0 , 1 ) at fixed values Γ . The vortex N-gon changes the rotation direction at a fixed Γ > 3 N 1 2 with the growth of parameter q on the curve
Γ N ( q ) = 1 2 3 N 1 2 N 1 q N .
In the case, Γ < 3 N 1 2 the rotation direction with the growth of q does not change. The graphs of Γ N ( q ) are shown in Figure 2.
The problem of stability of stationary rotation (4) is studied in the following sections.

3. Stability of a Regular Vortex N -gon for Γ = 0

We consider the case when circulation around the boundary is zero: Γ = 0 .
In this case, the stability problem of the stationary rotation of the Thomson vortex N-gon was solved in the linear setting by Havelock [13]. He found that the corresponding linearized system has exponentially growing solutions for N = 2 , , 6 , if the parameter q = R 2 / R 0 2 is greater than the critical value q 1 N : q 1 N < q < 1 , and for N 7 the instability occurs for all q ( 0 , 1 ) .
In the papers [17,18] for odd N = 3 , 5 and in [16] for even N = 2 , 4 , 6 , it is shown with the help of analysis of the quadratic part of the Hamiltonian that for q ( 0 , q 0 N ) there is orbital stability in the exact nonlinear setting. In the case of even N the value q 0 N = q 1 N . The critical values of q 0 N and q 1 N are shown in the Table 1.
Nonlinear analysis was performed using KAM-theory methods (see, for example [23,24]). It was held at a point q = q 1 N in the case of even N = 2 , 4 , 6 [16] and on the segment q [ q 0 N , q 1 N ] for odd N = 3 , 5 [17,18]. In particular, all resonances encountered in the problem were found and studied. It turned out that only two of them lead to instability. These are the resonance of the double zero at N = 3 , if q = q 03 , and the resonance 1 : 2 at N = 5 in the point q = q 1 : 2 = 0.333377 .

4. Linear Analysis of Stability of the Vortex N -gon in the Case of Arbitrary Γ

4.1. The Perturbations Equation

The change of variables
z k ( t ) = e i Ω N t v k ( t ) , v k = R 0 1 + 2 r k e i ( 2 π N ( k 1 ) + θ k ) , k = 1 , 2 , , N
and scaling of time t t / R 0 2 transforms the system (2) and (3) to the perturbations equations
r ˙ k = E θ k ( ρ ) , θ ˙ k = E r k ( ρ ) , ρ = ( r 1 , , r N , θ 1 , , θ N ) R 2 N ,
with the Hamiltonian
E ( ρ ) = H ( z ( ρ ) , z ¯ ( ρ ) ) + Ω N R 0 2 2 k = 1 N ( 1 + 2 r k ) .
The Taylor series expansion of the Hamiltonian E ( ρ ) in powers of ρ in the neighborhood of the zero solution is
E ( ρ ) = 1 4 π ( E 0 + E 2 ( ρ ) + E 3 ( ρ ) + ) .
The terms of degree higher than three are omitted. The linear terms is zero: E 1 = 0 . The quadratic form E 2 can be represented as
E 2 = S ρ , ρ , S = F 1 1 2 G 0 1 2 G 0 F 2 ,
where · , · is the scalar product of two vectors.
The linearization matrix of the system (8) in the neighborhood of the zero solution has the form
L = G 0 2 F 2 2 F 1 G 0 .
The matrices F 1 , F 2 , G 0 are circulants of size N × N , moreover F 1 , F 2 are symmetric matrices, and G 0 is a skew-symmetric one. They are polynomials of the cyclic matrix C (see, for example Appendix A of [10]):
F m = d e f f m 0 I + j = 1 N 1 f m j C j , m = 1 , 2 , G 0 = d e f j = 1 N 1 g 0 j C j ,
C = 0 1 0 0 0 0 0 1 0 0 0 0 0 0 1 1 0 0 0 0 ,
where I is unit matrix.
Only the coefficient f 10 depends on the parameter Γ :
f 10 ( q , Γ ) = f 10 0 ( q ) + 2 Γ ,
where
f 10 0 ( q ) = N 2 1 12 + ( N 1 ) 2 N q N 1 q N N 2 q N ( 1 q N ) 2 q ( 1 q ) 2 .
The coefficients f 10 0 , f 1 j , f 20 , f 2 j and g 0 j , j = 1 , , N 1 do not depend on Γ , and therefore, they match the corresponding coefficients written out in the papers [13] for the case of Γ = 0 :
f 1 j ( q ) = 1 2 1 1 cos 2 π j N q ( ( 1 + q 2 ) cos 2 π j N 2 q ) ( 1 2 q cos 2 π j N + q 2 ) 2 , f 20 ( q ) = N 2 1 12 + N 2 q N ( 1 q N ) 2 q ( 1 q ) 2 , f 2 j ( q ) = 1 2 1 1 cos 2 π j N q ( ( 1 + q 2 ) cos 2 π j N 2 q ) ( 1 2 q cos 2 π j N + q 2 ) 2 , g 0 j ( q ) = 2 q ( 1 q 2 ) sin 2 π j N ( 1 2 q cos 2 π j N + q 2 ) 2 .
The eigenvalues λ m k , m = 1 , 2 and i λ 0 k of matrices F 1 , F 2 and G 0 are calculated with the help of the formulas
λ m k = j = 0 N 1 f m j e i 2 π j k N , m = 1 , 2 , i λ 0 k = j = 1 N 1 g 0 j e i 2 π j k N , k = 1 , , N .
The value λ 1 k has the form
λ 1 k ( q , Γ ) = λ 1 k 0 ( q ) + 2 Γ ,
where
λ 1 k 0 ( q ) = k ( N k ) 2 + ( N 1 ) 2 N q N 1 q N N 2 q N k ( 1 + q k ) 2 2 ( 1 q N ) 2 N k ( q k q N k ) 2 ( 1 q N ) .
The values λ 1 k 0 and the eigenvalues λ 2 k , i λ 0 k are the same as in the papers [13]:
λ 2 k ( q ) = 1 2 k ( N k ) N k ( q k q N k ) 2 ( 1 q N ) N 2 q N k ( 1 q k ) 2 2 ( 1 q N ) 2 , λ 0 k ( q ) = N k ( q k + q N k ) 1 q N + N 2 q N k ( 1 q 2 k ) ( 1 q N ) 2 .
The matrices F 1 , F 2 have the same eigenvectors
h = ( 1 , cos ( α ) , , cos ( ( N 1 ) α ) ) T , h N = ( 0 , sin ( α ) , , sin ( ( N 1 ) α ) ) T , h N = ( 1 , 1 , , 1 ) T , α = 2 π N , = 1 , , N 1 2 .
In the case of even N = 2 n , there is another eigenvector
h N 2 = ( 1 , 1 , , 1 , 1 ) T .
For the matrices F 1 , F 2 and G 0 , the relations are performed:
F m h k = λ m k h k , m = 1 , 2 , k = 1 , , N , G 0 h k = λ 0 k h N k , G 0 h N k = λ 0 h k , G 0 h N = 0 .
The eigenvalues of matrix S are the roots of the polynomials (see [14], Section 3.1 of [27]):
P ( N , k , Λ ) = Λ 2 ( λ 1 k + λ 2 k ) Λ + λ 1 k λ 2 k 1 4 λ 0 k 2 , k = 1 , , N .
The eigenvalues of the linearization matrix (12) are calculated with the help of the formulas [13]
σ k ± = i λ 0 k ± 2 λ 1 k λ 2 k , k = 1 , , N .

4.2. Orbital Stability and Instability of Two-Dimension Invariant Set: General Definitions and Theorems

The stationary rotation regime (4) is Lyapunov unstable for any N 2 . Indeed, the angular speed of rotation Ω N depends on the radius R 0 . Therefore, for an arbitrarily small perturbation of the radius R 0 , the perturbed Thomson N-gon will move away from the undisturbed one by an order of diameter 2 R 0 . Thus, we need to clarify our definitions of stability and instability.
In this paper, as in [27], stability of a stationary solution (4) is understood as orbital stability or stability of the continuous family of equilibria
C = { ρ R 2 N : r 1 = = r N = 0 , θ 1 = = θ N } .
of the perturbations Equations (8) and (9).
The instability is instability of the invariant set C 2 of the system (8) and (9):
C 2 = { ρ R 2 N : r 1 = = r N , θ 1 = = θ N } ,
formed by all the orbits of a continuous family (4) depending on the parameter R 0 > R (For more information, see Section 3 of [27]).
The matrices S and L have the same structure as in [27]. Therefore, repeating verbatim the arguments about the orbital stability of the equilibrium family C and the instability of the invariant set C 2 , given in Sections 3.2 and 3.3 of [27], we obtain the following statements.
Proposition 1
(Section 3.2 of [27]). Let roots of the polynomials P ( N , k , Λ ) , k = 1 , , N 2 and the value λ 1 N , given by Formulas (23) and (18)–(20), have identical sign. Then,
the eigenvalues of the matrix S , given by (23), have identical sign except the simple zero;
the quadratic form E 2 = S ρ , ρ and the Hamiltonian (9) reach an extremum on the equilibria family (25);
the stationary rotation (4) and (5) N-gon of the system (2) with the Hamiltonian (3) is orbitally stable.
Proposition 2
(Section 3.3 of [27]). Let the inequality
λ 1 k λ 2 k < 0 ,
be satisfied for at least one k = 1 , , N 2 . The values λ 1 k and λ 2 k are given by Formulas (18)–(20).
Then the linearization matrix L has at least one eigenvalue (24) with a positive real part and the invariant set C 2 is unstable.
Proposition 3
(Section 3.3 of [27]). Let inequalities
λ 1 k λ 2 k 0
are valid for all k = 1 , , N 2 , then all eigenvalues of the matrix L lie on the imaginary axis, and stability by linear approximation takes place.
If, in addition, the conditions of Proposition 1 are violated, then the linear approximation is not enough to analyze the stability.
According to the properties 1 ° and 3 ° of the Proposition A1 of Appendix A, the value λ 2 k is positive for any N and k = 1 , , N 1 :
λ 2 k ( q ) > 0 , 0 < q < 1
and the value
λ 1 N 0 ( q ) > λ 1 N 2 0 ( q ) , 0 < q < 1 .
Taking into account these inequalities, we reformulate the Propositions 1–3 as follows.
Theorem 1.
The stationary rotation (4) of the Thomson vortex N-gon is
orbitally stable in an exact nonlinear setting if the inequalities
λ 1 k λ 2 k λ 0 k 2 4 = ( λ 1 k 0 + 2 Γ ) λ 2 k λ 0 k 2 4 > 0 , k = 1 , , N 2 ;
are valid;
unstable if there is at least one value k = 1 , , N 2 such that
λ 1 k = λ 1 k 0 + 2 Γ < 0 .
If conditions 1 ° and 2 ° are not valid, then the solution (4) is stable in a linear setting, and nonlinear analysis is required to solve the stability problem.

4.3. Stability Diagram of a Regular Vortex N-gon

Let’s introduce the following notation
Γ ˜ 0 k = d e f 1 2 λ 1 k 0 + λ 0 k 2 8 λ 2 k .
It is numerically established at N 100 that the following conditions
min 1 k N 2 λ 1 k 0 = λ 1 N 2 0 , 0 < q < 1 ,
max 1 k N 2 Γ ˜ 0 k ( q ) = Γ ˜ 0 N 2
are valid.
Conditions (33) and (34) have been checked for a number of values N > 100 and no violations of them have been found anywhere. It is possible that these conditions are also true for arbitrary N.
In similar stability problems for other models of vortex dynamics the conditions (29), (33) and (34) may not be valid, for example, for a two-fluid plasma model [28].
Let the values Γ 0 N and Γ 1 N be given by the formulas:
Γ 1 N ( q ) = 1 2 λ 1 N 2 0 , Γ 0 N ( q ) = 1 2 λ 1 0 + λ 0 2 8 λ 2 , in the odd case N = 2 + 1 ; Γ 1 N ( q ) , in the even case N = 2 n .
If the relations (33) and (34) are fulfilled, then the condition (31) is
Γ > Γ 0 N ,
and the condition (32) has the form
Γ < Γ 1 N .
Applying the Theorem 1 leads to results on the stability of stationary rotation of a Thomson vortex N-gon, schematically shown in Figure 3. The curves Γ 0 N and Γ 1 N divide the parameter space of the problem into the areas with different types of stability. Gray color indicates the area where condition 1 ° of the Theorem 1 is valid, the vertical shading is area, where the condition 2 ° is valid, and white color is one, where both conditions 1 ° and 2 ° are not valid.
The value Γ N is
Γ N = Γ 0 N ( 0 ) = Γ 1 N ( 0 ) = N 2 8 N + 7 16 , in the case of odd N = 2 + 1 ; N 2 8 N + 8 16 , in the case of even N = 2 n .
The values of Γ N for N = 2 , , 10 are shown in the Table 2. For N < 7 the values Γ N < 0 .
In the case of odd N = 2 + 1 , a linear stability domain, which requires nonlinear analysis, exists for all q ( 0 , 1 ) , since according to the property 5 ° of the Proposition A1, the inequality Γ 1 N ( q ) < Γ 0 N ( q ) is valid.
Note that instability occurs for all values of 0 < q < 1 , if Γ < Γ N . If Γ > Γ N , all three areas always exist. In the case of even N, the domain, where nonlinear analysis is required, consists only of a curve Γ 1 N . If Γ = 0 , the stability interval exists only for N < 7 [13]. In the case N = 7 the linear stability domain consists of one point q = 0 , that corresponds the problem of stability of a Thomson heptagon on a plane [10].
According to the general theory (see, for example [23,24]) in the white domain in the Figure 3 and in the neutral curves Γ 0 N and Γ 1 N instability is possible when the eigenvalues of the linearization matrix (12) are in resonant relations. Nonlinear terms of the Taylor series of the Hamiltonian (9), starting from the third degree, must satisfy some additional conditions. In [17,18] for Γ = 0 in the segment q 0 N q q 1 N two resonances leading to instability were found in the case N = 3 , 5 : in the point q = q 03 = 0.262542 in the case N = 3 and in the point q = q 5 = 0.333377 in the case N = 5 . These are the double-zero resonance (diagonalizable case) and the 1:2 resonance, respectively. In the following sections, they are studied for arbitrary values Γ .

5. The Double-Zero Resonance in the Stability Problem of the Thomson Vortex Triangle

Let N = 3 and a point ( q , Γ ) lies on the resonance curve Γ = Γ 03 ( q ) given by Equation (35). There is a four-fold zero eigenvalue among the eigenvalues (24) of the linearization matrix L . In addition to σ 3 ± = 0 , two other null eigenvalues appear: σ 2 + and σ 1 . In this section, it will be shown that in this case the invariant set C 2 given by the Formula (26) is unstable. To prove it, we construct a reduced system following the paper [17].
We introduce the change of variables [17,23]
ρ = A 3 ξ ζ , ξ = ( ξ 1 , ξ 2 , ξ 3 ) T , ζ = ( ζ 1 , ζ 2 , ζ 3 ) T ,
where A 3 is a symplectic matrix [29]:
A 3 = 1 3 ν h 2 ν h 2 h 3 2 λ 13 ν h 1 ν h 1 h 0 h 1 ν h 1 ν h 0 h 2 ν h 2 ν 2 λ 13 h 3
or
A 3 = A 31 , A 32 , A 33 , A 34 , A 35 , A 36 .
Here, A 3 k R 6 is k-th column of matrix A 3 , h k , k = 1 , 2 , 3 are basis vectors defined by equalities (21) for N = 3 :
h 1 = ( 1 , 1 2 , 1 2 ) T , h 2 = ( 0 , 3 2 , 3 2 ) T , h 3 = ( 1 , 1 , 1 ) T ,
h 0 is zero vector, ν = | λ 21 / λ 11 | 4 . According to the Formula (35) and the Proposition A1, values λ 11 ( q , Γ 03 ( q ) ) and λ 13 ( q , Γ 03 ( q ) ) are positive since Γ 03 > Γ 13 .
In new variables ( ξ , ζ ) the Hamiltonian E ( ρ ( ξ , ζ ) ) given by the Formula (9) has a cyclic variable ζ 3 , since the Hamiltonian (9) is invariant with respect to change of variables θ k θ k + η , η R , which in the vector form is written as
ρ = ρ + 3 2 λ 13 η A 36 .
The vectors A 36 and A 33 are the eigenvector and adjoined vector of a linearization matrix L . They correspond to the double zero eigenvalue σ 3 ± = 0 :
L A 36 = 0 , L A 33 = A 36 .
The cyclic variable ζ 3 and the variable ξ 3 correspond to the vectors A 36 and A 33 , respectively.
Methods for investigating equilibrium stability for systems with cyclic variables were developed by Routh [30]. Assuming ξ 3 = 0 , we obtain a reduced system with the reduced Hamiltonian
W ( ξ 1 , ξ 2 , ζ 1 , ζ 2 ) = E ( ρ ( ξ 1 , ξ 2 , 0 , ζ 1 , ζ 2 , 0 ) ) .
The invariant set C 2 is the subspace formed by the vectors A 36 and A 33 . A reduced system with the Hamiltonian (43) is defined in a subspace orthogonal to the invariant set C 2 . Hence, the instability of the zero equilibrium of a reduced system is the instability of the invariant set (26).
The Taylor series expansion of Hamiltonian (43) in the resonant curve Γ = Γ 03 ( q ) has the form
W ( ξ 1 , ξ 2 , ζ 1 , ζ 2 ) = 1 4 π ( W 0 + W 2 + W 3 + ) ,
W 2 = ω 1 ( q ) ( ξ 1 2 + ζ 1 2 ) , ω 1 = 1 2 σ 1 + ( q , Γ 03 ( q ) ) .
The terms of degree higher than three are omitted.
According to [31], the zero equilibrium is unstable if ω 1 ( q ) 0 and the third-degree form W 3 is different from the identity zero on the set of zeros of the quadratic part W 2 .
In the Figure 4 it is shown that coefficient ω 1 ( q ) is positive for all 0 < q < 1 . The quadratic form W 2 vanishes on the set
{ ( ξ 1 , ξ 2 , ζ 1 , ζ 2 ) R 4 | ξ 1 = ζ 1 = 0 } .
The calculations showed that the third-degree form W 3 on this set has the form
W 3 ( 0 , ξ 2 , 0 , ζ 2 ) = a 3 ( q ) ( ζ 2 3 3 ζ 2 ξ 2 2 )
As you can see from Figure 5 the coefficient a 3 ( q ) is not zero.
Thus, the stationary rotation (4) of the Thomson vortex triangle is unstable at all points of the resonance curve Γ 03 ( q ) at the interval 0 < q < 1 .
The function Γ 03 ( q ) turns to zero at the point q = q 03 = 0.262542 , that is consistent with the result of [17].

6. The Resonance 1:2 in the Stability Problem of the Thomson Vortex Pentagon

Let N = 5 and the condition
Γ 15 ( q ) < Γ < Γ 05 ( q ) , 0 < q < 1
is valid. Here functions Γ 05 and Γ 15 are given by Formula (35). Then the points ( q , Γ ) lie in the white domain in the Figure 3, where linear stability takes place and the stability problem requires nonlinear analysis.
Following the paper [18], we introduce change of variables
ρ = A 5 ξ ζ , ξ = ( ξ 1 , , ξ 5 ) T , ζ = ( ζ 1 , , ζ 5 ) T ,
where A 5 is symplectic matrix
A 5 = [ A 51 ; A 52 ] , A 51 = 1 5 ν 1 h 4 ν 2 h 3 ν 2 h 3 ν 1 h 4 h 5 2 λ 15 h 1 ν 1 h 2 ν 2 h 2 ν 2 h 1 ν 1 h 0 , A 52 = 1 5 ν 1 h 1 ν 2 h 2 ν 2 h 2 ν 1 h 1 h 0 h 4 ν 1 h 3 ν 2 h 3 ν 2 h 4 ν 1 2 λ 15 h 5 .
Here h k R 5 , k = 1 , , 5 are basis vectors given by equalities (21) for N = 5 , h 0 is zero vector,
ν m = λ 2 m λ 1 m 4 , m = 1 , 2 .
In the considered condition (47) the statement 3 ° of the Theorem 1 is performed and the inequality λ 1 k ( q , Γ ) > 0 is valid for any k = 1 , , 5 .
In the new variables ( ξ , ζ ) , the Hamiltonian E ( ρ ( ξ , ζ ) ) defined by Formula (9) has cyclic variable ζ 5 . Assuming ξ 5 = 0 , we obtain reduced system.
We introduce the complex variables Z 1 , , Z 4 :
ξ k = 1 2 ( Z k + Z ¯ k ) , ζ k = 1 2 i ( Z k Z ¯ k ) , k = 1 , , 4 ,
and obtain reduced system with the Hamiltonian W. The calculations showed that the Taylor series expansion of this Hamiltonian in a neighborhood of the zero equilibrium has the form
W = 1 4 π ( W 0 + W 2 + Im ( a 5 Z 3 2 Z 4 ) + ) .
The terms of the third degree and higher except for the written one are omitted. The quadratic part W 2 of Hamiltonian W has a normal form:
W 2 = ω 1 | Z 1 | 2 + ω 2 | Z 2 | 2 ω 3 | Z 3 | 2 + ω 4 | Z 4 | 2 .
Values ω j ( j = 1 , , 4 ) are calculated with the help of the formulas
ω 3 = d e f 1 2 Im σ 3 + , ω k = d e f 1 2 Im σ k + , k = 1 , 2 , 4 ,
where σ k ± defined by equality (24).
According to the paper [18] in the case Γ = 0 the resonance ω 4 = 2 ω 3 leads to instability for q = q 1 : 2 = 0.333377 .
We study the resonance
ω 4 ( q , Γ ) = 2 ω 3 ( q , Γ )
in the case of the arbitrary Γ . The Equation (54) gives the curve Γ = Γ 1 : 2 ( q ) where the resonance 1:2 occurs:
Γ 1 : 2 ( q ) = B + B 2 A C A , A = 64 16 λ 24 λ 23 λ 24 2 16 λ 23 2 , B = 8 ( 2 λ 03 + λ 04 ) 2 ( λ 24 + 4 λ 23 ) + 256 λ 23 λ 24 ( λ 14 0 + λ 13 0 ) , C = 256 λ 14 0 λ 24 λ 13 0 λ 23 ( 2 λ 03 + λ 04 ) 2 4 λ 14 0 λ 24 16 λ 13 0 λ 23 2 .
Here the values λ 1 k 0 , λ 2 k and λ 0 k are given by Formulas (19) and (20). The curve Γ 1 : 2 ( q ) is defined on the interval q 1 : 2 q < 1 , q 1 : 2 = 0.2576828 . Its schematic image is shown in the Figure 6. The graphs of coefficient ω k ( q , Γ 1 : 2 ( q ) ) are shown in Figure 7 for q 1 : 2 q < 1 .
According to [23,32], in the absence of other resonances of the form
ω j = ω 3 , ω j = 2 ω 3 , 2 ω j = ω 3 , ω j + ω k = ω 3 , j , k = 1 , 2 , 4 ,
the zero equlibrium of the system with Hamiltonian (51)–(54) is unstable, if the inequality
a 5 ( q , Γ 1 : 2 ( q ) ) 0
is valid, that is shown in the Figure 8.
As seen in Figure 7, in addition to the equality (54) on the curve Γ 1 : 2 , other resonant relations of the form (56) take place only at two points. The equality ω 2 = ω 3 is performed at the point ( q 1 : 2 , Γ 1 : 2 ) , Γ 1 : 2 = Γ 1 : 2 ( q 1 : 2 ) = 0.2490584 , and the equality ω 2 = 2 ω 3 is valid at the point ( q 2 : 3 , Γ 2 : 3 ) . Here q 2 : 3 = 0.2956825 , Γ 2 : 3 = Γ 1 : 2 ( q 2 : 3 ) = 0.1410773 . The stability problem at these points requires additional nonlinear analysis.
Thus, the invariant set C 2 is unstable on the curve Γ 1 : 2 ( q ) when q 1 : 2 < q < 1 and q q 2 : 3 . The domain 0 < q < q 1 : 2 and Γ 15 < Γ < Γ 05 does not contain the points of resonance 1:2 ( ω 4 = 2 ω 3 ).

7. Conclusions

This paper is a natural continuation of the series of articles [13,16,17,18,19]. The novelty of the result is that the stability problem of the Thomson vortex N-gon outside circular domain is considered in the case of an arbitrary circulation Γ around the boundary. Previously in the papers [13,16,17,18,19] it was assumed that Γ = 0 .
The stability analysis was based on the Routh theory for the study of Hamiltonian systems with cyclic variables [30] and on a generalization of the approach developed in Section 3 of [27] for the study of a Thomson vortex N-gon in the case of an arbitrary Hamiltonian that depends only on the distances between point vortices.
The results of the study of the linearization matrix and the quadratic part of the Hamiltonian are presented as stability diagrams in Figure 3 if the conditions (33) and (34) are valid. For N 100 and in some cases N > 100 , we checked their execution numerically. For N = 3 a nonlinear stability analysis has been performed in the critical case of two-fold zero on the corresponding resonant curve Γ 03 ( q ) , q ( 0 , 1 ) . It is shown that the stationary rotation of the vortex triangle is unstable. For a vortex pentagon N = 5 , a nonlinear analysis has been performed in the critical case of the resonance 1:2 on the curve Γ 1 : 2 ( q ) , q ( q 1 : 2 , 1 ) \ { q 2 : 3 } , shown schematically in Figure 6. Instability has been proved. The resonance curves Γ 1 : 2 and Γ 15 intersect in the point ( q 1 : 2 , Γ 1 : 2 ) . There are no points of resonance 1:2 in domain 0 < q < q 1 : 2 and Γ 15 ( q ) < Γ < Γ 1 : 2 , although it requires nonlinear analysis.
Possible directions of development of this paper are related to the nonlinear analysis of the stability problem in the white domain in Figure 3 at Γ 1 N Γ Γ 0 N and on its boundaries. To do this, we need to apply the KAM-theory methods (see, for example [23]), in particular, list and study all resonances encountered in this problem. Among them we will consider the points found in the Section 6. These are the point ( q 1 : 2 , Γ 1 : 2 ) where two resonances ω 4 = 2 ω 3 and ω 2 = ω 3 intersect, and the point ( q 2 : 3 , Γ 2 : 3 ) , where the resonances are ω 4 = 2 ω 3 and ω 2 = 2 ω 3 .
The method of nonlinear analysis was previously tested for a model of vortices inside a circular domain in [14,15] and outside the circle in the case of zero circulation Γ = 0 in [16,17,18,19].
Of particular interest is the study of the boundary Γ 1 N of the instability domain. The answer to this question is important in order to find out what this border is: “dangerous” or “safe” according to Bautin [24,33]? In other words, there is a hard or soft loss of stability of the Thomson polygon when the parameter q increases and passes through the critical value q 1 N at a fixed Γ ?
The results of the study can also be generalized to the case of a moving cylinder interacting with the configuration of N point vortices [26], the model of Thomson N-gon in two-layer fluid [27,34] and others.

Author Contributions

Conceptualization, L.K. and I.O.; funding acquisition, L.K. and I.O; investigation, L.K. and I.O; methodology, L.K. and I.O; writing, original draft, L.K. and I.O; writing, review and editing, L.K. and I.O. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by Russian Foundation for Basic Research (Projects No. 20-55-10001).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Properties of Eigenvalues λ1k, λ2k and 0k of Matrices F1, F2 and G0

The following proposition sets the properties of quantities λ 1 k , λ 2 k and λ 0 k given by equalities (18)–(20).
Proposition A1.
The following statements are valid.
1 °
The function λ 2 k ( q ) is positive in the interval 0 < q < 1 :
λ 2 k ( q ) > 0 , 0 < q < 1 , k = 1 , , N 1 .
2 °
The function λ 0 k ( q ) does not have zeros in the interval 0 < q < 1 :
λ 0 k ( q ) 0 , 0 < q < 1 , k = 1 , , N 1 , k N 2 .
3 °
The function λ 1 N 0 ( q ) given by Formula (19) satisfies inequality
λ 1 N 0 ( q ) > λ 1 N 2 0 ( q )
in the interval 0 < q < 1 .
4 °
In the case Γ > Γ 1 N ( q ) the value λ 1 N is positive:
λ 1 N ( q , Γ ) > 0 ,
Γ 1 N is defined by Formula (35).
5 °
In the case of odd N = 2 + 1 for any q ( 0 , 1 ) the inequality
Γ 1 N ( q ) < Γ 0 N ( q )
is valid.
Proof. 
From the equalities (18)–(20), it follows that
λ m k = λ m , N k , m = 1 , 2 , λ 0 k = λ 0 , N k , k = 1 , , N 1 .
Thus, it is sufficient to consider the cases 1 k N 2 .
To prove the statement 1 ° , we write the value λ 2 k as
λ 2 k = Q 2 ( q ) 2 ( 1 q N ) 2 ,
where
Q 2 ( q ) = k ( N k ) q 2 N k N q 2 N k N ( N k ) q N + k + 2 ( N 2 + k 2 k N ) q N N ( N k ) q N k k N q k + k ( N k ) .
According to Descartes’ rule of signs [35], the number of positive roots of the polynomial Q 2 ( k , N ) is at most the number of sign changes in the sequence of its coefficients, and that the difference between these two numbers is always even. Zero coefficients are not taken into account when calculating sign changes.
The series of the polynomial Q 2 ( q ) coefficients has four sign changes, and the polynomial Q 2 ( q ) has a four-fold root q = 1 . Indeed, the polynomial Q 2 ( q ) and its first three derivatives Q 2 ( s ) ( q ) , s = 1 , 2 , 3 are
Q 2 ( q ) = 2 N k ( N k ) q 2 N 1 k N ( 2 N k ) q 2 N k 1 N ( N 2 k 2 ) q k + N 1 + + 2 N ( N 2 + k 2 k N ) q N 1 N ( N k ) 2 q N k 1 N k 2 q k 1
Q 2 ( q ) = 2 N k ( N k ) ( 2 N 1 ) q 2 N 2 N k ( 2 N k ) ( 2 N k 1 ) q 2 N k 2 N ( N 2 k 2 ) ( k + N 1 ) q N + k 2 + 2 N ( N 1 ) ( N 2 + k 2 k N ) q N 2 N ( N k ) 2 ( N k 1 ) q N k 2 N k 2 ( k 1 ) q k 2
Q 2 ( q ) = 2 N k ( N k ) ( 2 N 1 ) ( 2 N 2 ) q 2 N 3 N k ( 2 N k ) ( 2 N k 1 ) ( 2 N k 2 ) q 2 N k 3 N ( N 2 k 2 ) ( k + N 1 ) ( N + k 2 ) q N + k 3 + + 2 N ( N 1 ) ( N 2 ) ( N 2 + k 2 k N ) q N 3 N ( N k ) 2 ( N k 1 ) ( N k 2 ) q N k 3 N k 2 ( k 1 ) ( k 2 ) q k 3 .
They vanish in the point q = 1 . The four derivative at this point is different from zero:
Q 2 ( 4 ) ( 1 ) = 6 k 2 N 2 N k 2 0 .
Therefore, the polynomial Q 2 ( q ) does not have roots in the interval 0 < q < 1 . Since the value Q 2 ( 0 ) = k ( N k ) > 0 , then Q 2 ( q ) > 0 for all q ( 0 , 1 ) . The statement 1 ° is proved.
Similarly, the statement 2 ° is proved. The function λ 0 k is rewritten as
λ 0 k = Q 0 ( q ) ( 1 q N ) 2 ,
where
Q 0 ( q ) = N k q 2 N k N ( N k ) q k + N + N ( N k ) q N k N k q k .
The series of the polynomial Q 0 ( q ) coefficients has three sign changes, and the polynomial itself has a three-fold root q = 1 . According to Descartes’ rule of signs, there are no other positive roots of the polynomial Q 0 ( q ) and, accordingly, the function λ 0 k .
To prove the inequality (A3) in the statement 3 ° , we present the difference λ 1 N 0 λ 1 N 2 0 as:
λ 1 N 0 λ 1 N 2 0 = Q 1 ( q ) 2 ( 1 q N ) 2 , Q 1 ( q ) = k N k q 2 N + k N q 2 N k + N N k q N + k 2 k N + N 2 k 2 q N + + N N k q N k + k N q k + k N k , k = N 2 .
The polynomial Q 1 ( q ) has double root q = 1 and two sign changes in the series of coefficients. Therefore, there are no other positive roots. Since the value Q 1 ( 0 ) > 0 is positive, the inequalities Q 1 ( q ) > 0 and (A3) are valid for all q ( 0 , 1 ) and any N.
The property 4 ° follows from the statement 3 ° and the equality (35). In the case of Γ > Γ 1 N ( q ) we have:
λ 1 N ( q , Γ ) = λ 1 N 0 ( q ) + 2 Γ > λ 1 N 2 0 ( q ) + 2 Γ 1 N ( q ) 0 .
The validity of the property 5 ° follows from the equality (35) and the statements 1 ° and 2 ° . Indeed, for N = 2 + 1
Γ 0 N ( q ) = 1 2 λ 1 0 ( q ) + λ 0 ( q ) 2 λ 2 ( q ) > 1 2 λ 1 0 ( q ) = Γ 1 N ( q ) .
The Proposition A1 is proved. □

References

  1. Aref, H. Point vortex dynamics: A classical mathematics playground. J. Math. Phys. 2007, 48, 065401. [Google Scholar] [CrossRef] [Green Version]
  2. Aref, H.; Newton, P.K.; Stremler, M.A.; Tokieda, T.; Vainchtein, D.L. Vortex Crystals. Adv. Appl. Mech. 2003, 39, 1–79. [Google Scholar]
  3. Newton, P.K. The n-Vortex Problem: Analytical Techniques; Appl. Math. Sci.; Springer: New York, NY, USA, 2001; Volume 145. [Google Scholar]
  4. Borisov, A.V.; Mamaev, I.S. Mathematical Methods in the Dynamics of Vortex Structures; Institute of Computer Sciences: Moscow/Izhevsk, Russia, 2005. [Google Scholar]
  5. Kilin, A.A.; Borisov, A.V.; Mamaev, I.S. Dynamics of point vortices inside and outside of a circular domain. In Fundamental and Applied Problems in Vortex Theory; Institute of Computer Sciences: Moscow/Izhevsk, Russia, 2003; pp. 414–440. [Google Scholar]
  6. Saffman, P.G. Vortex Dynamics; Series Cambridge Monogr. Mech. and Appl. Math. Cambridge University Press: Cambridge, UK, 1992. [Google Scholar]
  7. Kozlov, V.V. General Theory of Vortices; Dynamical systems X., Encyclopaedia Math. Sci., 67; Springer: Berlin, Germany, 2003. [Google Scholar]
  8. Yarmchuk, E.; Gordon, M.; Packard, R. Observation of stationary vortex array in rotating superfluid helium. Phys. Rev. Lett. 1979, 43, 214–217. [Google Scholar] [CrossRef]
  9. Fine, K.; Cass, A.; Flynn, W.; Dryscoll, C. Relaxation of 2D Turbulence to Vortex Crystal. Phys. Rev. Lett. 1995, 75, 3277–3280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Kurakin, L.G.; Yudovich, V.I. The stability of stationary rotation of a regular vortex polygon. Chaos 2002, 12, 574–595. [Google Scholar] [CrossRef] [PubMed]
  11. Lin, C.C. On the motion of vortices in two dimensions. I, II. Proc. Natl. Acad. Sci. USA 1941, 27, 575–577. [Google Scholar] [CrossRef] [Green Version]
  12. Milne-Thomson, L.M. Theoretical Hydrodynamics; Macmillan: London, UK, 1968. [Google Scholar]
  13. Havelock, T.H. The stability of motion of rectilinear vortices in ring formation. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1931, 11, 617–633. [Google Scholar] [CrossRef]
  14. Kurakin, L.G. Stability, resonances, and instability of the regular vortex polygons in the circular domain. Dokl. Phys. 2004, 49, 658–661. [Google Scholar] [CrossRef]
  15. Kurakin, L.G. On the Stability of Thomson’s Vortex Pentagon Inside a Circular Domain. Regul. Chaotic Dyn. 2012, 17, 150–169. [Google Scholar] [CrossRef]
  16. Kurakin, L.G.; Ostrovskaya, I.V. Stability of the Thomson Vortex Polygon with Evenly Many Vortices outside a Circular Domain. Sib. Math. J. 2010, 51, 463–474. [Google Scholar] [CrossRef]
  17. Kurakin, L.G. The stability of the steady rotation of a system of three equidistant vortices outside a circle. J. Appl. Math. Mech. 2011, 75, 227–234. [Google Scholar] [CrossRef]
  18. Kurakin, L.G.; Ostrovskaya, I.V. Nonlinear stability analysis of a regular vortex pentagon outside a circle. Regul. Chaotic Dyn. 2012, 17, 385–396. [Google Scholar] [CrossRef]
  19. Kurakin, L.G.; Melekhov, A.P.; Ostrovskaya, I.V. A survey of the stability criteria of Thomson’s vortex polygons outside a circular domain. Bol. Soc. Mat. Mex. 2016, 22, 733–744. [Google Scholar] [CrossRef]
  20. Kolmogorov, A.N. On conservation of conditionally periodic motions for a small change in hamilton’s function. Dokl. Akad. Nauk SSSR 1954, 98, 527–530. [Google Scholar]
  21. Arnold, V.I. Small denominators and problems of stability of motion in classical and celestial mechanics. Russ. Math. Surv. 1963, 18, 85–191. [Google Scholar] [CrossRef]
  22. Moser, J. Lectures on Hamiltonian Systems; Memoirs of the American Mathematical Society: Providence, RI, USA, 1968; Volume 81. [Google Scholar]
  23. Markeev, A.P. Libration Points in Celestial Mechanics and Space Dynamics; Nauka: Moscow, Russia, 1978. [Google Scholar]
  24. Khazin, L.G.; Shnol, E.E. Stability of Critical Equilibrium States; Manchester University Press: Manchester, NH, USA, 1991. [Google Scholar]
  25. Koshel, K.V.; Ryzhov, E.A.; Carton, X.J. Vortex Interactions Subjected to Deformation Flows: A Review. Fluids 2019, 4, 14. [Google Scholar] [CrossRef] [Green Version]
  26. Borisov, A.V.; Mamaev, I.S.; Romodanov, S.M. Dynamic interaction of point vortices and a two-dimentional cylinder. J. Math. Phys. 2007, 48, 065403. [Google Scholar] [CrossRef] [Green Version]
  27. Kurakin, L.G.; Lysenko, I.A.; Ostrovskaya, I.V.; Sokolovskiy, M.A. On stability of the Thomson’s vortex n-gon in the geostrophic model of the point vortices in two-layer fluid. J. Nonlinear Sci. 2019, 29, 1659–1700. [Google Scholar] [CrossRef]
  28. Kurakin, L.G.; Lysenko, I.A. On the stability of the orbit and the invariant set of Thomson’s vortex polygon in a two-fluid plasma. Rus. J. Nonlin. Dyn. 2020, 16, 3–11. [Google Scholar]
  29. Gantmacher, F.R. Lectures on Analytical Mechanics; Glavnaya Redaktsiya Fiziko-Matematicheskoj Literatury, Nauka: Moscow, Russia, 1966. [Google Scholar]
  30. Routh, E.J. A Treatise on the Stability of a Given State Motion; Macmillan: London, UK, 1877; 108p. [Google Scholar]
  31. Sokol’sky, A.G. On Stability of an Autonomous Hamiltonian System with Two Degrees of Freedom under First-Order Resonance. J. Appl. Math. Mech. 1977, 41, 20–28. [Google Scholar] [CrossRef]
  32. Markeev, A.P. Stability of a canonical system with two degrees of freedom in the presence of resonance. J. Appl. Math. Mech. 1968, 32, 766–772. [Google Scholar] [CrossRef]
  33. Bautin, N.N. Behavior of Dynamics Systems Near the Boundary of the Stability Region; Gostekhizdat: Moscow, Russia, 1984. [Google Scholar]
  34. Kurakin, L.G.; Ostrovskaya, I.V.; Sokolovskiy, M.A. On the stability of discrete tripole, quadrupole, Thomson’ vortex triangle and square in a two-layer/homogeneous rotating fluid. Regul. Chaotic Dyn. 2016, 21, 291–334. [Google Scholar] [CrossRef]
  35. Kurosh, A.G. Course of Higher Algebra; Nauka: Moscow, Russia, 1962. [Google Scholar]
Figure 1. Thomson vortex N-gon outside a circle.
Figure 1. Thomson vortex N-gon outside a circle.
Mathematics 08 01033 g001
Figure 2. The graphs of the curves Γ N ( q ) in the interval 0 < q < 1 such that the angular velocity Ω N ( q , Γ N ( q ) ) 0 . The graphs are shown from top to bottom in order of increasing N = 2 , , 10 .
Figure 2. The graphs of the curves Γ N ( q ) in the interval 0 < q < 1 such that the angular velocity Ω N ( q , Γ N ( q ) ) 0 . The graphs are shown from top to bottom in order of increasing N = 2 , , 10 .
Mathematics 08 01033 g002
Figure 3. Schematic stability diagram of the Thomson vortex N-gon: (a) The case of odd N = 2 + 1 ; (b) The case of even N = 2 n . Gray is the area of orbital stability in an exact nonlinear setting, the vertical shading is the area of instability, white color and neutral curves Γ 0 N , Γ 1 N are the area of linear stability, where is required nonlinear analysis. The diagram is true for N 100 and for N > 100 if the conditions (33) and (34) are valid.
Figure 3. Schematic stability diagram of the Thomson vortex N-gon: (a) The case of odd N = 2 + 1 ; (b) The case of even N = 2 n . Gray is the area of orbital stability in an exact nonlinear setting, the vertical shading is the area of instability, white color and neutral curves Γ 0 N , Γ 1 N are the area of linear stability, where is required nonlinear analysis. The diagram is true for N 100 and for N > 100 if the conditions (33) and (34) are valid.
Mathematics 08 01033 g003
Figure 4. The graph of the coefficient ω 1 ( q ) given by equality (45) in the interval 0 < q < 1 .
Figure 4. The graph of the coefficient ω 1 ( q ) given by equality (45) in the interval 0 < q < 1 .
Mathematics 08 01033 g004
Figure 5. The graph of the coefficient a 3 ( q ) given by Formula (46) in the interval 0 < q < 1 .
Figure 5. The graph of the coefficient a 3 ( q ) given by Formula (46) in the interval 0 < q < 1 .
Mathematics 08 01033 g005
Figure 6. Schematic stability diagram of the Thomson vortex pentagon. Grey color is the domain of orbital stability in an exact nonlinear setting, the vertical shading and the curve Γ 1 : 2 are instability domain, white color and neutral curves Γ 05 , Γ 15 are linear stability domains, which need nonlinear analysis.
Figure 6. Schematic stability diagram of the Thomson vortex pentagon. Grey color is the domain of orbital stability in an exact nonlinear setting, the vertical shading and the curve Γ 1 : 2 are instability domain, white color and neutral curves Γ 05 , Γ 15 are linear stability domains, which need nonlinear analysis.
Mathematics 08 01033 g006
Figure 7. The graphs of coefficients ω j ( q , Γ ) , j = 1 , , 4 on the resonance curve Γ = Γ 1 : 2 ( q ) for q 1 : 2 < q < 1 , ω 4 = 2 ω 3 .
Figure 7. The graphs of coefficients ω j ( q , Γ ) , j = 1 , , 4 on the resonance curve Γ = Γ 1 : 2 ( q ) for q 1 : 2 < q < 1 , ω 4 = 2 ω 3 .
Mathematics 08 01033 g007
Figure 8. The resonance 1:2. The graph of coefficient a 5 ( q , Γ 1 : 2 ( q ) ) in the interval q 1 : 2 < q < 1 .
Figure 8. The resonance 1:2. The graph of coefficient a 5 ( q , Γ 1 : 2 ( q ) ) in the interval q 1 : 2 < q < 1 .
Mathematics 08 01033 g008
Table 1. The critical values q 0 N and q 1 N for N = 2 , , 6 in the case Γ = 0 .
Table 1. The critical values q 0 N and q 1 N for N = 2 , , 6 in the case Γ = 0 .
N q 0 N q 1 N q 1 N q 0 N
2 q 02 = q 12 q 12 = 0.148536 0
3 q 03 = 0.262542 q 13 = 0.273695 0.011152
4 q 04 = q 14 q 14 = 0.308125 0
5 q 05 = 0.330399 q 15 = 0.334596 0.004197
6 q 06 = q 16 q 16 = 0.295985 0
Table 2. The critical values of Γ N for N = 2 , , 11 .
Table 2. The critical values of Γ N for N = 2 , , 11 .
N234567891011
Γ N 0.25 0.5 0.5 0.5 0.25 0 0.5 1 1.75 2.5

Share and Cite

MDPI and ACS Style

Kurakin, L.; Ostrovskaya, I. On the Effects of Circulation around a Circle on the Stability of a Thomson Vortex N-gon. Mathematics 2020, 8, 1033. https://doi.org/10.3390/math8061033

AMA Style

Kurakin L, Ostrovskaya I. On the Effects of Circulation around a Circle on the Stability of a Thomson Vortex N-gon. Mathematics. 2020; 8(6):1033. https://doi.org/10.3390/math8061033

Chicago/Turabian Style

Kurakin, Leonid, and Irina Ostrovskaya. 2020. "On the Effects of Circulation around a Circle on the Stability of a Thomson Vortex N-gon" Mathematics 8, no. 6: 1033. https://doi.org/10.3390/math8061033

APA Style

Kurakin, L., & Ostrovskaya, I. (2020). On the Effects of Circulation around a Circle on the Stability of a Thomson Vortex N-gon. Mathematics, 8(6), 1033. https://doi.org/10.3390/math8061033

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop