Asymptotically balanced schemes for non-homogeneous hyperbolic systems – application to the Shallow Water equations

In this work we introduce a class of balanced numerical schemes, up to second order, for the solution of general non-homogeneous hyperbolic systems of conservation laws. We give a general technique to build such schemes. We also prove that they balance up to second order a large class of steady solutions in the whole domain but some subset whose measure tends to zero as the grid size decreases to zero. We finally present an application to Shallow Water equations that exhibit the good performances of some of the schemes introduced. To cite this article: T. Chacón Rebollo et al., C. R. Acad. Sci. Paris, Ser. I 338 (2004).     

A Survey of High Order Schemes for the Shallow Water Equations

In this paper, we survey our recent work on designing high order positivity-preserving well-balanced finite difference and finite volume WENO (weighted essentially non-oscillatory) schemes, and discontinuous Galerkin finite element schemes for solving the shallow water equations with a non-flat bottom topography. These schemes are genuinely high order accurate in smooth regions for general solutions, are essentially non-oscillatory for general solutions with discontinuities, and at the same time they preserve exactly the water at rest or the more general moving water steady state solutions. A simple positivity-preserving limiter, valid under suitable CFL condition, has been introduced in one dimension and reformulated to two dimensions with triangular meshes, and we prove that the resulting schemes guarantee the positivity of the water depth.     

Local boundedness of weak solutions to the Diffusive Wave Approximation of the Shallow Water equations

In this paper we prove that weak solutions to the Diffusive Wave Approximation of the Shallow Water equations

${?}_{t}u?\mathrm{?}?({(u?z)}^{\alpha }|\mathrm{?}u{|}^{\gamma ?1}\mathrm{?}u)=f$

are locally bounded. Here, u describes the height of the water, z is a given function that represents the land elevation and f is a source term accounting for evaporation, infiltration or rainfall.     

Remarks on the derivation of the hydrostatic Euler equations

The motion of an inviscid incompressible fluid between two horizontal plates is studied in the limit when the plates are infinitesimally close. The convergence of the solutions of the Euler equations to those of their formal ‘hydrostatic’ limit can be established in the case when the initial velocity field satisfies a local Rayleigh conditions. This result, originally obtained by Grenier through weighted energy estimates based on Arnold's stability analysis of the Euler equations, is proven here by a more straightforward method even closer to Arnold's method.     

浅水方程组的数值方法

构造了浅水方程组的二阶精度的TVD格式。格式由简单的TVD Runge-Kutta型时间离散和有坡度限制的空间对称离散格式组成。数值耗散项用局部棱柱化河道流的特征变量构造。格式的主要优点是能够计算天然河道中浅水方程组的弱解并且构造简单。格式能够求出天然河道或非平底部渠道中的精确静水解。给出了渠道溃坝问题数值解与解析解的比较,验证格式精度高。实际天然河道型梯级水库溃坝的数值实验表明格式稳定,适应性强。     

A Posterior Estimate of the Particle Method Error for Shallow Water Models

The article examines the application of the particle method to numerical solution of the Cauchy problem for a quasi-linear system of first-order partial differential equations with discontinuous piecewise-constant initial conditions. A posterior estimate of the particle method error is derived for some shallow-water models.     

A monotone iterative technique for solving the bending elastic beam equations

This paper discusses the solvability of the fourth-order boundary value problem

     

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

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

Shallow water equations for power law and Bingham fluids

In this note,we provide a consistant thin layer theory for power law and Bingham incompressible fluids flowing down an inclined plane under the effect of gravity.The derivation of such equations is based on formal asymptotic expansions of solutions of Cauchy momentum equations in the shallow water scaling and in the neighbourhood of steady solutions so that we can close the average equations on the fluid height h and the total discharge rate q.     

Local existence and uniqueness for the hydrostatic Euler equations on a bounded domain

We address the question of well-posedness in spaces of analytic functions for the Cauchy problem for the hydrostatic incompressible Euler equations (inviscid primitive equations) on domains with boundary. By a suitable extension of the Cauchy-Kowalewski theorem we construct a locally in time, unique, real-analytic solution and give an explicit rate of decay of the radius of real-analyticity.     

Fast Wave Averaging for the Equatorial Shallow Water Equations

The equatorial shallow water equations in a suitable limit are shown to reduce to zonal jets as the Froude number tends to zero. This is a theorem of a singular limit with a fast variable coefficient due to the vanishing of the Coriolis force at the equator. Although it is not possible to get uniform estimates in classical Sobolev spaces (other than L²) by differentiating the system, a new method exploiting the particular structure of the fast coefficient leads to uniform estimates in slightly different functional spaces. The computation of resonances shows that fast waves may interact with a strong external forcing, introduced to mimic the effects of moisture, to create zonal jets.     

Unconditionally Stable Splitting Methods for the Shallow Water Equations

The front-tracking method for hyperbolic conservation laws is combined with operator splitting to study the shallow water equations. Furthermore, the method includes adaptive grid refinement in multidimensions and shock tracking in one dimension. The front-tracking method is unconditionally stable, but for practical computations feasible CFL numbers are moderately above unity (typically between 1 and 5). The method resolves shocks sharply and is highly efficient. The numerical technique is applied to four test cases, the first being an expanding bore with rotational symmetry. The second problem addresses the question of describing the time development of two constant water levels separated by a dam that breaks instantaneously. The third problem compares the front-tracking method with an explicit analytic solution of water waves rotating over a parabolic bottom profile. Finally, we study flow over an obstacle in one dimension.     

A reliable technique for solving the weakly singular second-kind Volterra-type integral equations

In this paper, the Adomian decomposition method and the phenomenon of the self-canceling “noise” terms are used for solving the weakly singular Volterra-type, linear and nonlinear, integral equations. The solution obtained is in the form of a convergent power series with elegantly computable terms. Comparing this scheme with many collocation-type methods that use the nonpolynomial basis shows that the present approach is effective and powerful. Many test modeling problems from mathematical physics, linear and nonlinear, are discussed to illustrate the effectiveness and the performance of the decomposition method.     

A Least‐Squares Formulation for the Shallow Water Equations with Small Viscosity

A least‐squares finite element method for the shallow water equations with viscosity parameter μ > 0 is proposed and studied. The shallow water equations are reformulated as a first order system by adding a new variable for the velocity flux. The reformulated first order system is combined with a characteristic‐based time discretization and a least squares approach. For the correct boundary treatment in the limit case μ → 0, a trace theorem is presented. For the numerical computation of the velocity, the finite element spaces introduced recently by Mardal, Tai and Winther (SIAM Journal on Numerical Analysis 40, pp. 1605–1631) are used. The degrees of freedom in these spaces can be identified with the normal and tangential components, respectively. Numerical results for some test examples are shown. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Wiener-Hopf equations technique for variational inequalities

In recent years, the theory of Wiener-Hopf equations has emerged as a novel and innovative technique for developing efficient and powerful numerical methods for solving variational inequalities and complementarity problems. In this paper, we provide an account of some of the fundamental aspects of the Wiener-Hopf equations with major emphasis on the formulation, computational algorithms, various generalizations and their applications. We also suggest some open problems for further research with sufficient information and references.     

Resolvent equations technique for variational inequalities

In this paper, we establish the equivalence between the general resolvent equations and variational inequalities. This equivalence is used to suggest and analyze a number of iterative algorithms for solving variational inclusions. We also study the convergence criteria of the iterative algorithms. Our results include several previously known results as special cases.     

The Cauchy Problem for the Fifth Order Shallow Water Equation

The local well-posedness of the Cauchy problem for the fifth order shallow water equation δtu＋αδx^5u＋βδx^3u＋rδxu＋μuδru=0,x,t∈R, is established for low regularity data in Sobolev spaces H^s（s≥-3/8） by the Fourier restriction norm method. Moreover, the global well-posedness for L^2 data follows from the local well-posedness and the conserved quantity. For data in H^s（s〉0）, the global well-posedness is also proved, where the main idea is to use the generalized bilinear estimates associated with the Fourier restriction norm method to prove that the existence time of the solution only depends on the L^2 norm of initial data.     

Interaction of Shallow Water Waves

In this paper, we consider the Riemann problem and interaction of elementary waves for a nonlinear hyperbolic system of conservation laws that arises in shallow water theory. This class of equations includes as a special case the equations of classical shallow water equations. We study the bore and dilatation waves and their properties, and show the existence and uniqueness of the solution to the Riemann problem. Towards the end, we discuss numerical results for different initial data along with all possible interactions of elementary waves. It is noticed that in contrast to the p -system, the Riemann problem is solvable for arbitrary initial data, and its solution does not contain vacuum state.