An implicit time‐marching algorithm for shallow water models based on the generalized wave continuity equation

Wave equation models currently discretize the generalized wave continuity equation with a three‐time‐level scheme centered at k and the momentum equation with a two‐time‐level scheme centered at k+1/2; non‐linear terms are evaluated explicitly. However in highly non‐linear applications, the algorithm becomes unstable at even moderate Courant numbers. This paper examines an implicit treatment of the non‐linear terms using an iterative time‐marching algorithm. Depending on the domain, results from one‐dimensional experiments show up to a tenfold increase in stability and temporal accuracy. The sensitivity of stability to variations in the G‐parameter (a numerical weighting parameter in the generalized wave continuity equation) was examined; results show that the greatest increase in stability occurs with G/τ=2–50. In the one‐dimensional experiments, three different types of node spacing techniques—constant, variable, and LTEA (Localized Truncation Error Analysis)—were examined; stability is positively correlated to the uniformity of the node spacing. Lastly, a scaling analysis demonstrates that the magnitudes of the non‐linear terms are positively correlated to those that most influence stability, particularly the term containing the G‐parameter. It is evident that the new algorithm improves stability and temporal accuracy in a cost‐effective manner. Copyright © 2001 John Wiley & Sons, Ltd.     

Similarities between the quasi‐bubble and the generalized wave continuity equation solutions to the shallow water equations

Two common strategies for solving the shallow water equations in the finite element community are the generalized wave continuity equation (GWCE) reformulation and the quasi‐bubble velocity approximation. The GWCE approach has been widely analysed in the literature. In this work, the quasi‐bubble equations are analysed and comparisons are made between the quasi‐bubble approximation of the primitive form of the shallow water equations and a linear finite element approximation of the GWCE reformulation of the shallow water equations. The discrete condensed quasi‐bubble continuity equation is shown to be identical to a discrete wave equation for a specific GWCE weighting parameter value. The discrete momentum equations are slightly different due to the bubble function. In addition, the dispersion relationships are shown to be almost identical and numerical experiments confirm that the two schemes compute almost identical results. Analysis of the quasi‐bubble formulation suggests a relationship that may guide selection of the optimal GWCE weighting parameter. Copyright © 2004 John Wiley & Sons, Ltd.     

A comparison of the GWCE and mixed P − P1 formulations in finite‐element linearized shallow‐water models

The appearance of spurious pressure modes in early shallow‐water (SW) models has resulted in two common strategies in the finite element (FE) community: using mixed primitive variable and generalized wave continuity equation (GWCE) formulations of the SW equations. One FE scheme in particular, the P ? P₁ pair, combined with the primitive equations may be advantageously compared with the wave equation formulations and both schemes have similar data structures. Our focus here is on comparing these two approaches for a number of measures including stability, accuracy, efficiency, conservation properties, and consistency. The main part of the analysis centres on stability and accuracy results via Fourier‐based dispersion analyses in the context of the linear SW equations. The numerical solutions of test problems are found to be in good agreement with the analytical results. Copyright © 2011 John Wiley & Sons, Ltd.     

Semi‐Lagrangian semi‐implicit time‐splitting scheme for the shallow water equations

Time‐splitting technique applied in the context of the semi‐Lagrangian semi‐implicit method allows the use of extended time steps mainly based on physical considerations and reduces the number of numerical operations at each time step such that it is approximately proportional to the number of the points of spatial grid. To control time growth of the additional truncation errors, the standard stabilizing correction method is modified with no penalty for accuracy and efficiency of the algorithm. A linear analysis shows that constructed scheme is stable for time steps up to 2h. Numerical integrations with actual atmospheric fields of pressure and wind confirm computational efficiency, extended stability and accuracy of the proposed scheme. Copyright © 2006 John Wiley & Sons, Ltd.     

High‐order semi‐implicit time‐integrators for a triangular discontinuous Galerkin oceanic shallow water model

We extend the explicit in time high‐order triangular discontinuous Galerkin (DG) method to semi‐implicit (SI) and then apply the algorithm to the two‐dimensional oceanic shallow water equations; we implement high‐order SI time‐integrators using the backward difference formulas from orders one to six. The reason for changing the time‐integration method from explicit to SI is that explicit methods require a very small time step in order to maintain stability, especially for high‐order DG methods. Changing the time‐integration method to SI allows one to circumvent the stability criterion due to the gravity waves, which for most shallow water applications are the fastest waves in the system (the exception being supercritical flow where the Froude number is greater than one). The challenge of constructing a SI method for a DG model is that the DG machinery requires not only the standard finite element‐type area integrals, but also the finite volume‐type boundary integrals as well. These boundary integrals pose the biggest challenge in a SI discretization because they require the construction of a Riemann solver that is the true linear representation of the nonlinear Riemann problem; if this condition is not satisfied then the resulting numerical method will not be consistent with the continuous equations. In this paper we couple the SI time‐integrators with the DG method while maintaining most of the usual attributes associated with DG methods such as: high‐order accuracy (in both space and time), parallel efficiency, excellent stability, and conservation. The only property lost is that of a compact communication stencil typical of time‐explicit DG methods; implicit methods will always require a much larger communication stencil. We apply the new high‐order SI DG method to the shallow water equations and show results for many standard test cases of oceanic interest such as: standing, Kelvin and Rossby soliton waves, and the Stommel problem. The results show that the new high‐order SI DG model, that has already been shown to yield exponentially convergent solutions in space for smooth problems, results in a more efficient model than its explicit counterpart. Furthermore, for those problems where the spatial resolution is sufficiently high compared with the length scales of the flow, the capacity to use high‐order (HO) time‐integrators is a necessary complement to the employment of HO space discretizations, since the total numerical error would be otherwise dominated by the time discretization error. In fact, in the limit of increasing spatial resolution, it makes little sense to use HO spatial discretizations coupled with low‐order time discretizations. Published in 2009 by John Wiley & Sons, Ltd.     

Efficient approximation of Sparse Jacobians for time‐implicit reduced order models

This paper introduces a sparse matrix discrete interpolation method to effectively compute matrix approximations in the reduced order modeling framework. The sparse algorithm developed herein relies on the discrete empirical interpolation method and uses only samples of the nonzero entries of the matrix series. The proposed approach can approximate very large matrices, unlike the current matrix discrete empirical interpolation method, which is limited by its large computational memory requirements. The empirical interpolation indices obtained by the sparse algorithm slightly differ from the ones computed by the matrix discrete empirical interpolation method as a consequence of the singular vectors round‐off errors introduced by the economy or full singular value decomposition (SVD) algorithms when applied to the full matrix snapshots. When appropriately padded with zeros, the economy SVD factorization of the nonzero elements of the snapshots matrix is a valid economy SVD for the full snapshots matrix. Numerical experiments are performed with the 1D Burgers and 2D shallow water equations test problems where the quadratic reduced nonlinearities are computed via tensorial calculus. The sparse matrix approximation strategy is compared against five existing methods for computing reduced Jacobians: (i) matrix discrete empirical interpolation method, (ii) discrete empirical interpolation method, (iii) tensorial calculus, (iv) full Jacobian projection onto the reduced basis subspace, and (v) directional derivatives of the model along the reduced basis functions. The sparse matrix method outperforms all other algorithms. The use of traditional matrix discrete empirical interpolation method is not possible for very large dimensions because of its excessive memory requirements. Copyright © 2016 John Wiley & Sons, Ltd.     

用离散速度方法计算浅水长波方程

用离散速度法计算浅水波方程，将空气动力学方程和浅水波方程作了比较，用Nadiga提出的近平衡流动方法模拟浅水波方程的连续和间断解。计算了一维的溃坝波问题和Thacker提出的连续解问题，结果与精确解作了比较，并且计算了水流跃过障碍物的问题。     

Impact of mass lumping on gravity and Rossby waves in 2D finite‐element shallow‐water models

The goal of this study is to evaluate the effect of mass lumping on the dispersion properties of four finite‐element velocity/surface‐elevation pairs that are used to approximate the linear shallow‐water equations. For each pair, the dispersion relation, obtained using the mass lumping technique, is computed and analysed for both gravity and Rossby waves. The dispersion relations are compared with those obtained for the consistent schemes (without lumping) and the continuous case. The P₀?P₁, RT₀ and P?P₁ pairs are shown to preserve good dispersive properties when the mass matrix is lumped. Test problems to simulate fast gravity and slow Rossby waves are in good agreement with the analytical results. Copyright © 2008 John Wiley & Sons, Ltd.     

求解浅水波方程的并行物理信息神经网络算法

双曲守恒律方程是一类比较特殊的偏微分方程,其数值求解方法的研究一直是一个热点问题,一个显著特性是即使初始条件是光滑的,其解也可能会发展成间断。浅水波方程作为非线性双曲守恒律方程,由于间断解的存在,其精确求解存在很大困难。针对浅水波方程数值求解问题,本文基于PINN(Physics informed neural networks)反问题网络结构构造新的网络,构造的网络结构包括两个并行的神经网络,其中一个网络与已知状态数据(熵稳定格式加密求出)相关,另一个网络与方程本身相关。利用已知速度数据结合浅水波方程本身求解未知水深,最终通过一些数值算例验证网络的可行性。结果表明,新的网络结构可用于浅水波方程求解,利用速度数据可以较为精确地推算出水深。     

A compact numerical algorithm for solving the time‐dependent mild slope equation

The mild slope equation has been widely used to describe combined wave refraction and diffraction. In this study, a new numerical algorithm is developed to solve the time‐dependent mild slope equation in a second‐order hyperbolic form. The numerical algorithm is based on a compact and explicit finite difference method that is second‐order accurate in both time and space. The algorithm has the similar structure to the leap‐frog method but is constructed on three time levels for the second‐order time derivative term. The numerical model has the capability of simulating transient wave motion by correctly predicting the speed of wave energy propagation, which is important for the real‐time forecast of the arrival time of storm waves generated in the far field. The model is validated against analytical solution for wave shoaling and experimental data for combined wave refraction and diffraction over a submerged elliptic shoal on a slope (Coastal Eng. 1982; 6 :255). Lastly, the realistic scale Homma's island (Geophys. Mag. 1950; 21 :199) is studied with the use of various wave periods of T = 720s, T = 120 s, and T = 24 s. These wave periods correspond to long, intermediate, and short waves for the given topography, respectively. Comparisons are made between numerical results and existing analytical solutions in terms of the wave amplification around the island, which serves as the indicator for the potential wave runup. Excellent agreements are obtained. The model runs on a PC (Pentium IV 1.8GHz) and the computer capacity allows the computation of a mesh system up to 3000 × 3000, which is equivalent to about 150 × 150 waves or a large area of 540km × 540km for a wave train with the period of T = 60 s. Copyright 2004 John Wiley & Sons, Ltd.     

An unconditionally stable semi‐Lagrangian method for the spherical atmospherical shallow water equations

A semi‐implicit, semi‐Lagrangian, mixed finite difference–finite volume model for the shallow water equations on a rotating sphere is introduced and discussed. Its main features are the vectorial treatment of the momentum equation and the finite volume approach for the continuity equation. Pressure and Coriolis terms in the momentum equation and velocity in the continuity equation are treated semi‐implicitly. Moreover, a splitting technique is introduced to preserve symmetry of the numerical scheme. An alternative asymmetric scheme (without splitting) is also introduced and the efficiency of both is discussed. The model is shown to be conservative in geopotential height and unconditionally stable for 0.5≤θ≤1. Numerical experiments on two standard test problems confirm the performance of the model. Copyright © 2000 John Wiley & Sons, Ltd.     

BOUNDARY CONDITIONS IN SHALLOW WATER MODELS—AN ALTERNATIVE IMPLEMENTATION FOR FINITE ELEMENT CODES

Finite element solution of the shallow water wave equations has found increasing use by researchers and practitioners in the modelling of oceans and coastal areas. Wave equation models, most of which use equal-orderC⁰ interpolants for both the velocity and the surface elevation, do not introduce spurious oscillation modes, hence avoiding the need for artificial or numerical damping. An important question for both primitive equation and wave equation models is the interpretation of boundary conditions. Analysis of the characteristics of the governing equations shows that for most geophysical flows a single condition at each boundary is sufficient, yet there is not a consensus in the literature as to what that boundary condition must be or how it should be implemented in a finite element code. Traditionally (partly because of limited data), surface elevation is specified at open ocean boundaries while the normal flux is specified as zero at land boundaries. In most finite element wave equation models both of these boundary conditions are implemented as essential conditions. Our recent work focuses on alternative ways to numerically implement normal flow boundary conditions with an eye towards improving the mass-conserving properties of wave equation models. A unique finite element formulation using generalized functions demonstrates that boundary conditions should be implemented by treating normal fluxes as natural conditions with the flux interpreted as external to the computational domain. Results from extensive numerical experiments show that the scheme does conserve mass for all parameter values. Furthermore, convergence studies demonstrate that the algorithm is consistent, as residual errors at the boundary diminish as the grid is refined.     

An upstream flux‐splitting finite‐volume scheme for 2D shallow water equations

An upstream flux‐splitting finite‐volume (UFF) scheme is proposed for the solutions of the 2D shallow water equations. In the framework of the finite‐volume method, the artificially upstream flux vector splitting method is employed to establish the numerical flux function for the local Riemann problem. Based on this algorithm, an UFF scheme without Jacobian matrix operation is developed. The proposed scheme satisfying entropy condition is extended to be second‐order‐accurate using the MUSCL approach. The proposed UFF scheme and its second‐order extension are verified through the simulations of four shallow water problems, including the 1D idealized dam breaking, the oblique hydraulic jump, the circular dam breaking, and the dam‐break experiment with 45° bend channel. Meanwhile, the numerical performance of the UFF scheme is compared with those of three well‐known upwind schemes, namely the Osher, Roe, and HLL schemes. It is demonstrated that the proposed scheme performs remarkably well for shallow water flows. The simulated results also show that the UFF scheme has superior overall numerical performances among the schemes tested. Copyright © 2005 John Wiley & Sons, Ltd.     

An unstructured grid,three‐dimensional model based on the shallow water equations

A semi‐implicit finite difference model based on the three‐dimensional shallow water equations is modified to use unstructured grids. There are obvious advantages in using unstructured grids in problems with a complicated geometry. In this development, the concept of unstructured orthogonal grids is introduced and applied to this model. The governing differential equations are discretized by means of a semi‐implicit algorithm that is robust, stable and very efficient. The resulting model is relatively simple, conserves mass, can fit complicated boundaries and yet is sufficiently flexible to permit local mesh refinements in areas of interest. Moreover, the simulation of the flooding and drying is included in a natural and straightforward manner. These features are illustrated by a test case for studies of convergence rates and by examples of flooding on a river plain and flow in a shallow estuary. Copyright © 2000 John Wiley & Sons, Ltd.     

A Godunov‐type fractional semi‐implicit method based on staggered grid for dam‐break modeling

A semi‐implicit finite volume model based upon staggered grid is presented for solving shallow water equation. The model employs a time‐splitting scheme that uses a predictor–corrector method for the advection term. The fluxes are calculated based on a Riemann solver in the prediction step and a downwind scheme in the correction step. A simple TVD scheme is employed for shock capturing purposes in which the Minmond limiter is used for flux functions. As a consequence of using staggered grid, an ADI method is adopted for solving the discretized equations for 2‐D problems. Several 1‐D and 2‐D flows have been modeled with satisfactory results when compared with analytical and experimental test cases. The model is also capable of simulating supercritical as well as subcritical flow. Copyright © 2010 John Wiley & Sons, Ltd.     

Baroclinic stability for a family of two‐level,semi‐implicit numerical methods for the 3D shallow water equations

The baroclinic stability of a family of two time‐level, semi‐implicit schemes for the 3D hydrostatic, Boussinesq Navier–Stokes equations (i.e. the shallow water equations), which originate from the TRIM model of Casulli and Cheng (Int. J. Numer. Methods Fluids 1992; 15 :629–648), is examined in a simple 2D horizontal–vertical domain. It is demonstrated that existing mass‐conservative low‐dissipation semi‐implicit methods, which are unconditionally stable in the inviscid limit for barotropic flows, are unstable in the same limit for baroclinic flows. Such methods can be made baroclinically stable when the integrated continuity equation is discretized with a barotropically dissipative backwards Euler scheme. A general family of two‐step predictor‐corrector schemes is proposed that have better theoretical characteristics than existing single‐step schemes. Copyright © 2006 John Wiley & Sons, Ltd.     

On unconditionally positive implicit time integration for the DG scheme applied to shallow water flows

We present a new unconditionally positivity‐preserving (PP) implicit time integration method for the DG scheme applied to shallow water flows. This novel time discretization enhances the currently used PP DG schemes, because in the majority of previous work, explicit time stepping is implemented to deal with wetting and drying. However, for explicit time integration, linear stability requires very small time steps. Especially for locally refined grids, the stiff system resulting from space discretization makes implicit or partially implicit time stepping absolutely necessary. As implicit schemes require a lot of computational time solving large systems of nonlinear equations, a much larger time step is necessary to beat explicit time stepping in terms of CPU time. Unfortunately, the current PP implicit schemes are subject to time step restrictions due to a so‐called strong stability preserving constraint. In this work, we hence give a novel approach to positivity preservation including its theoretical background. The 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. In the DG context, we prove consistency of the discretization as well as a truncation error of the third order away from the wet–dry transition. Because of the proposed modification, the implicit scheme can take full advantage of larger time steps and is able to beat explicit time stepping in terms of CPU time. The performance and accuracy of this new method are studied for several classical test cases. Copyright © 2014 John Wiley & Sons, Ltd.     

Discontinuous boundary implementation for the shallow water equations

Quasi‐bubble finite element approximations to the shallow water equations are investigated focusing on implementations of the surface elevation boundary condition. We first demonstrate by numerical results that the conventional implementation of the boundary condition degrades the accuracy of the velocity solution. It is also shown that the degraded velocity leads to a critical instability if the advection term is present in the momentum equation. Then we propose an alternative implementation for the boundary condition. We refer to this alternative implementation as a discontinuous boundary (DB) implementation because it introduces at each boundary node two independent mass–flux values that result in a discontinuity at the boundary. Numerical results show that the proposed DB implementation is consistent, stabilizes the quasi‐bubble scheme, and leads to second‐order accuracy at the surface elevation specified boundary. Copyright © 2004 John Wiley & Sons, Ltd.