Two-dimensional numerical modelling of shallow water flows over multilayer movable beds

The two-dimensional modelling of shallow water flows over multi-sediment erodible beds is presented. A novel approach is developed for the treatment of multiple sediment types in morphodynamics. The governing equations include the two-dimensional shallow water equations for hydrodynamics, an Exner-type equation for morphodynamics, a two-dimensional transport equation for the suspended sediments, and a set of empirical equations for entrainment and deposition. Multilayer sedimentary beds are formed of different erodible soils with sediment properties and new exchange conditions between the bed layers are developed for the model. The coupled equations yield a hyperbolic system of balance laws with source terms. As a numerical solver for the system, we implement a fast finite volume characteristics method. The numerical fluxes are reconstructed using the method of characteristics which employs projection techniques. The proposed finite volume solver is simple to implement, satisfies the conservation property and can be used for two-dimensional sediment transport problems in non-homogeneous isotropic beds without need of complicated three-dimensional equations. To assess the performance of the proposed models, we present numerical results for a wide variety of shallow water flows over sedimentary layers. Comparisons to experimental data for dam-break problems over movable beds are also included in this study.     

Regularized Equations for Numerical Simulation of Flows in the Two-Layer Shallow Water Approximation

Regularized equations describing hydrodynamic flows in the two-layer shallow water approximation are constructed. A conditionally stable finite-difference scheme based on the finitevolume method is proposed for the numerical solution of these equations. The scheme is tested using several well-known one-dimensional benchmark problems, including Riemann problems.     

二相流计算的一种差分算法

考虑两相流的力学行为,忽略相间的耗散作用,建立了Euler型的基本控制方程．状态方程采用刚性状态方程．基于Abgrall提出的准则,在流动区域内,对可压两相流提出了一个精度较高的Euler型数值方法,数值格式是Godunov型格式,对守恒型和非守恒型方程采用HLLC型和Lax-Friedrichs型近似Riemann解算器,引入了速度驰豫和压强驰豫过程来代替两相间的相互作用．在一维情形下给出数值算例,并且和Saurel的算例进行了比较,结果表明该算法不但精确而且稳定,且在间断处没有数值振荡．     

A robust central scheme for the shallow water flows with an abrupt topography based on modified hydrostatic reconstructions

We study a second-order central scheme for the shallow water flows with a discontinuous bottom topography based on modified hydrostatic reconstructions (HRs). The first HR scheme was proposed in Audusse et al, which may be missing the effect of the large discontinuous bottom topography. We introduce a modified HR method to cope with this numerical difficulty. The new scheme is well-balanced for still water solutions and can guarantee the positivity of the water depth. Finally, several numerical results of classical problems of the shallow water equations confirmed these properties of the new scheme. Especially, the new scheme yields superior results for the shallow water downhill flow over a step.     

A conservative extension of the method of characteristics for 1-D shallow flows

The method of characteristics (MOC) has been used for a long time in open channels and pipes flows. It is based on non-conservative equations, and hence it cannot be used directly for solving discontinuous shallow flows. In this paper we develop a conservative version of the MOC scheme for 1-D shallow flows by imposing the conservation law at the interpolation step. The conservation property of the scheme ensures the production of an accurate shock modeling and enables the MOC scheme to simulate dam-break type flows. By using a proper interpolation function, the proposed method can also produce quite accurate low-oscillatory results. A number of challenging test cases show considerable improvement compared to the traditional non-conservative MOC scheme in the case of dam-break type and trans-critical flow simulations.     

On unconditionally positivity preserving DG schemes for shallow water flows with shock capturing by adaptive filtering procedures

In this work, we present an unconditionally positivity preserving implicit time integration scheme for the DG method applied to shallow water flows. For locally refined grids with very small elements, the ODE resulting from space discretization is stiff and requires implicit or partially implicit time stepping. However, for simulations including wetting and drying or regions with small water height, severe time step restrictions have to be imposed due to positivity preservation. Nevertheless, as implicit time stepping demands a significant amount of computational time in order to solve a large system of nonlinear equations for each time step we need large time steps to obtain an efficient scheme. In this context, we propose a novel approach to the strategy of positivity preservation. This new technique is based on the so-called Patankar trick and guarantees non-negativity of the water height for any time step size while still preserving conservativity. Due to this modification, the implicit scheme can take full advantage of larger time steps and is therefore able to beat explicit time stepping in terms of CPU time. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonhydraulic Effects in Particle-Driven Gravity Currents in Deep Surroundings

In this article, we present an approach to modeling the flow of particle-driven gravity currents produced by the sudden release of well-mixed, fixed-volume suspensions into deep surroundings. Our model accounts for the initial turbulent energy of mixing in the release volume, characteristic of the classical lock–release experiments, as well as the spatiotemporal variability in the driving buoyancy forces attributable to particle settling. We show that, in contrast to compositionally driven flows, particle-driven flows cannot be described consistently in terms of shallow water theory. Specifically, we show that the presence of particles in the flow dynamics produces significant horizontal velocity shear, thereby changing the flow configuration in important ways from flows assumed to be governed by the shallow water equations. These new flow properties are calculated and contrasted with flow properties derived on the basis of the shallow water equations to show that the shallow water analysis misses dynamical features of the flow. We also show that our model provides significant improvement over the previous shallow water-based models in predicting the experimentally determined deposition patterns associated with the lock–release experiments.     

正变坐标系下完全深度平均湍流方程组

本文针对宽浅型水域，对三维湍流时均方程组逐项进行深度平均，推导出包含自由水面和地形影响的深度平均流动控制方程组．本文还同时获得了深度平均形式的ｋ－ε湍流模型方程组．因计入了水流的三维效应，该模型称为完全深度平均模型．考虑到天然水域几何边界复杂，本文运用较简便的方法，将上述模型方程组交换至正交坐标系下．所得控制方程组可以直接运用于对实际问题的数值模拟．     

A two-dimensional artificial viscosity technique for modelling discontinuity in shallow water flows

In this study, a two-dimensional cell-centred finite volume scheme is used to simulate discontinuity in shallow water flows. Instead of using a Riemann solver, an artificial viscosity technique is developed to minimise unphysical oscillations. This is constructed from a combination of a Laplacian and a biharmonic operator using a maximum eigenvalue of the Jacobian matrix. In order to achieve high-order accuracy in time, we use the fourth-order Runge–Kutta method. A hybrid formulation is then proposed to reduce computational time, in which the artificial viscosity technique is only performed once per time step. The convective flux of the shallow water equations is still re-evaluated four times, but only by averaging left and right states, thus making the computation much cheaper. A comparison of analytical and laboratory results shows that this method is highly accurate for dealing with discontinuous flows. As such, this artificial viscosity technique could become a promising method for solving the shallow water equations.     

Solving shallow-water systems in 2D domains using Finite Volume methods and multimedia SSE instructions

The goal of this paper is to construct efficient parallel solvers for 2D hyperbolic systems of conservation laws with source terms and nonconservative products. The method of lines is applied: at every intercell a projected Riemann problem along the normal direction is considered which is discretized by means of well-balanced Roe methods. The resulting 2D numerical scheme is explicit and first-order accurate. In [M.J. Castro, J.A. García, J.M. González, C. Pares, A parallel 2D Finite Volume scheme for solving systems of balance laws with nonconservative products: Application to shallow flows, Comput. Methods Appl. Mech. Engrg. 196 (2006) 2788–2815] a domain decomposition method was used to parallelize the resulting numerical scheme, which was implemented in a PC cluster by means of MPI techniques.     

Stationary Expansion Shocks for a Regularized Boussinesq System

Stationary expansion shocks have been identified recently as a new type of solution to hyperbolic conservation laws regularized by nonlocal dispersive terms that naturally arise in shallow‐water theory. These expansion shocks were studied previously for the Benjamin‐Bona‐Mahony (BBM) equation using matched asymptotic expansions. In this paper, we extend the BBM analysis to the regularized Boussinesq system by using Riemann invariants of the underlying dispersionless shallow‐water equations. The extension for a system is nontrivial, requiring a combination of small amplitude, long‐wave expansions with high order matched asymptotics. The constructed asymptotic solution is shown to be in excellent agreement with accurate numerical simulations of the Boussinesq system for a range of appropriately smoothed Riemann data.     

The structure of roll waves in two-layer flows

The structure of non-linear waves in a two-layer flow of an incompressible fluid in extended channels is investigated. Periodic discontinuous solutions, describing roll waves of finite amplitude, are constructred for the equations of two-layer shallow water. “Anomalous” waves of limited amplitude are found which correspond to the transition from stratified to slug flow conditions.     

The problem of the failure of a dam in two-layer shallow water (linear approximation)

The flows that occur in the solution, in the linear approximation, of the problem of the failure of a dam for the model of two-layer shallow water with a free boundary are analysed qualitatively. It is shown that, apart from symmetry, four basic processes and four transients exist. They are distinguished by the type of jumps in the levels on the lines of discontinuity and by the direction of the velocities in the liquid layers. Examples of the profiles of all these flows are presented.     

基于非结构自适应网格的复合有限体积法

利用文献[1]中将Lax-Wendroff格式和Lax-Friedrichs格式整体复合作用构成二维无结构网格上的复合型有限体积法,同时利用Delaunay方法,根据流场流动特性变化的梯度值为指示器对网格进行加密和粗化,实现自适应,并将此方法应用到二维浅水波方程的求解上,进行了二维部分溃坝,倾斜水跃的数值实验.结果表明,该方法是一个计算稳定、能适应复杂的求解域、能很好地捕捉激波、且计算速度快的算法.     

Relaxation approximation to bed-load sediment transport

In this work we propose and apply a numerical method based on finite volume relaxation approximation for computing the bed-load sediment transport in shallow water flows, in one and two space dimensions. The water flow is modeled by the well-known nonlinear shallow water equations which are coupled with a bed updating equation. Using a relaxation approximation, the nonlinear set of equations (and for two different formulations) is transformed to a semilinear diagonalizable problem with linear characteristic variables. A second order MUSCL-TVD method is used for the advection stage while an implicit–explicit Runge–Kutta scheme solves the relaxation stage. The main advantages of this approach are that neither Riemann problem solvers nor nonlinear iterations are required during the solution process. For the two different formulations, the applicability and effectiveness of the presented scheme is verified by comparing numerical results obtained for several benchmark test problems.     

A gas-kinetic scheme for shallow-water equations with source terms

In this paper, the Kinetic Flux Vector Splitting (KFVS) scheme is extended to solving the shallow water equations with source terms. To develop a well-balanced scheme between the source term and the flow convection, the source term effect is accounted in the flux evaluation across cell interfaces. This leads to a modified gas-kinetic scheme with particular application to the shallow water equations with bottom topography. Numerical experiments show better resolution of the unsteady solution than conventional finite difference method and KFVS method with little additional cost. Moreover, some positivity properties of the gas-kinetic scheme is established.     

Self-Similarity and the Evaluation of Long-Time Nonhydrostatic Effects in Compositionally Driven Gravity Flows in Deep Surroundings

In the study of compositionally driven gravity currents involving one or more homogeneous fluid layers, it has been customary to adopt the hydrostatic assumption for the pressure field in each layer which, in turn, leads to a depth‐independent horizontal velocity field in each of these layers and significant simplifications to the governing equations. Under this hydrostatic paradigm, each layer will then have its motion governed by the well‐known reduced dimension shallow‐water equations. For the so‐called ‐layer or reduced gravity shallow‐water equations, similarity solutions for fixed volume gravity currents released in rectangular geometry have been found. Very few attempts have been made to evaluate contributions arising from the possible loss of hydrostatic balance in the context of the problems treated using the classic shallow‐water approach. Where such attempts have been pursued, they have usually been carried out in a time‐independent context or using layer‐averaged equations and very small amplitude disturbances. The vast majority of these studies into nonhydrostatic effects do not include any relevant numerical work to assess these effects. In this paper, we develop an approach for evaluating nonhydrostatic contributions to the flow field for bottom gravity currents in deep surroundings and rectangular geometry. Our approach makes no assumptions on the amplitudes of the disturbances and does not depend on layer‐averaging in the governing equations. We seek asymptotic expansions of the solutions to the Euler equations for a shallow fluid by using the small parameter δ², where δ is the aspect ratio of the flow regime. At leading order the equations enforce hydrostatic balance while those obtained at first order retain certain nonhydrostatic effects which we evaluate. Our method for evaluation of these first‐order contributions employs the self‐similar nature of the solution to the leading‐order equations in the new first‐order equations without any vertical averaging procedures being employed.     

An efficient splitting technique for two‐layer shallow‐water model

We consider the numerical approximation of the weak solutions of the two‐layer shallow‐water equations. The model under consideration is made of two usual one‐layer shallow‐water model coupled by nonconservative products. Because of the nonconservative products of the system, which couple both one‐layer shallow‐water subsystems, the usual numerical methods have to consider the full model. Of course, uncoupled numerical techniques, just involving finite volume schemes for the basic shallow‐water equations, are very attractive since they are very easy to implement and they are costless. Recently, a stable layer splitting technique was introduced [Bouchut and Morales de Luna, M2AN Math Model Numer Anal 42 (2008), 683–698]. In the same spirit, we exhibit new splitting technique, which is proved to be well balanced and non‐negative preserving. The main benefit issuing from the here derived uncoupled method is the ability to correctly approximate the solution of very severe benchmarks. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1396–1423, 2015     

Regularized shallow water equations and an efficient method for numerical simulation of shallow water flows

Regularized shallow water equations are derived as based on a regularization of the Navier-Stokes equations in the form of quasi-gasdynamic and quasi-hydrodynamic equations. Efficient finite-difference algorithms based on the regularized shallow water equations are proposed for the numerical simulation of shallow water flows. The capabilities of the model are examined by computing a test Riemann problem, the flow over an obstacle, and asymmetric dam break.     

Mixing Closures for Conservation Laws in Stratified Flows

A closure for shocks involving the mixing of the fluids in two-layer stratified flows is proposed. The closure maximizes the rate of mixing, treating the dynamical hydraulic equations and entropy conditions as constraints. This closure may also be viewed as yielding an upper bound on the mixing rate by internal shocks. It is shown that the maximal mixing rate is accomplished by a shock moving at the fastest allowable speed against the upstream flow. Depending on whether the active constraint limiting this speed is the Lax entropy condition or the positive dissipation of energy, we distinguish precisely between internal hydraulic jumps and bores. Maximizing entrainment is shown to be equivalent to maximizing a suitable entropy associated to mixing. By using the latter, one can describe the flow globally by an optimization procedure, without treating the shocks separately. A general mathematical framework is formulated that can be applied whenever an insufficient number of conservation laws is supplemented by a maximization principle.