Classical plume theory: 1937-2010 and beyond

Developing a theoretical description of turbulent plumes, the likes of which may be seen rising above industrial chimneys, is a daunting thought. Plumes are ubiquitous on a wide range of scales in both the natural and the man-made environments. Examples that immediately come to mind are the vapour plumes above industrial smoke stacks or the ash plumes forming particle-laden clouds above an erupting volcano. However, plumes also occur where they are less visually apparent, such as the rising stream of warmair above a domestic radiator, of oil from a subsea blowout or, at a larger scale, of air above the so-called urban heat island. In many instances, not only the plume itself is of interest but also its influence on the environment as a whole through the process of entrainment. Zeldovich (1937, The asymptotic laws of freely-ascending convective flows. Zh. Eksp. Teor. Fiz., 7, 1463-1465 (in Russian)), Batchelor (1954, Heat convection and buoyancy effects in fluids. Q. J. R. Meteor. Soc., 80, 339-358) and Morton et al. (1956, Turbulent gravitational convection from maintained and instantaneous sources. Proc. R. Soc. Lond. A, 234, 1-23) laid the foundations for classical plume theory, a theoretical description that is elegant in its simplicity and yet encapsulates the complex turbulent engulfment of ambient fluid into the plume. Testament to the insight and approach developed in these early models of plumes is that the essential theory remains unchanged and is widely applied today. We describe the foundations of plume theory and link the theoretical developments with the measurements made in experiments necessary to close these models before discussing some recent developments in plume theory, including an approach which generalizes results obtained separately for the Boussinesq and the non-Boussinesq plume cases. The theory presented - despite its simplicity - has been very successful at describing and explaining the behaviour of plumes across the wide range of scales they are observed. We present solutions to the coupled set of ordinary differential equations (the plume conservation equations) that Morton et al. (1956) derived from the Navier-Stokes equations which govern fluid motion. In order to describe and contrast the bulk behaviour of rising plumes from general area sources, we present closed-form solutions to the plume conservation equations that were achieved by solving for the variation with height of Morton's non-dimensional flux parameter Γ - this single flux parameter gives a unique representation of the behaviour of steady plumes and enables a characterization of the different types of plume. We discuss advantages of solutions in this form before describing extensions to plume theory and suggesting directions for new research. © 2010 The Author. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.

Lyapunov theory for zeno stability

Zeno behavior is a dynamic phenomenon unique to hybrid systems in which an infinite number of discrete transitions occurs in a finite amount of time. This behavior commonly arises in mechanical systems undergoing impacts and optimal control problems, but its characterization for general hybrid systems is not completely understood. The goal of this paper is to develop a stability theory for Zeno hybrid systems that parallels classical Lyapunov theory; that is, we present Lyapunov-like sufficient conditions for Zeno behavior obtained by mapping solutions of complex hybrid systems to solutions of simpler Zeno hybrid systems defined on the first quadrant of the plane. These conditions are applied to Lagrangian hybrid systems, which model mechanical systems undergoing impacts, yielding simple sufficient conditions for Zeno behavior. Finally, the results are applied to robotic bipedal walking. © 2012 IEEE.

Transient Ventilation Dynamics Following a Change in Strength of a Point Source of Heat

We investigate the transient ventilation flow within a confined ventilated space, with high- and low-level openings, when the strength of a low-level point source of heat is changed instantaneously. The steady-flow regime in the space involves a turbulent buoyant plume, which rises from the point source to a well-mixed warm upper layer. The steady-state height of the interface between this layer and the lower layer of exterior fluid is independent of the heat flux, but the upper layer becomes progressively warmer with heat flux. New analogue laboratory experiments of the transient adjustment between steady states identify that if the heat flux is increased, the continuing plume propagates to the top of the room forming a new, warmer layer. This layer gradually deepens, and as the turbulent plume entrains fluid from the original warm layer, the original layer is gradually depleted and disappears, and a new steady state is established. In contrast, if the source buoyancy flux is decreased, the continuing plume is cooler than the original plume, so that on reaching the interface it is of intermediate density between the original warm layer and the external fluid. The plume supplies a new intermediate layer, which gradually deepens with the continuing flow. In turn, the original upper layer becomes depleted, both as a result of being vented through the upper opening of the space, but also due to some penetrative entrainment of this layer by the plume, as the plume overshoots the interface before falling back to supply the new intermediate layer. We develop quantitative models which are in good accord with our experimental data, by combining classical plume theory with models of the penetrative entrainment for the case of a decrease in heating. Typically, we find that the effect of penetrative entrainment on the density of the intruding layer is relatively weak, provided the change in source strength is sufficiently large. However, penetrative entrainment measurably increases the rate at which the depth of the draining layer decreases. We conclude with a discussion of the importance of these results for the control of naturally ventilated spaces.

Development of a classical force field for the oxidized Si surface: application to hydrophilic wafer bonding.

We have developed a classical two- and three-body interaction potential to simulate the hydroxylated, natively oxidized Si surface in contact with water solutions, based on the combination and extension of the Stillinger-Weber potential and of a potential originally developed to simulate SiO(2) polymorphs. The potential parameters are chosen to reproduce the structure, charge distribution, tensile surface stress, and interactions with single water molecules of a natively oxidized Si surface model previously obtained by means of accurate density functional theory simulations. We have applied the potential to the case of hydrophilic silicon wafer bonding at room temperature, revealing maximum room temperature work of adhesion values for natively oxidized and amorphous silica surfaces of 97 and 90 mJm(2), respectively, at a water adsorption coverage of approximately 1 ML. The difference arises from the stronger interaction of the natively oxidized surface with liquid water, resulting in a higher heat of immersion (203 vs 166 mJm(2)), and may be explained in terms of the more pronounced water structuring close to the surface in alternating layers of larger and smaller densities with respect to the liquid bulk. The computed force-displacement bonding curves may be a useful input for cohesive zone models where both the topographic details of the surfaces and the dependence of the attractive force on the initial surface separation and wetting can be taken into account.

Analysis of CFRP strap-strengthened reinforced concrete beams using the modified compression field theory (MCFT)

External, prestressed carbon fiber reinforced polymer (CFRP) straps can be used to enhance the shear strength of existing reinforced concrete beams. In order to effectively design a strengthening system, a rational predictive theory is required. The current work investigates the ability of the modified compression field theory (MCFT) to predict the behavior of rectangular strap strengthened beams where the discrete CFRP strap forces are approximated as a uniform vertical stress. An unstrengthened control beam and two strengthened beams were tested to verify the predictions. The experimental results suggest that the MCFT could predict the general response of a strengthened beam with a uniform strap spacing < 0.9d. However, whereas the strengthened beams failed in shear, the MCFT predicted flexural failures. It is proposed that a different compression softening model or the inclusion of a crack width limit is required to reflect the onset of shear failures in the strengthened beams.