Efficient non-hydrostatic modelling of flow and bed shear stress in a pier scour hole

2014 ◽  
Vol 41 (5) ◽  
pp. 450-460 ◽  
Author(s):  
S. Pournazeri ◽  
S.S. Li ◽  
F. Haghighat

Predicting 3-D flow in a pier scour hole and the associated bed shear stress τb is important for the safe and economical design of bridge piers. This paper combines layered, hydrostatic hydrodynamic computations with non-hydrostatic pressure corrections, exploring a new modelling approach for efficient and reliable predictions of 3-D flow velocity. The law of the wall method is used for estimating τb. Its suitability for incorporation into layered models for bedload transport and pier scour simulations is also discussed. The predicted flow shows realistic features: strong downward flow adjacent to the upstream nose of a circular pier, vortex motions in the vertical and horizontal direction, and meandering flow wakes. The velocity results compare well with available experimental data. In the approach region, τb is uniform. It attains a local maximum immediately before flow enters the scour hole and then drops non-linearly in the scour-hole region toward the pier. In the wake region, τb has very low values. The τb predictions are consistent with the experimental data. In multi-layer models, when applying the law of wall method, one should use near-bed velocities as opposed to bottom-layer velocities to obtain more reliable τb estimates and avoid noisy results, which can cause a numerical instability problem in bedload transport simulations.

1989 ◽  
Vol 111 (2) ◽  
pp. 160-164 ◽  
Author(s):  
R. J. Kind ◽  
F. M. Yowakim ◽  
S. A. Sjolander

Expressions for the logarithmic portion of the law of the wall are derived for the axial and tangential velocity components of swirling flow in annular ducts. These expressions involve new shear-velocity scales and curvature terms. They are shown to agree well with experiment over a substantial portion of the flow near both walls of an annulus. The resultant velocity data also agree with the law of the wall. The success of the proposed logarithmic expressions implies that the mixing-length model used in deriving them correctly describes flow-velocity behavior. This model indicates that the velocity gradient at any height y in the near-wall region is determined by the wall shear stress, not by the local shear stress. This suggests that the influence of wall shear stress is dominant and that it determines the near-wall wall flow even in flows with curvature and pressure gradient. A physical explanation is suggested for this.


Author(s):  
Dipankar Biswas ◽  
Steven A. Lottes ◽  
Pradip Majumdar ◽  
Milivoje Kostic

Bridges are a significant component of the ground transportation infrastructure in the United States. With about sixty percent of bridge failures due to hydraulic causes, primarily scour, application of computational fluid dynamics (CFD) analysis techniques to the assessment of risk of bridge failure under flood conditions can provide increased accuracy in scour risk assessment at a relatively low cost. The analysis can be used to make optimum use of limited federal and state funds available to maintain and replace bridges and ensure public safety while traveling on the nation’s roads and highways during and after floods. Scour is the erosion of riverbed material during high flow conditions, such as floods. When scouring of the supporting soil around the piers and abutments of bridges takes place, risk of bridge failure increases. A simulation methodology to conservatively predict equilibrium shape and size of the scour hole under pressure flow conditions for flooded bridge decks using commercial CFD software was developed. The computational methodology has been developed using C++ to compute changes in the bed contour outside of the CFD software and generate a re-meshing script to change the bed boundary contour. STAR-CD was used to run the hydrodynamic analysis to obtain bed shear stress, and a BASH script was developed to automate cycling between computing bed shear stress with the CFD software and computing changes in the bed contour due to scour predicted using the computed shear stress for the current bed contour. A single-phase moving boundary formulation has been developed to compute the equilibrium scour hole contour that proceeds through a series of quasi-steady CFD computations. It is based on CFD analysis of the flow fields around the flooded bridge deck and shear stress computed at the bed modeled as a rough wall. A high Reynolds number k-ε turbulence model with standard wall functions, based on a Reynolds-Averaged Navier-Stokes (RANS) turbulence model, was used to compute bed shear stress. The scour sites on the bed were identified as those sites where the computed shear stress exceeded the critical shear stress computed from a published correlation for flat bed conditions. Comparison with experimental data obtained from the Turner-Fairbank Highway Research Center (TFHRC), McLean, VA, USA, revealed larger discrepancies than anticipated between the bridge inundation ratio and the scour hole depth. Although scour hole slopes were small for the cases tested, a correction to critical shear stress to account for bed slope was also tested. It did not significantly improve the correlation between CFD prediction and experimental observations. These results may be a consequence of using only excess shear stress above critical as a criteria for scour when other physical mechanisms also contribute to the initiation of scour. Prediction of scour depth using federal guidelines over predicts scour depth by as much as an order of magnitude in some cases. Over prediction is acceptable for purposes of ensuring bridge safety. CFD methods for scour prediction can be a significant improvement of current methods as long as under prediction of scour depth is avoided. Conservative scour prediction using CFD methods can be achieved by using conservative values of parameters such as critical shear stress and effective bed roughness.


2018 ◽  
Vol 66 (1) ◽  
pp. 23-31 ◽  
Author(s):  
Magdalena M. Mrokowska ◽  
Paweł M. Rowiński ◽  
Leszek Książek ◽  
Andrzej Strużyński ◽  
Maciej Wyrębek ◽  
...  

Abstract Two sets of triangular hydrographs were generated in a 12-m-long laboratory flume for two sets of initial bed conditions: intact and water-worked gravel bed. Flowrate ranging from 0.0013 m3 s-1 to 0.0456 m3 s-1, water level ranging from 0.02 m to 0.11 m, and cumulative mass of transported sediment ranging from 4.5 kg to 14.2 kg were measured. Then, bedload transport rate, water surface slope, bed shear stress, and stream power were evaluated. The results indicated the impact of initial bed conditions and flow unsteadiness on bedload transport rate and total sediment yield. Difference in ratio between the amount of supplied sediment and total sediment yield for tests with different initial conditions was observed. Bedload rate, bed shear stress, and stream power demonstrated clock-wise hysteretic relation with flowrate. The study revealed practical aspects of experimental design, performance, and data analysis. Water surface slope evaluation based on spatial water depth data was discussed. It was shown that for certain conditions stream power was more adequate for the analysis of sediment transport dynamics than the bed shear stress. The relations between bedload transport dynamics, and flow and sediment parameters obtained by dimensional and multiple regression analysis were presented.


2007 ◽  
Vol 34 (10) ◽  
pp. 1312-1323 ◽  
Author(s):  
Bahram Gharabaghi ◽  
Chris Inkratas ◽  
Spyros Beltaos ◽  
Bommanna Krishnappan

The Mackenzie River has several anomalous deep scour holes in a number of river channels in its delta. Proposed gas pipeline crossings have renewed interest in studying the stability of these scour holes. The main goal of this research project was to study flow velocity and bed shear stress distributions for a 30 m deep hole in the East Channel of the Mackenzie Delta as a first step toward assessing the stability of the scour hole and the risk of its migration during various flow conditions. In this study, a three-dimensional (3D) finite element flow model, FLUENT, using the renormalization group (RNG) k-ε turbulence model (where k is the turbulent kinetic energy and ε is the turbulence dissipation rate) was set up for the scour hole and calibrated using detailed measurements of 3D flow velocities, obtained with an acoustic doppler current profiler. The numerical model was then applied to predict flow velocity and bed shear stress distributions in and around the scour hole for three flow conditions (720, 1000, and 1400 m3/s). Results indicate that two vortices are formed in the river elbow above the scour hole. As the flow rate changed, the sizes of the vortices varied. The region upstream of the hole experienced the greatest magnitudes of bed shear stress.Key words: computational fluid dynamics, finite element, bed shear stress, deep hole, flow reversal.


2018 ◽  
Vol 40 ◽  
pp. 04012 ◽  
Author(s):  
Renske C. Terwisscha van Scheltinga ◽  
Heide Friedrich ◽  
Giovanni Coco

Traditional sediment transport equations calculate sediment flux from bed shear stress and the equations predict that transport increases nonlinearly with an increase in flow velocity. In a dune field, the dune geometry affects the flow velocity causing accelerating flow over the dune crest and de- and reattachment of the flow downstream of the dune crest. Sediment flux predicted from the reach-averaged bed shear stress gives fairly good results for dune fields, though their simplification is discordant for the complexity of the processes involved. Measurements of the displacement of sand particles over the dune bed were derived from highfrequency image capturing. The two main methods to measure particle velocities from images are particle tracking velocimetry (PTV) and particle image velocimetry (PIV). We compare individual particle tracking with a PIV-based correlation method. The PIV-based method promises to be a more efficient and effective approach to track particle motion. It is more suitable for the conditions of high bedload transport, as present in our experiments. The PIV-based method is based on using images of difference (IoD) and is fully automated and identifies spatial gradients at a support scale in the order of centimetres. Findings align with our general knowledge of accelerating flow over the dune crest. The mean streamwise particle velocity and activity over a dune stoss slope increase. At the scale of 0.026 m the observed particle velocity variability can be explained in the context of general onset and cessation of sediment transport, the effect of the reattachment zone and observed sweep/burst events. By decreasing the streamwise distance between cross-sections, the variations in mean particle velocity induced by superimposed bed defects are distinguished as well. The maximum particle velocity and activity occurred at the same location and consequently the location of the maximum transport over the dune crest was identified. The measurements bridge the gap between individual particle motion studies and (non-local) sediment transport flux measurements.


2015 ◽  
Vol 42 (6) ◽  
pp. 367-376
Author(s):  
S. Attar ◽  
S. Samuel Li

River dunes are important bedforms. Problems associated with the development and evolution of dune bedforms include increased flood risks, channel erosion, and damage to fish habitats. The purpose of this paper is to investigate the near-bed flow structure and bedload transport along gavel dune-like bedforms. The velocity field is computed using a relatively simple multi-layer hydrodynamic model, with a parameterization of flow separation in the leeside of dunes. The computation is of high efficiency and avoids uncertainties caused by flow separation. Fractional transport rates for a sediment mixture of sands and gravel are calculated using surface-based techniques. The computed flow velocities and bed shear stresses are in good comparison with acoustic Doppler velocimeter measurements. Bedload transport is shown to increase non-linearly with distance toward the dune crest and reach the maximum at the crest. This implies that dune-length averaged bed shear stress is not suitable for bedload calculations. At low discharges, the bed shear stress is the limiting factor, resulting in insignificant bedload. At high discharges when the bed shear stress exceeds a threshold, the effect of sediment-grain hiding and sediment-size availability are important for bedload calculations. The discharge–transport relationship is highly non-linear. This paper has demonstrated selective transport and potential dune surface coarsening. The simplified modelling approach has a good potential for application to field conditions.


2016 ◽  
Vol 11 (6) ◽  
pp. 1073-1081 ◽  
Author(s):  
Robin K. Biswas ◽  
◽  
Atsuhiro Yorozuya ◽  
Shinji Egashira ◽  
◽  
...  

A method is proposed to predict bank erosion and sand bar migration in river reaches where suspended sediment transport is dominant. The method focuses on the influence of the lateral bed slope on the erosion and deposition rate of suspended sediment, as well as on the profile of lateral bedload transport, assuming that geometric similarity holds in the bank region. In the proposed model, the erosion and deposition rate can be evaluated using either the bed shear stress at a reference location or the average bed shear stress in the bank region. In order to simulate bank erosion and associated bank shifting with a depth-integrated-base treatment, stretchable grids were added to the conventional coarse grid system near the bank. The proposed method, including the bank erosion model, is applied to the lower reach of the Brahmaputra River, which is ∼90 km long and ∼12.50 km wide. The computed results on bank shifting, sand bar migration, and sediment transport rates are compared with data obtained from field investigations and remote sensing. These results suggest that the proposed method is applicable for predicting sediment issues in river reaches dominated by suspended sediment.


Sign in / Sign up

Export Citation Format

Share Document