首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
In many atmospheric flows, a dispersed phase is actively suspended by turbulence, whose competition with gravitational settling ultimately dictates its vertical distribution. Examples of dispersed phases include snow, sea-spray droplets, dust, or sand, where individual elements of much larger density than the surrounding air are carried by turbulent motions after emission from the surface. In cases where the particle is assumed to deviate from local fluid motions only by its gravitational settling (i.e., they are inertialess), traditional flux balances predict a power-law dependence of particle concentration with height. It is unclear, however, how particle inertia influences this relationship, and this question is the focus of this work. Direct numerical simulations are conducted of turbulent open-channel flow, laden with Lagrangian particles of specified inertia; in this way the study focuses on the turbulent transport which occurs in the lowest few meters of the planetary boundary layer, in regions critical for connecting emission fluxes to the fluxes felt by the full-scale boundary layer. Simulations over a wide range of particle Stokes number, while holding the dimensionless settling velocity constant, are performed to understand the role of particle inertia on vertical dispersion. It is found that particles deviate from their inertialess behaviour in ways that are not easily captured by traditional theory; concentrations are reduced with increasing Stokes number. Furthermore, a similarity-based eddy diffusivity for particle concentration fails as particles experience inertial acceleration, precluding a closed-form solution for particle concentration as in the case of inertialess particles. The primary consequence of this result is that typical flux parametrizations connecting surface emission models (e.g., saltation models or sea-spray generation functions) to elevated boundary conditions may overestimate particle concentrations due to the reduced vertical transport caused by inertia in between; likewise particle emission may be underestimated if inferred from concentration measurements aloft.  相似文献   

2.
An atmospheric surface-layer model is used to investigate the interactionbetween suspended snow particles and the near-surface flow. Themodel incorporates the effects of upward diffusion, gravitational settling and sublimation of snow particles in 48 size classes, the effects of snowdrift sublimation on the heat and moisture budget of the surface layer, and the buoyancy destruction of turbulent kinetic energy (TKE) caused by the presence of suspended particles. A new term in the E- closure model representing the buoyancy destruction due to suspended particles is included in the prognostic equation for TKE. Generally, model results indicate that the presence of suspended particles causes significant decreases in TKE, the dissipation rate, turbulent length scales and eddy exchange coefficients (up to 40%). It is found that the reduction in the eddy exchangecoefficients is due mainly to reductions in turbulent length scales. Theassociated particle Richardson number peaks near the saltation-suspensioninterface, but at higher levels in the surface layer the particle-induced buoyancy can also significantly affect the flow. A detailed analysis of the various snowdrift quantities, the TKE budget and the particle buoyancy effects on the flow is presented.  相似文献   

3.
Snowdrift is one of the manymanifestations of two-phase flow, in which theinteraction between suspended particles and theambient fluid brings about some interesting features.Specifically, the drag required to keep particles insuspension against the downward gravitational pullrequires expenditure of turbulent kinetic energy(TKE). Other effects include the increased density of theair-snow mixture and the stable thermal stratificationcaused by the snowdrift sublimation-induced cooling.An atmospheric surface-layer model that includes snowdriftsuspension is described that includes the effects ofupward diffusion, gravitational settling andsublimation of snow particles in 48 size classes, theeffects of snowdrift sublimation on the heat andmoisture budget of the surface layer and the dampingof turbulence in the presence of suspended particles. Thewell-known E- closure model is applied toevaluate the eddy exchange coefficient, with a newterm representing buoyancy reduction induced by thestably stratified suspended particle profile includedin the prognostic equation for TKE.  相似文献   

4.
An expression for the vertical equilibrium concentration profile of heavy particles, including the effects of canopy on the eddy diffusivity as well as corrections for atmospheric stability, is proposed. This expression is validated against measurements of vertical concentration profiles of corn pollen above a corn field. The fitted theoretical profiles show that particle settling is correctly accounted for. The sensitivity to variations in the turbulent Schmidt number, settling velocity and stability corrections are explicitly characterized. The importance of independent measurements of the surface flux of pollen in future experiments is noted.  相似文献   

5.
A laboratory experiment was performed to investigate mixing across a density interface which separates two turbulent fluid layers and coexists with a stabilizing buoyancy flux. It was found that the buoyancy flux (q0) across the interface and through the turbulent layers (of depth D) becomes steady and constant in magnitude in the vertical direction, only when , where u is the horizontal r.m.s. velocity at the base of the mixed layers. The results suggest that mixing across the density interface is controlled by a dynamically important buoyancy gradient induced in the turbulent layers and that parameters such as the bulk Richardson number, , where Δb is the interfacial buoyancy jump, are of secondary importance. Measurements are used to infer the mixing mechanism at the interface, the mixing efficiency of stratified fluids and the entrainment law. Some geophysical applications of the results are also discussed.  相似文献   

6.
Backward Lagrangian stochastic models calculate particle trajectories in the atmosphere upstream of an observation point. A feature of these models is the ease with which the vertical flux contribution from surface area sources can be calculated. The flux contribution at an observation point P is found by summing the ratio of the particle vertical velocity at P to the touchdown velocity for particles which impact the ground within the source boundary.  相似文献   

7.
Dust particles frequently become mixtures of mineral dust and sea salt during their transport in the marine boundary layer, consequently growing in size, which causes changes in their settling velocities. In this study, the effect of sea salt on the gravitational settling of dust particles is investigated. Results show that the adhering of sea salt to dust particles can dramatically increase the gravitational settling of the particles, in particular if the particles become larger than 3–4 μm. Estimates with the observational data from six dust events in southwestern Japan revealed that, due to sea salt adhering, the gravitational settling flux of mineral dust increased approximately 14–17% in well-mixed events and 4–6% in less-mixed events, indicating a potential significant effect of sea salt on dust settling and the importance of considering this effect in the schemata of particle gravitational settling when mapping dust flux to the ocean.  相似文献   

8.
Flux densities of carbon dioxide were measured over an arid, vegetation-free surface by eddy covariance techniques and by a heat budget-profile method, in which CO2 concentration gradients were specified in terms of mixing ratios. This method showed negligible fluxes of CO2, consistent with the bareness of the experimental site, whereas the eddy covariance measurements indicated large downward fluxes of CO2. These apparently conflicting observations are in quantitative agreement with the results of a recent theory which predicts that whenever there are vertical fluxes of sensible or latent heat, a mean vertical velocity is developed. This velocity causes a mean vertical convective mass flux (= cw for CO2, in standard notation). The eddy covariance technique neglects this mean convective flux and measures only the turbulent flux c w. Thus, when the net flux of CO2 is zero, the eddy covariance method indicates an apparent flux which is equal and opposite to the mean convective flux, i.e., c w = – c w. Corrections for the mean convective flux are particularly significant for CO2 because cw and c w are often of similar magnitude. The correct measurement of the net CO2 flux by eddy covariance techniques requires that the fluxes of sensible and latent heat be measured as well.  相似文献   

9.
Reliable predictions of the daytime dispersal of heavy particles in the unstable atmospheric boundary layer are important in a variety of disciplines. For many applications, particles disperse from area sources near the ground, and corresponding theoretical solutions are desired to reveal insight into the physical processes. Here, theoretical solutions recently developed for neutral conditions are modified to include the effects of atmospheric instability. The Obukhov length L O and convection velocity w ? are introduced to characterize the patterns of particle dispersion, in additional to friction velocity u ? and settling velocity w s used in the neutral case. The major effects of atmospheric instability are accounted for by modifying the vertical velocity variance profile and considering the ratio of velocity scales w ?/u ?. Theoretical predictions including the mean concentration profile, plume height, and horizontal transport above the source, and ground deposition flux downwind from the source agree well with large-eddy simulation results while the particle plume is within the atmospheric surface layer. The deposition curve is characterized by a power-law decay whose exponent depends on u ?, w s, and w ?. A second steeper power-law develops once the plume extends into the mixed layer. This effect is enhanced with increasing atmospheric instability, implying that particles disperse farther from the source.  相似文献   

10.
The main goals of this work are climatological analysis of characteristics of vertical wind in the stratosphere and estimation of potential opportunities of its influence on stratospheric aerosol particles. High-altitude, temporal, and latitude dependences of zonal mean vertical wind velocity for the period of 1992?C2006 from the UKMO atmospheric general circulation model are analyzed. It is shown that monthly averaged amplitudes of the vertical wind are approximately ±5?mm/s, while annual averaged ones are ±1?mm/s. The upward wind can provide the vertical lifting against gravity for sufficiently large (up to 3?C5???m) aerosol particles with a density up to 1.0?C1.5?g/cm3 at stratospheric and mesospheric altitudes. The vertical wind, probably, is a substantial factor for particle motion up to altitudes of 30?C40?km and can change essentially the sedimentation velocities and the residence times of stratospheric aerosols. The structure of the averaged fields of vertical wind supposes the opportunity of formation of dynamically stable aerosol layers in the middle stratosphere. With the problem regarding the action of a permanent source of monodisperse particles near the stratopause taken as an example, it is shown that if the action of the averaged vertical component is taken into account along with the gravitational sedimentation and turbulent diffusion, the standard vertical profiles of the relative concentration of particles change cardinally. Estimations for the levitation heights for particles of different densities and sizes in the stratosphere under action of gravity and vertical wind pressure are presented.  相似文献   

11.
A three-dimensional (3-D) inertial particle – Lagrangian stochastic model for heavy particles in turbulent flows has been constructed. In this model, particle velocities are computed by adopting a non-linear drag law, while fluid velocity in the vicinity of a particle is calculated using a 3-D Langevin equation. Our model results have shown that the inclusion of the horizontal fluid velocity fluctuation computations and a non-linear drag law have an impact on the statistics of both fluid and particles when compared with our earlier one-dimensional (1-D) model with a linear drag law. Model results are compared and contrasted with Businger’s 1965 theory in terms of effective settling velocity.  相似文献   

12.
The derivation of the Panofsky–Dutton internal boundary-layer(IBL) height formula has been revisited. We propose that the upwindroughness length (rather than downwind) should be used in theformula and that a turbulent vertical velocity (w) ratherthan the surface friction velocity (u*) should be considered asthe appropriate scaling for the rate of propagation ofdisturbances into the turbulent flow. A published set ofwind-tunnel and atmospheric data for neutral stratification hasbeen used to investigate the influence of the magnitude ofroughness change on the IBL height.  相似文献   

13.
Two almost identical eddy covariance measurement setups were used to measure the fogwater fluxes to a forest ecosystem in the “Fichtelgebirge” mountains (Waldstein research site, 786 m a.s.l.) in Germany. During the first experiment, an intercomparison was carried out with both setups running simultaneously at the same measuring height on a meteorological tower, 12.5 m above the forest canopy. The results confirmed a close agreement of the turbulent fluxes between the two setups, and allowed to intercalibrate liquid water content (LWC) and gravitational fluxes. During the second experiment, the setups were mounted at a height of 12.5 and 3 m above the canopy, respectively. For the 22 fog events, a persistent negative flux divergence was observed with a greater downward flux at the upper level. To extrapolate the turbulent liquid water fluxes measured at height z to the canopy of height hc, a conversion factor 1/[1+0.116(zhc)] was determined. For the fluxes of nonvolatile ions, no such correction is necessary since the net evaporation of the fog droplets appears to be the primary cause of the vertical flux divergence. Although the net evaporation reduces the liquid water flux reaching the canopy, it is not expected to change the absolute amount of ions dissolved in fogwater.  相似文献   

14.
Turbulent fluctuations in CO2 concentrations over a paddy field are measured by a fastresponse device with an open sensing path. This IR device coupled with a sonic anemometer constitutes an eddy correlation instrument to measure CO2 fluxes. Three experiments were conducted in the surface layer over paddy 90 cm high. The stability (z – d)/L ranged from -0.14 to 0.20, where L denotes the Monin-Obukhov length.CO2 power spectra show the range of applicability of the -2/3 power law to be between f = 0.2 and f = 2, where f is the frequency normalized by wind speed and height. The cospectral estimate between CO2 and vertical component of wind speed ranging from f = 0.005 to f = 2 shows a peak at about f = 0.15 under near-neutral stratification.Hourly means of CO2 flux measured by the eddy correlation method increase with intensity of net radiation. The maximum value of downward flux of CO2 rises to 0.6 mg cm-2 hr-1 over the paddy field at the stage of ear emergence.Some turbulence statistics relating to the CO2 transport are evaluated: the correlation coefficient between CO2 and vertical velocity is about -0.3, and that between CO2 and humidity attains -0.7 ~ -0.8 under unstable stratification; nondimensional gradients c for CO2 and m for wind speed are 0.89 and 0.99, respectively.  相似文献   

15.
Wind-tunnel experiments were used to investigate the ground deposition of fine airborne particles in an array of idealized tree crowns. The particle ground deposition was modelled with a gaseous tracer instead of solid particles, which is an approach for very fine particles. A chemical method based on the reaction of ammonia and manganese chloride was used to quantify the mass transfer from the simulated atmospheric boundary-layer flow to the surface. Using a tracer gas instead of solid particles can be considered only if turbulent diffusion is the decisive deposition mechanism and effects of sedimentation, impaction, interception or molecular diffusion can be approximately ignored. These constraints are necessary due to scaling problems concerning particle modelling in the small-scale experiment. The intention was to determine the obstacle arrangement density in which the mean ground deposition is maximized for a defined crown form. A deposition amplification factor was defined as the quotient of deposition efficiencies for an area with tree crowns and an open ground with identical similarity parameters. Based on this calculation an increase of the ground deposition by up to 60% should be realistic through a favourable arrangement of tree crowns and tree number density. An increase in turbulence intensity in the flow leads to a significant amplification of the mean ground deposition.  相似文献   

16.
Statistics of atmospheric turbulence within and above a corn canopy   总被引:1,自引:2,他引:1  
Two three-dimensional split-film anemometers were used to measure turbulence statistics within and above a corn canopy. Normalised profiles of mean windspeed, root-mean-square velocity, momentum flux, and heat flux were constructed from half-hourly averages by dividing within-canopy measurements by the simultaneous canopy-top measurement. With the exception of the heat flux, these profiles showed consistent shape from day to day. Time series of the three velocity components were recorded on magnetic tape and subsequently analysed to obtain Eulerian time and length scales and the power spectrum of each component at several heights. The timescale was found to have a local minimum value at the top of the canopy. However the length scale L wformed from the timescale and the root-mean-square vertical velocity varied with height as L w 0.1 z. The power-spectra were non-dimensionalised to facilitate comparison of spectra at different heights and times. All spectra had -5/3 regions spanning at least two decades in frequency.  相似文献   

17.
This is the first of a series of three papers describing experiments on the dispersion of trace heat from elevated line and plane sources within a model plant canopy in a wind tunnel. Here we consider the wind field and turbulence structure. The model canopy consisted of bluff elements 60 mm high and 10 mm wide in a diamond array with frontal area index 0.23; streamwise and vertical velocity components were measured with a special three-hot-wire anemometer designed for optimum performance in flows of high turbulence intensity. We found that:
  1. The momentum flux due to spatial correlations between time-averaged streamwise and vertical velocity components (the dispersive flux) was negligible, at heights near and above the top of the canopy.
  2. In the turbulent energy budget, turbulent transport was a major loss (of about one-third of local production) near the top of the canopy, and was the principal gain mechanism lower down. Wake production was greater than shear production throughout the canopy. Pressure transport just above the canopy, inferred by difference, appeared to be a gain in approximate balance with the turbulent transport loss.
  3. In the shear stress budget, wake production was negligible. The role of turbulent transport was equivalent to that in the turbulent energy budget, though smaller.
  4. Velocity spectra above and within the canopy showed the dominance of large eddies occupying much of the boundary layer and moving downstream with a height-independent convection velocity. Within the canopy, much of the vertical but relatively little of the streamwise variance occurred at frequencies characteristic of wake turbulence.
  5. Quadrant analysis of the shear stress showed only a slight excess of sweeps over ejections near the top of the canopy, in contrast with previous studies. This is a result of improved measurement techniques; it suggests some reappraisal of inferences previously drawn from quadrant analysis.
  相似文献   

18.
Observations of vegetation–atmosphere exchange of carbon dioxide (CO2) by the eddy covariance (EC) technique are limited by difficult conditions such as nighttime and heterogeneous terrain. Thus, advective flux components are included into the net ecosystem exchange (NEE) budget. However, advection measurements are experimentally challenging and do not always help to solve the night flux problem of the EC technique. This study investigates alternative methods for the observation of horizontal advection, in particular horizontal concentration gradients, as well as different approaches to coordinate rotation and vertical advection. Continuous high-frequency measurements of the horizontal CO2 concentration field are employed and compared to the often used discontinuous sequential sampling. Significant differences were found in the case of 30-min mean concentration values between the conventional discontinuous sampling approach and the complete observation of the time series by continuous sampling. Estimates of vertical advection rely on accurate estimates of vertical wind velocity ( $\emph{w}$ ). Therefore, different approaches to the planar fit coordinate rotation have been investigated. Sector-wise rotation was able to eliminate directional dependencies of mean $\emph{w}$ . Furthermore, the effect of the data set length used for rotation (window length) was investigated and was found to have significant impact on estimates of vertical advection, with larger window lengths yielding about 50% larger vertical advection. A sequential planar fit with controlled window length is proposed to give reproducible results. The different approaches to the measurement and calculation of horizontal and vertical advection presented are applied to data obtained during the exchange processes in mountainous region experiment at the FLUXNET site Waldstein–Weidenbrunnen (DE-Bay). Estimates of NEE including advection are compared to NEE from turbulent and storage flux alone without advection. NEE including vertical advection with sector-wise planar fit rotation and controlled window length and including horizontal advection from continuous gradient measurements, which were comprehensively bias corrected by a new approach, did compare well with the expected night flux error, with meteorological drivers of the fluxes and with soil chamber measurements. Unrealistically large and noisy values of horizontal advection from the conventional discontinuous sampling approach, which lead to unrealistic values of NEE, could be eliminated by the alternative approaches presented. We therefore suggest the further testing of those approaches at other sites in order to improve the accuracy of advection measurements and, subsequently, estimates of NEE.  相似文献   

19.
During spring and autumn, many lakes in temperate latitudes experience intensive convective mixing in the vertical, which leads to almost isothermal conditions with depth. Thus the regime of turbulence appears to be similar with that characteristic of convective boundary layers in the atmosphere. In the present paper a simple analytical approach, based on boundary-layer theory, is applied to convective conditions in lakes. The aims of the paper are firstly to analyze in detail the temperature distribution during these periods, and secondly to investigate the current system, created by the horizontal temperature gradient and wind action. For these purposes, simple analytical solutions for the current velocities are derived under the assumption of depth-constant temperatures. The density-induced current velocities are shown to be small, in the order of a few mm/sec. The analytical model of wind-driven currents is compared with field data. The solution is in good qualitative agreement with observed current velocities under the condition that the wind field is steady for a relatively long time and that residual effects from former wind events are negligible.The effect of the current system on an approximately depth-constant temperature distribution is then checked by using the obtained current velocity fields in the heat transfer equation and deriving an analytical solution for the corrected temperature field. These temperature corrections are shown to be small, which indicates that it is reasonable to describe the temperature distribution with vertical isotherms.Notation T temperature - t time - x, y, z cartesian coordinates - molecular viscosity - h , v horizontal and vertical turbulent viscosity - K h ,K v horizontal and vertical turbulent conductivity - Q heat flux through the water surface - D depth - u, v, w average current velocity components inx, y andz directions - f Coriolis parameter - p pressure - density - g gravity acceleration - a constant in the freshwater state equation - h s deviation from the average water surface elevation - L *,H * length and depth scale - U *,W * horizontal and vertical velocity scale - T temperature difference scale - bottom slope - u * friction velocity at the water surface - von Karman constant - L Monin-Obukhov length scale - buoyancy parameter - l turbulence length scale - C 1,C 2,C 3 dimensionless constants in the expressions for the vertical turbulent viscosity - , dimensionless vertical coordinate and dimensionless local depth - angle between surface stress direction andx-axis - T bx ,T by bottom stress components - c bottom drag coefficient  相似文献   

20.
The existence of universal power laws at low wavenumbers (K) in the energy spectrum (Eu) of the turbulent longitudinal velocity (u) is examined theoretically and experimentally for the near-neutral atmospheric surface layer. Newly derived power-law solutions to Tchen's approximate integral spectral budget equation are tested for strong- and weak-interaction cases between the mean flow and turbulent vorticity fields. To verify whether these solutions reproduce the measured Eu at low wavenumbers, velocity measurements were collected in the dynamic sublayer of the atmosphere at three sites and in the inner region of a laboratory open channel. The atmospheric surface layer measurements were carried out using triaxial sonic anemometers over tall corn, short grass, and smooth desert-like sandy soil. The open channel measurements were performed using a two-dimensional boundary-layer probe above a smooth stainless steel bed. Comparisons between the proposed analytical solution for Eu, the dimensional analysis by Kader and Yaglom, and the measured Haar wavelet Eu spectra are presented. It is shown that when strong interaction between the mean flow and turbulent vorticity field occurs, wavelet spectra measurements, predictions by the analytical solution, and predictions by the dimensional analysis of Kader-Yaglom (KY) are all in good agreement and confirm the existence of a -1 power law in Eu(= Cuuu2 * K-1, where Cuu is a constant and u* is the friction velocity). The normalized upper wavenumber limit of the -1 power law (Kz = 1, where z is the height above the zero-plane displacement) is estimated using two separate approaches and compared to the open channel and atmospheric surface-layer measurements. It is demonstrated that the measured upper wavenumber limit is consistent with Tchen's budget but not with the KY assumptions. The constraints as to whether the mean flow and turbulent vorticity strongly interact are considered using a proposed analysis by Panchev. It is demonstrated that the arguments by Panchev cannot be consistent with surface-layer turbulence. Using dimensional analysis and Heisenberg's turbulent viscosity model, new constraints are proposed. The new constraints agree with the open channel and atmospheric surface-layer measurements, Townsend's inactive eddy motion hypothesis, and the Perry et al. analysis.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号