A grid nesting method for large-eddy simulation of planetary boundary-layer flows

A method for performing nested grid calculations with a large-eddy simulation code is described. A common numerical method is used for all meshes, and the grid architecture consists of a single outer or coarse grid, and nested or fine grids, which overlap in some common region. Inter-grid communication matches the velocity, pressure and potential temperature fields in the overlap region. Resolved and sub-grid scale (SGS) turbulent fluxes and kinetic energy on the fine grid are averaged to the coarse grid using a conservation rule equivalent to Germano's identity used to develop dynamic SGS models.Simulations of a slightly convective, strong shear planetary boundary layer were carried out with varying surface-layer resolutions. Grid refinements in the (x, y, z) directions of up to (5, 5, 2) times were employed. Two-way interaction solutions on the coarse and fine meshes are successfully matched in the overlap region on an instantaneous basis, and the turbulent motions on the fine grid blend smoothly into the coarse grid across the grid interface. With surface-layer grid nesting, significant increases in resolved eddy fluxes and variances are found. The energy-scale content of the vertical velocity, and hence vertical turbulent fluxes, appear to be most influenced by increased grid resolution. Vertical velocity spectra show that the dominant scale shifts towards higher wavenumbers (smaller scales) and the magnitude of the peak energy is increased by more than a factor of 3 with finer resolution. Outside of the nested region the average heat and momentum fluxes and spectra are slightly influenced by the fine resolution in the surface layer. From these results we conclude that fine resolution is required to resolve the details of the turbulent motions in the surface layer. At the same time, however, increased resolution in the surface layer does not appreciably alter the ensemble statistics of the resolved and SGS motions outside of the nested region.     

Large-Eddy Simulation for the Mechanism of Pollutant Removal from a Two-Dimensional Street Canyon

Large-eddy simulation (LES) is conducted to investigate the mechanism of pollutant removal from a two-dimensional street canyon with a building-height to street-width (aspect) ratio of 1. A pollutant is released as a ground-level line source at the centre of the canyon floor. The mean velocities, turbulent fluctuations, and mean pollutant concentration estimated by LES are in good agreement with those obtained by wind-tunnel experiments. Pollutant removal from the canyon is mainly determined by turbulent motions, except in the adjacent area to the windward wall. The turbulent motions are composed of small vortices and small-scale coherent structures of low-momentum fluid generated close to the plane of the roof. Although both small vortices and small-scale coherent structures affect pollutant removal, the pollutant is largely emitted from the canyon by ejection of low-momentum fluid when the small-scale coherent structures appear just above the canyon where the pollutant is retained. Large-scale coherent structures also develop above the canyon, but they do not always affect pollutant removal.     

Characteristics of turbulent structures in the unstable atmospheric surface layer

An atmospheric surface-layer (ASL) experiment conducted at a meteorological site in the Oostelijk-Flevoland polder of the Netherlands is described. Turbulent fluctuations of wind velocity, air temperature and static pressure were measured, using three 10 m towers.Simultaneous turbulent signals at several heights on the towers were used to investigate the properties of the turbulent structures which contribute most significantly to the turbulent vertical transports in the unstable ASL. These turbulent structures produce between 30 and 50% of the mean turbulent vertical transport of horizontal alongwind momentum and they contribute to between 40 and 50% of the mean turbulent vertical heat transport; in both cases this occurs during 15 to 20% of the total observation time.The translation speed of the turbulent structures equals the wind speed averaged over the depth of the ASL, which scales on the surface friction velocity. The inclination angle of the temperature interface at the upstream edge of the turbulent structures to the surface is significantly smaller than that of the internal shear layer, which is associated with the temperature interface. The turbulent structures in the unstable ASL are determined by a large-scale temperature field: Convective motions, which encompass the whole depth of the planetary boundary layer (PBL), penetrate into the ASL. The curvature of the vertical profile of mean horizontal alongwind velocity forces the alignment of the convective cells in the flow direction (Kuettner, 1971), which have an average length of several hundreds of metres and an average width of a few tens of metres. This mechanism leads to the formation of turbulent structures, which extend throughout the depth of the ASL.     

Large-Scale Turbulence Structures and Their Contributions to the Momentum Flux and Turbulence in the Near-Neutral Atmospheric Boundary Layer Observed from a 213-m Tall Meteorological Tower

Large-scale turbulence structures in the near-neutral atmospheric boundary layer (ABL) are investigated on the basis of observations made from the 213-m tall meteorological tower at Tsukuba, Japan. Vertical profiles of wind speed and turbulent fluxes in the ABL were obtained with sonic anemometer-thermometers at six levels of the tower. From the archived data, 31 near-neutral cases are selected for the analysis of turbulence structures. For the typical case, event detection by the integral wavelet transform with a large time scale (180 s) from the streamwise velocity component (u) at the highest level (200 m) reveals a descending high-speed structure with a time scale of approximately 100 s (a spatial scale of 1 km at the 200-m height). By applying the wavelet transform to the u velocity component at each level, the intermittent appearance of large-scale high-speed structures extending also in the vertical is detected. These structures usually make a large contribution to the downward momentum transfer and induce the enhancement of turbulent kinetic energy. This behaviour is like that of “active” turbulent motions. From the analysis of the two-point space–time correlation of wavelet coefficients for the u velocity component, the vertical extent and the downward influence of large-scale structures are examined. Large fluctuations in the large-scale range (wavelet variance at the selected time scale) at the 200-m level tend to induce the large correlation between the higher and lower levels.     

An analysis of wind velocity fluctuations in the atmospheric surface layer using an orthonormal wavelet transform

A wavelet analysis can supply information of both the location (time) and the scale of fluctuations. This method is applied to the fluctuations of the natural wind and the turbulent transport of momentum in the atmospheric surface layer. The shapes of both the wavelet spectra and the Fourier spectra of the three components of the wind velocity fluctuations are similar to each other. The quadrant representation of momentum transport shows the scale difference of the transport. The large-scale fluctuations mainly contribute to the downward transport of momentum.     

大涡模式垂直分辨率在夹卷及示踪物垂直传输模拟结果分析中的作用

基于敦煌野外观测资料和大涡模式,研究了垂直方向不同尺度湍涡对夹卷及示踪物垂直传输的影响,明确了模式垂直分辨率在模拟结果分析中的作用。结果表明:垂直方向上小尺度湍涡对夹卷作用贡献更大,小尺度湍涡较多时夹卷层相对更暖,而夹卷层厚度、夹卷强度和风速变化受垂直方向湍涡尺度影响较小。当垂直分辨率为50 m时,越往夹卷层上部,上升气流和下沉气流分布较多且强度较大;分辨率为10、20和30 m时,夹卷层各高度垂直速度、位温和示踪物浓度分布较接近。另外,垂直方向湍涡尺度对示踪物垂直传输高度影响不大,而对示踪物的空间分布有一定作用。当大尺度湍涡较多且强度较强时,越有利于将高浓度的示踪物向上传输。综合考虑到模式采用较高分辨率模拟时产生的噪音及计算时间等问题,认为模式采用30 m的垂直分辨率,既能较好地模拟出夹卷层平均结构特征,又能模拟出夹卷层湍流的精细分布,是较为理想的选择。     

Effect of Incoming Turbulent Structure on Pollutant Removal from Two-Dimensional Street Canyon

Large-eddy simulations are conducted to investigate the effects of the incoming turbulent structure of the flow on pollutant removal from an ideal canyon. The target canyon is a two-dimensional street canyon with an aspect ratio of 1.0 (building height to street width). Three turbulent flows upwind of the street canyon are generated by using different block configurations, and a tracer gas is released as a ground-level line source at the centre of the canyon floor. Mean velocity profiles for the three flows are similar, except near the roof. However, the root-mean-square values of the velocity fluctuations and the Reynolds shear stress increase with the friction velocity of the incoming turbulent flow. The spatially-averaged concentration within the canyon decreases with increasing friction velocity. Coherent structures of low-momentum fluid, generated above the upwind block configurations, contribute to pollutant removal, and the amount of pollutant removal is directly related to the size of the coherent structure.     

Characteristics of vertical turbulent velocities in the urban convective boundary layer

Turbulence measurements of the vertical velocity component were obtained by an instrumented aircraft under fair weather conditions in the St. Louis, Missouri, metropolitan area. Time series of vertical velocity fluctuations from horizontal flight segments made in the lower part of and near the middle of the convective boundary layer (CBL) over the urban area and surrounding region were subjected to various statistical and objective analyses. Higher order vertical velocity moments, and positive and negative velocity statistics, were computed. The horizontal dimensions of updrafts and downdrafts, and related properties of these turbulent eddies were derived by conditional sampling analysis. Emphasis is on a comparison of the results from urban and selected rural measurements from the lower part of the CBL.The vertical velocity probability density distribution for each segment was positively skewed and the mode was negative. The means and standard deviations of positive and negative velocity fluctuations were greater over the urban area. The urban vertical velocity variance was 50% greater than rural values, and power spectra revealed greater production of vertical turbulent energy in the urban area over a wide frequency range.The mean and maximum widths of downdrafts were generally larger than the corresponding values for updrafts. Differences between urban and rural eddy sizes were not statistically significant. The widths of the largest updraft and downdraft are comparable to the boundary-layer depth, Z _i, and the mean value of the ratio of spectral peak wavelength to Z _iwas about 1.3 and 1.1 for urban and rural areas, respectively. Convective similarity scaling parameters appeared to order both the urban and rural measurements.On assignment from the National Oceanic and Atmospheric Administration, U.S. Dept. of Commerce.     

Fluctuations of line integrated concentrations across plume diffusion in grid generated turbulence and in shear flows

The fluctuations of the instantaneous values of line integrated concentrations across plumes from point sources diffusing in turbulent shear flows, and in grid generated turbulence, have been studied experimentally using a fast response system which measured the attenuation of the intensity of an infrared beam crossing the plume. Analysis of the measurements show that the dimensionless statistical properties of the fluctuations at different distances from the source at each flow are approximately similar, in the sense that they depend primarily on the relative off-center location of the line of integration and almost independent of the distance from the source and the nature of the turbulence in the flows, as long as the characteristic length of the mean plume is not large compared to the size of the large eddies. The characteristic time of the fluctuations, on the other hand, was found to grow with the distance from the source and the autocorrelations of the fluctuations, particularly in the case of a plume diffusing in grid generated turbulence, were it found to be proportional to the lateral size of the mean plume. A—5/3 decay law of the power spectrum of the fluctuations was observed in the low frequency range which corresponds to the scale of the large eddies. The decay of the fluctuations caused by smaller eddies was much faster, as expected.     

A spectral and lag-correlation analysis of turbulence in a deciduous forest canopy

The processes influencing turbulence in a deciduous forest and the relevant length and time scales are investigated with spectral and cross-correlation analysis. Wind velocity power spectra were computed from three-dimensional wind velocity measurements made at six levels inside the plant canopy and at one level above the canopy. Velocity spectra measured within the plant canopy differ from those measured in the surface boundary layer. Noted features associated with the within-canopy turbulence spectra are: (a) power spectra measured in the canopy crown peak at higher wavenumbers than do those measured in the subcanopy trunkspace and above the canopy; (b) peak spectral values collapse to a relatively universal value when scaled according to a non-dimensional frequency comprised of the product of the natural frequency and the Eulerian time scale for vertical velocity; (c) at wavenumbers exceeding the spectral peak, the slopes of the power spectra are more negative than those observed in the surface boundary layer; (d) Eulerian length scales decrease with depth into the canopy crown, then increase with further depth into the canopy; (e) turbulent events below crown closure are more correlated with turbulent events above the canopy than are those occurring in the canopy crown; and (f) Taylor's frozen eddy hypothesis is not valid in a plant canopy. Interactions between plant elements and the mean wind and turbulence alter the processes that produce, transport and remove turbulent kinetic energy and account for the noted observations.     

Simplified expressions for adjusting higher-order turbulent statistics obtained from open path gas analyzers

When density fluctuations of scalars such as CO₂ are measured with open-path gas analyzers, the measured vertical turbulent flux must be adjusted to take into account fluctuations induced by ‘external effects’ such as temperature and water vapour. These adjustments are needed to separate the effects of surface fluxes responsible for ‘natural’ fluctuations in CO₂ concentration from these external effects. Analogous to vertical fluxes, simplified expressions for separating the ‘external effects’ from higher-order scalar density turbulence statistics are derived. The level of complexity in terms of input to these expressions are analogous to that of the Webb–Pearman–Leuning (WPL), and are shown to be consistent with the conservation of dry air. It is demonstrated that both higher-order turbulent moments such as the scalar variances, the mixed velocity-scalar covariances, and the two-scalar covariance require significant adjustments due to ‘external effects’. The impact of these adjustments on the turbulent CO₂ spectra, probability density function, and dimensionless similarity functions derived from flux-variance relationships are also discussed.     

Power Spectra and Cospectra for Wind and Scalars in a Disturbed Surface Layer at the Base of an Advective Inversion

This paper reports power spectra and cospectra of windspeed and several scalars measured at two heights nearthe base of an advective inversion. The inversion hadformed over a paddy field downwind of an extensive dryregion. Winds over the paddy field were variable instrength and direction, as a result of convectivemotions in the atmospheric boundary layer passing overfrom the dry region upwind. Fetch over the rice waslarge enough that advective effects on the transportprocesses were small at the upper level and negligibleat the lower level. Results from the lower level areinterpreted in terms of a horizontally homogeneous,but disturbed, surface layer.Power spectra of longitudinal and lateral velocitywere substantially enhanced at low frequencies. Theresulting vertical motions added only a small amountto the spectrum of vertical velocity but this stronglyaffected scalar power spectra and cospectra. Thesewere all substantially enhanced over a range of lowfrequencies. We also found that differences in lowerboundary conditions cause differences among scalarspectra at low frequencies.Our analysis shows that the spectra and cospectra havethree components, characterized by different scalingregimes. We call these the ILS (inner-layer scaling),OLS (outer-layer scaling) and CS (combined scaling)components. Of these, the CS component had notpreviously been identified. We identify CS componentsof spectra by their independence of height andfrequency. Spectra with these characteristics had beenpredicted by Kader and Yaglom for a layer of theatmosphere where spectral matching between ILS and OLSwas proposed. However, we find that the velocity andscalar scales used by Kader and Yaglom do not fit ourresults and that their concept of a matching layer isincompatible with our application. An alternativebasis for this behaviour and alternative scales areproposed.We compare our decomposition of spectra into ILS, CSand OLS components with an extended form of Townsend'shypothesis, in which wind and scalar fluctuations aredivided into active and inactive components. Wefind the schemes are compatible if we identify all OLSspectral components as inactive, and all CS and ILScomponents as active.By extending the implications of our results toordinary unstable daytime conditions,we predict that classical Monin–Obukhovsimilarity theory should be modified. We find that theheight of the convective boundary layer is animportant parameter when describing transportprocesses near the ground, and that the scalar scalein the ILS part of the spectrum, which includes theinertial subrange, is proportional to observationheight times the local mean scalar gradient, and notthe Monin–Obukhov scalar scale parameter. The formerdepends on two stability parameters: the Monin–Obukhovstability parameter and the ratio of the inner-layerand outer-layer velocity scales. The outer-layer scalecan reflect disturbances by topographically-inducededdying as well as by convective motions.     

Turbulent Structures in a Pine Forest with a Deep and Sparse Trunk Space: Stand and Edge Regions

Forested landscapes often exhibit large spatial variability in vertical and horizontal foliage distributions. This variability may affect canopy-atmosphere exchanges through its action on the development of turbulent structures. Here we investigate in neutral stratification the turbulent structures encountered in a maritime pine forest characterized by a high, dense foliated layer associated with a deep and sparse trunk space. Both stand and edge regions are considered. In situ measurements and the results of large-eddy simulations are used and analyzed together. In stand conditions, far from the edge, canopy-top structures appear strongly damped by the dense crown layer. Turbulent wind fluctuations within the trunk space, where the momentum flux vanishes, are closely related to these canopy-top structures through pressure diffusion. Consequently, autocorrelation and spectral analyses are not quite appropriate to characterize the vertical scale of coherent structures in this type of canopy, as pressure diffusion enhances the actual scale of structures. At frequencies higher than those associated with canopy-top structures, wind fluctuations related to wake structures developing behind tree stems are observed within the trunk space. They manifest themselves in wind velocity spectra as secondary peaks in the inertial subrange region, confirming the hypothesis of spectral short-cuts in vegetation canopies. In the edge region specific turbulent structures develop just below the crown layer, in addition to canopy-top structures. They are generated by the wind shear induced by the sub-canopy wind jet that forms at the edge. These structures provide a momentum exchange mechanism similar to that observed at the canopy top but in the opposite direction and with a lower magnitude. They may develop as in plane mixing-layer flows, with some perturbations induced by canopy-top structures. Wake structures are also observed within the trunk space in the edge region.     

General motors sulfate dispersion experiment

The mean wind field and the wind fluctuations observed near a controlled, heavily traveled level roadway were investigated. The mean wind field was studied in terms of velocity defects relative to the ambient mean wind. The wind fluctuations were studied in terms of energy spectra, velocity cross-spectra, excesses of velocity variances and cross-covariances over their ambient values. Two-point cross-spectra and correlations were also used to study the propagation of wind fluctuations. The influence of traffic on the ambient wind field extended vertically up to at least 4.5m above ground, and horizontally to at least 30 m downwind of the road. The vertical extension was more gradual at high cross-road winds than at low cross-road winds. The wind tended to move up as it approached the road and down as it left the road, except when the cross-road wind was low, where the downward motion was reduced or reversed. This is likely due to the buoyancy effect of the heated exhaust. The traffic effects propagated outward at about 1 m s^-1. Thus at cross-road winds of less than 1 ms^-1, the influence of traffic could be observed upwind of the road. When intense shear existed between the local wind field and an opposing traffic direction, large turbulent energy production resulted, generating large eddies, in additions to the high frequency (> 0.1 Hz) eddies generated by vehicular motions in the absence of intense shear. The large eddies also increased the mixing volume, so that when the ambient wind opposed the traffic direction on the upwind lanes, the pollutant concentration upwind of the road could be higher than that downwind. While the cross-covariance excesses generally agreed with the gradient transport assumption, disagreements also existed, which may not be surprising in view of the participation of large eddies in the momentum transport.     

Propagation and Effects of a Mesoscale Gravity Wave Over a Weakly-Stratified Nocturnal Boundary Layer During the SABLES2006 Field Campaign

During the SABLES2006 (Stable Atmospheric Boundary Layer Experiment in Spain 2006) field campaign, a gravity-wave episode was observed on the night of July 11 by the microbarometers deployed at the surface and on the 100-m tower. The high-amplitude, low-frequency periodic pressure fluctuations were very well correlated with the wind speed and direction. Data from neighbouring automatic stations showed that the gravity wave was not local, but long-lived and mesoscale. The propagation of the wave over the experimental site had significant effects on the structure of the weakly-stratified nocturnal boundary layer that developed that night: the stability increased, turbulent vertical motions were suppressed, the nocturnal low-level jet was disrupted, and periodic temperature fluctuations of amplitude up to 3–4 K were observed. In this work we analyse the different available data sources (tower data, RASS-SODAR, microbarometric, satellite imagery, automatic stations) to describe the phenomena in depth and to find a suitable explanation for the generation and propagation of the wave. The linear wave theory explains remarkably well most of the observations, and the wave parameters could be estimated by applying a wavelet-based technique to surface microbarometric measurements. We also analyse the vertical structure of the wave and find wave ducting conditions above the surface. Finally, by means of the multi-resolution flux decomposition, we analyse in detail the changes in vertical turbulent fluxes and the spectra of turbulent motions produced by the interaction between the gravity wave and the local flow.     

Entrainment Processes in the Convective Boundary Layer with Varying Wind Shear

Large-eddy simulations (LES) are performed to investigate the entrainment andthe structure of the inversion layer of the convective boundary layer (CBL) withvarying wind shears. Three CBLs are generated with the constant surface kinematicheat flux of 0.05 K m s^-1 and varying geostrophic wind speeds from 5 to 15m s^-1. Heat flux profiles show that the maximum entrainment heat flux as afraction of the surface heat flux increases from 0.13 to 0.30 in magnitude withincreasing wind shear. The thickness of the entrainment layer, relative to the depthof the well-mixed layer, increases substantially from 0.36 to 0.73 with increasingwind shear. The identification of vortices and extensive flow visualizations nearthe entrainment layer show that concentrated vortices perpendicular to the meanboundary-layer wind direction are identified in the capping inversion layer for thecase of strong wind shear. These vortices are found to develop along the mean winddirections over strong updrafts, which are generated by convective rolls and to appearas large-scale wavy motions similar to billows generated by the Kelvin–Helmholtzinstability. Quadrant analysis of the heat flux shows that in the case of strong windshear, large fluctuations of temperature and vertical velocity generated by largeamplitude wavy motions result in greater heat flux at each quadrant than that inthe weak wind shear case.     

Comparison of Kansas Data with High-resolution Large-eddy Simulation Fields

The surface layer of an atmospheric boundary layer (ABL) is most accessible to field measurements and hence its ensemble-mean structure has been well established. The Kansas field measurements were the first detailed study of this layer, providing numerous benchmark statistical profiles for a wide range of stability states. Large-eddy simulation (LES), in contrast, is most suitable for studying the mixed layer of the ABL where the energy-containing range of the vertical velocity field is well resolved. In the surface layer, typical large-eddy simulations barely resolve the energy-containing vertical-velocity fields and hence do not provide sufficient data for a detailed analysis.We carried out a nested-mesh simulation of a moderately convective ABL (-z_i/L = 8) in which the lower 6% of the boundary layer had an effective grid resolution of 512³. We analyze the LES fields above the 6th vertical grid level (z = 23 m) where the vertical velocity field has a well formed inertial subrange, for a detailed comparison with the Kansas results. Various terms in the budgets of turbulent kinetic energy, temperature variance, Reynolds stress, temperature flux, and some higher order moments are compared. The agreement is generally quite good; however, we do observe certain discrepancies, particularly in the terms involving pressure fluctuations.     

Some evidence of organized flow over natural waves

Measurements of the flow characteristics at 2 m over unobstructed wave surfaces on Lake Michigan were made using an anemometer-bivane as a velocity sensor. During one 40-min period of measurement, significant energy concentration was observed at the frequency of dominant surface waves in the vertical and cross wind spectra. Cross spectra between the surface elevation and vertical motions in the flow indicate that the surface lags the vertical motions by about 55 ° at the frequency of dominant waves.     

Three-dimensional air pollutant modeling in the lower atmosphere

A steady-state, three-dimensional turbulent diffusion equation describing the concentration distribution of an air pollutant from an elevated point source in the lower atmosphere is solved analytically. The same formulation can be used to obtain solutions from line, area or other kinds of sources. The solution is developed for the cases in which the velocity, vertical and lateral diffusivities are given by the power law. The model preserves the beauty of analytical solution without sacrificing much on the accuracy of approximating the velocity and eddy diffusivities. Methods of evaluating the parameters, which are required for the model applications, are discussed. Results indicate that the ratio of the plume width to the plume length increases with decreasing stability and with increasing source height. These consequences are in response to the variations of the size of eddies in the vertical direction.     

Effects of Fetch on Turbulent Flow and Pollutant Dispersion Within a Cubical Canopy

The effects of fetch on turbulent flow and pollutant dispersion within a canopy formed by regularly-spaced cubical objects is investigated using large-eddy simulation. Six tracer gases are simultaneously released from a ground-level continuous pollutant line source placed parallel to the spanwise axis at the first, second, third, fifth, seventh and tenth rows. Beyond the seventh row, the standard deviations of the fluctuations in the velocity components and the Reynolds shear stresses reach nearly equivalent states. Low-frequency turbulent flow is generated near the bottom surface around the first row and develops as the fetch increases. The turbulent flow eventually passes through the canopy at a near-constant interval. The mean concentration within the canopy reaches a near-constant value beyond the seventh row. In the first and second rows, narrow coherent structures frequently affect the pollutant escape from the top of the canopy. These structures increase in width as the fetch increases, and they mainly affect the removal of pollutants from the canopy.