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.     

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

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

适合复杂地形的高阶Boussinesq水波方程

针对海底坡度较大(量阶为O(1))或海底曲率较大的复杂地形,建立了一个新型高阶Boussinesq水波方程.该方程可用于研究海底存在若干相互平行沙坝引起的Bragg反射问题.方程的水平速度沿水深的分布为四次多项式,色散性和变浅作用性能的精度比经典Boussinesq方程高了一阶.方程在浅水水域可以是完全非线性的.     

Extended Boussinesq equations for rapidly varying topography

We developed a new Boussinesq-type model which extends the equations of Madsen and Sørensen [1992. A new form of the Boussinesq equations with improved linear dispersion characteristics. Part 2. A slowly varying bathymetry. Coastal Engineering 18, 183-204.] by including both bottom curvature and squared bottom slope terms. Numerical experiments were conducted for wave reflection from the Booij's [1983. A note on the accuracy of the mild-slope equation. Coastal Engineering 7, 191-203] planar slope with different wave frequencies using several types of Boussinesq equations. Madsen and Sørensen's model results are accurate in the whole slopes in shallow waters, but inaccurate in intermediate water depths. Nwogu's [1993. Alternative form of Boussinesq equation for nearshore wave propagation. Journal of Waterway, Port, Coastal and Ocean Engineering 119, 618-638] model results are accurate up to 1:1 (V:H) slope, but significantly inaccurate for steep slopes. The present model results are accurate up to the slope of 1:1, but somewhat inaccurate for very steep slopes. Further, numerical experiments were conducted for wave reflections from a ripple patch and also a Gaussian-shaped trench. For the two cases, the results of Nwogu's model and the present model are accurate, because these models include the bottom curvature term which is important for the cases. However, Madsen and Sørensen's model results are inaccurate, because this model neglects the bottom curvature term.     

A composite numerical model for wave diffraction in a harbor with varying water depth

A composite numerical model is presented for computing the wave field in a harbor. The mild slope equation is discretized by a finite element method in the domain concerned. Out of the computational domain, the water depth is assumed to be constant. The boundary element method is applied to the outer boundary for dealing with the infinite boundary condition. Because the model satisfies strictly the infinite boundary condition, more accurate results can be obtained. The model is firstly applied to compute the wave diffraction in a narrow rectangular bay and the wave diffraction from a porous cylinder. The numerical results are compared with the analytical solution, experimental data and other numerical results. Good agreements are obtained. Then the model is applied to computing the wave diffraction in a square harbor with varying water depth. The effects of the water depth in the harbor and the incoming wave direction on the wave height distribution are discussed.     

A composite numerical model for wave diffraction in a harbor with varying water depth

A composite numerical model is presented for computing the wave field in a harbor. The mild slope equation is discretized by a finite element method in the domain concerned. Out of the computational domain, the water depth is assumed to be constant. The boundary element method is applied to the outer boundary for dealing with the infinite boundary condition. Because the model satisfies strictly the infinite boundary condition, more accurate results can be obtained. The model is firstly applied to compute the wave diffraction in a narrow rectangular bay and the wave diffraction from a porous cylinder. The numerical results are compared with the analytical solution, experimental data and other numerical results. Good agreements are obtained. Then the model is applied to computing the wave diffraction in a square harbor with varying water depth. The effects of the water depth in the harbor and the incoming wave direction on the wave height distribution are discussed.     

One-dimensional numerical models of higher-order Boussinesq equations with high dispersion accuracy

Abstract-Nonlinear water wave propagation passing a submerged shelf is studied experimentally andnumerically. The applicability of the wave propagation model of higher-order Boussinesq equations de-rived by Zou(2000, Ocean Engneering, 27, 557～575) is investigated. Physical experiments areconducted; three different front slopes (1:10, 1:5 and 1:2) of the shelf are set-up in the experimentand their effects on the wave propagation are investigated. Comparisons of the numerical results withtest data are made and the higher-order Boussinesq equations agree well with the measurements since thedispersion of the model is of high accuracy. The numerical results show that the good results can also beobtained for the steep-slope cases although the mild-slope assumption is employed in the derivation of thehigher-order terms in the higher-order Boussinesq equations.     

基于四阶完全非线性Boussinesq水波方程二维波浪传播数值模型

建立基于四阶完全非线性Boussinesq水波方程的二维波浪传播数值模型。采用Kennedy等提出的涡粘方法模拟波浪破碎。在矩形网格上对控制方程进行离散,采用高精度的数值格式对离散方程进行数值求解。对规则波在具有三维特征地形上的传播过程进行了数值模拟,通过数值模拟结果与实验结果的对比,对所建立的波浪传播模型进行了验证。同时,为了考察非线性对波浪传播的影响,给出和上述模型具有同阶色散性、变浅作用性能但仅具有二阶完全非线性特征的波浪模型的数值结果。通过对比两个模型的数值结果以及实验数据,讨论非线性在波浪传播过程中的作用。研究结果表明,所建立的Boussinesq水波方程在深水范围内不但具有较精确的色散性和变浅作用性能,而且具有四阶完全非线性特征,适合模拟波浪在近岸水域的非线性运动。     

A new form of higher order Boussinesq equations

On the basis of the higher order Boussinesq equations derived by the author (1999), a new form of higher order Boussinesq equations is developed through replacing the depth-averaged velocity vector by a new velocity vector in the equations in order to increase the accuracy of the linear dispersion, shoaling property and nonlinear characteristics of the equations. The dispersion of the new equations is accurate to a [4/4] Pade expansion in kh. Compared to the previous higher order Boussinesq equations, the accuracy of quadratic transfer functions is improved and the shoaling property of the equations have higher accuracy from shallow water to deep water.     

A temporally second-order accurate Godunov-type scheme for solving the extended Boussinesq equations

A numerical scheme for solving the class of extended Boussinesq equations is presented. Unlike previous schemes, where the governing equations are integrated through time using a fourth-order method, a second-order Godunov-type scheme is used thus saving storage and computational resources. The spatial derivatives are discretised using a combination of finite-volume and finite-difference methods. A fourth-order MUSCL reconstruction technique is used to compute the values at the cell interfaces for use in the local Riemann problems, whilst the bed source and dispersion terms are discretised using centred finite-differences of up to fourth-order accuracy. Numerical results show that the class of extended Boussinesq equations can be accurately solved without the need for a fourth-order time discretisation, thus improving the computational speed of Boussinesq-type numerical models. The numerical scheme has been applied to model a number of standard test cases for the extended Boussinesq equations and comparisons made to physical wave flume experiments.     

Internal generation of waves for extended Boussinesq equations

It is studied whether the mass transport or energy transport is the proper viewpoint for internally generating waves in the extended Boussinesq equations of Nwogu [J. Waterw., Port, Coastal Ocean Eng. 119 (1993) 618–638]. Numerical solutions of the Boussinesq equations with the internal generation of sinusoidal waves show that the energy transport approach yields the required wave amplitude properly while the mass transport approach yields wave amplitude different from the required one by the ratio of phase velocity to energy velocity. The waves which pass through the wave generation point do not cause any numerical distortion while the incident waves are generated. The technique of internal generation of waves shows its capability of generating nonlinear cnoidal waves as well as linear sinusoidal waves.     

Comparison between characteristics of mild slope equations and Boussinesq equations

Boussinesq-type equations and mild-slope equations are compared in terms of their basic forms and characteristics. It is concluded that linear mild-slope equations on dispersion relation are better than non-linear Boussinesq equations. In addition, Berkhoffexperiments are computed and compared by the two models, and agreement between model results and available experimental data is found to be quite reasonable, which demonstrates the two models＇ capacity to simulate wave transformation. However they can deal with different physical processes respectively, and they have their own characteristics.     

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

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

基于双层Boussinesq方程的三维波浪数值模型及其验证

Liu等给出的最高导数为2的双层Boussinesq水波方程具有较好的色散性和非线性,基于该方程建立了有限差分法的三维波浪数值模型。在矩形网格上对方程进行了空间离散,采用高阶导数近似方程中的时、空项,时间积分采用混合4阶Adams-Bashforth-Moulton的预报—校正格式。模拟了深水条件下的规则波传播过程,计算波面与解析结果吻合较好,反映出数值模型能很好地刻画波面过程及波面处的速度变化;在kh=2π条件下可较为准确获得沿水深分布的水平和垂向速度,这与理论分析结果一致。最后,利用数值模型计算了规则波在三维特征地形上的传播变形,数值结果和试验数据吻合较好;高阶非线性项会对波浪数值结果产生一定的影响,当波浪非线性增强,水深减少将产生更多的高次谐波。建立的双层Boussinesq模型对强非线性波浪的演化具有较好的模拟精度。     

A general theory of three-dimensional wave groups Part I: The formal derivation

If we know that a wave of given exceptionally large crest-to-trough height occurs at a fixed point at an instant t_o in a random wind-generated sea state, we can predict what happens with a very high probability before and after t_o in an area surrounding . The expressions of the surface displacement and velocity potential in this area are obtained in closed form. They are exact to the first order in a Stokes expansion and hold for nearly arbitrary bandwidth and solid boundary. It will be shown in Part II that these expressions represent either the evolution of a single three-dimensional wave group or the collision of two wave groups, according to the configuration of the solid boundary. The theory was developed in a series of papers starting on 1981. This paper presents the whole theory in a compact form thanks to a radical simplification of the mathematical proof.     

On the internal wave generation in Boussinesq and mild-slope equations

The literature contains empirical knowledge on whether the wave celerity or the group velocity should be used in the line source function for internal wave generation for at given set of Boussinesq or mild-slope equations. Theoretical derivations that confirm and explain these empirical findings are devised. For Boussinesq equations with, e.g. Padé[2,2]-type of dispersion relation some procedures for internal wave generation are affected by their excitation of an evanescent mode. This has some undesirable consequences, but the evanescent-mode excitation can be avoided by the use of an “internal flux boundary”.     

A note on linear dispersion and shoaling properties in extended Boussinesq equations

A set of optimum parameter α is obtained to evaluate the linear dispersion and shoaling properties in the extended Boussinesq equations of [Madsen and Sorensen, 1992 and Nwogu, 1993], and [Chen and Liu, 1995]. Optimum α values are determined to produce minimal errors in each wave property of phase velocity, group velocity, or shoaling coefficient relative to the analytical one given by the Stokes wave theory. Comparisons are made of the percent errors in phase velocity, group velocity, and shoaling coefficient produced by the Boussinesq equations with a different set of optimum α values. The case with a fixed value of α = −0.4 is also presented in the comparison. The comparisons reveal that the optimum α value tuned for a particular wave property gives in general poor results for other properties. Considering all the properties simultaneously, the fixed value of α = −0.4 may give overall accuracies in phase velocity and shoaling coefficient for all the types of Boussinesq equations selected in this study.     

Alternative forms of the higher-order Boussinesq equations: Derivations and validations

An alternative form of the Boussinesq equations is developed, creating a model which is fully nonlinear up to O(μ⁴) (μ is the ratio of water depth to wavelength) and has dispersion accurate to the Padé [4,4] approximation. No limitation is imposed on the bottom slope; the variable distance between free surface and sea bottom is accounted for by a σ-transformation. Two reduced forms of the model are also presented, which simplify O(μ⁴) terms using the assumption ε = O(μ^2/3) (ε is the ratio of wave height to water depth). These can be seen as extensions of Serre's equations, with dispersions given by the Padé [2,2] and Padé [4,4] approximations. The third-order nonlinear characteristics of these three models are discussed using Fourier analysis, and compared to other high-order formulations of the Boussinesq equations. The models are validated against experimental measurements of wave propagation over a submerged breakwater. Finally, the nonlinear evolution of wave groups along a horizontal flume is simulated and compared to experimental data in order to investigate the effects of the amplitude dispersion and the four-wave resonant interaction.     

An integral contravariant formulation of the fully non-linear Boussinesq equations

In this paper we propose an integral form of the fully non-linear Boussinesq equations in contravariant formulation, in which Christoffel symbols are avoided, in order to simulate wave transformation phenomena, wave breaking and nearshore currents in computational domains representing the complex morphology of real coastal regions. Following the approach proposed by Chen (2006), the motion equations retain the term related to the approximation to the second order of the vertical vorticity. A new Upwind Weighted Essentially Non-Oscillatory scheme for the solution of the fully non-linear Boussinesq equations on generalised curvilinear coordinate systems is proposed. The equations are rearranged in order to solve them by a high resolution hybrid finite volume–finite difference scheme. The conservative part of the above-mentioned equations, consisting of the convective terms and the terms related to the free surface elevation, is discretised by a high-order shock-capturing finite volume scheme in which an exact Riemann solver is involved; dispersive terms and the term related to the approximation to the second order of the vertical vorticity are discretised by a cell-centred finite difference scheme. The shock-capturing method makes it possible to intrinsically model the wave breaking, therefore no additional terms are needed to take into account the breaking related energy dissipation in the surf zone. The model is verified against several benchmark tests, and the results are compared with experimental, theoretical and alternative numerical solutions.     

一维全非线性Green-Naghdi水波方程数值模型

建立了求解一维全非线性Green-Naghdi水波方程的中心有限体积/有限差分混合数值格式。采用结构化网格对守恒形式的控制方程进行离散和积分,界面数值通量采用有限体积法计算,剩余项则采用中心有限差分格式求解。其中,采用中心迎风有限体积格式计算控制体界面数值通量,并结合界面变量的线性重构方法,使其在空间上具有四阶精度,通过引入静压重构技术和波浪破碎指标使模型具备处理海岸水-陆动边界及波浪破碎的能力。时间积分则采用具有总时间变差减小(Total Variation Diminishing,TVD)性质的三阶龙格-库塔法进行。应用该模型对孤立波在常水深和斜坡海岸上的传播过程及规则波跨越潜堤传播的实验进行了数值模型研究,数值计算同解析解及实验数据吻合良好。