Large-Eddy Simulation of Turbulent Flow Over a Rough Surface

A family of wall models is proposed that exhibits moresatisfactory performance than previousmodels for the large-eddy simulation (LES) of the turbulentboundary layer over a rough surface.The time and horizontally averaged statistics such asmean vertical profiles of windvelocity, Reynolds stress, turbulent intensities, turbulentkinetic energy and alsospectra are compared with wind-tunnel experimental data.The purpose of the present study is to obtain simulatedturbulent flows that are comparable with wind-tunnelmeasurements for use as the wind environment for thenumerical prediction by LES of source dispersion in theneutral atmospheric boundary layer.     

Numerical Study of the Response of an Atmospheric Surface Layer to a Spatially Nonuniform Plant Canopy

High-accuracy large-eddy simulations of neutral atmospheric surface-layer flow over a gapped plant canopy strip have been performed. Subgrid-scale (SGS) motions are parameterized by the Sagaut mixed length SGS model, with a modification to compute the SGS characteristic length self-adaptively. Shaw’s plant canopy model, taking the vertical variation of leaf area density into account, is applied to study the response of the atmospheric surface layer to the gapped dense forest strip. Differences in the region far away from the gap and in the middle of the gap are investigated, according to the instantaneous velocity magnitude, the zero-plane displacement, the potential temperature and the streamlines. The large-scale vortex structure, in the form of a roll vortex, is revealed in the region far away from the gap. The nonuniform spatial distribution of plants appears to cause the formation of the coherent structure. The roll vortex starts in the wake of the canopy, and results in strong fluctuations throughout the entire canopy region. Wind sweeps and ejections in the plant canopy are also attributed to the large vortex structure.     

Subfilter-scale Fluxes over a Surface Roughness Transition. Part I: Measured Fluxes and Energy Transfer Rates

A wind-tunnel experiment was designed and carried out to study the effect of a surface roughness transition on subfilter-scale (SFS) physics in a turbulent boundary layer. Specifically, subfilter-scale stresses are evaluated that require parameterizations and are key to improving the accuracy of large-eddy simulations of the atmospheric boundary layer. The surface transition considered in this study consists of a sharp change from a rough, wire-mesh covered surface to a smooth surface. The resulting magnitude jump in aerodynamic roughnesses, M = ln(z ₀₁/z ₀₂), where z ₀₁ and z ₀₂ are the upwind and downwind aerodynamic surface roughnesses respectively, is similar to that of past experimental studies in the atmospheric boundary layer. The two-dimensional velocity fields used in this study are measured using particle image velocimetry and are acquired at several positions downwind of the roughness transition as well as over a homogeneous smooth surface. Results show that the SFS stress, resolved strain rate and SFS transfer rate of resolved kinetic energy are dependent on the position within the boundary layer relative to the surface roughness transition. A mismatch is found in the downwind trend of the SFS stress and resolved strain rate with distance from the transition. This difference of behaviour may not be captured by some eddy-viscosity type models that parameterize the SFS stress tensor as proportional to the resolved strain rate tensor. These results can be used as a benchmark to test the ability of existing and new SFS models to capture the spatial variability SFS physics associated with surface roughness heterogeneities.     

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.     

An Analysis Of Secondary Circulations And Their Effects Caused By Small-Scale Surface Inhomogeneities Using Large-Eddy Simulation

A parallelized large-eddy simulation model has been used to investigate the effects of two-dimensional, discontinuous, small-scale surface heterogeneities on the turbulence structure of the convective boundary layer.Heterogeneities had a typical size of about the boundary-layer heightz_i. They were produced by a surface sensible heat flux pattern ofchessboard-type and of strong amplitude as typical, e.g., for the marginalice zone. The major objectives of this study were to determinethe effects of such strong amplitude heat flux variations and to specify theinfluence of different speeds and directions of the background wind.Special emphasis has been given to investigate the secondary circulations induced by the heterogeneities by means of three-dimensional phase averages.Compared with earlier studies of continuous inhomogeneities, the same sizeddiscontinuous inhomogeneities in this study show similar but stronger effects.Significant changes compared with uniform surface heating are only observedwhen the scale of the inhomogeneities is increased to z_i. Especially the vertical energy transport is much more vigorous and even the mean emperature profile shows a positive lapse rate within the whole mixed layer. However, the effects are not directly caused by the different shape of the inhomogeneities but can mainly be attributed to the large amplitude of the imposed heat flux,as it is typical for the partially ice covered sea during cold air outbreaks.The structure of the secondary flow is found to be very sensitive to the wavelength and shape of the inhomogeneities as well as to the heatflux amplitude, wind speed and wind direction. The main controlling parameter is the near-surface temperature distribution and the related horizontal pressure gradient perpendicular to the main flow direction. The secondary flow varies from a direct circulation with updraughts mainly above the centre of the heated regions to a more indirect circulation with updraughts beneath the centre and downdraughts above it. For background winds larger than 2.5 m s^–1 a roll-like circulation pattern is observed.From previous findings it has often been stated that moderate backgroundwinds of 5 m s^–1 eliminate all impacts of surface inhomogeneitiesthat could potentially be produced in realistic landscapes. However, this studyshows that the effects caused by increasing the wind speed stronglydepend on the wind direction relative to the orientation of theinhomogeneities. Secondary circulations remain strong, even for abackground wind of 7.5 m s^–1, when the wind direction is orientatedalong one of the two diagonals of the chessboard pattern. On the otherhand, the effects of inhomogeneities are considerably reduced, even undera modest background wind of 2.5 m s^–1, if the wind direction isturned by 45°. Mechanisms for the different flow regimesare discussed.     

Optimized Computation of Transfer Coefficients in Surface Layer with Different Momentum and Heat Roughness Lengths

A new method for computing the surface transfer coefficients is proposed, based on state-of-the-art empirical flux-profile relationships. The influence of the roughness length ratio is first demonstrated with the classical iterative calculation method. Then a non-iterative algorithm is developed, taking into account the difference between momentum and heat roughness lengths.The new method is validated by comparison with the reference iterative computation. The large gain-in computer processing time (CPU) time gain for the calculation of surface fluxes in Eulerian grid models is finally assessed.     

Application of a Large-Eddy Simulation Database to Optimisation of First-Order Closures for Neutral and Stably Stratified Boundary Layers

Large-eddy simulation (LES) is a well-established numerical technique, resolving the most energetic turbulent fluctuations in the planetary boundary layer. By averaging these fluctuations, high-quality profiles of mean quantities and turbulence statistics can be obtained in experiments with well-defined initial and boundary conditions. Hence, LES data can be beneficial for assessment and optimisation of turbulence closure schemes. A database of 80 LES runs (DATABASE64) for neutral and stably stratified planetary boundary layers (PBLs) is applied in this study to optimize first-order turbulence closure (FOC). Approximations for the mixing length scale and stability correction functions have been made to minimise a relative root-mean-square error over the entire database. New stability functions have correct asymptotes describing regimes of strong and weak mixing found in theoretical approaches, atmospheric observations and LES. The correct asymptotes exclude the need for a critical Richardson number in the FOC formulation. Further, we analysed the FOC quality as functions of the integral PBL stability and the vertical model resolution. We show that the FOC is never perfect because the turbulence in the upper half of the PBL is not generated by the local vertical gradients. Accordingly, the parameterised and LES-based fluxes decorrelate in the upper PBL. With this imperfection in mind, we show that there is no systematic quality deterioration of the FOC in the strongly stable PBL provided that the vertical model resolution is better than 10 levels within the PBL. In agreement with previous studies, we found that the quality improves slowly with the vertical resolution refinement, though it is generally wise not to overstretch the mesh in the lowest 500 m of the atmosphere where the observed, simulated and theoretically predicted stably stratified PBL is mostly located. The submission to a special issue of the “Boundary-Layer Meteorology” devoted to the NATO advanced research workshop “Atmospheric Boundary Layers: Modelling and Applications for Environmental Security”.     

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.     

Geometric Alignments of the Subgrid-Scale Force in the Atmospheric Boundary Layer

In recent years field experiments have been undertaken in the lower atmosphere to perform a priori tests of subgrid-scale (SGS) models for large-eddy simulations (LES). The experimental arrangements and data collected have facilitated studies of variables such as the filtered strain rate, SGS stress and dissipation, and the eddy viscosity coefficient. However, the experimental set-ups did not permit analysis of the divergence of the SGS stress (the SGS force vector), which is the term that enters directly in the LES momentum balance equations. Data from a field experiment (SGS2002) in the west desert of Utah, allows the calculation of the SGS force due to the unique 4 × 4 sonic anemometer array. The vector alignment of the SGS force is investigated under a range of atmospheric stabilities. The eddy viscosity model is likely aligned with the measured SGS force under near-neutral and unstable conditions, while its performance is unsatisfactory under stable conditions.     

A Simple Model Of The Convective Internal Boundary Layer And Its Application To Surface Heat Flux Estimates Within Polynyas

A simple model of the convective (thermal) internalboundary layer has been developed for climatologicalstudies of air-sea-ice interaction, where in situobservations are scarce and first-order estimates ofsurface heat fluxes are required. It is amixed-layer slab model, based on a steady-statesolution of the conservation of potentialtemperature equation, assuming a balance betweenadvection and turbulent heat-flux convergence. Boththe potential temperature and the surface heat fluxare allowed to vary with fetch, so the subsequentboundary-layer modification alters the fluxconvergence and thus the boundary-layer growth rate.For simplicity, microphysical and radiativeprocesses are neglected.The model is validated using several case studies.For a clear-sky cold-air outbreak over a coastalpolynya the observed boundary-layer heights,mixed-layer potential temperatures and surface heatfluxes are all well reproduced. In other cases,where clouds are present, the model still capturesmost of the observed boundary-layer modification,although there are increasing discrepancies withfetch, due to the neglected microphysical andradiative processes. The application of the model toclimatological studies of air-sea interaction withincoastal polynyas is discussed.     

A Simple And General Subgrid Model Suitable Both For Surface Layer And Free-Stream Turbulence

A new and general approach is presented to allow standard subgrid schemes to besuitable both for surface layer and free-stream turbulence. Simple modificationsto subgrid schemes are proposed and derived for any vertical stabilityconditions. They are simple to implement in models and are suitable for morecomplicated simulations such as large-eddy simulation with inhomogeneoussurface conditions or complex topography. They are also applicable to mesoscaleand large-scale models. These modifications are physically justified by recentmeasurements of spectra close to the ground. The spectral analysis presentedshows how the energy deficit of blocked turbulence for a given dissipation(`anomalous dissipation') dramatically affects the coefficients to be used insubgrid schemes. As shown for neutral and convective cases with wind shear,these changes allow us to substantially improve the prediction of profiles for themean quantities in the surface layer. Agreement with similarity laws in the unstablecase is found up to about 0.2z_i, for simulated shear, stabilityprofiles and dissipation rates of turbulent kinetic energy.     

A Homogeneous Langevin Equation Model, Part i: Simulation of Particle Trajectories in Turbulence with a Skewed Velocity Distribution

A Lagrangian stochastic model for the time evolution of the velocity of a fluid particle is presented. This model is based on a one-dimensional generalized Langevin equation, and assumes the velocity probability distribution of the turbulent fluid is skewed and spatially homogeneous. This has been shown to be an effective approach to simulating vertical dispersion in the convective boundary layer. We use a form of the Langevin equation that has a linear (in velocity) deterministic acceleration and a random acceleration that is a non-Gaussian, skewed process. For the case of homogeneous fluid velocity statistics, this 'linear-skewed' Langevin equation can be integrated explicitly, resulting in an efficient numerical simulation method. Model simulations were tested using cases for which exact, analytic statistical properties of particle velocity are known. Results of these tests show that, for homogeneous turbulence, a linear-skewed Langevin equation model can overcome the difficulties encountered in applying a Langevin equation with a skewed random acceleration. The linear-skewed Langevin equation model results are compared to results of a 'nonlinear-Gaussian' Langevin equation model, and show that the linear-skewed model is significantly more efficient.     

Spatial representativeness of single tower measurements and the imbalance problem with eddy-covariance fluxes: results of a large-eddy simulation study

A large-eddy simulation (LES) study is presented that investigates the spatial variability of temporal eddy covariance fluxes and the systematic underestimation of representative fluxes linked to them. It extends a prior numerical study by performing high resolution simulations that allow for virtual measurements down to 20 m in a convective boundary layer, so that conditions for small tower measurement sites can be analysed. It accounts for different convective regimes as the wind speed and the near-surface heat flux are varied. Moreover, it is the first LES imbalance study that extends to the stable boundary layer. It reveals shortcomings of single site measurements and the necessity of using horizontally-distributed observation networks. The imbalances in the convective case are attributed to a locally non-vanishing mean vertical advection due to turbulent organised structures (TOS). The strength of the TOS and thus the imbalance magnitude depends on height, the horizontal mean wind and the convection type. Contrary to the results of a prior study, TOS cannot generally be responsible for large energy imbalances: at low observation heights (corresponding to small towers and near-surface energy balance stations) the TOS related imbalances are generally about one order of magnitude smaller than those in field experiments. However, TOS may cause large imbalances at large towers not only in the case of cellular convection and low wind speeds, as found in the previous study, but also in the case of roll convection at large wind speeds. In the stably stratified boundary layer for all observation heights neither TOS nor significant imbalances are observed. Attempting to reduce imbalances in convective situations by applying the conventional linear detrending method increases the systematic flux underestimation. Thus, a new filter method is proposed.     

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.     

Forward-in-time and Backward-in-time Dispersion in the Convective Boundary Layer: the Concentration Footprint

The turbulence field obtained using a large-eddy simulation model is used to simulate particle dispersion in the convective boundary layer with both forward-in-time and backward-in-time modes. A Lagrangian stochastic model is used to treat subgrid-scale turbulence. Results of forward dispersion match both laboratory experiments and previous numerical studies for different release heights in the convective boundary layer. Results obtained from backward dispersion show obvious asymmetry when directly compared to results from forward dispersion. However, a direct comparison of forward and backward dispersion has no apparent physical meaning and might be misleading. Results of backward dispersion can be interpreted as three-dimensional or generalized concentration footprints, which indicate that sources in the entire boundary layer, not only sources at the surface, may influence a concentration measurement at a point. Footprints at four source heights in the convective boundary layer corresponding to four receptors are derived using forward and backward dispersion methods. The agreement among footprints derived with forward and backward methods illustrates the equivalence between both approaches. The paper shows explicitly that Lagrangian simulations can yield identical footprints using forward and backward methods in horizontally homogeneous turbulence.     

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.     

Model evaluation of the atmospheric boundary layer and mixed-layer evolution

In the present study, an attempt is made to assess the atmospheric boundary-layer (ABL) depth over an urban area, as derived from different ABL schemes employed by the mesoscale model MM5. Furthermore, the relationship of the mixing height, as depicted by the measurements, to the calculated ABL depth or other features of the ABL structure, is also examined. In particular, the diurnal evolution of ABL depth is examined over the greater Athens area, employing four different ABL schemes plus a modified version, whereby urban features are considered. Measurements for two selected days, when convective conditions prevailed and a strong sea-breeze cell developed, were used for comparison. It was found that the calculated eddy viscosity profile seems to better indicate the mixing height in both cases, where either a deep convective boundary layer develops, or a more confined internal boundary layer is formed. For the urban scheme, the incorporation of both anthropogenic and storage heat release provides promising results for urban applications.     

A Case Study on the Effects of Heterogeneous Soil Moisture on Mesoscale Boundary-Layer Structure in the Southern Great Plains, U.S.A. Part I: Simple Prognostic Model

The atmospheric boundary-layer (ABL) depth was observed by airborne lidar and balloon soundings during the Southern Great Plains 1997 field study (SGP97). This paper is Part I of a two-part case study examining the relationship of surface heterogeneity to observed ABL structure. Part I focuses on observations. During two days (12–13 July 1997) following rain, midday convective ABL depth varied by as much as 1.5 km across 400 km, even with moderate winds. Variability in ABL depth was driven primarily by the spatial variation in surface buoyancy flux as measured from short towers and aircraft within the SGP97 domain. Strong correlation was found between time-integrated buoyancy flux and airborne remotely sensed surface soil moisture for the two case-study days, but only a weak correlation was found between surface energy fluxes and vegetation greenness as measured by satellite. A simple prognostic one-dimensional ABL model was applied to test to what extent the soil moisture spatial heterogeneity explained the variation in north–south ABL depth across the SGP97 domain. The model was able to better predict mean ABL depth and variations on horizontal scales of approximately 100 km using observed soil moisture instead of constant soil moisture. Subsidence, advection, convergence/divergence and spatial variability of temperature inversion strength also contributed to ABL depth variations. In Part II, assimilation of high-resolution soil moisture into a three-dimensional mesoscale model (MM5) is discussed and shown to improve predictions of ABL structure. These results have implications for ABL models and the influence of soil moisture on mesoscale meteorology