Advances in Mathematical Physics

Volume 2015, Article ID 613683, 9 pages

http://dx.doi.org/10.1155/2015/613683

## The Ellipsoidal Vortex: A Novel Approach to Geophysical Turbulence

ISMAR-CNR, Arsenale-Tesa 104, Castello 2737/F, 30122 Venice, Italy

Received 8 October 2014; Accepted 14 January 2015

Academic Editor: Touvia Miloh

Copyright © 2015 William J. McKiver. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We review the development of the ellipsoidal vortex model within the field of geophysical fluid dynamics. This vortex model is built on the classical potential theory of ellipsoids and applies to large-scale fluid flows, such as those found in the atmosphere and oceans, where the dynamics are strongly affected by the Earth's rotation. In this large-scale limit the governing equations reduce to the quasi-geostrophic system, where all the dynamics depends on a single scalar field, the potential vorticity, which is a dynamical marker for vortices. The solution of this system is achieved by the inversion of a Poisson equation, that in the case of an ellipsoidal vortex can be solved exactly. From this ellipsoidal solution equilibria have been determined and their stability properties have been studied. Many studies have shown that this ellipsoidal vortex model, while being conceptually simple, is an extremely powerful tool in eliciting some of the fundamental characteristics of turbulent geophysical flows.

#### 1. Introduction

In turbulent flow there exist structures whose collective interactions are at the very heart of the dynamics of turbulence—vortices. The vortex, a word signifying a type of swirling or spinning fluid, has become a familiar concept in human society, being associated with extreme and highly energetic weather phenomena such as typhoons, tornados, and hurricanes. The relative infrequency of these most dramatic manifestations of vortices could lead one to think that the concept of the vortex itself is limited only to these extreme events, but, on the contrary, vortices are present in almost all fluid phenomena, from the smallest scale flow of water in pipes to the largest scales within the atmosphere and the oceans. Both observations [1, 2] and numerical simulations [3–5] have revealed the widespread appearance and dominance of vortices in turbulent geophysical flows. The prolific nature of vortices has led to their study becoming a field in its own right as a means to better understand and resolve the great unsolved problem of classical physics—turbulence.

While vorticity has a clear mathematical definition, being defined as a vector field which is the curl of the velocity flow field (), the definition of a vortex itself is far more elusive. This is due to the fact that a vortex does not have a particular shape. Ideas of tornadoes and hurricanes usually conjure up an image of circular or cylindrical type shapes for the vortex, but real vortices are not constrained to these perfect forms. However, when approaching them mathematically it is useful to apply a shape that is amenable to simple mathematical definition. One of the earliest incarnations of such an approach was that of Kirchhoff [6], who studied an isolated two-dimensional patch of uniform vorticity bounded within an ellipse. He showed that this is an exact analytical solution to the Euler equations which rotates with a constant angular velocity. Further work by Love [7] showed that there existed stable equilibria when this vortex is subjected to disturbances.

The extension of these early works to the three-dimensional case of an ellipsoidal vortex had to wait almost a century for the works of Zhmur and Pankratov [8, 9] and Meacham [10] who extended this approach to the quasi-geostrophic system [11], a set of equations relevant to large-scale geophysical fluid dynamics. In fact, the ellipsoidal form has a long history stretching back (at least) to the seminal works of MacClaurin [12] and Laplace [13] on gravitating bodies (see [14, 15] for a review of this subject). As the ellipsoidal form has an exact and simple mathematical definition that yet goes beyond the absolute symmetry of the sphere, it is a powerful tool in determining shape, structure, and orientation of “bodies” concerned in many different fields of application (see [16] for a review of the many applications). In the case of geophysical flows, it gives a means to determine the characteristics of vortices (such as typical height-to-width ratios, orientations, and rotation rates) that are the most robust and hence the most likely to occur in real turbulent flows.

Here we will review the recent history of the development of the ellipsoidal vortex in quasi-geostrophic dynamics. The beginning is in Section 2 where we give a brief overview of the governing system of equations, the quasi-geostrophic model [11]. Following this we review the ellipsoidal vortex model, starting in Section 3 with the case of the isolated ellipsoid as originally derived by Meacham [10]. In Section 4 we consider the case of an ellipsoidal vortex immersed in a background shear flow, following the matrix formulation approach of McKiver and Dritschel [17]. In Section 5 we discuss the results of studies of the equilibria and stability properties of this vortex system and what insights they provide for the full solution to the equations in turbulence simulations. In Section 6 we discuss the generalization of the single vortex ellipsoidal theory to the multiple vortex case, in particular outlining the Hamiltonian method of Dritschel et al. [18]. Finally we conclude in Section 7, with an overview of the consequences the ellipsoidal model has for the theory of turbulence and discuss future avenues of research in this field.

#### 2. The Quasi-Geostrophic System

Here we summarize the quasi-geostrophic model (for a detailed derivation see [19]). Quasi-geostrophic (hereafter QG) theory is one of the simplest geophysical fluid dynamical models that incorporates two crucial features underpinning large-scale fluid dynamics: the effect of the Earth’s rotation and the effects of strong density stratification (lighter fluid lying on top of denser fluid). As a consequence of these features certain terms in the full equations of motion (Navier-Stokes equations) tend to dominate over others. Strong rotation leads to what is termed “geostrophic balance,” which is a balance between the horizontal pressure gradient and the term associated with the Earth’s rotation (the Coriolis term): where is the local vertical component of the planetary rotation rate, is a vertical unit vector, is the horizontal velocity field, is pressure, and is the basic-state density distribution ( is nearly constant in the oceans and exponential, proportional to with km, in the atmosphere). Strong stratification leads to “hydrostatic balance,” which is a balance between the vertical pressure gradient and buoyancy: where is the acceleration due to gravity. The QG system, which incorporates both hydrostatic and geostrophic balance at leading order, is obtained from the full “primitive” equations by an asymptotic expansion in three small nondimensional parameters: the height-width aspect ratio of characteristic motions (constrained to be small due to the shallow geometry in the oceans and atmosphere); the Rossby number ; and the Froude number , where is a characteristic relative horizontal fluid speed, is the horizontal scale of the flow, is the average fluid depth, and is the buoyancy frequency (the oscillation frequency that a small fluid volume would exhibit if it were displaced by a small vertical distance). The QG model is obtained assuming that and . Also in this model we consider the limit of infinite Reynolds number, ; that is, the effect of viscous dissipation, , is negligible (in the real atmosphere, ). As done in many previous works, one can consider the constant-coefficient (“oceanic”) system when , , and are all constant. For an inviscid, adiabatic fluid, the QG equations may be expressed in terms of the conservation of a single material scalar, the “potential vorticity,” ; that is, The potential vorticity (hereafter PV) is proportional to the component of vorticity perpendicular to stratification surfaces (surfaces of constant density ). Thus vortices are identified as contiguous regions of PV. In the QG model, the vertical velocity to this order of approximation, so is advected in a layerwise manner. However the horizontal velocity depends on all three spatial coordinates, and this is recovered from via the inversion of for the stream function , followed by the incompressibility relations Note that in this form the vertical coordinate is rescaled by the factor . The three-dimensional distribution of PV generates the velocity field through (4) and (5), which carries the PV, conservatively, to the next instant of time. While (4) is an isotropic relation between and , the advection remains anisotropic as there is no vertical motion. The QG equations are widely used in studies of atmospheric and oceanic phenomena [20–23], in part because they are far easier to deal with both theoretically and numerically than the complete equations. But, also, the QG equations capture many qualitative features of atmospheric and oceanic dynamics.

#### 3. Potential Theory of an Isolated Ellipsoidal Vortex

Now we consider the case of an isolated ellipsoidal vortex, with uniform potential vorticity, as first considered in the works of Zhmur and Pankratov [8, 9] (see also [24]) and independently by Meacham [10]. We will review the full solution to the QG ellipsoid as derived by Meacham [10]. We consider the ellipsoid of PV having semiaxes lengths , , and , where . Without loss of generality we assume the vortex axes are aligned along the coordinate axis. Thus the boundary of the ellipsoid is defined by the equation Inside the ellipsoid we have from (4) where is the uniform potential vorticity. A general solution to this equation can be written as where , , , and are unknown coefficients to be determined. Substituting this into (7) gives the condition Outside the ellipsoid the potential vorticity is zero and the external potential satisfies With such a problem the natural approach is to move from Cartesian coordinates to an ellipsoidal coordinate system using the relations [25] where the variables , , and represent the roots to the cubic equation where . The -surfaces are ellipsoids, the -surfaces are hyperboloids of one sheet, and the -surfaces are hyperboloids of two sheets [16]. The largest root, , is analogous to the radial variable in spherical coordinates with each representing a family of concentric ellipsoids, with being the value on the surface of the original ellipsoid with semiaxes , , and . The boundary conditions on the surface of the ellipsoid are the continuity of the normal and tangential velocity components, which can be expressed in ellipsoidal coordinates as Substituting (11a), (11b), and (11c) into the expression for the internal field (8) we obtain A general solution for the external flow field can be written as an expansion in ellipsoidal harmonics [16, 25, 26] where denote Lame functions of the first kind, while are Lame functions of the second kind defined by The dependence on the coordinates and for the inner field in (14) implies that the exterior solution needed in order to match the interior solution at the boundary has the form where , , and are unknown coefficients and the second-order Lame functions of the first kind are where Substituting these Lame functions into the external field gives us Substituting (14) and (20) into the boundary conditions equations (13a) and (13b) provides 6 conditions, since there are terms proportional to and as well as the constant terms. These 6 conditions along with the condition given by (9) provide 7 equations in the 7 unknown coefficients and so can be solved as was shown by Meacham [10]. Solving these simultaneous equations we can write the interior stream function as where , , and are constant coefficients given by where is the vortex “strength” and and are elliptic integrals of the first and second kind, respectively, defined as This solution was originally determined by Laplace [13] for the case of a gravitating ellipsoid of uniform mass [14]. The exterior stream function is where , , and and now we have While the interior solution has quadratic dependence on spatial coordinates, the exterior field has a more complicated dependence, due to the coefficients being functions of the variable . However, the flow field that acts on the vortex itself is the interior solution, and so one can compute the flow field acting on the vortex itself from (5). This gives rise to a self-induced velocity field that is linear in spatial coordinates.

#### 4. Evolution of a Vortex in a Background Shear Flow

We now consider the case of an ellipsoidal vortex immersed in a background shear flow. This background shear flow idealizes the effects of other vortices and so is the first step towards understanding interactions between vortices. This case was first considered by Meacham et al. [27] for a number of specific types of background flow and later generalized by Meacham et al. [28] where it was recast within a Hamiltonian framework. Hashimoto et al. [29] extended the analytical solutions of Meacham [10] for the isolated vortex by solving for exact equilibrium solutions for the case of the ellipsoid in an external straining flow. McKiver and Dritschel [17] introduced a reformulation of the system that greatly simplified the equations to a simple matrix problem. In this section we review their derivation of this matrix formulation. This approach is built on the simple fact that the stream function for the ellipsoid has quadratic dependence on spatial coordinates, and so the velocity field has a linear dependence. Naturally any linear flow field imposed on the ellipsoidal vortex, while possibly altering the size (, , and ) and orientation of the vortex, preserves the ellipsoidal form. Thus we consider a general velocity field of the form
where is a matrix. Within this flow, we place a single ellipsoidal vortex at the origin and require only that its velocity field has the same form as (26); that is, . The total velocity field felt by the ellipsoid is
where , and this depends only on time . We refer to as the* flow matrix* subsequently.

The ellipsoid is specified by its axis half lengths and the unit vectors , **,** and directed along these axes. This geometric information may be encapsulated in a single symmetric matrix defined by
where the superscript denotes transpose, is an orthonormal rotation matrix () defined by
and is a diagonal matrix whose terms are the inverse squared semiaxes lengths:
The equation describing the surface of the ellipsoid is conveniently written as
By taking a time derivative of this equation and using , as well as the implied relation , we obtain
Since this must be true for all points on the boundary of the ellipsoid, it follows that
which is simply the equation for the evolution of in an arbitrary linear background flow. Note that remains a symmetric matrix for all time since is automatically symmetric if is.

It proves convenient to rewrite (33) to evolve not but its inverse . Using the fact that the product is the identity matrix, it follows that giving us the evolution equation (as originally derived in [17]) where In order to determine the evolution of the ellipsoid it is necessary to be able to calculate the semiaxes lengths, , , and , and the orientation vectors , , and of the ellipsoid from the matrix . But (36), right multiplied by , together with (37) implies Hence the semiaxes lengths and the orientation vectors of the ellipsoid can be found directly from the matrix by solving a simple eigenvalue problem. Finally one must specify the flow matrix . The self-induced part is derived from the stream function solution to the inner potential problem, . This solution can be written more generally for an arbitrarily orientated ellipsoid as where and is the diagonal matrix formed from the solution coefficients: From (5) it follows that where is defined as

For the background flow matrix one can consider the effect of a single distant vortex, of strength , centered at . Let be the distance between the two vortices; then, to leading order in the vortices rotate about each other at a rate At this order of approximation the vortices appear as points, with no shape or internal structure. If we adopt a frame of reference rotating with the -axis passing through the joint centre of the two vortices, so that their positions and are fixed, and assuming the original vortex is located at the origin, , then in the vicinity of the origin the background stream function takes the form (Taylor expanding) where notably the linear terms are absent and the matrix is given by The corresponding velocity field is given by , where as before the flow matrix is , leading to where , hereinafter referred to as “the strain rate,” and is a parameter depending only on the ratio of the vortex strengths. It can be simplified further by choosing either or to be zero, leaving only three nonzero components. For example, if , where , , and the range of the angle is restricted to because of symmetry. This form for the background flow matrix was originally considered by Meacham et al. [27] but using a different notation.

The complete model now consists of the evolution equation (35) for the matrix together with expressions (42) and (50) for the flow matrices and . Note that in QG flow—this follows because there is no vertical velocity. The element is the squared half-height of the vortex. Note also that the determinant is proportional to the squared vortex volume and is also invariant (because QG flow is incompressible).

#### 5. Equilibria and Stability

Given the equations governing the evolution of an ellipsoidal vortex in a background shear flow one can then determine equilibrium vortices from (35), where they satisfy This was first exploited by Reinaud et al. [30] using an iterative linear method to determine equilibria. These equilibria have important consequences for vortex interactions. As a vortex evolves over its lifetime it may come into close proximity with other vortices. When this happens the strain exerted on one vortex by the others increases, as the magnitude of this strain scales with the inverse cube of the separation distance between vortices (see (48)). At some “critical” value of strain the vortex will begin to deform significantly and will undergo a strong interaction. If the vortices are close enough they may partially or fully coalesce or “merge” into one, forming a larger vortex along with filamentary debris. Often the onset of such strong interactions occurs at the point beyond which equilibria can no longer exist for strains greater than this critical value, or beyond which equilibria are unstable to disturbances. Thus the stability analysis of a vortex in equilibrium with a background shear is an important indicator of when strong interactions may occur. Next we give a brief summary of the findings of a number of studies on this topic.

The stability of a freely rotating ellipsoid of potential vorticity was first considered by Meacham [10] who found instabilities over a large range of the parameter space characterising the ellipsoid. However, a more recent study by Dritschel et al. [31] found that one of the unstable modes Meacham encountered is in fact stable, and overall their results indicate that the freely rotating ellipsoid is widely stable.

Miyazaki et al. [32] studied the stability of a spheroid of uniform PV tilted by some inclination angle from the vertical axis. They found that highly prolate spheroids are unstable if the inclination angle is large, while oblate spheroids are unstable even if the inclination angle is very small. In both cases the instability is nonellipsoidal; that is, it destroys the ellipsoidal form of the vortex.

Hashimoto et al. [29] derived the equations for the linear stability of an ellipsoid in a 2D strain field by expanding the disturbing surface in Lame functions. Using these they studied the cases of a pure strain field (when there is no background rotation) and a simple shear flow (when the strain rate is equal to the background rotation). They found that, in a pure strain field, highly elongated ellipsoids are unstable to modes whose order is greater than , that is, to nonellipsoidal modes. In a simple shear flow they found that a highly elongated ellipsoid whose major axis is perpendicular to the flow direction is unstable, whereas any ellipsoidal vortex seems to be stable if the major axis is parallel to the flow direction.

Meacham et al. [27] derived equations for the evolution of an ellipsoid of uniform PV in a background flow with both horizontal strain and vertical shear. They examined the stability for some special background flows, either pure horizontal strain or pure vertical shear (in the absence of horizontal strain). A full linear stability analysis for the case of a general linear background flow (50) was conducted by McKiver and Dritschel [33]. They characterized the system by four parameters: the height-to-width aspect ratio of the vortex , the strength ratio , the angle , and the strain rate . Over nearly the entire parameter space examined they found that the ellipsoidal steady states are stable, and where instability does occur the dominant modes encountered are ellipsoidal. Thus the ellipsoidal vortex in a linear background flow is a fairly robust model. They also found that for vortices with aspect ratios greater than the most unstable vortices (i.e., the ones which destabilize at the smallest strain values ) are found in the ranges and , for opposite- and like-signed vortices, respectively. Because the strain is proportional to the inverse cube of the separation distances between two vortices, this implies that vortices which are vertically offset by these values destabilize from the greatest separation distances; that is, they are least resilient. This agrees with the findings of Reinaud and Dritschel [34] in which they examined the merger of two identical uniform PV vortices which are offset both horizontally and vertically and found that vortices which are moderately offset in the vertical merge from a greater separation distance than vortices that are not offset vertically. In McKiver and Dritschel [33], for opposite-signed vortices, they found that the most stable vortices have an aspect ratio between and , whereas for the like-signed vortices the most stable vortices have an aspect ratio near . This agrees with the findings in a study by Reinaud et al. [30] where they performed simulations of a turbulent flow comprising hundreds of QG vortices in order to determine the typical shape of vortices in QG turbulence. Remarkably their study also found that typically vortices are slightly oblate having a height-to-width aspect ratio of .

#### 6. Multiple Ellipsoidal Vortices

While the single vortex case provides great insight into vortex behaviour, in order to understand some of the dynamics that precipitate strong interactions and vortex mergers, one must consider a multivortex scenario. There are a number of works that have generalized the single ellipsoidal vortex solution to a multivortex scenario [35–38]. However, because of the form of the external stream function (24) generated by a vortex, in general it can induce nonellipsoidal disturbances on external vortices. This means that assumptions must be made in order to retain the ellipsoidal form. Here we will briefly summarize an accurate approach introduced by Dritschel et al. [18]. This approach, referred to as ELM, represents vortices as PV ellipsoids and filters higher-order nonellipsoidal deformations simply by preserving only the linear part of the exterior field. Following the matrix formulation in Section 4 each vortex can be specified by the value of its potential vorticity, , the symmetric matrix , and the location of the center of the vortex , in terms of which the surface of each ellipsoid can be written as where , where is the number of vortices. If we consider that the velocity field for a particular ellipsoid has a linear dependence on spatial coordinates, that is, where is an arbitrary velocity field and is the flow matrix, then this velocity field is the most general form which exactly preserves the ellipsoidal form of the vortex. The governing equation for ELM [18] is a finite Hamiltonian system where the equations of motion for the th ellipsoid are where is the vortex strength, is the Hamiltonian of the system defined as the total energy of the system divided by , and the flow matrix is decomposed into two parts, , the self-induced part and the part induced by all other ellipsoids . For efficiency of the method this background flow field is modelled for each ellipsoid as a finite sum of singular point vortices. The position and strength of the singularities are obtained such that the approximate stream function induced by them matches exactly the function at a given order of accuracy in , being the distance from to the point of evaluation. Dritschel et al. [18] find that high accuracy () can typically be achieved by using 7 point singularities for each ellipsoid. Comparisons of the ELM with full numerical simulations (Figure 1) of the QG vortices show remarkable agreement (see http://www-vortex.mcs.st-and.ac.uk/~dgd/ELM/ to see more comparisons).