有限体积KFVS方法在二维溃坝中的应用

本文采用了基于KFVS格式的有限体积方法 (FVM)求解了控制水流运动的二维浅水方程 ,建立了二维水坝瞬间溃坝的洪水演进模型 .并应用此模型模拟了二维非对称溃坝和对称溃坝情形下坝左下角有障碍物时的洪水波演进过程 .模拟结果表明该数学模型对二维浅水运动的模拟很有效 .     

动力学的流矢量分裂法模拟浅水波方程

Based on the analogy to gas dynamics, the kinetic flux vector splitting (KFVS) method is used to stimulate the shallow water wave equations. The flux vectors of the equations are split on the basis of the local equilibrium Maxwell-Boltzmann distribution. One dimensional examples including a dam breaking wave and flows over a ridge are calculated. The solutions exhibit second-order accuracy with no spurious oscillation.     

Well-balanced unstaggered central schemes for one and two-dimensional shallow water equation systems

We propose a new well-balanced unstaggered central finite volume scheme for hyperbolic balance laws with geometrical source terms. In particular we construct a new one and two-dimensional finite volume method for the numerical solution of shallow water equations on flat/variable bottom topographies. The proposed scheme evolves a non-oscillatory numerical solution on a single grid, avoids the time consuming process of solving Riemann problems arising at the cell interfaces, and is second-order accurate both in space and time. Furthermore, the numerical scheme follows a well-balanced discretization that first discretizes the geometrical source term according to the discretization of the flux terms, and then mimics the surface gradient method and discretizes the water height according to the discretization of the water level. The resulting scheme exactly satisfies the C-property at the discrete level. The proposed scheme is then applied and classical one and two-dimensional shallow water equation problems with flat or variable bottom topographies are successfully solved. The obtained numerical results are in good agreement with corresponding ones appearing in the recent literature, thus confirming the potential and efficiency of the proposed method.     

Seismic-generated unsteady motions in shallow basins and channels. Part II: Numerical modelling

In this paper, unsteady motions generated by seismic-type excitation are simulated by a 2D depth-averaged mathematical model based on the classic shallow water approximation. A suitable time-dependent forcing term is added in the governing equations, and these are solved by a MUSCL-type shock-capturing finite volume scheme with a splitting treatment of the source term. The HLL approximate Riemann solver is used to estimate the numerical fluxes. The accuracy of the numerical scheme is assessed by comparison with novel exact solutions of test cases concerning sinusoidally-generated sloshing in a prismatic tank, a rectangular open channel, and a parabolic basin. A sensitivity analysis is performed on the influence of the relevant dimensionless parameters. Moreover, numerical results are validated against experimental data available in literature concerning shallow water sloshing in a swaying tank. Finally, real‐scale applications to a reservoir created by a dam and an urban water-supply storage tank are presented. The results show that the model provides accurate solutions of the shallow water equations with a seismic-type source term and can be effectively adopted to predict the main flow features of the unsteady motion induced by horizontal seismic acceleration when the long wave assumption is valid.     

Well-balanced central finite volume methods for the Ripa system

We propose a new well-balanced central finite volume scheme for the Ripa system both in one and two space dimensions. The Ripa system is a nonhomogeneous hyperbolic system with a non-zero source term that is obtained from the shallow water equations system by incorporating horizontal temperature gradients. The proposed numerical scheme is a second-order accurate finite volume method that evolves a non-oscillatory numerical solution on a single grid, avoids the process of solving Riemann problems arising at the cell interfaces, and follows a well-balanced discretization that ensures the steady state requirement by discretizing the geometrical source term according to the discretization of the flux terms. Furthermore the proposed scheme mimics the surface gradient method and discretizes the water height according to the discretization of the water level. The proposed scheme is then applied and classical one and two-dimensional Ripa problems with flat or variable bottom topographies are successfully solved. The obtained numerical results are in good agreement with corresponding ones appearing in the recent literature, thus confirming the potential and efficiency of the proposed method.     

The Runge‐Kutta DG finite element method and the KFVS scheme for compressible flow simulations

This article presents a new type of second‐order scheme for solving the system of Euler equations, which combines the Runge‐Kutta discontinuous Galerkin (DG) finite element method and the kinetic flux vector splitting (KFVS) scheme. We first discretize the Euler equations in space with the DG method and then the resulting system from the method‐of‐lines approach will be discretized using a Runge‐Kutta method. Finally, a second‐order KFVS method is used to construct the numerical flux. The proposed scheme preserves the main advantages of the DG finite element method including its flexibility in handling irregular solution domains and in parallelization. The efficiency and effectiveness of the proposed method are illustrated by several numerical examples in one‐ and two‐dimensional spaces. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A study of splitting scheme for hyperbolic conservation laws with source terms

This study deals with the convergence of a numerical scheme for conservation laws including source terms. A splitting method for source term integration is presented. More precisely, the convergence of the numerical solution towards the entropy solution is proved in the scalar case. Because of the effect of source term, the constructed scheme is total variation bounded. Numerical experiments for one-dimensional shallow water equation are presented to demonstrate the performance of the scheme.     

求解浅水波方程的熵相容格式

提出了一种求解浅水波方程组的熵相容格式．在熵稳定通量中添加特征速度差分绝对值的项来抵消解在跨过激波时所产生的熵增,从而实现熵相容．新的数值差分格式具有形式简单、计算效率高、无需添加任何的人工数值粘性的特点．数值算例充分说明了其显著的优点．利用新格式成功地模拟了不同类型溃坝问题的激波、稀疏波传播及溃坝两侧旋涡的形成,是求解浅水波方程组较为理想的方法.     

Numerical solution of the two-dimensional shallow water equations by the application of relaxation methods

A generalization and extension of a finite difference method for calculating numerical solutions of the two dimensional shallow water system of equations is investigated. A previously developed non-oscillatory relaxation scheme is generalized as to included problems with source terms in two dimensions, with emphasis given to the bed topography, resulting to a class of methods of first- and second-order in space and time. The methods are based on classical relaxation models combined with TVD Runge–Kutta time stepping mechanisms where neither Riemann solvers nor characteristic decompositions are needed. Numerical results are presented for several test problems with or without the source term present. The wetting and drying process is also considered. The presented schemes are verified by comparing the results with documented ones.     

Robust a Simulation for Shallow Flows with Friction on Rough Topography

In this paper, we propose a robust finite volume scheme to numerically solve the shallow water equations on complex rough topography. The major difficulty of this problem is introduced by the stiff friction force term and the wet/dry interface tracking. An analytical integration method is presented for the friction force term to remove the stiffness. In the vicinity of wet/dry interface, the numerical stability can be attained by introducing an empirical parameter, the water depth tolerance, as extensively adopted in literatures. We propose a problem independent formulation for this parameter, which provides a stable scheme and preserves the overall truncation error of $\mathbb{O}$∆$x^3$. The method is applied to solve problems with complex rough topography, coupled with $h$-adaptive mesh techniques to demonstrate its robustness and efficiency.     

Seismic-generated unsteady motions in shallow basins and channels. Part I: Smooth analytical solutions

In this paper time-dependent water motions generated by seismic-type horizontal excitation in shallow basins and channels are modelled by the two-dimensional depth-averaged shallow water equations in which a specific source term is added in order to include an earthquake-induced forcing effect. Sinusoidal excitation is considered as a first approximation, and the response of shallow basins and channels to this simple external forcing is characterized. The nondimensional form of the governing equations shows that the Strouhal number and a ratio representing the amplitude of the forcing acceleration are the influential dimensionless parameters. Novel exact solutions of sinusoidally-forced smooth waves in a prismatic tank, a rectangular open channel, and a parabolic basin are presented. In the first two cases, a sway motion occurs, and reflections take place at the side walls. In the last case, the water sloshes back and forth flowing up the sloping sides of the basin; the free surface remains planar and a moving circular shoreline is present. These analytical solutions provide useful standards for assessing the accuracy of the numerical models used to solve the two-dimensional shallow water equations with source terms.     

A discontinuous Galerkin method for a new class of Green–Naghdi equations on simplicial unstructured meshes

In this paper, we introduce a discontinuous Finite Element formulation on simplicial unstructured meshes for the study of free surface flows based on the fully nonlinear and weakly dispersive Green-Naghdi equations. Working with a new class of asymptotically equivalent equations, which have a simplified analytical structure, we consider a decoupling strategy: we approximate the solutions of the classical shallow water equations supplemented with a source term globally accounting for the non-hydrostatic effects and we show that this source term can be computed through the resolution of scalar elliptic second-order sub-problems. The assets of the proposed discrete formulation are: (i) the handling of arbitrary unstructured simplicial meshes, (ii) an arbitrary order of approximation in space, (iii) the exact preservation of the motionless steady states, (iv) the preservation of the water height positivity, (v) a simple way to enhance any numerical code based on the nonlinear shallow water equations. The resulting numerical model is validated through several benchmarks involving nonlinear wave transformations and run-up over complex topographies.     

A new composite scheme for two-layer shallow water flows with shocks

This paper is devoted to solve the system of partial differential equations governing the flow of two superposed immiscible layers of shallow water flows. The system contains source terms due to bottom topography, wind stresses, and nonconservative products describing momentum exchange between the layers. The presence of these terms in the flow model forms a nonconservative system which is only conditionally hyperbolic. In addition, two-layer shallow water flows are often accompanied with moving discontinuities and shocks. Developing stable numerical methods for this class of problems presents a challenge in the field of computational hydraulics. To overcome these difficulties, a new composite scheme is proposed. The scheme consists of a time-splitting operator where in the first step the homogeneous system of the governing equations is solved using an approximate Riemann solver. In the second step a finite volume method is used to update the solution. To remove the non-physical oscillations in the vicinity of shocks a nonlinear filter is applied. The method is well-balanced, non-oscillatory and it is suitable for both low and high values of the density ratio between the two layers. Several standard test examples for two-layer shallow water flows are used to verify high accuracy and good resolution properties for smooth and discontinuous solutions.     

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.     

求解二维浅水波方程的旋转混合格式北大核心CSCD

针对二维浅水波方程数值求解问题,构造了一种旋转通量混合格式.空间方向上,该算法利用浅水波方程通量函数的旋转不变性,在单元界面法线方向及单元界面切线方向上采用可消除红斑现象的HLL与满足热力学第二定律的熵稳定加权混合数值通量函数,时间方向上采用三阶强稳定Runge-Kutta法.数值结果表明,该混合格式对于二维浅水波方程数值求解具有分辨率高的良好特性.     

一类通矢量分裂方法的保正性研究Ⅰ.显式格式

This paper is about the positivity analysis of a class of flux-vector splitting (FVS) methods for the compressible Euler equations, which include gas-kinetic Beam scheme[8], Steger-Warming FVS method[9], and Lax-Friedrichs scheme. It shows that the density and the internal energy could keep non-negative values under the CFL condition for all above three schemes once the initial gas stays in a physically realizable state. The proof of positivity is closely related to the pseudo-particle representation of FVS schemes.     

一类通矢量分裂方法的保正性研究Ⅰ.显式格式

１．引 言 近几年,流体力学数值方法的保正性一直受到人们的关注,特别是在高速流或低密度流体的数值模拟中,具有保正性的数值方法就显得更重要,因为用它们计算的密度或内能不出现负值．Ｅｉｎｆｅｌｄｔ等人［１］首先研究了Ｇｏｄｕｎｏｖ－型方法在密度较低时的特性．他们证明了Ｇｏｄｕｎｏｖ格式［２］是正守恒的,而 Ｒｏｅ的近似 Ｒｉｅｍａｎｎ方法［７］却不是．他们也证明了Ｈａｒｔｅｎ等人的近似Ｒｉｅｍａｎｎ方法［４］也是正守恒的,如果最大和最小波速的绝对值满足一定的稳定性界．最近,Ｔａｎｇ和 Ｘｕ［１１］分析了…     

Balance-characteristic scheme as applied to the shallow water equations over a rough bottom

The CABARET scheme is used for the numerical solution of the one-dimensional shallow water equations over a rough bottom. The scheme involves conservative and flux variables, whose values at a new time level are calculated by applying the characteristic properties of the shallow water equations. The scheme is verified using a series of test and model problems.     

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.