Pier scour within long contraction in cohesive sediment bed

In the present study, an attempt has been made to study the bridge pier scour embedded long contraction for the bed sediment prepared by a mixture of clay and fine sand with varying proportions having a specific range of antecedent moisture content. Results particularly focused on the clay–sand mixed cohesive bed at varying clay fractions of the sediment bed, approach flow velocity, contraction ratio and different pier shapes, on maximum equilibrium scour depth for pier scour within long contraction. Further, regression based equations for the estimation of non-dimensional maximum scour depth for piers within long contraction in clay–sand mixed cohesive bed embedded were proposed.     

Simulation-based optimization of in–stream structures design: bendway weirs

Bendway weirs are one of the most practical in–stream rock structures utilized to protect the outer bend of meandering streams and rivers from erosion. We present development of a simulation-based paradigm for effective design of bendway weir structures to enhance meandering stream bank stability and control lateral migration. To do so, we employ the St. Anthony Falls Laboratory Virtual StreamLab (VSL3D) code to elucidate the flow and sediment transport phenomena induced by interaction of flow, mobile bed, and in–stream structures in large rivers under prototype conditions. We carried out numerous numerical experiments to systematically simulate various arrangements of bendway weir in two river test-beds and gaining insights into the physical mechanisms via which such bendway weirs modify turbulent flow, sediment transport and scour processes. The so-gained physical insights are then taken into account to develop a set of practical physics-based design criteria for optimal placement of bendway weirs in large rivers.     

Three-dimensional analysis of coherent turbulent flow structure around a single circular bridge pier

The coherent turbulent flow around a single circular bridge pier and its effects on the bed scouring pattern is investigated in this study. The coherent turbulent flow and associated shear stresses play a major role in sediment entrainment from the bed particularly around a bridge pier where complex vortex structures exist. The conventional two-dimensional quadrant analysis of the bursting process is unable to define sediment entrainment, particularly where fully three-dimensional flow structures exist. In this paper, three-dimensional octant analysis was used to improve understanding of the role of bursting events in the process of particle entrainment. In this study, the three-dimensional velocity of flow was measured at 102 points near the bed of an open channel using an Acoustic Doppler Velocity meter (Micro-ADV). The pattern of bed scouring was measured during the experiment. The velocity data were analysed using the Markov process to investigate the sequential occurrence of bursting events and to determine the transition probability of the bursting events. The results showed that external sweep and internal ejection events were an effective mechanism for sediment entrainment around a single circular bridge pier. The results are useful in understanding scour patterns around bridge piers.     

Hydrodynamics and contaminant transport on a degraded bed at a 90-degree channel confluence

Channel confluences at which two channels merge have an important effect on momentum exchange and contaminant diffusion in both natural rivers and artificial canals. In this study, a three-dimensional numerical model, which is based on the Reynolds Averaged Navier–Stokes equations and Reynolds Stress Turbulence model, is applied to simulate and compare flow patterns and contaminant transport processes for different bed morphologies. The results clearly show that the distribution of contaminant concentrations is mainly controlled by the shear layer and two counter-rotating helical cells, which in turn are affected by the discharge ratio and the bed morphology. As the discharge ratio increases, the shear flow moves to the outer bank and the counter-clockwise tributary helical cell caused by flow deflection is enlarged, leading the mixing happens near the outer bank and the mixing layer distorted. The bed morphology can induce shrinkage of the separation zone and increase of the clockwise main channel helical cell, which is initiated by the interaction between the tributary helical cell and the main channel flow and strengthened by the deep scour hole. The bed morphology can also affect the distortion direction of the mixing layer. Both a large discharge ratio and the bed morphology could lead to an increase in mixing intensity.     

Estimation of maximum scour depths at upstream of front and rear piers for two in-line circular columns

Previous investigations indicate that scour around bridge piers is one of the most important factors for the failure of waterway bridges. Hence, it is essential to determine the accurate scour depth around the bridge piers. Most of the previous studies were based on scour around a single pier; however, in practice, new bridges are usually wide and then piers comprise two circular piers aligned in the flow direction that together support the loading of the structure. In this study, the effect on maximum scour depth of the spacing between two piers aligned in the flow direction was investigated experimentally under clear water scour conditions. The results show that the maximum scour depth at upstream of the front pier occurs when the spacing between the two piers is 2.5 times the diameter of the pier. Two semi empirical equations have been developed to predict the maximum scour depth at upstream of both front and rear piers as a function of the spacing between the piers, in terms of a pier-spacing factor. If the new equations for the pier-spacing factor are used with some of the existing equations for scour at a single pier, the predicted scouring depths are in good agreement with observed results. The S/M equation exhibited the best performance among the various equations tested and was recommended for use in prediction of the equilibrium scour depth. The findings of this study can be used to facilitate the positioning of piers when scouring is a design concern.     

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.     

Evaluation of existing equations for temporal scour depth around circular bridge piers

Scour is defined as the processes of removal of sediment particles from water stream bed by the erosive action of activated water, and also carries sediment away from the hydraulic structures. Scour is the main cause of pier failure. Numerous equations are available for estimating temporal and equilibrium scour depth. The present study describes the phenomenon of temporal scour depth variation at bridge piers and deals with the methods for its estimation. The accuracy of six temporal scour depth equations are also checked in this study. After graphical and statistical analysis, it was found that the relationship proposed by Oliveto and Hager (J Hydraul Eng (ASCE) 128(9):811–520, 2002) predicts temporal scour depth better than other equations. Three equations of equilibrium time of scour are also used for computing equilibrium time. Equilibrium time equation proposed by Choi and Choi (Water Environ J 30(1–2):14–21, 2016) gives better agreements with observed values.     

Adsorptive removal of Drimarine Red HF-3D dye from aqueous solution using low-cost agricultural waste: batch and column study

This study involves the utilisation of peanut husk for the removal of Drimarine Red HF-3D dye from aqueous solutions. Batch study experiments were conducted with native, HNO₃-treated and Na-alginate-immobilised peanut husk biomass. Maximum dye removal (95.24 mg/g) was obtained with HNO₃-treated biomass. The experimental data were successfully explained with a pseudo-second-order kinetic model for all types of biosorbents. The equilibrium data fitted well to the Freundlich adsorption isotherm model. A thermodynamic study was also carried out to check the nature of the adsorption process. A fixed-bed column study for Drimarine Red HF-3D was carried out to optimise the effect of bed height, flow rate and initial dye concentration using peanut husk biomass. The column study showed that biosorption capacity increased with the increase in initial dye concentration and bed height, but decreased with increased flow rate. Data for Drimarine Red HF-3D were in very good agreement with the bed depth service time model. Fourier transform infrared analysis demonstrated the involvement of different functional groups in dye biosorption. These results showed that peanut husk biomass possessed good potential for the removal of Drimarine Red HF-3D from aqueous solution.     

Steady Shallow-Water Current and Solute Transport around a Semi-Conical Headland

Laboratory experiments have been carried out to investigate the effects of a sloping wall headland on the flow characteristics and the associated concentration distributions from a point source around the headland. A semi-conical headland with a slope of 1:2 was set up in a flow basin, 4.8 m long and 3.8 m wide. In this paper, the experimental results of a steady shallow-water current are reported. Three dimensional flow velocities in the basin were measured using Sontek-ADV instrument. The dye concentration levels in the basin were measured by two fluorometers. The experimental results showed a large-scale re-circulation region behind the semi-conical headland. The peak turbulence energy, at about 53% of the local kinetic flow energy, coincides with the region of high velocity gradient. Significant vertical flows were observed around the area near the downhill slope of the headland, with a maximum ratio of vertical to horizontal velocities being about 22%. Such relatively significant vertical scouring velocities, coupled with strong turbulence energy and high horizontal velocity gradients in the same region, could cause severe bed erosion. The experimental results have also been compared with the predicted results of a depth-averaged numerical model. The predicted eddy structure and the concentration distribution in the re-circulation area were found to compare favourably with the experimental results. However, the discrepancies in the flow velocities and the concentration levels near the headland were apparent. It was observed that the dye concentration continued to spread in the cross-stream direction after passing the headland, whereas only a limited extent of the lateral spreading was predicted by the numerical model further downstream of the headland.     

Fixed-bed adsorption of Malachite Green onto binary solid mixture of adsorbents: seashells and eggshells

In this work, continuous removal of Malachite Green from aqueous solution by adsorption was investigated using a laboratory scale fixed-bed column packed with binary solid mixture of seashell powder (SSP) and eggshell powder (ESP). The effects of initial dye concentration, feed flow rate, bed height, and the SSP?:?ESP ratio in the binary solid mixture on the breakthrough characteristics were studied. The breakthrough time decreased with increase in the flow rate and the initial dye concentration. The breakthrough time also varied greatly with change in bed height. The breakthrough time increased from 17.2?h to 34.5?h with increase in bed height from 6?cm to 12?cm. The breakthrough time was significantly affected by the change in the SSP?:?ESP ratio in the binary mixture. The breakthrough time decreased with increase in the SSP ratio in the binary mixture. In order to determine the column kinetic parameters, three different column kinetic models, namely Adams–Bohart, Bed Depth Service Time (BDST), and Yoon–Nelson models were fitted to the dynamic flow experimental data. The BDST model showed good agreement with the experimental results at all the process parameters studied. Results suggest that this novel approach of simultaneous utilization of adsorbents can be adopted for the treatment of dye-bearing effluents.     

A new model for soil–water drainage problems

Seepage flow is an agent related to the transport and dispersion of contamination in groundwater. Steady two-dimensional seepage flow is governed by Laplace’s equation, for which several solution techniques are available. Because computations are complex from a practical point of view, simplified models encompass the Dupuit-Forchheimer approach assuming a horizontal flow. However this approach is inaccurate in seepage problems involving steep drawdowns. In this research, a new theoretical model for 2D seepage flow is proposed based on Fawer’s theory for curved flows Castro-Orgaz (Environ Fluid Mech 10(3):2971–2310, 2010), from which a second-order equation results describing the seepage surface. From this development, a numerical solution for the rectangular dam problem based on the second-order model is presented, whereas a simple first-order equation is found to describe flow to drains under a uniform rainfall. The results of this new model are compared with the full 2D solution of Laplace’s equation for typical test cases, resulting in an excellent agreement.     

Numerical simulation of scour around pipelines using an Euler–Euler coupled two-phase model

The scour around a long fixed pipeline placed just above a non-cohesive sandy bed is numerically simulated using an Eulerian two-phase model that implements Euler–Euler coupled governing equations for fluid and solid phases and a modified k−ɛ turbulence closure for the fluid phase, the modeling system being a part of the CFD software package FLUENT. Both flow–particle and particle–particle interactions are considered in the model. During the simulations, the interface between sand and water is specified using a threshold volume fraction of sand, and the evolution of the bedforms is studied in detail. The predictions of bedform evolution are in good agreement with previous laboratory measurements. Investigations into the mechanisms of scour reveal that three sediment transport modes (bed-load, suspended-load and laminated-load) are associated with the scour development. While some previously proposed scour development formulae for cylindrical objects are in good general agreement with the simulations, scour prediction based on a commonly used operational mine-burial model (DRAMBUIE) shows disparities with present simulations.     

2D numerical analysis of the influence of near-bank vegetation patches on the bed morphological adjustment

This paper investigates the influence of near-bank vegetation patches on the bed morphological adjustment in open channel flow systems. The 2D depth-averaged hydro-morphological model is adopted for this investigation, which is first validated by laboratory experimental data measured in an open channel with a single near-bank vegetation patch. The validated model is then applied for extensive numerical simulations, with the aim of conducting a systematic analysis of the influence of different geometric controlling parameters on the bed morphological evolution. The controlling parameters taken into account for numerical analysis include the angle of repose value (RAV) of sediment, vegetation density (VD), patch length (PL) and patch width (PW). The numerical results and analysis show that: (1) the RAV of sediment with slope-failure parametrization only influences the shape of the transverse bed topography in the junction region; (2) increase in VD, PL and PW that substantially enhances flow blockage effect encourages the growth of the pool adjacent to the patch in three dimensions; (3) increase in VD, PL and PW produces analogous retrogressive erosion (erosion toward the upstream) in the pool region, presumably due to the increase in flow resistance. Additional numerical experiments suggest that the staggered-order distribution of multiple patches might be an optimal choice for channel restoration and conservation since pools and riffles with larger scales can be produced.

     

Toxiwasp∗

TOXIWASP combines most of the kinetic structure of EXAMS 2 with the transport capabilities of WASP (Water Analysis Simulation Program). TOXIWASP uses variable chemical degradation rates from chemical properties and the environmental conditions of the aquatic ecosystem. These rates are reduced from pseudo first‐order rates to first‐order rates including the processes hydrolysis, biotransfor‐mation, phototransformation, oxidation, and volatilisation. Assuming ultimate local equilibrium, and using a chemical dependent partition coefficient as well as spatially varying environmental carbon fractions, sorption onto sediments and biomass is calculated. Environmental alternations could be specified in any time scale by providing monitoring data.

TOXIWASP generates total sediment and chemical concentrations every time step in every segment, including surface water, subsurface water, surface bed and subsurface bed. Advection, dispersion, mass loading, sedimentation, and scour affect sediment concentration in the water column and in the bed sediment concentrations depend on burial and erosion. In addition chemical concentrations are influenced by degradation, sediment‐water dispersion, and percolation. Lateral transport of chemical within the bed is neglected and transport data are not calculated in the program. TOXIWASP is developed to model stratified lakes, reservoirs, large rivers, estuaries, and coastal waters. As for EXAMS 2 (Burns et al.²) the TOXIWASP user has to accept the model's inability to connect the water body to a chemically contaminated atmosphere.     

RANS-based simulation of transverse turbulent mixing in a 2D geometry

Although transverse mixing is a significant process in river engineering when dealing with the discharge of pollutants from point sources or the mixing of tributary inflows, no theoretical basis exists for the prediction of its rate, which is indeed based upon the results of experimental works carried on in laboratory channels or in streams and rivers. The paper presents the preliminary results of a numerical study undertaken to simulate the transverse mixing of a steady-state point source of a tracer in a two-dimensional rectangular geometry, which is expected to reproduce a shallow flow. This geometry is that of Lau and Krishnappan (J Hydraul Div 13(HY10):1173–1189, 1977), who collected turbulent mixing data for a shallow flow. In the numerical study an approach based on the Reynolds Averaged Navier–Stokes (RANS) equations was applied, where the closure problem was solved by using turbulent viscosity concept. Particularly, the classical two-equations k–? model was used. Two methods were applied to the model results to evaluate the turbulent transverse mixing coefficient. The effect on transverse mixing of a grid located upstream the tracer source was also studied. Numerical results were generally higher than the experimental data. This overestimation could be explained considering the hypothesis of isotropic turbulence underlying the k–? model, which can lead to large turbulent viscosities and rate of mixing. However, RANS-based results may still be considered acceptable also providing the large uncertainties associated with literature predictive equations.     

Prediction of scour depth and dune morphology around circular bridge piers in seepage affected alluvial channels

Prediction of scour depth is one of the most significant problems in designing of bridges. Due to complexity of scour phenomenon, available empirical equations do not always offer accurate scour depth prediction. In this work, experiments were conducted on non-uniform sands using two piers of different diameter with and without seepage conditions. It has been found that, the scour depth at upstream of piers is decreased with application of downward seepage. The observed scour depths are compared with scour depth predicted by various prediction methods and it has been observed that the available equations are not suitable for seepage conditions. The present work thus introduces a new empirical equation for prediction of scour depth at piers with inclusion of seepage. The features of migrating dune like bedforms at downstream of piers due to deposition of scoured bed material are also explored. Height of deposition is found to be increased with downward seepage. The empirical equation describing morphology of dunes behind piers is also developed by incorporating downward seepage parameter.     

One-dimensional modelling of sediment deposition in reservoirs with sinking streams

It is common in karst regions that rivers are occasionally cut by mountains, resulting in the alternate appearances of open channel flow and pressurized flow. With more and more reservoirs being built in this region, the complicated sediment transport processes of such mixture flow are urgently needed to study. In this paper a one-dimensional numerical model with free surface and pressurized flows coupled together is presented. The simulated results are analyzed to explore deposition process in reservoirs with sinking streams; impacts of different hydraulic conditions on the sedimentation are also studied. To verify the computed results, a flume experiment is also conducted. The results show that deposition of sediments mainly occurs in open channel reaches where the longitudinal profile is similar to that of conventional reservoirs, i.e. typical delta has formed, demonstrating characteristics of deltaic deposition morphology in that the crest of delta moves to the downstream direction over time. The model provided by this paper is not only proved to simulate the characteristics of deposition in karst reservoirs successfully, but also reveal the impacts of hydraulic conditions in such circumstances.     

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

为探究和量化在水流冲刷下,排水管道中不同粒径颗粒物冲刷沉积的过程,本文模拟排水管道内沉积颗粒的冲刷过程.冲刷过程中,粒径较小的悬移质颗粒(小于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 of dam/levee-break hydrodynamics with a three-dimensional implicit unstructured-mesh finite element model

Dam failures usually cause huge economic and life losses , especially in urban areas where there is a high concentration of inhabitants and economic actors. In order to understand the physical mechanisms of the formation and development of dam-break flooding, lots of efforts have been put into different types of modelling techniques. However, most of existing models are 1D (one-dimensional) or 2D models based on the shallow water equations. In this paper, we present a 3D numerical modelling investigation of dam-break flow hydrodynamics in an open L-shape channel. A newly developed 3D unstructured mesh finite element model is used here. An absorption-like term is introduced to the Navier–Stokes equations in order to control the conditioning of the matrix equation in the numerical solution process and thus improve the stability. A wetting and drying algorithm is used here to allow the free surface height to be treated with a high level of implicitness and stability. The 3D model has been validated by comparing the results with the published experimental data. Good agreement has been achieved at six selected locations. This study shows that the 3D unstructured mesh model is capable of capturing the 3D hydraulic aspects and complicated local flows around structures in simulation of dam-break flows.