考虑动量和浮力通量耗散的烟气抬升解析模式

环流湍流对烟流抬升有重要的作用。本文考虑了环境湍流引起的烟流的动量和浮力通量耗散，导出了烟流抬升的控制方程，给出了烟流的轨迹方程和中性条件下热浮力烟流的终极抬升高度。在烟囱附近，轨迹方程接近2/3次律；当下风距离增大时，与2/3次律有明显偏离并逐渐变平。与外场试验资料的比较表明，本模式能较好地模拟浮力烟流的轨迹和终极高度。     

On the dispersion parameters of plumes from tall stacks in a shoreline environment

Experimental data for buoyant plumes released from high sources into layers having little ambient turbulence show that plume dispersion parameters vary in a manner similar to that during initial plume rise. This is consistent with general plume rise theory. Dispersion of plumes from tall stacks in a shoreline environment where a thermal internal boundary layer is formed often demonstrates this behaviour.     

Unification Of Non-Dimensional Solutions To Asymptotic Equations For Plumes Of Different Shape

A model of buoyant fluid rise is developed in aLagrangian framework. Results from the model arevalidated by comparison with laboratory and fieldexperiments. The model is sufficiently general toinclude geometries of any cross-sectional shape. Exact solutions to the asymptotic equations forcontinuous discharges from line and point sources andfor an instantaneous discharge from a point source areconsidered. Prismatic, cylindrical andspherical shapes approximate these three geometries,respectively. Accommodation of these shapes withinthe same general model allows for direct comparison. It is shown that, for discharges into a cross-flowthat may be stratified or unstratified, thenon-dimensional trajectory, thickness and dilution canbe uniquely specified using three parameters. Theseare the non-dimensional size of the source, therelative importance of the initial fluxes of momentumand buoyancy and the number of orthogonal axes throughwhich entrainment can occur. Such non-dimensionalresults are particularly useful for examining thosescenarios for which there are limited experimentaldata.     

Penetrative convection in rapidly rotating flows: preliminary results from numerical simulation

Turbulent convection forced by a surface heat flux into a stably stratified region is a feature of both the atmospheric and oceanic planetary boundary layers. Of particular interest is the interface between the convective layer and the stable stratification, where the entrainment of fluid into the convective layer by penetrating plumes may lead to a reverse buoyancy flux, and an enhancement of the stable stratification. Whereas in the atmosphere the influence of rotation on this penetrative convection is negligible, oceanic convection may be subjected to lower Rossby numbers and hence greater rotational influence. To isolate the effects of rotation, we present three numerical solutions for turbulent penetrative convection, characterised by different rotation rates, with all other parameters being held constant. Our results indicate that at lower Rossby numbers the lateral scale of the plumes is reduced, whereas the vertical vorticity of the plumes is much enhanced. Vertical transports of buoyancy and kinetic energy across the convective layer are reduced, leading to less efficient penetration at the interface with the stratified layer, and hence less reverse buoyancy flux in this region.     

Interaction of a negatively buoyant line plume with a density interface

Laboratory experiments were carried out to investigate the interaction between turbulent line buoyant plumes and sharp density interfaces, with the aim of using the results to interpret oceanic observations pertinent to crack openings in the polar ice-cap (leads). These openings take the form of long narrow channels, and are often modeled as line bouyant plumes. The plumes descend as in a homogenoous fluid, impinge on the density interface, and then spread horizontally as gravity currents. Depending on the Richardson number , where Δb is the buoyancy jump across the interface, l_D is the half-width of the plume before the impingement and q₀is the buoyancy flux per unit length of the source, different flow patterns were identified. When Ri < 0.5, the plumes penetrate deep into the bottom layer, deflect horizontally and then spread while showing little vertical rise. When 0.6 < Ri < 5, the penetration is significant, but the fluid bounces back after entraining heavy fluid from the lower layer and then spreads horizontally above the interface as a gravity current. Appreciable mixing between this current and the lower layer was detected when Ri <1. When Ri > 10, the penetration was small and a sharp-nosed gravity current emerged some time after the impact. Measurements were made on the penetration depth, the velocities of the gravity current and the subsurface flow towards the plume, the entrainment rate and other wave parameters. Possible implications of the results for oceanic cases are also discussed.     

The Entrainment Rate for Buoyant Plumes in a Crossflow

We consider large-eddy simulations (LES) of buoyant plumes from a circular source with initial buoyancy flux F ₀ released into a stratified environment with constant buoyancy frequency N and a uniform crossflow with velocity U. We make a systematic comparison of the LES results with the mathematical theory of plumes in a crossflow. We pay particular attention to the limits [(U)\tilde] << 1{\tilde{U}\ll1} and [(U)\tilde] >> 1{\tilde{U}\gg 1}, where [(U)\tilde]=U/(F₀ N)^1/4{\tilde{U}=U/(F_0 N)^{1/4}}, for which analytical results are possible. For [(U)\tilde] >> 1{\tilde{U}\gg 1}, the LES results show good agreement with the well-known two-thirds law for the rise in height of the plume. Sufficiently far above the source, the centreline vertical velocity of the LES plumes is consistent with the analytical z ^−1/3 and z ^−1/2 scalings for respectively [(U)\tilde] << 1{\tilde{U}\ll 1} and [(U)\tilde] >> 1{\tilde{U}\gg 1}. In the general case, where the entrainment is assumed to be the sum of the contributions from the horizontal and vertical velocity components, we find that the discrepancy between the LES data and numerical solutions of the plume equations is largest for [(U)\tilde]=O(1){\tilde{U}=O(1)}. We propose a modified additive entrainment assumption in which the contributions from the horizontal and vertical velocity components are not equally weighted. We test this against observations of the plume generated by the Buncefield fire in the U.K. in December 2005 and find that the results compare favourably. We also show that the oscillations of the plume as it settles down to its final rise height may be attenuated by the radiation of gravity waves. For [(U)\tilde] << 1{\tilde{U}\ll 1} the oscillations decay rapidly due to the transport of energy away from the plume by gravity waves. For ${\tilde{U}>rsim 1}${\tilde{U}>rsim 1} the gravity waves travel in the same direction and at the same speed as the flow. In this case, the oscillations of the plume do not decay greatly by radiation of gravity waves.     

The Fluxes and Behaviour of Plumes Inferred from Measurements of Coherent Structures within Images of the Bulk Flow

This paper describes how measurements of the movement of identifiable features at the edge of a turbulent plume can be interpreted to determine the properties of the mean flow and consequently, using plume theory, can be used to make estimates of the fluxes of volume (mass), momentum, and buoyancy in a plume. This means that video recordings of smoke rising from a chimney or buoyant material from a source on the sea bed can be used to make accurate estimates of the source conditions for the plume. At best we can estimate the volume flux and buoyancy flux to within about 5% and 15% of the actual values, respectively. Although this is restricted to the case of a plume rising in a stationary and unstratified environment, we show that the results may be of practical use in other more complex situations. In addition, we demonstrate that large-scale (turbulent) coherent structures at the plume edge form on a scale approximately 40% of the local (mean) plume half-width and travel at almost 60% of the average local (mean) velocity in the plume.     

ON THE EFFECT OF AMBIENT TURBULENCE ON THE BUOYANT PLUME RISE

Considering the effects of both entrainment and intensification of the exchanges of momentum and heat between plume and ambient air,we have derived the trajectory equation of buoyant plume under neutral conditions,and the final rise prediction formula theoretically without any hypotheses on the rise termination.Compared with the field experiments,the final rise formula simulates the observed final rise heights well.     

A one-dimensional model for the parameterization of deep convection in the ocean

A one-dimensional penetrative plume model has been constructed to parameterize the process of deep convection in ocean general circulation models (OGCMs). This research is motivated by the need for OGCMs to better model the production of deep and intermediate water masses. The parameterization scheme takes the temperature and salinity profiles of OGCM grid boxes and simulates the subgrid-scale effects of convection using a one-dimensional parcel model. The model moves water parcels from the surface layer down to their level of neutral buoyancy, simulating the effect of convective plumes. While in transit, the plumes exchange water with the surrounding environment; however, the bulk of the plume water mass is deposited at e level of neutral buoyancy. Weak upwelling around the plumes is included to maintain an overall mass balance. The process continues until the negative buoyant energy of the one-dimensional vertical column is minimized. The parameterized plume entrainment rate, which plays a central role in the parameterization, is calculated using modified equations based on the physics of entraining buoyant plumes. This scheme differs from the convective adjustment techniques currently used in OGCMs, because the parcels penetrate downward with the appropriate degree of mixing until they reach their level of neutral stability.     

考虑环境湍流作用的烟气上升路径方程

本文利用我国的观测资料对中性层结时的烟气拾升路径作了细致的分析。发现抬升高度与距离之间符合幂次关系,但与不考虑环境湍流作用的2/3次律有系统偏离;这一偏离与烟流半径和高度的非线性关系相一致。上述事实与国外的一些观测结果均表明,环境湍流对烟气抬升的累积作用不容忽视。为此,提出了环境湍流与自生湍流的联合作用模式,由此导得的抬升路径方程可以满意地解释上述观测事实。模式还显示,卷挟速度除取决于切变速度及湍流强度以外,还是烟流半径的函数,因而在物理上更加合理。由联合作用模式导出的终极抬升公式比不考虑环境湍流累积作用的公式更符合观测结果。     

Rise of Buoyant Emissions from Low-Level Sources in the Presence of Upstream and Downstream Obstacles

Field and laboratory studies have been conducted to investigate the effect of surrounding buildings on the plume rise from low-level buoyant sources, such as distributed power generators. The field experiments were conducted in Palm Springs, California, USA in November 2010 and plume rise from a 9.3 m stack was measured. In addition to the field study, a laboratory study was conducted in a water channel to investigate the effects of surrounding buildings on plume rise under relatively high wind-speed conditions. Different building geometries and source conditions were tested. The experiments revealed that plume rise from low-level buoyant sources is highly affected by the complex flows induced by buildings stationed upstream and downstream of the source. The laboratory results were compared with predictions from a newly developed numerical plume-rise model. Using the flow measurements associated with each building configuration, the numerical model accurately predicted plume rise from low-level buoyant sources that are influenced by buildings. This numerical plume rise model can be used as a part of a computational fluid dynamics model.     

论环境湍流对烟气抬升的作用

考虑环境湍流对烟气抬升具有增强烟气夹卷和促进烟气与环境空气间动、热量交换的双重作用，导出了中性层结条件下热烟气浮升的轨迹方程，并且无需对抬升的截止作任何假设，导得了理论的终极抬升公式。经与国内外烟气抬升实测资料对比，发现这些结果与实况基本一致。     

Dispersion of spores released from an elevated line source within a wheat canopy

Turbulent dispersion of spores was studied near a source located inside a wheat canopy. Two colors of Lycopodium spores were released simultaneously at a steady rate from line sources at two heights (0.4–0.5 m and 0.7–0.8 m) in a 0.8 to 1.0 m tall crop. The number of spores of each color released was estimated by weighing the sources before and after each release. Aerial spore concentrations were measured at 2 and 4 m downwind of the sources using rotorods placed at four heights above the canopy and small suction traps at two heights inside the canopy. Concentrations near the ground were estimated from deposits on sticky glass microscope slides placed on the ground. Experiments were conducted on six different days. Friction velocities ranged from about 0.3 to 0.5 m s^–1. The average horizontal fluxes of spores were calculated as the product of the observed concentrations and average wind speeds. At a distance of 2 m downwind from the sources, more than 16 to 44% of the flux of spores released from the lower source and more than 41 to 50% of the flux of spores released from the upper source were estimated to be above the canopy. These fluxes were compared with fluxes calculated using both a K-theory model and a random-flight-fluid-parcel-trajectory simulation model. The fluxes predicted by the models were generally considerably less than the values determined experimentally.     

Reynolds-Number Dependence of Gas Dispersion Over a Wavy Wall

Direct numerical simulations of an Ekman layer are performed to study flow evolution during the response of an initially neutral boundary layer to stable stratification. The Obukhov length, L, is varied among cases by imposing a range of stable buoyancy fluxes at the surface to mimic ground cooling. The imposition of constant surface buoyancy flux , i.e. constant-flux stability, leads to a buoyancy difference between the ground and background that tends to increase with time, unlike the constant-temperature stability case where a constant surface temperature is imposed. The initial collapse of turbulence in the surface layer owing to surface cooling that occurs over a time scale proportional to \(L/u_*\), where \(u_*\) is the friction velocity, is followed by turbulence recovery. The flow accelerates, and a “low-level jet” (LLJ) with inertial oscillations forms during the turbulence collapse. Turbulence statistics and budgets are examined to understand the recovery of turbulence. Vertical turbulence exchange, primarily by pressure transport, is found to initiate fluctuations in the surface layer and there is rebirth of turbulence through enhanced turbulence production as the LLJ shear increases. The turbulence recovery is not monotonic and exhibits temporal intermittency with several collapse/rebirth episodes. The boundary layer adjusts to an increase in the surface buoyancy flux by increased super-geostrophic velocity and surface stress such that the Obukhov length becomes similar among the cases and sufficiently large to allow fluctuations with sustained momentum and heat fluxes. The eventual state of fluctuations, achieved after about two inertial periods (\(ft \approx 4\pi \)), corresponds to global intermittency with turbulent patches in an otherwise quiescent background. Our simplified configuration is sufficient to identify turbulence collapse and rebirth, global and temporal intermittency, as well as formation of low-level jets, as in observations of the stratified atmospheric boundary layer.     

Nonuniform Distribution Of High-Frequency Turbulence In The Unstable Boundary Layer

Turbulent flow data of wind velocity and temperature in the unstably stratifiedatmospheric boundary layer, derived from steel tower observations in the field and wind-tunnel experiments were used to study the relationship between the plumes and the small-scale eddies in the inertial subrange. Flow visualisation experiments in the wind tunnel were also conducted to observe the structure of the flow in the plumes, and time series data were analysed by using wavelet transforms. The results show that variances of velocity and temperature due to the small-scale eddies are large in the plumes and small outside of the plumes, and that the momentum and heat fluxes due to the small-scale eddies follow the same tendency as found in the variances. The ratios of the variances caused by the small-scale eddies in the plumes to the whole of the variances caused by the small-scale eddies in and out the plumes increase with non-dimensional height -z/L in which L is the local Obukhov length. Similar ratios of the fluxes caused by the small-scale eddies also show the same tendency. These ratios can be expressed as functions of -z/L for results based on field observationand the wind tunnel experiments. This relation hardly changes even if the wavelet function is changed. The flow visualisation experiments show that the plumes have a complicated structure in which mushroom type flows are stacked on top of each other. This characteristic structure seems to increase the energy of the small-scale eddies in the plumes.     

Spectral Analysis of Turbulent Aerosol Fluxes by Fourier Transform,Wavelet Analysis,and Multiresolution Decomposition

Aerosol number fluxes are spectrally analyzed using fast Fourier transform analysis, wavelet analysis and multiresolution decomposition. All three methods yield similar spectral features in general, although a detailed evaluation of the cospectra shows some differences, e.g. due to different resolutions in the time and frequency domains. Wavelet analysis yields aerosol flux estimates with a high time resolution that can be used to assess the flux variability. Multiresolution decomposition has been applied successfully to evaluate cospectra of the aerosol number flux, the buoyancy flux and the momentum flux of three 1-day datasets from diverse environments. For all scalars and all environments, the dimensionless frequency (f) of the cospectral peak was found between $f = 0.1$ and 0.2. In addition, the cospectral gap time scale of the aerosol number flux was found between 100 and 1,000 s. Thus, in this study several spectral features such as the dominant time scale and the cospectral gap time scale of aerosol number fluxes are similar to buoyancy fluxes. However, the shape of aerosol number flux cospectra often deviates from buoyancy and momentum flux cospectra, especially at very small and at very large time scales.     

A statistical representation of sub-grid variation in mixing-length models of the boundary layer

In conventional gradient-transfer models of the atmospheric boundary layer, the fluxes are calculated using the gradients of buoyancy and momentum over the grid interval. When using finite-differences, these gradients are approximated using what are essentially the grid-volume means of the quantities involved. This is unrealistic, as within a large volume of the atmosphere fluxes are supported by an ensemble of small-scale motions, which may be organized by the presence of large coherent eddies. In this paper, a new method of computing the flux is developed, in which the gradients of grid-volume mean buoyancy and momentum are replaced by grid-volume ensembles. The distributional parameters of the ensembles are derived using dimensional arguments, augmented by the results of a large-eddy simulation.     

The vertical profile of concentration fluctuations in near-surface plumes

Field experiments on concentration fluctuations have frequently measured horizontal cross-sections of fluctuation statistics through plumes at fixed heights near the surface, but have not considered the effect of height above the ground in any detail. A set of tracer experiments designed to measure vertical profiles of concentration fluctuations in plumes, with a range of source heights, is described, and profiles of statistics are presented. Considerable variation of the statistics with both source and detector height is found. Near the surface, fluctuation intensity is a minimum and the time and length scales of the fluctuations are greatly increased. Profiles are consistent with the idea that concentration fluctuations near the surface are like those higher up at a greater distance from the source. Lowering the source height reduces the fluctuation intensity at all heights, and also alters the form of the concentration PDF. Results may be explained by the reduced length scale of sheargenerated turbulence near the surface causing enhanced small-scale mixing, which rapidly smooths out much of the fine structure with the plume.     

Methods to Estimate Surface Fluxes of Momentum and Heat from Routine Weather Observations for Dispersion Applications under Stable Stratification

We examine the efficacy of two methods commonly used to estimate the vertical turbulent fluxes of momentum and sensible heat from routinely observed mean quantities in the surface layer under stable stratification. The single-level method uses mean wind speed and temperature measurements at a single height, whereas the two-level method uses mean wind speed measurements at a single height and mean temperature measurements at two heights. These methods are used in popular meteorological processors such as the U.S. Environmental Protection Agency approved AERMET and CALMET to generate inputs for dispersion simulations. We use data from a flux station of the U.K. Met Office at Cardington for comparison. It is found that the single-level method does not describe the flux variation in the weakly stable regime at all, because of its assumption that the temperature scale, i.e. the ratio of the kinematic sensible heat flux to the friction velocity, is constant, which is plausible only under strongly stable conditions. On the other hand, the two-level method provides a physically realistic variation of the fluxes with stability, but the required temperature measurements at two levels are usually not available on a routine basis. If measurements of the standard deviation of temperature are also available, in addition to the mean temperature at a single level, then they can be usefully employed in a third (single-level) method, with the consequence that the computed fluxes are very similar to those obtained from the two-level method. An improvement to the original single-level method is considered, and flux calculations under low wind conditions are also discussed.     

Flux Corrections Revisited

In a horizontally homogeneous atmosphere, the calculation of thevertical fluxes of mass, momentum, and heat from point measurementsare to lowest order calculated as the standard eddy covariance.Additional terms primarily related to the mean vertical velocity areadded to improve the accuracy of the calculation. A review ofprevious work is undertaken. The fluxes are reformulated beginningwith the equations of conservation of mass, momentum, and energy. Aformula for the mean vertical velocity is found that is accurate tothe ideal gas approximation. Relations are developed between fluxesat different heights. In addition to the relations for fluxes ofmass, sensible heat, and latent heat, the fluxes of momentum andturbulent kinetic energy are also considered.