Coupled 2D Hydrodynamic and Sediment Transport Modeling of Megaflood due to Glacier Dam-break in Altai Mountains,Southern Siberia

One of the largest known megafloods on earth resulted from a glacier dam-break,which occurred during the Late Quaternary in the Altai Mountains in Southern Siberia.Computational modeling is one of the viable approaches to enhancing the understanding of the flood events.The computational domain of this flood is over 9460 km2 and about 3.784 × 106 cells are involved as a 50 m × 50 m mesh is used,which necessitates a computationally efficient model.Here the Open MP(Open Multiprocessing) technique is adopted to parallelize the code of a coupled 2D hydrodynamic and sediment transport model.It is shown that the computational efficiency is enhanced by over 80% due to the parallelization.The floods over both fixed and mobile beds are well reproduced with specified discharge hydrographs at the dam site.Qualitatively,backwater effects during the flood are resolved at the bifurcation between the Chuja and Katun rivers.Quantitatively,the computed maximum stage and thalweg are physically consistent with the field data of the bars and deposits.The effects of sediment transport and morphological evolution on the flood are considerable.Sensitivity analyses indicate that the impact of the peak discharge is significant,whilst those of the Manningroughness,medium sediment size and shape of the inlet discharge hydrograph are marginal.     

Positivity-preserving high order well-balanced discontinuous Galerkin methods for the shallow water equations

Shallow water equations with a non-flat bottom topography have been widely used to model flows in rivers and coastal areas. An important difficulty arising in these simulations is the appearance of dry areas where no water is present, as standard numerical methods may fail in the presence of these areas. These equations also have still water steady state solutions in which the flux gradients are nonzero but exactly balanced by the source term. In this paper we propose a high order discontinuous Galerkin method which can maintain the still water steady state exactly, and at the same time can preserve the non-negativity of the water height without loss of mass conservation. A simple positivity-preserving limiter, valid under suitable CFL condition, will be introduced in one dimension and then extended to two dimensions with rectangular meshes. Numerical tests are performed to verify the positivity-preserving property, well-balanced property, high order accuracy, and good resolution for smooth and discontinuous solutions.     

High-order finite volume WENO schemes for the shallow water equations with dry states

The shallow water equations are used to model flows in rivers and coastal areas, and have wide applications in ocean, hydraulic engineering, and atmospheric modeling. These equations have still water steady state solutions in which the flux gradients are balanced by the source term. It is desirable to develop numerical methods which preserve exactly these steady state solutions. Another main difficulty usually arising from the simulation of dam breaks and flood waves flows is the appearance of dry areas where no water is present. If no special attention is paid, standard numerical methods may fail near dry/wet front and produce non-physical negative water height. A high-order accurate finite volume weighted essentially non-oscillatory (WENO) scheme is proposed in this paper to address these difficulties and to provide an efficient and robust method for solving the shallow water equations. A simple, easy-to-implement positivity-preserving limiter is introduced. One- and two-dimensional numerical examples are provided to verify the positivity-preserving property, well-balanced property, high-order accuracy, and good resolution for smooth and discontinuous solutions.     

Numerical solution of the discontinuous-bottom Shallow-water Equations with hydrostatic pressure distribution at the step

Free-surface flows are usually modelled by means of the Shallow-water Equations: this system of hyperbolic equations exhibits a source term which is proportional to the product of the water depth by the bed slope, and which takes into account the effect of gravity onto fluid mass. Recently, much attention has been paid to the case in which bottom discontinuities are present in the physical domain to be represented: in this case, it is difficult to define the non-conservative product in the distributional sense. Here, the discontinuous-bottom Shallow-water Equations with hydrostatic pressure distribution at the bed step (Bernetti et al., 2006) are discussed in the context of the theory of Dal Maso et al. (1995) [9]; finally, a first-order numerical scheme is presented, which is consistent for regular solutions, and which is able to capture contact discontinuities at bottom steps. Numerous tests are presented to show the feasibility of the scheme and its ability to converge to the exact solution in the cases of smooth as well as discontinuous bed profiles.     

Balancing the source terms in a SPH model for solving the shallow water equations

A shallow flow generally features complex hydrodynamics induced by complicated domain topography and geometry. A numerical scheme with well-balanced flux and source term gradients is therefore essential before a shallow flow model can be applied to simulate real-world problems. The issue of source term balancing has been exhaustively investigated in grid-based numerical approaches, e.g. discontinuous Galerkin finite element methods and finite volume Godunov-type methods. In recent years, a relatively new computational method, smooth particle hydrodynamics (SPH), has started to gain popularity in solving the shallow water equations (SWEs). However, the well-balanced problem has not been fully investigated and resolved in the context of SPH. This work aims to discuss the well-balanced problem caused by a standard SPH discretization to the SWEs with slope source terms and derive a corrected SPH algorithm that is able to preserve the solution of lake at rest. In order to enhance the shock capturing capability of the resulting SPH model, the Monotone Upwind-centered Scheme for Conservation Laws (MUSCL) is also explored and applied to enable Riemann solver based artificial viscosity. The new SPH model is validated against several idealized benchmark tests and a real-world dam-break case and promising results are obtained.     

Numerical resolution of well-balanced shallow water equations with complex source terms

This paper presents a well-balanced numerical scheme for simulating frictional shallow flows over complex domains involving wetting and drying. The proposed scheme solves, in a finite volume Godunov-type framework, a set of pre-balanced shallow water equations derived by considering pressure balancing. Non-negative reconstruction of Riemann states and compatible discretization of slope source term produce stable and well-balanced solutions to shallow flow hydrodynamics over complex topography. The friction source term is discretized using a splitting implicit scheme. Limiting value of the friction force is derived to ensure stability. This new numerical scheme is validated against four theoretical benchmark tests and then applied to reproduce a laboratory dam break over a domain with irregular bed profile.