A scale-dependent dynamic model for large-eddy simulation: application to a neutral atmospheric boundary layer

A scale-dependent dynamic subgrid-scale model for large-eddy simulation of turbulent flows is proposed. Unlike the traditional dynamic model, it does not rely on the assumption that the model coefficient is scale invariant. The model is based on a second test-filtering operation which allows us to determine from the simulation how the coefficient varies with scale. The scale-dependent model is tested in simulations of a neutral atmospheric boundary layer. In this application, near the ground the grid scale is by necessity comparable to the local integral scale (of the order of the distance to the wall). With the grid scale and/or the test-filter scale being outside the inertial range, scale invariance is broken. The results are compared with those from (a) the traditional Smagorinsky model that requires specification of the coefficient and of a wall damping function, and (b) the standard dynamic model that assumes scale invariance of the coefficient. In the near-surface region the traditional Smagorinsky and standard dynamic models are too dissipative and not dissipative enough, respectively. Simulations with the scale-dependent dynamic model yield the expected trends of the coefficient as a function of scale and give improved predictions of velocity spectra at different heights from the ground. Consistent with the improved dissipation characteristics, the scale-dependent model also yields improved mean velocity profiles

Scale dependence of subgrid-scale model coefficients: An a priori study

Dynamic subgrid-scale models require an a priori assumption about the variation in the model coefficients with filter scale. The standard dynamic model assumes independence of scale while the scale dependent model assumes power-law dependence. In this paper, we use field experimental data to investigate the dependence of model coefficients on filter scale for the Smagorinsky and the nonlinear models. The results indicate that the assumption of a power-law dependence, which is often used in scale dependent dynamic models, holds very well for the Smagorinsky model. For the nonlinear model, the power-law assumption seems less robust but still adequat

The influence of local meteorological conditions on the circadian rhythm of corn (Zea mays L.) pollen emission

Field experiments were performed to study the diurnal cycle of corn pollen emission and its relation to local meteorological conditions, including temperature, relative humidity, solar radiation, mean wind speed, and turbulence quantities. Pollen concentrations were measured from canopy height to twice this level using four Rotorod samplers located on a pole in the middle of the corn field. The measured pollen concentration at canopy height was used as a surrogate for the pollen source strength while the concentrations above canopy represented the pollen transported upwards from canopy height. At twice the canopy height, the pollen concentration decreased considerably to about 30% of the canopy height values. During the mornings, pollen was emitted in large quantities while during the afternoons, airborne pollen concentrations decreased and no significant atmospheric pollen was measured from about 2 h prior to sunset until sunrise the next morning. The actual time that airborne pollen was first recorded differed from day-to-day and depended on the time required for the anthers to dry and open, as well as there being sufficiently strong winds to entrain the pollen away from the plants. On four consecutive mornings diurnal atmospheric pollen concentration distributions were bi-modal in time. The first pollen concentration peak happened shortly before the direct irradiance peak on the anthers suggesting that direct solar irradiation might be important for drying the anthers. The subsequent dip in pollen concentration seemed to be linked to a lull in mean and turbulent wind conditions. Analysis of the vertical velocity fluctuations, sw, showed that the fraction of pollen transported upwards from canopy height increased with increasing sw. In addition, Quadrant– Hole analysis applied to the turbulence data sets of 2 days suggested that low values of ejection duration fractions were associated with low values of pollen concentration, while high values were associated with high concentration values. The diurnal pattern of ejection duration fractions was similar to the pattern of the fraction of pollen that reached twice the canopy height, indicatin

PIV Measurements in the Atmospheric Boundary Layer within and above a Mature Corn Canopy. Part I: Statistics and Energy Flux

Particle image velocimetry (PIV) measurements just within and above a mature corn canopy have been performed to clarify the small-scale spatial structure of the turbulence. The smallest resolved scales are about 15 times the Kolmogorov length scale ($\nu \approx$ 0.4 mm), the Taylor microscales are about $100\nu$, and the Taylor scale Reynolds numbers range between $R_{\lambda} =2000$ and 3000. The vertical profiles of mean flow and turbulence parameters match those found in previous studies. Frequency spectra, obtained using the data as time series, are combined with instantaneous spatial spectra to resolve more than five orders of magnitude of length scales. They display an inertial range spanning three decades. However, the small-scale turbulence in the dissipation range exhibits anisotropy at all measurement heights, in spite of apparent agreement with conditions for reaching local isotropy, following a high-Reynolds-number wind tunnel study. Directly calculated subgrid-scale (SGS) energy flux, determined by spatially filtering the PIV data, increases significantly with decreasing filter size, providing support for the existence of a spectral shortcut that bypasses the cascading process and injects energy directly into small scales. The highest measured SGS flux is about 40% of the estimated energy cascading rate as determined from a -5/3 fit to the spectra. Terms appearing in the turbulent kinetic energy budget that can be calculated from the PIV data are in agreement with previous results. Evidence of a very strong correlation between dissipation rate and out-of-plane component of the vorticity is demonstrated by a striking similarity between their time series. 1. Introductio

Evaporation from three water bodies of different sizes and climates: Measurements and scaling analysis

Evaporation from small reservoirs, wetlands, and lakes continues to be a theoretical and practical problem in surface hydrology and micrometeorology because atmospheric flows above such systems can rarely be approximated as stationary and planar-homogeneous with no mean subsidence (hereafter referred to as idealized flow state). Here, the turbulence statistics of temperature(T)and water vapor (q)most pertinent to lake evaporation measurementsover three water bodies differing in climate, thermal inertia and degree of advective conditions are explored. The three systems included Lac Le´man in Switzerland (high thermal inertia, near homogeneous conditions with no appreciable advection due to long upwind fetch), Eshkol reservoir in Israel (intermediate thermal inertia, frequent strong advective conditions) and Tilopozo wetland in Chile (low thermal inertia, frequent but moderate advection). The data analysis focused on how similarity constants for the flux-variance approach, CT/Cq, and relative transport efficiencies RwT/Rwq, are perturbed from unity with increased advection or the active role of temperature. When advection is small and thermal inertia is large, CT/Cq 1)primarily due to the active role of temperature, which is consistent with a large number of studies conducted over bare soil and vegetated surfaces. However, when advection is significantly large, then CT/Cq >1 (orRwT/Rwq < 1). When advection is moderate and thermal inertia is low, then CT/Cq �1. This latter equality, while consistent with Monin–Obukhov similarity theory (MOST), is due to the fact that advection tends to increase CT/Cq above unity while the active role of temperature tends to decrease CT/Cq below unity. A simplified scaling analysis derived from the scalar variance budget equation, explained qualitatively how advection could per¬turb MOST scaling (assumed to represent the idealized flow state)

Evaporation from a small water reservoir: Direct measurements and estimates

Summary Knowing the rate of evaporation from surface water resources such as chan¬nels and reservoirs is essential for precise management of the water balance. However, evaporation is difficult to measure experimentally over water surfaces and several tech¬niques and models have been suggested and used in the past for its determination. In this research, evaporation from a small water reservoir in northern Israel was measured and estimated using several experimental techniques and models during the rainless summer. Evaporation was measured with an eddy covariance (EC) system consisting of a three-dimensional sonic anemometer and a Krypton hygrometer. Measurements of net radia¬tion, air temperature and humidity, and water temperature enabled estimation of other energy balance components. Several models and energy balance closure were evaluated. In addition, evaporation from a class-A pan was measured at the site. EC evaporation measurements for 21 days averaged 5.48 mm day�1. Best model predictions were obtained with two combined flux-gradient and energy balance models (Penman–Mon¬teith–Unsworth and Penman–Brutsaert), which with the water heat flux term, gave sim¬ilar daily average evaporation rates, that were up to 3% smaller than the corresponding EC values. The ratio between daily pan and EC evaporation varied from 0.96 to 1.94. The bulk mass transfer coefficient was estimated using a model based on measurements of water surface temperature, evaporation rate and absolute humidity at 0.9 and 2.9 m above the water surface, and using two theoretical approaches. The bulk transfer coefficient

Baltimore Supersite: Highly time- and size-resolved concentrations of urban PM2.5 and its constituents for resolution of sources and immune responses

Protection of public health from the effects of air particulate matter (PM) requires measurements and methods that assess the PM chemical constituents, physical properties, and their sources. Sampling was conducted at three sites in the Baltimore area: a source-oriented (industrial) area in south Baltimore (FMC site), and two receptor area sites (Clifton Park and Ponca Street). FMC measurements were made for the initial 1-month of the project; Clifton measurements lasted for about 2 months, while measurements at Ponca Street lasted for about 9.5 months. Pollutant samples were collected at intervals ranging from 5 min to 1 h using semi-continuous monitors for PM2.5 mass, sulfate, nitrate, elemental and organic carbon, particle number size distributions (10–20,000 nm), CO, NOx, O3, 11 metals, and mass spectra of individual particles, throughout the project. In addition to standard meteorological measurements, a 3D-sonic anemometer and a LIDAR system were operated during selected periods as were a rotating drum impactor with 3- to 6-h resolution and a filter/PUF sampler for 3-h measurements of organic compounds. Standard speciation and FRM mass measurements were also made. This report describes the types of measurements that were made at the various sites of the Baltimore Supersite program as well as presents the summary statistics for some of the PM measurements that have been made. The measurements of aerosol mass, major components, and size distribution data for the three sites are compared. Results show comparable PM concentrations at Ponca Street and Clifton Park. Increased variability was observed at Ponca Street