Simulation of the meander cut-off by 2D hydrodynamic model for erodible bed

2021 ◽  
Author(s):  
Tatiana Fedorova ◽  
Vitaly Belikov ◽  
Andrei Alabyan
Keyword(s):  
Bottom Topography ◽  
Gravitational Effect ◽  
High Water ◽  
Computational Domain ◽  
River Channels ◽  
Time Step ◽  
Channel Deformations ◽  
Computational Mathematics ◽  
Secondary Currents ◽  
Meandering Channel

<p>The retrospective simulation of the Pyoza river (Arkhangelsk region, Russia) meander cut-off in 2003-2008 has been undertaken. As a result of the river bend straightening two large villages were cut off from the road network of the region.</p><p>The initial data for modeling were obtained by analyses of archive satellite images for the period from 1997 together with the runoff data, as well as by the field survey of September 2019. The simulation was performed by the latest version of the STREAM_2D CUDA software, using a new method for the numerical solution of two-dimensional Saint-Venant equations [1]. It was adapted for the complicated bottom topography typical for a wide floodplain with a meandering channel flooded in high water stage.</p><p>The mass-exchange equations for three layers of sediment over the unerodible bed were solved together with the hydrodynamic equations. When calculating channel deformations, the gravitational effect due to bottom slope and the influence of secondary currents on the sediment shift were taken into account [2].</p><p>The Pyoza river is the lowest large tributary of the Mezen’ river flowing into the White sea. It is distinguished by a typical alluvial channel, meandering along wide floodplain composed by sands and sandy loams. The maximum runoff usually corresponds to spring snow-melting and can reach 1500-2000 m<sup>3</sup>/s.</p><p>To schematize the computational domain of the Pyoza river section of 13 km long, a hybrid grid of irregular structure was constructed, consisting of 37 329 cells of a quadrangular shape for the channel and a triangular one for the floodplain.</p><p>The simulation started at the year 1997 when where was no any rill across the meander neck. The time step of calculation was taken to be one day.</p><p>Modeling made it possible to simulate realistically the essential steps and mechanisms of the meander cut-off: the development of a pioneer straightening rill, its widening and deepening, as well as blocking of the old channel by a point bar in its upper reaches, as well as its further silting and aggradation.</p><p>1. Aleksyuk A.I., Belikov V.V. (2017): Simulation of shallow water flows with shoaling areas and bottom discontinuities. Computational Mathematics and Mathematical Physics, Volume 57, issue 2, pp. 318–339. https://doi.org/10.1134/S0965542517020026</p><p>2. Aleksyuk А. I., Belikov V. V., Borisova N. M., Fedorova T. A. (2018): Numerical modeling of non-uniform sediment transport in river channels. Water Resources, Volume 45, Special Issue S1, pp. 11–17. http://dx.doi.org/10.1134/S0097807818050275</p>

Transient bottom topography changes in alluvial streams

2000 ◽  
Vol 27 (2) ◽  
pp. 309-318
Author(s):  
Y R Fares
Keyword(s):  
Bottom Topography ◽  
Computational Modelling ◽  
Field Measurements ◽  
Computational Procedure ◽  
Bed Load ◽  
River Channels ◽  
Total Load ◽  
Secondary Currents ◽  
Sediment Bed ◽  
Transient Variations

The development of a quasi two-dimensional computational model for simulating the transient variations of bed topography profiles in alluvial river channels is reported. The formulation of the model is based on combining the longitudinal flow momentum with the continuity principle of the sediment bed load. The Engelund-Hansen formula is employed in estimating the total sediment bed load along the reach of a river channel. The lateral bed load contribution from the total load is calculated in the same way as in calculating the lateral secondary currents from the main flow velocities. The numerical scheme and the computational procedure used in the study are described in detail. The simulated bed level profiles are verified through comparisons with experimental and field measurements taken from case studies in the literature for different flow conditions, channel characteristics, and sediment properties. The correlation between flow discharge, bed load, boundary friction, and channel slope is discussed. On the basis of the reasonably good comparisons with field data, it may be deduced that the model can be used for predicting the bottom topography variations in river channels.Key words: meandering rivers, bottom topography, sediment transport, bed load, boundary roughness, field measurements, experimental data, computational modelling, finite difference method.

Application of the method of mathematical modeling for forecasting channel deformations at the underwater crossing of the main pipeline

2020 ◽  
Vol 10 (6) ◽  
pp. 599-609
Author(s):  
Valery А. Gruzdev ◽  
◽  
Georgy V. Mosolov ◽  
Ekaterina A. Sabayda ◽  
◽  
...  
Keyword(s):  
Mathematical Modeling ◽  
Mathematical Model ◽  
Roughness Coefficient ◽  
Computational Domain ◽  
Channel Deformations ◽  
Geological Surveys ◽  
Kuban River ◽  
The Stability ◽  
Protective Structures

In order to determine the possibility of using the method of mathematical modeling for making long-term forecasts of channel deformations of trunk line underwater crossing (TLUC) through water obstacles, a methodology for performing and analyzing the results of mathematical modeling of channel deformations in the TLUC zone across the Kuban River is considered. Within the framework of the work, the following tasks were solved: 1) the format and composition of the initial data necessary for mathematical modeling were determined; 2) the procedure for assigning the boundaries of the computational domain of the model was considered, the computational domain was broken down into the computational grid, the zoning of the computational domain was performed by the value of the roughness coefficient; 3) the analysis of the results of modeling the water flow was carried out without taking the bottom deformations into account, as well as modeling the bottom deformations, the specifics of the verification and calibration calculations were determined to build a reliable mathematical model; 4) considered the possibility of using the method of mathematical modeling to check the stability of the bottom in the area of TLUC in the presence of man-made dumping or protective structure. It has been established that modeling the flow hydraulics and structure of currents, making short-term forecasts of local high-altitude reshaping of the bottom, determining the tendencies of erosion and accumulation of sediments upstream and downstream of protective structures are applicable for predicting channel deformations in the zone of the TLUC. In all these cases, it is mandatory to have materials from engineering-hydro-meteorological and engineering-geological surveys in an amount sufficient to compile a reliable mathematical model.

scholarly journals Sediment Transport and Water Flow Resistance in Alluvial River Channels: Modified Model of Transport of Non-Uniform Grain-Size Sediments

Water ◽  
2021 ◽  
Vol 13 (15) ◽  
pp. 2038
Author(s):  
Gennady Gladkov ◽  
Michał Habel ◽  
Zygmunt Babiński ◽  
Pakhom Belyakov
Keyword(s):  
Sediment Transport ◽  
Water Flow ◽  
Transport Model ◽  
Flow Characteristics ◽  
Field Measurements ◽  
Field Investigation ◽  
Empirical Modeling ◽  
River Channels ◽  
Channel Deformations ◽  
East European

The paper presents recommendations for using the results obtained in sediment transport simulation and modeling of channel deformations in rivers. This work relates to the issues of empirical modeling of the water flow characteristics in natural riverbeds with a movable bottom (alluvial channels) which are extremely complex. The study shows that in the simulation of sediment transport and calculation of channel deformations in the rivers, it is expedient to use the calculation dependences of Chézy’s coefficient for assessing the roughness of the bottom sediment mixture, or the dependences of the form based on the field investigation data. Three models are most commonly used and based on the original formulas of Meyer-Peter and Müller (1948), Einstein (1950) and van Rijn (1984). This work deals with assessing the hydraulic resistance of the channel and improving the river sediment transport model in a simulation of riverbed transformation on the basis of previous research to verify it based on 296 field measurements on the Central-East European lowland rivers. The performed test calculations show that the modified van Rijn formula gives the best results from all the considered variants.

scholarly journals Numerical Simulation of Effects of River Reconstruction on Flooding: A Case Study of the Ba River, China

2021 ◽  
Vol 2021 ◽  
pp. 1-10
Author(s):  
Haixiao Jing ◽  
Yongbiao Lang ◽  
Xinhong Wang ◽  
Mingyang Yang ◽  
Zongxiao Zhang
Keyword(s):  
Flow Velocity ◽  
Eddy Currents ◽  
Human Life ◽  
River Channel ◽  
Computational Domain ◽  
River Channels ◽  
One Dimensional ◽  
Local Reconstruction ◽  
Before And After ◽  
Flood Flow

The local reconstruction of river channels may pose obstacles of flood flow, local eddy currents, or high flow velocity which pose potential threats to human life and infrastructures nearby. In the design of such projects, the effects of local reconstruction of the river channel on flooding are often evaluated by the one-dimensional method, which is based on the formula of one-dimensional nonuniform flow. In this study, a two-dimensional hydrodynamic model based on shallow water equations is employed to investigate the impacts of river reconstruction on flooding in the Ba River, China. The finite volume method and an unstructured triangular mesh are used to solve the governing equations numerically. The numerical model is validated by comparison with the results of a physical model of 1 : 120 scale. The backwater effects and impacts of flood flow fields under two flood frequencies are analyzed by comparing the numerical results before and after local reconstruction. The results show that the backwater length under both 10-year and 100-year floods can be reached up to the upstream boundary of the computational domain. However, the maximum water level rises are limited, and the levees in this river channel are safe enough. The flow velocity fields under both floods are changed obviously after local reconstruction in the Ba River. Areas with the potential for scour and deposition of the river bed are also pointed out. The findings of this study are helpful for the evaluation of flood risks of the river.

scholarly journals Hydro-morphodynamics 2D modelling using a discontinuous Galerkin discretisation

2020 ◽  
Author(s):  
Mariana Clare ◽  
James Percival ◽  
Stephan Kramer ◽  
Athanasios Angeloudis ◽  
Colin Cotter ◽  
...  
Keyword(s):  
Sediment Transport ◽  
Discontinuous Galerkin ◽  
Transport Model ◽  
Ocean Model ◽  
Suspended Sediment Transport ◽  
Secondary Currents ◽  
Meandering Channel ◽  
Complex Process ◽  
Coastal Ocean Model ◽  
High Degree

<p>The development of morphodynamic models to simulate sediment transport accurately is a challenging and highly complex process given the non-linear and coupled nature of the sediment transport problem. We implement a new depth-averaged coupled hydrodynamic and sediment transport model within the coastal ocean model Thetis, built using the code generating framework Firedrake which facilitates code flexibility and optimisation benefits. To the best of our knowledge, this represents the first full morphodynamic model using a discontinuous Galerkin based finite element discretisation, to include both bedload and suspended sediment transport. We apply our model to problems with non-cohesive sediment and account for effects of gravity and helical flow by adding slope gradient terms and parametrising secondary currents. For validation purposes and to demonstrate model capability, we present results from the common test cases of a migrating trench and a meandering channel comparing against experimental data and the widely used model Telemac-Mascaret.</p><p>There is a high degree of uncertainty associated with morphodynamic models, in part due to incomplete knowledge of various physical, empirical and numerical closure related parameters in both the hydrodynamic and morphodynamic solvers. We therefore also present examples of how an adjoint model can be used to calibrate or invert for the values of these parameters from either experimental results or real-world erosion profiles.</p>

scholarly journals Algorithms for Solving Darcian Flow in Structured Porous Media

10.14311/1829 ◽  
2013 ◽  
Vol 53 (4) ◽  
Author(s):  
Michal Kuráž ◽  
Petr Mayer
Keyword(s):  
Proper Time ◽  
Time Integration ◽  
Mass Term ◽  
Step Length ◽  
Richards Equation ◽  
Parabolic Problems ◽  
Computational Domain ◽  
Wetting Front ◽  
Time Step ◽  
The Matrix

This paper presents several algorithms that were implemented in DRUtES [1], a new open source project. DRUtES is a finite element solver for coupled nonlinear parabolic problems, namely the Richards equation with the dual porosity approach (modeling the flow of liquids in a porous medium). Mass balance consistency is crucial in any hydrological balance and contaminant transportation evaluations. An incorrect approximation of the mass term greatly depreciates the results that are obtained. An algorithm for automatic time step selection is presented, as the proper time step length is crucial for achieving accuracy of the Euler time integration method. Various problems arise with poor conditioning of the Richards equation: the computational domain is clustered into subregions separated by a wetting front, and the nonlinear constitutive functions cover a high range of values, while a very simple diagonal preconditioning method greatly improves the matrix properties. The results are presented here, together with an analysis.

A Meshfree-Based Lattice Boltzmann Approach for Simulation of Fluid Flows Within Complex Geometries

2018 ◽  
pp. 188-222 ◽  
Author(s):  
Sonam Tanwar
Keyword(s):  
Lattice Boltzmann Method ◽  
Lattice Boltzmann ◽  
Moving Boundary ◽  
Fluid Flows ◽  
Computational Domain ◽  
Time Step ◽  
Boundary Problems ◽  
Complex Process ◽  
Element Free ◽  
Boltzmann Method

This chapter develops a meshless formulation of lattice Boltzmann method for simulation of fluid flows within complex and irregular geometries. The meshless feature of proposed technique will improve the accuracy of standard lattice Boltzmann method within complicated fluid domains. Discretization of such domains itself may introduce significant numerical errors into the solution. Specifically, in phase transition or moving boundary problems, discretization of the domain is a time-consuming and complex process. In these problems, at each time step, the computational domain may change its shape and need to be re-meshed accordingly for the purpose of accuracy and stability of the solution. The author proposes to combine lattice Boltzmann method with a Galerkin meshfree technique popularly known as element-free Galerkin method in this chapter to remove the difficulties associated with traditional grid-based methods.

scholarly journals Domain Decomposition Modified with Characteristic Mixed Finite Element and Numerical Analysis for Three-Dimensional Slightly Compressible Oil-Water Seepage Displacement

2017 ◽  
Vol 9 (1) ◽  
pp. 143
Author(s):  
Yirang Yuan ◽  
Luo Chang ◽  
Changfeng Li ◽  
Tongjun Sun
Keyword(s):  
Finite Element ◽  
Domain Decomposition ◽  
Method Of Characteristics ◽  
Mixed Finite Element ◽  
Three Dimensional ◽  
Optimal Error Estimate ◽  
Computational Domain ◽  
Time Step ◽  
Slightly Compressible ◽  
The Method Of Characteristics

A parallel algorithm is presented to solve three-dimensional slightly compressible seepage displacement where domain decomposition and characteristics-mixed finite element are combined. Decomposing the computational domain into several subdomains, we define a special function to approximate the derivative at interior boundary explicitly and obtain numerical solutions of the saturation implicitly on subdomains in parallel. The method of characteristics can confirm strong stability at the fronts, and can avoid numerical dispersion and nonphysical oscillation. It can adopt large-time step but can obtain small time truncation error. So a characteristic domain decomposition finite element scheme is put forward to compute the saturation. The flow equation is computed by the method of mixed finite element and numerical accuracy of Darcy velocity is improved one order. For a model problem we apply some techniques such as variation form, domain decomposition, the method of characteristics, the principle of energy, negative norm estimates, induction hypothesis, and the theory of priori estimates of differential equations to derive optimal error estimate in $l^2$ norm. Numerical example is given to testify theoretical analysis and numerical data show that this method is effective in solving actual applications. Then it can solve the well-known problem.

Development of a Three-Dimensional Iterative Methodology Using a Commercial CFD Code for Flow Scouring Around Bridge Piers

2012 ◽  
Author(s):  
Phani Ganesh Elapolu ◽  
Pradip Majumdar ◽  
Steven A. Lottes ◽  
Milivoje Kostic
Keyword(s):  
Shear Stress ◽  
Three Dimensional ◽  
Critical Shear Stress ◽  
Shear Stresses ◽  
Computational Domain ◽  
Time Step ◽  
Mesh Morphing ◽  
River Bed ◽  
Scour Hole ◽  
The Stability

One of the major concerns affecting the safety of bridges with foundation supports in river-beds is the scouring of river-bed material from bridge supports during floods. Scour is the engineering term for the erosion caused by water around bridge elements such as piers, monopiles, or abutments. Scour holes around a monopile can jeopardize the stability of the whole structure and will require deeper piling or local armoring of the river-bed. About 500,000 bridges in the National Bridge Registry are over waterways. Many of these are considered as vulnerable to scour, about five percent are classified as scour critical, and over the last 30 years bridge failures caused by foundation scour have averaged about one every two weeks. Therefore it is of great importance to predict the correct scour development for a given bridge and flood conditions. Apart from saving time and money, integrity of bridges are important in ensuring public safety. Recent advances in computing boundary motion in combination with mesh morphing to maintain mesh quality in computational fluid dynamic analysis can be applied to predict the scour hole development, analyze the local scour phenomenon, and predict the scour hole shape and size around a pier. The main objective of the present study was to develop and implement a three dimensional iterative procedure to predict the scour hole formation around a cylindrical pier using the mesh morphing capabilities in the STARCCM+ commercial CFD code. A computational methodology has been developed using Python and Java Macros and implemented using a Bash script on a LINUX high performance computer cluster. An implicit unsteady approach was used to obtain the bed shear stresses. The mesh was iteratively deformed towards the equilibrium scour position based on the excess shear stress above the critical shear stress (supercritical shear stress). The model solves the flow field using Reynolds Averaged Navier-Stokes (RANS) equations, and the standard k–ε turbulence model. The iterative process involves stretching (morphing) a meshed domain after every time step, away from the bottom where scouring flow parameters are supercritical, and remeshing the relevant computational domain after a certain number of time steps when the morphed mesh compromises the stability of further simulation. The simulation model was validated by comparing results with limited experimental data available in the literature.

Numerical Simulation of Nonlinear Water Wave Propagation Over Rippled Bed

2007 ◽  
Author(s):  
Gerasimos A. Kolokythas ◽  
Athanassios A. Dimas
Keyword(s):  
Free Surface ◽  
Boundary Conditions ◽  
Water Waves ◽  
Stokes Equations ◽  
Surface Flow ◽  
Surface Boundary ◽  
Computational Domain ◽  
Time Step ◽  
Outflow Boundary ◽  
Rippled Bed

In the present study, numerical simulations of the free-surface flow, developing by the propagation of nonlinear water waves over a rippled bottom, are performed assuming that the corresponding flow is two-dimensional, incompressible and viscous. The simulations are based on the numerical solution of the Navier-Stokes equations subject to the fully-nonlinear free-surface boundary conditions and the suitable bottom, inflow and outflow boundary conditions. The equations are properly transformed so that the computational domain becomes time-independent. For the spatial discretization, a hybrid scheme with finite-differences and Chebyshev polynomials is applied, while a fractional time-step scheme is used for the temporal discretization. A wave absorption zone is placed at the outflow region in order to efficiently minimize reflection of waves by the outflow boundary. The numerical model is validated by comparison to the analytical solution for the laminar, oscillatory, current flow which develops a uniform boundary layer over a horizontal bottom. For the propagation of finite-amplitude waves over a rigid rippled bed, the case with wavelength to water depth ratio λ/d0 = 6 and wave height to wavelength ratio H0/λ = 0.05 is considered. The ripples have parabolic shape, while their dimensions — length and height — are chosen accordingly to fit laboratory and field data. Results indicate that the wall shear stress over the ripples and the form drag forces on the ripples increase with increasing ripple height, while the corresponding friction force is insensitive to this increase. Therefore, the percentage of friction in the total drag force decreases with increasing ripple height.

Sign in / Sign up

Export Citation Format

Share Document