A numerical study of rotating annulus flows using a modified Galerkin method

Summary Physical phenomena fundamental to rotating, baroclinically driven flows are studied with reference to results of numerical simulation of rotating annulus flows, using a modified Galerkin Model. Both local and global effects of sources, sinks, and transports of heat and momentum are discussed. A convenient energy exchange diagram reveals detailed information that is used to analyze nonlinear equilibration and amplitude vacillation of quasi-geostrophic baroclinic eddies. Transient inertial oscillations, sidewall boundary layers, and internal boundary layers are also discussed.A detailed study of symmetric flows is made, eleven of which are tested numerically for stability with respect to three-dimensional disturbances of a given zonal wave number. Two of the four unstable cases are integrated to a numerical steady state with finite-amplitude, quasi-geostrophic baroclinic waves. With the rigid-lid geometry assumed, the average zonal velocity is zero, resulting in zero phase velocity of the waves. The structure of the thermal wave is nearly coherent in the vertical. These numerical results are consistent with laboratory observations.The eddy flow is quasi-geostrophic except in horizontal boundary layers, where the flow is driven toward low pressure. A small cross-isotherm advection is sufficient to maintain the temperature wave against diffusion and vertical advection. The eddy flow adjusts spontaneously toward the form of the fastest growing or slowest decaying disturbance representable by the truncated space resolution. The eddy flow feeds energy into the mean zonal flow in barotropic-type interactions reflected mainly by the familiar tilted trough. During equilibration, the eddy flow alters the mean zonal flow in such a way that eddy energy sources are reduced relative to energy sinks. However, this adjustment is small compared to the change of total flow, which reflects a relatively large change of eddy amplitude. This suggests that small errors in the mean zonal flow representation can lead to relatively large errors in total flow representation.In most flows studied the kinetic energy dissipation is concentrated in thin boundary layers. In spite of this thinness, the basically laminar character of these dissipative boundary layers allows accurate and economical numerical simulation through the use of characteristic functions, which is a natural refinernent of the basic Galerkin method used. (In this prototype study, only moderately characteristic functions are used, thus sacrificing numerical economy while simplifying the programming.) Similarly, the generation of potential energy, which is transformed into the kinetic energy of the flow, is accurately simulated. In most cases studied, this generation is also concentrated in thin boundary layers where thermal energy is extracted from cold fluid and added to warm fluid.Contribution number 76 of the Geophysical Fluid Dynamics Institute, Florida State University, USA.     

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.     

Two-dimensional numerical modeling of dam-break flows over natural terrain using a central explicit scheme

A two-dimensional (2D) numerical model has been developed to solve shallow water equations for simulation of dam-break flows. The spatial derivatives are discretized using a well-balanced explicit central upwind conservative scheme. The scheme is Riemann solver free and guarantees the positivity of the flow depth over complex topography if the Courant number is kept less than 0.25. The time integration is performed by Euler’s scheme. The model is verified against analytical results for water surface elevation and discharge for three benchmark test cases. A good agreement between analytical solutions and computed results is observed. The property of well-balancing in still water over an uneven bottom is also confirmed. The model is then validated by simulating a laboratory experiment in which a dam break flow propagates over a triangular obstacle. The model performance was found to be satisfactory. A dam break laboratory experimental test case on a frictionless horizontal bottom is also simulated for 2D validation of the model, and good agreement between simulation and the experimental data is observed. The suitability of the proposed model for real life applications is demonstrated by simulating the Malpasset dam-break event, which occurred in 1959 in France. The computed arrival time of the flood wave front and the maximum flow depths at various observation points matched well with the measurements on a 1/400 scale physical model. The overall performance indicates that this model can be applied for simulation of dam-break waves in real life cases.     

Processes of dike-break induced flows:A combined experimental and numerical model study

Dike breaking is a disaster that could cause extensive damage. It could lead to flood flows outside the dike and induce water level fluctuations in the main channel. Numerical models are increasingly used to simulate flood flows due to dike-break, because direct observations from field surveys and physical models are rather limited. Existing knowledge concerning dam-break flows cannot be applied directly to dike-break flows because the effect of channel discharge cannot be neglected in the latter. In this study,physical experiments are done in a large laboratory flume to simulate the process of dike-break induced flood wave propagation in the floodplain and flow fluctuations in the main channel. The variations of water levels and velocities are measured and recorded using an array of pressure sensors and two acoustic Doppler velocimetry devices. A numerical model has been set up according to the experimental layout. The experiments have high repeatability and the numerical model predictions agree closely with the physical model data. The experimental results provide reliable information for improving the understanding of dike-break flow dynamics and for the verification of numerical models.     

Nonstationary pearl pulsations as a signature of magnetospheric disturbances

We analyse long-lasting (several hours) Pc1 pearl pulsations with decreasing, increasing or constant central frequencies. We show that nonstationary pearl events (those with either decreasing or increasing central frequency) are observed simultaneously with increasing auroral magnetic activity at the nightside magnetosphere while the stationary events (constant central frequency) correspond to quiet magnetic conditions. Events with decreasing central frequency are observed mostly in the late morning and daytime whereas events with increasing central frequency appear either early in the morning or in the afternoon. We explain the diurnal distribution of the nonstationary pearl pulsations in terms of proton drifts depending on magnetic activity, and evaluate the magnetospheric electric field based on the variation of the central frequency of pearl pulsations.     

Bubble coalescence in basalt flows: comparison of a numerical model with natural examples

Gas accumulation in magma may be aided by coalescence of bubbles because large coalesced bubbles rise faster than small bubbles. The observed size distribution of gas bubbles (vesicles) in lava flows supports the concept of post-eruptive coalescence. A numerical model predicts the effects of rise and coalescence consistent with observed features. The model uses given values for flow thickness, viscosity, volume percentage of gas bubbles, and an initial size distribution of bubbles together with a gravitational collection kernel to numerically integrate the stochastic collection equation and thereby compute a new size spectrum of bubbles after each time increment of conductive cooling of the flow. Bubbles rise and coalesce within a fluid interior sandwiched between fronts of solidification that advance inward with time from top and bottom. Bubbles that are overtaken by the solidification fronts cease to migrate. The model predicts the formation of upper and lower vesicle-rich zones separated by a vesicle-poor interior. The upper zone is broader, more vesicular, and has larger bubbles than the lower zone. Basaltic lava flows in northern California exhibit the predicted zonation of vesicularity and size distribution of vesicles as determined by an impregnation technique. In particular, the size distribution at the tops and bottoms of flows is essentially the same as the initial distribution, reflecting the rapid initial solidification at the bases and tops of the flows. Many large vesicles are present in the upper vesicular zones, consistent with expected formation as a result of bubble coalescence during solidification of the lava flows. Both the rocks and model show a bimodal or trimodal size distribution for the upper vesicular zone. This polymodality is explained by preferential coalescence of larger bubbles with subequal sizes. Vesicularity and vesicle size distribution are sensitive to atmospheric pressure because bubbles expand as they decompress during rise through the flow. The ratio of vesicularity in the upper to that in the lower part of a flow therefore depends not only on bubble rise and coalescence, but also on flow thickness and atmospheric pressure. Application of simple theory to the natural basalts suggests solidification of the basalts at 1.0±0.2 atm, consistent with the present atmospheric pressure. Paleobathymetry and paleoaltimetry are possible in view of the sensitivity of vesicle size distributions to atmospheric pressure. Thus, vesicular lava flows can be used to crudely estimate ancient elevations and/or sea level air pressure.     

下击暴流风的数值仿真研究

作为自然界的一种强风荷载,雷暴天气产生的下击暴流风在世界范围内造成了大量工程结构物的破坏。研究证实,下击暴流风具有同常规的大气边界层近地风完全不同的风场特征。利用通用商业CFD(Computational fluid dynamics)软件FLUENT完成了几何缩尺比为1:2000的下击暴流过程的数值仿真。非定常的数值模拟将有助于更好地理解下击暴流风的动力特征。数值结果表明,基于雷诺平均的Navier-Stokes方程进行的非定常数值模拟结果同根据实测数据提出的三种经验风剖面吻合良好,在目前实测资料稀少的情况下,有助于更好地研究下击暴流风的基本特征。     

A numerical study of the string function using a primitive equation ocean model

Abstract

We use results from a primitive-equation ocean numerical model (SCRUM) to test a theoretical 'string function' formulation put forward by Tyler and Käse in another article in this issue. The string function acts as a stream function for the large-scale potential energy flow under the combined beta and topographic effects. The model results verify that large-scale anomalies propagate along the string function contours with a speed correctly given by the cross-string gradient. For anomalies having a scale similar to the Rossby radius, material rates of change in the layer mass following the string velocity are balanced by material rates of change in relative vorticity following the flow velocity. It is shown that large-amplitude anomalies can be generated when wind stress is resonant with the string function configuration.     

Interactions between river flows and colonizing vegetation on a braided river: exploring spatial and temporal dynamics in riparian vegetation cover using satellite data

Field, laboratory, and numerical modelling research are increasingly demonstrating the potential of riparian tree colonization and growth to influence fluvial dynamics and the evolution of fluvial landforms. This paper jointly analyses multi‐temporal, multispectral ASTER data, continuous river stage and discharge data, and field observations of the growth rates of the dominant riparian tree species (Populus nigra) along a 21 km reach of the Tagliamento River, Italy. Research focuses on the period 2004–2009, during which there was a bankfull flood on 24 October 2004, followed by 2 years with low water levels, nearly 2 years with only modest flow pulses, and then a final period from 15 August 2008 that included several intermediate to bankfull flow events. This study period of increasing flow disturbance allows the exploration of vegetation dynamics within the river's active corridor under changing flow conditions. The analysis demonstrates the utility of ASTER data for investigating vegetation dynamics along large fluvial corridors and reveals both spatial and temporal variations in the expansion, coalescence, and erosion of vegetated patches within the study reach. Changes in the extent of the vegetated area and its dynamics vary along the study reach. In sub‐reaches where riparian tree growth is vigorous, the vegetated area expands rapidly during time periods without channel‐shaping flows, and is subsequently able to resist erosion by bankfull floods. In contrast, in sub‐reaches where tree growth is less vigorous, the vegetated area expands at a slower rate and is more readily re‐set by bankfull flood events. This illustrates that the rate of growth of riparian trees is crucial to their ability to contribute actively to river corridor dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

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 α².     

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     

Analytical study of the transverse distribution of along-channel and transverse residual flows in tidal estuaries

We present an analytical model to decompose complex along-channel and transverse residual flows into components induced by individual mechanisms. The model describes the transverse distribution of residual flows in tidally dominated estuaries. Scaling and perturbation techniques are used to obtain analytical solutions for residual flows over arbitrary across-channel bed profiles. The flows are induced by horizontal density gradients, tidal rectification processes, river discharge, wind, channel curvature and the earth's rotation. These rectification processes induce residual flows that are up-estuary to the right and down-estuary to the left of an estuarine channel (looking up-estuary in the northern hemisphere). The tidal rectification processes fundamentally change the transverse structure of along-channel residual flows in many tidal estuaries, as these processes cause the flows to be internally asymmetric about the mid-axis of the channel for relatively large tidal velocities, steep channels or narrow estuaries. In addition, velocity scales are derived from the analytical solutions to estimate the relative importance of the various residual flow mechanisms from estuarine parameters. A case study of a transect across the Upper Chesapeake Bay showed that important features of the residual flow observed in that transect are reproduced and explained by the analytical model. The velocity scales were able to identify the relevant residual flow mechanisms as well. The tidal rectification processes considered here result from advection of along-channel tidal momentum by Coriolis-induced transverse tidal currents.     

Water storage dynamics in a till hillslope: the foundation for modeling flows and turnover times

Studies on hydrology, biogeochemistry, or mineral weathering often rely on assumptions about flow paths, water storage dynamics, and transit times. Testing these assumptions requires detailed hydrometric data that are usually unavailable at the catchment scale. Hillslope studies provide an alternative for obtaining a better understanding, but even on such well‐defined and delimited scales, it is rare to have a comprehensive set of hydrometric observations from the water divide down to the stream that can constrain efforts to quantify water storage, movement, and turnover time. Here, we quantified water storage with daily resolution in a hillslope during the course of almost an entire year using hydrological measurements at the study site and an extended version of the vertical equilibrium model. We used an exponential function to simulate the relationship between hillslope discharge and water table; this was used to derive transmissivity profiles along the hillslope and map mean pore water velocities in the saturated zone. Based on the transmissivity profiles, the soil layer transmitting 99% of lateral flow to the stream had a depth that ranged from 8.9 m at the water divide to under 1 m closer to the stream. During the study period, the total storage of this layer varied from 1189 to 1485 mm, resulting in a turnover time of 2172 days. From the pore water velocities, we mapped the time it would take a water particle situated at any point of the saturated zone anywhere along the hillslope to exit as runoff. Our calculations point to the strengths as well as limitations of simple hydrometric data for inferring hydrological properties and water travel times in the subsurface.     

Baroclinic dynamics of wind-driven circulation in a stratified bay: A numerical study using models of varying complexity

The baroclinic response of a stratified coastal embayment (Lunenburg Bay of Nova Scotia) to the observed wind forcing is examined using two numerical models. A linear baroclinic model based on the normal mode approach shows skill at reproducing the observed isotherm movements and sub-surface currents during a time of strong stratification in the bay. The linear model also shows that the isotherm movement in Lunenburg Bay is influenced by the wind forcing and propagation of baroclinic Kelvin waves from neighbouring Mahone Bay. The effects of nonlinearity and topography are investigated using a three-dimensional nonlinear coastal circulation model. The nonlinear model results demonstrate that the nonlinear advection terms generate a gyre circulation at the entrance of Lunenburg Bay, and the slope bottom topography at the mouth of the bay strengthens the sub-surface time-mean inflow on the southern side of the bay. A comparison of model-calculated currents in different numerical experiments clearly shows that baroclinicity plays a dominant role in the dynamics of wind-driven circulation in Lunenburg Bay.     

A quasi-geostrophic numerical model of a rotating internally heated fluid

Abstract

A quasi-geostrophic numerical model of flow in a rotating channel is integrated under conditions typical of laboratory experiments with an internally heated annulus system. Compared to a laboratory experiment, or a full Navier-Stokes simulation, the quasi geostrophic numerical model is a simple system. It includes nonlinear interactions, dissipation via conventional parameterizations of Ekman layers and internal diffusion, and a steady forcing term which represents heating near the centre of the channel and cooling near both sides. Explicit boundary layers, cylindrical geometry effects, horizontal variations in static stability and variations in conductivity and diffusivity with temperature are all absent, and ageostrophic advection is incompletely represented. Nevertheless, over a range of parameters, flows are produced which strongly resemble those seen in the laboratory thus suggesting that the most important physical processes are represented. The numerical model is used to map out a regime diagram which includes examples of steady flows, flows with periodic time dependence (wavenumber vacillations) and flows which are irregularly time dependent.     

A slow mode wave as a possible source of Pi 2 and associated particle precipitation: a case study

An intensification of auroral luminosity referred to as an auroral break-up often accompanies the onset of geomagnetic pulsation (Pi 2) at the dip-equator. One such auroral break-up occurred at 2239 UT on 16 June, 1986, being accompanied by weak substorm activity (AE≈50 nT) which was recorded in all-sky image of Syowa Station, Antarctica (66.2°S, 71.8°E in geomagnetic coordinates). The associated Pi 2 magnetic pulsation was detected by a fluxgate magnetometer in the afternoon sector at the dip-equator (Huancayo, Peru; 1.44°N, 355.9° in geomagnetic coordinates; 12.1°S, 75.2°W in geographic coordinates; L = 1.00). In spite of the large separation of the two stations in longitude and latitude, the auroral break-up and subsequent luminosity modulation were seen to be correlated with the wave form of the ground Pi 2 pulsation. This occurred in such a way that the luminosity maximum was seen to occur at the phase of maximum amplitudes of Pi 2 wave form. We argue that the observed correlation could be interpreted as indicating a Pi 2-modulation of a field-aligned acceleration of the low energy electrons that may occur near the equator of the midnight magnetosphere.     

Effective velocities in fractured media: a numerical study using the rotated staggered finite‐difference grid

The modelling of elastic waves in fractured media with an explicit finite‐difference scheme causes instability problems on a staggered grid when the medium possesses high‐contrast discontinuities (strong heterogeneities). For the present study we apply the rotated staggered grid. Using this modified grid it is possible to simulate the propagation of elastic waves in a 2D or 3D medium containing cracks, pores or free surfaces without hard‐coded boundary conditions. Therefore it allows an efficient and precise numerical study of effective velocities in fractured structures. We model the propagation of plane waves through a set of different, randomly cracked media. In these numerical experiments we vary the wavelength of the plane waves, the crack porosity and the crack density. The synthetic results are compared with several static theories that predict the effective P‐ and S‐wave velocities in fractured materials in the long wavelength limit. For randomly distributed and randomly orientated, rectilinear, non‐intersecting, thin, dry cracks, the numerical simulations of velocities of P‐, SV‐ and SH‐waves are in excellent agreement with the results of the modified (or differential) self‐consistent theory. On the other hand for intersecting cracks, the critical crack‐density (porosity) concept must be taken into account. To describe the wave velocities in media with intersecting cracks, we propose introducing the critical crack‐density concept into the modified self‐consistent theory. Numerical simulations show that this new formulation predicts effective elastic properties accurately for such a case.     

Parametric decay of beam-driven Langmuir wave and enhanced ion-acoustic fluctuations in the ionosphere: a weak turbulence approach

We present a model that describes the decay of beam generated Langmuir waves into ion-acoustic waves in the topside ionosphere. This calculation is done within the frame of the weak turbulence approximation. We study the spectral signature of such a process as seen by a VHF incoherent scatter radar. An incoherent scatter (IS) spectrum is characterized by two maxima at k_radar and −k_radar, the right and left ion lines respectively. It is shown that, for reasonable beam parameters, the parametric decay of beam-generated Langmuir waves can enhance either the right, the left or both ion lines simultaneously. The shape of the spectrum can change drastically on time scale of about 0.1 to 1 s. The role of the beam parameter as well as the ionospheric parameters is also investigated. For a given beam number density, the beam energy or the background density are important to trigger either the left or the right ion line. A large energy spread of the beam or low electron collision frequencies can explain the simultaneous observations of the left and the right ion line. The importance of the electron collision frequency can explain the altitude distribution of the coherent echoes observed by incoherent scatter radars.     

Statistical evaluation of compositional changes in volcanic gas chemistry: a case study

The geochemical analysis of fumarolic gases collected at quiescent and active volcanic systems over time is one of the main tools to understand changes in the state of activity for surveillance and risk assessment. The continuous output of chemical species through fumarolic activity, which characterizes the inter-eruptive intervals, has also a major and general influence on the environment. The mobilization of chemical species due to weathering of volcanic rocks, or the input of gaseous components from fumarolic activity, results in some kind of modification of the environment affecting, in particular, water, soils, and the consequent growth of the plants present in these areas. In this paper, an investigation on the chemical composition of fumarolic gases collected at Vulcano island (Sicily, southern Italy) is performed, with the aim to discover how data changes during the monitored period of time, and to design a strategy for the environmental surveillance of volcanic systems taking into account the nature of the analyzed data. In order to summarize the contribution of all the components that can affect the chemical composition of volcanic gases, a multivariate statistical approach appears to be suitable. Since many of those methods assume independent observations, the possible presence of time-dependent structures should be carefully verified. In this framework, given the compositional nature of geochemical data, we have applied recent theoretical and practical developments in the field of compositional data analysis to work in the correct sample space and to isolate groups of parts responsible for significant changes in the gas chemistry. The proposed approach can be generalized to the investigation of complex environmental systems.     

Contemporary hydrodynamics and morphological change of a microtidal estuary: a numerical modelling study

Contemporary hydrodynamics and morphological change are examined in a shallow microtidal estuary, located on a wave-dominated coast (Port Stephens, NSW, Australia). Process-based numerical modelling is undertaken by combining modules for hydrodynamics, waves, sediment transport and bathymetry updates. Model results suggest that the complex estuarine bathymetry and geometry give rise to spatial variations in the tidal currents and a marked asymmetry between ebb and flood flows. Sediment transport paths correspond with tidal asymmetry patterns. The SE storms significantly enhance the quantities of sediment transport, while locally generated waves by the westerly strong winds also are capable of causing sediment entrainment and contribute to the delta morphological change. The wave/wind-induced currents are not uniform with flow over shoals driven in the same direction as waves/winds while a reverse flow occurring in the adjacent channel. The conceptual sediment transport model developed in this study shows flood-directed transport occurs on the flood ramp while ebb-directed net transport occurs in the tidal channels and at the estuary entrance. Accretion of the intertidal sand shoals and deepening of tidal channels, as revealed by the model, suggest that sediment-infilling becomes advanced, which may lead to an ebb-dominated estuary. It is likely that a switch from flood- to ebb-dominance occurs during the estuary evolution, and the present-day estuary acts as a sediment source rather than sediment sink to the coastal system. This is conflictive to the expectation drawn from the estuarine morphology; however, it is consistent with previous research suggesting that, in an infilling estuary, an increase in build-up of intertidal flats/shoals can eventually shift an estuary towards ebb dominance. Thus, field data are needed to validate the result presented here, and further study is required to investigate a variety of estuaries in the Australian area.