A preliminary study on sea wave packet equations on slowly varying topography

There is a common hypothesis for the presently popular mild-slope equations that wave particle motion is irrotational. In this paper, an attempt is made to abandon the irrotational assumption and to set up new sea wave packet equations on slowly varying topography by use of the WKBJ method. To simplify the deduction, the two-dimensional shallow water equations are used to describe the sea wave particle motion in the very shallow nearshore area. The established equations can give some characteristics of wave propagation near shore.     

High resolution Godunov‐type schemes with small stencils

Higher‐order Godunov‐type schemes have to cope with the following two problems: (i) the increase in the size of the stencil that make the scheme computationally expensive, and (ii) the monotony‐preserving treatments (limiters) that must be implemented to avoid oscillations, leading to strong damping of the solution, in particular linear waves (e.g. acoustic waves). When too compressive, limiting procedures may also trigger the instability of oscillatory numerical solutions (e.g. in advection–dispersion phenomena) via the artificial amplification of the shorter modes. The present paper proposes a new approach to carry out the reconstruction. In this approach, the values of the flow variable at the edges of the computational cells are obtained directly from the reconstruction within these cells. This method is applied to the MUSCL and DPM schemes for the solution of the linear advection equation. The modified DPM scheme can capture contact discontinuities within one computational cell, even after millions of time steps at Courant numbers ranging from 1 to values as low as 10^‐4. Linear waves are subject to negligible damping. Application of the method to the DPM for one‐dimensional advection–dispersion problems shows that the numerical instability of oscillatory solutions caused by the over compressive, original DPM limiter is eliminated. One‐ and two‐dimensional shallow water simulations show an improvement over classical methods, in particular for two‐dimensional problems with strongly distorted meshes. The quality of the computational solution in the two‐dimensional case remains acceptable even for mesh aspect ratios Δx/Δy as large as 10. The method can be extend to the discretization of higher‐order PDEs, allowing third‐order space derivatives to be discretized using only two cells in space. Copyright © 2004 John Wiley & Sons, Ltd.     

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

This paper builds upon earlier work that developed and evaluated a 1D predictor–corrector time‐marching algorithm for wave equation models and extends it to 2D. Typically, the generalized wave continuity equation (GWCE) utilizes a three time‐level semi‐implicit scheme centred at k, and the momentum equation uses a two time‐level scheme centred at k+12. It has been shown that in highly non‐linear applications, the algorithm becomes unstable at even moderate Courant numbers. This work implements and analyses an implicit treatment of the non‐linear terms through the use of an iterative time‐marching algorithm in the two‐dimensional framework. Stability results show at least an eight‐fold increase in the maximum time step, depending on the domain. Studies also examined the sensitivity of the G parameter (a numerical weighting parameter in the GWCE) with results showing the greatest increase in stability occurs when 1?G/τ_max?10, a range that coincides with the recommended range to minimize errors. Convergence studies indicate an increase in temporal accuracy from first order to second order, while overall error is less than the original algorithm, even at higher time steps. Finally, a parallel implementation of the new algorithm shows that it scales well. Copyright © 2004 John Wiley & Sons, Ltd.     

Development and testing of a simple 2D finite volume model of sub‐critical shallow water flow

Many environmental applications of shallow water flow modelling can be characterized as only slowly varying and everywhere sub‐critical. A simplified finite volume model is therefore developed that is capable of describing pertinent shallow water flow processes more efficiently than the usual Godunov/ Riemann characteristics approaches. The model is tested against a number of analytical and numerical solutions to the governing equations. The model reproduces accurately flow round a circular bend, flow over topography, flow up an initially dry beach and floodwave propagation down a meandering river reach, with mass conservative solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Quasi-potential Method to Study Nonlinear Surface Shallow Water Waves

By using the potential method and the perturbation method under the condition of small amplitude and shallow water waves, we analytically get the KdV-type equation for a viscous shallow water. It indicates that for one soliton-like solution, its amplitude will decrease as it propagates away due to the viscous effects of water.     

Extension of an explicit finite volume method to large time steps (CFL>1): application to shallow water flows

In this work, the explicit first order upwind scheme is presented under a formalism that enables the extension of the methodology to large time steps. The number of cells in the stencil of the numerical scheme is related to the allowable size of the CFL number for numerical stability. It is shown how to increase both at the same time. The basic idea is proposed for a 1D scalar equation and extended to 1D and 2D non‐linear systems with source terms. The importance of the kind of grid used is highlighted and the method is outlined for irregular grids. The good quality of the results is illustrated by means of several examples including shallow water flow test cases. The bed slope source terms are involved in the method through an upwind discretization. Copyright © 2005 John Wiley & Sons, Ltd.     

High‐order boundary conditions for linearized shallow water equations with stratification,dispersion and advection

The two‐dimensional linearized shallow water equations are considered in unbounded domains with density stratification. Wave dispersion and advection effects are also taken into account. The infinite domain is truncated via a rectangular artificial boundary ??, and a high‐order open boundary condition (OBC) is imposed on ??. Then the problem is solved numerically in the finite domain bounded by ??. A recently developed boundary scheme is employed, which is based on a reformulation of the sequence of OBCs originally proposed by Higdon. The OBCs can easily be used up to any desired order. They are incorporated here in a finite difference scheme. Numerical examples are used to demonstrate the performance and advantages of the computational method, with an emphasis is on the effect of stratification. Published in 2004 by John Wiley & Sons, Ltd.     

变系数浅水波方程的精确解

均衡作用法给出了一种求非线性发展方程孤波解的有效方法.利用该方法,运用计算机符号计算,求出了变系数的一般浅水波方程的孤子解.     

SWAMP: A two‐dimensional hydrodynamic and quality modeling platform for shallow waters

Water quality two‐dimensional models are often partitioned into separate modules with separate hydraulic and biological units. In most cases this approach results in poor flexibility whenever the biological dynamics has to be adapted to a specific situation. Conversely, an integrated approach is pursued in this article, producing a two‐dimensional hydraulic‐water quality model, named Shallow Water Analysis and Modeling Program (SWAMP) designed for shallow water bodies. The major objective of the work is to create a comprehensive two‐dimensional water quality assessment tool, based on an open framework and combining easy programming of additional procedures with a user‐friendly interface. The model is based on the numerical solution of the partial differential equations describing advection‐diffusion and biological processes on a two‐dimensional rectangular finite elements mesh. The hydraulics and advection‐diffusion modules model were validated both with experimental tracer data collected at a constructed wetland site and a comparison with a commercial hydrodynamic software, showing good agreement in both cases. Moreover, the model was tested in critical conditions for mass conservation, such as time‐varying wet boundary, showing a considerable numerical robustness. In the last part of the article water quality simulations are presented, though validation data are not yet available. Nevertheless, the observed model response demonstrates general consistency with expected results and the advantages of integrating the hydraulic and quality modules. The interactive graphical user interface (GUI) is also shown to represent a simple and effective connective tool to the integrated package. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 663–687, 2002; DOI 10.1002/num.10014     

A stabilized SPH method for inviscid shallow water flows

In this paper, the smoothed particle hydrodynamics (SPH) method is applied to the solution of shallow water equations. A brief review of the method in its standard form is first described then a variational formulation using SPH interpolation is discussed. A new technique based on the Riemann solver is introduced to improve the stability of the method. This technique leads to better results. The treatment of solid boundary conditions is discussed but remains an open problem for general geometries. The dam‐break problem with a flat bed is used as a benchmark test. Copyright © 2004 John Wiley & Sons, Ltd.