Advanced Numerical Solver for Dam-Break Flow Application

In this paper, a HLL (Harten Lax van Leer) approximate Riemann solver with MUSCL scheme (Monotonic Upwind Schemes for Conservative Laws) is implemented in the presented FV (Finite Volume) model. The presented model is used to simulate different dam-break flow events to verify its capability. Four test cases are presented in this paper. In the first test case, a 1-Dimensional (1D) dambreak flow is simulated over a rectangular channel with different slope limiters of the FV model (namely Godunov, Superbee, Minmod, van Leer, and van Albada). The second test case consists of a simulation of shallow water discontinuous dam-break flow over a dry-downstream bed channel. The third test simulates the shallow water dam-break flow with the existence of bed slope and bed shear stress. Finally, in the last test, the HLL-MUSCL model used in this paper and some other solver models used in literature are compared against the referred exact solution in dam-break flow application. The presented HLL-MUSCL scheme is found to give the best agreement to the exact solution

Numerical and experimental turbulence studies on shallow open channel flows

YesBased on the previous studies, the shallow water equations (SWEs) model was proven to be insufficient to consider the flow turbulence due to its simplified Reynolds-averaged form. In this study, the k-ε model was used to improve the ability of the SWEs model to capture the flow turbulence. In terms of the numerical source terms modelling, the combined k-ε SWEs model was improved by a recently proposed surface gradient upwind method (SGUM) to facilitate the extra turbulent kinetic energy (TKE) source terms in the simulation. The laboratory experiments on both the smooth and rough bed flows were also conducted under the uniform and non-uniform flow conditions for the validation of the proposed numerical model. The numerical simulations were compared to the measured data in the flow velocity, TKE and power spectrum. In the power spectrum comparisons, a well-studied Kolmogorov’s rule was also employed to complement both the numerical and experimental results and to demonstrate that the energy cascade trend was well-held by the investigated flows.The Major State Basic Research Development Program (973 program) of China (Grant Number 2013CB036402). Open Fund from the State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, China (Grant Number SKLH-OF-1103)

Shallow sediment transport flow computation using time-varying sediment adaptation length

YesBased on the common approach, the adaptation length in sediment transport is normally estimated in the temporal independence. However, this approach might not be theoretically justified as the process of reaching of the sediment transport equilibrium stage is affected by the flow conditions in time, especially for those fast sediment moving flows, such as scour-hole developing flow. In this study, the 2D shallow water formulation together with a sediment continuity-concentration (SCC) model were applied to flow with mobile sediment boundary. A time-varying approach was proposed to determine the sediment transport adaptation length to treat the flow sediment erosion-deposition rate. The proposed computational model was based on the Finite Volume (FV) method. The Monotone Upwind Scheme of Conservative Laws (MUSCL)-Hancock scheme was used with the Harten Lax van Leer-contact (HLLC) approximate Riemann solver to discretize the FV model. In the flow applications of this paper, a highly discontinuous dam-break fast sediment transport flow was used to calibrate the proposed time-varying sediment adaptation length model. Then the calibrated model was further applied to two separate experimental sediment transport flow applications documented in literature, i.e. a highly concentrated sediment transport flow in a wide alluvial channel and a sediment aggradation flow. Good agreements with the experimental data were presented by the proposed model simulations. The tests prove that the proposed model, which was calibrated by the discontinuous dam-break bed scouring flow, also performed well to represent the rapid bed change and the steady sediment mobility conditions.The National Natural Science Foundation of China NSFC (Grant Number 20101311246), Major State Basic Research Development Program (973 program) of China (Grant Number 2013CB036402) and Open Fund of the State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University of China (Grant Number SKLH-OF-1103)

Dominant features in three-dimensional turbulence structure: comparison of non-uniform accelerating and decelerating flows

The results are presented from an experimental study to investigate three-dimensional turbulence structure profiles, including turbulence intensity and Reynolds stress, of different non-uniform open channel flows over smooth bed in subcritical flow regime. In the analysis, the uniform flow profiles have been used to compare with those of the non-uniform flows to investigate their time-averaged spatial flow turbulence structure characteristics. The measured non-uniform velocity profiles are used to verify the von Karman constant κ and to estimate sets of log-law integration constant B r and wake parameter П, where their findings are also compared with values from previous studies. From κ, B r and П findings, it has been found that the log-wake law can sufficiently represent the non-uniform flow in its non-modified form, and all κ, B r and П follow universal rules for different bed roughness conditions. The non-uniform flow experiments also show that both the turbulence intensity and Reynolds stress are governed well by exponential pressure gradient parameter β equations. Their exponential constants are described by quadratic functions in the investigated β range. Through this experimental study, it has been observed that the decelerating flow shows higher empirical constants, in both the turbulence intensity and Reynolds stress compared to the accelerating flow. The decelerating flow also has stronger dominance to determine the flow non-uniformity, because it presents higher Reynolds stress profile than uniform flow, whereas the accelerating flow does not

Efficient numerical computation and experimental study of temporally long equilibrium scour development around abutment

YesFor the abutment bed scour to reach its equilibrium state, a long flow time is needed. Hence, the employment of usual strategy of simulating such scouring event using the 3D numerical model is very time consuming and less practical. In order to develop an applicable model to consider temporally long abutment scouring process, this study modifies the common approach of 2D shallow water equations (SWEs) model to account for the sediment transport and turbulence, and provides a realistic approach to simulate the long scouring process to reach the full scour equilibrium. Due to the high demand of the 2D SWEs numerical scheme performance to simulate the abutment bed scouring, a recently proposed surface gradient upwind method (SGUM) was also used to improve the simulation of the numerical source terms. The abutment scour experiments of this study were conducted using the facility of Hydraulics Laboratory at Nanyang Technological University, Singapore to compare with the presented 2D SGUM-SWEs model. Fifteen experiments were conducted over a total period of 3059.7 hours experimental time (over 4.2 months). The comparison shows that the 2D SGUM-SWEs model gives good representation to the experimental results with the practical advantage

Numerical modeling of the Sakuma Dam reservoir sedimentation

YesThe present study attempts to predict the reservoir sedimentation in 32 km region of the Tenryu River between the Hiraoka and Sakuma Dams in Japan. For numerical simulations of the reservoir sedimentation, the one-dimensional model of the Hydrologic Engineering Centre-River Analysis System (HEC-RAS) is used together with the inclusion of channel geometry, bed gradation curve, Exner-5 bed sorting mechanisms, fall velocity of the particle, and flow and sediment boundary conditions pertaining to modeling region. The modeling region of the Tenryu River is divided into 48 river stations with 47 reaches in the numerical simulations. The numerical model is calibrated using the available data for 48 years from 1957 to 2004. The formulae of sediment transport function, Manning’s roughness coefficient, computational increment and fall velocity have been identified for getting the best estimation of the Sakuma Dam reservoir sedimentation. Combination of obtained sensitive parameters and erodible limits of 2 m gave the best comparison with the measured bed profile. The computed results follow the trend of measured data with a small underestimation. Although Manning’s roughness coefficient has an effect on the sedimentation, no direct relation is found between the Manning’s roughness coefficient and reservoir sedimentation. It is found that the temperature of water has no effect on the reservoir sedimentation