Next Article in Journal
Anatomical Considerations and Study of the Fractal Dimension around the Posterior Superior Alveolar Artery
Next Article in Special Issue
Quadratic Stability of Non-Linear Systems Modeled with Norm Bounded Linear Differential Inclusions
Previous Article in Journal / Special Issue
Analysis of Heat and Mass Transfer for Second-Order Slip Flow on a Thin Needle Using a Two-Phase Nanofluid Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Motion of a Point Vortex in Multiply-Connected Polygonal Domains

by
El Mostafa Kalmoun
*,
Mohamed M. S. Nasser
and
Khalifa A. Hazaa
Department of Mathematics, Statistics and Physics, Qatar University, Doha 2713, Qatar
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(7), 1175; https://doi.org/10.3390/sym12071175
Submission received: 22 December 2019 / Revised: 12 January 2020 / Accepted: 13 January 2020 / Published: 16 July 2020
(This article belongs to the Special Issue Symmetry in Numerical Analysis and Numerical Methods)

Abstract

:
We study the motion of a single point vortex in simply- and multiply-connected polygonal domains. In the case of multiply-connected domains, the polygonal obstacles can be viewed as the cross-sections of 3D polygonal cylinders. First, we utilize conformal mappings to transfer the polygonal domains onto circular domains. Then, we employ the Schottky-Klein prime function to compute the Hamiltonian governing the point vortex motion in circular domains. We compare between the topological structures of the contour lines of the Hamiltonian in symmetric and asymmetric domains. Special attention is paid to the interaction of point vortex trajectories with the polygonal obstacles. In this context, we discuss the effect of symmetry breaking, and obstacle location and shape on the behavior of vortex motion.

1. Introduction

Since the seminal work by Helmholtz [1], 2D point vortex motion has sparked intense research due to its important applications in fluid dynamics and geophysics [2,3,4,5]. A point vortex is a model of potential flow in which vorticity of the surrounding flow is present only at a single point. In addition to the various physical applications, the 2D point vortex model, despite its simplicity, provides a mathematics playground in the words of Aref [2] that links a variety of mathematical areas such as dynamical systems, ordinary and partial differential equations, Hamiltonian dynamics, theory of polynomials, and elliptic functions, to name a few.
The dynamics of point vortices is governed by a Hamiltonian system, which can be derived from the Euler equations describing the velocity field and pressure of a 2D incompressible and inviscid flow. If the fluid domain has no boundaries, e.g., the whole complex plane, the Hamiltonian of the governing dynamical system is written in terms of the complex Green function [5]. The presence of solid boundaries might break symmetries that would otherwise exist and thus the Hamiltonian system is constrained by the fluid non-penetration condition expressed by zero normal velocity on all solid boundaries. For simple geometries containing special symmetries such as the upper-half plane, a channel between two infinite walls or the exterior of a circular cylinder, an explicit formula of the Hamiltonian of a point vortex can be obtained by the classical method of images in which solid boundaries are treated as streamlines of a flow field. This is usually combined with conformal mapping ideas, as can be seen in [4,5] for a variety of examples. Similarly, in bounded and simply-connected domains, basic results concerning single vortex motion are well established [6]. In this case, the generalized Hamiltonian is known as the Kirchhoff–Routh path function [7,8].
As many fluid domains are not simply-connected (e.g., in aerodynamics), a growing interest has centered around studying N-point vortex dynamics in multiply-connected domains. Lin [9] early proved the existence and uniqueness of the Kirchhoff–Routh path function in multiply-connected domains but no explicit formula was given. Recently, Johnson and McDonald [10] studied the motion of one point vortex in doubly- and triply-connected domains exterior to either one or two circular cylinders using elliptic functions, provided that zero circulations are imposed around the circles. Under the same boundary conditions, Crowdy and Marshall [11] presented an analytical formula for the Hamiltonian in multiply-connected circular domains by applying methods of complex function theory. In particular, their analysis makes use of a special transcendental function known as the Schottky-Klein prime function. The formula can be employed for a general multiply-connected domain as long as the conformal mapping onto the canonical circular domain is known. Indeed, this formula was already applied to study point vortex motion around an arbitrary finite number of circular cylinders [12] and through gaps in walls [13]. When the multiply-connected domain is symmetric with respect to the real axis, Sakajo [14] showed that the two point vortices’ motion is reduced to that of a single point vortex in a multiply-connected semicircle. Thus, the Crowdy and Marshall analytic formula was again utilized to plot the trajectories of two point vortices for several circular domains.
In this paper, we study the motion of a single point vortex in simply- and multiply-connected polygonal domains. In the case of multiply-connected domains, the polygonal obstacles can be viewed as the cross-sections of polygonal cylinders in 3D. Our method first makes use of conformal mappings to transfer the polygonal domain onto the circular domain where the Crowdy and Marshall analytic formula can be employed. The objective of our paper is to compare the topological structures of the contour lines of the Hamiltonian in circular and polygonal domains. We also aim to emphasize how these trajectories interact with the polygonal obstacles. In this context, we discuss the effect of symmetry breaking, and obstacle location and shape on the behavior of vortex motion.

2. Computing the Hamiltonian of a Point Vortex Motion

We present in this section a numerical method for computing the Hamiltonian of a point vortex motion in multiply-connected polygonal domains. The method is based on using conformal mappings to map polygonal domains onto circular ones, and then appealing to the analytic formula presented in [11] for the Hamiltonian in multiply-connected circular domains.
Let G denote the multiply-connected domain exterior to m non-overlapping polygons Γ 1 , , Γ m (clockwise orientated) and interior to a polygon Γ m + 1 (in the counterclockwise orientation). Such a domain is called a polygonal domain. We assume that no corner of these polygons is a cusp. Let us suppose that the interior of G is occupied with incompressible fluid. The flow is assumed to undergo irrotational motion except for a single point vortex. Therefore, the circulations around the polygonal obstacles Γ 1 , , Γ m are zero, and the point vortex has a given circulation χ .

2.1. The Conformal Mapping

It is known that there exists a conformal mapping ζ = f ( z ) from the polygonal domain G onto a multiply-connected circular domain D interior to the unit circle and exterior to m non-overlapping circles [15]. We denote the inner m circles by C k , k = 1 , 2 , , m and the unit circle by C m + 1 . By Carathéodory’s theorem ([15] p. 111), the mapping function f can be extended to G ¯ = G G in such a way that the polygon Γ k is mapped onto the circle C k , for  k = 1 , 2 , , m + 1 . Note that G denotes the boundary of the domain G, which consists of m + 1 piecewise smooth Jordan curves. For the uniqueness of the conformal mapping, we fix a point α in G, and we assume that f ( α ) = 0 and f ( α ) > 0 .
The centers and radii of the inner circles are unknown and uniquely determined by the polygonal domain G. Thus, computing the conformal mapping f requires also computing the centers and radii in the circular domain. In this study, the conformal mapping f was computed using the MATLAB toolbox (PlgCirMap) presented in [16].

2.2. The Hamiltonian for Circular Domains

The trajectories of a point vortex in the multiply-connected circular domain D can be written as the contours of the Hamiltonian, also referred to as the Kirchhoff–Routh path function [11]. The existence and uniqueness of such a function was established by Lin [9].
Let H D ( ζ , ζ ¯ ) denote the Hamiltonian for the motion of a single point vortex in the circular domain D. An explicit formula for H D ( ζ , ζ ¯ ) was derived by Crowdy and Marshall [11],
H D ( ζ , ζ ¯ ) = χ 2 8 π log 1 ζ 2 ω ^ ( ζ , ζ ) ω ^ ¯ ( 1 / ζ , 1 / ζ ) ω ( ζ , 1 / ζ ¯ ) ω ^ ¯ ( 1 / ζ , ζ ¯ ) ,
where ζ ¯ is the complex conjugate of ζ , ω ( ζ , ξ ) is the Schottky-Klein prime function (prime function, for short) associated with the circular domain D, and the function ω ^ ( ζ , ξ ) is defined through
ω ( ζ , ξ ) = ( ζ ξ ) ω ^ ( ζ , ξ ) , ζ , ξ D ,
and ω ^ ¯ ( ζ , ξ ) = ω ^ ( ζ ¯ , ξ ¯ ) ¯ .
The Schottky-Klein prime function is an important special transcendental function. Due to the works of Crowdy with various collaborators [17] in the last two decades, the prime function has become a key tool to construct analytic formulas in solving several problems involving multiply-connected circular domains. For more details on the properties of the prime function, its applications and numerical computation, the reader is referred to the review article [18] and the recent monograph [17].
A simplified form of Equation (1) was proposed by Sakajo [14],
H D ( ζ , ζ ¯ ) = χ 2 4 π log ω ^ ( ζ , ζ ) ζ ω ( ζ , 1 / ζ ¯ ) .
In the two numerical methods presented in [18], the prime function is normalized with w ( ζ , ζ ) = 1 , which yields
H D ( ζ , ζ ¯ ) = χ 2 4 π log ζ ω ( ζ , 1 / ζ ¯ ) .
By making use of an important property of the prime function [17] that states
ω ( ζ , 1 / ζ ¯ ) = 0 for ζ k = 1 m + 1 C k ,
it follows from Equation (2) that H D ( ζ , ζ ¯ ) as ζ approaches D .

2.3. The Hamiltonian for Polygonal Domains

By the help of conformal mappings, the above exact formula in Equation (2) can be extended to the general multiply-connected domains [11]. In fact, using the conformal mapping ζ = f ( z ) to map the polygonal domain G onto the circular domain D, the Hamiltonian H G ( z , z ¯ ) of a single vortex motion of circulation χ in the bounded polygonal domain G is given by [11]
H G ( z , z ¯ ) = χ 2 4 π log ζ ω ( ζ , 1 / ζ ¯ ) + χ 2 4 π log ( f 1 ) ( ζ ) ,
which can be simplified as
H G ( z , z ¯ ) = χ 2 4 π log | f ( z ) | ω f ( z ) , 1 / f ( z ) ¯ f ( z ) .
Equation (3) implies that H G ( z , z ¯ ) as z approaches G as well. In the sequel, we assume that the point vortex has a circulation χ = 1 .

2.4. Numerical Implementation

Recently, a MATLAB toolbox (PlgCirMap) for computing the conformal mapping ζ = f ( z ) , from the polygonal domain G onto the circular domain D, and its inverse z = f 1 ( ζ ) was presented [16]. The method used in [16] is based on a fast implementation of Koebe’s iterative method [19] using the boundary integral equation with the generalized Neumann kernel [20] and fast multipole method [21]. The toolbox also works if G is simply-connected ( m = 1 ), and hence D is the unit disk. On the other hand, to compute the prime function, two numerical methods are presented in [18] with MATLAB codes available online. Their first procedure is a spectral method based on global relations, while the second depends on the boundary integral equation with the generalized Neumann kernel.
In this study, the Hamiltonian in Equation (4) was computed by a combination of the MATLAB toolbox PlgCirMap and one of the two numerical methods presented in [18]. The outline of our procedure is summarized in Algorithm 1.
Algorithm 1: Computing the Hamiltonian H G for the polygonal domain G.
1:
Define the vertices of the polygons Γ k , k = 1 , , m + 1 , as a cell array ver{k}.
2:
Choose an auxiliary point alpha in G.
3:
Use the MATLAB function plgcirmap to compute
f = plgcirmap(ver, alpha)
where f is a MATLAB struct with several fields.
4:
Extract the centers and radii of the inner circles in D from f by f.cent(1:m) and f.rad(1:m).
5:
Discretize the domain G by a matrix of points Z.
6:
Call the functions evalu and evalud from the toolbox PlgCirMap to compute the values of the conformal mapping fZ and its derivative dfZ on the meshgrid Z:
fZ = evalu(f,Z,’d’);  dfZ = evalud(f,Z,’d’)
7:
Compute the prime function by applying one of the two methods presented in [18].
8:
Use Equation (4) to compute the Hamiltonian.
Note that the point α in Step 2 of Algorithm 1 has no effect on the numerical computation of the Hamiltonian as long as it is sufficiently far from the boundary of G.

3. Simply-Connected Domains

When D is the unit disk, the associated prime function is simply
ω ( ζ , α ) = ζ α .
Thus, the Hamiltonian in Equation (2) of a single vortex motion of circulation χ in the unit disk D reduces to [11]
H D ( ζ , ζ ¯ ) = 1 4 π log | 1 ζ ζ ¯ | = 1 4 π log ( 1 | ζ | 2 ) .
Hence, the contour lines of the Hamiltonian are circles with a center point at the origin. It is clear that 0 < 1 | ζ | 2 1 when ζ is in the unit disk D and 1 | ζ | 2 0 + as ζ approaches D . Hence H D ( ζ , ζ ¯ ) as ζ approaches D and H D ( 0 , 0 ) = 0 is the unique global maximum.
Similarly, the Hamiltonian in Equation (4) in the bounded simply-connected polygonal domain G is given by [11]
H G ( z , z ¯ ) = 1 4 π log | f ( z ) | f ( z ) 1 / f ( z ) ¯ f ( z ) = 1 4 π log 1 | f ( z ) | 2 | f ( z ) | ,
where f is the above described conformal mapping from G onto D. Note also that H G ( z , z ¯ ) as z approaches G .
For convex domains, Gustafsson [22] proved that the Hamiltonian is strictly concave, and therefore the point vortex has a unique equilibrium point, which is the global maximum of H G . In this case, contour plots of the Hamiltonian H G , which correspond to vortex trajectories, span the domain G by encircling the unique critical point in a similar topological way as in the disk D. When passing to non convex domains, the uniqueness of critical points does not hold in general, as noticed by Gustafsson [23], who provided a counter-example for star-shaped domains. In the next subsection, we present a numerical example that illustrates this fact.

3.1. A Star-Shaped Domain

The vertices of the star-shaped domain in Figure 1 are 1 + 0.5 i , 1 + i , 0.5 + i , a i , 0.5 + i , 1 + i , 1 + 0.5 i , b , 1 0.5 i , 1 i , 0.5 i , a i , 0.5 i , 1 i , 1 0.5 i , and b for three different values of ( a , b ) . The X-shaped polygon in Figure 1a has ( a , b ) = ( 0.5 , 0.5 ) , which makes the four concave corners located on the same circle of center zero and radius 0.5 . Similar to the unit disk case, the Hamiltonian in this polygonal domain has a unique critical point, which is a center at the origin surrounded by closed orbits.
If the two vertices on the imaginary axis are displaced by 0.25 unit toward each other, so that they are not longer located on the same circle with the other two vertices on the real axis, then we have more than one critical point (see Figure 1b). Indeed, the neutrally stable center point becomes a saddle, which is connected to itself by a pair of symmetric homoclinic orbits forming together a separatrix. Two new vortex centers appear on the real axis at the same distance from the saddle point zero. The Hamiltonian has four other centers where each point is surrounded by a homoclinic orbit connecting a saddle point. These four hyperbolic points are themselves connected by a heteroclinic loop. Observe also how all equilibria of this system are distributed in a way that reflects line symmetries in the domain.
Another worth commenting dynamical feature in this example is the reappearance of the center point at the origin instead of the homoclinic structure when the four concave corners are placed again on a circle of center zero and radius 0.25 . However, we do not recover the uniqueness of critical points since the previous outer homoclinic-heteroclinic pattern still occurs (see Figure 1c).

3.2. Non-Convex Polygonal Domains

We consider in Figure 2, Figure 3 and Figure 4 three other examples of general bounded simply-connected polygonal domains that are non convex. The C-shaped polygon in Figure 2 has eight vertices at 1 + i , 1 + i , 1 i , 1 i , 1 0.5 i , a 0.5 i , a + 0.5 i , 1 + 0.5 i ; the H-shaped polygon in Figure 3 has twelve vertices at 1 + i , 0.5 + i , 0.5 + a i , 0.5 + a i , 0.5 + i , 1 + i , 1 i , 0.5 i , 0.5 a i , 0.5 a i , 0.5 i , 1 i ; and the twelve vertices of the E-shaped polygon in Figure 4 are 1 i , 1 0.6 i , 0.6 i , 0.2 i , 1 0.2 i , 1 + 0.2 i , 0.2 i , 0.6 i , 1 + 0.6 i , 1 + i , a + i , a i .
The common dynamical feature between the three shapes is that the unique stable equilibrium point is always located at the center of the maximal rectangle that has antiknobs in the boundary. When this maximal rectangle becomes narrower in either the C-shaped polygon or the H-shaped domain, the center point changes to a saddle connected by two homoclinic trajectories that surround two new vortex centers. Observe that this not the case in the E-shaped domain, as this center point still occurs but it is now surrounded by a heteroclinic orbit. In this example, two new saddle points and two new vortex centers are born making each group of elliptic and hyperbolic points vertically aligned.

4. Doubly-Connected Domains

Sakajo [14] studied the motion of a single point vortex in doubly-connected circular domains. If the center of the inner circle is not the origin, then the Hamiltonian always has one elliptic center point and one saddle point connected to itself by two homoclinic vortex trajectories (see Figure 5b). In this section, we discuss what happens to the dynamics of point vortex when we replace circular domains with polygonal domains.
We assume the external boundary component of the domain G is composed of the square with vertices ± 1 ± i . For the inner obstacle shape, we consider four cases: a square, a rotated square, an equilateral triangle, and a hexagon. In Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11, we show vortex motion in these doubly-connected polygonal domains by plotting some contours of the Hamiltonian for different positions of the polygonal obstacle center.

4.1. Concentric Domains

We first consider the case of concentric domains in Figure 6a, Figure 7a, Figure 8a, and Figure 11a. Unlike the annulus domain (Figure 5a) where there is no center or saddle, the polygonal structure gives rise to several equilibria of this type. The number of vortex centers of the Hamiltonian is the same as that of saddle points, which is again in accordance with Morse theory. The elliptic equilibria are surrounded by closed vortex trajectories and the hyperbolic ones are connected by either homoclinic or heteroclinic orbits.
In Table 1, we show the number of saddle points and how they are linked by vortex trajectories in concentric polygonal doubly-connected domains. By looking at these numbers, there are aspects to be highlighted regarding the topology of vortex trajectories near the obstacle boundaries. For a square obstacle (standard and rotated), we observe a pure heteroclinic structure, which is not the case for triangle and hexagon obstacles. Near the square obstacle in Figure 6a, there occur eight saddle points that are connected to each other by eight heteroclinic loops. At first, this number seems to result from the total number of reflectional and rotational symmetries or the number of corner points in the domain. However, when we rotate the square obstacle by an angle of π / 4 in Figure 7a so that its sides become perpendicular to the diagonals of the outer square, this number is reduced by half, although the total symmetries and corner points in the domain remain the same.
If we change the obstacle shape to an equilateral triangle as in Figure 8a, there are only three saddle points that are connected by four homoclinic and one heteroclinic loops. For the case of a regular hexagon obstacle in Figure 11a, the number of saddles increases again to four points connected by four homoclinic and two heteroclinic loops. This suggests that the number of equilibria and topology pattern of vortex motion might be more related to the interaction between the vertices of the polygon obstacle and the outer square boundaries.
In the next subsections, we are interested in the manner the vortex motion changes due to obstacle movement.

4.2. A Square Obstacle

After breaking the symmetry with respect to the outer square in Figure 6b by a perturbation of the center p of the square obstacle with a translation of 0.01 unit along the imaginary axis, the topology pattern of vortex motion changes dramatically. Only two saddle points which are situated on the side of the obstacle movement still occur. The saddle connection is made by two heteroclinic loops encircling two center points, which are symmetrically located on the imaginary axis with one above the obstacle and the other below it. The same configuration appears if we continue moving the square obstacle northward, as shown in Figure 6c,d.
On the other hand, if we choose to break the symmetry by shifting the obstacle center to the upper right corner of the outer square, as in Figure 6e, then a small perturbation changes again the topology pattern, but in a slightly different way. Besides the heteroclinic cycle containing two saddles in the direction of the obstacle movement as in the previous case, the other elliptic point, which is situated on the left lower side of the diagonal line y = x , is now the center of a new heteroclinic loop itself joining two new saddles. Each of the latter hyperbolic points belongs to one homoclinic orbit surrounding a vortex center, which gives a total of eight critical points.
When displacing the center of the obstacle to p = 0.11 + 0.11 i , as shown in Figure 6f, the vortex motion loses all its rest points except one saddle and one center, which are, respectively, located on the upper right and lower left sides of the diagonal line y = x . Observe also that a topological configuration similar to the vortex motion in Figure 6c,d appears but on the diagonal instead of the imaginary axis when p = 0.65 + 0.65 i (see Figure 6g). After the separation of the opposite corners in the outer and inner squares becomes smaller, as shown in Figure 6h, the motion dynamics consists of four saddles and three center points in the separation region plus a fourth elliptic point situated near the center of the outer square. It is worth emphasizing that homoclinic orbits are only present in the phase portraits in Figure 6e,f while all other cases contain pure heteroclinic cycles.

4.3. A Rotated Square Obstacle

We now examine the effect of symmetry breaking when the obstacle is the rotated square in Figure 7. The small displacement of the center along the imaginary axis in Figure 7b does not reduce the number of equilibria as in the square obstacle. However, this breakdown destroys the heteroclinic structure shown in Figure 7a, which confirms the structurally instability of heteroclinic cycles. Indeed, none of the four saddles is still an intersection of two heteroclinic loops. However, each saddle gives rise to the bifurcation of two homoclinic orbits that form a separatrix.
On the other hand, a perturbation of 2 × 0.01 toward the top right corner of the outer square, as done in Figure 7c, seems to damage less the heteroclinic structure since we still have two heteroclinic loops surrounding the two centers on the main diagonal line. The two other heteroclinic loops change to homoclinic orbits as they become separated from the two saddles near the diagonal line y = x .
Figure 7d reveals that when we move the obstacle center to p = 0.11 i , the vortex motion only involves homoclinic trajectories that intersect at two hyperbolic fixed points situated on the imaginary axis. Observe also that the two bottom homoclinic orbits encircle two elliptic points.
We next consider the point vortex motion when the rotated square obstacle interacts more closely with the outer square boundaries. As the obstacle gets closer to the upper outer boundary in Figure 7e, the lower saddle point changes to a center yielding a configuration similar to the circular domain in Figure 5b, while the homoclinic pattern is maintained. This seems related to the interaction between the obstacle top vertex and the nearest outer square side. Figure 7f confirms this observation as two saddle points occur when two obstacle vertices become closer to the outer boundary. Note also that the saddle connections are of a heteroclinic nature in this case.

4.4. An Equilateral Triangle Obstacle

In Figure 8, we show the effect of obstacle displacement on vortex motion for the case of the equilateral triangle. We note first that the imaginary axis is the unique line of reflectional symmetry in the starting configuration of this example; that is, when the obstacle center is at p = 0 (see Figure 8a). This configuration seems structurally more stable to obstacle movement to the upper side external boundary. Indeed, a small perturbation by a 0.01 unit does not change at all the topology pattern of vortex trajectories.
In Figure 9, we illustrate how a small perturbation of the triangle obstacle location in different directions impacts the topology pattern of vortex trajectories. We can see that the configuration does not change for a vertical translation of 0.05 while it does with the same translation on the opposite direction. On the other hand, we observe a change to the same configuration when a small perturbation is applied along the horizontal and diagonal directions (see Figure 9d–f).
Next, we explore closely how the topological pattern of vortex motion changes after displacement of the obstacle triangle. First, for the vertical translation to the upper side boundary, the bifurcation from the two lower saddles and three centers to a unique center point does not happen until an approximate displacement of 0.102 (see Figure 10a). Afterwards, the same configuration continues to hold when moving the obstacle close to the boundary, as shown in Figure 8c.
An interesting dynamical feature in this example arises if we move the obstacle on the other side of the imaginary axis. A merging of two pairs of homoclinic orbits and one heteroclinic loop to a new structure composed of three heteroclinic connections and one triheteroclinic loop occurs at a very small bifurcation location of the triangle center equal to p = 0.004985 (see Figure 10b). This configuration changes immediately after this value to the topology pattern shown in Figure 9c. Indeed, the upper saddle point gets separated from the two lower saddles and becomes an intersection of a pair of homoclinic orbits. The same topological pattern appears also for a displacement of 0.11 unit in Figure 8e. The upper saddle changes to a center point for an approximate bifurcation displacement of 0.145 yielding the pattern given in Figure 10c. This configuration stays the same all the way as the obstacle approaches the lower boundary (see Figure 8f).
Now, when we break the line symmetry with respect to the imaginary axis by moving the obstacle triangle with 0.01 unit eastward in Figure 8g, the heteroclinic loop disappears as each of the three saddles becomes an intersection point of homoclinic orbits. The inner saddle point disappears first at an approximate bifurcation displacement of 0.063 (Figure 10d), and then the other inner saddle changes to a center point at almost p = 0.102 (Figure 10e). Afterwards, the vortex motion maintains a unique unstable rest point until the obstacle comes close the right boundary (see Figure 8h–i).
Note also that, as mentioned above, a small perturbation along the line y = x yields a topologically equivalent pattern to the horizontal translation case. However, when the obstacle triangle approaches the upper right corner of the outer square, we get a different configuration to the horizontal movement case as the vortex motion keeps two unstable equilibria (see Figure 8j–k). It is worth mentioning that this does not happen when the obstacle comes close to the lower right corner. The configuration in this case is topologically equivalent to the horizontal translation, as shown in Figure 8I.

4.5. A Hexagon Obstacle

In this example, we consider how the displacement of the hexagon obstacle in Figure 11 impacts the point vortex motion.
First, the upper heteroclinic loop connecting the two saddles on the real axis splits into two homoclinic orbits after a small perturbation of the center p of the hexagon by 0.01 unit along the imaginary axis (see Figure 11b). The upper saddle on the imaginary axis becomes unique after the three other points disappear at a vertical displacement of 0.11 unit. When approaching the upper boundary side, this saddle changes to a center point of a heteroclinic loop linking two new saddles.
When we displace the center of the triangle by 0.01 unit eastward in Figure 11e, a topologically equivalent pattern to the ( 6 , 1 ) homoclinic-heteroclinic structure we got after a small vertical perturbation of p appears again. The only difference in the horizontal displacement is that the pattern appearing at p = 0.11 stays the same even after approaching the right side boundary.
On the other hand, a diagonal displacement of the hexagon obstacle changes the dynamics of the point vortex in a different manner. A small perturbation as in Figure 11h destroys the heteroclinic structure completely, and all four saddles become intersection points of homoclinic orbits. At p = 0.11 + 0.11 i , the pair of saddle points on the side of the obstacle movement are maintained with the same homoclinic structure, but the other pair of saddles and pair of centers disappear except one vortex center on the main diagonal line. This pattern remains topologically invariant as the obstacle approaches the upper right corner.

5. Multiply-Connected Domains

In his paper [14], Sakajo studied also the motion of a single point vortex in multiply-connected circular domains. The centers of the inner circles are assumed to be on the real axis. The contour plot of the Hamiltonian for these domains is given in Figure 12. Crowdy and Marshal [11] also considered several circular domains.
Here, we investigate vortex motion in several multiply-connected polygonal domains. First, we replace the circles in Figure 12 by squares where the outer one has vertices at ± 1 ± i . The contour plot of the Hamiltonian in this case is then shown in the top row of Figure 13. The topology pattern of vortex motion looks similar to the circular domain case. Looking at how the point vortex trajectories interact with the polygonal obstacles, we observe that this interaction depends on the distance between these obstacles. When they are close to each other, as in Figure 13a, we have two saddle points between them connected by a heteroclinic loop, compared to one saddle for the circular case. As the separation between the two square obstacles increases, we get a vortex merging from two saddles and one center to one saddle at which two homoclinic orbits intersect and surround the previous upper and lower centers separately (see Figure 13b). Increasing further the separation distance, for example to 0.8 unit in Figure 13c, we get another vortex merging from the last homoclinic structure to a vortex center located at zero. This configuration continues to hold when each of the two square obstacles approaches the external eastern and western boundaries, as shown in Figure 13d. The only difference to the circular domain case is that we have two saddle points to the right and two to the left of the square obstacles instead of one to the right and one to the left.
We consider also the case when the centers of the inner squares are on the diagonal line y = x , as shown in the bottom row of Figure 13. When these centers are located at ± p for p = 0.2 + 0.2 i or p = 0.3 + 0.3 i in Figure 13e,f, the same dynamical configuration occurs, which is topologically equivalent to the pattern in Figure 13b. Observe that all three saddle points are located on the main diagonal line. After displacing the square obstacles toward the top right and bottom left outer corners in Figure 13g, the homoclinic structure at the center disappears, and a unique vortex center is born in a similar way to the horizontal separation. As expected, when the obstacles gets closer to the corners in Figure 13h, the phase portrait between each obstacle and the external boundary becomes similar to the case of one square obstacle. The Hamiltonian has four saddles near each side and all eight points are joined by heteroclinic connections.
We examine now the effect of reducing reflectional symmetry to a unique line by placing the obstacle at different locations, as shown in the top row of Figure 14. We first place one obstacle centered at zero and the other once in the north and then in the right east (see Figure 14a,b). In these two examples, the dynamics restricted near the lower side of the line of reflectional symmetry consists of a ( 2 , 1 ) homoclinic–heteroclinic structure similar to the one square obstacle case. The Hamiltonian also has one saddle between the two obstacles and two saddles connected by a heteroclinic loop between the obstacle and the external boundary. In Figure 14c,d, we consider two cases where the two obstacles are placed at symmetric locations with respect to the line y = x and the imaginary axis, respectively. The behavior of vortex trajectories in these two examples is the same. There is a vortex center located almost at zero as well as two saddles between each obstacle and the external boundary. Every two points of the four saddles are connected by a heteroclinic loop.
In the bottom row of Figure 14, we demonstrate the effect of symmetry breaking on vortex motion. As the patterns in Figure 14a,b are topologically equivalent, we break the vertical line symmetry only in the first example. First, after a small horizontal perturbation of 0.01 unit on the center of the lower obstacle, as shown in Figure 14e, the ( 2 , 1 ) homoclinic–heteroclinic pattern in the lower side of this obstacle is no longer present, and a new vortex center is born near the lower left corner of the same obstacle. The other center and three saddles are still present in the phase space with the same patterns. However, a similar horizontal perturbation of the upper obstacle does not cause the two lower saddles to disappear, but it does split their heteroclinic connections. Observe also that the upper critical points still occur (see Figure 14f). Looking at Figure 14g,h, we see the effect of breaking the line symmetry in the two examples shown in Figure 14c,d. Indeed, if we move one obstacle and keep the other at the same location, then the two saddle points near one obstacle get disconnected from the two saddles near the other.
In Figure 15, we show the effect of increasing the size of the square obstacles. Since the squares become too close, a new saddle point appears between them in Figure 15a,b. The dynamics of the point vortex motion is more remarkable in the case of Figure 15c. Two new saddles connected by a heteroclinic loop appear near the top right corner of the outer square. Moreover, on the other side, the ( 2 , 1 ) homoclinic–heteroclinic pattern in Figure 14b becomes an inner structure within a heteroclinic loop in a new ( 2 , 1 ) pattern.
Finally, when we consider a quintuply-connected polygonal domain with four square obstacles in Figure 16, we notice the same dynamical feature as in the corresponding case of triply-connected domains. It is worth mentioning that the presented method can also be efficiently used to compute vortex trajectories in multiply-connected polygonal domains of high connectivity, as illustrated in Figure 17.

6. Conclusions

We have numerically computed the Hamiltonian of a point vortex motion in simply- and multiply-connected bounded polygonal domains. We have first conformally mapped the polygonal domain onto a circular domain, and then applied the analytic formula presented in [11] for the Hamiltonian in multiply-connected circular domains.
For simply-connected domains, we have considered four polygonal shapes that are non convex. In these four examples, we have studied first the case when the system is topologically equivalent to the unit disk where the Hamiltonian has a unique critical point. Then, keeping the same general shape and changing the distance between the center point location and boundaries has yielded a lost of this uniqueness property. In all cases, the number of critical points satisfies the Morse Theorem; that is, the number of centers is equal to that of saddles plus one.
We have then examined the case of doubly-connected polygonal domains with four examples of an obstacle in a square. For concentric domains, we have noticed a pure heteroclinic structure of vortex trajectories around the standard and rotated squares obstacle, which was not the case for triangle and hexagon obstacles. Afterwards, we have concerned ourselves with the manner the vortex motion changes due to small perturbations of the obstacle location and its displacement toward the external boundary. Through various examples, the effect of symmetry breaking, and obstacle shape and location on vortex trajectories has been well demonstrated.
In triply-connected polygonal domains, we have looked at how vortex trajectories interact with the square obstacles and the outer square boundaries. We have observed that this interaction depends on the distance between these obstacles, and between each obstacle and the external boundary. The common dynamical feature at close boundaries interaction involves two saddles for side to side and one saddle for corner to side or corner to corner. As the separation between the two obstacles or between one obstacle and the outer boundaries increases, a vortex merging has been noticed. Another worth commenting dynamical feature is that heteroclinic connections between saddle points are quickly broken by lost of symmetry, which confirms their structural instability.
We have also demonstrated that our method can be efficiently used to study single point vortex motion in multiply-connected polygonal domains of high connectivity. Indeed, we have shown vortex trajectories in a square domain with four square obstacles, and a rectangle with fourteen rectangular obstacles.
Finally, as an outlook on future perspectives, we mention three directions. We have seen in this paper that heteroclinic loops is generally associated with the presence of symmetries in the model. Breaking such symmetries often breaks heteroclinicity and produces new homoclinic orbits. The formation of homoclinic orbits might be important for the mixing of passive tracers [24]. Future work could also focus on the case in which the Hamiltonian is not a log function but follows a power law, which are an important new frontier for applied mathematics and physics [25]. On the other hand, in this paper, we have considered the multiply-connected polygonal domains. With the help of the method presented in [19,20,26], our method can be extended to multiply-connected domains other than polygonal domains such as domains exterior to rectilinear slits.

Author Contributions

Conceptualization, E.M.K., M.M.S.N. and K.A.H.; methodology, E.M.K. and M.M.S.N.; software, E.M.K. and M.M.S.N.; validation, E.M.K., M.M.S.N., and K.A.H.; formal analysis, E.M.K., M.M.S.N., and K.A.H.; investigation, E.M.K. and K.A.H.; resources, E.M.K. and M.M.S.N.; writing—original draft preparation, E.M.K.; writing—review and editing, E.M.K., M.M.S.N., and K.A.H.; and visualization, E.M.K. and M.M.S.N. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by Qatar National Library.

Acknowledgments

The authors would like to thank the three anonymous reviewers for their careful reading and valuable suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Helmholtz, H. Über integrale der hydrodynamischen gleichungen, welche den Wirbelbewegungen entsprechen. J. Reine Angew. Math. 1858, 55, 25–55. [Google Scholar]
  2. Aref, H. Point vortex dynamics: A classical mathematics playground. J. Math. Phys. 2007, 48, 065401. [Google Scholar] [CrossRef] [Green Version]
  3. Boatto, S.; Crowdy, D.G. Point vortex dynamics. In Encyclopedia of Mathematical Physics; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  4. Newton, P. The N-Vortex Problem: Analytical Techniques; Springer: New York, NY, USA, 2002. [Google Scholar]
  5. Saffman, P. Vortex Dynamics; Cambridge University Press: Cambridge, UK, 1992. [Google Scholar]
  6. Flucher, M.; Gustafsson, B. Vortex motion in two-dimensional hydrodynamics (TRITA-MAT-97-MA-02). Royal Inst. Techn. Stockh. 1997. [Google Scholar]
  7. Kirchhoff, G. Vorlesungen über Mathematische Physik; Teubner: Leipzig, Germany, 1876. [Google Scholar]
  8. Routh, E.J. Some applications of conjugate functions. Proc. Lond. Mat. Soc. 1881, 12, 73–89. [Google Scholar] [CrossRef]
  9. Lin, C.C. On the motion of vortices in two dimensions. Proc. Nat. Acad. Sci. USA 1941, 27, 570–577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Johnson, E.R.; McDonald, N.R. The motion of a vortex near two circular cylinders. Proc. R. Soc. A 2004, 460, 939–954. [Google Scholar] [CrossRef]
  11. Crowdy, D.; Marshall, J. Analytical formulae for the Kirchhoff-Routh path function in multiply connected domains. Proc. R. Soc. A 2005, 461, 2477–2501. [Google Scholar] [CrossRef]
  12. Crowdy, D.; Marshall, J. The motion of a point vortex around multiple circular islands. Phys. Fluids 2005, 17, 056602. [Google Scholar] [CrossRef] [Green Version]
  13. Crowdy, D.; Marshall, J. The motion of a point vortex through gaps in walls. J. Fluid Mech. 2006, 551, 31–48. [Google Scholar] [CrossRef] [Green Version]
  14. Sakajo, T. Equation of motion for point vortices in multiply connected circular domains. Proc. R. Soc. A 2009, 465, 2589–2611. [Google Scholar] [CrossRef]
  15. Krantz, S.G. Geometric Function Theory: Explorations in Complex Analysis; Birkhäuser: Boston, MA, USA, 2006. [Google Scholar]
  16. Nasser, M.M.S. PlgCirMap: A MATLAB toolbox for computing the conformal mapping from polygonal multiply connected domains onto circular domains. arXiv 2019, arXiv:1911.01787. [Google Scholar]
  17. Crowdy, D. Solving Problems in Multiply Connected Domains. In Proceedings of the SIAM, CBMS-NSF Regional Conference Series in Applied Mathematics, Irvine, CA, USA, 18–22 June 2018. [Google Scholar]
  18. Crowdy, D.G.; Kropf, E.H.; Green, C.C.; Nasser, M.M.S. The Schottky-Klein prime function: A theoretical and computational tool for applications. IMA J. Appl. Math. 2016, 81, 589–628. [Google Scholar] [CrossRef] [Green Version]
  19. Nasser, M.M.S. Fast Computation of the Circular Map. Comput. Methods Funct. Theory 2015, 15, 187–223. [Google Scholar] [CrossRef] [Green Version]
  20. Nasser, M.M.S. Fast solution of boundary integral equations with the generalized Neumann kernel. Electron. Trans. Numer. Anal. 2015, 44, 189–229. [Google Scholar]
  21. Greengard, L.; Gimbutas, Z. FMMLIB2D: A MATLAB tOolbox for Fast Multipole Method in Two Dimensions. Version 1.2. Available online: http://www.cims.nyu.edu/cmcl/fmm2dlib/fmm2dlib.html (accessed on 1 January 2018).
  22. Gustafsson, B. On the motion of a vortex in two-dimensional flow of an ideal fluid in simply and multiply connected domains (RITA-MAT-1979-7). Royal Inst. Techn. Stockh. 1979. [Google Scholar]
  23. Gustafsson, B. On the convexity of a solution of Liouville’s equation. Duke Math. J. 1990, 60, 303–311. [Google Scholar] [CrossRef]
  24. Scott, R.K. Local and nonlocal advection of a passive scalar. Phys. Fluids 2006, 18, 116601. [Google Scholar] [CrossRef]
  25. Badin, G.; Barry, A.M. Collapse of generalized Euler and surface quasigeostrophic point vortices. Phys. Rev. E 2018, 98, 023110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Nasser, M.M.S.; Green, C.C. A fast numerical method for ideal fluid flow in domains with multiple stirrers. Nonlinearity 2018, 31, 815–837. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Vortex trajectories in a star-shaped domain. The vertices of the polygon are: 1 + 0.5 i , 1 + i , 0.5 + i , a i , 0.5 + i , 1 + i , 1 + 0.5 i , b , 1 0.5 i , 1 i , 0.5 i , a i , 0.5 i , 1 i , 1 0.5 i , and b.
Figure 1. Vortex trajectories in a star-shaped domain. The vertices of the polygon are: 1 + 0.5 i , 1 + i , 0.5 + i , a i , 0.5 + i , 1 + i , 1 + 0.5 i , b , 1 0.5 i , 1 i , 0.5 i , a i , 0.5 i , 1 i , 1 0.5 i , and b.
Symmetry 12 01175 g001
Figure 2. Vortex trajectories in a C-shaped polygon of vertices: 1 + i , 1 + i , 1 i , 1 i , 1 0.5 i , a 0.5 i , a + 0.5 i , 1 + 0.5 i .
Figure 2. Vortex trajectories in a C-shaped polygon of vertices: 1 + i , 1 + i , 1 i , 1 i , 1 0.5 i , a 0.5 i , a + 0.5 i , 1 + 0.5 i .
Symmetry 12 01175 g002
Figure 3. Vortex trajectories in a H-shaped polygon of vertices: 1 + i , 0.5 + i , 0.5 + a i , 0.5 + a i , 0.5 + i , 1 + i , 1 i , 0.5 i , 0.5 a i , 0.5 a i , 0.5 i , 1 i .
Figure 3. Vortex trajectories in a H-shaped polygon of vertices: 1 + i , 0.5 + i , 0.5 + a i , 0.5 + a i , 0.5 + i , 1 + i , 1 i , 0.5 i , 0.5 a i , 0.5 a i , 0.5 i , 1 i .
Symmetry 12 01175 g003
Figure 4. Vortex trajectories in a E-shaped polygon of vertices: 1 i , 1 0.6 i , 0.6 i , 0.2 i , 1 0.2 i , 1 + 0.2 i , 0.2 i , 0.6 i , 1 + 0.6 i , 1 + i , a + i , a i .
Figure 4. Vortex trajectories in a E-shaped polygon of vertices: 1 i , 1 0.6 i , 0.6 i , 0.2 i , 1 0.2 i , 1 + 0.2 i , 0.2 i , 0.6 i , 1 + 0.6 i , 1 + i , a + i , a i .
Symmetry 12 01175 g004
Figure 5. Vortex trajectories in doubly-connected circular domains. The outer circle is the unit circle | ζ | = 1 and the inner circle is | ζ p | = 0.1 .
Figure 5. Vortex trajectories in doubly-connected circular domains. The outer circle is the unit circle | ζ | = 1 and the inner circle is | ζ p | = 0.1 .
Symmetry 12 01175 g005
Figure 6. Vortex trajectories in doubly-connected square domains where the square obstacle has a center p and a horizontal radius of 0.1 .
Figure 6. Vortex trajectories in doubly-connected square domains where the square obstacle has a center p and a horizontal radius of 0.1 .
Symmetry 12 01175 g006
Figure 7. Vortex trajectories in doubly-connected polygonal domains where the vertices of the rotated square obstacle are on the circle of center p and radius 0.1 .
Figure 7. Vortex trajectories in doubly-connected polygonal domains where the vertices of the rotated square obstacle are on the circle of center p and radius 0.1 .
Symmetry 12 01175 g007
Figure 8. Vortex trajectories in doubly-connected square domains where the vertices of the triangle obstacle are: p + 0.1 e π i / 2 , p + 0.1 e π i / 6 , p + 0.1 e 5 π i / 6 .
Figure 8. Vortex trajectories in doubly-connected square domains where the vertices of the triangle obstacle are: p + 0.1 e π i / 2 , p + 0.1 e π i / 6 , p + 0.1 e 5 π i / 6 .
Symmetry 12 01175 g008
Figure 9. Critical vortex trajectories after a small displacement ( β = 0.05 ) of the triangle obstacle in different directions. The vertices of the triangle are: p + 0.1 e π i / 2 , p + 0.1 e π i / 6 , p + 0.1 e 5 π i / 6 .
Figure 9. Critical vortex trajectories after a small displacement ( β = 0.05 ) of the triangle obstacle in different directions. The vertices of the triangle are: p + 0.1 e π i / 2 , p + 0.1 e π i / 6 , p + 0.1 e 5 π i / 6 .
Symmetry 12 01175 g009
Figure 10. Critical vortex trajectories for bifurcation locations of the triangle obstacle. The vertices of the triangle are: p + 0.1 e π i / 2 , p + 0.1 e π i / 6 , p + 0.1 e 5 π i / 6 .
Figure 10. Critical vortex trajectories for bifurcation locations of the triangle obstacle. The vertices of the triangle are: p + 0.1 e π i / 2 , p + 0.1 e π i / 6 , p + 0.1 e 5 π i / 6 .
Symmetry 12 01175 g010
Figure 11. Vortex trajectories in doubly-connected polygonal domains where the vertices of the hexagonal obstacle are: p + 0.1 , p + 0.1 e π i / 3 , p + 0.1 e 2 π i / 3 , p 0.1 , p + 0.1 e 4 π i / 3 , p + 0.1 e 5 π i / 3 .
Figure 11. Vortex trajectories in doubly-connected polygonal domains where the vertices of the hexagonal obstacle are: p + 0.1 , p + 0.1 e π i / 3 , p + 0.1 e 2 π i / 3 , p 0.1 , p + 0.1 e 4 π i / 3 , p + 0.1 e 5 π i / 3 .
Symmetry 12 01175 g011
Figure 12. Vortex trajectories in triply-connected circular domains. The outer circle is the unit circle | ζ | = 1 and the inner circles are | ζ ± p | = 0.1 .
Figure 12. Vortex trajectories in triply-connected circular domains. The outer circle is the unit circle | ζ | = 1 and the inner circles are | ζ ± p | = 0.1 .
Symmetry 12 01175 g012
Figure 13. Vortex trajectories in triply-connected polygonal domains where the inner squares have side length equal to 0.2 and centers at p and p .
Figure 13. Vortex trajectories in triply-connected polygonal domains where the inner squares have side length equal to 0.2 and centers at p and p .
Symmetry 12 01175 g013
Figure 14. Vortex trajectories in triply-connected polygonal domains where the inner squares have side length equal to 0.2 , and centers at p and q, respectively.
Figure 14. Vortex trajectories in triply-connected polygonal domains where the inner squares have side length equal to 0.2 , and centers at p and q, respectively.
Symmetry 12 01175 g014
Figure 15. Vortex trajectories in triply-connected polygonal domains where the large and small inner squares have side length equal to 0.5 and 0.4 , and centers at p and q, respectively.
Figure 15. Vortex trajectories in triply-connected polygonal domains where the large and small inner squares have side length equal to 0.5 and 0.4 , and centers at p and q, respectively.
Symmetry 12 01175 g015
Figure 16. Vortex trajectories in quintuply-connected polygonal domains where the inner squares have side length equal to 0.2 and centers at ± p and ± q .
Figure 16. Vortex trajectories in quintuply-connected polygonal domains where the inner squares have side length equal to 0.2 and centers at ± p and ± q .
Symmetry 12 01175 g016
Figure 17. Vortex trajectories in a multiply-connected polygonal domain with 15 obstacles.
Figure 17. Vortex trajectories in a multiply-connected polygonal domain with 15 obstacles.
Symmetry 12 01175 g017
Table 1. Number of saddle points and their connecting critical loops in doubly-connected concentric polygonal domains with different obstacles shapes.
Table 1. Number of saddle points and their connecting critical loops in doubly-connected concentric polygonal domains with different obstacles shapes.
Obstacle ShapeSaddleHomoclinicHeteroclinic
Square808
Rotated square404
Triangle341
Hexagon442

Share and Cite

MDPI and ACS Style

Kalmoun, E.M.; Nasser, M.M.S.; Hazaa, K.A. The Motion of a Point Vortex in Multiply-Connected Polygonal Domains. Symmetry 2020, 12, 1175. https://doi.org/10.3390/sym12071175

AMA Style

Kalmoun EM, Nasser MMS, Hazaa KA. The Motion of a Point Vortex in Multiply-Connected Polygonal Domains. Symmetry. 2020; 12(7):1175. https://doi.org/10.3390/sym12071175

Chicago/Turabian Style

Kalmoun, El Mostafa, Mohamed M. S. Nasser, and Khalifa A. Hazaa. 2020. "The Motion of a Point Vortex in Multiply-Connected Polygonal Domains" Symmetry 12, no. 7: 1175. https://doi.org/10.3390/sym12071175

APA Style

Kalmoun, E. M., Nasser, M. M. S., & Hazaa, K. A. (2020). The Motion of a Point Vortex in Multiply-Connected Polygonal Domains. Symmetry, 12(7), 1175. https://doi.org/10.3390/sym12071175

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