Numerical and wind tunnel modeling on the windbreak effectiveness to control the aeolian erosion of conical stockpiles

There are many industrial sites where open aggregate conical piles exist to store granular materials, like coal, industrial residuals, or other minerals. Usually these storage piles are placed in open areas, making them susceptible to wind erosion, which can create health, environmental, and/or economical concerns. It is common to minimize the dust emission through the placement of windbreaks in the vicinity of the storage piles, which reduces the wind speed in the vicinity of the pile’s surface. In this work, some experimental results from a wind tunnel study on the erosion of a conical sand pile, exposed or protected by a fence with porosities of 0, 70, and 83%, are shown. For the sheltered cases, the windbreak was placed at several distances from the pile’s leading edge, ranging from H to 4H, where H is the initial height of the non-eroded pile. The evolution of the shape of the sand pile is shown, at different instances in time, and the pile deformation quantified, using a novel experimental setup developed for wind erosion studies. This information might be regarded as a useful dataset for the benchmark of computational models aiming to produce the transient simulation of the aeolian erosion of stockpiles. The CFD results are comprised of the modeling of several experimental scenarios. The computational results for the surface wind velocity show a good correlation with the initial deformation of the pile. Based on the results, the isocontours of (u _s /u _r ) presented might be regarded as a good basis for the estimation of the pile shear velocity.     

Effect of the topography of an industrial site on dust emissions from open storage yards

The erosion, transport and subsequent deposition of particles by wind has been identified as an important environmental problem worldwide. Investigations assessing fugitive dust emissions represent a relevant subject for various applications such as prevention of desert expansion, farmland erosion, quantification of dust emissions due to wind erosion from mining zones or open storage yards of granular materials. The present study, within this latter framework, aims to improve the methods of quantification of fugitive emissions from granular materials stockpiles on industrial sites. The methodology mostly used by the industrial operators to quantify dust emissions from their plants is based on the wind flow analysis and application of emission factors into such pattern. In this paper, three-dimensional numerical simulations were done to simulate wind flow over a real configuration of a power plant to highlight the influence of the surroundings buildings on the real exposure of the granular materials stockpiles. This study emphasizes the necessity to take into consideration the presence of surrounding buildings for the calculation of the particle emissions rate on stockpiles in order to obtain a more accurate and relevant evaluation of fugitive dust emissions from open storage systems on industrial sites.     

Experimental surface flow visualization and numerical investigation of flow structure around an oblong stockpile

Emission factors are largely used to quantify particle emissions from industrial open storage piles. These factors are based on the knowledge of velocity distribution and flow patterns over the stockpile surface which still requires further research. The aim of the present work is to investigate the airflow characteristics over a single typical oblong pile and in its near-ground surroundings for various wind flow directions. Wind tunnel experiments using an oil-film surface coating technique were carried out for near-wall flow visualization. Numerical simulation results, favorably compared to PIV measurements, were used to allow comparison analysis of flow features. For the stockpile oriented 90° to the wind main direction, typical topology of flow around wall-mounted obstacles were observed, notably a wake zone downstream the pile including two main counter-rotating vortices. Further analysis of numerical wall shear stress distribution and streamlines indicates that two complex three-dimensional vortical flow structures develop downstream the pile. For other incoming wind flow directions (30 and 60°), the flow characteristics over the storage pile greatly differ as a single helical main vortex develops from the pile’s crest. Corresponding high values of wall shear stress are noticed downstream the storage pile. For each configuration studied, downwash and upwash zones are induced by the vortical structures developed. This near-wall flow topology combined with areas of high friction levels may be linked to potential dust emission from the ground surface surrounding industrial stockpiles.     

Aeolian erosion of storage piles yards: contribution of the surrounding areas

Dust emissions from stockpiles surfaces are often estimated applying mathematical models such as the widely used model proposed by the USEPA. It employs specific emission factors, which are based on the fluid flow patterns over the near surface. But, some of the emitted dust particles settle downstream the pile and can usually be re-emitted which creates a secondary source. The emission from the ground surface around a pile is actually not accounted for by the USEPA model but the method, based on the wind exposure and a reconstruction from different sources defined by the same wind exposure, is relevant. This work aims to quantify the contribution of dust re-emission from the areas surrounding the piles in the total emission of an open storage yard. Three angles of incidence of the incoming wind flow are investigated ( $30^{\circ }, 60^{\circ }$ and $90^{\circ }$ ). Results of friction velocity from numerical modelling of fluid dynamics were used in the USEPA model to determine dust emission. It was found that as the wind velocity increases, the contribution of particles re-emission from the ground area around the pile in the total emission also increases. The dust emission from the pile surface is higher for piles oriented $30^{\circ }$ to the wind direction. On the other hand, considering the ground area around the pile, the $60^{\circ }$ configuration is responsible for higher emission rates (up to 67 %). The global emissions assumed a minimum value for the piles oriented perpendicular to the wind direction for all wind velocity investigated.     

Development and testing of a micro wind tunnel for on-site wind erosion simulations

Wind erosion processes affect soil surfaces across all land uses worldwide. Understanding the spatial and temporal scales of wind erosion is a challenging undertaking because these processes are diverse and highly variable. Wind tunnels provide a useful tool as they can be used to simulate erosion at small spatial scales. Portable wind tunnels are particularly valued because erosion can be simulated on undisturbed soil surfaces in the field. There has been a long history of use of large portable wind tunnels, with consensus that these wind erosion simulation tools can meet real world aerodynamic criteria. However, one consequence of striving to meet aerodynamic reality is that the size of the tunnels has increased, making them logistically difficult to work with in the field and resulting in a tendency to homogenise naturally complex soil surfaces. This homogenisation is at odds with an increasing awareness of the importance that small scale processes have in wind erosion. To address these logistical and surface homogenisation issues we present here the development and testing of a micro wind tunnel (MWT) designed to simulate wind erosion processes at high spatial resolution. The MWT is a duct-type design—0.05 m tall 0.1 m wide and with a 1.0 m working section. The tunnel uses a centrifugal motor to suck air through a flow‐conditioning section, over the working section and then through a sediment collection trap. Simulated wind velocities range from 5 to 18 m s^?1, with high reproducibility. Wind speeds are laterally uniform and values of u _* at the tunnel bed (calculated by measuring the pressure gradients within the MWT) are comparable with those of larger tunnels in which logarithmic profiles can be developed. Saltation sediment can be added. The tunnel can be deployed by a single person and operated on slopes ranging from 0 to 10°. Evidence is presented here that the MWT provides new and useful understanding of the erodibility of rangelands, claypans and ore stockpiles.     

Design and Testing of a Turbulence Probe for Harsh Flows

The force of wind on the ground created by turbulent eddies is commonly used to describe the horizontal flux of material during wind erosion. Here we present the Murdoch Turbulence Probe, an instrument for use in both clean and eroding flows which uses pressure differences to measure the three components of wind velocity. Correlation techniques calculate the forces near the ground and turbulence statistics in nearly real time, including turbulent velocity fluctuations from less than 0.1 Hz to 200 Hz, mean flow velocities, Reynolds stresses as well as the integral length and time scales. In the portable wind-tunnel used by Agriculture Western Australia, turbulence statistics were recorded over stable surfaces and in blowing sand from the initiation of erosion up to the time the sand supply was exhausted. Estimates of the friction velocity derived from the turbulence probe were compared with estimates obtained from the wind speed profile measured with a rake of pitot and static tubes. The Murdoch Turbulence Probe appears to work well in sandblasting conditions. Relative turbulence intensities ranged from 0.11 to 0.2 and are in close agreement with values in the literature. The ratio of the turbulence to the friction velocity (3 to 3.2) is at the high end of the reported range. The Reynolds stress measurements agree closely with predictions of the threshold friction velocities of the sand and estimates from the wind speed profile with a von Kármán constant of about 0.3, lower than the commonly accepted value of 0.4. We suggest that the wind-tunnel profile represents the `outer layer' of the boundary-layer that may best be described by a `Wake Law' or `Defect Law'. At about 54 mm above the surface, the friction velocity decreases from 0.64 m/s to 0.39 m/s and the mean velocity increases from 9.6 m/s to 11.6 m/s as the supply of sand is depleted. In addition to the friction velocity, other scales may be necessary to characterise the overriding effect of the wind and in extending wind-tunnel results to the field.     

Mixing at a sediment concentration interface in turbulent open channel flow

In this work we address the role of turbulence on mixing of clear layer of fluid with sediment-laden layer of fluid at a sediment concentration interface. This process can be conceived as the entrainment of sediment-free fluid into the sediment-laden layer, or alternatively, as the transport of sediment into the top sediment-free flow. This process is governed by four parameters—Reynolds number of the flow \(Re_\tau\), non-dimensional settling velocity of the sediment (proxy for sediment size) \(\tilde{V}\), Richardson number \(Ri_\tau\) and Schmidt number Sc. For this work we have performed direct numerical simulations for fixed Reynolds and Schmidt numbers while varying the values of Richardson number and particle settling velocity. In the simple model considered here, the flow’s momentum and turbulence pre-exists over the entire layer of fluid, while the sediment is initially confined to a layer close to the bed. Mixing of sediment-free fluid with the sediment-laden layer is associated primarily with upward transport of sediment and buoyancy. There is no simultaneous upward transport of fluid momentum and turbulence into the sediment-free fluid layer, which is already in motion and turbulent. The analysis performed shows that the ability of the flow to transport a given sediment size decreases with the distance from the bottom, and thus only fine enough sediment particles are transported across the sediment concentration interface. For these cases, the concentration profiles evolve to a final steady state in good agreement with the well-known Rouse profile. The approach towards the Rouse profile happens through a transient self-similar state. This behavior of the flow is not seen for larger particles. Detailed analysis of the three dimensional structure of the sediment concentration interface shows the mechanisms by which sediment particles are lifted up by tongues of sediment-laden fluid with positive correlation between vertical velocity and sediment concentration. Finally, the mixing ability of the flow is addressed by monitoring the time evolution of the center of mass of the sediment-laden layer and the vertical location of the sediment-free/sediment-laden interface.     

Analytical model for vertical velocity profiles in flows with submerged shrub-like vegetation

An analytical solution for the vertical profiles of the horizontal velocity of channel flow with submerged shrub-like vegetation is investigated in this paper. At first, a shape function is proposed to fit the diameter change of different types of shrub-like vegetation. Using the momentum theorem and the mixing-length turbulence model, an analytical solution for the vertical profile of the horizontal velocity within the vegetation is obtained. The velocity distribution of the whole column is determined in tandem with the logarithmic velocity profile above the vegetation. The solution is compared with experimental data in excellent agreement. The results show that the flow above the vegetation has a logarithmic velocity profile while the flow within the vegetation is impacted greatly by the shape and density of vegetation. The flows within shrub-like vegetations are non-uniform and vary inversely with the shrub diameter.     

Double-average mean flow and local turbulence intensity profiles from PIV measurements for an open channel flow with rigid vegetation

This contribution presents particle image velocimetry measurements for an open channel stationary uniform and fully developed flow of water over a horizontal flat bed of uniform glass beads in presence of a staggered array of vertical cylindrical stems. The main objective was to explore and quantify the influence of the stems-to-flow relative submergence, h _v /h, over the mean flow and local turbulence intensities. A comparison with measurements for the non-vegetated flow over the same granular bed is presented. Results indicate a remarkable influence of h _v /h over the whole flow field. The time-average mean flow presents a strong spatial variation in the layer of the flow occupied by the stems. The local velocity fluctuations are strongly affected by the presence of the stems, with regions in between the stems where they reach peaks that are several times larger than those encountered in the flow in absence of vegetation. The turbulence intensity profiles are noticeably different when compared to those measured in the non-vegetated flow conditions. From previous works it was possible to derive an equation for the mean velocity, U _v , of the flow through the vegetated layer of height h _v . The prediction of this equation is in good agreement with the uniform value for the double-average longitudinal velocity profile in this layer. A final brief discussion about the possible impact of these vegetated-flow features on the sediment transport is presented.     

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.     

Characteristics of flow structure of free-surface flow in a partly obstructed open channel with vegetation patch

Free-surface flows over patchy vegetation are common in aquatic environments. In this study, the hydrodynamics of free-surface flow in a rectangular channel with a bed of rigid vegetation-like cylinders occupying half of the channel bed was investigated and interpreted by means of laboratory experiments and numerical simulations. The channel configurations have low width-to-depth aspect ratio (1.235 and 2.153). Experimental results show that the adjustment length for the flow to be fully developed through the vegetation patch in the present study is shorter than observed for large-aspect-ratio channels in other studies. Outside the lateral edge of the vegetation patch, negative velocity gradient (\(\partial \overline{u}/\partial z < 0\)) and a local velocity maximum are observed in the vertical profile of the longitudinal velocity in the near-bed region, corresponding to the negative Reynolds stress (\(- \overline{{u^{\prime}w^{\prime}}} < 0\)) at the same location. Assuming coherent vortices to be the dominant factor influencing the mean flow field, an improved Spalart–Allmaras turbulence model is developed. The model improvement is based on an enhanced turbulence length scale accounting for coherent vortices due to the effect of the porous vegetation canopy and channel bed. This particular flow characteristic is more profound in the case of high vegetation density due to the stronger momentum exchange of horizontal coherent vortices. Numerical simulations confirmed the local maximum velocity and negative gradient in the velocity profile due to the presence of vegetation and bed friction. This in turn supports the physical interpretation of the flow processes in the partly obstructed channel with vegetation patch. In addition, the vertical profile of the longitudinal velocity can also be explained by the vertical behavior of the horizontal coherent vortices based on a theoretical argument.     

Mass-based depth and velocity scales for gravity currents and related flows

Gravity driven flows on inclines can be caused by cold, saline or turbid inflows into water bodies. Another example are cold downslope winds, which are caused by cooling of the atmosphere at the lower boundary. In a well-known contribution, Ellison and Turner (ET) investigated such flows by making use of earlier work on free shear flows by Morton, Taylor and Turner (MTT). Their entrainment relation is compared here with a spread relation based on a diffusion model for jets by Prandtl. This diffusion approach is suitable for forced plumes on an incline, but only when the channel topography is uniform, and the flow remains supercritical. A second aspect considered here is that the structure of ET’s entrainment relation, and their shallow water equations, agrees with the one for open channel flows, but their depth and velocity scales are those for free shear flows, and derived from the velocity field. Conversely, the depth of an open channel flow is the vertical extent of the excess mass of the liquid phase, and the average velocity is the (known) discharge divided by the depth. As an alternative to ET’s parameterization, two sets of flow scales similar to those of open channel flows are outlined for gravity currents in unstratified environments. The common feature of the two sets is that the velocity scale is derived by dividing the buoyancy flux by the excess pressure at the bottom. The difference between them is the way the volume flux is accounted for, which—unlike in open channel flows—generally increases in the streamwise direction. The relations between the three sets of scales are established here for gravity currents by allowing for a constant co-flow in the upper layer. The actual ratios of the three width, velocity, and buoyancy scales are evaluated from available experimental data on gravity currents, and from field data on katabatic winds. A corresponding study for free shear flows is referred to. Finally, a comparison of mass-based scales with a number of other flow scales is carried out for available data on a two-layer flow over an obstacle. Mass-based flow scales can also be used for other types of flows, such as self-aerated flows on spillways, water jets in air, or bubble plumes.     

Dust emission calculations in open storage piles protected by means of barriers, CFD and experimental tests

The importance of open mineral storage piles in bulk solids port terminals have increased considerably in recent years in Europe and USA (in Spain, great extensions of transoceanic ports are being made) to address the increasing demand of raw material importation due to local mining operation closure. These storage piles are affected by the climatologically processes and the dust emission to the atmosphere, causing repercussions to health and environmental, which is intimately related to the air velocity, according to USEPA studies, maximum ratio of dust emission emitted from CFD and software Ansys CFX 10.0 is determined. In addition, an emission variability study based on the pile distance of different height solid barriers is developed, reducing the emission to 66%. The studies have been contrasted by published investigations and industrial measurement.     

Dressler’s theory for curved topography flows: iterative derivation,transcritical flow solutions and higher-order wave-type equations

The Dressler equations are a system of two non-linear partial differential equations for shallow fluid flows over curved topography. The theory originated from an asymptotic stretching method formulating the equations of motion in terrain-fitted curvilinear coordinates. Apparently, these equations failed to produce a transcritical flow profile changing from sub- to supercritical flow conditions. Further, wave-like motions over a flat bottom are excluded because the bed-normal velocity component is not accounted for. However, the theory was found relevant for several environmental flow problems including density currents over mountains and valleys, seepage flow in hillslope hydrology, the development of antidunes, the formation of geological deposits from hyper-concentrated flows, and shallow-water flow modeling in hydraulics. In this work, Dressler’s theory is developed in an alternative way by a systematic iteration of the stream and potential functions in terrain-fitted coordinates. The first iteration was found to be the former Dressler’s theory, whereas a second iteration of the governing equations results in velocity components generalizing Dressler’s theory to wave-like motion. Dressler’s first-order theory produces a transcritical flow solution over topography only if the total head is fixed by a minimum value of the specific energy at the transition point. However, the theory deviates from measurements under subcritical flow conditions, given that the bed-normal velocity component is significant. A second iteration to the velocity field was used to produce a second-order differential equation that resembles the cnoidal-wave theory. It accurately produces flow over an obstacle including the critical point and the minimum specific energy as part of the numerical solution. The new cnoidal-wave model compares well with the theory of a Cosserat surface for directed fluid sheets, whereas the Saint-Venant theory appears to be poor.     

Turbulent velocity profile in fully-developed open channel flows

The determination of velocity profile in turbulent narrow open channels is a difficult task due to the significant effects of the anisotropic turbulence that involve the Prandtl’s second type of secondary flow occurring in the cross section. With these currents the maximum velocity appears below the free surface that is called dip phenomenon. The well-known logarithmic law describes the velocity distribution in the inner region of the turbulent boundary layer but it is not adapted to define the velocity profile in the outer region of narrow channels. This paper relies on an analysis of the Navier–Stokes equations and yields a new formulation of the vertical velocity profile in the center region of steady, fully developed turbulent flows in open channels. This formulation is able to predict time averaged primary velocity in the outer region of the turbulent boundary layer for both narrow and wide open channels. The proposed law is based on the knowledge of the aspect ratio and involves a parameter C_Ar depending on the position of the maximum velocity (ξ_dip). ξ_dip may be derived, either from measurements or from an empirical equation given in this paper. A wide range of longitudinal velocity profile data for narrow open channels has been used for validating the model. The agreement between the measured and the computed velocities is rather good, despite the simplification used.     

Prediction of particulate air pollution from a landfill site using CFD and LIDAR techniques

Landfill sites have been the most common way of eliminating solid urban waste, as well as that of public and mining wastes. Landfill sites are a constant source of environmental pollution and wind is the most important contributing factor to air pollution, due to the erosion which it produces over the landfill site surfaces, transporting dust away from the source point. This causes serious risks for human health and general dirt in the surrounding areas of the landfill site. The result of dust air pollution produced in a landfill site is analysed through CFD3D (Computational Fluid Dynamics) by joining the LIDAR (Light Detection and Ranging) technique and Ansys CFX 10.0 software. The CFD simulations determine the wind velocity distribution on the landfill site surface and the different particle threshold friction velocities which determine the dust emission in multiphase simulations (air-particles). These simulations are validated from field data obtained in three measurement programmes for each type of landfill site surface treatment which has been studied. It was determined that the superficial landfill site treatment with the lowest air pollution is tall grass and bushes. The methodology used can be applied to the dust emission calculation in the design or evaluation of other landfill sites.     

Coastal erosion and shore protection: A brief historical analysis

Since the Neolithic, humans have gathered along coastal plains, where they had to face sea level rise and subsidence without the technology to oppose these processes. When sea level stabilized, approx. 6.000 yr. B.P., coastal colonization was allowed, but where mountain deforestation was carried out river sediment input increased tremendously: settlements were disconnected from the shore and harbour siltation occurred. Shore erosion was a limited process at the time and local solutions were adopted: clay dikes, wood piles, fascinates and rock revetments. Along the Mediterranean, in China and Japan the construction of more complex structures has been documented since the Middle Ages. Further human development, with river bed quarrying, wetland reclamation, dam construction and mountain re-afforestation, favoured a coastal erosion that nowadays threatens most world coasts. From the Venetian “Murazzi” to the recent 114-km-long concrete element defence at the Yellow river delta, shore protection structures evolved following the different needs (protect coastal communication routes, urban and industrial settlements, tourist resorts). Beach nourishment, previously performed with inland quarried materials, is now intensively carried out with marine aggregates. The vernacular solution, left to undeveloped countries, has been recently revaluated by “green engineering”. However, with Sea Level Rise, the debate of whether to defend, accommodate or retreat is open.     

华南红壤坡面侵蚀水动力学机制试验研究

华南红壤区是我国发生潜在水土流失危害的重点区域之一,红壤坡地的水土流失现象较为严重。华南红壤地区气候温暖,雨量丰沛,降雨特点通常表现为强度大、历时长,暴雨或大暴雨事件频发,因此暴雨驱动因子影响作用更为显著,研究探讨华南红壤坡面侵蚀过程及动力特征有着重要的现实意义。通过室内35场不同雨强（30、60、90、120、180、210、270 mm·h-1）、不同坡度（5°、10°、15°、20°、25°）的人工模拟降雨试验,研究了华南红壤侵蚀的水动力学机制。主要取得以下研究结果：（1）红壤坡面侵蚀水流的流态大部分均处于＂层流—缓流区＂,除平均流速随着模拟雨强的增加而增加外,其他水动力参数（坡面侵蚀水流的阻力、曼宁糙率、径流剪切力、以及水流功）均会受到雨强和坡度改变的影响,但是不呈现特定的函数关系,规律不明显;（2）从水动力学角度解析了红壤坡面水蚀过程,得出了相对水深和曼宁糙率系数两种水动力因子共同组成的复合水动力特征参数,可作为表征不同雨强及坡度下的坡面水蚀侵蚀产沙的特征水动力参数指标。