Numerical modelling of cohesive sediment transport in rivers

2004 ◽  
Vol 31 (5) ◽  
pp. 749-758 ◽  
Author(s):  
David H Willis ◽  
B G Krishnappan

Techniques available to practicing civil engineers for numerically modelling cohesive mud in rivers and estuaries are reviewed. Coupled models, treating water and sediment as a single process, remain research tools but are usually not three-dimensional. The decoupled approach, which separates water and sediment computations at each model time step, allows the three-dimensional representation of at least the bed and the use of well-proven, commercial, numerical, hydrodynamic models. Most hydrodynamic models compute sediment transport in suspension but may require modification of the dispersion coefficients to account for the presence of sediment. The sediment model deals with the sediment exchange between the water column and the bed using existing equations for erosion and deposition. Both equations relate the sediment exchange rates to the shear stress in the bottom boundary layer. In real rivers and estuaries, a depositional bed layer is associated with a period of low flow and shear, at slack tide for example, whereas in numerical models a layer is defined by the model time step. The sediment model keeps track of the uppermost layers at each model grid point, including consolidation and strengthening. Although numerical hydrodynamic models are based strongly on physics, sediment models are only numerical frameworks for interpolating and extrapolating full-scale field or laboratory measurements of "hydraulic sediment parameters," such as threshold shear stresses. Calibration and verification of models against measurement are therefore of prime importance.Key words: cohesive sediment, mathematical modelling, settling velocity, erosion, resuspension, deposition, fluid mud, bed layers.

2008 ◽  
Vol 31 (7) ◽  
pp. 1004-1014 ◽  
Author(s):  
Xiaobo Chao ◽  
Yafei Jia ◽  
F. Douglas Shields ◽  
Sam S.Y. Wang ◽  
Charles M. Cooper

Author(s):  
Phani Ganesh Elapolu ◽  
Pradip Majumdar ◽  
Steven A. Lottes ◽  
Milivoje Kostic

One of the major concerns affecting the safety of bridges with foundation supports in river-beds is the scouring of river-bed material from bridge supports during floods. Scour is the engineering term for the erosion caused by water around bridge elements such as piers, monopiles, or abutments. Scour holes around a monopile can jeopardize the stability of the whole structure and will require deeper piling or local armoring of the river-bed. About 500,000 bridges in the National Bridge Registry are over waterways. Many of these are considered as vulnerable to scour, about five percent are classified as scour critical, and over the last 30 years bridge failures caused by foundation scour have averaged about one every two weeks. Therefore it is of great importance to predict the correct scour development for a given bridge and flood conditions. Apart from saving time and money, integrity of bridges are important in ensuring public safety. Recent advances in computing boundary motion in combination with mesh morphing to maintain mesh quality in computational fluid dynamic analysis can be applied to predict the scour hole development, analyze the local scour phenomenon, and predict the scour hole shape and size around a pier. The main objective of the present study was to develop and implement a three dimensional iterative procedure to predict the scour hole formation around a cylindrical pier using the mesh morphing capabilities in the STARCCM+ commercial CFD code. A computational methodology has been developed using Python and Java Macros and implemented using a Bash script on a LINUX high performance computer cluster. An implicit unsteady approach was used to obtain the bed shear stresses. The mesh was iteratively deformed towards the equilibrium scour position based on the excess shear stress above the critical shear stress (supercritical shear stress). The model solves the flow field using Reynolds Averaged Navier-Stokes (RANS) equations, and the standard k–ε turbulence model. The iterative process involves stretching (morphing) a meshed domain after every time step, away from the bottom where scouring flow parameters are supercritical, and remeshing the relevant computational domain after a certain number of time steps when the morphed mesh compromises the stability of further simulation. The simulation model was validated by comparing results with limited experimental data available in the literature.


2017 ◽  
Vol 34 (7) ◽  
pp. 1469-1482 ◽  
Author(s):  
Daosheng Wang ◽  
Jicai Zhang ◽  
Ya Ping Wang ◽  
Xianqing Lv ◽  
Yang Yang ◽  
...  

AbstractThe model parameters in the suspended cohesive sediment transport model are quite important for the accurate simulation of suspended sediment concentrations (SSCs). Based on a three-dimensional cohesive sediment transport model and its adjoint model, the in situ observed SSCs at four stations are assimilated to simulate the SSCs and to estimate the parameters in Hangzhou Bay in China. Numerical experimental results show that the adjoint method can efficiently improve the simulation results, which can benefit the prediction of SSCs. The time series of the modeled SSCs present a clear semidiurnal variation, in which the maximal SSCs occur during the flood tide and near the high water level due to the large current speeds. Sensitivity experiments prove that the estimated results of the settling velocity and resuspension rate, especially the temporal variations, are robust to the model settings. The temporal variations of the estimated settling velocity are negatively correlated with the tidal elevation. The main reason is that the mean size of the suspended sediments can be reduced during the flood tide, which consequently decreases the settling velocity according to Stokes’s law, and it is opposite in the ebb tide. The temporal variations of the estimated resuspension rate and the current speeds have a significantly positive correlation, which accords with the dynamics of the resuspension rate. The temporal variations of the settling velocity and resuspension rate are reasonable from the viewpoint of physics, indicating the adjoint method can be an effective tool for estimating the parameters in the sediment transport models.


2005 ◽  
Vol 127 (4) ◽  
pp. 637-647 ◽  
Author(s):  
Yong He ◽  
Nandini Duraiswamy ◽  
Andreas O. Frank ◽  
James E. Moore

Background: Restenosis after stent implantation varies with stent design. Alterations in secondary flow patterns and wall shear stress (WSS) can modulate intimal hyperplasia via their effects on platelet and inflammatory cell transport toward the wall, as well as direct effects on the endothelium. Method of Approach: Detailed flow characteristics were compared by estimating the WSS in the near-strut region of realistic stent designs using three-dimensional computational fluid dynamics (CFD), under pulsatile high and low flow conditions. The stent geometry employed was characterized by three geometric parameters (axial strut pitch, strut amplitude, and radius of curvature), and by the presence or lack of the longitudinal connector. Results: Stagnation regions were localized around stent struts. The regions of low WSS are larger distal to the strut. Under low flow conditions, the percentage restoration of mean axial WSS between struts was lower than that for the high flow by 10–12%. The largest mean transverse shear stresses were 30–50% of the largest mean axial shear stresses. The percentage restoration in WSS in the models without the longitudinal connector was as much as 11% larger than with the connector. The mean axial WSS restoration between the struts was larger for the stent model with larger interstrut spacing. Conclusion: The results indicate that stent design is crucial in determining the fluid mechanical environment in an artery. The sensitivity of flow characteristics to strut configuration could be partially responsible for the dependence of restenosis on stent design. From a fluid dynamics point of view, interstrut spacing should be larger in order to restore the disturbed flow; struts should be oriented to the flow direction in order to reduce the area of flow recirculation. Longitudinal connectors should be used only as necessary, and should be parallel to the axis. These results could guide future stent designs toward reducing restenosis.


1984 ◽  
Vol 1 (19) ◽  
pp. 199
Author(s):  
E.J. Hayter ◽  
A.J. Mehta

Cohesive sediment related problems in estuaries include shoaling in navigable waterways and water pollution. A two-dimensional, depth averaged, finite element cohesive sediment transport model, CSTM-H, has been developed and may be used to assist in predicting the fate of sorbed pollutants and the frequency and quantity of dredging required to maintain navigable depths. Algorithms which describe the transport processes of redispersion, resuspenslon, dispersive transport, settling, deposition, bed formation and bed consolidation are incorporated in CSTM-H. The Galerkin weighted residual method is used to solve the advection-dispersion equation with appropriate source/sink terms at each time step for the nodal suspended sediment concentrations. The model yields stable and converging solutions. Verification was carried out against a series of erosion-deposition experiments in the laboratory using kaolinite and a natural mud as sediment. A model application under prototype conditions is described.


Sign in / Sign up

Export Citation Format

Share Document