Run-up heights of nearshore tsunamis based on quadtree grid system

To investigate the run-up heights of nearshore tsunamis in the vicinity of a circular island, a numerical model has been developed based on quadtree grids. The governing equations of the model are the nonlinear shallow-water equations. The governing equations are discretized explicitly by using the leap-frog scheme on adaptive hierarchical quadtree grids. The refined quadtree grids are generated around a circular island on a combined domain of rectangular and circular grids. The predicted numerical results have been verified by comparing to available laboratory measurements. A good agreement has been observed.     

植被斜坡岸滩海啸波消减数值模拟研究

An explicit one-dimensional model based on the shallow water equations(SWEs) was established in this work to simulate tsunami wave propagation on a vegetated beach. This model adopted the finite-volume method(FVM)for maintaining the mass balance of these equations. The resistance force caused by vegetation was taken into account as a source term in the momentum equation. The Harten–Lax–van Leer(HLL) approximate Riemann solver was applied to evaluate the interface fluxes for tracing the wet/dry transition boundary. This proposed model was used to simulate solitary wave run-up and long-periodic wave propagation on a sloping beach. The calibration process suitably compared the calculated results with the measured data. The tsunami waves were also simulated to discuss the water depth, tsunami force, as well as the current speed in absence of and in presence of forest domain. The results indicated that forest growth at the beach reduced wave energy loss caused by tsunamis. A series of sensitivity analyses were conducted with respect to variable parameters(such as vegetation densities, wave heights, wave periods, bed resistance, and beach slopes) to identify important influences on mitigating tsunami damage on coastal forest beach.     

Shock-capturing Boussinesq-type model for nearshore wave processes

This paper describes the formulation and validation of a nearshore wave model for tropical coastal environment. The governing Boussinesq-type equations include the conservative form of the nonlinear shallow-water equations for shock capturing. A Riemann solver supplies the inter-cell flux and bathymetry source term, while a Godunov-type scheme integrates the evolution variables in time. The model handles wave breaking through momentum conservation with energy dissipation based on an eddy viscosity concept. The computed results show very good agreement with laboratory data for wave propagation over a submerged bar, wave breaking and runup on plane beaches as well as wave transformation over fringing reefs. The model accurately describes transition between supercritical and subcritical flows as well as development of dispersive waves in the processes.     

Simulation of nearshore wave processes by a depth-integrated non-hydrostatic finite element model

This paper presents CCHE2D-NHWAVE, a depth-integrated non-hydrostatic finite element model for simulating nearshore wave processes. The governing equations are a depth-integrated vertical momentum equation and the shallow water equations including extra non-hydrostatic pressure terms, which enable the model to simulate relatively short wave motions, where both frequency dispersion and nonlinear effects play important roles. A special type of finite element method, which was previously developed for a well-validated depth-integrated free surface flow model CCHE2D, is used to solve the governing equations on a partially staggered grid using a pressure projection method. To resolve discontinuous flows, involving breaking waves and hydraulic jumps, a momentum conservation advection scheme is developed based on the partially staggered grid. In addition, a simple and efficient wetting and drying algorithm is implemented to deal with the moving shoreline. The model is first verified by analytical solutions, and then validated by a series of laboratory experiments. The comparison shows that the developed wave model without the use of any empirical parameters is capable of accurately simulating a wide range of nearshore wave processes, including propagation, breaking, and run-up of nonlinear dispersive waves and transformation and inundation of tsunami waves.     

Interaction of dispersive water waves with weakly sheared currents of arbitrary profile

A set of depth-integrated equations describing combined wave–current flows is derived and validated. To account for the effect of turbulence induced by interactions between waves and currents with arbitrary horizontal vorticity, new additional stress terms are introduced. These stresses are functions of a parameter b that relates the relative importance of wave radiation stress and bottom friction stress to the wave–current interaction. To solve the equations, a fourth-order MUSCL-TVD scheme with an approximate Riemann solver is adopted. As a first-order check of the model, the Doppler shift effect and wave dispersion over linearly sheared currents are analytically shown to be retained appropriately in the equation set. The model results are then validated through comparisons with three experimental data sets. First, based on the experiments of Kemp and Simons (1982, 1983), a reasonable functional form of b is estimated. Second, simulations examining the propagation of a weakly dispersive wave over a depth-uniform or linearly sheared current are performed. Finally, the model is applied to a more complex configuration where bichromatic waves interact with spatially varying currents. Simulated results indicate that the model is capable of predicting nearshore interactions of waves with currents of arbitrary vertical structure. One of the unique properties of the developed model is its ability to assimilate an external current field from any source, be it from a circulation model or an observation, and predict the interaction of a nonlinear and dispersive wave field with that current.     

Mathematical analysis of long-wave breaking on open channels with bottom friction

In this paper we study the breaking of long waves propagating along an open channel with linear friction on the bottom. The equations governing the wave propagation consist of a pair of first-order nonlinear hyperbolic partial differential equations (PDEs). We first transformed the PDEs into a pair of ordinary differential equations (ODEs) along the characteristic directions by means of a pair of Riemann invariants. By analyzing the ODEs, we found that the breaking of waves can be identified by the singularity of the derivative of the Riemann invariants. Thus, we derived an analytical solution for the derivative of the Riemann invariants. Then, a breaking criterion and an analytical formula for the estimation of breaking time were developed and validated through numerical experiments. It is also shown in the paper that the present model includes the previous model neglecting bottom friction as a special case.     

Benchmark computations of wave run-up on single cylinder and four cylinders by naoe-FOAM-SJTU solver

The benchmark simulations of wave run-up on a fixed single truncated circular cylinder and four circular cylinders are presented in this paper. Our in-house CFD solver naoe-FOAM-SJTU is adopted which is an unsteady two-phase CFD code based on the open source package OpenFOAM. The Navier-Stokes equations are employed as the governing equations, and the volume of fluid (VOF) method is applied for capturing the free surface. Monochromatic incident waves with the specified wave period and wave height are simulated and wave run-up heights around the cylinder are computed and recorded with numerical virtual wave probes. The relationship between the wave run-up heights and the incident wave parameters are analyzed. The numerical results indicate that the presented naoe-FOAM-SJTU solver can provide accurate predictions for the wave run-up on one fixed cylinder and four cylinders, which has been proved by the comparison of simulated results with experimental data.     

Tsunami generation,propagation, and run-up with a high-order Boussinesq model

In this work we extend a high-order Boussinesq-type (finite difference) model, capable of simulating waves out to wavenumber times depth kh < 25, to include a moving sea-bed, for the simulation of earthquake- and landslide-induced tsunamis. The extension is straight forward, requiring only an additional term within the kinematic bottom condition. As first test cases we simulate linear and nonlinear surface waves generated from both positive and negative impulsive bottom movements. The computed results compare well against earlier theoretical, numerical, and experimental values. Additionally, we show that the long-time (fully nonlinear) evolution of waves resulting from an upthrusted bottom can eventually result in true solitary waves, consistent with theoretical predictions. It is stressed, however, that the nonlinearity used far exceeds that typical of geophysical tsunamis in the open ocean. The Boussinesq-type model is then used to simulate numerous tsunami-type events generated from submerged landslides, in both one and two horizontal dimensions. The results again compare well against previous experiments and/or numerical simulations. The new extension compliments recently developed run-up capabilities within this approach, and as demonstrated, the model can therefore treat tsunami events from their initial generation, through their later propagation, and final run-up phases. The developed model is shown to maintain reasonable computational efficiency, and is therefore attractive for the simulation of such events, especially in cases where dispersion is important.     

片状分布植被对孤立波爬高影响的数值模拟研究

近岸植被对波浪爬坡具有一定的衰减作用。在自然界中,由于植被的死亡、再生或人为破坏等原因,近岸植被通常呈片状分布,且其内部分布也是不均匀的。本文以完全非线性Boussinesq方程为基础,引入植被作用项,建立了模拟近岸植被区波浪传播的数值模型,验证了模型可靠性,进而采用该模型模拟分析了片状分布植被对孤立波爬高的影响。数值模拟结果表明,片状分布植被能有效减小孤立波爬高;对于均匀分布的片状植被,高密度片状植被对孤立波爬高的消减效果优于低密度片状植被;对于相同密度、不同分布形式的片状植被,均匀分布的片状植被对孤立波的消减效果优于不均匀分布的片状植被;对于不均匀分布的片状植被,前密后疏的片状植被对孤立波的消减效果优于前疏后密的片状植被。     

基于全和谐型浅水方程的有限体积海啸数值模型开发与应用

数值模拟作为海啸预报的主要研究方法在海啸预警中起着关键作用。本文采用Godunov格式的有限体积方法,使用MUSCL-Hancock格式,并利用HLLC Riemann近似求解器计算单元界面上的流体通量,建立了球坐标系下二阶精度的海啸数值模型。模型所基于的全和谐型浅水方程保证了数值的稳定性,而地形重构方法实现了干湿边界的精准模拟。本文模拟了2015年9月16日智利M_w8.3级地震海啸,通过与智利近岸14个测站和环太平洋20个DART浮标实测数据比较,验证了模型对实际越洋海啸模拟预报的能力。     

A higher-order σ-coordinate non-hydrostatic model for nonlinear surface waves

A higher-order non-hydrostatic model in a σ-coordinate system is developed. The model uses an implicit finite difference scheme on a staggered grid to simultaneously solve the unsteady Navier-Stokes equations (NSE) with the free-surface boundary conditions. An integral method is applied to resolve the top-layer non-hydrostatic pressure, allowing for accurately resolving free-surface wave propagation. In contrast to the previous work, a higher-order spatial discretization is utilized to approximate the large horizontal pressure gradient due to steep surface waves or rapidly varying topographies. An efficient direct solver is developed to solve the resulting block hepta-diagonal matrix system. Accuracy of the new model is validated by linear and nonlinear standing waves and progressive waves. The model is then used to examine freak (extreme) waves. Features of downshifting focusing location and wave asymmetry characteristics are predicted on the temporal and spatial domains of a freak wave.     

A spectral model for unidirectional nonlinear wave propagation over arbitrary depths

A weakly-nonlinear and dispersive wave equation recently developed by the authors is used for formulating a spectral-type unidirectional wave propagation model describing spectral transformations of narrow-band waves travelling over arbitrary depths. The essential characteristics of the model equation are recapitulated first and then the spectral domain representation in terms of spatially varying harmonic amplitudes is presented. The resulting evolution equations are used to simulate the experiments concerning harmonic generation in shallow water and nonlinear random wave transformations over a submerged bar. Furthermore, the spectral model predictions are compared with the field measurements in nearshore with satisfactory results.     

A Numerical Model of Nonlinear Wave Propagation in Curvilinear Coordinates

For the simulation of the nonlinear wave propagation in coastal areas with complex boundaries,a numerical model is developed in curvilinear coordinates. In the model,the Boussinesq-type equations including the dissipation terms are employed as the governing equations. In the present model,the dependent variables of the transformed equations are the free surface elevation and the utility velocity variables,instead of the usual primitive velocity variables. The introduction of utility velocity variables which...     

The application of a Godunov-type shock capturing scheme for the simulation of waves from deep water up to the swash zone

A two-dimensional vertical (2DV) non-hydrostatic boundary fitted model based on a Godunov-type shock-capturing scheme is introduced and applied to the simulation of waves from deep water up to the swash zone. The effects of shoaling, breaking, surf zone dissipation and swash motions are considered. The application of a Godunov-type shock-capturing algorithm together with an implicit solver on a standard staggered grid is proposed as a new approach in the 2DV simulation of large gradient problems such as wave breaking and hydraulic jumps. The complete form of conservative Reynolds averaged Navier–Stokes (RANS) equations are solved using an implicit finite volume method with a pressure correction technique. The horizontal advection of the horizontal velocity is solved by an explicit predictor–corrector method. Fluxes are predicted by an exact Riemann solver and corrected by a downwind scheme. A simple total variation diminishing (TVD) method with a monotonic upstream-centered scheme for conservation laws (MUSCL) limiter function is employed to eliminate undesirable oscillations across discontinuities. Validation of the model is carried out by comparing the results of the simulations with several experimental test cases of wave breaking and run-up and the analytical solution to linear short waves in deep water. Promising performance of the model has been observed.     

Hybrid finite volume – finite difference scheme for 2DH improved Boussinesq equations

In this paper, a hybrid scheme based on a set of 2DH extended Boussinesq equations for slowly varying bathymetries is introduced. The numerical code combines the finite volume technique, applied to solve the advective part of the equations, with the finite difference method, used to discretize dispersive and source terms. Time integration is performed using the fourth-order Adams–Bashforth–Moulton predictor–corrector method; the Riemann problem is solved employing an approximate HLL solver, a fourth-order MUSCL-TVD technique is applied. Five test cases, for non-breaking and breaking waves, are reproduced to verify the model comparing its results to laboratory data or analytical solutions.     

波浪与多孔介质结构相互作用SPH模拟

应用经CSPM法和黎曼解修正后的光滑粒子流体动力学(SPH)方法,建立了主动吸收无反射数值波浪水槽,研究波浪作用下多孔介质结构的水动力特性。流体运动控制方程采用N-S方程,多孔渗水结构内流体的运动控制方程考虑渗流力的影响。数值计算结果给出了水槽内不同位置测点的波面历时曲线和越浪量随时间变化曲线,并同试验结果和Philip Liu的数值计算结果进行了比较。并对一个波浪周期内斜坡堤多孔介质结构内外的速度场和压力分布进行了讨论分析。计算分析表明,数值计算波面较Philip Liu的计算结果与试验结果吻合更好。说明应用SPH方法建立的二维数学模型能够较好地模拟破碎波在多孔渗水斜坡上的爬坡和越浪。     

Numerical study of vegetation damping effects on solitary wave run-up using the nonlinear shallow water equations

Vegetation damping effects on propagating water waves have been investigated by many researchers. This paper investigates the effects of damping due to vegetation on solitary water wave run-up via numerical simulation. The numerical model is based on an implementation of Morison's formulation for vegetation induced inertia and drag stresses in the nonlinear shallow water equations. The numerical model is solved via a finite volume method on a Cartesian cut cell mesh. The accuracy of the numerical scheme and the effects of the vegetation terms in the present model are validated by comparison with experiment results. The model is then applied to simulate a solitary wave propagating on a plane slope with vegetation. The sensitivity of solitary wave run-up to plant height, diameter and stem density is investigated by comparison of the numerical results for different patterns of vegetation. The numerical results show that vegetation can effectively reduce solitary wave propagation velocity and that solitary wave run-up is decreased with increase of plant height in water and also diameter and stem density.     

Maximum run-up,breaking conditions and dynamical forces in the swash zone: a boundary value approach

On the basis of the approximate analytical solution for the nonlinear shallow water equations of Antuono and Brocchini [M. Antuono & M. Brocchini, The boundary value problem for the nonlinear shallow water equation, Stud. Appl. Maths, 119, 71–91 (2007).], we propose useful regression curves for the prediction of maximum run-up and dynamical forces in the swash zone on a frictionless, uniformly sloping beach. For the first time the dependence of the results on both the wave height and the wave steepness is analyzed in detail providing formulae able to describe a wide class of wave inputs. Finally, the regression formulae are validated through comparison with maximum run-up laws and breaking conditions already available in the literature, the present model results appearing to better account for nonlinear effects.     

Boundary integral equation solutions for solitary wave generation, propagation and run-up

The boundary integral equation method (BIEM) is developed as a tool for studying two-dimensional, nonlinear water wave problems, including the phenomena of wave generation, propagation and run-up. The wave motions are described by a potential flow theory. Nonlinear free-surface boundary conditions are incorporated in the numerical formulation. Examples are given for either a solitary wave or two successive solitary waves. Special treatment is developed to trace the run-up and run-down along a shoreline. The accuracy of the present scheme is verified by comparing numerical results with experimental data of maximum run-up.     

Practical modified scheme of linear shallow-water equations for distant propagation of tsunamis

A simple but practical numerical model describing a distant propagation of tsunamis is newly proposed by introducing an additional term to the existing modified scheme. The numerical dispersion of the proposed model is manipulated to replace the physical dispersion of the linear Boussinesq equations without any limitation. The new model developed in this study is applied to propagation of a Gaussian hump over a constant water depth and the predicted free surface displacements are compared with available analytical solutions. A very reasonable agreement is observed.