Numerical Simulation of Bed Load and Suspended Load Sediment Transport Using Well-Balanced Numerical Schemes

Sediment transport can be modelled using hydrodynamic models based on shallow water equations coupled with the sediment concentration conservation equation and the bed con-servation equation.The complete system of equations is made up of the energy balance law and the Exner equations.The numerical solution for this complete system is done in a seg-regated manner.First,the hyperbolic part of the system of balance laws is solved using a finite volume scheme.Three ways to compute the numerical flux have been considered,the Q-scheme of van Leer,the HLLCS approximate Riemann solver,and the last one takes into account the presence of non-conservative products in the model.The discretisation of the source terms is carried out according to the numerical flux chosen.In the second stage,the bed conservation equation is solved by using the approximation computed for the system of balance laws.The numerical schemes have been validated making comparisons between the obtained numerical results and the experimental data for some physical experiments.The numerical results show a good agreement with the experimental data.

The Bedload Movement in the Changjiang Estuary

Sandwaves in the Changjiang estuary were measured with a shallow sediment profiler and an echosounder from 1978 to 1988. The data, together with grain size and bedform of sediment indicates that the bedload movement by rolling and saltation is of great significance to sediment transport and is the principal factor responsible for sandwave and sandbody development in the estuary. The sandwaves were found well-developed, which is related to the tidal range and the velocity of ebb current. However, the further growth is restricted by strong flood current prevailing in the estuary. Because of the significant bedload, the sandbodies shift obviously and frequently, and sometimes the exchange of position occurs between the sandbodies and tidal channels. As a result, ships are regularly forced to change their navigation course.

Two-Dimensional Mathematical Model of Tidal Current and Sediment for Oujiang Estuary and Wenzhou Bay

A 2-D mathematical model of tidal current and sediment has been developed for the Oujiang Estuary and the Wenzhou Bay. This model accomodates complicated features including multiple islands, existence of turbidity, and significant differ-ence in size distribution of bed material. The governing equations for non-uniform suspended load and bed load transport are presented in a boundary-fitted orthogonal curvilinear coordinate system. The numerical solution procedures along with their initial conditions, boundary conditions, and movable boundary technique are presented. Strategies for computation of the critical condition of deposition or erosion, sediment transport capacity, non-uniform bed load discharge, etc. are suggested. The model verification computation shows that, the tidal levels computed from the model are in good agreement with the field data at the 18 tidal gauge stations. The computed velocities and flow directions also agree well with the values measured along the totally 52 synchronously observed verticals distributed over 8 cross sections. The coraputed tidal water throughputs through the Huangda'ao cross section are close to the measured data. And the computed values of bed deformation from Yangfushan to the estuary outfall and in the outer-sea area are in good agreement with the data observed from 1986 to 1992. The changes of tidal volumes through the estuary, velocities in different channels and the bed form due to the influence of the reclamation project on the Wenzhou shoal are predicted by means of this model.

A numerical simulation study on mechanical behaviour of coal with bedding planes under coupled static and dynamic load

To investigate the bedding influence on coal mechanical behaviour in underground environments such as coal or rock burst, simulations of dynamic SHPB tests of pre-stressed coal specimens with different bedding angles were carried out using a particle flow code 2-dimensional(PFC2D). Three impact velocities of 4, 8 and 12 m/s were selected to study dynamic behaviours of coal containing bedding planes under different dynamic loads. The simulation results showed that the existence of bedding planes leads to the degradation of the mechanical properties and their weakening effect significantly depends on the angle h between the bedding planes and load direction. With h increaseing from 0° to 90°, the strength first decreased and subsequently increased and specimens became most vulnerable when h was 30° or 45°.Five failure modes were observed in the specimens in the context of macro-cracks. Furthermore, energy characteristics combined with ultimate failure patterns revealed that maximum accumulated energy and failure intensity have a positive relation with the strength of specimen. When bedding planes were parallel or perpendicular to loading direction, specimens absorbed more energy and experienced more violent failure with increased number of cracks. In contrast, bedding planes with h of 30° or 45° reduced the specimens' ability of storing strain energy to the lowest with fewer cracks observed after failure.

Experiment and simulation on the pyrite removal from the recirculating load of pulverizer with a dilute phase gas-solid fluidized bed

In order to reduce the energy consumption and subsequent air pollution of coal-fired power station, based on the analysis to size and density distribution of particles from the recirculating load of the classifier of pulverizer, the separation experiment on sampling material from power plant with a dilute phase fluidized bed to remove pyrite and other minerals and numerical simulation on the separation process were done. The results show that the minimum fluidization velocity is 1.62 cm/s. Pyrite and other minerals in the material are separated. Ash of the upper and bottom layer material account for 33.34% and 73.42% respectively and sulfur content occupy 1.12% and 8.96% respectively. Scanning electron microscopy and spectroscopy tests show that sulfur in the bottom material exist in the form of pyrite. Numerical simulation on the flow field form of the dilute phase separation bed with gas-solid two phase and particle motion verifies the experimental results.

Modeling on Flash Flood Disaster Induced by Bed Load

Flash floods result from a complex interaction among hydro-meteorological, hydrological, and hydraulic processes across various spatial and temporal scales. Sichuan Province suffers flash floods frequently owing to mountain weather and topography. A flash flood and gravel bed load transport are two key relative problems in mountain river engineering. Bed materials are often encountered in alternate scouring and deposition in mountain fluvial processes during a flash flood. In this circumstance, CRS-1 bed load numerical model jointly with scale physical model is employed to predict water level and gravel bed scour and deposition for design of flood control dykes and flash flood disaster mitigation. A case study on the mechanism of a flash flood disaster induced by bed load transport for a hydropower station in Sichuan Province is conducted. Finally, suggestions to protect the hydropower station are proposed.

Concept of Bed Roughness Boundary Layer and Its Application to Bed Load Transport in Flow with Non-Submerged Vegetation

Ecosystem conservation has become one of the purposes in river management as well as flood mitigation and water resources management, and understanding of river flow and morphology in a stream with vegetation becomes important. Recently 2D depth averaged analysis is familiar even in a stream with vegetation by taking account of form drag due to vegetation. However, the shear stress in vegetated area is not properly described because the resistance law due to bed roughness is not reasonably modified in vegetated area. In this study, we discussed the bed roughness boundary layer in flow with non-submerged vegetation to deduce a reasonable relation between U and u* in vegetated area toward improving the analysis of sediment transport. The results show that the modification of resistance law using by thickness, velocity distribution in that layer was found to bring significant improvement of accurate estimation of shear velocity and subsequently the sediment transport. The proposed modification is improved by 2D depth averaged analysis based on this concept, and its application is certificated through flume experiment.

Unsymmetrical load effect of geologically inclined bedding strata on tunnels of passenger dedicated lines

In order to study the unsymmetrical load effect in geological bedding strata for the Muzhailing tunnel on the Lanzhou-Chongqing passenger dedicated line in China, we investigated the deformation, mechanical response and pressure of the surrounding rock and the mechanical characteristics of bolts of the tunnel. The results suggest that open zones appear at arch and invert where joints open up, when layered stratum is horizontal, or when the dip angle of in- clined bedding is small. Open zones occur perpendicular to a joint. The failure mode is bending disjunction at the arch tain shear displacement, and lead to obvious geological bedding unsymmetrical load. The failure mode is shear damage. For the joint dip angle in the range of 75-90°, the failure mode is flexural crushing at the wall and vertical shear rup- ture above the arch. The restraining effect of two sides weakens for vertical dip. On the whole, shear failure instabilitytrend would occur and the tunnel collapses evenly. When the angle between the bolt and structure plane is greater than 23°, bolts can enhance the shearing stiffness of joint plane. Unfortunately, in the general purpose graph of tunnel for 250 km/h of passenger dedicated lines, the bolts have equal length and spacing. The rationale behind this is worthy offurther study. For inclined bedding, the surrounding rock pressure at the left wall is more than that at the right wall. In addition, lining is likely to be damaged at left shoulder and side wall. With the dip angle increasing, the unsymmetrical load gradually achieves symmetry. Asymmetry design for support is recommended to reduce the unsymmetrical load on lining disturbed by excavation.

Non-capacity transport of non-uniform bed load sediment in alluvial rivers

Rivers often witness non-uniform bed load sedim ent transport. For a long tim e, non-uniform bed load transport has been assum ed to be at capacity regime determined exclusively by local flow. Yet whether the capacity assumption for non-uniform bed load transport is justified remains poorly understood. Here, the relative time scale of non-uniform bed load transport is evaluated and non-capacity and capacity models are compared for both aggradation and degradation cases with observed data. As characterized by its relative time scale, the adaptation of non-uniform bed load to capacity regime should be fulfilled quickly. However, changes in the flow and sedim ent inputs from upstream or tributaries hinder the adaptation. Also, the adaptation to capacity regime is size dependent, the finer the sediment size the slower the adaptation is, and vice versa. It is shown that the capacity model may entail considerable errors compared to the non-capacity model. For modelling of non-uniform bed load, non-capacity modelling is recommended, in which the temporal and spatial scales required for adaptation are explicitly appreciated.

Experimental study on bed-load sediment transport under irregular wave and current combined flow

Using an irregularly oscillating tray and flume, a series of experiments are completed to evaluate bed-load sediment transport rate under irregular wave -current coexistent field. Testing conditions include three interaction angles 0', 45', 90' and two kinds of median sizes (0.38 and 1.10 mm). The results of transport rate show that the net sediment transport rate can be expressed approximately as the function of the maximum bottom shear stress of waves, mean shear stress of current and the grain size.

The Effect of Grain Properties for Artificial Bedload Supply

The German Armed Forces University in Munich has conducted experiments in a laboratory flume to determine the influence of roundness on bed load transport.The investigations were assigned by the German Federal Institute of Hydrology (BIG),with a focus on incipient motion,transport velocity and the depth of mixing into the riverbed. The results of the experiments show that the transport velocity of angular graim is lower and the critical shear stress for incipient motion is slightly higher than of well-ro...

Comparison of depth-averaged concentration and bed load flux sediment transport models of dam-break flow

This paper presents numerical simulations of dam-break flow over a movable bed. Two different mathematical models were compared: a fully coupled formulation of shallow water equations with erosion and deposition terms(a depth-averaged concentration flux model), and shallow water equations with a fully coupled Exner equation(a bed load flux model). Both models were discretized using the cell-centered finite volume method, and a second-order Godunov-type scheme was used to solve the equations. The numerical flux was calculated using a Harten, Lax, and van Leer approximate Riemann solver with the contact wave restored(HLLC). A novel slope source term treatment that considers the density change was introduced to the depth-averaged concentration flux model to obtain higher-order accuracy. A source term that accounts for the sediment flux was added to the bed load flux model to reflect the influence of sediment movement on the momentum of the water. In a onedimensional test case, a sensitivity study on different model parameters was carried out. For the depth-averaged concentration flux model,Manning's coefficient and sediment porosity values showed an almost linear relationship with the bottom change, and for the bed load flux model, the sediment porosity was identified as the most sensitive parameter. The capabilities and limitations of both model concepts are demonstrated in a benchmark experimental test case dealing with dam-break flow over variable bed topography.

Incipient Motion of Bed Load Sediment and Its Application in the Three Gorges Hydro Project

A formula for the calculation of critical shear stress of non-uniform sediment was developed by on river bed has been established while considering the relative degree expose of grains.The values of the Coeffi- cients in the formula have been obtained.The formula has been verified with prototype observation data, flumes data of both nature and lightweight sand.The computed results fully reflect the characteristics of thresh- old motion of non-uniform sediment.The result have been used in some physical mo...

An Investigation of Stochastic Nature of Bed Load Motion in Chuanjiang River

Bed Load Motion and its transport rate is one of the basic issues in river dynamics.In this paper, the authors discussed the stochastic nature of bed load motion in Chuanjiang River in details.Chungjiang lies in the upstream reach of Yangtze River.Its stochastic nature is shown in the following four aspects.Firstly, even though all the conditions are the same,due to the fluctuation of the flow,the bed load discharge and the location and width of sediment transport belts are different.Secondly,during the ...

Evaluation and Improvement of Bed Load Formula Using Tapi River Data, India

The effect of non uniformity of bed material on the sediment transport has been studied by various investi-gators in the past. In the present paper the bed load transport rate has been estimated for non uniform bed material considering the various variables like discharge, hydraulic mean depth, flow velocity, bed slope, average diameter of particle etc. by collecting field data of Tapi River. The majority of the bed load formulae represent a functional relationship between bed load discharge and shear stress. This study focuses on evalu-ating the bed load using Einstein’s formulae. The bed load of pre monsoon season is estimated using various field parameters. The mathematical model has been developed using effective shear stress and bed load dis-charge. The statistical analysis, multiple regression and curve fitting (by nonlinear square fitter) is carried out using allometric function of Micro cal Origin 7.5. The proposed model has been tested using five years field data of Tapi River other than that used for the development of model. The value of rmse is close to zero in-dicates a perfect fit between measured and predicted values. The inequality coefficient is close to 0.50 sug-gest moderate relationship between estimated and computed bed load.

Development of All-Weather and Real-Time Bottom-Mounted Monitor of Bed Load Quantity

Quantity of bed load is an important physical parameter in sediment transport research. Aiming at the difficulties in the bed load measurement, this paper develops a bottom-mounted monitor to measure the bed load transport rate by adopting the sedimentation pit method and resolving such key problems as weighing and desilting, which can achieve long-time, all-weather and real-time telemeasurement of the bed load transport rate of plain rivers, estuaries and coasts. Both laboratory and field tests show that this monitor is reasonable in design, stable in properties and convenient in measurement, and it can be used to monitor the bed load transport rate in practical projects.

An Experimental Study on the Response of Model Footings Resting on Reinforced Flyash Beds under Repeated Load

In the present, investigation an attempt has been made to study the behavior of reinforced flyash beds under repeated loads by conducting carefully designed experiments. The main objective of the experiments conducted is to evaluate the beneficial effects of using flyash as a backfill material in the reinforced beds to resist the repeated loads. As no unique standard equipment is available for the application of repeated loads and to measure the response, different researchers have designed and fabricated different types of equipments for the testing. The repeated load of known intensity with waveform type and frequency is applied on the surface and embedded footing in unreinforced and reinforced flyash beds. The response of the flyash beds, in the form of settlement is measured using linear variable differential transducers. The experimental results clearly demonstrated that the provision of reinforcement in the flyash beds is effective in improving the performance of both surface and embedded footing under repeated loads.

Bed Load Transport Rate in Scouring and Armoring Process of Non-uniform Sediment River Bed

Flume experiments were carried out to study bed load transport rate during rive bed scouring and ar- moring.A theoretical differential equation linking the transport rate to the probability of incipient motion of non-uniform sediment is solved.The transport rate is shown to decrease exponentially with time,according to the theory,which is in good agreement with the experiment data.

ANALYSIS OF SALTATION CHARACTERISTICS OF BED-LOAD TRANSPORT IN FLOWING WATER

The study of bed-load transport is of great significance both in theory and in practice.This paper discusses the saltation of bed-load solid grains in flowing water.Experiments and theoretic analysis have been made by means of high-speed photographing and advanced data processing technique with a combined method based on mechanical and statistical theories.It indicates that the saltation is the main form of the bed-load transport under ordinary flowing conditions.In the meantime,taking suecessive saltation as the mod- el of bed-load transport,systematic analysis has been made with regard to the kinematic properties and mechanism of saltation.The statistical analysis shows that the probability density functions of the relative height and length of saltation are in conformity with Γ-type distribution,while the probability density func- tions of the relative velocities of saltation are in conformity with the Gaussian distribution.

各向异性层状千枚岩渗透率及有效孔隙率试验研究

采用岩石全自动三轴伺服仪和气体渗透装置,对一种致密的各向异性层状千枚岩开展了气体渗透率及有效孔隙率试验,研究常规三轴压缩和围压循环加卸载2种应力路径下,气体渗透率、有效孔隙率随层理倾角及偏应力的演化规律。结果表明:围压相同时,岩样的初始气体渗透率k0随着层理倾角β的增大呈“W”型变化,在β=45°时取最大值;在围压循环加卸载过程中,气体渗透率先随围压的加载而减小,后随围压的卸载而增大,卸载时的气体渗透率小于加载时的渗透率;层状千枚岩有效孔隙率和气体渗透率呈指数关系;平行于层理方向的气体渗透率远大于垂直于层理方向的气体渗透率;岩样有效孔隙率和气体渗透率随偏应力变化经历初始压密阶段、线弹性阶段和塑性变形阶段,随着偏应力的增大,岩样有效孔隙率和气体渗透率先减小,接着保持稳定,最后快速增大,并在岩样应力-应变曲线斜率接近于0时达到最大值。