Modelling of three-dimensional flow velocities in a deep hole in the East Channel of the Mackenzie Delta, Northwest Territories

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.

2009 ◽  
Vol 36 (5) ◽  
pp. 791-800 ◽  
Author(s):  
Chris Inkratas ◽  
Bahram Gharabaghi ◽  
Spyros Beltaos ◽  
Bommanna Krishnappan

The discovery of vast quantities of hydrocarbons during exploration of the Mackenzie Delta has led to the detection of numerous anomalous deep holes in various delta channels in the past four decades. The proposed pipeline crossing from the Mackenzie Delta to northern Alberta has renewed interest in studying the stability of these deep holes. The main goal of this research project was to study ice-covered flow velocity and bed shear stress distributions using a previously calibrated computational fluid dynamics (CFD) model for a 30 m deep hole in the East Channel of the Mackenzie Delta to provide an initial assessment of its stability. In this study, a previously developed 3D CFD model using the FLUENT code with the renormalization group theory (RNG) k-ε turbulence closure model with nonequilibrium wall functions was adapted to represent ice-covered flow conditions for the study reach. The numerical model was applied to simulate four flow conditions (500, 720, and 1000 m3/s along with a moving ice cover during a flow rate of 720 m3/s). Results indicated that a single vortex was formed near the inner bend above the hole and deposition may be a key process occurring during ice-covered flow conditions.


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.


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.


2021 ◽  
Vol 11 (15) ◽  
pp. 6898
Author(s):  
Ahmed M. Helmi ◽  
Ahmed H. Shehata

In the current study, a three-dimensional CFD model is utilized to investigate the variation of the flow structure and bed shear stress at a single cylindrical pier during scour development. The scour development is presented by seven solidified geometries of the scour hole, collected during previous experimental work at different scour stages. Different turbulence models are evaluated and the (k-ω) model is chosen due to its relative accuracy in capturing the flow oscillation and vortex shedding at the pier downstream side with personal computer computational and storage resources. The numerical results are verified against dimensionless parameters from different previous experimental works. This research describes in detail the flow structure and bed shear stress variations through seven stages of the scour hole development. The dimensionless area-averaged circulation coefficient (Ψi) is developed to evaluate the changes in the vortex strength through the scouring process by eliminating the calculation area effect. It was concluded that the circulation in the (Y) direction is the main driving factor in the development of the scour hole more than the circulation in the (X) direction. The ratio between the horseshoe vortex (HV) mean size and the scouring depth (DV/dS) in addition to the location of the maximum bed shear stress are investigated during different stages of the scour development.


2007 ◽  
Vol 19 (02) ◽  
pp. 117-127 ◽  
Author(s):  
Yang-Yao Niu ◽  
Ding-Yu Chang

In this work, a preliminary numerical simulation of the lower urinary system using Computational Fluid Dynamics (CFD) is performed. Very few studies have been done on the simulation of three-dimensional urine through the lower urinary system. In this study, a simplified lower urinary model with rigid body assumption is proposed. The distributions of urine flow velocity, wall pressure and shear stress along the urethra are simulated based on MRI scanned uroflowmetry of a normal female. Numerical results show that violent secondary flows appear on the cross surface near the end of the urethra when the inflow rate is increased. The oscillative variation of pressure and shear stress distributions are found around the beginning section of the urethra when flow rate is at the peak value.


PLoS ONE ◽  
2021 ◽  
Vol 16 (1) ◽  
pp. e0245245
Author(s):  
Yun-Feng Liu ◽  
Ke Gu ◽  
Yi-Ming Shu ◽  
Xian-Lei Zhang ◽  
Xin-Xin Liu ◽  
...  

As a type of flexible impermeable material, a PVC geomembrane must be cooperatively used with cushion materials. The contact interface between a PVC geomembrane and cushion easily loses stability. In this present paper, we analyzed the shear models and parameters of the interface to study the stability. Two different cushion materials were used: the common extrusion sidewall and non-fines concrete. To simulate real working conditions, flexible silicone cushions were added under the loading plates to simulate hydraulic pressure loading, and the loading effect of flexible silicone cushions was demonstrated by measuring the actual contact areas under different normal pressures between the geomembrane and cushion using the thin-film pressure sensor. According to elastomer shear stress, there are two main types of shear stress between the PVC geomembrane and the cushion: viscous shear stress and hysteresis shear stress. The viscous shear stress between the geomembrane and the cement grout was measured using a dry, smooth concrete sample, then the precise formula parameters of the viscous shear stress and viscous friction coefficient were obtained. The hysteresis shear stress between the geomembrane and the cushion was calculated by subtracting the viscous shear stress from the total shear stress. The formula parameters of the hysteresis shear stress and hysteresis friction coefficient were calculated. The three-dimensional box-counting dimensions of the cushion surface were calculated, and the formula parameters of the hysteresis friction were positively correlated with the three-dimensional box dimensions.


Water ◽  
2020 ◽  
Vol 12 (5) ◽  
pp. 1423
Author(s):  
Amir Golpira ◽  
Fengbin Huang ◽  
Abul B.M. Baki

This study experimentally investigated the effect of boulder spacing and boulder submergence ratio on the near-bed shear stress in a single array of boulders in a gravel bed open channel flume. An acoustic Doppler velocimeter (ADV) was used to measure the instantaneous three-dimensional velocity components. Four methods of estimating near-bed shear stress were compared. The results suggested a significant effect of boulder spacing and boulder submergence ratio on the near-bed shear stress estimations and their spatial distributions. It was found that at unsubmerged condition, the turbulent kinetic energy (TKE) and modified TKE methods can be used interchangeably to estimate the near-bed shear stress. At both submerged and unsubmerged conditions, the Reynolds method performed differently from the other point-methods. Moreover, a quadrant analysis was performed to examine the turbulent events and their contribution to the near-bed Reynolds shear stress with the effect of boulder spacing. Generally, the burst events (ejections and sweeps) were reduced in the presence of boulders. This study may improve the understanding of the effect of the boulder spacing and boulder submergence ratio on the near-bed shear stress estimations of stream restoration practices.


Water ◽  
2020 ◽  
Vol 12 (4) ◽  
pp. 964 ◽  
Author(s):  
Shampa ◽  
Yuji Hasegawa ◽  
Hajime Nakagawa ◽  
Hiroshi Takebayashi ◽  
Kenji Kawaike

This paper focuses on finding efficient solutions for the design of a highly permeable pile spur (or slit type) dike field used in morphologically dynamic alluvial rivers. To test the suitability of different arrangements of this type of permeable pile spur dike field, laboratory experiments were conducted, and a three-dimensional multiphase numerical model was developed and applied, based on the experimental conditions. Three different angles to the approach flow and two types of individual pile position arrangements were tested. The results show that by using a series of slit-type spurs, the approach velocity of the flow can be considerably reduced within the spur dike zone. Using different sets of angles and installation positions, this type of permeable spur dike can be used more efficiently than traditional dikes. Notably, this type of spur dike can reduce the longitudinal velocity, turbulence intensity, and bed shear stress in the near-bank area. Additionally, the deflection of the permeable spur produces more transverse flow to the opposite bank. Arranging the piles in staggered grid positions among different spurs in a spur dike field improves functionality in terms of creating a quasi-uniform turbulence zone while simultaneously reducing the bed shear stress. Finally, the efficacy of the slit-type permeable spur dike field as a solution to the riverbank erosion problem is numerically tested in a reach of a braided river, the Brahmaputra–Jamuna River, and a comparison is made with a conventional spur dike field. The results indicate that the proposed structure ensures the smooth passing of flow compared with that for the conventional impermeable spur structure by producing a lower level of scouring (low bed shear stress) and flow intensification.


Water ◽  
2020 ◽  
Vol 12 (10) ◽  
pp. 2753
Author(s):  
Liyuan Zhang ◽  
Faxing Zhang ◽  
Ailing Cai ◽  
Zhaoming Song ◽  
Shilin Tong

Bed shear stress is closely related to sediment transport in rivers. Bed shear stress estimation is very difficult, especially for complex flow fields. In this study, complex flow field measurement experiments in a 60° bend with a groyne were performed. The feasibility and reliability of bed shear stress estimations using the log-law method in a complex flow field were analyzed and compared with those associated with the Reynolds, Turbulent Kinetic Energy (TKE), and TKE-w′ methods. The results show that the TKE, Reynolds, and log-law methods produced similar bed shear stress estimates, while the TKE-w′ method produced larger estimates than the other methods. The TKE-w′ method was found to be more suitable for bed shear stress estimation than the TKE method, but the value of its constant C2 needed to be re-estimated. In a complex, strong, three-dimensional flow field, the height of the measurement point (relative or absolute) should be re-estimated when a single point measurement is used to estimate the bed shear stress. The results of this study provide guidance for experimental measurement of bed shear stress in a complex flow field.


Water ◽  
2019 ◽  
Vol 11 (11) ◽  
pp. 2346 ◽  
Author(s):  
Kiraga ◽  
Popek

Numerous approaches in sediment mobility studies highlighted the key meaning of channel roughness, which results not only from bed material granulation but also from various bed forms presence, caused by continuous sediment transport. Those forms are strictly connected with the intensity of particle transport, and they eventuate from bed shear stress. The present paper comprised of local scours geometric dimensions research in three variants of lengthwise development of laboratory flume in various hydraulic properties, both in “clear-water” and “live-bed” conditions of sediment movement. Lots of measurements of the bed conformation were executed using the LiDAR device, marked by a very precise three-dimensional shape description. The influence of the bed shear stress downstream model on scours hole dimensions of water structure was investigated as one of the key factors that impact the sediment transport intensity. A significant database of 39 experimental series, lasting averagely 8 hours, was a foundation for delineating functional correlations between bed shear stress-and-critical shear stress ratio and geometry properties of local scours in various flume development cases. In the scope of mutual influence of bed shear stress and water depth, high correlation coefficients were attained, indicating very good and good functional correlations. Also, the influence of bed shear stress and the total length of the scour demonstrated a high correlation coefficient.


Sign in / Sign up

Export Citation Format

Share Document