Advection of passive scalars induced by a bay-trapped nonstationary vortex

A simple model of fluid particle advection induced by the interaction of a point vortex and incident plane flow occurring near a curved boundary is analyzed. The use of the curved boundary in this case is aimed at mimicking the geometry of an isolated bay of a circular shape. An introduction of such a boundary to the model results in the appearance of retention zones, where the vortex can be permanently trapped being either stationary or periodically oscillating. When stationary, it induces a steady velocity field that in turn ensures regular advection of nearby fluid particles. When the vortex oscillates periodically, the induced velocity field turns unsteady leading to the manifestation of chaotic advection of fluid particles. We show that the size of the fluid region engaged into chaotic advection increases almost monotonically with the increased magnitude of the vortex oscillations provided the magnitude remains relatively small. The monotonicity is accounted for the fact that the frequency of the vortex oscillations incommensurable with the proper frequency of fluid particle rotations in the steady state. Another point of interest is that it is demonstrated that bounded regions, in which the vortex may be trapped, can appear even at a significant distance from the bay. Making use of a Lagrangian indicator, examples of fluid particle advection induced by the periodic motion of the vortex inside the bay are adduced.     

利用GS流场重构方法研究磁尾等离子体片涡流

2000年9月30日Geotail卫星分别于17∶54∶36~18∶09∶00UT和18∶59∶00~19∶30∶00UT在磁尾晨侧等离子体片内(n≈0.4 cm^-3,T≈6 keV)观测到等离子体涡流事件.本文采用Grad-Shafranov (GS)流场重构技术再现了这些涡流的二维速度场、离子数密度和离子温度的分布图像.结果显示:从地心太阳磁层坐标系(GSM)赤道面上面看, 涡流的尺度约为5000 km×1400 km , 朝地球的运动速度约为15~25 km/s.所有5个涡流的旋转方向都为顺时针方向,旋转周期约为6~11 min.相邻涡流的相互作用导致它们之间的磁场强度增强.考察观测数据发现,涡流内不仅包含等离子体片热等离子体成分,也包含较大通量的类似源自磁鞘的冷等离子体成分(T<1 keV).这与观测到涡流等离子体的平均温度(T≈4 keV)较磁尾等离子体片等离子体的典型温度(T≈6 keV)明显偏低的事实是一致的.不仅如此,离子数密度和温度在结构内的分布也不均匀,数密度在涡流内部偏离中心的位置比较低而在每个涡流的边缘位置比较高,温度的分布大体上与密度相反.分析认为观测到的磁尾等离子体涡流事件可能由发生在低纬边界层的Kelvin-Helmholtz不稳定性引起,涡流结构内的冷等离子体可能来自磁层顶外部的磁鞘.     

Hydrodynamical Modeling Of Oceanic Vortices

Mesoscale coherent vortices are numerous in the ocean.Though they possess various structures in temperature and salinity,they are all long-lived, fairly intense and mostly circular. Thephysical variable which best describes the rotation and the density anomaly associated with coherent vortices is potential vorticity. It is diagnostically related to velocity and pressure, when the vortex is stationary. Stationary vortices can be monopolar (circular or elliptical) or multipolar; their stability analysis shows thattransitions between the various stationary shapes are possible when they become unstable. But stable vortices can also undergo unsteady evolutions when perturbed by environmental effects, likelarge-scale shear or strain fields, -effect or topography. Changes in vortex shapes can also result from vortex interactions. such as the pairing, merger or vertical alignment of two vortices, which depend on their relative polarities and depths. Such interactions transfer energy and enstrophy between scales, and are essential in two-dimensional and in geostrophic turbulence. Finally, in relation with the observations, we describe a few mechanisms of vortex generation.     

Successive formation of planetary lenses in an intermediate layer

Abstract

We study the formation of lenses of the ocean's intermediate water using a 2.5-layerβ-plane primitive equation model with localized injection of water mass. For the injecting rate of 1.0 Sv, we have observed that strong vortices are shed regularly. These vortices propagate westward much faster than the second baroclinic long Rossby wave. They are totally isolated from each other and show strong baroclinicity as well. Moreover, they remain stable over a sufficiently long period of time. Regular formation of such strong vortices in the intermediate layer has not been reported previously. The translation speed is explained using the Euler's momentum integral theorem for the nonlinear baroclinic vortex on the β-plane. We have demonstrated that coupling between the primary motion in the intermediate layer and the secondary motion in the upper layer with a meridional shift is crucial to the fast westward translation of the intense vortices. A simple dispersion formula relating the zonal translation speed with the vortex radius is also derived under the assumption of quasi-geostrophy. It has turned out that the analytical relation explains the numerical results surprisingly well despite the limitation of its derivation.     

Wave-induced steady streaming,mass transport and net sediment transport in rough turbulent ocean bottom boundary layers

The interaction between two important mechanisms which causes streaming has been investigated by numerical simulations of the seabed boundary layer beneath both sinusoidal waves and Stokes second order waves, as well as horizontally uniform bottom boundary layers with asymmetric forcing. These two mechanisms are streaming caused by turbulence asymmetry in successive wave half-cycles (beneath asymmetric forcing), and streaming caused by the presence of a vertical wave velocity within the seabed boundary layer as earlier explained by Longuet-Higgins. The effect of wave asymmetry, wave length to water depth ratio, and bottom roughness have been investigated for realistic physical situations. The streaming induced sediment dynamics near the ocean bottom has been investigated; both the resulting suspended load and bedload are presented. Finally, the mass transport (wave-averaged Lagrangian velocity) has been studied for a range of wave conditions. The streaming velocities beneath sinusoidal waves (Longuet-Higgins streaming) is always in the direction of wave propagation, while the streaming velocities in horizontally uniform boundary layers with asymmetric forcing are always negative. Thus the effect of asymmetry in second order Stokes waves is either to reduce the streaming velocity in the direction of wave propagation, or, for long waves relative to the water depth, to induce a streaming velocity against the direction of wave propagation. It appears that the Longuet-Higgins streaming decreases as the wave length increases for a given water depth, and the effect of wave asymmetry can dominate, leading to a steady streaming against the wave propagation. Furthermore, the asymmetry of second order Stokes waves reduces the mass transport (wave-averaged Lagrangian velocity) as compared with sinusoidal waves. The boundary layer streaming leads to a wave-averaged transport of suspended sediments and bedload in the direction of wave propagation.     

Point vortex dynamics for coupled surface/interior QG and propagating heton clusters in models for ocean convection

Abstract

The dynamic behavior of baroclinic point vortices in two-layer quasi-geostrophic flow provides a compact model for studying the transport of heat in a variety of geophysical flows including recent heton models for open ocean convection as a response to spatially localized intense surface cooling. In such heton models, the exchange of heat with the region external to the compact cooling region reaches a statistical equilibrium through the propagation of tilted heton clusters. Such tilted heton clusters are aggregates of cyclonic vortices in the upper layer and anti-cyclonic vortices in the lower layer which collectively propagate almost as an elementary tilted heton pair even though the individual vortices undergo shifts in their relative locations. One main result in this paper is a mathematical theorem demonstrating the existence of large families of long-lived propagating heton clusters for the two-layer model in a fashion compatible to a remarkable degree with the earlier numerical simulations. Two-layer quasi-geostrophic flow is an idealization of coupled surface/interior quasi-geostrophic flow. The second family of results in this paper involves the systematic development of Hamiltonian point vortex dynamics for coupled surface/interior QG with an emphasis on propagating solutions that transport heat. These are novel vortex systems of mixed species where surface heat particles interact with quasi-geostrophic point vortices. The variety of elementary two-vortex exact solutions that transport heat include two surface heat particles of opposite strength, tilted pairs of a surface heat particle coupled to an interior vortex of opposite strength and two interior tilted vortices of opposite strength at different depths. The propagation speeds of the tilted elementary hetons in the coupled surface/interior QG model are compared and contrasted with those in the simpler two-layer heton models. Finally, mathematical theorems are presented for the existence of large families of propagating long-lived tilted heton clusters for point vortex solutions in coupled surface/interior QG flow.     

Calculating the macrodispersion coefficient of the ensemble averaged solute transport equation in the discrete domain

Recognizing the spatial heterogeneity of hydraulic parameters, many researchers have studied the solute transport by both groundwater and channel flow in a stochastic framework. One of the methodologies used to up‐scale the stochastic solute transport equation, from a point‐location scale to a grid scale, is the cumulant expansion method combined with the calculus for the time‐ordered exponential and the calculus for the Lie operator. When the point‐location scale transport equation is scaled up to the grid scale, using the cumulant expansion method, a new dispersion coefficient emerges in the dispersive term of the solute transport equation in addition to the molecular dispersion coefficient. This velocity driven dispersion is called ‘macrodispersion’. The macrodispersion coefficient is the integral function of the time‐ordered covariance of the random velocity field. The integral is calculated over a Lagrangian trajectory of the flow. The Lagrangian trajectory depends on the following: (i) the spatial origin of the particle; (ii) the time when the macrodispersion is calculated; and (iii) the mean velocity field along the trajectory itself. The Lagrangian trajectory is a recursive function of time because the location of the particle along the trajectory at a particular time depends on the location of the particle at the previous time. This recursive functional form of the Lagrangian trajectory makes the calculation of the macrodispersion coefficient difficult. Especially for the unsteady, spatially non‐stationary, non‐uniform flow field, the macrodispersion coefficient is a highly complex expression and, so far, calculated using numerical methods in the discrete domains. Here, an analytical method was introduced to calculate the macrodispersion coefficient in the discrete domain for the unsteady and steady, spatially non‐stationary flow cases accurately and efficiently. This study can fill the gap between the theory of the ensemble averaged solute transport model and its numerical implementations. Copyright © 2011 John Wiley & Sons, Ltd.     

The formation of new quasi-stationary vortex patterns from the interaction of two identical vortices in a rotating fluid

Within the framework of the quasi-geostrophic approximation, the interactions of two identical initially circular vortex patches are studied using the contour dynamics/surgery method. The cases of barotropic vortices and of vortices in the upper layer of a two-layer fluid are considered. Diagrams showing the end states of vortex interactions and, in particular, the new regime of vortex triplet formation are constructed for a wide range of external parameters. This paper shows that, in the nonlinear evolution of two such (like-signed) vortices, the filaments and vorticity fragments surrounding the merged vortex often collapse into satellite vortices. Therefore, the conditions for the formation and the quasi-steady motions of a new type of triplet-shaped vortex structure are obtained.     

The interaction of two co-rotating quasi-geostrophic vortices in the vicinity of a surface buoyancy filament

In this paper, we investigate the interaction between two like-signed quasi-geostrophic uniform potential vorticity internal vortices in the vicinity of a surface buoyancy anomaly filament in a three dimensional, stably stratified and rapidly rotating fluid. The surface buoyancy distribution locally modifies the pressure fields and generates a shear flow. We start the study by first considering the effects of a uniform linear horizontal shear on the binary vortex interaction. We confirm that a cooperative shear facilitates the merger of a pair of vortices while an adverse shear has the opposite effect. We next investigate the binary vortex interaction in the vicinity of the surface buoyancy filament explicitly. Here, not only the filament generates a shear flow, but it also responds dynamically to the forcing by the vortex pair. The filament destabilises and forms buoyancy billows at the surface. These billows interact with the internal vortices. In particular, a surface billow may pair with one of the internal vortices. In such cases, the like-signed internal vortex pair may separate if they are initially moderately distant from each other.     

Some interactions between small numbers of baroclinic,geostrophic vortices

Abstract

Various interactions between small numbers (two and four) of baroclinic, geostrophic point vortices in a two-layer system are studied with attention to the qualitative changes in behavior which occur as size of the deformation radius is varied.

A particularly interesting interaction, which illustrates the richness of baroclinic vortex dynamics, is a collision between two hetons. (A heton is a vortex pair in which the constituent vortices have opposite signs and are in opposite layers. The “breadth” of a heton is the distance between its constituent vortices. A translating heton transports heat.) When two hetons, which initially have different breadths, collide, the result is either an exchange of partners, or a “slip-through” collision in which the initial structures are preserved. It is shown here that the outcome is always an exchange, provided the deformation radius is sufficiently small. This strongly contrasts with a collision between pairs of classical, one-layer vortices in which no exchange occurs if the initial ratio of the breadths is sufficiently extreme.

Finally the transport of passive fluid by a translating baroclinic pair is investigated. A pair of vortices in the top layer transports no lower layer fluid if the distance between the vortices is less than 1.72 deformation radii. By contrast, the size of the region trapped by a heton increases without bound as the spacing between the vortices increases.     

基于非结构化网格的边界积分方程方法的断层自发破裂模拟

地震自发破裂模拟是震源动力学研究的重要内容,了解复杂的断层动力学破裂过程对深入认识震源特征和解释运动学反演结果具有重要意义.基于边界积分方程方法的破裂模拟已经被广泛使用,大多采用的是平面断层模型的结构化网格划分.由于实际的断层往往具有较为复杂的几何特征,为了更为灵活地刻画断层几何复杂性,我们建立断层模型的三角形网格离散方案,通过精确的解析解形式来计算断层各个单元之间的应力格林函数,联立滑动弱化摩擦准则和非奇异边界积分方程,对断层的自发破裂过程进行了模拟.在简单的平面断层模型下,将计算结果与前人的结果进行了对比,验证了方法的正确性与有效性.对于几种常见的复杂断层模型,例如弯折、阶跃、含障碍体断层等,我们模拟了其破裂过程并对计算结果进行了比较与分析.模拟结果表明,非结构化网格划分的边界积分方程方法能够很好地模拟平面矩形断层或由其组成的规则断层,同时也能成功地模拟具有复杂几何形状的不规则断层上的动力学破裂过程.本研究的结果显示了边界积分方程方法在模拟复杂断层系统的动力学破裂问题上具有较广阔的应用前景.     

Buoyancy-driven perturbations in a rapidly rotating,electrically conducting fluid: part I – flow and magnetic field

The steady velocity, perturbation pressure and perturbation magnetic field, driven by an isolated buoyant parcel of Gaussian shape in a rapidly rotating, unconfined, incompressible electrically conducting fluid in the presence of an imposed uniform magnetic field, are obtained by means of the Fourier transform in the limit of small Ekman number. Lorentz and inertial forces are neglected. The solution requires at most evaluation of a single integral and is found in closed form in some spatial regions. The solution has structure on two disparate scales: on the scale of the buoyant parcel and on the scale of the Taylor column, which is elongated in the direction of the rotation axis. The detailed structures of the flow and pressure depend linearly on the relative orientation of gravity and rotation, with the solution for arbitrary orientation being a linear combination of two limiting cases in which these vectors are colinear (polar case) and perpendicular (equatorial case). The perturbation magnetic field depends additionally on the relative orientation of the imposed magnetic field, and three limiting cases of interest are presented in which gravity and rotation are colinear (polar–toroidal case), gravity and imposed field are colinear (equatorial–radial case) and all three are mutually perpendicular (equatorial–toroidal case). Visualization and analysis of the velocity and perturbation magnetic field vectors are facilitated by dividing these vector fields into geostrophic and ageostrophic protions. In all cases, the geostrophic and ageostrophic portions have different structure on the Taylor-column scale. The buoyancy force is balanced by a pressure force in the polar case and by a flux of momentum in the equatorial case. The pressure force and momentum flux do not decay in strength with increasing axial distance. Far from the parcel, the axial mass flux varies as the inverse one-third power of distance from the parcel. The velocity has a single geostrophic vortex in the polar case and two vortices in the equatorial case. The perturbation magnetic field has two, four and one geostrophic vortices in the polar–toroidal, equatorial–radial and equatorial–toroidal cases, respectively. To facilitate comparison of the present results with numerical simulations carried out in a finite domain, a set of boundary conditions are developed, with may be applied at a finite distance from the parcel.     

On the Structures Observed in Thin Rotating Layers of a Conductive Fluid and the Anomalies of the Geomagnetic Field

The results of the laboratory and numerical experiments in circular rotating trays with thin layers of a conductive fluid under the MHD generation of small-scale velocity fields are presented. The configurations of constant magnets for MHD generation were determined based on the numerical calculations with shallow water equations. Both the laboratory and numerical experiments with rotating trays demonstrate the emergence of nonaxisymmetric structures and large-scale near-circular vortices caused by the energy transfer from the system of the externally generated small-scale vortices to the large-scale velocity fields under the action of the Coriolis force. The near-circular vortex has areas with differential rotation when the angular velocity of rotation decreases with the radius. The single large-scale vortices and wide jet flows arise in the regimes of subrotation and superrotation relative to the external rotation depending on its angular velocity. The emergence of the flow structures with the azimuthal wave number m = 2 is demonstrated, and their probable relation to the anomalies of the geomagnetic field observed on the Earth’s surface is considered.     

On the exact shape of the horizontal profile of a topographically rectified tidal flow

Abstract

Starting from the nonlinear shallow water equations of a homogeneous rotating fluid we derive the equation describing the evolution of vorticity by a fluctuating bottom topography of small amplitude, using a multiple scale expansion in a small parameter, which is the topographic length scale relative to the tidal wave length. The exact response functions of residual vorticity for a sinusoidal bottom topography are compared with those obtained by a primitive perturbation series and by harmonic truncation, showing the former to be invalid for small topographic length scales and the latter to be only a fair approximation for vorticity produced by planetary vortex stretching. In deriving the exact shape of the horizontal residual velocity profile at a step-like break in the bottom topography, it is shown that the Lagrangian profile only exists in a strip having the width of the amplitude of the tidal excursion at both sides of the break, and that it vanishes outside that interval. Moreover, in the limit of small amplitude topography at least, it vanishes altogether for the generation mechanism by means of planetary vortex stretching. The Eulerian profile is shown to extend over twice the interval of the Lagrangian profile both for production by vortex stretching and by differential bottom friction. These finite intervals over which the residual velocity profiles exist for a step-like topography are not reproduced by harmonic truncation of the basic equation. This method gives exponentially decaying profiles, indicating spurious horizontal diffusion of vorticity. In terms of orders of magnitude, the method of harmonic truncation is reliable for residual velocity produced by vortex stretching but it overestimates the residual velocity produced by differential bottom friction by a factor 2.     

Dynamics of intrathermocline vortices in a gyre flow over a seamount chain

The interaction of meddies with a complex distribution of seamounts is studied in a three-layer quasi-geostrophic model on the f-plane. This study aims at understanding if and how this seamount chain can represent a barrier to the propagation of these eddies and how it can be involved in their decay. The eddies are idealized as vortex patches in the middle layer, interacting with a regional cyclonic current and with ten idealized seamounts. The numerical code is based on the contour surgery technique. The initial position, radius, shape, number and polarity of the eddies are varied. The main results are the following: (1) Though they do not describe the unsteady flow, the streamlines of the regional and topographic flow provide a useful estimate of the vortex trajectories, in particular towards the major seamounts, where stronger velocity shears are expected. (2) The tallest and widest seamounts which have the largest vorticity reservoir are able to considerably erode the vortices, but also to draw anticyclones towards the seamount top. The ability of narrower seamounts to erode vortices is related to their multiplicity. (3) Only 1/3 of the anticyclones with about 30-km radius reach the southern boundary of the seamount chain, and their erosion is larger than 50 %. The other anticyclones are either completely eroded or trapped over a wide seamount top. Cyclones are less affected by seamounts because they oppose the topographic draft towards the seamount top and they drift along the side of the seamount. (4) Large vortices resist topographic erosion more efficiently. The rate of erosion grows from a few percent to about 35–50 % as the vortex radius decreases from about 60 to 30 km. Small cyclones are not eroded, contrary to small anticyclones (which completely decay), in relation with the different trajectories of these eddies in the vicinity of the seamounts. (5) The detailed vortex shape does not appear critical for their evolution, if they are close enough to the seamount chain initially. The interaction between a group of vortices initially north of the seamount chain can modify their trajectory to such an extent that they finally avoid collision with seamounts. (6) Finally, meddy trajectories across the Horseshoe Seamounts (data from the AMUSE experiment) show qualitative similarity with the vortex paths in the model. Several events of vortex decay also occur at comparable locations (in particular over the wide and tall seamounts) in the model and observations.     

On the nonlinear evolution of columnar vortices in a rotating environment

We examine the three-dimensional, nonlinear evolution of columnar vortices in a rotating environment. As the initial vorticity distribution, a wavetrain of finite amplitude Kelvin-Helmholtz vortices in shear is employed. Through direct numerical simulation of the Navier-Stokes equations we seek to better understand the process of maturation of the various three-dimensional modes of instability to which such vortical flows are subject, especially those which exist as a consequence of the action of the Coriolis force. In the absence of rotational influence, we thereby demonstrate that the nonlinear evolution of columnar vortices is most strongly controlled by one or the other of two mechanisms. One mechanism of instability is identifiable as a so-called elliptical instability, which promotes the initial bending of vortex tubes in a sinusoidal fashion, while the other is a hyperbolic mode, which is responsible for the development of streamwise vortex streaks in the "braids" between adjacent vortex cores. In the rotating case, anticyclonic vortices are strongly destabilized by weak background rotation, while rapid rotation stabilizes both the cyclones and anticyclones. The strong anticyclones are subject to two distinct forms of instability, namely a Coriolis force modified elliptical instability and an inertial (centrifugal) instability. The former instability is very similar to the nonrotating form of the elliptical instability as it promotes bending of vortex tubes, while the latter instability grows on the edge of the vortex core and generates streaks of vorticity, which surround the vortex core itself. These results of direct numerical simulation fully verify the results of previous linear stability analyses. Taken together, they provide a simple explanation for the broken symmetry that is often observed to be characteristic of the von Karman vortex streets that develop in the atmospheric lee of oceanic islands.     

A model of secondary upwelling over the shelf break

Motivated by our puzzling high-resolution radar observations of surface vortices in the nearly rectangular Gulf of Eilat/Aqaba, northern Red Sea, we propose and explore the driven cavity approach to this geophysical phenomenon. While the lid-driven cavity has long been considered a benchmark problem in computational fluid dynamics, its oceanographic context has not been considered. Despite the additional effects of rotation and stratification, our modeling demonstrates that when the fluid within a cavity geometrically similar to the Gulf of Eilat is driven by the external current, an interior vortex can develop as in our observations. Furthermore, the Eilat vortices appear only under relatively calm conditions, adding evidence to the intriguing possibility of their simple shear-driven origin.