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.     

Characteristics of clustered particles in skimming flows on a stepped spillway

Air–water flows at hydraulic structures are commonly observed and called white waters. The free-surface aeration is characterised by some intense exchanges of air and water leading to complex air–water structures including some clustering. The number and properties of clusters may provide some measure of the level of particle-turbulence and particle–particle interactions in the high-velocity air–water flows. Herein a re-analysis of air–water clusters was applied to a highly aerated free-surface flow data set (Chanson and Carosi, Exp Fluids 42:385–401, 2007). A two-dimensional cluster analysis was introduced combining a longitudinal clustering criterion based on near-wake effect and a side-by-side particle detection method. The results highlighted a significant number of clustered particles in the high-velocity free-surface flows. The number of bubble/droplet clusters per second and the percentage of clustered particles were significantly larger using the two-dimensional cluster analysis than those derived from earlier longitudinal detection techniques only. A number of large cluster structures were further detected. The results illustrated the complex interactions between entrained air and turbulent structures in skimming flow on a stepped spillway, and the cluster detection method may apply to other highly aerated free-surface flows.     

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.     

Subgrid-Scale Modeling of Reacting Scalar Fluxes in Large-Eddy Simulations of Atmospheric Boundary Layers

In large-eddy simulations of atmospheric boundary layer turbulence, the lumped coefficient in the eddy-diffusion subgrid-scale (SGS) model is known to depend on scale for the case of inert scalars. This scale dependence is predominant near the surface. In this paper, a scale-dependent dynamic SGS model for the turbulent transport of reacting scalars is implemented in large-eddy simulations of a neutral boundary layer. Since the model coefficient is computed dynamically from the dynamics of the resolved scales, the simulations are free from any parameter tuning. A set of chemical cases representative of various turbulent reacting flow regimes is examined. The reactants are involved in a first-order reaction and are injected in the atmospheric boundary layer with a constant and uniform surface flux. Emphasis is placed on studying the combined effects of resolution and chemical regime on the performance of the SGS model. Simulations with the scale-dependent dynamic model yield the expected trends of the coefficients as function of resolution, position in the flow and chemical regime, leading to resolution-independent turbulent reactant fluxes.     

Embedded large eddy simulation approach for pollutant dispersion around a model building in atmospheric boundary layer

In the present article, the potential of embedded large eddy simulation (ELES) approach to reliably predict pollutant dispersion around a model building in atmospheric boundary layer is assessed. The performance of ELES in comparison with large eddy simulation (LES) is evaluated in several ways. These include a number of qualitative and quantitative comparisons of time-averaged and instantaneous results with wind tunnel measurements supplemented by statistical data analyses using scatter plots and standard evaluation metrics. Results obtained by both LES and ELES approaches show very good agreement with the experiment. However, addition of turbulence to mean flow at Reynolds averaged Navier–Stokes (RANS)–LES interface in ELES approach not only increases the turbulence intensity, it also results in larger values of turbulent kinetic energy (TKE) as well as a shorter reattachment length in the wake region. Accordingly, higher levels of TKE predicted by ELES increase the local intensity of concentration leading to shorter plume shapes as compared with LES. In general, ELES shows better agreement with experiment on the surfaces of model building and also in the downstream wake region. In terms of computational costs, the CPU time required to obtain statistical values in ELES is about 49 % lower than that of LES and the number of iterations per time step is also reduced by 55 % as compared with LES.     

Faster phototransformation of the formate (terrestrial) versus methanesulphonate (marine) markers of airborne particles: implications for modelling climate change

Atmospheric particulate matter is altering climate. For instance marine biogenic particles are cooling climate. Organic markers are major tools to elucidate the sources of atmospheric particulate matter. Formate is commonly used as a marker of continental aerosols, whereas methanesulphonate is used as tracer of biogenic marine aerosols. However, transformation processes during aerosol transport may modify their relative concentrations and, in turn, introduce a bias in the assessment of particle sources. Actually very little is known about the transformation of formate and methanesulphonate in aerosols. Therefore, we irradiated formate and methanesulphonate in the presence of nitrate and haematite. Nitrate and haematite are aerosol photosensitisers, producing reactive species that degrade organic compounds. The time evolution of formate and methanesulphonate was monitored by ion chromatography. Our results show that formate is transformed from 1.6 to 4.1 times faster than methanesulphonate. This trend is partly due to higher reactivity with the hydroxyl radical and partly due to additional reaction with other transients such as nitrogen dioxide. Such results strongly suggest faster formate transformation during particle transport. Therefore, when formate and methanesulphonate are used as particle tracers, an overestimation of marine biogenic versus continental particle sources is expected. This bias has major implications for climate prediction models, because marine biogenic particles have a cooling effect on climate.     

Experimental assessment of characteristic turbulent scales in two-phase flow of hydraulic jump: from bottom to free surface

A hydraulic jump is a turbulent shear flow with a free-surface roller. The turbulent flow pattern is characterised by the development of instantaneous three-dimensional turbulent structures throughout the air–water column up to the free surface. The length and time scales of the turbulent structures are key information to describe the turbulent processes, which is of significant importance for the improvement of numerical models and physical measurement techniques. However, few physical data are available so far due to the complexity of the measurement. This paper presents an investigation of a series of characteristic turbulent scales for hydraulic jumps, covering the length and time scales of turbulent flow structures in bubbly flow, on free surface and at the impingement point. The bubbly-flow turbulent scales are obtained for Fr = 7.5 with 3.4 × 10⁴ < Re < 1.4 × 10⁵ in both longitudinal and transverse directions, and are compared with the free-surface scales. The results highlight three-dimensional flow patterns with anisotropic turbulence field. The turbulent structures are observed with different length and time scales respectively in the shear flow region and free-surface recirculation region. The bubbly structures next to the roller surface and the free-surface fluctuation structures show comparable length and time scales, both larger than the scales of vortical structures in the shear flow and smaller than the scales of impingement perimeter at the jump toe. A decomposition of physical signals indicates that the large turbulent scales are related to the unsteady motion of the flow in the upper part of the roller, while the high-frequency velocity turbulence dominates in the lower part of the roller. Scale effects cannot be ignored for Reynolds number smaller than 4 × 10⁴, mainly linked to the formation of large eddies in the shear layer. The present study provides a comprehensive assessment of turbulent scales in hydraulic jump, including the analyses of first data set of longitudinal bubbly-flow integral scales and transverse jump toe perimeter integral scales.     

A modification of the classical Ekman model accounting for the Stokes drift and stratification effects

A modification of the classical Ekman model of oceanic wind-driven currents including the Stokes drift and stratification effects is discussed. The modification is formulated as an application of turbulence mechanics accounting for the curvature effect of velocity fluctuation streamlines. It is shown that similar to the Stokes drift effect, the presence of a density jump layer (pycnocline) decreases the veering of the flow velocity vector at the surface from the direction of the wind stress. It is shown also that in the pycnocline the decrease of the norm of the velocity vector as well as its rotation with depth is smaller than in the regions adjacent to the pycnocline. If the Stokes drift and stratification effects are neglected, the model reduces to the classical Ekman solution with the coefficient of the turbulent shear viscosity replaced by an effective viscosity coefficient. The vertical distributions of velocity predicted by the modified model are compared with the velocity data measured in the Drake Passage and within the Long-Term Upper Ocean Study (LOTUS) in the North Atlantic.     

Small-scale entrainment in inclined gravity currents

We investigate the effect of buoyancy on the small-scale aspects of turbulent entrainment by performing direct numerical simulation of a gravity current and a wall jet. In both flows, we detect the turbulent/nonturbulent interface separating turbulent from irrotational ambient flow regions using a range of enstrophy iso-levels spanning many orders of magnitude. Conform to expectation, the relative enstrophy isosurface velocity \(v_n\) in the viscous superlayer scales with the Kolmogorov velocity for both flow cases. We connect the integral entrainment coefficient E to the small-scale entrainment and observe excellent agreement between the two estimates throughout the viscous superlayer. The contribution of baroclinic torque to \(v_n\) is negligible, and we show that the primary reason for reduced entrainment in the gravity current as compared to the wall-jet are 1) the reduction of \(v_n\) relative to the integral velocity scale \(u_T\); and 2) the reduction in the surface area of the isosurfaces.     

Hydrodynamic characteristics and related mass-transfer properties in open-channel flows with rectangular embayment zone

Consecutive groynes and embayments form dead water zones, where sedimentation and high concentrations of pollutants are often observed. It is thus very important to understand the mass and momentum exchange between the main channel and side cavities in rivers and hydraulic engineering structures. The spanwise gradient of the streamwise velocity near the junction produces small-scale turbulent vortices because of shear instability. Furthermore, large-scale horizontal circulation is also generated in the cavity zone. These coherent turbulent structures play a significant role in mass and sediment transfer at the boundary between the mainstream and embayment. However, the relation between turbulence and mass transfer is poorly understood. In this study, we performed particle image velocity and laser-induced fluorescence experiments using a laboratory flume, laser light sheets and a high-speed CMOS camera. We examined the exchange properties of a dye as a function of bed configuration and sedimentation effect. Both primary and secondary gyres were observed in the flat bed and downward-sloping bed, whereas the primary gyre was prevalent in the upward-sloping bed. Moreover, the horizontal circulation strongly affected the mass-transfer properties between the mainstream and side cavity.     

大气气溶胶研究新动向

本文概述了近年大气气溶胶研究的四个主要方面:大气气溶胶的表征、大气化学过程与气候变化、健康效应等的概况.阐明了当今大气气溶胶研究的趋向,主要着重于PM_(10)和PM_(2.5)细颗粒(可吸入颗粒物或二次颗粒物)的物理化学特性与环境行为、生态效应,更重视气溶胶的非均相化学反应过程;其研究范围,从平流层向对流层发展,并密切结合气候变化、健康影响等有关的一些实际问题,进行深层次的综合性研究.大气气溶胶化学已成为当今大气化学研究最前沿的领域.     

Roof-level large- and small-scale coherent structures in a street canyon flow

Environmental Fluid Mechanics - The characteristics of large- and small-scale turbulent motions at roof-level in a street canyon flow were experimentally investigated along with their...     

烟炱气溶胶的光学辐射特性研究

烟炱气溶胶是大气气溶胶的重要组成部分,它能直接或间接地影响地-气辐射平衡,进而对全球气候产生重大影响.目前,关于烟炱气溶胶光学性质的研究已经成为气溶胶科学中的前沿和热点课题,对其进行系统的实验室研究可以为更好地分析和认识外场观测结果提供基础,并为数值模式模拟研究提供准确的基本参数.本文简要介绍了烟炱气溶胶的来源和分类,...     

Turbulence and aeration in hydraulic jumps: free-surface fluctuation and integral turbulent scale measurements

In an open channel, a change from a supercritical to subcritical flow is a strong dissipative process called a hydraulic jump. Herein some new measurements of free-surface fluctuations of the impingement perimeter and integral turbulent time and length scales in the roller are presented with a focus on turbulence in hydraulic jumps with a marked roller. The observations highlighted the fluctuating nature of the impingement perimeter in terms of both longitudinal and transverse locations. The results showed further the close link between the production and detachment of large eddies in jump shear layer, and the longitudinal fluctuations of the jump toe. They highlighted the importance of the impingement perimeter as the origin of the developing shear layer and a source of vorticity. The air–water flow measurements emphasised the intense flow aeration. The turbulent velocity distributions presented a shape similar to a wall jet solution with a marked shear layer downstream of the impingement point. The integral turbulent length scale distributions exhibited a monotonic increase with increasing vertical elevation within 0.2 < L_z/d₁ < 0.8 in the shear layer, where L_z is the integral turbulent length scale and d₁ the inflow depth, while the integral turbulent time scales were about two orders of magnitude smaller than the period of impingement position longitudinal oscillations.     

A Fine-Scale <Emphasis Type="Italic">k</Emphasis>−ε Model for Atmospheric Flow over Heterogeneous Landscapes

A multi-purpose model for small-scale atmospheric flows over heterogeneous landscapes is being developed. The aim of this research is to build a tool able to predict the dynamical (wind, turbulence) and diffusive (gases, particles) fields over landscapes characterised by heterogeneous plant cover. In its present stage of development the model is based on the numerical integration of neutral atmospheric flow equations, using an energy-dissipation closure scheme and over a domain that may include vegetation layers. Three validation cases of the model are presented: (i) response of the airflow to a change in surface roughness; (ii) airflow within and above a horizontally homogeneous plant canopy; (iii) airflow over two complex forest-to-clearing and clearing-to-forest transitions. All simulations provide results in good agreement with the experimental data, except for turbulent kinetic energy just after a clearing-to-forest transition. This result is not surprising for a statistical k– model in a flow region characterised by strong distorsion and intermittent turbulence. However the overall good performance of the model is promising for environmental research at fine scales over heterogeneous landscapes.     

Analytical expression for predicting the reduced settling velocity of small particles in turbulence

The influence of turbulence on the settling velocity of small particles remains an inconclusive research subject. Both enhanced and retarded particle settling compared to quiescent settling have been reported in previous literature, and several theories have been proposed. Among the mechanisms that account for reduced settling velocities in turbulence, the loitering effect is an important one that generally exists in various conditions. This study focuses on an analytical prediction of reduced particle settling velocities due to the loitering effect. By considering the velocity autocorrelation function as a step function analogous to the free path theory, and using a modified integral time scale to capture the essence of the loitering effect, an analytical expression is derived for predicting the reduction in particle settling velocities in turbulent flows. Calculation results of the expression are then examined by comparing with the results of a random walk model, direct numerical simulations that conditionally captured the reduction of particle settling velocities due to the loitering effect, and several representative experiments. Major possible influencing factors on the prediction of the analytical expression and applicable conditions of the expression are then further discussed. The proposed analytical expression is shown to be suitable for predicting the reduced settling velocities of small particles with relatively weak inertia in turbulent flows and could provide a reasonable explanation for reported cases in which small to moderate reductions in particle settling velocities were observed.

     

Observations on the small-scale distribution of estuarine phytoplankton

Analyses of phytoplankton collections from Narragansett Bay, Rhode Island and the Damariscotta River, Maine, USA showed that statistically significant differences in population density can usually be demonstrated within any series of collections when the collecting interval is more than 10 cm. Individual species in any series are in general over-dispersed. Phytoplankton patches ranging in size from <1 m to 12 m were tentatively identified. No correlation was found between temperature or salinity and the small-scale distribution of individual species. Possible causes of the inhomogeneities observed in phytoplankton distributions are considered.Contribution No. 17 from the Ira C. Darling Center for Research, Teaching, and Service, University of Maine.     

Wind-induced pressure at a tunnel portal

In order to properly size the mechanical ventilation system of a tunnel, it is essential to estimate the wind-driven pressure difference that might rise between its two portals. In this respect, we explore here the pressure distribution over a tunnel portal under the influence of an incident atmospheric boundary layer and, in particular, its dependency on wind direction and on tunnel geometry. Reduced scale models of generic configurations of a tunnel portal are studied in an atmospheric wind tunnel. Pressure distributions over the front section of different open cavities are measured with surface taps, which allows us to infer the influence of the tunnel aspect ratio and wind direction on a pressure coefficient \(C_{P}\), defined as a spatially and time averaged non-dimensional pressure. Experiments reveal that the magnitude of the coefficient \(C_{P}\), as a function of the wind direction, is significantly influenced by the portal height-to-width ratio and almost insensitive to its length. The experimental data set is completed by hot-wire anemometry measurements providing vertical distribution of velocity statistics. The same configurations are simulated by numerically solving the Reynolds-averaged Navier–Stokes equations, adopting the standard \(k - \varepsilon\) turbulence model. Despite some discrepancies between numerical and experimental estimates of some flow parameters (namely the turbulent kinetic energy field), the numerical estimates of the pressure coefficients \(C_{P}\) show very good agreement with experimental data. The latter is also compared to the predictions of an analytical model, based on the estimate of a spatially averaged velocity within an infinitely long street canyon. The results of the model, which takes into account varying canyon aspect ratios, are in reasonable agreement with experimental data for all cases studied. Notably, its predictions are significantly better than those provided by the simple analytical relations usually adopted as a reference in tunnel ventilation studies.     

A CFD study of the transport and fate of airborne droplets in a ventilated office: The role of droplet−droplet interactions

• Coulomb and Lennard−Jones forces were considered for droplet interactions. • The net droplet interactions were repulsive. • Repulsive droplet interactions increased the transport of droplets. • Repulsive droplet interactions significantly modified the fate of droplets. Previous studies reported that specially designed ventilation systems provide good air quality and safe environment by removing airborne droplets that contain viruses expelled by infected people. These water droplets can be stable in the environment and remain suspended in air for prolonged periods. Encounters between droplets may occur and droplet interactions should be considered. However, the previous studies focused on other physical phenomena (air flow, drag force, evaporation) for droplet transport and neglected droplet interactions. In this work, we used computational fluid dynamics (CFD) to simulate the transport and fate of airborne droplets expelled by an asymptomatic person and considered droplet interactions. Droplet drag with turbulence for prediction of transport and fate of droplets indicated that the turbulence increased the transport of 1 μm droplets, whereas it decreased the transport of 50 μm droplets. In contrast to only considering drag and turbulence, consideration of droplet interactions tended to increase both the transport and fate. Although the length scale of the office is much larger than the droplet sizes, the droplet interactions, which occurred at the initial stages of release when droplet separation distances were shorter, had a significant effect in droplet fate by considerably manipulating the final locations on surfaces where droplets adhered. Therefore, it is proposed that when an exact prediction of transport and fate is required, especially for high droplet concentrations, the effects of droplet interactions should not be ignored.     

Stability and Mixing of a Vertical Axisymmetric Buoyant Jet in Shallow Water

The stability, mixing and effect of downstream control on axisymmetric turbulent buoyant jets discharging vertically into shallow stagnant water is studied using 3D Reynolds-averaged Navier–Stokes equations (RANS) combined with a buoyancy-extended k –ε model. The steady axisymmetric turbulent flow, temperature (or tracer concentration) and turbulence fields are computed using the finite volume method on a high resolution grid. The numerical predictions demonstrate two generic flow patterns for different turbulent heated jet discharges and environmental parameters (i) a stable buoyant discharge with the mixed fluid leaving the vertical jet region in a surface warm water layer; and (ii) an unstable buoyant discharge with flow recirculation and re-entrainment of heated water. A stratified counterflow region always appears in the far-field for both stable and unstable buoyant discharges. Provided that the domain radius L exceeds about 6H, the near field interaction and hence discharge stability is governed chiefly by the jet momentum length scale to depth ratio l_M/H, regardless of downstream control. The near field jet stability criterion is determined to be l_M/H = 3.5. A radial internal hydraulic jump always exists beyond the surface impingement region, with a 3- to 6-fold increase in dilution across the jump compared with vertical buoyant jet mixing. The predicted stability category, velocity and temperature/concentration fields are well-supported by experiments of all previous investigators.