Effect of Roughness on Surface Boundary Conditions for Large-Eddy Simulation

An important parameterization in large-eddy simulations (LESs) of high- Reynolds-number boundary layers, such as the atmospheric boundary layer, is the specification of the surface boundary condition. Typical boundary conditions compute the fluctuating surface shear stress as a function of the resolved (filtered) velocity at the lowest grid points based on similarity theory. However, these approaches are questionable because they use instantaneous (filtered) variables, while similarity theory is only valid for mean quantities. Three of these formulations are implemented in simulations of a neutral atmospheric boundary layer with different aerodynamic surface roughness. Our results show unrealistic influence of surface roughness on the mean profile, variance and spectra of the resolved velocity near the ground, in contradiction of similarity theory. In addition to similarity-based surface boundary conditions, a recent model developed from an a priori experimental study is tested and it is shown to yield more realistic independence of the results to changes in surface roughness. The optimum value of the model parameter found in our simulations matches well the value reported in the a priori wind-tunnel study.     

A Variable Mesh Spacing for Large-Eddy Simulation Models in the Convective Boundary Layer

A variable vertical mesh spacing for large-eddy simulation (LES) models in a convective boundary layer (CBL) is proposed. The argument is based on the fact that in the vertical direction the turbulence near the surface in a CBL is inhomogeneous and therefore the subfilter-scale effects depend on the relative location between the spectral peak of the vertical velocity and the filter cut-off wavelength. From the physical point of view, this lack of homogeneity makes the vertical mesh spacing the principal length scale and, as a consequence, the LES filter cut-off wavenumber is expressed in terms of this characteristic length scale. Assuming that the inertial subrange initial frequency is equal to the LES filter cut-off frequency and employing fitting expressions that describe the observed convective turbulent energy one-dimensional spectra, it is feasible to derive a relation to calculate the variable vertical mesh spacing. The incorporation of this variable vertical grid within a LES model shows that both the mean quantities (and their gradients) and the turbulent statistics quantities are well described near to the ground level, where the LES predictions are known to be a challenging task.     

Large-Eddy Simulation Of The Stably Stratified Planetary Boundary Layer

In this work, we study the characteristics of a stably stratifiedatmospheric boundary layer using large-eddy simulation (LES).In order to simulate the stable planetary boundary layer, wedeveloped a modified version of the two-part subgrid-scalemodel of Sullivan et al. This improved version of themodel is used to simulate a highly cooled yet fairly windy stableboundary layer with a surface heat flux of(W)_o = -0.05 m K s^-1and a geostrophic wind speed of U_g = 15 m s^-1.Flow visualization and evaluation of the turbulencestatistics from this case reveal the development ofa continuously turbulent boundary layer with small-scalestructures. The stability of the boundary layercoupled with the presence of a strong capping inversionresults in the development of a dominant gravity wave atthe top of the stable boundary layer that appears to be relatedto the most unstable wave predicted by the Taylor–Goldsteinequation. As a result of the decay of turbulence aloft,a strong-low level jet forms above the boundary layer.The time dependent behaviour of the jet is compared with Blackadar'sinertial oscillation analysis.     

Footprints in Homogeneously and Heterogeneously Driven Boundary Layers Derived from a Lagrangian Stochastic Particle Model Embedded into Large-Eddy Simulation

A Lagrangian stochastic (LS) model, which is embedded into a parallelised large-eddy simulation (LES) model, is used for dispersion and footprint evaluations. For the first time an online coupling between LES and LS models is applied. The new model reproduces concentration patterns, which were obtained in prior studies, provided that subgrid-scale turbulence is included in the LS model. Comparisons with prior studies show that the model evaluates footprints successfully. Streamwise dispersion leads to footprint maxima that are situated less far upstream than previously reported. Negative flux footprints are detected in the convective boundary layer (CBL). The wide range of applicability of the model is shown by applying it under neutral and stable stratification. It is pointed out that the turning of the wind direction with height leads to a considerable dependency of source areas on height. First results of an application to a heterogeneously heated CBL are presented, which emphasize that footprints are severely affected by the inhomogeneity.     

A Theoretical Model for the Study of Convective Turbulence Decay and Comparison with Large-Eddy Simulation Data

The development of a theoretical model fora decaying convective boundary layeris considered. The model relies on thedynamical energy spectrumequation in which the buoyancy andinertial transfer terms are retained,and a closure assumptionmade for both. The parameterization for thebuoyancy term is given providing a factorizationbetween the energy source termand its temporal decay. Regarding the inertialtransfer term a hypothesis ofsuperposition is used to describe theconvective energy source and time variationof velocity correlation separately.The solution of the budget equation for theturbulent kinetic energy spectrum is possible,given the three-dimensional initial energyspectrum. This is doneutilizing a version of the Kristensen et al.(see Boundary-Layer Meteorol. 47, 149–193)model valid for non-isotropic turbulence. During thedecay the locus of the spectralpeak remains at about the sameposition as the heat flux decreases.Comparison of the theoretical modelis performed against large-eddy simulationdata for a decaying convectiveboundary layer.     

A Large-Eddy Simulation and Lagrangian Stochastic Study of Heavy Particle Dispersion in the Convective Boundary Layer

Large-eddy simulation and Lagrangian stochastic dispersion models were used to study heavy particle dispersion in the convective boundary layer (CBL). The effects of various geostrophic winds, particle diameters, and subgrid-scale (SGS) turbulence were investigated. Results showed an obvious depression in the vertical dispersion of heavy particles in the CBL and major vertical stratification in the distribution of particle concentrations, relative to the passive dispersion. Stronger geostrophic winds tended to increase the dispersion of heavy particles in the lower CBL. The SGS turbulence, particularly near the surface, markedly influenced the dispersion of heavy particles in the CBL. For reference, simulations using passive particles were also conducted; these simulation results agreed well with results from previous convective tank experiments and numerical simulations.     

Large-Eddy Simulation Of A Nocturnal Stratocumulus-Topped Marine Atmospheric Boundary Layer: An Uncertainty Analysis

A large-eddy simulation (LES) model has been used to study a nocturnalstratocumulus-topped marine atmospheric boundary layer. The main objectivesof our study have been first to investigate the statistical significance of LES-derived data products. Second, to test the sensitivity of our LES results with respect to the representation of subgrid-scale mixing and microphysical processes, and third to evaluate and to quantify the parametric uncertainty arising from the incomplete knowledge of the environmental parameters that are required to specify the initial and boundary conditions of a particular case study. Model simulations were compared with observations obtained in solid stratocumulus during the third flight of the first 'Lagrangian' experiment of the Atlantic Stratocumulus Transition Experiment (ASTEX). Based on these simulations the following conclusions could be drawn. Resolution(50 × 50 × 25 m³) and domain size (3.2 × 3.2 × 1.5 km³) of the LES calculations were adequate from a numerical point of view to represent the essential features of the stratocumulus-topped boundary layer. However, the ensemble runs performed in our study to investigate the statistical significance of LES-derived data products demonstrate that the area-time averaging procedure for the second-order moments produces only a low degree of statistical reliability in the model results. This illustratesthe necessity of having LES model results that are not only of adequate resolution but also of sufficiently large domain. The impact of different subgrid schemes was small, but the primary effects of drizzle were found to influence the boundary-layer structure in a climatologically significant way. The parametric uncertainty analysis revealed that the largest contribution to the variance of the LES-derived data products is due to theuncertainties in the cloud-top jump of total water mixing ratio and the net radiative forcing. The differences between the model and measurements for most of the simulated quantities were within the modelling uncertainties, but the calculated precipitation rate was found to differ significantly from that derived in the observations.     

Two-Point Correlation Analysis Of Neutrally Stratified Flow Within And Above A Forest From Large-Eddy Simulation

Two-point space-time correlations ofvelocities, a passive scalar and static pressure arecalculated using the resolvable flow fields computedby large-eddy simulation (LES) of neutrally stratifiedflow within and above a sparse forest. Zero-time-lagspatial auto-correlation contours in thestreamwise-vertical cross-section for longitudinal andlateral velocities and for a scalar are tilted fromthe vertical in the downstream direction, as istypical in near-wall sheared flow. On the other hand,auto-correlations of vertical velocity and of staticpressure are vertically coherent. Zero-time-lagspatial auto-correlations in the spanwise-verticalcross-section show no distinct tilt, and those forboth longitudinal and vertical velocities demonstratedistinct negative side lobes in the middle forest andabove, while longitudinal velocity in the subcrowntrunk space is laterally in-phase. Static pressureperturbations appear to be spatially coherent in thespanwise direction at all heights, especially insidethe forest. Near the forest floor, longitudinalvelocity is found to be in-phase with static pressureperturbation and to be closely linked to theinstantaneous streamwise pressure gradient, supportinga previous proposal that longitudinal velocity in thisregion is dominantly modulated by the pressurepatterns associated with the coherent sweep/ejectionevents. Near treetop height, a lack of linkage betweenthe pressure gradient and the local time derivative ofthe longitudinal velocity supports the hypothesis ofadvection dominating turbulent flow.The major phase characteristics of the two-pointcorrelations essentially remained the same from fourLES runs with different domain size and/or gridresolution. A larger LES domain yielded betteragreement with field observations in a real forest onboth the magnitudes of the correlations and thesingle-point integral time scales. A finer gridresolution in the LES led to a faster rate of decreaseof correlation with increasing separation in space ortime, as did the higher frequency fluctuations in theturbulent records from field measurements. Convectivevelocities estimated from the lagged two-pointauto-correlations of the calculated flow fields werecompared with similar calculations from wind-tunnelstudies. At the canopy top, estimates from thecorrelation analyses agree with the translationvelocity estimated from instantaneous snapshots of ascalar microfront using both LES and field data. Thistranslation velocity is somewhat higher than the localmean wind speed. Convective velocities estimated fromlagged correlations increase with height above thecanopy. It is suggested that an appropriate filteringprocedure may be necessary to reduce the effects ofsmall-scale random turbulence, as was reported in astudy over an orchard canopy. The mean longitudinalvelocity near the treetops is found to be moreappropriate than the local mean longitudinal velocityat each height to link single-point integral timescales with directly calculated spatial integralstreamwise length scales.     

Coherent Structures near the Surface in a Strongly Sheared Convective Boundary Layer Generated by Large-Eddy Simulation

Well-developed low speed and high temperature streaks in association with the alignment of convection cells are observed in a large-eddy-simulation (LES) generated strongly sheared convective boundary-layer flow, which is driven by a geostrophic wind speed of 15 m s^-1 and a surface kinematic heat flux of 0.05 K m s^-1. Vortices that drive streaky structures are identified through an eigenvalue method (lambda;₂method) near the surface. These vortices are highly elongated along the quasi-streamwise direction alternating sign of the x-component of vorticity (_x). By conditional sampling of fully developed vortices, a statistically significant coherent structure is educed. The educed vortex is elongated to the streamwise direction with the elevation angle of about 17° above the horizontal surface. However, the horizontal tilting is not clearly demonstrated in the present simulation. Fluctuation fields in the domain of the educed vortex show the existence of a low speed and high temperature streak as a direct consequence of momentum and heat transport by vortical motions. The strong ejection(upward transport of low momentum or high temperature)occurring at the higher level than that of the strong sweep (downward transport of high momentum and low temperature) can be explained from the spatial distribution of the fluctuationfields of velocity and temperature. The contribution of ejection to the Reynolds stress at z/h₁ = 0.18 is about 75%, which is slightly greater than that (70% at z/h₁ = 0.173) for the neutrally stratified atmospheric boundary layer. Ejection is also found to be dominant for the turbulent heat flux.     

The Effects of Vegetation Density on Coherent Turbulent Structures within the Canopy Sublayer: A Large-Eddy Simulation Study

Large-eddy simulation has become an important tool for the study of the atmospheric boundary layer. However, since large-eddy simulation does not simulate small scales, which do interact to some degree with large scales, and does not explicitly resolve the viscous sublayer, it is reasonable to ask if these limitations affect significantly the ability of large-eddy simulation to simulate large-scale coherent structures. This issue is investigated here through the analysis of simulated coherent structures with the proper orthogonal decomposition technique. We compare large-eddy simulation of the atmospheric boundary layer with direct numerical simulation of channel flow. Despite the differences of the two flow types it is expected that the atmospheric boundary layer should exhibit similar structures as those in the channel flow, since these large-scale coherent structures arise from the same primary instability generated by the interaction of the mean flow with the wall surface in both flows. It is shown here that several important similarities are present in the two simulations: (i) coherent structures in the spanwise-vertical plane consist of a strong ejection between a pair of counter-rotating vortices; (ii) each vortex in the pair is inclined from the wall in the spanwise direction with a tilt angle of approximately 45°; (iii) the vortex pair curves up in the streamwise direction. Overall, this comparison adds further confidence in the ability of large-eddy simulation to produce large-scale structures even when wall models are used. Truncated reconstruction of instantaneous turbulent fields is carried out, testing the ability of the proper orthogonal decomposition technique to approximate the original turbulent field with only a few of the most important eigenmodes. It is observed that the proper orthogonal decomposition reconstructs the turbulent kinetic energy more efficiently than the vorticity.     

Exploring the Effects of Microscale Structural Heterogeneity of Forest Canopies Using Large-Eddy Simulations

The Regional Atmospheric Modeling System (RAMS)-based Forest Large-Eddy Simulation (RAFLES), developed and evaluated here, is used to explore the effects of three-dimensional canopy heterogeneity, at the individual tree scale, on the statistical properties of turbulence most pertinent to mass and momentum transfer. In RAFLES, the canopy interacts with air by exerting a drag force, by restricting the open volume and apertures available for flow (i.e. finite porosity), and by acting as a heterogeneous source of heat and moisture. The first and second statistical moments of the velocity and flux profiles computed by RAFLES are compared with turbulent velocity and scalar flux measurements collected during spring and winter days. The observations were made at a meteorological tower situated within a southern hardwood canopy at the Duke Forest site, near Durham, North Carolina, U.S.A. Each of the days analyzed is characterized by distinct regimes of atmospheric stability and canopy foliage distribution conditions. RAFLES results agreed with the 30-min averaged flow statistics profiles measured at this single tower. Following this intercomparison, two case studies are numerically considered representing end-members of foliage and midday atmospheric stability conditions: one representing the winter season with strong winds above a sparse canopy and a slightly unstable boundary layer; the other representing the spring season with a dense canopy, calm conditions, and a strongly convective boundary layer. In each case, results from the control canopy, simulating the observed heterogeneous canopy structure at the Duke Forest hardwood stand, are compared with a test case that also includes heterogeneity commensurate in scale to tree-fall gaps. The effects of such tree-scale canopy heterogeneity on the flow are explored at three levels pertinent to biosphere-atmosphere exchange. The first level (zero-dimensional) considers the effects of such heterogeneity on the common representation of the canopy via length scales such as the zero-plane displacement, the aerodynamic roughness length, the surface-layer depth, and the eddy-penetration depth. The second level (one-dimensional) considers the normalized horizontally-averaged profiles of the first and second moments of the flow to assess how tree-scale heterogeneities disturb the entire planar-averaged profiles from their canonical (and well-studied planar-homogeneous) values inside the canopy and in the surface layer. The third level (three-dimensional) considers the effects of such tree-scale heterogeneities on the spatial variability of the ejection-sweep cycle and its propagation to momentum and mass fluxes. From these comparisons, it is shown that such microscale heterogeneity leads to increased spatial correlations between attributes of the ejection-sweep cycle and measures of canopy heterogeneity, resulting in correlated spatial heterogeneity in fluxes. This heterogeneity persisted up to four times the mean height of the canopy (h _c ) for some variables. Interestingly, this estimate is in agreement with the working definition of the thickness of the canopy roughness sublayer (2h _c –5h _c ).