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.

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.

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.

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.

Heat Transfer and Flow Characteristics Predictions with a Refined k-ε-f_u Turbulent Model in Impinging Jet

Local heat transfer and flow characteristics in a round turbulent impinging jet for Re≈23 000 is predicted numerically with the RANS approach and a k-ε-fu turbulence model. The heat transfer predictions and turbulence parameters are verified against the axis-symmetric free jet impingement measurements and compared with previous other turbulence models, and results show the k-ε-fu model has a good performance in predictions of the local wall heat transfer coefficient, and in agreement with measurements in mean velocity profiles at different radial positions as well. The numerical model is further used to examine the effect of the fully confined impingement jet on the local Nusselt number. Local Nusselt profiles in x and y-centerlines for the target plate over three separation distances are predicted. Compared with the experimental data, the numerical results are accurate in the central domain around the stagnation region and present a consistent structure distribution.

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.

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.

A three-dimensional model for the transportation of Fe and Mn in Arha Reservoir

Coupling with a three dimensional (3D) hydrodynamic model and a suspended solids model, a 3D model for the transport of Fe and Mn in Arha Reservoir, China, was developed. The 3D velocity fields for the flood season are computed to drive the 3D model of Fe and Mn in which the processes of advection, diffusion, redox, sorption, desorption, deposition, and re suspension are included. The model has been calibrated by matching observed fluid, suspended solids, and total concentrations of Fe and Mn in the water column and in the sediment, successively. The model simulated both horizontal and vertical gradients of Fe and Mn in Arha Reservoir. It was found that Fe and especially Mn stratify in accordance with the stratification of DO during summer. The redox cycles across the water sediment interface has a principal role in the rise of Fe and Mn concentrations in the overlying water. It was also found that Fe and Mn loadings from the tributaries have a carryover effect on the water quality through a secondary contamination in the reservoir.

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.

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.

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 simulation of flow around square cylinder using different low-Reynolds number turbulence models

ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged unsteady flow field.Meanwhile,drag and lift coefficients of the four different low-Reynolds number turbulence models were analyzed.The simulated results of YANG-SHIH model are close to the large eddy simulation results and experimental results,and they are significantly better than those of ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMR models.The modification of the generation of turbulence kinetic energy is the key factor to a successful simulation for YANG-SHIH model,while the correction of the turbulence near the wall has minor influence on the simulation results.For ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMA models satisfactory simulation results cannot be obtained due to lack of the modification of the generation of turbulence kinetic energy.With the joint force of wall function and the turbulence models with the adoption of corrected swirl stream,flow around a square cylinder can be fully simulated with less grids by the near-wall.

THE CATASTROPHE FORECAST OF THE ANNUAL SIDEMENT TRANSPORT BY DEBRIS FLOW IN JIANGJIA GULLY,YUNNAN PROVINCE, CHINA

On the basis of the observational data on the annual sediment transport by debris flow in recent 8 years, appling the catastrophe forecast method of Grey System Theory, this study has established the catastrophe model of the annual sediment transport by debris flow in Jiangjia Gully. It has forecasted the next potential catastrophic year in which the annual sediment transport will be over the catastrophic. threshold 2 million m3. Furthermore, it has introduced the 'equal dimension-new information model', which makes the forecast be done continuously.

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.

Development of River Meander Model

In the studies of open-channel flow with suspended sediments,used a constant of Von Karmanκin a model for velocity profile.The augmentation parameters have been added by various researchers in more recent development of the boundary-layer theory of meander development.In this research new parameters will be included because of the existence of the turbulent flow region in meandering channels because of boundary-layer theory.

Sedimentary Characteristics and Model of Gravity Flow Depositional System for the First Member of Upper Miocene Huangliu Formation in Dongfang Area, Yinggehai Basin, Northwestern South China Sea

The gravity flow deposit were mainly developed in the lowstand systems tract（LST） of the first member of Upper Miocene Huangliu Formation（Ehl1） in Dongfang area, Yinggehai Basin, has become a valuable target for gas exploration and production. The gravity flow sedimentary characteristics of lithofacies associations, sedimentary texture, seismic facies and logging facies were described in detail on the basis of integrated analysis of cores, logging and seismic data. The sedimentary microfacies types composed of neritic sandbar, continental shelf mud, main channel, bifurcated or cross-cutting distributary channel, overspill, and natural levee are revealed under the constraint of high resolution sequence stratigraphic framework in the Ehl1. The gravity flow deposit system in the LST is divided into three evolution stages corresponding to periods of three parasequence sets. The gravity flow deposit was induced in the early LST, expanded rapidly in the middle LST and decreased slightly in the late LST. But its developing scale decreased sharply in the transgression systems tract（TST） and finally vanished in the highstand systems tract（HST）. This spatial evolution rule is constrained by the integrated function of sediments supply of the Vietnam Blue River in the LST, the development of local gradient change in sea floor（micro-topography, i.e., flexure slope break）, and the fall in relative sea level. On the basics of the deep study of the coupling relationship among the three main control factors, the sedimentary model is established as an optimal component of ＂source-channel-sink＂ for shallow marine turbidite submarine fan.

旋流排沙装置排沙效率的数值模拟与试验研究

为提高水库中排沙洞的排沙效率,研发了一种新型旋流排沙装置。利用数值模拟试验和水工模型试验两种方法,研究了旋流排沙装置中叶轮组合、旋流场的流速分布特征、旋流排沙装置对排沙洞口泥沙的冲淤特性和排沙效率的影响。数值模拟结果表明,旋流排沙装置中单叶轮的流场范围随转速的增大而扩大,但当叶轮转速增加至20r/s时流场范围不再增加;旋流排沙装置中组合旋流场范围取决于叶轮的转速及其布置方式,当以三角形方式布置时,组合流场范围也随叶轮转速的增加而扩大,但叶轮间距布置过小会导致能量耗散过高,叶轮间距布置过大则不能形成有效的组合流场,旋流排沙装置的最优运行状态为叶轮距离为70mm、转速为20r/s。数值模拟和模型试验结果表明,旋流排沙装置可有效增大排沙洞的排沙效率和冲淤范围,当库区水位较深形成冲沙漏斗时,旋流排沙装置可使排沙洞的排沙效率提高71.5%;当库区水位较浅形成冲沙廊道时,旋流排沙装置能使排沙洞的排沙量增加4.02kg,说明该装置对库区水深较浅的淤积泥沙具有较好的排出效果,对解决甘肃河西水少沙多型水库内的泥沙淤积问题具有重要意义。