Russian Journal of Nonlinear Dynamics, 2019, vol. 15, no. 4, pp. 533-542. Full-texts are available at http://nd.ics.org.ru DOI: 10.20537/nd190412
MSC 2010: 76B47, 76E20, 93B18
On the Stability of Thomson's Vortex N-gon and a Vortex Tripole/Quadrupole in Geostrophic Models of Bessel Vortices and in a Two-Layer Rotating Fluid: a Review
L. G. Kurakin, I. V. Ostrovskaya
In this paper the two-layer geostrophic model of the rotating fluid and the model of Bessel vortices are considered. Kirchhoff's model of vortices in a homogeneous fluid is the limiting case of both of these models. Part of the study is performed for an arbitrary Hamiltonian depending on the distances between point vortices.
The review of the stability problem of stationary rotation of regular Thomson's vortex N-gon of identical vortices is given for N ^ 2. The stability problem of the vortex tripole/quad-rupole is also considered. This axisymmetric vortex structure consists of a central vortex of an arbitrary intensity and two/three identical peripheral vortices. In the model of a two-layer fluid, peripheral vortices belong to one of the layers and the central vortex can belong to either another layer or the same.
The stability of the stationary rotation is interpreted as orbital stability (the stability of a one-parameter orbit of a stationary rotation of a vortex system). The instability of the stationary rotation is instability of equilibrium of the reduced system. The quadratic part of the Hamiltonian and eigenvalues of the linearization matrix are studied.
Received June 20, 2019 Accepted November 05, 2019
L.G. Kurakin was supported by the Russian Foundation for Basic Research (Projects No. 20-55-10001) and I. V. Ostrovskaya was supported by Ministry of Education and Science of the Russian Federation, Southern Federal University (Projects No. 1.5169.2017/8.9)
Leonid G. Kurakin kurakin@math.rsu.ru
Southern Mathematical Institute, Vladikavkaz Scientific Center of RAS
ul. Markusa 22, Vladikavkaz, 362027 Russia
Water Problems Institute of RAS
ul. Gubkina 3, Moscow, 119333 Russia
Irina V. Ostrovskaya ivostrovskaya@sfedu.ru
Southern Federal University
ul. Milchakova 8a, Rostov-on-Don, 344090 Russia
The parameter space is divided into three parts: A is the domain of stability in an exact nonlinear setting, B is the linear stability domain, where the stability problem requires nonlinear analysis, and C is the instability domain.
In the stability problem of a vortex multipole, another definition of stability is used; it is the stability of an invariant three-parametric set of all trajectories of the families of stationary orbits. It is shown that in the case of non zero total intensity, the stability of the invariant set implies orbital stability.
Keywords: N-vortex problem, Thomson's vortex N-gon, point vortices, two-layer fluid, stability, Hamiltonian equation
1. The stability of the Thomson vortex N-gon
for a system of particles with potential of interaction depending
on the distances between them
We consider the motion of the system of N equal point vortices with intensity r = 0 described by the Hamiltonian depending on the distances between vortices
r
Kj<k^N
where zk = qk + ipk, (qk ,pk) are Cartesian coordinates of the kth vortex and V(£) is a smooth function defined on R+.
The form of function V(£) defines the vortex dynamics model under study. The formula
V (0 = ln £ (1.2)
defines Kirchhoff's model describing the motion of point vortices in a homogeneous fluid. In the
case
V (£) = -Ko(7£) (1.3)
we have a model of the geostrophic Bessel vortices [11, 14, 19, 20], and if
h
the geostrophic vortices are located in the two-layer rotating fluid [4, 18]. Here K0(£) is a modified Bessel function of zero order, 1/y > 0 is the internal Rossby deformation radius, and hi, h2 are the thickness of the top and bottom layer, respectively, hi+h2 = 1. The case 7 = 1 can be reduced to the case 7 = 1 by scaling. Note that function (1.2) is the limit for functions (1.3) and (1.4) as 7 ^ 0 and hi ^ 0, respectively.
We consider the stability problem of the stationary rotation of a system of N identical point vortices lying uniformly on a circle of radius R. In the two-layer model the vortices are located in a single layer of a two-layer fluid (for example, in the bottom layer). The system with Hamiltonian (1.1) has the exact solution
zk = eiUNtuk, uk = Re2ni(k-i)/N, R> 0, k = 1,...,N, (1.5)
corresponding to the stationary rotation of the configuration of N point vortices lying on a circle of radius R at the vertices of the regular N-gon. The angular velocity wN is given by the expression
„ N — 1
= E F (2Esm iv")sm iv"- d-6)
m=1
The stability of the stationary motion is interpreted as orbital stability.
The stability problem of the stationary rotation of the system of N identical point vortices located uniformly on a circle (Thomson's vortex N-gon) in a homogeneous fluid (1.2) was posed by Kelvin (W.Thomson) [15]. J.J.Thomson [16] and T. H. Havelock [5] investigated this problem completely in the linear setting. A long history of research on this problem in the exact nonlinear setting for Kirchhoff's model is described in [1, 8], where a review of experimental and theoretical works is given. The change of variables
zk(t) = e^vUt), vk = k = l,2,...,N (1.7)
leads the Hamiltonian (1.1) to perturbation equations with the Hamiltonian
N
E(p) = H(p) + fM(p), M = J2(R2 + 2 rk), (1.8)
k=1
P = (ri, ...,TN ,01,.. .,0N).
The equilibrium family C = {r1 = ... = rN = O,01 = ... = 0N} corresponds to stationary rotation (1.5). The stability of this family is the orbital stability of stationary rotation. The instability of the stationary rotation of an N-gon (1.5) is the instability of an invariant two-dimensional set C2 of all orbits of a continuous family of stationary rotations (1.5) depending on R:
C2 = {P e R2N : r1 = ... = rN, 01 = ... = 0N},
or Lyapunov instability of zero equilibrium of the reduced system.
The quadratic part of the Taylor series expansion of the Hamiltonian E in powers of p in the neighborhood of the zero solution has the form
( F1 Go\
E2(p) = (Sp,p), S = 1 G . (1.9)
The linearization matrix is
—Go F2
t ~ Go F2 \
L = 2 g0 * • (1.10)
y— Fi —Goy The matrices F1, F2 and Go are circulants
N-1 N-1
Fj = ^ fjmCm, G0 = ^ gomCm.
m=0 m=1
The analytical formulas of coefficients fjm and g0m for an arbitrary Hamiltonian are given by formulas (3.15), (3.16) of [12]. The cyclic matrix C = {cjk}NNk=1 has nontrivial coefficients, namely, only those above its principal diagonal c1>2 = ... = cN-1>N = 1 and in the left lower corner cN 1 = 1.
Analysis of eigenvalues of matrices S and L gives the following statements [12].
A The Thomson N-gon (1.5) is orbitally stable in an exact nonlinear setting if all eigenvalues of matrix S have identical sign, except a simple zero.
C The Thomson N-gon (1.5) is instable if the linearization matrix L has at least one eigenvalue in the right half-plane.
B The stability problem requires nonlinear analysis if conditions A and C are not valid.
In the case of the Kirchhoff model the rotation of the vortex N-gon is stable only if N ^ 7, whereas for N ^ 8 the motion is unstable (C) [5, 8]. In the case N ^ 6 a linear analysis is enough to conclude nonlinear stability (A). For N = 7, it is necessary to involve nonlinear terms for analysis of the stability [8].
1.1. Thomson's N-gon in the geostrophic model of the Bessel vortices
The motion of the system of N point geostrophic (Bessel) vortices is described by the Hamiltonian (1.1), (1.3). The parameter space is (N,R). Stewart [19, 20] performed a linear analysis of the stability of the regular vortex N-gon and came to the following conclusions: 1) Stability always occurs for N ^ 6 (there are no exponentially growing modes); 2) The vortex heptagon (N = 7) is stable if R > 71, and unstable for R < 71; 3) Instability always occurs for N ^ 8.
Morikawa and Swenson [14] specified that the vortex hexagon (N = 6) is linearly stable at R < Ro6 and exponentially unstable if R > R06, Ro6 ~ 1.288.
AT = 2,3,4
" N = 21 + 1^7
0 R '' C B
A B 0 —N = 5 R*{N) R
0 Ros R C
A C —»► N = 6 0 R
0 Roe R
N = 2n>8
Fig. 1. Stability diagrams of the vortex N-gon in the model of Bessel vortices: A — the stability in the case of exact nonlinear problem setting, B — the problem requires nonlinear analysis, and C — instability.
In [11] it is proved that the vortex N-gon (N = 2,3,4) is orbitally stable (A) in the exact nonlinear setting for R > 0 (see Fig. 1). The vortex pentagon (N = 5) is orbitally stable (A) if R < R05, R05 ~ 3.76 and the stability problem requires nonlinear analysis (B) if R ^ R05. The vortex hexagon (N = 6) is orbitally stable (A) if R < R06. In the case of odd N = 21+1 ^ 7 the vortex N-gon is unstable (C) if R < R*(N) and the stability problem requires nonlinear analysis (B) if R ^ R*(N). In the case of even N = 2m ^ 8 the vortex N-gon is unstable (C) for any R > 0. The asymptotic formulas for the value R*(N) are
R*(N) « Ra(N) =
2 sin
3 _7L
N
2cOS°f
1 + cos
n(N - 1) N~
1
+
4sin^
(1.11)
1
1
1.2. Thomson's N-gon in a two-layer fluid
The motion of N -point vortices in a two-layer fluid is described by the Hamiltonian (1.1), (1.4). The parameter space is (N,R,a). The parameter a = h2 — h\ is the difference between the thicknesses of the lower and upper layers, —1 < a < 1.
The case A takes place for N = 2, 3,4 for all R> 0 and —1 < a < 1 [10].
In the case of N = 5 we have two domains: A and B (see Fig. 2a). At N = 6 cases A and C take place (see Fig. 2b). The case B is the neutral curve a6 defined for R06 ^ R < ж>. In the case of N = 7 there are all three domains A, B and C (see Fig. 3). The instability C takes place always for even N = 2m ^ 8. In the case of odd N = 2i + 1 > 9 there are two domains B and C (see Fig. 4). The analytical formulas of neutral curves an have been obtained in [10].
0 2 4 6 „ 10 12 14 16 K
(b)
Fig. 2. Stability diagrams of a vortex N-gon: a) vortex pentagon; b) vortex hexagon.
Fig. 3. Stability diagram of a vortex heptagon: a) 0 < R < 100, b) \R — R*(7)| < 10"21, R*(7) « 71.52720252.
в 0
R*(N) oo
Fig. 4. The schematic diagram of stability for an N-gon for odd N = 2^+1 > 9. R* (N) is given by (1.11).
2. The stability problem for a vortex tripole and a vortex quadrupole
2.1. The model of a homogenous fluid (Kirchhoff's equation)
The motion of the system of one vortex with arbitrary intensity r = 0 and N vortices of unit intensity is described by the Hamiltonian
n = ~ , 4n
1 N
rXW(?o -fl02 + (^-pfc) + (2-1)
k=i
+ Yj ln \ '//■ )"' • il'j
The system has the first integrals: the Hamiltonian (2.1)
NN
Ii = rqo + Y Qk, I2 = Y Pk, k=l k=0 2 N
k=l
In complex variables zk = qk + ipk, the problem has the exact solution
zk{t) = eluJNtuk, «o = 0, uk = Rei^k~1\ k = l,2,...,N, (2.2)
N - 1 + 2r , x
= -Siff— (2-3)
corresponding to stationary rotation of configuration of a central vortex of intensity r and N vortices of equal intensity lying uniformly on a circle of radius R.
The orbital stability is the stability of a one-parameter orbit of a stationary rotation of a vortex multipole (2.2).
The vortex tripole (N = 2) is linearly stable if T < and unstable (C) if T > -| [13].
The vortex quadrupole (N = 3) is linearly stable if r < 1 and unstable (C) if r > 1 [3]. The vortex quadrupole (N = 3) is orbitally stable if 0 < r < 1, the problem requires nonlinear analysis if r < 0 [9, 10]. To our surprise, we have found no earlier references to nonlinear stability analysis. The proof and the nonlinear stability condition given in [2] are incorrect.
Let Mn be the invariant three-parametric set of all trajectories of the families of stationary orbits
zk(t) = ukeiMN 1 + n, k = 0,...,N, n e C. (2.4)
Z.Kizner [6, 7] used the following definition of stability: the vortex multipole is steady if, at a reasonably small initial perturbation of distances between its vortices, these distances remain weakly perturbed for all time. Such stability is stability of the invariant set MN.
The invariant set Mn is nonlinearly stable if the first integral of the system exists which has a local extremum on the set MN.
To study the stability of the invariant set MN, we introduce new variables Xl, X2, 62 for N = 2 or Xi, X2, X3, 82, 83 for N = 3 (see Fig. 5).
Fig. 5
If N = 2, the system integrals H and L in the new variables are rewritten as
U{XllX2l 62) = r(lnX\ + lnX2) + In yjxf + xf- 2XiX2cose2, (2.5)
l(Xi ,X2, e2) = r(X2 + X22 ) + X2 + X22 - 2X1 X2 cos e2. (2.6)
The point (R, R, n) in the new variables corresponds to the invariant set M2. The invariant set M2 is stable if the integral
i (Xi, X2, e2) = (H(Xi ,X2, e2) - H*)2 + (L(Xi ,X2, e2) - L* )2
has a strict local minimum at the point (R, R, n), H* = H(R, R, n), L* = L(R, R, n).
The invariant set WI2 is stable if F < — | (see [6, 7]), and the invariant set WI3 is stable if r < -3 and 0 < r < 1 [10].
2.2. The model of a two-layer rotating fluid
A set of N +1 point vortices in a two-layer fluid is considered for N = 2, 3: one vortex of arbitrary intensity k0 = 0 with coordinates (q0 ,p0) and N vortices with coordinates (qk,pk), k = 1,...,N with equal intensity k lying in one layer, for definiteness, in the lower layer (kh2 = 1). Two cases are distinguished: (I) the vortex of intensity k0 is located in the upper layer with thickness h1, r = k0h1; (II) all N + 1 vortices are located in the lower layer with thickness h2, r = k0h2. The motion of the vortex system is described by the Hamilto-nian H [4, 18]:
(2.7)
L k=i v v ' + E WVfe-^ + fe-»)2)
1- a i ln £ + K0 (£) incase (I),
1 + a W(£) in case (II).
A complete analysis of the stability problem of the stationary rotation (2.2) of a vortex tripole and a vortex quadrupole with constant angular velocity
= -¿x{rV'(.R) + W(2R)), = ^(rV\R) + V3W'{V3R))
is presented in [6, 7, 9, 10, 17]. Its results are presented in Figs. 6-9. The formulas of curves Tj, j = 1,..., 6 are obtained in [9, 10].
0.9 1.0 1.1 1.2
Fig. 6. Stability diagram of a vortex tripole. Case I. a = 0. A — dark domain, B — light domain, C dashed domain.
Fig. 7. Stability diagram of a vortex tripole. Case II: a) a = 0.8; b) a = 0; c) a = —0.8.
Fig. 8. Stability diagram of a vortex quadrupole. Case I: a-b) a = 0.8. The stability of the invariant set — grey domain.
2. 0.
-2. -4.
0.0 2.0 4.0 6.0
T = -3
M
2.0
<9-° -2.0
-4.0
0.0 2.0 4.0 6.0
T = -3
M
2.0 0.0 -2.0 -4.0
r
0.0 2.0 4.0 6.0 8.
r = -3
(c)
Fig. 9. Stability diagram of a vortex quadrupole. Case II: a) a = 0.8; b) a = 0; c) a = —0.8.
The invariant set M2 is stable everywhere in the orbital stability domain A and in domain B, where nonlinear analysis is required [6, 7].
The stability of the invariant set in the case of a two-layer fluid is studied in [10]. In the stability diagrams presented in Figs. 8, 9, the domain of stability of the invariant set consists of a dark domain of orbital stability A and a grey domain, where the invariant set is stable and condition B is valid (see Section 1).
2.3. On equivalence of stability of invariant set and orbital stability in case of non zero total intensity
Note that the integral K = If + on set MN has the value
K = (r + N )2|n|2. (2.9)
The integral
J = (H — H*)2 + (L — L*)2 + K, (2.10)
where H* = HIMn, L* = LIMn, has a strict local minimum at n = 0 on the orbit (2.2), (2.3) if r + N = 0.
Remark. In all the above cases the stability of the invariant set MN implies orbital stability of the stationary rotation of the vortex multipole.
References
[1] Borisov, A. V. and Mamaev, I.S., Mathematical Methods in the Dynamics of Vortex Structures, Izhevsk: R&C Dynamics, Institute of Computer Science, 2005 (Russian).
[2] Cabral, H. E. and Schmidt, D. S., Stability of Relative Equilibria in the Problem of N +1 Vortices, SIAM J. Math. Anal, 1999/2000, vol.31, no. 2, pp. 231-250.
[3] Campbell, L. J., Transverse Normal Modes of Finite Vortex Arrays, Phys. Rev. A, 1981, vol.24, no. 1, pp. 514-534.
[4] Gryanik, V. M., Dynamics of Singular Geostrophical Vortices in a 2-Level Model of the Atmosphere (Ocean), Izv. Atmos. Ocean Phys., 1983, vol. 19, no. 3, pp. 171-179; see also: Izv. Akad. Nauk SSSR. Fiz. Atmos. Okeana, 1983, vol. 19, no. 3, pp. 227-240.
[5] Havelock, T.H., The Stability of Motion of Rectilinear Vortices in Ring Formation, Philos. Mag., 1931, vol. 11, no. 70, pp. 617-633.
[6] Kizner, Z., Stability of Point-Vortex Multipoles Revisited, Phys. Fluids, 2001, vol.23, no. 6, 064104, 11pp.
[7] Kizner, Z., On the Stability of Two-Layer Geostrophic Point-Vortex Multipoles, Phys. Fluids, 2014, vol.26, no.4, 046602, 18pp.
[8] Kurakin, L. G. and Yudovich, V. I., The Stability of Stationary Rotation of a Regular Vortex Polygon, Chaos, 2002, vol. 12, no. 3, pp. 574-595.
[9] Kurakin, L.G., Ostrovskaya, I. V., and Sokolovskiy, M. A., Stability of Discrete Vortex Multipoles in Homogeneous and Two-Layer Rotating Fluid, Dokl. Phys., 2015, vol. 60, no. 5, pp. 217-223; see also: Dokl. Akad.. Nauk, 2015, vol.462, no. 2, pp. 161-167.
[10] Kurakin, L.G., Ostrovskaya, I.V., and Sokolovskiy, M. A., On the Stability of Discrete Tripole, Quadrupole, Thomson' Vortex Triangle and Square in a Two-Layer/Homogeneous Rotating Fluid, Regul. Chaotnc Dyn, 2016, vol.21, no. 3, pp. 291-334.
[11] Kurakin, L.G. and Ostrovskaya, I.V., On Stability of the Thomson's Vortex N-Gon in the Geostrophic Model of the Point Bessel Vortices, Regul. Chaotic Dyn., 2017, vol. 22, no. 7, pp. 865-879.
[12] Kurakin, L. G., Lysenko, I. A., Ostrovskaya, I. V., and 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 Set., 2019, vol.29, no.4, pp. 1659-1700.
[13] Mertz, G., Stability of Body-Centered Polygonal Configurations of Ideal Vortices, Phys. Fluids, 1978, vol.21, no. 7, pp. 1092-1095.
[14] Morikawa, G.K. and Swenson, E. V., Interacting Motion of Rectilinear Geostrophic Vortices, Phys. Fluids, 1971, vol. 14, no. 6, pp. 1058-1073.
[15] Thomson, W., Floating Magnets (Illustrating Vortex Systems), Nature, 1878, vol. 18, pp. 13-14. See also: Kelvin, W. T., Mathematical and Physical Papers: Vol. 4, Cambridge: Cambridge Univ. Press, 1910, pp. 162-164.
[16] Thomson, J. J., Treatise on the Motion of Vortex Rings, London: Macmillan, 1883, pp. 94-108.
[17] Sokolovskiy, M. A. and Verron, J., Some Properties of Motion of A +1 Vortices in a Two-Layer Rotating Fluid, Nelin. Dinam., 2006, vol.2, no. 1, pp. 27-54 (Russian).
[18] Sokolovskiy, M.A. and Verron, J., Dynamics of Vortex Structures in a Stratified Rotating Fluid, Atmos. Oceanogr. Sci. Libr., vol.47, Cham: Springer, 2014.
[19] Stewart, H. J., Periodic Properties of the Semi-Permanent Atmospheric Pressure Systems, Quart. Appl. Math, 1943, vol.1, pp. 262-267.
[20] Stewart, H. J., Hydrodynamic Problems Arising from the Investigation of the Transverse Circulation in the Atmosphere, Bull. Amer. Math. Soe., 1945, vol.51, pp. 781-799.