A new form of generalized Boussinesq equations for varying water depth

A new set of equations of motion for wave propagation in water with varying depth is derived in this study. The equations expressed by the velocity potentials and the wave surface elevations include first-order non-linearity of waves and have the same dispersion characteristic to the extended Boussinesq equations. Compared to the extended Boussinesq equations, the equations have only two unknown scalars and do not contain spatial derivatives with an order higher than 2. The wave equations are solved by a finite element method. Fourth-order predictor–corrector method is applied in the time integration and a damping layer is applied at the open boundary for absorbing the outgoing waves. The model is applied to several examples of wave propagation in variable water depth. The computational results are compared with experimental data and other numerical results available in literature. The comparison demonstrates that the new form of the equations is capable of calculating wave transformation from relative deep water to shallow water.     

An explicit finite difference model for simulating weakly nonlinear and weakly dispersive waves over slowly varying water depth

In this paper, a modified leap-frog finite difference (FD) scheme is developed to solve Non linear Shallow Water Equations (NSWE). By adjusting the FD mesh system and modifying the leap-frog algorithm, numerical dispersion is manipulated to mimic physical frequency dispersion for water wave propagation. The resulting numerical scheme is suitable for weakly nonlinear and weakly dispersive waves propagating over a slowly varying water depth. Numerical studies demonstrate that the results of the new numerical scheme agree well with those obtained by directly solving Boussinesq-type models for both long distance propagation, shoaling and re-fraction over a slowly varying bathymetry. Most importantly, the new algorithm is much more computationally efficient than existing Boussinesq-type models, making it an excellent alternative tool for simulating tsunami waves when the frequency dispersion needs to be considered.     

An Improved Coupling of Numerical and Physical Models for Simulating Wave Propagation

An improved coupling of numerical and physical models for simulating 2D wave propagation is developed in this paper. In the proposed model, an unstructured finite element model (FEM) based Boussinesq equations is applied for the numerical wave simulation, and a 2D piston-type wavemaker is used for the physical wave generation. An innovative scheme combining fourth-order Lagrange interpolation and Runge-Kutta scheme is described for solving the coupling equation. A Transfer function modulation method is presented to minimize the errors induced from the hydrodynamic invalidity of the coupling model and/or the mechanical capability of the wavemaker in area where nonlinearities or dispersion predominate. The overall performance and applicability of the coupling model has been experimentally validated by accounting for both regular and irregular waves and varying bathymetry. Experimental results show that the proposed numerical scheme and transfer function modulation method are efficient for the data transfer from the numerical model to the physical model up to a deterministic level.     

An Improved Coupling of Numerical and Physical Models for Simulating Wave Propagation简

An improved coupling of numerical and physical models for simulating 2D wave propagation is developed in this paper. In the proposed model, an unstructured finite element model （FEM） based Boussinesq equations is applied for the numerical wave simulation, and a 2D piston-type wavemaker is used for the physical wave generation. An innovative scheme combining fourth-order Lagrange interpolation and Runge-Kutta scheme is described for solving the coupling equation. A Transfer function modulation method is presented to minimize the errors induced from the hydrodynamic invalidity of the coupling model and/or the mechanical capability of the wavemaker in area where nonlinearities or dispersion predominate. The overall performance and applicability of the coupling model has been experimentally validated by accounting for both regular and irregular waves and varying bathymetry. Experimental results show that the proposed numerical scheme and transfer function modulation method are efficient for the data transfer from the numerical model to the physical model up to a deterministic level.     

非线性波传播的新型数值模拟模型及其实验验证

以一种新型的Boussinesq型方程为控制方程组,采用五阶Runge-Kutta-England格式离散时间积分,采用七点差分格式离散空间导数,并通过采用恰当的出流边界条件,从而建立了非线性波传播的新型数值模拟模型.通过对均匀水深水域内波浪传播的数值模拟说明,模型能较好地模拟大水深水域和强非线性波的传播.通过设置不同的入射波参数来进行潜堤地形上波浪传播的物理模型实验,并将数值解与物理模型实验结果进行了比较.     

加强的适合复杂地形的水波方程及其一维数值模型验证

在他人给出的方程的基础上,通过在其动量方程中引入含4个参数的公式,推导出了加强的适合复杂地形的水波方程,新方程的色散、变浅作用以及非线性均比原来适合复杂地形的方程有了改善:色散关系式与斯托克斯线性波的Padé(4,4)阶展开式一致;变浅作用在相对水深(波数乘水深)不大于6时与解析解符合较好;非线性在相对水深不大于1.05时保持在5%的误差之内.基于该方程,在非交错网格下建立的时间差分格式为混合4阶Adams-Bashforth-Moulton的一维数值模型,并在数值计算中利用了五对角宽带解法.数值模拟了潜堤上波浪传播变形,并将数值计算结果与实验结果进行了对比,验证了该数值模型是合理的.     

非线性波传播的新型数值模拟模型及其实验验证引入变换速度变量

以一种新型的含变换速度变量的Boussinesq型方程为控制方程组,采用五阶Runge-Kutta-England格式离散时间积分,采用七点差分格式离散空间导数,并采用恰当的出流边界条件,从而建立了非线性波传播的新型数值模拟模型.对均匀水深水域内波浪传播的数值模拟,说明在引入变换速度后进一步增大了模型的水深适用范围.对潜堤地形上波浪传播的数值模拟说明,在引入变换速度后进一步提高了模型的数值模拟精度.     

A depth-integrated model for weakly dispersive, turbulent, and rotational fluid flows

A set of weakly dispersive Boussinesq-type equations, derived to include viscosity and vorticity terms in a physically consistent manner, is presented in conservative form. The model includes the approximate effects of bottom-induced turbulence, in a depth-integrated sense, as a second-order correction. Associated with this turbulence, vertical and horizontal rotational effects are captured. While the turbulence and horizontal vorticity models are simplified, a model with known physical limitations has been derived that includes the quadratic bottom friction term commonly added in an ad hoc manner to the inviscid equations. An interesting result of this derivation is that one should take care when adding such ad hoc models; it is clear from this exercise that (1) it is not necessary to do so – the terms can be included through a consistent derivation from the viscous primitive equations – and (2) one cannot properly add the quadratic bottom friction term without also adding a number of additional terms in the integrated governing equations. To solve these equations numerically, a highly accurate and stable model is developed. The numerical method uses a fourth-order MUSCL-TVD scheme to solve the leading order (shallow water) terms. For the dispersive terms, a cell averaged finite volume method is implemented. To verify the derived equations and the numerical model, four cases of verifications are given. First, solitary wave propagation is examined as a basic, yet fundamental, test of the models ability to predict dispersive and nonlinear wave propagation with minimal numerical error. Vertical velocity distributions of spatially uniform flows are compared with existing theory to investigate the effects of the newly included horizontal vorticity terms. Other test cases include comparisons with experiments that generate strong vorticity by the change of bottom bathymetry as well as by tidal jets through inlet structures. Very reasonable agreements are observed for the four cases, and the results provide some information as to the importance of dispersion and horizontal vorticity.     

Wave Dispersion Study in the Indian Ocean-Tsunami of December 26, 2004

A numerical study which takes into account wave dispersion effects has been carried out in the Indian Ocean to reproduce the initial stage of wave propagation of the tsunami event that occurred on December 26, 2004. Three different numerical models have been used: the nonlinear shallow water (nondispersive), the nonlinear Boussinesq, and the full Navier-Stokes aided by the volume of fluid method to track the free surface. Numerical model results are compared against each other. General features of the wave propagation agreed very well in all numerical studies. However some important differences are observed in the wave patterns, i.e., the development in time of the wave front is shown to be strongly connected to the dispersion effects. Discussions and conclusions are made about the spatial and temporal distribution of the free surface reaffirming that the dispersion mechanism is important for tsunami hazard mitigation.     

Wave Dispersion Study in the Indian Ocean-Tsunami of December 26, 2004

A numerical study which takes into account wave dispersion effects has been carried out in the Indian Ocean to reproduce the initial stage of wave propagation of the tsunami event that occurred on December 26, 2004. Three different numerical models have been used: the nonlinear shallow water (nondispersive), the nonlinear Boussinesq, and the full Navier-Stokes aided by the volume of fluid method to track the free surface. Numerical model results are compared against each other. General features of the wave propagation agreed very well in all numerical studies. However some important differences are observed in the wave patterns, i.e., the development in time of the wave front is shown to be strongly connected to the dispersion effects. Discussions and conclusions are made about the spatial and temporal distribution of the free surface reaffirming that the dispersion mechanism is important for tsunami hazard mitigation.     

具有精确色散性的非线性波浪水平二维数学模型

建立具有色散性的水平二维非线性波浪方程,方程的非线性近似到了三阶。方程以波面升高和自由表面速度势表达的微分-积分型数学方程,给出方程的数值求解方法和算例,对方程积分项的处理给出了计算方法。计算结果与Boussinesq方程模型和缓坡方程模型的对应计算结果进行了对比。     

礁坪上波浪传播的数值模拟

本文基于具备间断捕捉能力的二阶全非线性Boussinesq数值模型,对规则波和随机波在礁坪地形上的传播变形进行了数值模拟。该模型采用高阶有限体积法和有限差分方法求解守恒格式的控制方程,将波浪破碎视为间断,同时采用静态重构技术处理了海岸动边界问题。重点针对礁坪上波浪传播过程中的波高空间分布和沿程衰减,礁坪上的平均水位变化,以及波浪能量频谱的移动和空间差异等典型水动力现象开展数值计算。将数值结果与实验结果对比,两者吻合情况良好,验证了模型具有良好的稳定性,具备模拟破碎波浪和海-岸动边界的能力,能较为准确地模拟波浪在礁坪地形上的传播过程中发生的各种水动力现象。     

Three-Dimensional Numerical Model for Tsunami Propagation Passing Through An Obstacle

1 .Introduction In the present numerical analysis of a tsunami ,atwo-dimensional numerical model based on non-linear shallowwater theoryis mainly used (Aburaya and Imamura ,2002 ;Imamura ,1995 ; Goto andOgawa ,1992) .Thoughthis model representstsunami hei…     

Time splitting and linear stability of the slow part of the barotropic component

In this paper, the authors analyse the stability of particular numerical schemes used in oceanic general circulation models to deal with the barotropic momentum advection term. It is shown that, when this term is integrated using time splitting, its stability properties can be drastically reduced in configurations where there exists shallow areas, where velocities become comparable to the propagation speed of external gravity waves. A simple alternative scheme with improved stability is proposed and discussed.     

Wave dispersion study for tsunami propagation in the Sea of Marmara

In this paper the aim is to investigate whether there are differences between the dispersion and non-dispersion solutions on tsunami propagation. For this purpose, two numerical models of tsunami propagation are compared. One of these numerical models is a nondispersive model that uses Saint Venant equations and the other is a dispersive model that uses Boussinesq equations. The tsunamis resulting from a submarine mass failure (SMF) which is settled at the bottom of the north eastern Sea of Marmara are examined. An analytical solution considering wave dispersion is developed for obtaining near-field tsunami amplitudes above the submarine mass failure. Numerical modeling is used at the sea surface from the common boundary called as liquid boundary with incident waves up to the coastal regions to get the tsunami amplitudes. The output of the analytical model is taken as the disturbances for the numerical method. In the numerical solutions TELEMAC-2D software system is used for both dispersive and nondispersive modeling. The results of the dispersive and nondispersive models are compared to each other. Both temporal and spatial differences in the amplitudes and wave shapes are examined. The obtained results demonstrate that there are no noticeable differences between the dispersion and non-dispersion solutions except some special cases and some special landslide velocities.     

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方程及其数值验证

通过改进二阶全非线性 Boussinesq 波浪方程中的色散项,得到了一组没有改变原方程的数学形式但适用于更大变化水深的新方程,其色散性能和变浅性能都比原方程有了很大改进,所适用的水深范围更大,能更好地描述从深水到近岸浅水处的波浪传播;并基于新方程建立了波浪数值模型,通过模拟波浪从浅水到深水的传播变形来验证新方程的有效性.     

关于Boussinesq型水波方程理论和应用研究的综述

Boussinesq型方程是研究水波传播与演化问题的重要工具之一,本文就1967-2018年常用的Boussinesq型水波方程从理论推导和数值应用两个方面进行了回顾,以期推动该类方程在海岸(海洋)工程波浪水动力方向的深入研究和应用。此类方程推导主要从欧拉方程或Laplace方程出发。在一定的非线性和缓坡假设等条件下,国内外学者建立了多个Boussinesq型水波方程,并以Stokes波的相关理论为依据,考察了这些方程在相速度、群速度、线性变浅梯度、二阶非线性、三阶非线性、波幅离散、速度沿水深分布以及和(差)频等多方面性能的精度。将Boussinesq型水波方程分为水平二维和三维两大类,并对主要Boussinesq型水波方程的特性进行了评述。进而又对适合渗透地形和存在流体分层情况下的Boussinesq型水波方程进行了简述与评论。最后对这些方程的应用进行了总结与分析。     

Propagation and run-up of nearshore tsunamis with HLLC approximate Riemann solver

A numerical model describing the propagation and run-up process of nearshore tsunamis in the vicinity of shorelines is developed based on an approximate Riemann solver. The governing equations of the model are the nonlinear shallow-water equations. The governing equations are discretized explicitly by using a finite volume method. The nonlinear terms in the momentum equations are solved with the Harten-Lax-van Leer-Contact (HLLC) approximate Riemann solver. The developed model is first applied to prediction of water motions in a parabolic basin, and propagation and subsequent run-up process of nearshore tsunamis around a circular island. Computed results are then compared with available analytical solutions and laboratory measurements. Very reasonable agreements are observed.     

高阶Boussinesq类方程数值求解及试验验证

基于二阶非线性与色散的Boussinesq类方程,采用改善的Crank-Nicolson方法对不同情况下淹没潜堤上的波浪传播进行数值模拟。高阶方程与传统、改进型的Boussinesq方程计算结果进行比较,高阶方程的计算结果与实验吻合得更好。表明该高阶Boussinesq方程能够精确预测变水深、强非线性的复杂波况,可用于实际近岸海域波浪问题的计算。