Numerical modelling of alternate bar formation,development and sediment sorting in straight channels

A two‐dimensional shallow water hydro‐sediment‐morphodynamic model is applied to investigate alternate bar formation, development and sediment sorting in straight channels. The model is coupled, explicitly incorporating the flow–sediment–bed interactions by using the full mass and momentum conservation equations, which are numerically solved by a well‐balanced version of the finite volume Slope Limiter Centred (SLIC) scheme. The model is first tested against a flume experiment on alternate bars formed over a uniform sediment bed, which clearly exhibits processes of bar formation, migrating and finally approaching an equilibrium state. Then it is applied to another flume experiment on alternate bars due to non‐uniform sediment transport. The computational results are evaluated, with a focus on the longitudinal and vertical sediment sorting. It is argued for the first time that the inconsistent sediment sorting patterns observed in previous studies are determined by different sediment transport conditions, i.e. full versus partial transport. When a condition of full transport is achieved, under which all size fractions are fully mobilized and transported, the longitudinal surface sediment shows a sorting pattern of coarse‐on‐head and fine‐in‐pool, and the vertical substrate sediment exhibits an immobile‐fine‐coarse structure upwards. In contrast, for a partial transport condition, under which only finer fraction participates in the transport process, an opposite longitudinal pattern (i.e. fine‐on‐head and coarse‐in‐pool) and a different vertical structure (i.e. immobile‐coarse‐fine) are observed. Concurrently, numerical experiments with specified conditions show that the critical aspect ratio for the formation of migrating alternate bars is approximately equal to 12. With the increase of the aspect ratio, the bar length grows gradually, while the bar height increases rapidly for moderate values of the aspect ratio and then keeps nearly stable. The bar celerity, however, is weakly sensitive to the variation of this ratio. Copyright © 2016 John Wiley & Sons, Ltd.     

2-D SIMULATION OF CHANNEL FLOWS WITH MOVEABLE BED

1 INTRODUCTIONFor many hydraulic engineering problems, the analysis of flow and bed level variations in openchannels is a fundamental prerequisite. forcal methOds fOr alluvial rivers are well develoPednowadays as far as onediInensional descriPtions are concemed. A cOmPrhensive analysis of Ihe wellknown models is Presented by Habersack(l998). HOwever, for a number of Problems such as channelwidening, flow pattem close to sPuds and etc. a more deailed knowledge of the bed level behavio…     

DEPTH-AVERAGED 2-D CALCULATION OF TIDAL FLOW,SALINITY AND COHESIVE SEDIMENT TRANSPORT IN ESTUARIES

A depth-averaged 2-D numerical model for unsteady flow, salinity and cohesive sediment transport in estuaries is established using the finite volume method on the non-staggered, curvilinear grid. The convection terms are discretized by upwind schemes, the diffusion terms are by the central difference scheme, and the time derivative terms are by the three-time-level implicit scheme. The coupling of flow velocity and water level in the 2-D shallow water equations is achieved by the SIMPLEC algorithm with the Rhie and Chow's momentum interpolation method. The sediment model calculates the settling, deposition, erosion and transport of cohesive sediment, taking into account the influence of sediment size, sediment concentration, salinity and turbulence intensity on the flocculation of cohesive sediment. The flow model is first tested against the measurement data in the Tokyo Bay and San Francisco Bay, showing good agreements. And then, the entire model of flow, salinity and sediment transport is verified in the Gironde Estuary. The water elevation, flow velocity, salinity and sediment concentration are well predicted.     

THREE-DIMENSIONAL NUMERICAL MODELING OF COHESIVE SEDIMENT TRANSPORT BY TIDAL CURRENT IN PEARL RIVER ESTUARV

IINTRODUCTIONEstUariesareprominentcoastalfeatUres.Estuariesareofgreateconomicssignificancetomankind.Attheseareas,manyharborsandwaterchannelshavetobebuiltforeconomicpurposes.ThedesignandconstrUctionofcoastalstrUctUresinestUariesrequireknowledgeofhydrodynamicsaswellassedimenttransportinsuchregions.ThenatUreofestuariesiscontrolledbyvariouscoastalhydrodynamicprocesses.Undertheactionofhydrodynamics,sedimentdepositionsorerosionswilloccurinestuariesornearcoastalstrUCtures.Tomaintainnavigati…     

2-D NUMERICAL MODEL FOR FLOOD FORECAST IN RIVERS WITH HEAVILY SEDIMENT LADEN FLOW

I.INTRODUCTIONTwo-dimensionalnumericalmodelisaPOwerfoltoolforengineersandriVermanagerstopredictfloodhydxaulics,identifyareasofinundation,anddesignoptionsforfloodcontrollingstructures.SomespecialproblemswithheavilysedimentladenflowriVershouldbecarefullyconsideredforthenumericalmodeldesigning;1.theplaneformofariVerisusuallybraidedanditsmainchannelshiftsoften.Themainchannelandbarreplaceeachotherseveraltimesinonefloodevent.Atagivenlocationthewaterdepthmaychangefromover10meterstoseveralcenhm…     

NUMERICAL SIMULATION OF SEDIMENT RELEASE FROM RESERVOIRS

1 INTRODUCTION Reservoirs sedimentation is a serious problem in many countries, including the I. R. of Iran. Accumulation of sediment deposits decreases worldwide reservoir storage capacity by one percent per year (Mahmood, 1987). The loss of reservoir st…     

A time-splitting pressure-correction projection method for complete two-fluid modeling of a local scour hole

A two-dimensional vertical (2DV), Eulerian two-phase model or complete two-fluid model of the free surface flow was developed to simulate water-sediment flow in a local scour hole. In the model, the complete forms of the vertical, two-dimensional, two-fluid Navier-Stokes equations were discretized using a finite volume scheme. This discretization was done based on a standard staggered grid system using a curvilinear network system in compliance with the bed boundaries and water level. At the beginning of the computational cycle, the equations governing the fluid phase were solved based on the two-step projection method with a pressure-correction technique. In the first step, the intermediate fluid velocities were obtained by solving different phases of the momentum equations of the fluid phase using the time-splitting technique. In the second step, pressure was obtained and fluid velocities were updated. In this step a simple discretization method was applied for decreasing the computational complexity. After obtaining all the fluid phase variables at a new time step, the sediment phase momentum equations were solved using the time-splitting technique and sediment velocities were obtained. Then, at the end of the computational cycle, the sediment phase mass equation was solved and the concentrations of both phases were updated. At last, the capacity of the model for simulating of the longitudinal fluid velocity and sediment concentration in a local scour hole was evaluated. Numerical results were found to be in good agreement with experimental data.     

Modelling dam-break flows over mobile beds using a 2D coupled approach

Dam-break flows usually propagate along rivers and floodplains, where the processes of fluid flow, sediment transport and bed evolution are closely linked. However, the majority of existing two-dimensional (2D) models used to simulate dam-break flows are only applicable to fixed beds. Details are given in this paper of the development of a 2D morphodynamic model for predicting dam-break flows over mobile beds. In this model, the common 2D shallow water equations are modified, so that the effects of sediment concentrations and bed evolution on the flood wave propagation can be considered. These equations are used together with the non-equilibrium transport equations for graded sediments and the equation of bed evolution. The governing equations are solved using a matrix method, thus the hydrodynamic, sediment transport and morphological processes can be jointly solved. The model employs an unstructured finite volume algorithm, with an approximate Riemann solver, based on the Roe-MUSCL scheme. A predictor–corrector scheme is used in time stepping, leading to a second-order accurate solution in both time and space. In addition, the model considers the adjustment process of bed material composition during the morphological evolution process. The model was first verified against results from existing numerical models and laboratory experiments. It was then used to simulate dam-break flows over a fixed bed and a mobile bed to examine the differences in the predicted flood wave speed and depth. The effects of bed material size distributions on the flood flow and bed evolution were also investigated. The results indicate that there is a great difference between the dam-break flow predictions made over a fixed bed and a mobile bed. At the initial stage of a dam-break flow, the rate of bed evolution could be comparable to that of water depth change. Therefore, it is often necessary to employ the turbid water governing equations using a coupled approach for simulating dam-break flows.     

MODELING UNSTEADY TURBIDITY CURRENTS IN A STEEP NARROW RESERVOIR

1.IN~DUCTIONTurbiditycurrentisoneclassofflowsnameddensitycurrentorgravitycurrent(therHunterRouse(Yih(1980)),whichdePictstheintmsionofheaVyfluidintoalighterone.Usually,thedensitydifferencebetWeentWonuidisrelativelysmallandmixingacrosstheimerfaceoccurs.ThedrivingforceofdensitycurrentsisnotdensitydifferenceitselfbutthedifferenceinspeCmcweights.Turbiditycurrentisnamedwhenthedensitydifferenceisespeciallycausedbysuspendedfinesedimentparticles.Sincesediment-ladenflowcaninteraCtwiththelowerbou…     

EFFECT OF TIDAL FLOW ON ESTUARINE SEDIMENTATION

1. INTRODUCTIONIn the regulation of an estuary and the maintenance of a tidal channel for navigation, it is highlyimportant to investigate the interplay between tide and river flow and to anticipate its effects on sedimentation. In this paper. sediment--laden flow in an estuary is treated as a case of two--phase flow. resulting in a set of four governing equations. A theoretical solution is then derived for a simplified casefor qualitative assessment. Numerical computations are also carri…     

Simulation of submerged granular flows by development of a two-phase incompressible explicit mesh-free method

The interaction between fluid and sediment particles is widely involved in hydraulic engineering problems. In the current study, an explicit incompressible mesh-free method in the framework of the Moving Particle Semi-implicit(MPS) method is proposed to simulate the interaction between the two phases in submerged conditions. The proposed method solves two sets of the continuity and momentum equations, respectively, for the fluid phase and the sediment phase according to the mixture theory. In th...     

A finite element approach to watershed runoff

A numerical finite element based implicit recurrence relationship is developed for the analysis of watershed direct runoff problems. The algorithm is presented in a two dimensional form for the full momentum and continuity equations and unidimensionally for the kinematic wave theory. Since the whole domain is represented and solved as a single set of matrix equations, the advent of shocks for domains where appreciable changes in slope and/or resistance to flow occur, is immediately distinguishable.     

Numerical investigation of the factors influencing the vertical profiles of current,salinity, and SSC within a turbidity maximum zone

The flow-sediment interaction plays a considerable role on the vertical (internal) profiles of current,salinity and suspended sediment concentration (SSC) within a turbidity maximum zone (TMZ).Numerical modeling provides valuable insights into the complex estuarine physical processes.By combining numerical modeling with field observations,the influencing factors of fine sediment dynamics within the TMZ of Yangtze Estuary have been explored in this study.Firstly,during the neap tide,the measured data present that the current is too weak to break the density stratification,and the vertical flow structure is effectively altered.Secondly,a three-dimensional numerical model based on the Delft3D has been developed and a range of numerical sensitivity analyses were carried out to distinguish the dominant mechanisms and physical processes responsible for the phenomena observed from the measurement data.The numerical investigation highlights the following findings.(1) The vertical profile of currents within the TMZ is largely affected by saltwater intrusion,especially during lower currents when the baroclinic pressure gradient can significantly reshape the local vertical profiles of velocity.(2) The baroclinic effects are primarily determined by the stratification of salinity.(3) In addition to salinity,SSC also influences the local density stratification when its contribution to fluid density is comparable to that of salinity.(4) The settling velocity determines the overall sediment distribution and vertical profiles of the SSC in the water column.The SSC-dependent settling velocity (including the flocculation-induced acceleration and hindered settling deceleration phases) affects the longitudinal movement of the sediments.(5) The vertical profiles of current,salinity and SSC within the TMZ are highly associated with the turbulence determined by the model.The approach to modulate the vertical eddy viscosity in the model,based on the empirical dependency between Rig and Prt,may lead to a numerical instability in the stratified flow.In order to improve the stratification of SSC,additional turbulence damping effect is suggested to be implemented in the model.     

Application of Hausdorff fractal derivative to the determination of the vertical sediment concentration distribution

The Rouse formula and its variants have been widely used to calculate the steady-state vertical concentration distribution for suspended sediment in steady sediment-laden flows, where the diffusive flux is assumed to be Fickian. Turbulent flow, however, exhibits fractal properties, leading to non-Fickian diffusive flux for sediment particles. To characterize non-Fickian dynamics of suspended sediment, the current study proposes a Hausdorff fractal derivative based advection-dispersion equation(H...     

TWO DIMENSIONAL NUMERICAL SIMULATION OF FLOW AND GEO-MORPHOLOGICAL PROCESSES NEAR HEADLANDS BY USING UNSTRUCTURED GRID

1 INTRODUCTION In recent years, due to the increase in population and industrial developments, mankind has faced manyproblems associated with rivers, coastal waters and reservoirs. Some of these problems are flood control,water supply, power generation, and irrigation. In addition, making new hydraulic structures changesnatural conditions. Prediction of these changes is necessary for designing such constructions. For solutionof these problems usually an assessment of flow pattern, sedim…     

Significance of particle interaction to the modelling of cohesive sediment transport in rivers

Aggregation processes of fine sediments have rarely been integrated in numerical simulations of cohesive sediment transport in riverine systems. These processes, however, can significantly alter the hydrodynamic characteristics of suspended particulate matter (SPM), modifying the particle settling velocity, which is one of the most important parameters in modelling suspended sediment dynamics. The present paper presents data from field measurements and an approach to integrate particle aggregation in a hydrodynamic sediment transport model. The aggregation term used represents the interaction of multiple sediment classes (fractions) with corresponding multiple deposition behaviour. The k–ε–turbulence model was used to calculate the coefficient of vertical turbulent mixing needed for the two‐dimensional vertical‐plane simulations. The model has been applied to transport and deposition of tracer particles and natural SPM in a lake‐outlet lowland river (Spree River, Germany). The results of simulations were evaluated by comparison with field data obtained for two levels of river discharge. Experimental data for both discharge levels showed that under the prevailing uniform hydraulic conditions along the river reach, the settling velocity distribution did not change significantly downstream, whereas the amount of SPM declined. It was also shown that higher flow velocities (higher fluid shear) resulted in higher proportions of fast settling SPM fractions. We conclude that in accordance with the respective prevailing turbulence structures, typical aggregation mechanisms occur that continuously generate similar distribution patterns, including particles that settle toward the river bed and thus mainly contribute to the observed decline in the total SPM concentration. In order to determine time‐scales of aggregation and related mass fluxes between the settling velocity fractions, results of model simulations were fitted to experimental data for total SPM concentration and of settling velocity frequency distributions. The comparison with simulations for the case of non‐interacting fractions clearly demonstrated the practical significance of particle interaction for a more realistic modelling of cohesive sediment and contaminant transport. Copyright © 2004 John Wiley & Sons, Ltd.     

MODELING OF THE HIGH CONCENTRATION LAYER OF COHESIVE SEDIMENT UNDER THE ACTION OF WAVES AND CURRENTS

MODELING OF THE HIGH CONCENTRATION LAYEROF COIIESIVE SEDMNT UNDER Tus ACTIONOF WAVES AND CURRENTSQinghe ZHANG', Yongsheng WU', Jiian LIAN1 and Pingxing DING3Abstract:High concenhation layer of cohesive sediment frequenhy occurs in muddy estUaries and coastalzones, and causes raPid siltation of the waterways. A one dimensional vertical coupled modeldescribing the interactions betWeen waves, currentS and suspended cohesive sediment is develoPed inthe pre…     

THREE DIMENSIONAL MATHEMATICAL MODEL OF SEDIMENT TRANSPORT IN ESTUARINE REGIONS -A CASE STUDY OF THE HAIHE RIVER MOUTH

1 INTRODUCTIONUnderstanding the flow characters in open channel or other water area, whether natufal water bodies(such as rivers and estuaries) or man-made strUctures (such as navigational channels, harbors), isimportant for addressing numerous hydraulic engineering problems. These include the selechons ofsuitable waste disposal sites, contndnant transport, sediment transport and other ecological problems.The major difficulty for solving these problems is that it often requires an optimal…     

A model based on Hirano-Exner equations for two-dimensional transient flows over heterogeneous erodible beds

In order to study the morphological evolution of river beds composed of heterogeneous material, the interaction among the different grain sizes must be taken into account. In this paper, these equations are combined with the two-dimensional shallow water equations to describe the flow field. The resulting system of equations can be solved in two ways: (i) in a coupled way, solving flow and sediment equations simultaneously at a given time-step or (ii) in an uncoupled manner by first solving the flow field and using the magnitudes obtained at each time-step to update the channel morphology (bed and surface composition). The coupled strategy is preferable when dealing with strong and quick interactions between the flow field, the bed evolution and the different particle sizes present on the bed surface. A number of numerical difficulties arise from solving the fully coupled system of equations. These problems are reduced by means of a weakly-coupled strategy to numerically estimate the wave celerities containing the information of the bed and the grain sizes present on the bed. Hence, a two-dimensional numerical scheme able to simulate in a self-stable way the unsteady morphological evolution of channels formed by cohesionless grain size mixtures is presented. The coupling technique is simplified without decreasing the number of waves involved in the numerical scheme but by simplifying their definitions. The numerical results are satisfactorily tested with synthetic cases and against experimental data.