Diffusive destabilisation of the baroclinic circular vortex

Abstract

It is shown that, even for vanishingly small diffusivities of momentum and heat, a rotating stratified zonal shear flow is more unstable to zonally symmetric disturbances than would be indicated by the classical inviscid adiabatic criterion, unless σ, the Prandtl number, = 1. Both monotonic instability, and growing oscillations ("overstability") are involved, the former determining the stability criterion and having the higher growth rates. The more σ differs from 1, the larger the region in parameter space for which the flow is stable by the classical criterion, but actually unstable.

If the baroclinity is sufficiently great for the classical criterion also to indicate instability, the corresponding inviscid adiabatic modes usually have the numerically highest growth rates. An exception is the case of small isotherm slope and small σ.

A single normal mode of the linearized theory is also, formally, a finite amplitude solution; however, no theoretical attempt is made to assess the effect of finite amplitude in general. But, in a following paper, viscous overturning (the mechanism giving rise to the sub‐classical monotonic instability when σ > 1) is shown to play an important role at finite amplitude in certain examples of nonlinear steady thermally‐driven axisymmetric flow of water in a rotating annulus. Irrespective of whether analogous mechanisms turn out to be identifiable and important in large‐scale nature, it appears then that a Prandtl‐type parameter should enter the discussion of any attempt to make laboratory or numerical models of zonally‐symmetric baroclinic geophysical or astrophysical flows.     

Instability of flows near the onset of convection in a rotating layer with stress-free horizontal boundaries

We investigate instability of convective flows of simple structure (rolls, standing and travelling waves) in a rotating layer with stress-free horizontal boundaries near the onset of convection. We show that the flows are always unstable to perturbations, which are linear combinations of large-scale modes and short-scale modes, whose wave numbers are close to those of the perturbed flows. Depending on asymptotic relations of small parameters α (the difference between the wave number of perturbed flows and the critical wave number for the onset of convection) and ε (ε² being the overcriticality and the perturbed flow amplitude being O(ε)), either small-angle or Eckhaus instability is prevailing. In the case of small-angle instability for rolls the largest growth rate scales as ε^8/5, in agreement with results of Cox and Matthews (Cox, S.M. and Matthews, P.C., Instability of rotating convection. J. Fluid. Mech., 2000, 403, 153–172) obtained for rolls with k = k _c . For waves, the largest growth rate is of the order ε^4/3. In the case of Eckhaus instability the growth rate is of the order of α².     

Finite amplitude holmboe waves

Abstract

We investigate the evolution of a parallel shear flow which has embedded within it a thin, symmetrically positioned layer of stable density stratification. The primary instability of this flow may deliver either Kelvin-Helmholtz waves or Holmboe waves, depending on the strength of the stratification. In this paper we describe a sequence of numerical simulations which reveal for the first time the behavior of the Holmboe wave at finite amplitude and clarify its structural relationship to the Kelvin-Helmholtz wave.

The flows investigated have initial profiles of horizontal velocity and Brunt-Vaisala frequency given in nondimensional form by U = tanhζ and N ²=J sech² RCζ, respectively, in which ζ is a nondimensional vertical coordinate, J is the value of the gradient Richardson number N ²/(dU/dζ)² at ζ=0, and R = 3. Linear stability theory predicts that the flow will develop Holmboe instability when J exceeds some critical value J_c' and Kelvin-Helmholtz instability when J is less than J_c; J_c being approximately equal to 0.25 when R=3. We simulate the evolution of flows with J=0.9, J=0.45, and J = 0.22, and find that the first two simulations yield Holmboe waves while the third yields a Kelvin-Helmholtz wave, as predicted.

The Holmboe wave is a superposition of two oppositely propagating disturbances, a right-going mode whose energy is concentrated in the region above the centre of the shear layer, and a left-going mode whose energy is concentrated below the centre of the shear layer. The horizontal speed of the modes varies periodically, and the variations are most pronounced at low values of J. If J ζ J_c' the minimum horizontal speed of the modes vanishes and the modes become phase-locked, whereupon they roll up to form a Kelvin-Helmholtz wave as predicted by Holmboe (1962). When J is moderately greater than J_c' the Holmboe wave ejects long, thin plumes of fluid into the regions above and below the shear layer, as has often been observed in laboratory experiments, and we examine in detail the mechanism by which this occurs.     

The instability of barotropic circular vortices

Abstract

The linear, normal mode instability of barotropic circular vortices with zero circulation is examined in the f-plane quasigeostrophic equations. Equivalents of Rayleigh's and Fjortoft's criteria and the semicircle theorem for parallel shear flow are given, and the energy equation shows the instability to be barotropic. A new result is that the fastest growing perturbation is often an internal instability, having a finite vertical scale, but may also be an external instability, having no vertical structure. For parallel shear flow the fastest growing perturbation is always an external instability; this is Squire's theorem. Whether the fastest growing perturbation is internal or external depends upon the profile: for mean flow streamfunction profiles which monotonically decrease with radius, the instability is internal for less steep profiles with a broad velocity extremum and external for steep profiles with a narrow velocity extremum. Finite amplitude, numerical model calculations show that this linear instability analysis is not valid very far into the finite amplitude range, and that a barotropic vortex, whose fastest growing perturbation is internal, is vertically fragmented by the instability.     

Dynamic simulation on hydraulic characteristic values of overland flow

The economic forest management is one of the main land use models on low hill gentle slope. In order to investigate the soil erosion properties of bare slope under economic forest, dynamic simulation on hydraulic characteristic values of overland flow was carried out under 0.5 mm min^?1, 1.2 mm min^?1 and 1.8 mm min^?1 rainfall intensities. Results indicated that runoff shear stress increased with increasing of slope length and their relationship can be described by quadratic equation. There were abnormal points at the length of 4 m and 5.5 m under rainfall intensity of 1.8 mm min^?1. The shallow flow was pseudo-laminar flow under 0.5 mm min^?1, 1.2 mm min^?1 and 1.8 mm min^?1 rainfall intensities, and the runoff at upslope was sluggish flow then changed to torrential flow at downslope with increasing of slope length. Critical Reynolds number varied from sluggish flow to torrential flow with 1.8 mm min^?1 rainfall intensity and was more than that under 0.5 mm min^?1. Reynolds number can be estimated by power function of slope length. And there was a positive correlation between runoff shear stress and both Froude number Fr and Reynolds number Re. We hope this study can provide scientific gist for soil erosion control under economic forest.     

An experimental study of global instabilities due to the tidal (elliptical) distortion of a rotating elastic cylinder

Abstract

Broad band secondary instability of elliptical vortex motion has been proposed as a principal source of shear-flow turbulence. Here experiments on such instability in an elliptical flow with no shear boundary layer are described. This is made possible by the mechanical distortion in the laboratory frame of a rotating fluid-filled elastic cylinder. One percent ellipticity of a 10 cm diameter cylinder rotating once each second can give rise to an exponentially-growing mode stationary in the laboratory frame. In first order this mode is a sub-harmonic parametric Faraday instability. The finite-amplitude equations represent angular momentum transfer on an inertial time scale due to Reynolds stresses. The growth of this mode is not limited by boundary friction but by detuning and centrifugal stabilization. On average, a generalized Richardson number achieves a marginal value through much of the evolved flow. However, the characteristic flow is intermittent with the cycle: rapid growth, stabilizing momentum transfer from the mean flow, interior re-spin up, and then again. Data is presented in which, at large Reynolds numbers, seven percent ellipticity causes a fifty percent reduction in the kinetic energy of the rotating fluid. In the geophysical setting, this tidal instability in the earth's interior could be inhibited by sub-adiabatic temperature gradients. A near adiabatic region greater than 10 km in height would permit the growth of tidally destabilized modes and the release of energy to three-dimensional disturbances. Such disturbances might play a central role in the geodynamo and add significantly to overall tidal dissipation.     

Axisymmetric annulus convection at unit Prandtl number

Abstract

We examine the role played in annulus flows by mechanisms dependent upon the Prandtl number, σ. Solutions are obtained at σ = 1 for both the real annulus system and for the hypothetical “free annulus” system (free slip lateral boundaries). These solutions are compared with previously obtained solutions at σ = 7.

In the free annulus, the solution at σ = 1 differs radically from that at σ = 7. The σ = 1 solution appears to be essentially a finite amplitude mode due to Solberg instability whereas the solution at σ = 7 manifests a flow caused by the diffusive overturning mechanism.

The variation with σ of the real annulus flow is not so fundamental but some differences in the dynamical structures are noted.     

The stability of dissipative magnetohydrodynamic shear flow in a parallel magnetic field

Abstract

Strong decay bounds are obtained for linearized perturbations to an unbounded, plane Couette flow in a parallel magnetic field. Finite conductivity and molecular viscosity are found to be stabilizing. Those modes decaying most slowly have the form of rolls aligned with the shear flow. The non-aligned rolls decay at an enhanced rate. Stability bounds at finite amplitude are obtained for flows bounded in one direction using energy methods.     

Three-dimensional primary instabilities of a stratified,dissipative, parallel flow

Abstract

We demonstrate the existence of a class of dissipative, stratified, parallel shear flows which, as a consequence of linear supercritical instability, evolve directly into three-dimensional flows without the requirement for an intermediate two-dimensional finite-amplitude state. This represents a counter-example to a common misinterpretation of Squire's theorm, namely that the fastest-growing unstable mode of a dissipative parallel shear flow must be two-dimensional.     

Dynamos on stream surfaces of a highly conducting fluid

Abstract

This paper discusses dynamo action in generalisations of the Ponomarenko dynamo at large magnetic Reynolds number. The original Ponomarenko dynamo consists of a spiralling flow in which the stream surfaces are concentric cylinders of circular cross section, and the flow depends only on distance from the axis in cylindrical polar coordinates.

In this study, the stream surfaces are allowed to be cylinders of arbitrary cross section, and the flow is only required to be independent of the coordinate along the cylinder axes. For smooth flows alpha and eddy diffusion effects are identified, in terms of the geometry of the stream surfaces, and asymptotic formulae for growth rates in the limit of large magnetic Reynolds number are obtained. Numerical support for these results is presented using direct simulation of dynamo action in selected flows at high conductivity. Finally the case is considered when in spherical polar coordinates the flow is independent of the azimuthal coordinate and the stream surfaces, which are tori, have arbitrary cross sections.     

The impact of flow discharge on the hydraulic characteristics of headcut erosion processes in the gully region of the Loess Plateau

Headcut erosion is associated with major hydraulic changes induced by the gully head of concentrated flow. However, the variation in the hydraulic characteristics of the headcut erosion process is still not clear in the gully region of the Loess Plateau. A series of rainfall combined scouring experiments (flow discharges ranging from 3.6 to 7.2 m³ hr⁻¹, with 0.8 mm min⁻¹ rainfall intensity) were conducted on experimental plots to clarify the variation in the hydraulic parameters induced by gully head and erosion processes under different flow discharges. The results showed that concentrated flows in the catchment area and gully bed were turbulent (Reynolds number ranging from 1,876 to 6,693) and transformed between supercritical and subcritical (Froude number ranging from 0.96 to 3.73). The hydraulic parameters, such as the flow velocity, Reynolds number, shear stress, stream power, Darcy–Weisbach friction factor, and unit stream power in the catchment area were 0.45–0.59 m s⁻¹, 2086–6693, 1.96–5.33 Pa, 0.89–2.86 W m⁻², 0.08–0.16, and 0.023–0.031 m s⁻¹, respectively. When the concentrated flows dropped from the gully head, the hydraulic parameters in the gully bed decreased by 3.39–26.07%, 1.49–29.99%, 65.19–67.14%, 67.25–74.96%, 28.53–61.31%, and 67.82–77.14%, respectively, which contributed to the flow energy consumption at the gully head. As flow discharge increased, Reynolds number, shear stress, and stream power increased, while flow velocity, Froude number, unit stream power, and Darcy–Weisbach friction factor did not. The flow energy consumption at the gully head was 9.66–10.13, 13.25–13.74, 15.68–16.41, and 19.28–20.25 J s⁻¹, respectively, under different flow discharges and accounted for 60.58–68.50% of the flow energy consumption of the experimental plots. Generally, the sediment discharges increased rapidly at the initial stage, then increased slowly, and finally reached a steady state condition, which showed a significant declining logarithmic trend with experimental duration (P<.01) and increased with increasing flow discharge. Accordingly, the flow energy consumption was significantly correlated with the sediment yield. These findings could improve our understanding of the hydraulic properties and flow energy characteristics of headcut erosion.     

Stability of buoyancy opposed mixed convection in a vertical channel and its dependence on permeability

Non-Darcy mixed convective flow of water due to external pressure gradient and buoyancy opposed forces are considered in a vertical channel filled with porous medium, which can be either isotropic or anisotropic. The linear theory of stability analysis has been used to numerically investigate the dependence of the transition behavior of the fully developed basic flow on the permeability of the medium. Numerical experiments indicate that mainly two main instability modes appear: Rayleigh–Taylor (R–T) and buoyant instability. For Darcy numbers (Da) ?10⁻⁹, R–T instability dominates within the entire Reynolds number (Re) range considered here. It was also found that for the same Re, the fully developed base flow is highly unstable (stable) for porous media with high (low) permeability. Further, it was seen that the disturbance isotherm cells migrate from the channel walls toward the centerline when permeability is reduced. Reducing the permeability by one order of magnitude (corresponding to a decrease of Darcy number from 10⁻⁶ to 10⁻⁷) increases base flow stability approximately 20-fold. For higher Reynolds numbers, buoyant, mixed and shear instability of the basic flow were found when Da was increased from 10⁻⁷ to 10⁻³. However, for cases in which permeability and porosity behaved as suggested by Carman–Kozeny relation (CKR), buoyant stability was the only mode of instability. Critical values of the Rayleigh (Ra) and Darcy (Da) numbers in the R–T mode of instability were related to each other by the hyperbolic function RaDa = −2.465.     

Hot springs mediate spatial exchange of heat and mass in the enclosed sediment domain: A stability perspective

In many natural environments, such as in underwater hot springs and hydrothermal vents, thermal gradients are accompanied with changes in the concentration of chemical compounds transported to the seawater, causing the so-called double-diffusive, mixed convection. To study the physical scenarios in such systems, a vertical channel filled with a porous medium saturated with saline water is considered. The motion in the sediment-filled channel is induced by two buoyancy forces and an external pressure gradient, similar to the situation in a vent with an upward flow direction. The fluid flow has been modeled by an extended Darcy model, and the flow instability mechanisms have been studied numerically. The linear stability analysis is performed considering a wide range of Darcy number (Da = 10⁻⁵ -10⁻⁸). The instability boundary curve showed three distinct dynamic regimes: (i) Rayleigh-Taylor (R-T), (ii) log-log non-linear variation, and (iii) log-log linear variation. The domain of different regimes were sensitive to external pressure gradient as well as permeability. Similar to cross-diffusive natural convection in pure viscous fluids, a linear relationship between logarithmic absolute values of critical thermal Rayleigh number (∣Ra_T∣) and solute Rayleigh number (Ra_C) is found in the third regime. Based on the permeability, for any solute Rayleigh number (Ra_C), there existed a minimum value of Reynolds number (Re), below which R-T type of instability appeared. Above this minimum value, the instability was due to two buoyancy forces, known as buoyant instability. Simulations of secondary flow via energy analysis demonstrated the development of complex dynamics at the critical state in all three regimes characterized by transition of multi to uni-cellular structures and vice verse.     

Shear flow instability of rotating hydromagnetic fluids

Abstract

The effect of an axial magnetic field on the linear stability of shear flows in rotating systems is examined by extending Busse's analysis of the nonmagnetic case to fluids of high magnetic diffusivity in the presence of a magnetic field. The shear is caused by differential rotation which creates slight deviations from a state of rigid rotation, corresponding to a small Rossby number. It is found that the Rossby number for the onset of instability is larger when a magnetic field is present than when it is absent.     

The nonlinear stabilization of a zonal shear flow instability

Abstract

The development of initially small perturbations in a weakly supercritical zonal shear flow on a β-plane is studied. Two different scenarios of evolution are possible. If the supercriticality is sufficiently small, the growth of a perturbation is stopped in the viscous critical layer regime; for this case the evolution equation (corrected by the inclusion of a quintic nonlinearity) is derived. At greater supercriticality the nonlinearity cannot stop the growth of the perturbation in a linear (viscous or unsteady) critical layer regime, and the evolution is more complicated. Transition to a nonlinear critical layer regime leads to a reduction in the growth rate and to a slowing (but not a stopping) of the increase in amplitude, A. These are connected to the formation of a plateau (S=constant) of width L=O(A ^½) in the profile of absolute vorticity, S. Careful analysis reveals that the growth in amplitude ceases only when the whole instability domain (where the slope of unperturbed S-profile is positive) becomes covered again by the plateau.     

The structure of separated flow past a circular cylinder in a rotating frame

Abstract

This paper examines the detailed E ^1/4-layer structure of separated flow past a circular cylinder in a low-Rossby-number rotating fluid as the Ekman number E tends to zero. This structure is based on an initial proposal by Page (1987) but with some modifications in response to further evidence, outlined both in this paper and elsewhere, on the behaviour of E ^1/4-layer flows in this context. Numerical calculations for flow in an E ^1/4 shear layer along the separated free streamline are described and the mass flux from this layer is then used to calculate the higher-order flow within the separation bubble. The flow structure is found to have two forms, depending on the value of the O(1) parameter λ, and these are compared with results from published “Navier-Stokes” type calculations for the flow at small but finite values of E.     

Barotropic instability of weakly non-parallel zonal flows

Abstract

Barotropic instability of weakly non-parallel zonal flows with localized intense shear regions is investigated numerically. The numerical integrations of the linear stability problem reveal the existence of unstable localized wave packets whose spatial structure and eigenfrequencies depend on two parameters which measure the degree of supercriticality and the zonal length-scale of the shear region. The results indicate that the structure of the instability is determined by conditions that ensure the decay of the wave packet at infinity and the transition from long to short waves across a turning point (critical layer) region which is controlled by non-parallel effects. The controlling influence exerted by the weak non-parallel effects on the evolution of the instability underlines the weakness of the parallel flow assumption which can be used locally, away from critical layers, as a diagnostic tool only.     

A laboratory study of baroclinic instability

Abstract

Experimental and theoretical results are presented for a simple system which exhibits baroclinic instability. We consider the motion of two immiscible fluids with densities ρ ₁ and ρ ₂ contained in a cylinder rotating with angular frequency ω. The motion is driven by a contact lid rotating with frequency ω + ω. In this paper ω, ω, 2(ρ ₂ – ρ ₁)/(ρ ₂ + ρ ₁), and the geometry are such that the interface does not intersect the “ground” (e.g. an almost horizontal boundary). The motions are described by two-layer quasi-geostrophic equations which are identical, except perhaps for the presence of interfacial friction and tension, with those used in meteorology and oceanography. For small enough internal Froude number F = 4ω² L ²/(g(Δρ/ρ)H) or small enough Rossby number ? = ω/2ω the flow is steady and axisymmetric, the velocity field in each layer being determined primarily by frictional effects in top, bottom, and interfacial Ekman layers. For certain (F, ?) the flow becomes non-axisymmetric. The transition points for the case where the basic potential vorticity gradient is due to interface slope alone have been carefully measured and are in very good agreement with a linear instability theory which neglects sidewall effects. Some preliminary observations of supercritical motion, which include repeatable amplitude and wavenumber vacillation, are reported.     

The instability of precessing flow

Abstract

An explanation is put forward for the instability observed within a precessing, rotating spheroidal container. The constant vorticity solution for the flow suggested by Poincaré is found to be inertially unstable through the parametric coupling of two inertial waves by the underlying constant strain field. Such resonant couplings are due either to the elliptical or shearing strains present which elliptically distort the circular streamlines and shear their centres respectively. For the precessing Earth's outer core, the shearing of the streamlines and the ensuing shearing instability are the dominant features. The instability of some exact, linear solutions for finite precessional rates is established and used to corroborate the asymptotic analysis. A complementary unbounded analysis of a precessing, rotating fluid is also presented and used to deduce a likely upperbound on the growth rate of a small disturbance. Connection is made with past experimental studies.