Comparison of dynamic subgrid-scale models for simulations of neutrally buoyant shear-driven atmospheric boundary layer flows

Several non-dynamic, scale-invariant, and scale-dependent dynamic subgrid-scale (SGS) models are utilized in large-eddy simulations of shear-driven neutral atmospheric boundary layer (ABL) flows. The popular Smagorinsky closure and an alternative closure based on Kolmogorov’s scaling hypothesis are used as SGS base models. Our results show that, in the context of neutral ABL regime, the dynamic modeling approach is extremely useful, and reproduces several establised results (e.g., the surface layer similarity theory) with fidelity. The scale-dependence framework, in general, improves the near-surface statistics from the Smagorinsky model-based simulations. We also note that the local averaging-based dynamic SGS models perform significantly better than their planar averaging-based counterparts. Lastly, we find more or less consistent superiority of the Smagorinsky-based SGS models (over the corresponding Kolmogorov’s scaling hypothesis-based SGS models) for predicting the inertial range scaling of spectra.     

Large-Eddy Simulations of the Atmospheric Boundary Layer Using a New Subgrid-Scale Model – I. Slightly unstable and neutral cases

Subgrid-scale (SGS) modeling is a long-standing problem and a critical component in the large-eddy simulation (LES) of atmospheric boundary layer. A variety of SGS models with different levels of sophistication have been proposed for different needs, such as Smagorinsky's (1963) eddy viscosity model, Mason and Thomson's (1992) stochastic backscattering model, and Sullivan et al.'s (1994) near surface model. A modified Smagorinsky SGS model has been used in the LES version of Terminal Area Simulation System (TASS-LES). It has successfully simulated the buoyancy-dominated, convective atmospheric boundary layer flows, while simulations of the shear-dominated, slightly unstable, neutral, and stably stratified boundary layer flows are not so good. For the later, we used a simpler version of Sullivan et al.'s subgrid-scale model in which turbulent kinetic energy equation is not included and the model is still the first-order closure. A momentum profile matching approach is adopted in the proposed model. A series of simulations for shear-dominated, slightly unstable and neutral boundary layers are performed using different subgrid-scale models and different grid resolutions. The results are compared with those of Sullivan et al. (1994) and with empirical similarity relations for the surface layer. The simulations with the new SGS model appear to be far more satisfactory than those with the modified Smagorinsky model.     

Large-eddy simulation of turbulent flow and dispersion over a complex urban street canyon

Turbulent flow and dispersion characteristics over a complex urban street canyon are investigated by large-eddy simulation using a modified version of the Fire Dynamics Simulator. Two kinds of subgrid scale (SGS) models, the constant coefficient Smagorinsky model and the Vreman model, are assessed. Turbulent statistics, particularly turbulent stresses and wake patterns, are compared between the two SGS models for three different wind directions. We found that while the role of the SGS model is small on average, the local or instantaneous contribution to total stress near the surface or edge of the buildings is not negligible. By yielding a smaller eddy viscosity near solid surfaces, the Vreman model appears to be more appropriate for the simulation of a flow in a complex urban street canyon. Depending on wind direction, wind fields, turbulence statistics, and dispersion patterns show very different characteristics. Particularly, tall buildings near the street canyon predominantly generate turbulence, leading to homogenization of the mean flow inside the street canyon. Furthermore, the release position of pollutants sensitively determines subsequent dispersion characteristics.     

Large-Eddy Simulations of the Atmospheric Boundary Layer Using a New Subgrid-Scale Model – II. Weakly and moderately stable cases

A series of simulations under weakly to moderately stable boundary layers (SBLs) have been performed using the proposed subgrid-scale (SGS) model implemented into the Terminal Area Simulation System (TASS). The proposed SGS model incorporates some aspects of the two-part eddy viscosity SGS model of Sullivan et al. (1994) and further refinements which include the dependence of SGS mixing length on stratification, two-part separation of the SGS eddy diffusivity of heat, and more realistic empirical forms of Monin–Obukhov similarity functions. The potential temperature profiles from simulations clearly show a three-layer structure: a stable surface layer of strong gradients, a middle layer of small gradients, and an inversion layer on the top. The wind speed profiles show the formation of low level jet (LLJ). However, the sub-layer structures under moderately SBLs differ from those under weakly SBLs. Both the momentum and heat fluxes decrease almost linearly in the lower part of the SBL. The near surface values of the normalized turbulent kinetic energy (TKE/u _* ²) in all simulations are about 4 which is much less than the typical value of 5.5 under the neutral condition. The decay of turbulence first occurs in the area with large values of Richardson number (R _i<0.2). Generally, instantaneous values of the TKE and R _i at the various grid points are negatively correlated, but there is not a unique relationship between the two parameters.     

Turbulence closure models and their application in RAMS

Turbulence closures are fundamental for modelling the atmospheric diffusion in numerical codes and the resulting eddy diffusivities are key parameters in describing the transport and dispersion in the boundary layer. In this work, four turbulence closure schemes have been applied for reproducing a neutral flow over schematic complex terrain using the meteorological model RAMS. Two of the closures, a one-equation (E-l) and a two-equations (E-) model, have been implemented in RAMS in alternative to the ones originally available. In these cases, an analytical method based on the similarity theory for the atmospheric surface layer and boundary layer is adopted to calculate the empirical constants of the turbulence closures. Some examples of numerical studies performed to simulate the flow and turbulence over a 3-D hill in wind-tunnel experiment in neutral stratification are presented and discussed. An intercomparison of simulations related to different closures is considered by analysing the main features of the flow over the hill and by comparing calculated vertical profiles of turbulent kinetic energy with measured ones.     

Simulation of Ekman Boundary Layers by Large Eddy Model with Dynamic Mixed Subfilter Closure

Theoretical analysis of boundary layer turbulence has suggested a feasibility of sufficiently accurate turbulence resolving simulations at relatively coarse meshes. However, large eddy simulation (LES) codes, which employ traditional eddy-viscosity turbulence closures, fail to provide adequate turbulence statistics at coarse meshes especially within a surface layer. Manual tuning of parameters in these turbulence closures may correct low order turbulence statistics but severely harms spectra of turbulence kinetic energy (TKE). For more than decade, engineering LES codes successfully employ dynamic turbulence closures. A dynamic Smagorinsky turbulence closure (DSM) has been already tried in environmental LES. The DSM is able to provide adequate turbulence statistics at coarse meshes but it is not completely consistent with the LES equations. This paper investigates applicability of an advanced dynamic mixed turbulence closure (DMM) to simulations of Ekman boundary layers of high Reynolds number flows. The DMM differs from the DSM by explicit calculation of the Leonard term in the turbulence stress tensor. The Horizontal Array Turbulence Study (HATS) field program has revealed that the Leonard term is indeed an important component of the real turbulence stress tensor. This paper presents validation of a new LES code LESNIC. The study shows that the LES code with the DMM provides rather accurate low order turbulence statistics and the TKE spectra at very coarse meshes. These coarse LES maintain more energetic small scale fluctuations of velocity especially within the surface layer. This is critically important for success of simulations. Accurate representation of higher order turbulence statistics, however, requires essentially better LES resolution. The study also shows that LES of the Ekman boundary layer cannot be directly compared with conventionally neutral atmospheric boundary layers. The depth of the boundary layer is an important scaling parameter for turbulence statistics.     

Monitoring the transport of biomass burning emissions in South America

The atmospheric transport of biomass burning emissions in the South American and African continents is being monitored annually using a numerical simulation of air mass motions; we use a tracer transport capability developed within RAMS (Regional Atmospheric Modeling System) coupled to an emission model. Mass conservation equations are solved for carbon monoxide (CO) and particulate material (PM2.5). Source emissions of trace gases and particles associated with biomass burning activities in tropical forest, savanna and pasture have been parameterized and introduced into the model. The sources are distributed spatially and temporally and assimilated daily using the biomass burning locations detected by remote sensing. Advection effects (at grid scale) and turbulent transport (at sub-grid scale) are provided by the RAMS parameterizations. A sub-grid transport parameterization associated with moist deep and shallow convection, not explicitly resolved by the model due to its low spatial resolution, has also been introduced. Sinks associated with the process of wet and dry removal of aerosol particles and chemical transformation of gases are parameterized and introduced in the mass conservation equation. An operational system has been implemented which produces daily 48-h numerical simulations (including 24-h forecasts) of CO and PM2.5, in addition to traditional meteorological fields. The good prediction skills of the model are demonstrated by comparisons with time series of PM2.5 measured at the surface.     

Numerical simulation of stratified flow around a tall building of a complex shape

Environmental Fluid Mechanics - Results of large-eddy simulations of stably stratified atmospheric flow around an isolated, complex-shaped tall building are presented. The study focuses on the...     

A stochastic Lagrangian micromixing model for the dispersion of reactive scalars in turbulent flows: role of concentration fluctuations and improvements to the conserved scalar theory under non-homogeneous conditions

Several reaction schemes, based on the conserved scalar theory, are implemented within a stochastic Lagrangian micromixing model to simulate the dispersion of reactive scalars in turbulent flows. In particular, the formulation of the reaction-dominated limit (RDL) reaction scheme is here extended to improve the model performance under non-homogeneous conditions (NHRDL scheme). The validation of the stochastic model is obtained by comparison with the available measurements of reactive pollutant concentrations in a grid-generated turbulent flow. This test case describes the dispersion of two atmospheric reactant species (NO and O₃) and their reaction product (NO₂) in an unbounded turbulent flow. Model inter-comparisons are also assessed, by considering the results of state-of-the-art models for pollutant dispersion. The present validation shows that RDL reaction scheme provides a systematic overestimation (relative error of ca. 85% around the centreline) in computing the local reactant consumption/production rate, whereas the NHRDL scheme drastically reduces this gap (relative error lower than 5% around the centreline). In terms of NO₂ production (or reactant consumption), neglecting concentration fluctuations determines overestimations of the product mean of around 100% and a NO₂ local production of one order of magnitude higher than the reference simulation. In terms of standard deviations, the concentration fluctuations of both the passive and reactive scalars are generally of the same order of magnitude or up to 1 or 2 orders of magnitudes higher than the corresponding ensemble mean values, except for the background reactant close to the plume edges. The study highlights the importance of modelling pollutant reactions depending on the instantaneous instead of the mean concentrations of the reactants, thus quantifying the role of the turbulent fluctuations of concentration, in terms of scalar statistics (mean, standard deviation, intensity of fluctuations, skewness and kurtosis of concentration, segregation coefficient, simulated reaction rate). This stochastic particle method represents an efficient numerical technique to solve the convection–diffusion equation for reactive scalars and involves several application fields: micro-scale air quality (urban and street-canyon scales), accidental releases, impact of odours, water quality and fluid flow industrial processes (e.g. combustion).     

Analytical evaluation of mesoscale fluxes and pressure field

Terrain in natural areas is never homogeneous: there may be a variety of vegetation types and patches of vegetated and unvegetated areas which can modify the mesoscale atmospheric flow. Moreover, horizontal thermal inhomogeneities in the planetary boundary layer are a well known source of mesoscale circulation systems such as land and sea breezes, mountain-valley winds, and urban heat island circulations. Since those phenomena are not resolved in regional scale numerical models, therefore an analytic procedure able to evaluate the relative importance of mesoscale and turbulent heat fluxes associated with surface thermal heterogeneities is of crucial importance in the optic of developing a parameterization of mesoscale effects generated by these heterogeneities for use in larger scale models. In the present paper we analyze how small a horizontal variation in surface heating can be and still produce a significant mesoscale circulation, how the heat and momentum fluxes associated to mesoscale flows can penetrate deeply into the mid-troposphere, and how they modify tropospheric relevant climate parameters, such as the atmospheric static stability. In addition, we evaluate the terms of the pressure gradient force, nonlinear and linear, non-hydrostatic and hydrostatic, as function of time and space scales of the mesoscale flow. The present paper is mainly a review of analytical results, the numerical comparison and verification using RAMS is in progress.     

A numerical (LES) investigation of a shallow-water, mid-latitude, tidally-driven boundary layer

We investigate turbulent mixing in a tidally driven, mid-latitude, shallow-water basin. The study is carried out numerically at a laboratory-scale, using large-eddy simulation. We compared the results of the simulation with those of a correspondent purely oscillatory flow (Stokes boundary layer). The effect of rotation on the flow dynamics is twofold. First, rotation gives rise to a mean spanwise flow that concurs to redistribute the turbulent energy among the Reynolds stresses, in particular between the horizontal directions, thus increasing the mixing across the water column and thickening the layer where developed turbulence is observable. Second, the presence of the horizontal component of the background vorticity (latitude effect) breaks the symmetry between the two semi-cycles of the oscillation, since turbulence results suppressed/enhanced during the first/second semi-cycle. These two effects significantly modify the turbulent characteristics with respect to the purely oscillating flow, although the mechanisms that generates turbulence present similar features. The qualitative agreement between our results and some measurements carried out in two sites with characteristics similar to the case analyzed suggests that the outcomes here provided may be of general use for the analysis of mid-latitude, neutrally stratified, shallow-water basins mainly driven by semi-diurnal tidal currents.     

Instantaneous transport of a passive scalar in a turbulent separated flow

The results of large-eddy simulations of flow and transient solute transport over a backward facing step and through a 180° bend are presented. The simulations are validated successfully in terms of hydrodynamics and tracer transport with experimental velocity data and measured residence time distribution curves confirming the accuracy of the method. The hydrodynamics are characterised by flow separation and subsequent recirculation in vertical and horizontal directions and the solute dispersion process is a direct response to the significant unsteadiness and turbulence in the flow. The turbulence in the system is analysed and quantified in terms of power density spectra and covariance of velocity fluctuations. The injection of an instantaneous passive tracer and its dispersion through the system is simulated. Large-eddy simulations enable the resolution of the instantaneous flow field and it is demonstrated that the instabilities of intermittent large-scale structures play a distinguished role in the solute transport. The advection and diffusion of the scalar is governed by the severe unsteadiness of the flow and this is visualised and quantified. The analysis of the scalar mass transport budget quantifies the mechanisms controlling the turbulent mixing and reveals that the mass flux is dominated by advection.     

A closure approach with multi-scale mixing length,and its application to wind and shear stress profiles within the atmospheric boundary layer

A method to determine flow specific first-order closure for the turbulent flux of momentum in the atmospheric boundary layer (ABL) is presented. This is based on the premise that eddy viscosity is a flow rather than a fluid property, and the physically more realistic assumption that the transfer of momentum and other scalar quantities in a turbulent flow takes place by a large, but finite number of length scales, than the often used single length scale, the ‘mixing length’. The resulting eddy viscosity is flow specific and when applied to the study of the ABL, yields the vertical profiles of shear stress and mean wind velocity in good agreement with observations. The method may be extended to other types of turbulent flows, however it should be recognized that each type of flow may yield a different eddy viscosity profile. Using the derived eddy viscosity the paper presents simple analytical solutions of the ABL equations to determine observationally consistent wind speed and shear stress profiles in the ABL for a variety of practical applications including air pollution modelling.     

Physical modeling of complex airflows developing above rural terrains

The rural atmospheric boundary layer (ABL) flow was reproduced in a wind tunnel at three different simulation length scales to investigate possible effects of the simulation length scale on flow characteristics. Performance of truncated vortex generators developed for part-depth ABL wind-tunnel simulations was tested in rural terrain exposure against the full-size Counihan vortex generators. A procedure to design the ABL developing above rural type terrain has been described. The 1:395 and 1:236 simulations were created as full-depth simulations, i.e., wind characteristics throughout an entire ABL were reproduced in the wind tunnel. The 1:208 simulation was a part-depth simulation, i.e., only a lower 70% of the ABL was experimentally modelled. The projected scaled-up ABL thicknesses are 395, 354, and 416 m full-scale in the 1:395, 1:236, and 1:208 simulations, respectively. Experimental results show similar trends in all three configurations not depending on the simulation length scale factor. This clearly indicates a possibility to physically, in the wind tunnel, reproduce the same rural atmospheric airflows at different simulation length scales.     

Similarity theory and parameterization of mixing in the upper ocean

Boundary layers with small thermal and mechanical inertia are close to steady-state conditions. This underlies the Monin-Obukhov similarity theory and explains why the surface values of the fluxes can be chosen as external parameters. For fluids with large thermal inertia, such as the ocean, the thermal time scale is relatively large, and the density flux is a complex function of depth; thus, the external thermal forcing is no longer a governing parameter. However, the mechanical inertia of the upper ocean is about three orders of magnitude smaller than the thermal inertia. Consequently, the upper ocean can be considered as steady-state in the dynamic sense, to any dynamic property depends primarily on the depth, the surface momentum flux, and the vertical density structure. This property allows us to suggest an alternative formulation of the similarity theory for the stratified boundary layers through specification of a new stratification parameter which characterizes the internal density structure instead of the external density flux. The turbulent mixing coefficient is derived as dependent on the stratification parameter. The latter includes the surface stress and the integral density deficit for the entire layer above. The general form and the asymptotic behavior of the nondimensional turbulent mixing coefficient as a function of the stratification parameter are obtained using dimensional considerations. Determination of numerical parameters is based on 8 years of temperature profiles acquired at the Ocean Weather Ship (OWS) PAPA. Finally, a method for calculating the profile of the turbulent mixing coefficient is obtained. This approach reproduces the 8-year evolution of the upper ocean with the maximum rms difference of approximately 1C and the bias of 1C over the depth range 0–150 m. Additional 1-year simulation of the upper ocean at OWS CHARLEY and 9-year simulation at OWS NOVEMBER confirms reasonable applicability of this approach. The proposed simple turbulent mixing scheme reproduces the evolution of the upper ocean with accuracies similar to those obtained using much more complicated models.     

Large eddy simulation of dispersion around an isolated cubic building: evaluation of localized dynamic k_\mathrm{SGS}-equation sub-grid scale model

In the present study, the prediction accuracy of a dynamic one-equation sub-grid scale model for the large eddy simulation of dispersion around an isolated cubic building is investigated. For this purpose, the localized dynamic $k_\mathrm{SGS} $ -equation model (LDKM) is employed and the results are compared with the available experimental data and two other classic sub-grid scale models, namely, standard Smagorinsky–Lilly model (SSLM) and dynamic Smagorinsky–Lilly model (DSLM). It is shown that the three SGS models give results in good agreement with experiment. However, near the ground level of the leeward wall, dimensionless time-averaged concentration, $\langle K\rangle $ , profile is not quite similar to the experimental data. It is also demonstrated that the LDKM predicts the values of $\langle K\rangle $ on the roof, leeward and side walls more acceptably than the SSLM and DSLM. Whereas, the streamwise elongation of time-averaged structures of the plume shape is more over-estimated with the LDKM than with the other two SGS models. In terms of numerical difficulty, the LDKM is found to be stable and computationally reasonable. In addition, it does not suffer from a flow dependent constant such as the Smagorinsky coefficient employed in the SSLM model.     

Sensitivity of atmospheric dispersion simulations by HYSPLIT to the meteorological predictions from a meso-scale model

Mesoscale transport and dispersion of air pollutants from a few major point sources in the Mississippi Gulf coastal region is calculated using a coupled modeling system consisting of the atmospheric dynamical model WRF and the lagrangian particle model HYSPLIT. The sensitivity of the dispersion model results to the meteorological fields is studied by conducting an ensemble of simulations using the WRF model for the same dispersion case. Several parameterization schemes for the physical processes of boundary layer turbulence and land surface temperature/moisture prediction in WRF are used in various combinations to produce different meteorological members which are then used for dispersion simulation. The uncertainty in the simulated concentration probabilities to the meteorological model configurations and the ensemble mean are presented. The parameters used for determining the uncertainties include the wind fields, temperature, area of concentration and the levels of concentration. The results indicate that dispersion model results are influenced by the choices made in respect of the planetary boundary layer and land surface schemes in the mesoscale model to produce the meteorological forecast thereby leading to certain amount of uncertainty in the resultant concentrations. Results show that the specific choices made about the atmospheric model configuration can significantly after the simulated concentrations.     

Direct numerical simulations of boundary condition effects on the propagation of density current in wall-bounded and open channels

The propagation of density current under different boundary conditions is investigated using high resolution direct numerical simulations (DNS). A revised Kleiser and Schumann influence-matrix method is used to treat the general Robin type velocity boundary conditions and the related “tau” error corrections in the numerical simulations. Comparison of the simulation results reveals that the boundary conditions change the turbulent flow field and therefore the propagation of the front. This paper mainly focuses on the effects of boundary conditions and initial depth of the dense fluid. The differences in energy dissipation and overall front development in wall-bounded and open channels are examined. Through DNS simulations, it is evident that with the decrease of initial release depth ratio ( $D/H$ ), the effect of the top boundary becomes less important. In wall-bounded channels, there are three distinctive layers in the vertical distribution of energy dissipation corresponding to the contributions from bottom wall, interface, and top wall, respectively. In open channels, there are only two layers with the top one missing due to the shear free nature of the boundary. It is found that the energy dissipation distribution in the bottom layer is similar for cases with the same $D/H$ ratio regardless the top boundary condition. The simulation results also reveal that for low Reynolds number cases, the energy change due to concentration diffusion cannot be neglected in the energy budget. To reflect the real dynamics of density current, the dimensionless Froude number and Reynolds number should be defined using the release depth $D$ as the length scale.     

Canopy-wake dynamics and wind sheltering effects on Earth surface fluxes

The atmospheric boundary layer adjustment at the abrupt transition from a canopy (forest) to a flat surface (land or water) is investigated in a wind tunnel experiment. Detailed measurements examining the effect of canopy turbulence on flow separation, reduced surface shear stress and wake recovery are compared to data for the classical case of a solid backward-facing step. Results provide new insights into the interpretation for flux estimation by eddy-covariance and flux gradient methods and for the assessment of surface boundary conditions in turbulence models of the atmospheric boundary layer in complex landscapes and over water bodies affected by canopy wakes. The wind tunnel results indicate that the wake of a forest canopy strongly affects surface momentum flux within a distance of 35–100 times the step or canopy height, and mean turbulence quantities require distances of at least 100 times the canopy height to adjust to the new surface. The near-surface mixing length in the wake exhibits characteristic length scales of canopy flows at the canopy edge, of the flow separation in the near wake and adjusts to surface layer scaling in the far wake. Components of the momentum budget are examined individually to determine the impact of the canopy wake. The results demonstrate why a constant flux layer does not form until far downwind in the wake. An empirical model for surface shear stress distribution from a forest canopy to a clearing or lake is proposed.     

Application of a three-dimensional model for assessing effects of small clear-cuttings on radiation and soil temperature

A three-dimensional model Mixfor-3D of soil–vegetation–atmosphere transfer (SVAT) was developed and applied to estimate possible effects of tree clear-cutting on radiation and soil temperature regimes of a forest ecosystem. The Mixfor-3D model consists of several closely coupled 3D sub-models describing: forest stand structure; radiative transfer in a forest canopy; turbulent transfer of sensible heat, H₂O and CO₂ between ground surface and the atmospheric surface layer; evapotranspiration of ground surface vegetation and soil; heat and moisture transfer in soil. The model operates with the horizontal grid resolution, 2 m × 2 m; vertical resolution, 1 m and primary time step, 1 h.