Internal hydraulic control in rotating fluids—applications to oceans

Abstract

Laboratory experiments and analysis of shallow water equations in a rotating fluid show that channel flow is governed by the ratio of the width of the channel to the Rossby radius of deformation R= √[g&Delta;ρH/ρf ²]. Flows through narrow ocean openings exhibit blocking and clear evidence of hydraulic control. These imply that formulae can be derived for width, volume flux, and velocity scales of the currents. A new version of the constant potential vorticity problem is solved, and it is shown to predict volume flux within 22% of the zero potential vorticity results. Next a systematic method of predicting volume flux through ocean passages is described. Some examples are given from the Denmark Straits overflow and the flow of Antarctic Bottom Water into the western Atlantic Ocean. Two-layer flows and counter-flows with rotation in a narrow passage, the so-called lock exchange flow problem, duplicate flows at a number of important straits and openings to bays. A potential vorticity formulation is reviewed. The flows in the mouths of various bays such as Funka Bay in Hokkaido, Japan, Spencer Gulf in South Australia, and Chesapeake Bay in the United States has R < width of the mouth, and the two currents are separated by a front. The width of the front and the density difference can be predicted with good results.     

Sill effects on physical dynamics in eastern Long Island Sound

This study investigates how Mattituck Sill influences circulation patterns and physical dynamics in eastern Long Island Sound, a major estuary on the U.S. east coast. Observations show there is pronounced across-estuary transport in the area and suggest there may be subtidal anticyclonic flow around the sill. Model runs, with and without sill bathymetry, exhibit this across-estuary transport and anticyclonic circulation. Comparison between these runs indicates that the sill intensifies the anticyclonic circulation. This study finds the sill does not exert internal hydraulic control during neap, mean, or spring tidal conditions. Nevertheless, along-estuary exchange is reduced over the sill and across-estuary fluxes are increased. The Connecticut River plume enters close to the estuary mouth. The sill deflects more of the plume waters towards the mouth, causing less freshwater to take the long looping route through the estuary. The subtidal circulation balance around the sill indicates a barotropic balance between the tidal advection of tidal vorticity and friction. The subtidal vorticity balance indicates the net effect of tidal advection of relative vorticity is balanced with frictional curl associated with lateral speed gradients and vorticity dissipation. Previously developed scalings based on the circulation balance (Nature 290:549–555, 1981), frictional vorticity generation mechanisms (Deep-Sea Res 28:195–212, 1981), and tidal diffusion of potential vorticity (J Phys Oceanogr 29:821–827, 1999) are applicable to Mattituck Sill and predict circulation with a similar magnitudes to model results.     

Instability of a geostrophic front and its energetics

Abstract

The instability of a current with a geostrophic surface density front is investigated by means of a reduced gravity model having a velocity profile with nearly uniform potential vorticity. It is shown that currents are unstable when the mean potential vorticity decreases toward the surface front at the critical point of the frontal trapped waves investigated by Paldor (1983). This instability is identical with that demonstrated by Killworth (1983) in the longwave limit.

The cross-stream component of mass flux and the rates of energy conversions among the five energy forms defined by Orlanski (1968) are also calculated. The main results are as follows, (a) The mass flux toward the surface front is positive near the front and negative around the critical point. The positive mass flux near the front does not vanish at the position of the undisturbed surface front, so that the mean position of the front moves outward and the region of the strong current spreads. (b) The potential energy of the mean flow integrated over the fluid is released through the work done by the force of the pressure gradient of the mean flow on the fluid, and is converted into the kinetic energy of the mean flow. (c) In the critical layer, the mean flow is rapidly accelerated with the growth of the unstable wave. This acceleration is caused by the rapid phase shift of the unstable wave in the critical layer.     

Rotating channel flow: Control and upstream currents

Abstract

Theory and experiments are presented for critically controlled flow of a layer of inviscid rotating fluid. Flow is controlled by a level passage. For a wide upstream channel of fixed depth (i.e. constant potential vorticity) the volume flux on the right-hand wall is unaffected by passage flow. This suggests that specifying Bernoulli potential on the right-hand passage wall produces a physically well-posed condition. The specification results in one less dimensionless number than was required by previous formulations to specify flow in the controlled passage. The upstream flow needs the same number as before, so that a range of upstream conditions produce exactly the same passage flow. A laboratory study is conducted using a thin layer of water under air. This is pumped in steadily at various locations in a deep rotating upstream basin, with fluid leaving through a level passage. All currents in the upstream basin cross to the left-hand wall as the current approaches the passage over a sloping bottom. The current crosses back to the right-hand wall within the passage. Velocity profiles of currents agree reasonably well with constant potential vorticity theory. To the right of the detached upstream current is a closed gyre that connects the upstream flows (that have different patterns depending on source location) with the unique passage flows. The results suggest that gyres upstream of critically controlling passages in the ocean might serve as adjustment regions between the relatively unconstrained upstream flows and the tightly controlled passage flows.     

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.     

On the coupling of fluid dynamics and electromagnetism at the top of the earth's core

Abstract

A kinematic approach to short-term geomagnetism has recently been based upon pre-Maxwell frozen-flux electromagnetism. A complete dynamic theory requires coupling fluid dynamics to electromagnetism.

A geophysically plausible simplifying assumption for the vertical vorticity balance, namely that the vertical Lorentz torque is negligible, is introduced and its consequences are developed. The simplified coupled magnetohydrodynamic system is shown to conserve a variety of magnetic and vorticity flux integrals. These provide costraints on eligible models for the geomagnetic main field, its secular variation, and the horizontal fluid motions at the top of the core, and so permit a number of tests of the underlying assumptions.     

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.     

Stratosphere-troposphere exchange in an extratropical cyclone, calculated with a Lagrangian method

A Lagrangian technique is developed and applied to calculate stratosphere-troposphere exchange in an extratropical cyclone. This exchange is computed from the potential vorticity or PV along trajectories, calculated from ECMWF circulation data. Special emphasis is put on the statistical significance of the results. The computed field of the cross-tropopause flux is dominated by elongated patterns of statistically significant large downward and small upward fluxes. The downward fluxes mainly occur in the lower part of the considered tropopause folds. The upward fluxes are found near the entrance of the folds, in the tropopause ridges. The ratio between the area averaged downward and upward cross-tropopause fluxes increases with increasing strength of the cyclone. Since the largest fluxes are shown to occur in the regions with the largest wind shear, where PV-mixing is thought to cause large cross-tropopause fluxes, the results are expected to be reliable, at least in a qualitative sense. The position of a tropopause fold along the northwest coast of Africa is confirmed by total ozone observations. The results indicate that the applied Lagrangian technique is an appropriate tool for diagnosing stratosphere-troposphere exchange.     

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.     

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.     

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.     

Usefulness of flux measurements during hydraulic tomographic survey for mapping hydraulic conductivity distribution in a fractured medium

Using the first-order analysis, we investigate the spatial cross-correlation between hydraulic conductivity variation and specific discharge (flux) as well as its components measured in a borehole under steady-state flow conditions during cross-hole pumping tests in heterogeneous aquifers. These spatial correlation patterns are found to be quite different from that between the hydraulic conductivity variation and the hydraulic head measurement in the same borehole. This finding suggests that a specific discharge measurement carries non-redundant information about the spatial distribution of heterogeneity, even this measurement is collected from the same location where the head measurement is taken. As such, specific discharge observations should be included in the analysis of hydraulic tomography to increase the resolution of estimated aquifer heterogeneity. Using numerical experiments, we demonstrate the effectiveness of the joint interpretation of both hydraulic heads and fluxes for mapping fracture distributions in a hypothetic geologic medium.     

Numerical simulation of groundwater flow patterns using flux as upper boundary

Since the 1960s, most of the studies on groundwater flow systems by analytical and numerical modelling have been based on given‐head upper boundaries. The disadvantage of the given‐head approach is that the recharge into and discharge from a basin vary with changes in hydraulic conductivity and/or basin geometry. Consequently, flow patterns simulated with given‐head boundaries but with different hydraulic conductivities and/or basin geometry may not reflect the effects of these variables. We conducted, therefore, numerical simulations of groundwater flow in theoretical drainage basins using flux as the upper boundary and realistically positioned fluid‐potential sinks while changing the infiltration intensity, hydraulic conductivities, and geometric configuration of the basin. The simulated results demonstrate that these variables are dominant factors controlling the flow pattern in a laterally closed drainage basin. The ratio of infiltration intensity to hydraulic conductivity (R_ic) has been shown to be an integrated pattern‐parameter in a basin with a given geometric configuration and possible fluid‐potential‐sink distribution. Successively, the changes in flow patterns induced by stepwise reductions in R_ic are identical, regardless of whether the reductions are due to a decrease in infiltration intensity or an increase in hydraulic conductivity. The calculated examples show five sequential flow patterns containing (i) only local, (ii) local–intermediate, (iii) local–intermediate–regional, (iv) local–regional, and (v) just regional flow systems. The R_ic was found to determine also whether a particular sink is active or not as a site of discharge. Flux upper boundary is preferable for numerical simulation when discussing the flow patterns affected by a change of infiltration, the hydraulic conductivity, or the geometry of a basin. Copyright © 2012 John Wiley & Sons, Ltd.     

Flows within the Ekman layer during spin-up of a thermally stratified fluid

Abstract

Finite-difference numerical solutions were obtained to present the flow and temperature field details within the transient Ekman layer during spin-up of a thermally stratified fluid in a cylinder. This complements the earlier studies on stratified spin-up which examined the flows in the interior core region. As the stratification increases, the following changes in the flow field are noticeable. The radial velocity in the Ekman layer decreases in magnitude. The azimuthal flows adjust smoothly from the interior region to the endwall boundary, and the Ekman layer in the azimuthal flow field fades. Vertical motions are inhibited, resulting in a weakened Ekman pumping. The axial vorticity field behaves similarly to the azimuthal flows. The temperature deviation from the equilibrium profile decreases, and the heat transfer flux from the endwall to the fluid decreases. The thickness of the thermal layer is larger than the velocity layer thickness. Illustrative comparisons of the relative sizes of the terms in the governing equations are conducted in order to assess the stratification effect in the adjustment process of the fluid.     

Non-hydrostatic and non-linear contributions to the internal wave energy flux in sill regions

A three-dimensional non-linear, non-hydrostatic model in cross-sectional form is used to determine the factors influencing the relative importance of the linear, non-hydrostatic and non-linear contributions to the internal wave energy flux in sill regions due to tidal forcing. The importance of the free surface elevation term is also considered. Idealised topography representing the sill at the entrance to Loch Etive, the site of a recent measurement programme, is used. Calculations show that the non-linear terms in the energy flux become increasingly important as the sill Froude Number (F _s) increases and the sill aspect ratio is increased. The vertical profile of the stratification, in particular its value close to the sill crest where internal waves are generated, has a significant influence on unsteady lee wave and mixed tidal–lee wave generation and the non-linear contribution to the energy flux. Calculations show that as F _s increases, the energy flux due to the non-linear and non-hydrostatic terms increases more rapidly than the linear term. The importance of the non-linear terms in the energy flux also increases as the sill aspect ratio is increased. Increasing the buoyancy frequency reduces the contribution of the non-hydrostatic and non-linear terms to the total energy flux. Also, as the buoyancy frequency is increased, this reduces unsteady lee wave and mixed tidal–lee wave generation. In essence, these calculations show that the energy flux due to the non-hydrostatic and non-linear terms is appreciable in sill regions.     

Numerical modelling of stratified tidal flow over a fjord sill

Field observations of tidally driven stratified flow in the sill area of Knight Inlet (British Columbia) revealed a very complicated structure, which includes solitary waves, upstream bifurcation, hydraulic jump and mixing processes. Recent observations suggest that the flow instabilities on the plunging pycnocline at the lee side of the sill may contribute to solitary wave generation through a subharmonic interaction. The present study reports on a series of numerical experiments of stratified tidal flow in Knight Inlet performed with the help of a fine resolution fully non-linear non-hydrostatic numerical model. The model reproduces all important stages of the baroclinic tidal dynamics observed in Knight Inlet. Results demonstrate that solitary waves are generated apart from the area of hydrodynamic instability. Accelerating tidal flux forms a baroclinic hydraulic jump just above the top of the sill, whereas the bifurcations and zones of shear instabilities are formed downstream of the sill. The first baroclinic mode having the largest velocity escapes from the generation area and propagates upstream, disintegrating further into a packet of solitary waves reviling the classical “non-subharmonic” mechanism of generation. The remaining part of the disturbance (slow baroclinic modes) is arrested by tidal flow and carried away to the lee side of the obstacle, where shear instability, billows and mixing processes are developed. Some sensitivity runs were performed for different value of tidal velocity.     

Advective flow in sediments under the influence of compaction

Abstract

Models on flow and transport in surface water sediments currently neglect compaction, although it is well understood that compaction is one of the major processes below the free fluid-sediment interface. Porosity changes in the sediment layers, as a result of compaction, are measured in almost all probes: porosity decreases with the distance from the surface water-sediment interface. This paper provides a rigorous derivation of basic flux terms for a frame of reference that is moving with the fluid-sediment interface. It is shown how burial rate, interface velocity, velocities of fluid and solid phase and porosity are connected—under steady-state conditions. It turns out that porosity and the velocities in a one-dimensional column can be directly computed from each other. These findings are important not only for the understanding of compaction-driven flow itself; they are crucial for all studies on storage and transport of chemical components in sediments. As mass fluxes across the sediment-water interface may be affected, there is an indirect link on surface water quality, making these findings relevant also for research on eutrophication of surface water bodies and/or on biogeochemical cycles.     

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.