Search

A framework to describe transport and deposition of ultrafine particles (UFP) within forests using large eddy simulation (LES) is presented. Comparison with measurements collected within and above a Scots pine stand in Southern Finland are used to explore the plausibility of the simulations. The numerical model is then employed to quantify the effects of canopy morphology (leaf area index and leaf area density), turbulence intensity, and particle size on the partitioning between upper canopy and subcanopy deposition and the overall deposition velocities. Results show a complex interplay between canopy morphology and turbulence, which is reflected on the particle flux profiles within the canopy. However, mean particle concentration profiles, total deposition, and deposition velocities at the canopy top are insensitive to the leaf area density profile but show dependence mostly on leaf area index, turbulence levels, and particle size. Finally, with the goal of understanding the sensitivity of the deposition velocity to all these parameters, an analytical model is developed that shows good agreement with LES results (within ± 20 % ) for all conditions simulated here.

In this study a stochastic wall model based on the “large-eddy” version of the One-Dimensional Turbulence (ODT) model was developed for Large-Eddy Simulation (LES) of smooth and rough channel flows, with the primary goal of providing a refined turbulent flow field near the wall. This LES-ODT coupling was tested with the dynamic Smagorinsky and the scale-dependent Lagrangian dynamic subgrid-scale models. When compared to the same LES with a wall model based on a local law-of-the-wall, LES-ODT improved the one-dimensional energy spectra for all three velocity components close to the wall for both subgrid-scale models tested. More importantly, improving the LES wall model had a more positive effect in the near-wall spectra than improving the subgrid-scale model from the traditional dynamic to the scale-dependent Lagrangian dynamic model. For smooth channels, LES-ODT results compared well with DNS of R e λ = 590 and 5200; however, the variance modeled by the ODT presents discrepancies for all three velocity components, an issue inherent to ODT. Finally, the simulation of a channel flow with additional roughness modeled by a drag force was compared to data of atmospheric flow through a maize field, providing evidence of the potential for this approach to directly simulate complex near-wall phenomena. Given its high computational cost, the main use of the LES-ODT coupling is in studies that require a refinement of the near-wall region without the need to refine the entire LES domain.

Author(s): Chamecki, M; Chor, T; Yang, D; Meneveau, C | Abstract: Material transport in the ocean mixed layer (OML) is an important component of natural processes such as gas and nutrient exchanges. It is also important in the context of pollution (oil droplets, microplastics, etc.). Observational studies of small-scale three-dimensional turbulence in the OML are difficult, especially if one aims at a systematic coverage of relevant parameters and their effects, under controlled conditions. Numerical studies are also challenging due to the large-scale separation between the physical processes dominating transport in the horizontal and vertical directions. Despite this difficulty, the application of large eddy simulation (LES) to study OML turbulence and, more specifically, its effects on material transport has resulted in major advances in the field in recent years. In this paper we review the use of LES to study material transport within the OML and then summarize and synthesize the advances it has enabled in the past decade or so. In the first part we describe the LES technique and the most common approaches when applying it in OML material transport investigations. In the second part we review recent results on material transport obtained using LES and comment on implications.

In the context of many applications of turbulent multi-phase flows, knowledge of the dispersed phase size distribution and its evolution is critical to predicting important macroscopic features. We develop a large eddy simulation (LES) model that can predict the turbulent transport and evolution of size distributions, for a specific subset of applications in which the dispersed phase can be assumed to consist of spherical droplets, and occurring at low volume fraction. We use a population dynamics model for polydisperse droplet distributions specifically adapted to a LES framework including a model for droplet breakup due to turbulence, neglecting coalescence consistent with the assumed small dispersed phase volume fractions. We model the number density fields using an Eulerian approach for each bin of the discretized droplet size distribution. Following earlier methods used in the Reynolds-averaged Navier–Stokes framework, the droplet breakup due to turbulent fluctuations is modelled by treating droplet–eddy collisions as in kinetic theory of gases. Existing models assume the scale of droplet–eddy collision to be in the inertial range of turbulence. In order to also model smaller droplets comparable to or smaller than the Kolmogorov scale we extend the breakup kernels using a structure function model that smoothly transitions from the inertial to the viscous range. The model includes a dimensionless coefficient that is fitted by comparing predictions in a one-dimensional version of the model with a laboratory experiment of oil droplet breakup below breaking waves. After initial comparisons of the one-dimensional model to measurements of oil droplets in an axisymmetric jet, it is then applied in a three-dimensional LES of a jet in cross-flow with large oil droplets of a single size being released at the source of the jet. We model the concentration fields using $N_{d}=15$ bins of discrete droplet sizes and solve scalar transport equations for each bin. The resulting droplet size distributions are compared with published experimental data, and good agreement for the relative size distribution is obtained. The LES results also enable us to quantify size distribution variability. We find that the probability distribution functions of key quantities such as the total surface area and the Sauter mean diameter of oil droplets are highly variable, some displaying strong non-Gaussian intermittent behaviour.

Inland deposition of sea salt aerosol (SSA) particles emitted over the ocean is studied via numerical and theoretical models. The focus is on the large particles that contribute most to the total mass deposition. Large eddy simulations of idealized sea wind are used to investigate the development of the particle plume over land for different particle sizes and to validate some of the assumptions in the theoretical model. An existing theoretical modeling framework for particle dispersion in the atmospheric boundary layer is adapted to the problem of SSA deposition and it is shown to be adequate for the large particles of interest here. The decay of monodisperse SSA particle deposition flux with distance from the shoreline is shown to have a power-law behavior far from the shoreline. A complete model for predicting mass deposition as a function of distance is formulated and shown to present reasonable agreement with existing data.

The transport in the ocean mixed layer (OML) of oil plumes originated from deepwater blowouts is studied using large eddy simulations. In particular, we focus on the effects of swell on the modulation of turbulence in the OML and its impact on oil transport. Results show that when the misalignment between the wind and the swell propagation is small, Langmuir cells develop and significantly enhance the vertical dilution of the oil plume. Conversely, when the misalignment is large, vertical dilution is suppressed when compared to the no-swell case. Due to the strong directional shear of the mean flow within the OML, plume depth significantly impacts mean transport direction. The size of oil droplets in the plume also plays an important role in vertical dilution and mean transport direction.

Large-eddy simulation (LES) has proven to be a valuable tool for high-fidelity modeling of environmental and geophysical turbulent flows. An important application of LES is to study the transport of effluents (e.g. oils from a subsea blowout) in the ocean mixed layer (OML). Oil plumes being transported in the OML experience the action of shear-generated turbulence, Langmuir circulations, Ekman transport and submesoscale quasi-geostrophic eddies. To resolve such turbulent processes, grid sizes of a few meters are desirable while horizontal domain sizes of LES are typically restricted from hundreds of meters to a few kilometers, for LES to remain practically affordable. Yet transported oil plumes evolve to large scales extending to tens or even hundreds of kilometers. In this study, the Extended Nonperiodic Domain LES for Scalar transport (ENDLESS) is proposed as a multi-scale approach to tackle this challenge while being computationally affordable. The basic idea is to simulate the shear turbulence and Langmuir circulations on a small horizontal domain with periodic boundary conditions while the resulting transport velocity field is replicated periodically following adaptively the large-scale plume as it evolves spatially towards much larger scales. This approach also permits the superposition of larger-scale quasi two-dimensional flow motions on the oil advection, allowing for coupling with regional circulation models. A validation case and two sample applications to oil plume evolution in the OML are presented in order to demonstrate key features and computational speedup associated with the ENDLESS method.

Characteristics of laboratory-scale bubble-driven buoyant plumes in a stably stratified quiescent fluid are studied using large-eddy simulation (LES). As a bubble plume entrains stratified ambient water, its net buoyancy decreases due to the increasing density difference between the entrained and ambient fluids. A large fraction of the entrained fluid eventually detrains and falls along an annular outer plume from a height of maximum rise (peel height) to a neutral buoyancy level (trap height), during which less buoyant scalars (e.g. small droplets) are trapped and dispersed horizontally, forming quasi-horizontal intrusion layers. The inner/outer double-plume structure and the peel/intrusion process are found to be more distinct for cases with small bubble rise velocity, while weak and unstable when the slip velocity is large. LES results are averaged to generate distributions of mean velocity and turbulent fluxes. These distributions provide data for assessing the performance of previously developed closures used in one-dimensional integral plume models. In particular, the various LES cases considered in this study yield consistent behaviour for the entrainment coefficients for various plume cases. Furthermore, a new continuous peeling model is derived based on the insights obtained from LES results. Comparing to previous peeling models, the new model behaves in a more self-consistent manner, and it is expected to provide more reliable performance when applied in integral plume models.

In this study, large eddy simulation (LES) is adopted to evaluate the aerodynamic and deposition effects of street trees on fine particulate matter (PM2.5) concentration within street canyons. The Extended Nonperiodic Domain LES for Scalar Transport (ENDLESS) is used to allow exploration of vegetation effects at a neighborhood scale (up to 100 canyons) while maintaining a reasonable resolution required to resolve flow patterns within each canyon. We investigate three emission scenarios: (i) only local traffic emissions within the canyon (the first canyon in the urban environment); (ii) only background pollution originating from the upwind canyons; and (iii) a combination of scenarios (i) and (ii). Numerical results show that the presence of trees has different effects on the PM2.5 level within canyons in different emission scenarios and at different spatial scales. At the street scale with only local traffic emissions, aerodynamic effect of trees results in an increase in the concentration near leeward walls and a decrease in the concentration near windward walls, which overwhelms the deposition effect. On the other hand, trees have a negligible impact on the transport of background pollution into the canyon or its distribution within the canyons. The deposition has beneficial effects that only manifest in a considerable way at the neighborhood scale. Finally, the effect of trees on a simple operational urban pollution model (OUPM) is investigated. After modifications for aerodynamic and deposition effect of trees, the predictions of the OUPM show good agreement with LES results.

A large eddy simulation is used to study the surface patterns of noninertial passive buoyant particles in the ocean mixed layer. A mechanism for preferential concentration that is independent of inertial effects is then identified and quantified as a function of the terminal rise velocity.

Once oil plumes such as those originating from underwater blowouts reach the ocean mixed layer (OML), their near-surface dispersion is influenced heavily by wind and wave-generated Langmuir turbulence. In this study, the complex oil spill dispersion process is modeled using large-eddy simulation (LES). The mean plume dispersion is characterized by performing statistical analysis of the resulting fields from the LES data. Although the instantaneous oil concentration exhibits high intermittency with complex spatial patterns such as Langmuir-induced striations, it is found that the time-averaged oil distribution can still be described quite well by smooth Gaussian-type plumes. LES results show that the competition between droplet rise velocity and vertical turbulent diffusion due to Langmuir turbulence is crucial in determining both the dilution rate and overall direction of transport of oil plumes in the OML. The smoothness of the mean plume makes it feasible to aim at modeling the oil dispersion using Reynolds-averaged type formulations, such as the K-profile parameterization (KPP) with sufficient vertical resolution to capture vertical profiles in the OML. Using LES data, we evaluate the eddy viscosity and eddy diffusivity following the KPP framework. We assess the performance of previous KPP models for pure shear turbulence and Langmuir turbulence by comparing them with the LES data. Based on the assessment a modified KPP model is proposed, which shows improved overall agreement with the LES results for both the eddy viscosity and the eddy diffusivity of the oil dispersion under a variety of flow conditions and droplet sizes.

Modelling the dispersion of small particles such as fungal spores, pollens and small seeds inside and above plant canopies is important for many applications. Transport of these particles is driven by strongly inhomogeneous and non-Gaussian turbulent flows inside the canopy roughness sublayer, the region that extends from the ground to approximately three canopy heights. A large-eddy simulation (LES) approach is refined to study particle dispersion within and above the canopy region. Effects of plant reconfiguration are parameterized through a velocity-dependent drag coefficient, which is shown to be critical for accurate reproduction of velocity statistics and mean spore concentrations. The model yields predictions of turbulence statistics that are in good agreement with measurements. This is particularly true of the stress fractions carried by strong events, as revealed by standard quadrant analysis of the resolved velocity fluctuations, which is a known weakness of earlier LES studies of canopy flow using a constant drag coefficient. Experimental data on spore dispersal inside and above a maize canopy are reproduced successfully as well. Characteristics of the particle plume are analysed using LES results, and a pre-existing theoretical framework is adapted to model particle dispersal above the canopy. The results suggest that the plume above the canopy can be approximated using a simple analytical solution if the fraction of spores that escape the canopy region is known. Source height and gravitational settling have strong effects on the plume inside the canopy region and consequently determine the escape fraction. These effects are parameterized in the theoretical model by using the escape fraction to rescale the source strength.

A numerical study of the effect of entrainment fluxes at the top of the atmospheric boundary layer (ABL) on dissimilarity between scalars within the mixed and surface layers is conducted. Simulation results clearly show that entrainment fluxes of opposite sign cause decorrelation between the scalars throughout the entire ABL. In the upper part of the mixed layer, this decorrelation is caused by changes in the covariance between the scalars and the scalar variance as well, and is distributed over the entire range of scales resolved in the simulation. Near the surface, the reduction in correlation coefficient originates from an increasing scalar variance, which is present exclusively in the large scales. These effects are noticeable on time scales of about 24 min or longer, and could be interpreted as nonstationarity for the typical 30-min periods used in surface-layer data processing. In addition, it is shown that, for the conditions studied here, the scalar correlation coefficient within the surface layer scales with the measurement height normalized by the ABL depth and not by the Obukhov length.

An analytical model for dispersion of biological particles such as spores and pollens released from finite area sources into the atmospheric boundary layer is proposed. The model is applicable to rectangular fields aligned with the mean wind and accounts for effects of turbulent dispersion in the vertical and crosswind direction, as well as depletion of airborne material due to dry deposition. Predictions for the mean particle concentration above the area source and ground deposition downwind from it are derived and validated against detailed numerical simulations of pollen dispersal using large eddy simulation. Effects of particle settling velocity and area source size and geometry on dispersal patterns are investigated. Results show that along the centerline above the source, effects of crosswind dispersion are negligible (except for very narrow area sources). Ground deposition fluxes downwind from the source are strongly affected by crosswind dispersion. Similarly to area sources infinitely long in the crosswind direction, deposition fluxes far from rectangular sources display a power-law decay with distance. The power-law exponent is larger than the one for infinitely wide fields and depends on the relative importance of gravitational settling and turbulent dispersion in the vertical and crosswind directions, but does not depend on source geometry. However, the distance from the source at which the power-law decay is first observed does depend on source geometry.

A scale-similarity model to estimate the subfilter-scale energy using the trace of the Leonard stress tensor is proposed and evaluated for large-eddy simulations of the atmospheric boundary layer (ABL). The model is derived from a stability-dependent model of the energy spectrum in the ABL, which accounts for the effects of buoyancy and mean shear as a function of z/L, the Monin–Obukhov stability variable. An a priori test using ABL turbulence data demonstrates that the model has accurate performance for dimensionless filter widths of Δ/z = 2, 1, and 0.5 for stabilities of −1 ≤ z/L ≤ 0.5, and improves considerably upon a similar model that is derived using an infinite κ −5/3 spectrum. This improvement is especially significant in the first several grid points near the surface in large-eddy simulations of the ABL, where Δ/z is necessarily large. The modelling procedure is then extended to develop a similarity model for the subfilter-scale scalar variance; it is shown to have robust performance for temperature.

This paper presents a framework to simulate pollen dispersal by the wind based on the large eddy simulation (LES) technique. Important phenomena such as the pollen emission by the plants and the ground deposition are parameterized by the lower boundary condition. The numerical model is validated against previously published experiments of point source releases of glass beads and pollen grains in the atmospheric boundary layer. The numerical model is used together with experimental data of pollen emission and downwind deposition from a natural field obtained near Washington, DC, in the summer of 2006. The combined analysis of experimental and numerical data allows to elucidate the emission/transport/deposition process in considerable detail. In particular, the relative fractions of pollen deposited inside the source field and airborne at the edge of the field can be quantified. The use of LES allows quantification of important intermittent deposition events far from the source field.

Pseudospectral methods are frequently used in the horizontal directions in large-eddy simulation of atmospheric flows. However, the same approach often creates unphysical oscillations for scalar fields if there are horizontal heterogeneities in the sources and/or sinks, as is usual in air pollution problems. A hybrid approach is developed to combine the use of pseudospectral representation of the velocity field and bounded finite-volumes for the scalar concentration. An interpolation scheme that yields a divergence-free interpolated velocity field is derived and implemented, and its importance is illustrated by two sample applications.

The mean drag–wind relationship is usually investigated assuming that field data are representative of spatially-averaged metrics of statistically stationary flow within and above a horizontally homogeneous canopy. Even if these conditions are satisfied, large-eddy simulation (LES) data suggest two major issues in the analysis of observational data. Firstly, the streamwise mean pressure gradient is usually neglected in the analysis of data from terrestrial canopies, which compromises the estimates of mean canopy drag and provides misleading information for the dependence of local mean drag coefficients on local velocity scales. Secondly, no standard approach has been proposed to investigate the instantaneous drag–wind relationship, a critical component of canopy representation in LES. Here, a practical approach is proposed to fit the streamwise mean pressure gradient using observed profiles of the mean vertical momentum flux within the canopy. Inclusion of the fitted mean pressure gradient enables reliable estimates of the mean drag–wind relationship. LES data show that a local mean drag coefficient that characterizes the relationship between mean canopy drag and the velocity scale associated with total kinetic energy can be used to identify the dependence of the local instantaneous drag coefficient on instantaneous velocity. Iterative approaches are proposed to fit specific models of velocity-dependent instantaneous drag coefficients that represent the effects of viscous drag and the reconfiguration of flexible canopy elements. LES data are used to verify the assumptions and algorithms employed by these new approaches. The relationship between mean canopy drag and mean velocity, which is needed in models based on the Reynolds-averaged Navier-Stokes equations, is parametrized to account for both the dependence on velocity and the contribution from velocity variances. Finally, velocity-dependent drag coefficients lead to significant variations of the calculated displacement height and roughness length with wind speed., National Science Foundation (U.S.) (Grant AGS1005363)

Pollen dispersion by the wind is a subject of interest in many fields which has gained renewed attention in the recent years due to the increasing level of landscape fragmentation and the introduction of genetically modified crops. Accurate quantitative estimates of the three-dimensional plume of airborne pollen grains and consequent ground deposition from an emitting field are usually required. A large number of field experiments has been conducted to investigate the problem. However, the results are expected to depend upon several factors, such as aerodynamic properties of the pollen grains, geometric properties of the field, and local vegetation and atmospheric conditions, limiting the validity of experimental results to specific field conditions. Up to date, simulations of pollen dispersal were restricted to the use of Lagrangian stochastic models [6] and two-dimensional RANS models [4]. In this work we describe tools and frameworks to simulate pollen dispersal using large eddy simulation (LES).

The performances of four models for the subgrid-scale heat flux under conditions of poor resolution typical of large-eddy simulation of atmospheric boundary layer flows are compared using observational data. It is argued that a key feature of a numerically stable model is to accurately predict the probability density function of the dissipation of resolved temperature variance (or at least not overpredict the amount of backscatter). The results show that the nonlinear model yields excessive backscatter in agreement with numerical instabilities observed in a posteriori implementations. It is also observed that the Daly–Harlow–Smagorinsky model performs much better, despite having a similar structure. The source of the excessive backscatter in the nonlinear model is tracked to the presence of the rotation component in the tensor eddy diffusivity. A modified version of the Daly–Harlow model is proposed on the basis of a closure for the subgrid-scale stress tensor using the nonlinear model after elimination o...

Phenomena such as large-scale shear, buoyancy, and the proximity to the ground surface significantly affect interactions among scales in atmospheric boundary layer turbulent flows. Hence, these phenomena impact parameters that enter subgrid-scale (SGS) parameterizations used in large eddy simulations (LES) of the atmospheric boundary layer. The effects of these phenomena upon SGS parameters have, to date, been studied mostly as functions of the global state of the flow. For instance, the Smagorinsky coefficient has been measured as a function of the mean shear and stability condition of the atmosphere as determined from the average surface heat and momentum fluxes. However, in LES the global average field values are often difficult to determine a priori and the SGS parameters ideally must be expressed as a function of local flow variables that characterize the instantaneous flow phenomena. With the goal of improving the Smagorinsky closure, in this study several dimensionless parameters characterizing the local structure and important dynamical characteristics of the flow are defined. These local parameters include enstrophy, vortex stretching, self-amplification of strain rate, and normalized temperature gradient and all are defined in such a way that they remain bounded under all circumstances. The dependence of the Smagorinsky coefficient on these local parameters is studied a priori from field data measured in the atmospheric surface layer and, as a reference point, from direct numerical simulation of neutrally buoyant, isotropic turbulence. To capture the local effects in a statistically meaningful fashion, conditional averaging is used. Results show various important and interrelated trends, such as significant increases of the coefficient in regions of large strain-rate self-amplification and vortex stretching. Results also show that the joint dependence on the parameters is rather complicated and cannot be described by products of functions that depend on single parameters. Dependence on locally defined parameters is expected to improve the SGS model by sensitizing it to local flow conditions and by enabling possible generalizations of the dynamic model based on conditional averaging.