Hamilton's principle and the vorticity laws for a relativistic perfect fluid

Abstract

The equations for a relativistic perfect fluid result from the requirement that the total mass-energy be stationary with respect to variations δx^α(a, b, c, s) in the space-time location of the fluid particle identified by Lagrangian labels (a, b, c) at the point s on its world-line. By considering variations of the Lagrangian labels that leave the specific volume and entropy unchanged, we obtain a general covariant statement of vorticity conservation. The conservation laws for circulation, potential vorticity, and helicity are simple corollaries. This Noether-theorem derivation shows that the vorticity laws have no analogues in particle mechanics, where the corresponding particle labels cannot be continuously vaned.     

The generation and decay of vorticity

Abstract

Vorticity, although not the primary variable of fluid dynamics, is an important derived variable playing both mathematical and physical roles in the solution and understanding of problems. The following treatment discusses the generation of vorticity at rigid boundaries and its subsequent decay. It is intended to provide a consistent and very broadly applicable framework within which a wide range of questions can be answered explicitly. The rate of generation of vorticity is shown to be the relative tangential acceleration of fluid and boundary without taking viscosity into account and the generating mechanism therefore involves the tangential pressure gradient within the fluid and the external acceleration of the boundary only. The mechanism is inviscid in nature and independent of the no-slip condition at the boundary, although viscous diffusion acts immediately after generation to spread vorticity outward from boundaries. Vorticity diffuses neither out of boundaries nor into them, and the only means of decay is by cross-diffusive annihilation within the fluid.     

Comment on rotating hydraulics

Abstract

Numerical solutions of the axisymmetric flows during the relatively early phase of spin-up from rest of a stratified fluid in a cylinder are presented. Detailed results are given for a cylinder of aspect ratio of O(l) and for a minute Ekman number, showing axisymmetric spin-up for three values of the stratification parameter. As the stratification increases, the meridional circulation is confined to a region closer to the Ekman layers. An axisymmetric shear wave propagates radially inward from the sidewall, but, unlike the strictly vertical front for a homogeneous fluid, the interface which separates rotating from nonrotating fluid is bow-shaped. For a stratified fluid, the axial vorticity distribution is nonuniform both in the vertical and in the radial directions. With increasing stratification, diffusive vorticity production near the sidewall is more pronounced. Axisymmetric flows in the early phase of spin-up of a stratified fluid are controlled by both the inviscid dynamic effect and the viscous diffusion effect. At a location close to the Ekman layers, the inviscid effect outweighs the viscous effect, in much the same way as in a homogeneous fluid. However, at a location close to mid-depth, the viscous diffusion effect, enhanced by substantial flow gradients in that region, is dominant. This points to the necessity of including the direct effect of viscous diffusion in the interior in formulating an analytical model of stratified spin-up problems.     

Rapid formation of taylor columns: Obstacles against sidewalls

Abstract

The problem of topographic forcing by an obstacle against the boundary of a rotating flow is considered in various parameter regimes. The timescale for the motion is the topographic vortex-stretching time, which is inversely proportional to the background rotation rate and the fractional height of the obstacle. For slow flows this time is short compared with the advection time and the governing equation of conservation of potential vorticity is linear. The final state satisfies the non-linear equation for the advection of potential vorticity, however, and so time dependence has given a specific solution to a non-linear problem. The presence of the sidewall causes a stagnant Taylor column to be set up far more rapidly than cases with no sidewall. It is shown that viscosity and mixing arrests the inviscid evolution at some stage, thus some fluid still crosses the obstacle in the steady state. These solutions suggest that experimental results on separation obtained by Griffiths and Linden (1983) can tentatively be ascribed to entrainment and expulsion of fluid through vertical shear layers at the edge of the topography.     

A note on the free-surface effect on the topographically induced vorticity field in a homogeneous rotating flow

Abstract

The vorticity field induced by a topography in a homogeneous rotating fluid bounded by a free-surface is investigated. It is found that the integrated strength of the topographically bound vortex is zero and that the volume of the free-surface perturbations is equal to the volume of the topography.     

On reduced-gravity flow through sills

Abstract

The hydraulic flow of a reduced-gravity fluid of non-negative potential vorticity through a sill is considered. It is shown that for any flow with a reversal of current, another, physically realisable, flow exists which is unidirectional and/or resting, and carries more flux than the original flow. Thus only non-negative flows need be considered when examining maximal hydraulic fluxes. Then, for a simple sill (one which slopes downward on the left and upward on the right, looking downstream), it is shown that zero potential vorticity flow, possibly modified by having a region of motionless fluid at its right, carries the maximum flux possible for that sill shape. This makes the calculation of maximal fluxes for a given sill considerably simpler, and examples of parabolic and V-shaped sills are computed.     

Superhelicity,helicity and potential vorticity

Abstract

‘‘Helicity'’ density H ≡ u · ω and other pseudo-scalar fields such as P ≡ ω · Vlnρ (which is related to Ertel potential vorticity) are useful quantities in theoretical fluid dynamics and magneto-fluid dynamics. Here u denotes the Eulerian flow velocity relative to the chosen frame of reference, ω ≡ V × u is the corresponding relative vorticity and ρ the mass density of the fluid. A general expression is readily obtained for ?H/?t (where t denotes time) in terms of P and the ‘‘superhelicity'’ density S ≡ ω · V × ω which, in fluids of low viscosity, has its highest values in boundary layers. One need for such a relationship became evident during an attempt to interpret the findings of laboratory experiments on thermal convection in rotating fluids in containers of various geometrical shapes and topological characteristics.

In electrodynamics an analogous expression can be found relating the time rate of change of ‘‘magnetic helicity'’ A · B to ‘‘magnetic superhelicity'’ B · ? × B (where B · ? × A is the magnetic field) and a scalar quantity analogous to P which involves non-Ohmic contributions to the relationship between the electric current density and the electric field.     

On the steady mass transport induced by wave motions in a viscous rotating fluid

Abstract

The steady second-order motion induced by a first-order wave motion in a homogeneous, viscous and rotating fluid is examined. If the wave motion produces a steady Ekman layer suction by non-linear interactions, this suction must induce a steady component of interior, relative vorticity parallel to the axis of rotation in order to conserve mass. A boundary value problem for the determination of the induced, steady interior mass transport velocity is presented. The mass transport induced by a Kelvin wave is examined as an illustration and possible application of the theory.     

Displacement and rectification of planetary fluids

Abstract

Laboratory experiments on the decay (spin-up) of fluid motion on the β-plane are compared with theory. Under weakly dissipative conditions, some particles conserve potential vorticity during the entire decay. We also study the rectified mean flow which is produced by the lateral Reynolds Stress when a low frequency force is applied to the planetary fluid. The possible connection of effects with oceanic phenomena is briefly discussed.     

Statistical properties of an enstrophy conserving finite element discretisation for the stochastic quasi-geostrophic equation

ABSTRACT

A framework of variational principles for stochastic fluid dynamics was presented by Holm, and these stochastic equations were also derived by Cotter, Gottwald and Holm. We present a conforming finite element discretisation for the stochastic quasi-geostrophic equation that was derived from this framework. The discretisation preserves the first two moments of potential vorticity, i.e. the mean potential vorticity and the enstrophy. Following the work of Dubinkina and Frank, who investigated the statistical mechanics of discretisations of the deterministic quasi-geostrophic equation, we investigate the statistical mechanics of our discretisation of the stochastic quasi-geostrophic equation. We compare the statistical properties of our discretisation with the Gibbs distribution under assumption of these conserved quantities, finding that there is an agreement between the statistics under a wide range of set-ups.     

Topographic eddies in temporally varying oceanic flows

Abstract

In this paper we examine the behaviour of oceanic unsteady flow impinging on isolated topography by means of numerical simulation. The ocean model is quasigeostrophic and forced by an oscillatory mean flow. The fluid domain is of the channel type and open-boundary numerical conditions are used to represent downstream and upstream flow.

In certain cases, vortex shedding, either cyclonic or anticyclonic, is observed in the lee of obstacles. Such shedding can be explained as the consequence of both an enhanced process of vorticity dissipation over the topography which locally affects the balance of potential vorticity on the advective timescale, and a periodic dominance of advective effects which sweep the fluid particles trapped on the seamount. For refined resolution and smallest viscosity the model will predict flows in which the shed eddies are coherent structures with closed streamlines.

The model suggests a mechanism by which topographically generated eddies may be swept away from a seamount in the ocean.     

Vortex evolution

Abstract

Friedmann's equation and the potential vorticity equation are generalised for turbulent motion. The generalised equations incorporate some new phenomena connected with turbulent transport of mass. It is proved that, if ?×[S×Ω+S(?·S)]≠0 where Ω is the absolute vorticity of the velocity and S is the turbulent density flux, then the Helmholtz-Kelvin theorem concerning the conservation of the velocity circulation around a closed path is violated and the potential vorticity is not a Lagrangian adiabatic invariant. The effects of this turbulent transport of mass on the creation or dissipation of vorticity discussed here is not equivalent to effects of baroclinicity or viscosity. Some possible implications of the new circulation theorem in geophysical and astrophysical fluid dynamics are discussed.     

Laboratory experiments on fronts

Abstract

Supercritically unstable density fronts near a vertical wall in a rotating, two-layer fluid were created on a laboratory turntable by withdrawing the outer wall of an annulus with a narrow gap, and allowing buoyant fluid from within the annulus to collapse toward a state of quasi-geostrophic balance. The resulting “coastal” current has a nearly uniform potential vorticity and is bounded by a front on which ageostrophic, wave-like disturbances grow. If the current width is comparable to the Rossby radius of deformation, the dominant length scale of disturbances is proportional to the width of the current. On the other hand, if the upper layer is much wider than the Rossby radius, then the observed length scale is a constant multiple of the Rossby radius. If the vertical boundary is omitted in the experiments, so that we are left with a circular anticyclonic vortex, the observed length scales and large-amplitude behaviour of disturbances are identical to those for the boundary currents, indicating that the wall has no significant influence on the flow.

At very large amplitude the growing waves lead to the formation of cyclone-anticyclone vortex pairs. For very wide currents, both the mean flow and the disturbances are first confined to a region within a few Rossby radii of the front. However, both the mean flow and the turbulent eddy motions slowly propagate into the previously stationary upper layer until, eventually, the whole of the upper layer is turbulent.     

Rectification of the wind-driven ocean circulation on the beta plane

Abstract

A nonlinear Stommel model of the ocean circulation on the beta plane, driven by a time periodic wind stress, is investigated in order to study symmetry properties of the observed time-mean ocean gyres. Due to the presence of vorticity advection terms the model will have a steady or rectified response to fluctuating wind fields. In this paper a small inverse Ekman number, “the small beta regime”, is considered. It is demonstrated that for this case all qualitative features of the residual circulation, obtained numerically by Veronis (1970). are reproduced in an analytical way. They include the dipole character of the gyre, its maximum symmetry breaking around the north-south axis for intermediate Reynolds numbers, measuring the ratio of vorticity forcing and dissipation, and the maximum residual response for intermediate forcing frequencies.     

Residual vorticity induced by the action of tidal currents in combination with bottom topography in the Southern Bight of the North Sea

Abstract

Results are presented of calculations on the generation of residual vorticity by tidal currents over the bottom topography of the Southern Bight of the North Sea. A typical order of magnitude is 10^?6 to 10^?7 s ^?1. This is compared with current measurements on calm days, when similar magnitudes are found. At windspeeds less than about 5 m/s tidal generation of residual vorticity is important; at higher windspeeds wind effects begin to dominate. Our results are relevant in understanding the spatial variability of residual currents, because a non-zero vorticity implies the existence of horizontal gradients in the residual current field.     

On the normal mode instability of modons and Wu-Verkley waves

Abstract

The normal mode instability of steady Wu-Verkley (1993) wave and modons by Verkley (1984, 1987, 1990) and Neven (1992) is considered. All these flows are solutions to the vorticity equation governing the motion of an ideal incompressible fluid on a rotating sphere. A conservation law for infinitesimal perturbations to each solution is derived and used to obtain a necessary condition for its exponential instability. By these conditions, Fjörtoft's (1953) average spectral number of the amplitude of an unstable mode must be equal to a specific number that depends on the degree of the solution in its inner and outer regions as well as on spectral distribution of the mode energy in these regions. Some properties of the conditions for different types of modons are discussed. The maximum growth (and decay) rate of the modes is estimated, and the orthogonality of the amplitude of each unstable, decaying, or non-stationary mode to the basic solution is shown in the energy inner product.

The new instability conditions confine the unstable disturbances of the WV wave and modon to a hypersurface in the perturbation space and allow interpretation of their energy structure. They are also useful both in estimating the maximum growth rate of unstable modes and in testing the numerical algorithms designed for the linear stability study.     

On the generation of baroclinic rossby waves in the presence of a semi-circular boundary

Abstract

An analytical model is constructed for the generation of baroclinic Rossby waves by a vorticity source in the presence of a semi-circular boundary. The vorticity source is used to represent the effect of the Agulhas retroflection to the south of Southern Africa. The displacement of the interface between the two layers of the model ocean consists of quantized waves near the coast and a train of Rossby waves drifting westward further offshore.     

Stability and stable modes of coastal fronts

Abstract

The stability of a single layer, geostrophic front of zero potential vorticity bounded by a vertical coast (wall) is investigated by means of a Rayleigh integral. It is proved that the flow of the density-driven current is stable at all wavenumbers provided the mean velocity of basic flow exceeds fL (where f is the Coriolis parameter and L is the distance between the wall and the free streamline). The frequency of the stable long waves is either zero or super-inertial.     

A note on the use of zero potential vorticity models

Abstract

The acceptability of zero potential vorticity models as approximations for natural systems of small, but finite, potential vorticity is studied for bounded frontal flows of arbitrary profile. It is demonstrated that all (infinitely) long-wave solutions of the zero potential vorticity front are asymptotic limits for some (not necessarily long-wave) solutions of the small potential vorticity front. In contrast, for downstream-varying solutions there is no simple way of demonstrating this property. These findings suggest that the use of zero potential vorticity models should be carefully examined in other, non-frontal, problems as well. Finally we show that the longwave solutions of the zero potential vorticity flow are at most neutral (quasi-stable).     

The nonlinear spin-up of a stratified ocean

Abstract

The adjustment of a nonlinear, quasigeostrophic, stratified ocean to an impulsively applied wind stress is investigated under the assumption that barotropic advection of vortex tube length is the most important nonlinearity. The present study complements the steady state theories which have recently appeared, and extends earlier, dissipationless, linear models.

In terms of Sverdrup transport, the equation for baroclinic evolution is a forced advection-diffusion equation. Solutions of this equation subject to a “tilted disk” Ekman divergence are obtained analytically for the case of no diffusion and numerically otherwise. The similarity between the present equation and that of a forced barotropic fluid with bottom topography is shown.

Barotropic flow, which is assumed to mature instantly, can reverse the tendency for westward propagation, and thus produce regions of closed geostrophic contours. Inside these regions, dissipation, or equivalently the eddy field, plays a central role. We assume that eddy mixing effects a lateral, down-gradient diffusion of potential vorticity; hence, within the closed geostrophic contours, our model approaches a state of uniform potential vorticity. The solutions also extend the steady-state theories, which require weak diffusion, by demonstrating that homogenization occurs for moderately strong diffusion.

The evoiution of potential vorticity and the thermocline are examined, and it is shown that the adjustment time of the model is governed by dissipation, rather than baroclinic wave propagation as in linear theories. If dissipation is weak, spin-up of a nonlinear ocean may take several times that predicted by linear models, which agrees with analyses of eddy-resolving general circulation models. The inclusion of a western boundary current may accelerate this process, although dissipation will still play a central role.