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.

The influence of geometry on jet plume development

Our recent efforts of using large-eddy simulation (LES) type methods to study complex and realistic geometry single stream and co-flow nozzle jets and acoustics are summarized in this paper. For the LES, since the solver being used tends towards having dissipative qualities, the subgrid scale (SGS) model is omitted, giving a numerical type LES (NLES). To overcome near wall streak resolution problems a near wall RANS (Reynolds averaged Navier-Stokes) model is smoothly blended in the LES making a hybrid RANS-NLES approach. Several complex nozzle geometries including the serrated (chevron) nozzle, realistic co-axial nozzles with eccentricity, pylon and wing-flap are discussed. The hybrid RANS-NLES simulations show encouraging predictions for the chevron jets. The chevrons are known to increase the high frequency noise at high polar angles, but decrease the low frequency noise at lower angles. The deflection effect of the potential core has an important mechanism of noise reduction. As for co-axial nozzles, the eccentricity, the pylon and the deployed wing-flap are shown to influence the flow development, especially the former to the length of potential core and the latter two having a significant impact on peak turbulence levels and spreading rates. The studies suggest that complex and real geometry effects are influential and should be taken into count when moving towards real engine simulations. © 2012 Elsevier Ltd. All rights reserved.

Transient natural ventilation of a space with localised heating

The natural ventilation of a well-mixed, pre-heated room with a point source of heating, and openings at the base and roof is investigated. The transient draining associated with the room being warmer than the exterior combined with the convective ow produced by the point source of heat leads to a fascinating series of transient ow regimes as the system evolves to the two-layer steady-state regime described by Linden, Lane-Ser_ and Smeed [1]. As the room begins to ventilate, a turbulent plume rises from the point source of heat to the ceiling, and typically forms a deepening layer of hot air. However, with a weak heat source, then at some point the ascending plume will intrude beneath the layer of original uid. Otherwise, the ascending plume always reaches the top of the room as the system evolves to a steady state. We develop a simpli_ed model of the transient evolution and test this with some new laboratory experiments. We conclude with a discussion of the implications of our results for real buildings.

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.

Non-aqueous phase liquid behavior in the subsurface: Transportation, source zone characterization and remediation

Contaminant behaviour in soils and fractured rock is very complex, not least because of the heterogeneity of the subsurface environment. For non-aqueous phase liquids (NAPLs), a liquid density contrast and interfacial tension between the contaminant and interstitial fluid adds to the complexity of behaviour, increasing the difficulty of predicting NAPL behaviour in the subsurface. This paper outlines the need for physical model tests that can improve fundamental understanding of NAPL behaviour in the subsurface, enhance risk assessments of NAPL contaminated sites, reduce uncertainty associated with NAPL source remediation and improve current technologies for NAPL plume remediation. Four case histories are presented to illustrate physical modelling approaches that have addressed problems associated with NAPL transport, remediation and source zone characterization. © 2006 Taylor & Francis Group, London.

Analysis of basic processes inside the keyhole during deep penetration ND-YAG CW laser welding

During laser welding, the keyhole is generated by the recoil pressure induced by the evaporation processes occurring mainly on the front keyhole wall (KW). In order to characterize the evaporation process, we have measured this recoil pressure by using a plume deflection technique, where the plume generated for static conditions (i. e. with no sample displacement) is deflected by a transverse side gas jet. From the measurement of the plume deflection angle, the recoil pressure can be determined as a function of incident intensity and sample material. From these data one can estimate the pressure generated on the front KW, during laser welding. Therefore, the corresponding dynamic pressure exerted by the vapor plume expansion on the rear KW, in contact with the melt pool, can be also estimated. These pressures appear to be in close agreement with those generated by an additional side jet that has been used in previous experiments, for stabilizing the observed melt pool oscillations or fluctuations.

Laser-vapour interaction in high-power cw Nd:YAG laser welding

During high-power cw Nd:YAG laser welding a vapour plume is formed containing vaporised material ejected from the keyhole. Spectroscopic studies of the vapour emission have demonstrated that the vapour can be considered as thermally excited gas with a stable temperature (less than 3000K), not as partially ionised plasma. In this paper, a review of temperatures in the vapour plume is presented. The difficulties in the analysis of the plume spectroscopic results are reviewed and explained. It is shown that particles present in the vapour interact with the laser beam, attenuating it. The attenuation can be calculated with Mie scattering theory, however, vaporisation and particle formation also both play a major role in this process. The laser beam is also defocused due to the scattering part of the attenuation mechanism, changing the energy density in the laser beam. Methods for mitigating the effects of the laser beam-vapour interaction, using control gases, are presented together with their advantages and disadvantages. This 'plume control' has two complementary roles: firstly, the gas must divert the vapour plume from out of the laser beam path, preventing the attenuation. Secondly, the gas has to stabilise the front wall of the keyhole, to prevent porosity formation.

Displacement and mixing ventilation driven by opposing wind and buoyancy

The effect of an opposing wind on the stratification and flow produced by a buoyant plume rising from a heat source on the floor of a ventilated enclosure is investigated. Ventilation openings located at high level on the windward side of the enclosure and at low level on the leeward side allow a wind-driven flow from high to low level, opposite to the buoyancy-driven flow. One of two stable steady flow regimes is established depending on a dimensionless parameter F that characterizes the relative magnitudes of the wind-driven and buoyancy-driven velocities within the enclosure, and on the time history of the flow. A third, unstable steady flow solution is identified. For small opposing winds (small F) a steady, two-layer stratification and displacement ventilation is established. Exterior fluid enters through the lower leeward openings and buoyant interior fluid leaves through the upper windward openings. As the wind speed increases, the opposing wind may cause a reversal in the flow direction. In this case, cool exterior fluid enters through the high windward openings and mixes the interior fluid, which exits through the leeward openings. There are now two possibilities. If the rate of heat input by the source exceeds the rate of heat loss through the leeward openings, the temperature of the interior increases and this flow reversal is only maintained temporarily. The buoyancy force increases with time, the flow reverts to its original direction, and steady two-layer displacement ventilation is re-established and maintained. In this regime, the increase in wind speed increases the depth and temperature of the warm upper layer, and reduces the ventilation flow rate. If, on the other hand, the heat loss exceeds the heat input, the interior cools and the buoyancy-driven flow decreases. The reversed flow is maintained, the stratification is destroyed and mixing ventilation occurs. Further increases in wind speed increase the ventilation rate and decrease the interior temperature. The transitions between the two ventilation flow patterns exhibit hysteresis. The change from displacement ventilation to mixing ventilation occurs at a higher F than the transition from mixing to displacement. Further, we find that the transition from mixing to displacement ventilation occurs at a fixed value of F, whereas the transition from displacement to mixing flow is dependent on the details of the time history of the flow and the geometry of the openings, and is not determined solely by the value of F. Theoretical models that predic t the steady stratification profiles and flow rates for the displacement and mixing ventilation, and the transitions between them, are presented and compared with measurements from laboratory experiments. The transition between these ventilation patterns completely changes the internal environment, and we discuss some of the implications for the natural ventilation of buildings. © 2004 Cambridge University Press.

Thermal stratification produced by plumes and jets in enclosed spaces

The airflow and thermal stratification produced by a localised heat source located at floor level in a closed room is of considerable practical interest and is commonly referred to as a 'filling box'. In rooms with low aspect ratios H/R ≲ 1 (room height H to characteristic horizontal dimension R) the thermal plume spreads laterally on reaching the ceiling and a descending horizontal 'front' forms separating a stably stratified, warm upper region from cooler air below. The stratification is well predicted for H/R ≲ 1 by the original filling box model of Baines and Turner (J. Fluid. Mech. 37 (1968) 51). This model represents a somewhat idealised situation of a plume rising from a point source of buoyancy alone-in particular the momentum flux at the source is zero. In practical situations, real sources of heating and cooling in a ventilation system often include initial fluxes of both buoyancy and momentum, e.g. where a heating system vents warm air into a space. This paper describes laboratory experiments to determine the dependence of the 'front' formation and stratification on the source momentum and buoyancy fluxes of a single source, and on the location and relative strengths of two sources from which momentum and buoyancy fluxes were supplied separately. For a single source with a non-zero input of momentum, the rate of descent of the front is more rapid than for the case of zero source momentum flux and increases with increasing momentum input. Increasing the source momentum flux effectively increases the height of the enclosure, and leads to enhanced overturning motions and finally to complete mixing for highly momentum-driven flows. Stratified flows may be maintained by reducing the aspect ratio of the enclosure. At these low aspect ratios different long-time behaviour is observed depending on the nature of the heat input. A constant heat flux always produces a stratified interior at large times. On the other hand, a constant temperature supply ultimately produces a well-mixed space at the supply temperature. For separate sources of momentum and buoyancy, the developing stratification is shown to be strongly dependent on the separation of the sources and their relative strengths. Even at small separation distances the stratification initially exhibits horizontal inhomogeneity with localised regions of warm fluid (from the buoyancy source) and cool fluid. This inhomogeneity is less pronounced as the strength of one source is increased relative to the other. Regardless of the strengths of the sources, a constant buoyancy flux source dominates after sufficiently large times, although the strength of the momentum source determines whether the enclosure is initially well mixed (strong momentum source) or stably stratified (weak momentum source). © 2001 Elsevier Science Ltd. All rights reserved.

Global civil aviation black carbon emissions.

Aircraft black carbon (BC) emissions contribute to climate forcing, but few estimates of BC emitted by aircraft at cruise exist. For the majority of aircraft engines the only BC-related measurement available is smoke number (SN)-a filter based optical method designed to measure near-ground plume visibility, not mass. While the first order approximation (FOA3) technique has been developed to estimate BC mass emissions normalized by fuel burn [EI(BC)] from SN, it is shown that it underestimates EI(BC) by >90% in 35% of directly measured cases (R(2) = -0.10). As there are no plans to measure BC emissions from all existing certified engines-which will be in service for several decades-it is necessary to estimate EI(BC) for existing aircraft on the ground and at cruise. An alternative method, called FOX, that is independent of the SN is developed to estimate BC emissions. Estimates of EI(BC) at ground level are significantly improved (R(2) = 0.68), whereas estimates at cruise are within 30% of measurements. Implementing this approach for global civil aviation estimated aircraft BC emissions are revised upward by a factor of ~3. Direct radiative forcing (RF) due to aviation BC emissions is estimated to be ~9.5 mW/m(2), equivalent to ~1/3 of the current RF due to aviation CO2 emissions.

Updated correlation between aircraft smoke number and black carbon concentration

Aircraft emissions of black carbon (BC) contribute to anthropogenic climate forcing and degrade air quality. The smoke number (SN) is the current regulatory measure of aircraft particulate matter emissions and quantifies exhaust plume visibility. Several correlations between SN and the exhaust mass concentration of BC (CBC) have been developed, based on measurements relevant to older aircraft engines. These form the basis of the current standard method used to estimate aircraft BC emissions (First Order Approximation version 3 [FOA3]) for the purposes of environmental impact analyses. In this study, BC with a geometric mean diameter (GMD) of 20, 30, and 60 nm and filter diameters of 19 and 35 mm are used to investigate the effect of particle size and sampling variability on SN measurements. For BC with 20 and 30 nm GMD, corresponding to BC emitted by modern aircraft engines, a smaller SN results from a given CBC than is the case for BC with 60 nm GMD, which is more typical of older engines. An updated correlation between CBC and SNthat accounts for typical size of BC emitted by modern aircraft is proposed. An uncertainty of ±25% accounts for variation in GMD in the range 20-30 nm and for the range of filter diameters. The SN-CBC correlation currently used in FOA3 underestimates by a factor of 2.5-3 for SN <15, implying that current estimates of aircraft BC emissions derived from SN are underestimated by the same factor. Copyright © American Association for Aerosol Research.

Nanoscale-accuracy transfer printing of ultra-thin AlInGaN light-emitting diodes onto mechanically flexible substrates

The transfer printing of 2 μm-thick aluminum indium gallium nitride (AlInGaN) micron-size light-emitting diodes with 150 nm (±14 nm) minimum spacing is reported. The thin AlInGaN structures were assembled onto mechanically flexible polyethyleneterephthalate/polydimethylsiloxane substrates in a representative 16 × 16 array format using a modified dip-pen nano-patterning system. Devices in the array were positioned using a pre-calculated set of coordinates to demonstrate an automated transfer printing process. Individual printed array elements showed blue emission centered at 486 nm with a forward-directed optical output power up to 80 μW (355 mW/cm 2) when operated at a current density of 20 A/cm2. © 2013 AIP Publishing LLC.