We consider flows on a spherical surface and use a transformation to transport some well-known periodic two-dimensional vortex streets to that spherical surface to arrive at some new expressions for vortex streets on a sphere.
1. Introduction
For fluid flow on a two-dimensional plane, the vorticity at a point is twice the angular rotation. A point vortex is a model of a flow in which the vorticity is zero except at the point itself where the vorticity is infinite, so that there is a nonzero circulation around the point. The study of point vortices on the plane, and other two-dimensional manifolds such as the cylinder, sphere, and torus, has a long history, dating back to the 19th century with Helmholtz [1] initiating the point vortex model and Kirchhoff [2] and Lin [3] formulating it as a Hamiltonian dynamical system. In this paper, we are concerned primarily with vortex streets, which consist of one or more periodic rows of point vortices, the simplest of which is a single infinite row of identical vortices [4]. These have important applications in engineering and geophysics, with a single row having been used to model the quasisteady large-scale vortices arising following the roll-up of a shear layer, and double rows, or von Kármán vortex streets [5], having been used to model the shedding of eddies behind a bluff body. An overview of two-dimensional vortex streets can be found in standard texts on hydrodynamics such as IN [6, 7].
In this paper, we are interested in transporting these well-known vortex streets from the plane to a curved two-dimensional manifold, the surface of a sphere. Flows on a sphere are important because of applications to planetary atmospheres. In his classic monograph, Lamb [6] briefly outlines a method of determining the motion of vortices on a curved manifold and discusses how some of the 19th century work on electrical conduction, such as charge-on-a-sphere problems, by Boltzmann, Kirchhoff, Töpler, and others could be applied to the problem of point vortices on the sphere although Gromeka [8] appears to have been the first to study vortices on a sphere specifically. More recently, the formulation of the motion of vortex streets on curved manifolds has been examined in more detail by Hally [9], with several subsequent studies [10–12] delving more deeply into the formulation of vortex motion on a sphere. A review of some of the work on point vortices on vortices on a sphere can be found in [13]. One interesting thread of research [14–16] has involved using numerical methods, such as contour surgery, to study the motion of vortex patches on a sphere, and although these results were numerical rather than the closed-form expressions sought in the current work, they serve to reveal the richness of vortex motion on a sphere.
The outline of the rest of the paper is as follows. In Section 2, we present our analysis, giving a brief overview of two-dimensional vortex streets in Section 2.1 and then transporting these streets to the sphere in Section 2.2, where we look at two vortex streets in detail, a row of corotating vortices and a row of counterrotating vortices. Finally, in Section 3, we make some closing remarks.
2. Analysis2.1. Two-Dimensional Vortex Streets
In plane two-dimensional hydrodynamics, the equation governing the motion of an inviscid incompressible fluid can be written in terms of a streamfunction ψ(x,y,t) as∂∂t∇2ψ-∂(ψ,∇2ψ)∂(x,y)=0,
where ∂(a,b)/∂(x,y)=(∂a/∂x)(∂b/∂y)-(∂a/∂y)(∂b/∂x) is a Jacobian and ∇2=(∂2/∂x2)+(∂2/∂y2) is the two-dimensional Laplacian. The streamfunction (2.1) admits steady-state solutions of the form∇2ψ=F(ψ),
for any function ℱ, and a number of solutions are known for two-dimensional hydrodynamics. If we set ℱ≡0, (2.2) becomes the two-dimensional Laplace equation ∇2ψ=0, the fundamental solution of which has the formψ=κln(x-x0)2+(y-y0)2,
where (x-x0)2+(y-y0)2 is the distance from the fixed point x=x0, y=y0, which represents a point vortex of strength κ at x=x0,y=y0, for which ∇2ψ=0 except at the point itself with a circulation around the vortex of κ. Because ∇2ψ=0 except at the point vortex itself, it is possible to construct steady-state solutions to (2.2) consisting of more than one vortex although for such a combination of vortices to be a steady-state solution, each vortex much be located at a point where the velocities induced by the other vortices sum to zero. A number of well-known steady-state vortex street solutions are given in Table 1, and an overview of vortex streets can be found in standard texts on hydrodynamics such as [6, 7]. The first solution in the table, which dates back to [4], is an infinite row of identical vortices, each of strength κ, located along y=0 at the points x=±π,±3π,±5π,… which is sketched in Figure 1(a). If we combine a row of vortices of strengths κ at x=±π,±3π,±5π,… with a row of vortices of strengths -κ at x=0,±2π,±4π,…, we have a row of vortices of alternating sign, which is sketched in Figure 1(b). There are also two combinations due to von Kármán [5] which propagate at constant speed in the negative x-direction: the symmetrical double row and the staggered double row, both of which consist of two rows of vortices, one along y=d and the other along y=-d, with the symmetrical row propagating at speed c=(κ/2π)cothd and the staggered row at speed c=(κ/2π)tanhd. It is straightforward to verify that the single rows in Table 1 are stationary and that the double rows propagate at constant speed. To calculate the velocity of a vortex in a vortex street, we must subtract the contribution for an isolated vortex from the streamfunction for the vortex street and then evaluate the resulting expression at the location of the vortex, so that, for example, for the vortex at (x,y)=(π,0) in the corotating single row,ux=κ2πlimx→π,y→0∂∂y(ln[coshy+cosx]-ln[y2+(x-π)2])=0,uy=-κ2πlimx→π,y→0∂∂x(ln[coshy+cosx]-ln[y2+(x-π)2])=0.
With this approach, it can be shown that the corotating and counterrotating rows are stationary, while the symmetrical and staggered double rows propagate in the x-direction with velocities -(κ/2π)cothd and -(κ/2π)tanhd, respectively.
Some two-dimensional vortex streets.
Single corotating row
ψ=κ2πln[coshy+cosx]
Single counterrotating row
ψ=κ2πln[coshy+cosxcoshy-cosx]
Symmetrical double row
ψ=κ2πln[cosh(y-d)+cos(x+ct)cosh(y+d)+cos(x+ct)]
Staggered double row
ψ=κ2πln[cosh(y-d)+cos(x+ct)cosh(y+d)-cos(x+ct)]
Point vortices on the plane: (a) a row of corotating vortices; (b) a row of counterrotating vortices.
For the plane case, the corotating single row was generalized to smooth finite-amplitude vortices satisfying Liouville’s equation by [17] and the counterrotating row to finite-amplitude vortices satisfying the sinh-Poisson equation by [18].
2.2. Vortex Streets on a Sphere
If we have a two-dimensional vortex street which is periodic in at least one direction, such as those in Table 1, it is possible to transport that street to the surface of a sphere, which to the best of our knowledge has not been done previously. In what follows, we will work in spherical polars (r,θ,ϕ) with associated velocities (ur,uθ,uϕ). If we consider only motions on the surface of a sphere of radius R, we will have r=R and the radial velocity will be zero, ur=0. If we introduce a streamfunction via uθ=(∂ψ/∂ϕ)cscθ and uϕ=-∂ψ/∂θ, the inviscid incompressible equations of motion become∂∂t∇2ψ-1Rsinθ∂(ψ,∇2ψ)∂(θ,ϕ)=0,
together with a radial pressure gradient ∂p/∂r=R-1(uθ2+uϕ2). In the above, ∇2≡R-2[csc2θ(∂2/∂ϕ2)+cotθ(∂/∂θ)+(∂2/∂θ2)] is the Laplacian on the sphere, known as the Laplace-Beltrami operator. The streamfunction (2.5) admits steady-state solutions of the form∇2ψ=F(ψ).
This formulation is of course not new and has been used by a number of authors to study vortices on a sphere in the past. At this point, however, we will depart from previous studies and make the change of variable η=lntan(θ/2), so that the interval 0≤θ≤π is mapped onto -∞<η<∞. With this change of variable, ∇2≡R-2cosh2η[(∂2/∂ϕ2)+(∂2/∂η2)]. Since we can write the surface of the sphere as 0≤θ≤π,0≤ϕ≤2π, or equivalently -∞<η<∞,0≤ϕ≤2π, it follows that if ψ=Ψ(x,y) is a point vortex solution on the plane which is 2π-periodic in x, then the transformation x→ϕ,y→η will produce a point vortex solution on the sphere, ψ=Ψ(ϕ,lntan(θ/2)). If we apply this transformation to the two-dimensional vortex streets in Table 1, we obtain the corresponding vortex streets on a sphere given in Table 2. We will examine a couple of these vortex streets in more detail.
It is straightforward to verify that the single rows in Table 2 are stationary and that the double rows propagate at constant speed. To calculate the velocity of a vortex in a vortex street, we must subtract the contribution of that vortex itself from the streamfunctions in Table 2 and then evaluate the resulting expression at the location of the vortex. Since dη/dθ=cscθ,uθ=(∂ψ/∂ϕ)cscθ and uϕ=-(∂ψ/∂η)cscθ, and therefore it follows from the velocities of the two-dimensional vortex streets in Section 2.1 that the corotating and counterrotating rows are stationary, while the symmetrical and staggered double rows propagate in the ϕ-direction with velocities κcothdcoshd/2π and κsinhd/2π, respectively, since dη/dθ|η=d=coshd and are therefore rotating steadily about the axis.
2.2.1. Single Corotating Row
The streamfunction and velocity components for a single row of identical vortices are given byψ=κ2πln[cosh(mlntanθ2)+cosmϕ]=κ2πln[12(tanmθ2+cotmθ2)+cosmϕ],(uθ,uϕ)=-m(sinmϕ,sinh(mlntan(θ/2)))sinθ[cosh(mlntan(θ/2))+cosmϕ]=-m(sinmϕ,(1/2)(tanm(θ/2)-cotm(θ/2)))sinθ[(1/2)(tanm(θ/2)+cotm(θ/2))+cosmϕ],
for m=1,2,3,…. This streamfunction is plotted in Figure 2 for m=4. When m=1, the streamfunction (2.7) reduces to ψ=ln[cosϕ+cscθ], which is part of a larger family of solutions ψ=ln[cosϕ+(1/2)(ctan(θ/2)+(1/c)cot(θ/2))], while when m=2, we have ψ=ln[2cot2θ+2cos2ϕ].
Corotating point vortices on a sphere with m=4 viewed from (a) positive y-axis; (b) positive z-axis; (c) in ϕ-θ space.
The points θ=0 and π, which are the poles of the sphere, correspond to y=-∞ and +∞ on the plane, and ψ→∞ as θ→0+ and π-. We can see from Figure 2(b) that the flow at each of the poles looks like a vortex, and (2.7) confirms that there are point vortices at the two poles, with ψ~-(κm/2π)lnθ as θ→0 and -(κm/2π)ln(π-θ) as θ→π. Close to the poles, we have uϕ~-mκ/2πθ as θ→0 and mκ/2π(π-θ) as θ→π, so that uϕ is singular at the poles, while uθ~-(mκθm-1/2mπ)sinmϕ as θ→0 and -(mκ(π-θ)m-1/2mπ)sinmϕ as θ→π, so that uθ→-(κ/2π)sinϕ as θ→0,π for m=1 and uθ→0 as θ→0,π for m>1.
2.2.2. Single Counterrotating Row
The streamfunction and velocity components for a single row of alternating vortices are given byψ=κ2πln[cosh(mlntan(θ/2))+cosmϕcosh(mlntan(θ/2))-cosmϕ]=κ2πln[(1/2)(tanm(θ/2)+cotm(θ/2))+cosmϕ(1/2)(tanm(θ/2)+cotm(θ/2))-cosmϕ],(uθ,uϕ)=4m(-sinmϕcosh(mlntan(θ/2)),cosmϕsinh(mlntan(θ/2)))sinθ[cosh(2mlntan(θ/2))-cos2mϕ]=2m(-sinmϕ(tanm(θ/2)+cotm(θ/2)),cosmϕ(tanm(θ/2)-cotm(θ/2)))sinθ[(1/2)(tan2m(θ/2)+cot2m(θ/2))-cos2mϕ],
for m=1,2,3,…. This streamfunction is plotted in Figure 3 for m=2. When m=1, the streamfunction (2.8) reduces to ψ=ln[(1+cosϕsinθ)/(1-cosϕsinθ)], which is the point vortex limit of a solution given in [19], while when m=2, we have ψ=ln[(cos2θ+cos2ϕsin2θ)/(1-cos2ϕsin2θ)].
Counterrotating point vortices on a sphere with m=2 viewed from (a) positive y-axis; (b) positive z-axis; (c) in ϕ-θ space.
The spherical arcs ϕ=π/2m,3π/2m,…,(2m-1)π/2m are streamlines since ψ=0 on these lines. The points θ=0 and π, which are the poles of the sphere, correspond to y=-∞ and +∞ on the plane, and ψ→0 as θ→0+ and π-. We can see from Figure 3(b) that the flow at each of the poles looks like flow in a corner of angle π/2m, with the streamlines ϕ=π/2m,3π/2m,…,(2m-1)π/2m acting like rigid walls, and (2.8) confirms that ψ~(κθm/π2m-1)cosmϕ as θ→0 and (κ(π-θ)m/π2m-1)cosmϕ as θ→π. Close to the poles, we have uϕ~-(mκθm-1/2m-1π)cosmϕ as θ→0 and (mκ(π-θ)m-1/2m-1π)cosmϕ as θ→π, while uθ~-(mκθm-1/2m-1π)sinmϕ as θ→0 and -(mκ(π-θ)m-1/2m-1π)sinmϕ as θ→π, so that for m=1,uϕ~-(κ/π)cosϕ as θ→0 and (κ/π)cosϕ as θ→π and uθ~-(κ/π)sinϕ as θ→0 or π, and for m>1, both uϕ and uθ→0 as θ→0 or π.
2.2.3. Symmetrical and Staggered Double Rows
Near to the poles, the streamfunction for the symmetrical double row, ψ~(κ/π)(md-21-msinhmdcosmϕθm) as θ→0 and -(κ/π)(md-21-msinhmdcosmϕ(π-θ)m) as θ→π, with uϕ~(mκθm-1/2m-1π)sinhmdcosmϕ as θ→0 and (mκ(π-θ)m-1/2m-1π)sinhmdcosmϕ as θ→π, while uθ~(mκθm-1/2m-1π)sinhmdsinmϕ as θ→0 and -(mκ(π-θ)m-1/2m-1π)sinhmdsinmϕ as θ→π, so that for m=1,uϕ~(κ/2π)sinhdcosϕ as θ→0 or π and uθ~-(κ/π)sinhdsinϕ as θ→0 and -(κ/π)sinhdsinϕ as θ→π, and for m>1, both uϕ and uθ→0 as θ→0 or π.
Near to the poles, the streamfunction for the staggered double row, ψ~(κ/π)(md+21-mcoshmdcosmϕθm) as θ→0 and -(κ/π)(md-21-mcoshmdcosmϕ(π-θ)m) as θ→π, with uϕ~-(mκθm-1/2m-1π)coshmdcosmϕ as θ→0 and (mκ(π-θ)m-1/2m-1π)coshmdcosmϕ as θ→π, while uθ~-(mκθm-1/2m-1π)coshmdsinmϕ as θ→0 and -(mκ(π-θ)m-1/2m-1π)coshmdsinmϕ as θ→π, so that for m=1,uϕ~-(κ/π)coshdcosϕ as θ→0 and (κ/π)coshdcosϕ as θ→π and uθ~-(κ/2π)coshdsinϕ as θ→0 or π, and for m>1, both uϕ and uθ→0 as θ→0 or π.
2.2.4. Gauss’ Constraint on the Vorticity
Since a sphere is a closed compact surface, it follows from Gauss’ theorem that the integral of the scalar vorticity field over the spherical surface must be zero. This is a global constraint on the vorticity distribution. In order to satisfy this constraint and simultaneously have an irrotational flow, each point vortex must be counterbalanced by another point vortex on the sphere. By inspection, the single counterrotating row and the symmetrical and staggered double rows satisfy this constraint because for these vortex streets the total vorticity on the spherical surface is the sum of the circulations of the point vortices, and for each point vortex, there is a point vortex of opposite strength, so that the circulations sum to zero. For the single corotating row (2.7), the situation is slightly more complex. There are m point vortices along the equator θ=π/2, and close to the vortex at θ=π/2 and ϕ=π/m, the streamfunction behaves like (κ/π)ln(θ-(π/2))2+(ϕ-(π/m))2, with similar behavior at the other vortices along the equator. There are also two point vortices at the poles: close to the vortex at θ=0, the streamfunction behaves like -(κm/2π)lnθ, while close to the vortex at θ=π, the streamfunction behaves like -(κm/2π)ln(π-θ), so that the circulations sum to zero and the constraint is satisfied for the single corotating row, and it would appear that the two polar vortices are generated to satisfy this constraint.
3. Discussion
In the previous section, we saw that if we had a vortex street ψ=Ψ(x,y) on the plane which was 2π-periodic in x, then the transformation x→ϕ, y→lntan(θ/2) would produce a vortex street on the sphere, ψ=Ψ(ϕ,lntan(θ/2)). We would mention that in [20], we used a much simpler transformation to transport vortices from the plane to the cylinder. Perhaps surprisingly, this very simple approach of transporting vortex streets from the plane to the sphere does not appear to have been used previously, with most previous studies having instead used a stereographic projection, an approach which Lamb [6] mentions was used by Kirchhoff as far back as 1875 to study electrical conduction in a spherical sheet. As an illustration, we applied this transformation to the two-dimensional vortex streets in Table 1 to obtain the corresponding vortex streets on a sphere given in Table 2, and the vortex streets corresponding to rows of corotating and counterrotating vortices were plotted in Figures 2 and 3. For the corotating row, two polar vortices are generated, presumably to satisfy the Gauss constraint on the vorticity.
As to the direction of possible future research, we mentioned in Section 1 that for the plane case, some of the vortex streets have been generalized to smooth finite amplitude vortices, resulting in Stuart vortices [17] and Mallier-Maslowe vortices [18]. We saw in Section 2.2 that when m=1, the streamfunction (2.8) for the counterrotating row reduced to the point vortex limit of a solution given in [19], where we presented some finite-amplitude vortex solutions for the sphere, so conversely it follows that the m=1 counterrotating case can be generalized to smooth finite-amplitude vortices. It would be interesting to see if the remaining vortex streets presented here could be generalized as well although we note that this is somewhat harder for the sphere, because the equation for nonlinear solutions on the sphere, [(∂2/∂ϕ2)+(∂2/∂η2)]ψ=r2sech2ηℱ(ψ), has a sech2η factor on the right hand side which is not present in the corresponding equation on the plane, [(∂2/∂x2)+(∂2/∂y2)]ψ=ℱ(ψ). Perhaps because of this additional factor, very few exact smooth solutions are known for flow on a sphere, and apart from [19], the only exact distributed vortex equilibria (as opposed to point vortices) on a sphere appear to be the exact solutions on a rotating sphere presented in [21, 22] and two studies by Crowdy, one involving a generalization of Stuart vortices to the sphere [23], the other involving a combination of vortex patches and point vortices [24].
We would not claim that the present study is exhaustive, and another possible direction for future research was suggested by an anonymous referee, who wondered whether any vortex street solutions might be found consisting of double rows with different strengths.
HelmholtzH.Über Integrale der hydrodynamishen Gleichgungen entsprechen1858552555KirchhoffG.1877Leipzig, GermanyTeubnerMechanik, Lecture 20LinC. C.1943Toronto, CanadaUniversity of Toronto Press39University of Toronto Studies, Applied Mathematics Series, no. 50008204MaxwellJ. C.1873Oxford, UKvon KármánT.Über den mechanismus des widerstandes, den ein bewegter Körper in einer Flussigkeit erfahrt1912547556LambH.1932Cambridge, UKDoverMilne-ThomsonL. M.1968London, UKMacMillanGromekaI. S.On vortex motions of liquid on a sphere. Uchenye Zapiski Imperatorskogo Kazanskogo Universiteta18853202236See also Collected Papers, Moscow Akademii Nauk, USSR, 1952HallyD.Stability of streets of vortices on surfaces of revolution with a reflection symmetry198021121121710.1063/1.524322556289ZBL0446.76027BogomolovV. A.Two-dimensional fluid dynamics on a sphere1979151822KimuraY.Vortex motion on surfaces with constant curvature1999455198124525910.1098/rspa.1999.03111700500KimuraY.OkamotoH.Vortex motion on a sphere198756124203420610.1143/JPSJ.56.4203927662ArefH.NewtonP.StremlerM.TokiedaT.VainchteinD. L.Vortex crystals20033917910.1016/S0065-2156(02)39001-XDibattistaM. T.PolvaniL. M.Barotropic vortex pairs on a rotating sphere199835810713310.1017/S00221120970081001616092ZBL0908.76096DritschelD. G.Contour dynamics/surgery on the sphere19887947748310.1016/0021-9991(88)90028-9ZBL0655.76018PolvaniL. M.DritschelD. G.Wave and vortex dynamics on the surface of a sphere1993255356410.1017/S00221120930023811244223ZBL0793.76022StuartJ. T.On finite amplitude oscillations in laminar mixing layers19672941744010.1017/S0022112067000941ZBL0152.45403MallierR.MasloweS. A.A row of counter-rotating vortices1993541074107510.1063/1.8586221208239ZBL0778.76022AlobaidiG.HaslamM. C.MallierR.Vortices on a sphere20061143573642285385ZBL1116.76016HaslamM. C.MallierR.Vortices on a cylinder20031572087208810.1063/1.15821811987480ZBL1186.76215VerkleyW. T. M.The construction of barotropic modons on a sphere198441162492250410.1175/1520-0469(1984)041<2492:TCOBMO>2.0.CO;2768284VerkleyW. T. M.Nonlinear structures with multivalued (q, θ) relationships—exact solutions of the barotropic vorticity equation on a sphere199369779410.1080/03091929308203575CrowdyD. G.Stuart vortices on a sphere200449838140210.1017/S00221120030070432256924ZBL1059.76012CrowdyD.ClokeM.Analytical solutions for distributed multipolar vortex equilibria on a sphere2003151223410.1063/1.15217271971791