LES-DEM simulations of sediment transport

In this work, a fully-coupled Computational Fluid Dynamics (CFD) model and Discrete Element Method (DEM) are used to simulate a unidirectional turbulent open-channel flow over the full range of sediment transport regimes. The fluid and particles are computed on separate grids using a dual-grid formulation to maintain consistency and avoid instability issues. The results of coupling the dispersed phase to a multiphase flow solver that uses volume-averaged Navier-Stokes equations are compared to those obtained from coupling through drag to a single flow solver. The current work also examines the applicability and limitations of lumping particles as a representative particle to reduce the cost of simulations. Insight to the impact of different turbulent events to the entrainment of particles is also given. The simulation results of sediment transport from both coupling techniques show good agreement with empirical formulas in the bedload regime, but under-predict sediment transport in the suspended load regime. In the suspended load regime, using partial coupling, the rate of sediment transport was found to be under-predicted as compared to full-coupling. The deviation in results in the suspended load regime was found to increase with increases in the applied shear stress. Both coupling methods revealed the same effect on the friction factor where friction increases in the bedload regime and decreases in the suspended load regime reaching a maximum at the transition between regimes. This result is contrary to past studies which have shown a discrete jump in the friction factor at the transition. Lumping particles as representative particles is shown to reduce the simulation cost by more than a factor of 5 when using a scaling factor of 2. By doing a quadrant analysis on information obtained from particle and flow field results, it was found that most of the particles are entrained by more frequent sweep events.     

An improved method for evaluating the seasonal variability of total suspended sediment flux field in the Yellow and East China Seas

The suspended sediment flux field in the Yellow and East China Seas(YECS) displays its seasonal variability.A new method is introduced in this paper to obtain the flux field via retrieval of ocean color remote sensing data,statistical analysis of historical suspended sediment concentration data,and numerical simulation of three-dimensional(3D) flow velocity.The components of the sediment flux field include(i) surface suspended sediment concentration inverted from ocean color remote sensing data;(ii) vertical distribution of suspended sediment concentration obtained by statistical analysis of historical observation data;and(iii) 3D flow field modeled by a numerical simulation.With the improved method,the 3D suspended sediment flux field in the YECS has been illustrated.By comparison with the suspended sediment flux field solely based on the numerical simulation of a suspended sediment transport model,the suspended sediment flux field obtained by the improved method is found to be more reliable.The 3D suspended sediment flux field from ocean colour remote sensing and in situ observation are more closer to the reality.Furthermore,by quantitatively analyzing the newly obtained suspended sediment flux field,the quantity of sediment erosion and deposition within the different regions can be evaluated.The sediment exchange between the Yellow Sea and the East China Sea can be evident.The mechanism of suspended sediment transport in the YECS can be better understood.In particular,it is suggested that the long-term transport of suspended sediment is controlled mainly by the circulation pattern,especially the current in winter.     

Sand settling through bedform-generated turbulence in rivers

Fluvial bedforms generate a turbulent wake that can impact suspended-sediment settling in the passing flow. This impact has implications for local suspended-sediment transport, bedform stability, and channel evolution; however, it is typically not well-considered in geomorphologic models. Our study uses a three-dimensional OpenFOAM hydrodynamic and particle-tracking model to investigate how turbulence generated from bedforms and the channel bed influences medium sand-sized particle settling, in terms of the distribution of suspended particles within the flow field and particle-settling velocities. The model resolved the effect of an engineered bedform, which altered the flow field in a manner similar to a natural dune. The modelling scenarios alternated bed morphology and the simulation of turbulence, using detached eddy simulation (DES), to differentiate the influence of bedform-generated turbulence relative to that of turbulence generated from the channel bed. The bedform generated a turbulent wake that was composed of eddies with significant anisotropic properties. The eddies and, to a lesser degree, turbulence arising from velocity shear at the bed substantially reduced settling velocities relative to the settling velocities predicted in the absence of turbulence. The eddies tended to advect sediment particles in their primary direction, diffuse particles throughout the flow column, and reduced settling likely due to production of a positively skewed vertical-velocity fluctuation distribution. Study results suggest that the bedform wake has a significant impact on particle-settling behaviour (up to a 50% reduction in settling velocity) at a scale capable of modulating local suspended transport rates and bedform dynamics. © 2020 John Wiley & Sons, Ltd.     

Flow analysis in a channel with flexible vegetation using double-averaging method

The paper addresses the problem of the resistance due to vegetation in an open channel flow, characterized by partially and fully submerged vegetation formed by colonies of bushes. The flow is characterized by significant spatial variations of velocity between vertical profiles that make the traditional approach based on time averaging of turbulent fluctuations inconvenient. A more useful procedure, based on time and spatial averaging (Double-Averaging Method) is applied for the flow field analysis and characterization. The vertical distribution of mean velocity and turbulent stresses at different spatial locations has been measured with a 3D Acoustic Doppler Velocimeter (ADV) for two different vegetation densities where fully submerged real bushes (salix pentandra) have been used. Velocity measurements were completed together with the measurements of drag exerted on the flow by bushes at different flow depths. The analysis of velocity measurements allows depicting the fundamental characteristics of both the mean flow field and turbulence. The experimental data show that the contribution of form-induced stresses to the momentum balance cannot be neglected. The mean velocity profiles and the spatially averaged turbulent intensity profiles allow inferring that the vegetation density is a driving parameter for the development of a mixing layer at the canopy top in the case of submerged vegetation. Moreover, the net upward turbulent momentum flux, evaluated with the methodology proposed by Lu and Willmarth (1973), appears to be damped for increased vegetation density; this finding can rationally explain the reduction of the suspended sediment transport capacity typically observed in free surface flows over a vegetated bed.     

What is suspended sediment?

Suspended sediment is conventionally regarded as that sediment transported by a fluid that it is fine enough for turbulent eddies to outweigh settling of the particles through the fluid. Early work in the fluvial field attributed suspension to turbulence, and led to the notion of a critical threshold for maintaining sediment in suspension. However, research on both turbulence structures and the interactions between suspended sediment and bedforms in rivers has shown a more complex story and, although there appear to have been no studies of the impact of bedforms on aeolian suspended sediment concentrations, turbulent flow structures and transport rates of saltating particles have been shown to be affected. This research indicates that suspended sediment neither travels with the same velocity as the flow in which it is suspended, nor is it likely to remain in suspension in perpetuity, even under conditions of steady flow or in unsteady flow the where dimensionless critical threshold is permanently exceeded. Rather, like bedload, it travels in a series of hops, and is repeatedly deposited on the bed where it remains until it is re‐entrained. Is there, therefore, a qualitative difference between suspended and saltating sediment, or is it just a quantitative difference in the size of the jump length and the frequency of re‐entrainment? It is our contention that the distinction of suspension as a separate class of sediment transport is both arbitrary and an unhelpful anthropocentric artefact. If we recognize that sediment transport is a continuum and applies to any fluid medium rather than split into different “processes” based on arbitrary thresholds and fluids, then recognizing the continuity will enable development of an holistic approach sediment transport, and thus sediment‐transport models that are likely to be viable across a wider range of conditions than hitherto. Copyright © 2015 John Wiley & Sons, Ltd.     

Application of Hausdorff fractal derivative to the determination of the vertical sediment concentration distribution

The Rouse formula and its variants have been widely used to calculate the steady-state vertical concentration distribution for suspended sediment in steady sediment-laden flows, where the diffusive flux is assumed to be Fickian. Turbulent flow, however, exhibits fractal properties, leading to non-Fickian diffusive flux for sediment particles. To characterize non-Fickian dynamics of suspended sediment, the current study proposes a Hausdorff fractal derivative based advection-dispersion equation(H...     

Multiscale modeling of the atmospheric environment over a forest canopy

Vegetation constitutes one of the fundamental types of land use on Earth. The presence of trees in urban areas can decrease local winds and exchange sensible and latent heat with the surrounding environments, thus exerting notable influences on the urban microenvironment. A better understanding of the turbulent transfer of momentum and scalars around vegetation canopy could significantly contribute to improvement of the urban environment. This work develops a large-eddy simulation(LES) method that is applicable to model the flow and scalar transport over the forest canopy. We study the atmospheric flow over complex forested areas under typical weather conditions by coupling LES to the mesoscale model. Models of radiation and energy balance have been developed with explicit treatment of the vegetation canopy. By examining the flow over a forest canopy under a range of stability conditions, we found that buoyancy enhances or suppresses turbulent mixing in unstable or stable atmosphere respectively, with decreasing or increasing wind shear, respectively. From the multiscale modeling of the Beijing Olympic Forest Park, the present coupling scheme proves to better resolve the diurnal variations in wind speed,temperature, and relative humidity over complex urban terrains. The coupling scheme is superior to the traditional mesoscale model in terms of wind field simulation. This is mainly because the coupling scheme not only takes the influences of external mesoscale flow into consideration, but also resolves the heterogeneous urban surface at a fine scale by downscaling, thus better reproducing the complex flow and turbulent transport in the urban roughness sublayer.     

Comparing the transitional behaviour of kaolinite and bentonite suspension flows

Past research has demonstrated the dramatic effects that variations in suspended clay can have on the properties of flow by producing a range of transitional flows between turbulent and laminar states, depending on clay concentration and fluid shear. Past studies have been restricted to kaolinite flows, a clay mineral that has relatively weak cohesive properties. This paper extends these studies to suspension flows of bentonite, a clay mineral that attains higher viscosities at far lower volumetric concentrations within a flow. The results show that the types of transitional flow behaviour recognized in past studies can also be found in bentonite suspension flows, but at lower suspended sediment concentrations, thus demonstrating an even more dramatic effect on flow properties, and potentially on sediment transport and resulting bed morphology, than kaolinite flows. The paper proposes new stability diagrams for the phase space of bentonite flows and compares these to past work on kaolinite suspension flows. These new data suggest that the transitional‐flow Reynolds number can be used to delineate the types of transitional flow across different clay types and assess modern and ancient clay‐suspension flows. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Finite element computation of turbulent flows

After reviewing briefly the major alternatives in the mathematical modelling of turbulent flows, this paper combines the kinetic energy-dissipation (k−ε) turbulent flow model and a mean flow model in a Galerkin finite element model applicable to incompressible turbulent two-dimensional flows with recirculation. Application of this model to examination of flow and particle settling in a rectangular sedimentation basin is presented; this approach completely determines the non-uniform velocity field and also resolves much uncertainty about the choice of values for both the eddy viscosity and the turbulent transport coefficients for settling. Finally, the expected directions of future research work in finite element simulation of turbulent flows are delineated.     

Burst-like sediment suspension events in a sand bed river

In a much quoted paper, Jackson (1976) hypothesized that turbulent [bursting] motions such as those documented in laboratory boundary layers play a major role in alluvial sediment suspension. To date, the hypothesis remained largely untested, due to difficulties in monitoring turbulent suspension in rivers. This study provides field data documenting burst-like turbulent motions over a sandy bed channel and quantifying the role of these motions in sand suspension. The data were collected in a 10 m deep channel of the Fraser River near Mission, British Columbia, Canada. Turbulent fluctuations of both flow components, downstream and normal to the bed, along with the output of an optical suspended sediment sensor, were monitored 1 m above the river bed. Typical flow velocities averaged 0·9 ms⁻¹ at the sensors, where mean suspended sediment concentrations were 500 mgl⁻¹; decimetre height small dunes on the backs of larger, half-metre amplitude dunes covered the channel bed in the area. Brief but intense, burst-like [ejection and inrush] events were identified in the flow records, where they are responsible for a high degree of [intermittency] in shear stress over the dunes: 80 per cent of the turbulent momentum exchange across the 1 m level can be ascribed to such brief (3-8 s duration) events, active under 12 per cent of the time. In addition, the record of fluctuating sediment concentrations reveals these burst-like motions to be highly effective in vertically mixing suspended sediment and thus, ultimately, in maintaining suspended sediment transport above the dune bed. The bulk (60 and 90 per cent in two deployments) of the vertical sediment mixing was accomplished by intense events active some 10 per cent of the time. No discrete recurrence timescale for these ‘burst-like’ mixing events is evident, however. Rather, a continuous variation of return periods was observed as a function of the magnitude of vertical mixing event considered. To that extent, conceptual models of sediment transport in terms of burst events with a predictable recurrence such as proposed by Jackson (1976) may be misleading.     

Large-eddy simulation of the turbulent structure in compound open-channel flows

A large-eddy simulation study has been undertaken to investigate the turbulent structure of open-channel flow in an asymmetric compound channel. The dynamic sub-grid scale model has been employed in the model, with the partial cell treatment being implemented using a Cartesian grid structure to deal with the floodplain. The numerical model was used to predict the: primary velocity and secondary currents, boundary shear stress, turbulence intensities, turbulent kinetic energy, and Reynolds stresses. These parameters were compared with experimental measurements published in the literature, with relatively close agreement being obtained between both sets of results. Furthermore, instantaneous flow fields and large-scale vortical structures were predicted and are presented herein. These vortical structures were found to be responsible for the significant lateral exchange of mass and momentum in compound channels.     

Modelling of hydrodynamics and cohesive sediment transport in Tanshui River estuarine system,Taiwan

A laterally averaged two-dimensional numerical model is used to simulate hydrodynamics and cohesive sediment transport in the Tanshui River estuarine system. The model handles tributaries as well as the main stem of the estuarine system. Observed time series of salinity data and tidally averaged salinity distributions have been compared with model results to calibrate the turbulent diffusion coefficients. The overall model verification is achieved with comparisons of residual currents and salinity distribution. The model reproduces the prototype water surface elevation, currents and salinity distributions. Comparisons of the suspended cohesive sediment concentrations calculated by the numerical model and the field data at various stations show good agreement. The validated model is applied to investigate the tidally averaged salinity distributions, residual circulation and suspended sediment concentration under low flow conditions in the Tanshui River estuarine system. The model results show that the limit of salt intrusion in the mainstem estuary is located at Hsin-Hai bridge in Tahan Stream, 26 km from the River mouth under Q75 flow. The null point is located at the head of salt intrusion, using 1 ppt isohaline as an indicator. The tidally averaged sediment concentration distribution exhibits a local maximum around the null point.     

Stochastic particle based models for suspended particle movement in surface flows

Modeling of suspended sediment particle movement in surface water can be achieved by stochastic particle tracking model approaches.In this paper,different mathematical forms of particle tracking models are introduced to describe particle movement under various flow conditions,i.e.,the stochastic diffusion process,stochastic jump process,and stochastic jump diffusion process.While the stochastic diffusion process can be used to represent the stochastic movement of suspended particles in turbulent flows,the stochastic jump and the stochastic jump diffusion processes can be used to describe suspended particle movement in the occurrences of a sequence of extreme flows.An extreme flow herein is defined as a hydrologic flow event or a hydrodynamic flow phenomenon with a low probability of occurrence and a high impact on its ambient flow environment.In this paper,the suspended sediment particle is assumed to immediately follow the extreme flows in the jump process（i.e.the time lag between the flow particle and the sediment particle in extreme flows is considered negligible）.In the proposed particle tracking models,a random term mainly caused by fluid eddy motions is modeled as a Wiener process,while the random occurrences of a sequence of extreme flows can be modeled as a Poisson process.The frequency of occurrence of the extreme flows in the proposed particle tracking model can be explicitly accounted for by the Poisson process when evaluating particle movement.The ensemble mean and variance of particle trajectory can be obtained from the proposed stochastic models via simulations.The ensemble mean and variance of particle velocity are verified with available data.Applicability of the proposed stochastic particle tracking models for sediment transport modeling is also discussed.     

Study on the simulation of transparency of Lake Taihu under different hydrodynamic conditions

It was indicated in this study that there were negative relations between the concentrations of suspended solid (SS) and transparency according to the analysis of measured data of Lake Taihu. Their relations in pervious studies were reviewed, which showed that the changes of transparency in Lake Taihu could be reflected by simulating suspended solid concentration (SSC). Measured data showed that the changes of SSC with wind speed were similar at different water depths. SSC increased with the increasing of wind speed. Both wave and lake current of Lake Taihu had positive relations with SSC. However, wave was the main factor affecting sediment suspension, while flow took the second place. In this study, a numerical model coupling lake current, wave and SSC of Lake Taihu was developed. In the SS model, the combined effects of wave and current were included. The amounts of suspended and deposited sediments near the lake bed surface layer were treated separately. The stochastic characteristics of turbulent flow pulsation near lake beds were also considered, and the start-up conditions of sediment suspension were introduced to the model. The model elucidated the mutual exchange processes between sediment particles in SS and active sediments within and on the bed surface layer. Simulated results showed that lake current had relatively significant effects on the SSC at littoral areas of Lake Taihu, while SSC at the central area of the lake was mainly influenced by wave. The changes of transparency with SSC were simulated for Lake Taihu using this model. Calculated results were validated by measured data with good fitness, which indicated that the model is basically suitable for the simulation and prediction of transparency of Lake Taihu.     

Study on the simulation of transparency of Lake Taihu under different hydrodynamic conditions

It was indicated in this study that there were negative relations between the concentrations of suspended solid (SS) and transparency according to the analysis of measured data of Lake Taihu. Their relations in pervious studies were reviewed, which showed that the changes of transparency in Lake Taihu could be reflected by simulating suspended solid concentration (SSC). Measured data showed that the changes of SSC with wind speed were similar at different water depths. SSC increased with the increasing of wind speed. Both wave and lake current of Lake Taihu had positive relations with SSC. However, wave was the main factor affecting sediment suspension, while flow took the second place. In this study, a numerical model coupling lake current, wave and SSC of Lake Taihu was developed. In the SS model, the combined effects of wave and current were included. The amounts of suspended and deposited sediments near the lake bed surface layer were treated separately. The stochastic characteristics of turbulent flow pulsation near lake beds were also considered, and the start-up conditions of sediment suspension were introduced to the model. The model elucidated the mutual exchange processes between sediment particles in SS and active sediments within and on the bed surface layer. Simulated results showed that lake current had relatively significant effects on the SSC at littoral areas of Lake Taihu, while SSC at the central area of the lake was mainly influenced by wave. The changes of transparency with SSC were simulated for Lake Taihu using this model. Calculated results were validated by measured data with good fitness, which indicated that the model is basically suitable for the simulation and prediction of transparency of Lake Taihu.     

Turbulence and suspended sediment processes in the Garonne River tidal bore in November 2016

A tidal bore is a water discontinuity at the leading edge of a ood tide wave in estuaries with a large tidal range and funneling topography. New measurements were done in the Garonne River tidal bore on 14 15 November 2016, at a site previously investigated between 2010 and 2015. The data focused on long, continuous, high-frequency records of instantaneous velocity and suspended sediment con- centration (SSC) estimate for several hours during the late ebb, tidal bore passage and ood tide. The bore passage drastically modi ed the ow eld, with very intense turbulent and sediment mixing. This was evidenced with large and rapid uctuations of both velocity and Reynolds stress, as well as large SSCs during the ood tide. Granulometry data indicated larger grain sizes of suspended sediment in water samples compared to sediment bed material, with a broader distribution, shortly after the tidal bore. The tidal bore induced a sudden suspended sediment ux reversal and a large increase in suspended sedi- ment ux magnitude. The time-variations of turbulent velocity and suspended sediment properties indicated large uctuations throughout the entire data set. The ratio of integral time scales of SSC to velocity in the x-direction was on average TE,SSC/TE,x 0.16 during the late ebb tide, compared to TE,SSC/ TE,x 0.09 during the late ood tide. The results imply different time scales between turbulent velocities and suspended sediment concentrations.     

Sediment transport in ice-covered channels

The existence of ice cover has important effects on sediment transport and channel morphology for rivers in areas with an annual occurrence of an ice season. The interaction of sediment transport and s...     

Experiment and simulation of turbulent flow in local scour around a spur dyke

The turbulent flow in the local scour hole around a single non-submerged spur dyke is investigated with both experimental and numerical methods. The experiments are conducted under clear-water scour regime with an impermeable spur dyke. The scour geometry and flow velocities are measured in details with a high-resolution laser displacement meter, electro-magnetic velocimetries and PIV （Particle image velocimetry）. A 3D non-linear k-ε model is developed to simulate the complex local flow field around the scour area. The numerical model is formulated using FVM （Finite volume method） on a collocated unstructured mesh, capable of resolving complex geometries and boundaries. It is found that the simulation results are reasonably consistent with those of the experimental measurements. Based on the study results, the nature of the flow structure around a spur dyke with local scour hole is analyzed.