Numerical modelling of horizontal sediment-laden jets

Sediment-laden turbulent flows are commonly encountered in natural and engineered environments. It is well known that turbulence generates fluctuations to the particle motion, resulting in modulation of the particle settling velocity. A novel stochastic particle tracking model is developed to predict the particle settling out and deposition from a sediment-laden jet. Particle velocity fluctuations in the jet flow are modelled from a Lagrangian velocity autocorrelation function that incorporates the physical mechanism leading to a reduction of settling velocity. The model is first applied to study the settling velocity modulation in a homogeneous turbulence field. Consistent with basic experiments using grid-generated turbulence and computational fluid dynamics (CFD) calculations, the model predicts that the apparent settling velocity can be reduced by as much as 30 % of the stillwater settling velocity. Using analytical solution for the jet mean flow and semi-empirical RMS turbulent velocity fluctuation and dissipation rate profiles derived from CFD predictions, model predictions of the sediment deposition and cross-sectional concentration profiles of horizontal sediment-laden jets are in excellent agreement with data. Unlike CFD calculations of sediment fall out and deposition from a jet flow, the present method does not require any a priori adjustment of particle settling velocity.     

Clustering of aerosols in atmospheric turbulent flow

A mechanism of formation of small-scale inhomogeneities in spatial distributions of aerosols and droplets associated with clustering instability in the atmospheric turbulent flow is discussed. The particle clustering is a consequence of a spontaneous breakdown of their homogeneous space distribution due to the clustering instability, and is caused by a combined effect of the particle inertia and a finite correlation time of the turbulent velocity field. In this paper a theoretical approach proposed in Elperin et al. (2002) Phys Rev E 66:036302 is further developed and applied to investigate the mechanisms of formation of small-scale aerosol inhomogeneities in the atmospheric turbulent flow. The theory of the particle clustering instability is extended to the case when the particle Stokes time is larger than the Kolmogorov time scale, but is much smaller than the correlation time at the integral scale of turbulence. We determined the criterion of the clustering instability for the Stokes number larger than 1. We discussed applications of the analyzed effects to the dynamics of aerosols and droplets in the atmospheric turbulent flow.     

Evaluation of CFD Simulation using RANS Turbulence Models for Building Effects on Pollutant Dispersion

CFD evaluations were performed to examine the applicability of the RANS methods in simulating pollutant dispersion near, within and over three typical building configurations: (1) an isolated building, (2) a building array and (3) an urban intersection. The CFD results are compared with values obtained from wind tunnel tests. In some situations major differences between the wind tunnel tests and the CFD results were observed. The main source of difference between the CFD and wind tunnel results was inadequate modelling of local flow patterns using the RANS turbulence models. Also inappropriate evaluation of high intermittent turbulent mixing in the RANS approach may lead to either over-prediction or under-prediction of the concentration level, by up to a factor of 10, depending on the case investigated.     

Deposition and resorption of organic aerosol constituents in the respiratory tract as calculated from particle size distribution data

Data on the particle size distributions of organic aerosol constituents were used as input for a study, designed to calculate the fractions of the particulate concentrations of these compounds, deposited into the respiratory tract. The known relation between the deposition probability and the particle size as described by the ICRP‐model was used. The organic constituents were from the classes of the aliphatic hydrocarbons, carboxylic acids, polycyclic and aza‐heterocyclic aromatic hydrocarbons. Aerosol samples were obtained by Hi‐Vol cascade impactor sampling at suburban, rural and sea shore background stations as well as in an industrial emission site (coke oven).

Our approach uses the measured concentrations, being average values within each impactor particle size interval, as well as the integrated average deposition probabilities. This procedure was validated experimentally for eight model distributions from the literature, for which an infinitesimal calculation of the deposited fractions was possible.

Dilution reduces total particulate concentrations in the remote areas and predominantly determines the total deposited pollutant concentrations. Of these, pulmonary and nasopharyngeal deposition are most significant and, as a first approximation, correspond largely to the relative importance of the accumulation and dispersion modes of the sampled aerosol. A particle size distribution shift toward larger particles within the accumulation mode occurs upon ageing of the aerosol and reduces the pulmonary deposited fraction of the measured compounds in the background sites, compared to the one in the suburb. The total deposited fraction, however, increases. The contributions of biogenic higher odd n‐alkanes and, to a lesser extent, of even carboxylic acids to the dispersion mode of the aerosol result in an increased nasopharyngeal deposition at the background sites mainly during summer.

Since little information on the bio‐availability of organic aerosol constituents is available in the literature, the fractions of the particulate pollutant concentrations, resorbed in the tissues from the deposited material, were calculated, assuming an average efficiency of 70% for pulmonary and of 10% for nasopharyngeal and tracheobronchial resorption. A nearly constant total resorbed fraction of 20±2% resulted, independent of the sampling station or the season chosen, in contrast with the total deposited fractions, for which significant differences were observed. The predominant pulmonary resorption as well as compensating effects of the nasopharyngeal resorption level out the relatively small differences in particle size distributions observed. Based on these data, a first estimate in nanogram of the daily intake by inhalation of the organic pollutants studied can be formulated as four times the particulate pollutant concentration, expressed in ngm^‐3.     

A Computational Fluid Dynamics Approach for Urban Area Transport and Dispersion Modeling

A simulation tool has been developed to model the wind fields, turbulence fields, and the dispersion of Chemical, Biological, Radiological and Nuclear (CBRN) substances in urban areas on the building to city blocks scale. A Computational Fluid Dynamics (CFD) approach has been taken that naturally accounts for critical flow and dispersion processes in urban areas, such as channeling, lofting, vertical mixing and turbulence, by solving the steady-state, Reynolds-Averaged Navier–Stokes (RANS) equations. Rapid generation of high quality cityscape volume meshes is attained by a unique voxel-based model generator that directly interfaces with common Geographic Information Systems (GIS) file formats. The flow and turbulence fields are obtained by solving the steady-state RANS equations using a collocated, pressure-based approach formulated for unstructured and polyhedral mesh elements. Turbulence modeling is based upon the Renormalization Group variant of the k–ε model (k–ε RNG). Neutrally buoyant simulations are made by prescribing velocity boundary condition profiles found by a power–law relationship, while turbulence quantities boundary conditions are defined by a prescribed mixing length in conjunction with the assumption of turbulence equilibrium. Dispersion fields are computed by solving an unsteady transport equation of a dilute gas, formulated in a Eulerian framework, using the velocity and turbulence fields found from the steady-state RANS solution. In this paper the model is explained and detailed comparisons of predicted to experimentally obtained velocity, turbulence and dispersion fields are made to neutrally stable wind tunnel and hydraulic flume experiments.     

Characterization of bedload intermittency near the threshold of motion using a Lagrangian sediment transport model

At the smallest scales of sediment transport in rivers, the coherent structures of the turbulent boundary layer constitute the fundamental mechanisms of bedload transport, locally increasing the instantaneous hydrodynamic forces acting on sediment particles, and mobilizing them downstream. Near the critical threshold for initiating sediment motion, the interactions of the particles with these unsteady coherent structures and with other sediment grains, produce localized transport events with brief episodes of collective motion occurring due to the near-bed velocity fluctuations. Simulations of these flows pose a significant challenge for numerical models aimed at capturing the physical processes and complex non-linear interactions that generate highly intermittent and self-similar bedload transport fluxes. In this investigation we carry out direct numerical simulations of the flow in a rectangular flat-bed channel, at a Reynolds number equal to Re = 3632, coupled with the discrete element method to simulate the dynamics of spherical particles near the bed. We perform two-way coupled Lagrangian simulations of 48,510 sediment particles, with 4851 fixed particles to account for bed roughness. Our simulations consider a total of eight different values of the non-dimensional Shields parameter to study the evolution of transport statistics. From the trajectory and velocity of each sediment particle, we compute the changes in the probability distribution functions of velocities, bed activity, and jump lengths as the Shields number increases. For the lower shear stresses, the intermittency of the global bedload transport flux is described by computing the singularity or multifr actal spectrum of transport, which also characterizes the widespread range of transport event magnitudes. These findings can help to identify the mechanisms of sediment transport at the particle scale. The statistical analysis can also be used as an ingredient to develop larger, upscaled models for predicting mean transport rates, considering the variability of entrainment and deposition that characterizes the transport near the threshold of motion.     

Prediction of the Transport and Dispersal of Oil in the South Caspian Sea Resulting from Blowouts

A 3-D hybrid flow/transport model is developed to predict the dispersal of oil pollution in coastal waters. The transport module of the model takes predetermined current and turbulent diffusivities and uses Lagrangian tracking to predict the motion of individual particles (droplets), the sum of which constitutes a hypothetical oil spill. Currents and turbulent diffusivities used in the model are generated by a numerical ocean circulation model (POM) implemented for the Caspian Sea. The basic processes affecting the fate of the oil spill are taken into account and parameterised in the transport model. The process of evaporation is modeled with the pseudo-component approach. The model is implemented for a simulated continuous release in the coastal waters of the south part of the Caspian Sea. Numerical experiments simulate 5- and 10-day blowout scenarios resulting from sources situated in areas were intensive and extensive development of oil deposits is expected soon. Oil slick movement and risk of coastline contamination by beaching of offshore oil spills are illustrated for different wind conditions.     

Deterministic particle tracking simulation of pollutant discharges in rivers and estuaries

A two-dimensional deterministic particle tracking model, in which the anisotropic-dispersive process is described by a particle strength exchange scheme, is established for the simulation of pollutant transport in vertically well-mixed rivers and estuaries. By simulating two benchmark problems with analytic solutions, the PSE scheme is shown to be accurate even if the anisotropic ratio of dispersion coefficients is very high. Further simulations of two specific problems concerning the optimal effluent discharge location and procedure are presented. The major conclusion is that in a tidal estuary with a relatively large fresh-water flow, setting the discharge position at the mixing center and making the discharge rate proportional to flow speed may minimize the peaks of concentration.     

排水管道内不同粒径沉积颗粒物冲刷率的分析与计算

为探究和量化在水流冲刷下,排水管道中不同粒径颗粒物冲刷沉积的过程,本文模拟排水管道内沉积颗粒的冲刷过程.冲刷过程中,粒径较小的悬移质颗粒(小于0.1 mm),在管道沿线取样测得水流中悬浮固体质量浓度(SS);粒径较大的推移质颗粒(0.1—2 mm),测得管道沿线不同位置沉积的颗粒质量.并建立两个数学模型分别计算排水管道内两类颗粒物的冲刷率.研究发现,悬移质颗粒,以水流中SS为依据,可计算出输送通量和冲刷率;推移质颗粒,以管道不同位置沉积下来的颗粒质量为依据,拟合出了管道中的沉积分布函数,是以e为底数的指数函数,可计算出冲刷量和冲刷率.模型计算出的冲刷率呈现出规律性:悬移质颗粒被冲刷成悬浮状态,随水流迁移过程中部分会再次沉积,使得悬移质颗粒的冲刷率从管道前段至中后段逐步降低,如初始沉积质量为100 g的0.045 mm悬移质颗粒,在0.30 m·s^-1的冲刷流速下,计算出管道前段冲刷率为78.94%,最终在管道后段降至13.89%;对于两类颗粒而言,颗粒物粒径越小,冲刷流速越大,初始沉积质量越小,冲刷率越高.     

Simulation-based optimization of in-stream structures design: rock vanes

We employ a three-dimensional coupled hydro-morphodynamic model, the Virtual Flow Simulator (VFS-Geophysics) in its Unsteady Reynolds Averaged Navier–Stokes mode closed with \(k-\omega\) model, to simulate the turbulent flow and sediment transport in large-scale sand and gravel bed waterways under prototype and live-bed conditions. The simulation results are used to carry out systematic numerical experiments to develop design guidelines for rock vane structures. The numerical model is based on the Curvilinear Immersed Boundary approach to simulate flow and sediment transport processes in arbitrarily complex rivers with embedded rock structures. Three validation test cases are conducted to examine the capability of the model in capturing turbulent flow and sediment transport in channels with mobile-bed. Transport of sediment materials is handled using the Exner equation coupled with a transport equation for suspended load. Two representative meandering rivers, with gravel and sand beds, respectively, are selected to serve as the virtual test-bed for developing design guidelines for rock vane structures. The characteristics of these rivers are selected based on available field data. Initially guided by existing design guidelines, we consider numerous arrangements of rock vane structures computationally to identify optimal structure design and placement characteristics for a given river system.     

Numerical study of aerosols motion and deposition outside and inside buildings with different geometries

In this study, motion and deposition of various sizes of aerosols in turbulent flow of air inside and outside of two-dimensional buildings with closed and open windows have been investigated, numerically. This simulation was based on FVM solving of RANS equations with k–ε model. Eulerian–Lagrangian method was used to simulate fluid and particles motions, respectively, and one-way coupling between them was considered. Effect of particles size on deposition efficiency has been calculated. The results shown that, the particles deposition outside and inside of the building with domical roof is less than triangular and flat roof buildings.     

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.     

Stormwater treatment: examples of computational fluid dynamics modeling

Control of rainfall-runoff particulate matter (PM) and PM-bound chemical loads is challenging; in part due to the wide gradation of PM complex geometries of many unit operations and variable flow rates. Such challenges and the expense associated with resolving such challenges have led to the relatively common examination of a spectrum of unit operations and processes. This study applies the principles of computational fluid dynamics (CFD) to predict the particle and pollutant clarification behavior of these systems subject to dilute multiphase flows, typical of rainfall-runoff, within computationally reasonable limits, to a scientifically acceptable degree of accuracy. The Navier-Stokes (NS) system of nonlinear partial differential equations for multiphase hydrodynamics and separation of entrained particles are solved numerically over the unit operation control volume with the boundary and initial conditions defined and then solved numerically until the desired convergence criteria are met. Flow rates examined are scaled based on sizing of common unit operations such as hydrodynamic separators (HS), wet basins, or filters, and are examined from 1 to 100 percent of the system maximum hydraulic operating flow rate. A standard turbulence model is used to resolve flow, and a discrete phase model (DPM) is utilized to examine the particle clarification response. CFD results closely follow physical model results across the entire range of flow rates. Post-processing the CFD predictions provides an in-depth insight into the mechanistic behavior of unit operations by means of three dimensional (3-D) hydraulic profiles and particle trajectories. Results demonstrate the role of scour in the rapid degradation of unit operations that are not maintained. Comparisons are provided between measured and CFD modeled results and a mass balance error is identified. CFD is arguably the most powerful tool available for our profession since continuous simulation modeling.     

Modelling of turbulent dispersion for numerical simulation of wind-driven rain on bridges

Wind-driven rain (WDR) is responsible for many potential negative effects on bridges, such as structural cracking, aggregate erosion, steel corrosion and storm water management problems and so on. Hence, accurate evaluations of the WDR effects on bridges are essential to provide solutions for preventing material degradation and improving durability capability of bridges. However, in most previous WDR numerical studies, the turbulent dispersion of raindrops was neglected. In this paper, the turbulent dispersion is integrated into Eulerian multiphase model to investigate the WDR effects on a bridge with rectangular cross-section. Especially, the influences of the turbulent dispersion are discussed in detail by comparing the WDR simulation results for the cases with and without consideration of the turbulent dispersion in terms of WDR flow fields, volume fraction, specific catch ratio, catch ratio, rain loads and aerostatic force coefficients. The results indicate that the turbulent dispersion for a certain range of raindrop size is needed to be taken into account for obtaining accurate WDR simulation results for bridges.     

Self-contamination of aquaculture cages in shallow water

Fish farms, which initially colonized quiet and protected natural coastal areas, are now frequently installed in open flow zones, due to the lack of space along coasts and to the emergence of new environmental constraints. For the past two decades, a salmon fish farm has been located inside the roadstead of Cherbourg (France) to benefit from both sea protection and tide currents which regularly refresh the water. In spite of these favourable environmental conditions, periods of non-negligible fish mortalities have been observed to occur without clear evidence of their origin. This motivated the turbidity measurements and the numerical simulations presented in this paper. Firstly, it is shown that high turbidities in the farm site under study are mainly due to the flow acceleration under the cages, which causes the re-suspension of sediments and bio-deposits. Secondly, particles which enter the fishnet can have different origins (external source, bottom, or the net itself). Numerical simulations, based on the Reynolds equations and on the discrete random walk model for particle dispersion, suggest that the rear area of the net can be reached by particles emerging from below the net. It is observed that turbulent dispersion is a key ingredient for such a behaviour, as it can lead particles towards a large recirculation cell behind the net. Dispersion by realistic unsteady vortices has also been analysed by means of a Lattice-Boltzmann model. Though these computations involve smaller Reynolds numbers, they confirm qualitatively the observations of the random walk model. In addition, they suggest that vortex shedding and unsteady recirculation cells near the bottom can force particles from the sand bed to be lifted up and reach the rear of the net.     

Two-phase flow insights into open-channel flows with suspended particles of different densities

The effect of particle density on the turbulent open-channel flow carrying dilute particle suspensions is investigated using two specific gravities and three concentrations of solid particles. The particles, identical in size and similar in shape, were natural sand and a neutrally buoyant plastic. The particles were fully suspended, and formed no particle streaks on the channel’s bed. Accordingly, the changes in the flow are attributed to the interactions between suspended particles and flow turbulence structures. Measurements were obtained by means of image velocimetry enabling simultaneous, but distinct, measurement of liquid and particle velocities. The experimental results show that, irrespective of particle specific gravity, particle suspension influences bulk velocity of flow and the Kármán coefficient, while friction velocity essentially remains constant. The results also show that particles in suspension modify local water turbulence over the flow depth, but in ways not accurately predicted using the customary parameters for characterizing turbulence modification.     

LES validation of urban flow,part II: eddy statistics and flow structures

Time-dependent three-dimensional numerical simulations such as large-eddy simulation (LES) play an important role in fundamental research and practical applications in meteorology and wind engineering. Whether these simulations provide a sufficiently accurate picture of the time-dependent structure of the flow, however, is often not determined in enough detail. We propose an application-specific validation procedure for LES that focuses on the time dependent nature of mechanically induced shear-layer turbulence to derive information about strengths and limitations of the model. The validation procedure is tested for LES of turbulent flow in a complex city, for which reference data from wind-tunnel experiments are available. An initial comparison of mean flow statistics and frequency distributions was presented in part I. Part II focuses on comparing eddy statistics and flow structures. Analyses of integral time scales and auto-spectral energy densities show that the tested LES reproduces the temporal characteristics of energy-dominant and flux-carrying eddies accurately. Quadrant analysis of the vertical turbulent momentum flux reveals strong similarities between instantaneous ejection-sweep patterns in the LES and the laboratory flow, also showing comparable occurrence statistics of rare but strong flux events. A further comparison of wavelet-coefficient frequency distributions and associated high-order statistics reveals a strong agreement of location-dependent intermittency patterns induced by resolved eddies in the energy-production range. The validation concept enables wide-ranging conclusions to be drawn about the skill of turbulence-resolving simulations than the traditional approach of comparing only mean flow and turbulence statistics. Based on the accuracy levels determined, it can be stated that the tested LES is sufficiently accurate for its purpose of generating realistic urban wind fields that can be used to drive simpler dispersion models.     

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.     

Kohäsive Feinpartikel in fluvialen Systemen: einige Gedanken zu Erkenntnissen und Forschungsdefiziten

Fluvial cohesive sediment is of fundamental environmental and multidisciplinary concern owing to its significant impact on energy, nutrient and trace element fluxes, river sediment budgets as well as habitat quality. Consequently, numerous studies in geomorphology, hydraulics, hydrology and river ecology accentuate the importance of fluvial fine grained sediment. Increased deposition of fine grained particles for instance can negatively effect benthic habitats. Settling velocity and deposition rates of fines are key terms in sediment transport models. Thus, the knowledge of transport and quantitative deposition dynamics are necessary prerequisites for a sustainable sediment and habitat management as well as a reliable, physically based modelling of fine sediment conveyance and contaminant transfer in streams. However, it is generally a challenging task to characterise and trace fine-grained particle transport, deposition and resuspension dynamics in a fluvial environment. One inherent complexity of fluvial cohesive sediment dynamics is for instance the relevance of interparticle forces and flocculation/aggregation processes. Owing to this complexity and the lack of standardized methods cohesive particle transport mechanisms are hardly predictable in a specific field situation. Consequently, interdisciplinary research approaches are needed. This article discusses specific findings and research gaps in the field of cohesive sediment research.     

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.