Barotropic unstable modes in zonal and meridional channel on the beta-plane

Abstract

The linear stability of a non-divergent barotropic parallel shear flow in a zonal and a non-zonal channel on the β plane was examined numerically. When the channel is non-zonal, the governing equation is slightly modified from the Orr-Sommerfeld equation. Numerical solutions were obtained by solving the discretized linear perturbation equation as an eigenvalue problem of a matrix. When the channel is zonal and lateral viscosity is neglected the problem is reduced to the ordinary barotropic instability problem described by Kuo's (1949) equation. The discrepancy between the stability properties of westward and eastward flows, which have been indicated by earlier studies, was reconfirmed. It has also been suggested that the unstable modes are closely related to the continuous modes discretized by a finite differential approximation. When the channel is non-zonal, the properties of unstable modes were quite different from those of the zonal problem in that: (1) The phase speed of the unstable modes can exceed the maximum value of the basic flow speed; (2) The unstable modes are not accompanied by their conjugate mode; and (3) The basic flow without an inflection point can be unstable. The dispersion relation and the spatial structure of the unstable modes suggested that, irrespective of the orientation of the channel, they have close relation to the neutral modes (Rossby channel modes) which are the solutions in the absence of a basic shear flow. The features mentioned above are not dependent on whether or not the flow velocity at the boundary is zero.     

Nonlinear dynamical modeling of solar cycles using dynamo formulation with turbulent magnetic helicity

Abstract

We describe a sequence of two-dimensional numerical simulations of inflection point instability in a stably stratified shear flow near the ground. The fastest growing Kelvin-Helmholtz modes are studied in detail; in particular we investigate the growth inhibiting effect of the ground which is predicted by linear theory and the Reynolds number dependence of the process of growth to finite amplitude. We consider flows which are both above and below the critical Reynolds number (Re = 300) which has been reported by Woods (1969) to mark the boundary between flows which have turbulent final states and those which do not. A global energy budget reveals a fundamental difference in character of the finite amplitude billows in these two Reynolds number regimes. However, for relatively high Reynolds numbers (Re = 10³) we do not find any explicit evidence for secondary instability. Above the transition Reynolds number the modified mean flow induced by wave growth is characterized by a splitting of the original shear layer and of the in version in which it is embedded.     

Improved bounds on linear instability of barotropic zonal flow within the shallow water equations

Here we develop mathematical results to describe the location of linear instability of a parallel mean flow within the framework of the shallow water equations; growth estimates of near neutral modes (for disturbances subcritical with respect to gravity wave speed) in the cases of non-rotating and rotating shallow water. The bottom topography is taken to be one-dimensional and the isobaths are parallel to the mean flow. In the case of a rotating fluid, the isobaths and the mean flow are assumed to be zonal. The flow is front-like: there is a monotonic increase of mean flow velocity. Our results show that for barotropic flows the location of instabilities will be a semi-ellipse region in the complex wave velocity plane, that is based on the wave-number, Froude number, and depth of the fluid layer. We also explore the instability region for the case of spatially unbounded mean velocity profiles for non-rotating shallow water.     

A Note on Standing Internal Inertial Gravity Waves of Finite Amplitude

The effects of finite amplitude are examined in two-dimensional, standing, internal gravity waves in a rectangular container which rotates about a vertical axis at frequency f/ 2. Expressions are given for the velocity components, density fluctuations and isopycnal displacements to second order in the wave steepness in fluids with buoyancy frequency, N , of general form, and the effect of finite amplitude on wave frequency is given in an expansion to third order. The first order solutions, and the solutions to second order in the absence of rotation, are shown to conserve energy during a wave cycle. Analytical solutions are found to second order for the first two modes in a deep fluid with N proportional to sech( az ), where z is the upward vertical coordinate and a is scaling factor. In the absence of rotation, results for the first mode in the latter stratification are found to be consistent with those for interfacial waves. An analytical solution to fourth order in a fluid with constant N is given and used to examine the effects of rotation on the development of static instability or of conditions in which shear instability may occur. As in progressive internal waves, an effect of rotation is to enhance the possibility of shear instability for waves with frequencies close to f . The analysis points to a significant difference between the dynamics of standing waves in containers of limited size and progressive internal waves in an unlimited fluid; the effect of boundaries on standing waves may inhibit the onset of instability. A possible application of the analysis is to transverse oscillations in long, narrow, steep-sided lakes such as Loch Ness, Scotland.     

The instability of internal gravity waves to localised disturbances

The instability of an internal gravity wave due to nonlinear wave-wave interaction is studied theoretically and numerically. Three different aspects of this phenomenon are examined. 1. The influence of dissipation on both the resonant and the nonresonant interactions is analysed using a normal mode expansion of the basic equations. In particular, the modifications induced in the interaction domain are calculated and as a result some modes are shown to be destabilised by dissipation. 2. The evolution of an initial unstable disturbance of finite vertical extent is described as the growth of two secondary wave packets travelling at the same group velocity. A quasi-linear correction to the basic primary wave is calculated, corresponding to a localised amplitude decrease due to the disturbance growth. 3. Numerical experiments are carried out to study the effect of a basic shear on wave instability. It appears that the growing secondary waves can have a frequency larger than that of the primary wave, provided that the shear is sufficient. The instability of waves with large amplitude and long period, such as tides or planetary waves, could therefore be invoked as a possible mechanism for the generation of gravity waves with shorter period in the middle atmosphere.     

切变基本纬向流中非线性赤道Rossby长波

为了解决观测和理论研究中的一些问题以及更好地了解热带大气动力学 ,有必要进一步研究基本气流的变化对大气中赤道Rossby波动的影响 .本文研究分析基本气流对赤道Rossby长波的影响 ,利用一个简单赤道 β平面浅水模式和摄动法 ,研究纬向基本气流切变中非线性赤道Rossby波 ,推导出在切变基本纬向流中赤道Rossby长波振幅演变所满足的非线性KdV方程并得到其孤立波解 .分析表明 ,孤立波存在的必要条件是基本气流有切变 ,而且基流切变不能太强 ,否则将产生正压不稳定 .     

A model of the wave and current boundary‐layer structure

A model of the wave and current boundary‐layer structure was developed using the k–ε turbulent closure model. The finite‐difference method was used to solve the governing equations. Vertical logarithmic grids and equal time steps were adopted. The following modelled simulations were obtained: (1) vertical profiles of wave velocity amplitude, eddy viscosity coefficient and turbulent kinetic energy with waves only; (2) vertical profiles of wave velocity amplitude, mean current velocity, eddy viscosity coefficient and turbulent kinetic energy with waves having a following current. To test the validity and the rationality of the present model, vertical profiles of modelled wave velocity amplitude and mean velocity were compared with corresponding experimental results available in the literature. Copyright © 2006 John Wiley & Sons, Ltd.     

Accelerating shear velocity in gravel-bed channels

Abstract

The behaviour of the shear velocity along a gravel-bed channel is investigated experimentally in the presence of a negative pressure gradient (accelerating flow). Different methods of estimation of the shear velocity, derived from vertical profiles of the mean longitudinal point velocity, are examined and a new method is proposed. Results show that the proposed method of estimation is comparable to the St Venant and Clauser's methods. At a specific cross section, for constant bottom slope and relative roughness, shear velocity increases with discharge.     

切变基本纬向流中非线性赤道Rossby长波

为了解决观测和理论研究中的一些问题以及更好地了解热带大气动力学 ,有必要进一步研究基本气流的变化对大气中赤道Rossby波动的影响 .本文研究分析基本气流对赤道Rossby长波的影响 ,利用一个简单赤道 β平面浅水模式和摄动法 ,研究纬向基本气流切变中非线性赤道Rossby波 ,推导出在切变基本纬向流中赤道Rossby长波振幅演变所满足的非线性KdV方程并得到其孤立波解 .分析表明 ,孤立波存在的必要条件是基本气流有切变 ,而且基流切变不能太强 ,否则将产生正压不稳定 .     

Instability of a thin magnetic tube in the solar atmosphere

Abstract

The static equilibrium of a thin vertical magnetic tube embedded in the solar atmosphere is shown to be dynamically unstable against the fundamental mode of perturbation having no nodes in the vertical displacement. The instability has its origin in the convection zone, and the eigenfunction is extended further up in the stable upper layers by the magnetic field which guides the displacement mainly in the longitudinal direction. It is suggested that the downdraft observed in the solar network structure is a finite amplitude consequence of this instability. The overtone modes are found to be stable.     

正应力扰动对断层滑动失稳影响的实验研究

利用双轴伺服控制加载装置,采用三块花岗闪长岩标本组成的含有两个滑动面的直剪结构,开展了摩擦滑动实验.实验中通过在垂直滑动面的载荷上叠加正弦波状和方波状的扰动,研究了正应力扰动对断层黏滑失稳的影响.研究表明,在恒定的正应力和位移速率下,标本表现为规则的黏滑,叠加正应力扰动后,随扰动振幅的增加黏滑发生时间与扰动的相关性增大,黏滑应力降和时间间隔的分布趋于离散.黏滑应力降和时间间隔的平均值随平均正应力的增加呈线性增长,扰动叠加后黏滑应力降的离散度随平均正应力的增加而增大;黏滑应力降和时间间隔主要受应力变化幅度的影响,而与应力变化的速率关系不大.剪应力和正应力扰动都会对断层黏滑失稳产生影响,而正应力扰动的影响更明显.这两种扰动对断层黏滑失稳影响的机制存在差异,剪应力扰动只是改变断层滑动的推动力,而正应力扰动则改变了断层面上凹凸体的接触状态.     

On the propagation of pressure-patterns: Part I

Summary The propagation speed of sinoidal troughs and wedges in a steady state flow is determined from consideration of the mass transport due to the bodily motion of the system. Fundamental propositions are established regarding the mutual motion of wind-, pressure-, temperature-, and density-fields.It is found that in a frictionless barotropic general flow, all perturbations are propagated with the same speed—the speed of the general current. In a baroclinic general flow a perturbation will only be propagated without dispersion if it has a specific (sinoidal) horizontal structure.When a sinoidal perturbation is embedded in a baroclinic general flow-field, it will be propagated as though by a barotropic flow with the sameeffective speed. The effective speed can be computed when the vertical structure of the perturbation and of the mean flow are known.It is frequently assumed that the speed of mean flow at some particular level (500 mb is often assumed) gives the «steering» of the surface perturbation by a baroclinic general flow, that is to say, a baroclinic flow steers a perturbation with the speed of an equivalent barotropic field. The present paper provides a rational basis for the concept of an equivalent barotropic flow, but it is to be remembered that the «steering level» does not depend uniquely on the vertical structure of the mean flow-field, but varies from perturbation to perturbation, being lower for shallow perturbations than for (vertically) deep ones.     

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.     

Barotropic local instability and severe storm process

By means of barotropic model, the characteristic and initial value problems are investigated to reveal the local two-dimensional barotropic instability of the nonuniform current to the dynamic mechanism of the formation of the Yangtze-Huaihe River severe storm in July 1991. Analytical theory and numerical experiment show that (i) the unstable developing modes are chiefly the two periods of about 44 d and 10 d, which are fundamentally consistent with that of the precipitation change of the Yangtze-Huaihe River. (ii) The growth rate of the local perturbation is dominated by the meridional wave numbern = 1–5 and zonal wave numberk = 1–12, i.e. the severe storm over the Yangtze-Huaihe River results from the interaction of the systems at different latitudes and waves of different scales, (iii) The perturbation over the Yangtze-Huaihe River possesses the property of local intensification, which slowly migrates westward over the lower and middle reaches of the Yangtze-Huaihe River. (iv) The growth rate of the instability, especially the propagation velocity of the perturbation, is sensitive to the external parameters ū and α. Project supported by the National Natural Science Foundation of China.     

斜压基本气流对东亚夏季风区气旋扰动低频发展影响的数值模拟研究

本文将夏季气候平均的基本气流分解为正压和斜压分量,使用一个线性斜压模式,研究了不同斜压基本气流对热带西北太平洋地区初始气旋性环流扰动低频发展演变的重要作用.其中,控制试验较好地模拟出初始气旋扰动向西北方向传播、在西北太平洋季风槽附近停滞增强、在东亚地区出现经向波列和在南海到海洋大陆地区形成西北—东南向波列等特征.改变斜压分量的敏感性试验结果表明,正压基流不能为西传的初始扰动供给足够的能量;海陆热力差异引起东亚地区的纬向温度梯度和北风垂直切变,是东亚太平洋型经向波列形成和维持的重要因素;当基本气流中的斜压纬向偏差部分线性增大时,扰动的能量会呈e指数迅速增强,提示在气候变化的背景下,基本气流微小的改变可能带来天气或季节内扰动强度的剧烈响应.     

Instability of planetary waves in a high vertical resolution model

Abstract

A high vertical resolution model is used to examine the instability of a baroclinic zonal flow and a finite amplitude topographically forced wave. Two families of unstable modes are found, consisting of zonally propagating most unstable modes, and stationary unstable modes. The former have time scale and spatial structure similar to baroclinic synoptic disturbances, but are localized in space due to interaction with the zonally asymmetric forcing. These modes transport heat efficiently in both the zonal and meridional directions. The second family of stationary unstable modes has characteristics of modes of low frequency variability of the atmosphere. They have time scales of 10 days and longer, and are of planetary scale with an equivalent barotropic vertical structure. The horizontal structure resembles blocking flows. They are maintained by available potential energy of the basic wave, and have large zonal heat fluxes. The results for both families of modes are interpreted in terms of an interaction between forcing and baroclinic instability to create favoured regions for eddy development. Applications to baroclinic planetary waves are also considered.     

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.     

风的垂直切变对滞弹性大气中静力适应过程的影响

本文基于描写滞弹性大气静力适应过程的线性方程组,从波动频散关系、气团运动规律和能量转换的角度出发,研究了水平基流及其垂直切变对该模式大气静力适应过程的影响.构造四种水平基流垂直分布模型进行比较,分别为常数型、线性切变型、反气旋切变型和气旋切变型,得到结论：（1）具有重力波性质的波动是滞弹性大气静力适应过程中扰动能量传播的方式,当垂直折射指数大于零时,基本气流及其垂直切变的存在,不仅改变了波动频率的大小,而且改变波动传播的方向;（2）在静力适应过程中气块的运动轨迹呈椭圆形,水平基流及其垂直切变影响椭圆的扁率,同时也影响扰动物理量之间的偏振关系;（3）水平基流的垂直切变是扰动能量和水平基流能量发生转换的媒介,当存在垂直向上的动量输送时,正的垂直风切变对应扰动能量的衰减,水平基流能量的增加,负的垂直风切变对应扰动能量的增加,水平基流能量的衰减;（4）不同的风的垂直切变模型,对静力适应过程的影响不同;对于水平基流呈反气旋切变型和气旋切变型,扰动发展的波动垂直结构为,急流轴上方波动等相位线自下而上向西倾斜,急流轴下方波动等相位线自下而上向东倾斜,反之亦然.

     

Magnitude of deviatoric terms in vertically averaged flow equations

The depth‐integrated momentum and kinetic energy equations contain velocity correlation terms that involve products of local deviations in velocity components about depth‐averaged values. Based on velocity data obtained from North Boulder Creek, Colorado, a simple scaling analysis suggests that certain of these terms, which normally can be neglected in the case of smooth channels, can be significant parts of the momentum and energy balances in steep, rough channels owing to the occurrence of non‐logarithmic velocity profiles. A linearized version of the kinetic energy equation suggests that, for flow accelerations over small‐amplitude bed forms, the energy of the mean motion is spatially partitioned between a form involving the depth‐averaged velocity and a form involving the deviatoric part of the velocity profile; this partitioning is associated with spatial variations in the uniformity of the vertical profile of the streamwise velocity. These points are consistent with published flume measurements involving flow over sand‐roughened dunes, and with published field measurements of flow over a gravel bar. Copyright © 2000 John Wiley & Sons, Ltd.     

Hydrological characteristics of vegetated river flows: a laboratory flume study

Abstract

Laboratory flume experiments were undertaken to measure the vertical profiles of mean flow velocity for three different flow discharges and four different stem densities of Hydrilla verticillata. The data were used to calculate three parameters, namely Manning's roughness coefficient, the Reynolds number and the Froude number. In addition, empirical equations were obtained for the vertical distribution of measured flow velocity within the transitional zone and above the plant canopy. The results show that: (a) the vertical distribution of measured flow velocity exhibits three zone profiles; (b) Manning's roughness coefficient decreases with increasing depth-averaged flow velocity; (c) the relationship between Manning's roughness coefficient and the depth-averaged flow velocity is within the smooth left inverse curve; (d) Manning's roughness coefficient significantly changes with increasing density of Hydrilla; (e) the Froude number is independent of the density of Hydrilla; and (f) both the Reynolds number and the Froude number increase with increasing depth-averaged flow velocity.

Citation Shi, J.Z., Li, Y.-H., Hughes, J.M.R., and Zhao, M., 2013. Hydrological characteristics of vegetated river flows: a laboratory flume study. Hydrological Sciences Journal, 58 (5), 1047–1058.

Editor Z.W. Kundzewicz