Coupled Model of Two-phase Debris Flow,Sediment Transport and Morphological Evolution

The volume fraction of the solid and liquid phase of debris flows, which evolves simultaneously across terrains, largely determines the dynamic property of debris flows. The entrainment process significantly influences the amplitude of the volume fraction. In this paper, we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth, the solid and fluid volume fractions and the bed morphological evolution. The model employs the Mohr–Coulomb plasticity for the solid stress, and the fluid stress is modeled as a Newtonian viscous stress. The interfacial momentum transfer includes viscous drag and buoyancy. A new extended entrainment rate formula that satisfies the boundary momentum jump condition （Iverson and Ouyang, 2015） is presented. In this formula, the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models. A finite volume method using Roe’s Riemann approximation is suggested to solve the equations. Three computational cases are conducted and compared with experiments or previous results. The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.

A Higher-Efficient Non-Hydrostatic Finite Volume Model for Strong Three-Dimensional Free Surface Flows and Sediment Transport

In order to accurately simulate strong three-dimensional （3-D） free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection-diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional （2-D） grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement.

THE STUDY ON THREE-DIMENSIONAL MATHEMATICAL MODEL OF RIVER BED EROSION FOR WATER-SEDIMENT TWO-PHASE FLOW

Based on the tensor analysis of water-sediment two-phase how, the basic model equations for clear water flow and sediment-laden flow are deduced in the general curve coordinates for natural water variable-density turbulent how. Furthermore, corresponding boundary conditions are also presented in connection with the composition and movement of non-uniform bed material. The theoretical results are applied to the calculation of the float open caisson in the construction period and good results are obtained.

Numerical Modeling of Sediment Transport and Its Effect on Algal Biomass Distribution in Lake Pontchartrain Due to Flood Release from Bonnet CarréSpillway

In order to protect the city of New Orleans from the Mississippi River flooding, the Bonnet Carré Spillway (BCS) was constructed from 1929 to 1936 to divert flood water from the river into Lake Pontchartrain and then into the Gulf of Mexico. During the BCS opening for flood release, large amounts of freshwater, nutrients, sediment, etc. were discharged into Lake Pontchartrain, and caused a lot of environmental problems. To evaluate the environmental impacts of the flood water on lake ecosystems, a two-dimensional numerical model was developed based on CCHE2D and applied to simulate the flow circulation, sediment transport and algal biomass distribution in Lake Pontchartrain. The effect of sediment concentration on the growth of algae was considered in the model. The numerical model was calibrated using field measured data provided by USGS, and then it was validated by the BCS Opening Event in 1997. The simulated results were generally in good agreement with filed data and satellite imagery. The field observation and numerical model show that during the spillway opening for flood release, the sediment concentration is very high, which greatly restricts the growth of algae, so there is no algal bloom observed in the lake. After the closure of BCS, the sediment concentration in the lake reduces gradually, and the nutrient concentration of the lake is still high. Under these conditions, numerical results and satellite imagery showed that the chlorophyll concentration was high and algal bloom might occur.

Application of the Sediment Mathematical Modelling on Planned Project of Lingdingyang Bay

Based on the combined hydraulic calculation for the eastern network region at the Pearl River estuary and several outlets to the Lingdingyang Bay, the sediment calculation modelling was introduced in the establishment of the sediment mathematical model for Lingdingyang Bay and the eastern region with one and two dimensional flow calculation. Model adjustment and verification were performed in conjunction with field data. The simulated results coincide well with measured data.In addition the model is applied to predict the shore-line planning scheme of Lingdingyang Bay.The theoretical criterion is provided for the shore line plan in the model.And a new mathematical simulated method is put out to research the planning engineering concerned with one-dimensional net rivers and two-dimensional estuary.

Numerical simulations of flow and sediment transport within the Ning-Meng reach of the Yellow River,northern China

Effective management of a river reach requires a sound understanding of flow and sediment transport generated by varying natural and artificial runoff conditions. Flow and sediment transport within the Ning-Meng reach of the Yellow River（NMRYR）, northern China are controlled by a complex set of factors/processes, mainly including four sets of factors：（1） aeolian sediments from deserts bordering the main stream;（2） inflow of water and sediment from numerous tributaries;（3） impoundment of water by reservoir/hydro-junction; and（4） complex diversion and return of irrigation water. In this study, the 1-D flow ＆ sediment transport model developed by the Yellow River Institute of Hydraulic Research was used to simulate the flow and sediment transport within the NMRYR from 2001 to 2012. All four sets of factors that primarily control the flow and sediment transport mentioned above were considered in this model. Compared to the measured data collected from the hydrological stations along the NMRYR, the simulated flow and sediment transport values were generally acceptable, with relative mean deviation between measured and simulated values of 〈15%. However, simulated sediment concentration and siltation values within two sub-reaches（i.e., Qingtongxia Reservoir to Bayan Gol Hydrological Station and Bayan Gol Hydrological Station to Toudaoguai Hydrological Station） for some periods exhibited relatively large errors（the relative mean deviations between measured and simulated values of 18% and 25%, respectively）. These errors are presumably related to the inability to accurately determine the quantity of aeolian sediment influx to the river reach and the inflow of water from the ten ephemeral tributaries. This study may provide some valuable insights into the numerical simulations of flow and sediment transport in large watersheds and also provide a useful model for the effective management of the NMRYR.

A Numerical Model for Natural Backfill of Pipeline Trenches Subjected to Unidirectional/Oscillatory Flows

A numerical model for the self-burial of a pipeline trench is developed. Morphological evolutions of a pipeline trench under steady-current or oscillatory-flow conditions are simulated with/without a pipeline inside the trench. The oscillatory flow in this study represents the action of waves. The two-dimensional Reynolds-averaged continuity and Navier-Stokes equations with the standard k-ε turbulence closure, as well as the sediment transport equations, are solved with the finite difference method in a curvilinear coordinate system. Both bed and suspended loads of sediment transport are included in the morphological model. Because of the lack of experimental data on the backfilling of pipeline trenches, the numerical model is firstly verified against three closely-relevant experiments available in literature. A detailed measurement of the channel migration phenomenon under steady currents is employed for the assessment of the integral performance of the model. The two experimental results from U-tube tests are used to validate the model's ability in predicting oscillatory flows. Different time-marching schemes are employed for the morphological computation under unidirectional and oscillatory conditions. It is found that vortex motions within the trench play an important role in the trench development.

Three-dimensional CFD simulation of inlet structure flow in pumping station based on Eulerian solid- liquid two-phase flow model

Sediment deposition in the pumping station has a huge negative impact on unit operation.The three-dimensional CFD method has been used to simulate inlet structure flow in pumping station based on the Eulerian solid- liquid two-phase flow model. The numerical results of the preliminary scheme show that sediment deposition occurs in the forebay of pumping station because of poor flow pattern therein. In order to improve hydraulic configuration in the forebay,one modified measure reconstructs water diversion weir shape,and another measure sets a water retaining sill in the approach channel. The simulation results of the modified scheme prove that back flow in the forebay has been eliminated and the sediment deposition region has also been reduced.

INFLUENCES OF THE FISH-MOUTH PROJECT AND THE GROINS ON THE FLOW AND SEDIMENT RATIO OF THE YANGTZE RIVER WATERWAY

A depth_integrated two_dimensional numerical model of current, salinity and sediment transport was proposed and calibrated by the observation data in the Yangtze River Estuary. It was then applied to investigate the flow and sediment ratio of the navigation channel, i.e. the North Channel of the Yangtze River Estuary,before and after the first phase waterway project is implemented. Particularly, the influences of groin length and the orientation of the submerged dam on the flow ratio and sediment load discharging into the North Channel were discussed. The numerical results demonstrate that less sediment load discharges into the navigation channel, which unburdens the waterway dredging, but in the meantime the flow ratio is also decreased. The flow and sediment ratio can be adjusted by changing layout and dimensions of the hydro_structures, such as the groin length, the top height, etc. The effect of the orientation of the submerged dam is more obvious than the groin lengh.

Water Flow and Sediment Flux Forecast in the ChókwèIrrigation Scheme, Mozambique

This study sought to forecast water flow and sediment flux in the scheme as potential contributions for improved management in the Chókwè Irrigation Scheme (CIS). Fieldwork data was collected during dry (DS) and wet (WS) seasons. Flow measurement was performed at 9 stations using a calibrated flow meter OTT C31. Water flow and sediment flux from 2004 to 2019 were used. Hydrodynamic forecast simulations were performed using Mann-Kendall test and ARIMA model for determination of temporal trends. Findings suggest higher values during DS for water discharge and sediment flux. Mann-Kendall test for sediment discharge trends was not significant at 95% significance level, except for the Offtake in WS. ARIMA test for the sediment discharges, at the Intake, for DS and WS, sediments were well described by the ARIMA model and gave a good result for the sediments. Good fit between the observed and the predicted ARIMA model was found. ARIMA model for sediment discharge at CIS based on AIC has a good fit for AR (p = 1), whereby, at the Intake the ARIMA p-value was 0.822 and 0.932, for WS and DS, respectively. Whilst in the Offtake, the ARIMA p-value was 0.877 and 0.893, respectively. These results can be used to improve the CIS management, both for water flow and sediment flux.

Morphodynamics and sediment connectivity index in an unmanaged,debris-flow prone catchment:a through time perspective

Torrential processes are among the main actors responsible for sediment production and mobility in mountain catchments.For this reason,the understanding of preferential pathways for sediment routing has become a priority in hazard assessment and mitigation.In this context,the sediment Connectivity Index(IC)enables to analyse the existing linkage between sediment sources and the selected target(channel network or catchment outlet).The IC is a grid-based index that allows fast computation of sediment connectivity based on landscape information derived from a single Digital Terrain Model(DTM).The index computation is based on the log-ratio between an upslope and a downslope component,including information about drainage area,slope,terrain roughness,and distance to the analysis target(e.g.outlet).The output is a map that highlights the degree of structural connectivity of sediment pathways over analysed catchments.Until now,these maps are however rarely used to help defining debris-flow hazard maps,notably due to a lack of guidelines to interpret the IC spatial distribution.This paper proposes an exploitation procedure along profiles to extract more information from the analysis of mapped IC values.The methodology relies on the analysis of the IC and its component variables along the main channel profile,integrated with information about sediment budgeting derived from Difference of DEMs(DoD).The study of connectivity was applied in the unmanaged sub-catchment(without torrent control works)of the Rio Soial(Autonomous Province of Trento–NE Italy)to understanding the geomorphic evolution of the area after five debris flows(in ten years)and the related changes of sediment connectivity.Using a recent DTM as validation,we demonstrated how an IC analysis over the older DTM can help predicting geomorphic changes and associated hazards.The results show an IC aptitude to capture geomorphic trajectories,anticipate debris flow deposits in a specific channel location,and depict preferential routing pathways.

A High-Order Conservative Semi-Lagrangian Solver for 3D Free Surface Flows with Sediment Transport on Voronoi Meshes

In this paper,we present a conservative semi-Lagrangian scheme designed for the numeri-cal solution of 3D hydrostatic free surface flows involving sediment transport on unstruc-tured Voronoi meshes.A high-order reconstruction procedure is employed for obtaining a piecewise polynomial representation of the velocity field and sediment concentration within each control volume.This is subsequently exploited for the numerical integration of the Lagrangian trajectories needed for the discretization of the nonlinear convective and viscous terms.The presented method is fully conservative by construction,since the transported quantity or the vector field is integrated for each cell over the deformed vol-ume obtained at the foot of the characteristics that arises from all the vertexes defining the computational element.The semi-Lagrangian approach allows the numerical scheme to be unconditionally stable for what concerns the advection part of the governing equations.Furthermore,a semi-implicit discretization permits to relax the time step restriction due to the acoustic impedance,hence yielding a stability condition which depends only on the explicit discretization of the viscous terms.A decoupled approach is then employed for the hydrostatic fluid solver and the transport of suspended sediment,which is assumed to be passive.The accuracy and the robustness of the resulting conservative semi-Lagrangian scheme are assessed through a suite of test cases and compared against the analytical solu-tion whenever is known.The new numerical scheme can reach up to fourth order of accu-racy on general orthogonal meshes composed by Voronoi polygons.

基于Flow-3D的多阶明渠工程数学模型研究

[目的]研究灌区多级阶梯明渠水流的水动力特性及泥沙淤积问题,并评价其对多阶明渠工程的影响.[方法]基于VOF方法,在Flow-3D平台上建立了山东省德州市牛角店输水工程的三维数学模型,通过模拟不同工况下的三维多阶明渠水流,对比分析水流流速分布、水面线变化及渠底泥沙淤积情况.[结果]水流流态平稳后,最大流速出现在第二级阶梯上,挟沙水流流速较清水流速大0.50~0.86 m/s;在挑坎及第二级阶梯上,水面线变化幅度较大;随着时间推移,泥沙主要淤积于挑坎前及第三级阶梯底部,200 s时淤积厚度最高达0.143 mm.[结论]在牛角店输水工程施工及后期运行中,适当加固第二级阶梯渠底,提升第二级侧墙高度,并通过调控下游水位以保证灌渠的安全运行.

泥沙冲刷模拟的FLOW-3D数值方法研究

泥沙冲刷是很多水利工程中常见的现象,采用FLOW-3D的数值模拟手段来研究泥沙冲刷应用广泛。从水射流冲刷泥沙角度,对FLOW-3D泥沙冲刷时网格质量和泥沙模型参数的选择对于冲刷效果的影响进行研究,提出了一种网格的设置方法,确定了泥沙模型中影响射流冲刷泥沙过程的关键参数。结果表明:泥沙模型参数中临界体积分数越大,悬浮的泥沙量越少;水下休止角需要通过试验测量验证;挟带系数越大,随冲刷起动的沙更多;推移质系数越大,沉积物沿河床的传输量越大。研究成果可为FLOW-3D在泥沙冲刷的应用提供一定的参考。

APPLICATION OF SWAT MODEL IN THE UPSTREAM WATERSHED OF THE LUOHE RIVER

In the Huanghe (Yellow) River basin,soil erosion is a serious problem,while runoff and sediment yield simulation has not been extensively studied on the basis of GIS (Geographic Information System) and distributed hydrological model. GIS-based SWAT (Soil and Water Assessment Tool) model was used to simulate runoff and sediment in the Huanghe River basin. The objective of this paper is to examine the applicability of SWAT model in a large river basin with high sediment runoff modulus,which could reach 770t/(km2·a). A two-stage "Brute Force" optimization procedure was used to calibrate the parameters with the observed monthly flow and sediment data from 1992 to 1997,and with input parameters set during the calibration process without any change the model was validated with 1998-1999’s observed data. Coefficient of examination (R2) and Nash-Suttcliffe simulation efficiency (Ens) were used to evaluate model prediction. The evaluation coefficients for simulated flow and sediment,and observed flow and sediment were all above 0.7,which shows that SWAT model could be a useful tool for water resources and soil conservation planning in the Huanghe River basin.

Numerical modeling of sediment transport based on unsteady and steady flows by incompressible smoothed particle hydrodynamics method

The purpose of the present paper is to introduce a simple two-part multi-phase model for the sediment transport problems based on the incompressible smoothed particle hydrodynamics（ISPH） method. The proposed model simulates the movement of sediment particles in two parts. The sediment particles are classified into three categories, including the motionless particles, moving particles behave like a rigid body, and moving particles with a pseudo fluid behavior. The criterion for the classification of sediment particles is the Bingham rheological model. Verification of the present model is performed by simulation of the dam break waves on movable beds with different conditions and the bed scouring under steady flow condition. Comparison of the present model results, the experimental data and available numerical results show that it has good ability to simulate flow pattern and sediment transport.

Identifying the driving factors of sediment delivery ratio on individual flood events in a long-term monitoring headwater basin

The Sediment Delivery Ratio(SDR) has multi-fold environmental implications both in evaluating the soil and water losses and the effectiveness of conservation measures in watersheds. Various factors, including hydrological regime and watershed properties, may influence the SDR at interannual timescales. However, the effect of certain important dynamic factors, such as rainfall peak distribution, runoff erosion power and sediment bulk density, on the sediment delivery ratio of single flood events(SDRe) has received little attention. The Qiaogou headwater basin is in the hilly-gully region of the Chinese Loess Plateau, and it encompasses a 0.45 km^2 catchment. Three large-scale field runoff plots at different geomorphological positions were chosen to obtain the observation data, and the 20-year period between 1986 and 2005 is presented. The results showed that the SDRe of the Qiaogou headwaters varied from 0.49 to 2.77. Among the numerous influential factors, rainfall and runoff were the driving factors causing slope erosion and sediment transport. The rainfall erosivity had a significant positive relationship with the sediment transport modulus(R^2=0.85, P<0.01) but had no significant relationship with SDRe. The rainfall peak coefficient was significantly positively correlated with the SDRe(R^2=0.64, P<0.05), indicating the influence of rainfall energy distribution on the SDRe. The runoff erosion power index was not only significantly related to the sediment transport modulus(R^2=0.84, P<0.01) but also significantly related to the SDRe(R^2=0.57, P<0.01). In addition, the relative bulk density was significantly related to the SDRe, indicating that hyper-concentrated flow characteristics contributed to more transported sediment in the catchment. Thus, the rainfall peak coefficient, runoff erosion power and sediment relative bulk density could be used as dynamic indexes to predict the SDRe in the hilly areas of the Chinese Loess Plateau.

1-D coupled non-equilibrium sediment transport modeling for unsteady flows in the discontinuous Galerkin framework

A high-resolution, 1-D numerical model has been developed in the discontinuous Galerkin framework to simulate 1-D flow behavior, sediment transport, and morphological evaluation under unsteady flow conditions. The flow and sediment concentration variables are computed based on the one-dimensional shallow water flow equations, while empirical equations are used for entrainment and deposition processes. The sediment transport model includes the bed load and suspended load components. New formulations for Harten-Lax-van Leer （HLL） and Harten-Lax-van Contact （HLLC） are presented for shallow water flow equations that include the bed load and suspended load fluxes. The computational results for the flow and morphological changes after two dam break events are compared with the physical model tests. Results show that the modified HLL and HLLC formulations are robust and can accurately predict morphological changes in highly unsteady flows.

Application of numerical model in evaluating the effects of a planned reservoir on reach-scale bankfull discharge of the lower Wei River

Predicting the responses of an alluvial channel to changes in flow and sediment supply is essential for engineering design. Many methods have been developed in the last few decades to describe sectional bankfull characteristics(elevation and discharge); however, studies on long-term reach-scale bankfull discharge are still limited. In this study, a hydraulic model is built to calculate the reach-scale bankfull discharge, and the effects of reservoir building on downstream bankfull discharges are discussed. The studied river reach is located at the lower Wei River(WR), where the planned Dongzhuang Reservoir would be built on its largest tributary, the Jing River. A quasi-two-dimensional numerical model coupled with a bankfull discharge estimating method is put forward to calculate the reach-scale bankfull discharge. The soundness of the model is verified. Results show that the temporal variation of reach-scale bankfull discharge of the lower reach of the WR would be highly influenced by the planned reservoir, especially during the first 20 years of operation. The effect of the planned reservoir on bankfull discharge may reach its maximum when the total trapped sediment load reaches approximately 75% of the reservoir capacity. Our results show that after the first 17 years of operation,the effect of the planned reservoir on bankfull discharge of the river reach may decrease gradually.The soundness and predictive capability of the coupled model have also been calibrated by comparing with existing reservoirs. All analyses indicate that the numerical model can be used to predict the changed reach-scale bankfull discharge of the lower WR.

A numerical study for boundary layer current and sheet flow transport induced by a skewed asymmetric wave

An analytical model with essential parameters given by a two-phase numerical model is utilized to study the net boundary layer current and sediment transport under skewed asymmetric oscillatory sheet flows. The analytical model is the first instantaneous type model that can consider phase-lag and asymmetric boundary layer development. The two-phase model supplies the essential phase-lead, instantaneous erosion depth and boundary layer development for the analytical model to enhance the understanding of velocity skewness and acceleration skewness in sediment flux and transport rate. The sediment transport difference between onshore and offshore stages caused by velocity skewness or acceleration skewness is shown to illustrate the determination of net sediment transport by the analytical model. In previous studies about sediment transport in skewed asymmetric sheet flows, the generation of net sediment transport is mainly concluded to the phase-lag effect.However, the phase-lag effect is shown important but not enough for the net sediment transport, while the skewed asymmetric boundary layer development generated net boundary layer current and mobile bed effect are key important in the transport process.