强非线性和色散性Boussinesq方程数值模型检验

采用同位网格有限差分法,建立了强非线性和色散性Boussinesq方程数值计算模型。以稳恒波Fourier近似解给定入射波边界条件,对均匀水深深水和浅水域不同非线性的行进波、缓坡地形上深水至浅水域的浅水变形波、以及缓坡和陡坡地形上的波浪水槽实验进行了数值计算,并将计算结果与解析解、解析数值解以及实验值进行了较为详细的比较,从而检验了模型的色散性、非线性以及不同底坡下非线性波的浅水变形性能。     

A GPU accelerated Boussinesq-type model for coastal waves

This study presents an efficient Boussinesq-type wave model accelerated by a single Graphics Processing Unit (GPU). The model uses the hybrid finite volume and finite difference method to solve weakly dispersive and nonlinear Boussinesq equations in the horizontal plane, enabling the model to have the shock-capturing ability to deal with breaking waves and moving shoreline properly. The code is written in CUDA C. To achieve better performance, the model uses cyclic reduction technique to solve massive tridiagonal linear systems and overlapped tiling/shared memory to reduce global memory access and enhance data reuse. Four numerical tests are conducted to validate the GPU implementation. The performance of the GPU model is evaluated by running a series of numerical simulations on two GPU platforms with different hardware configurations. Compared with the CPU version, the maximum speedup ratios for single-precision and double-precision calculations are 55.56 and 32.57, respectively.     

Velocity predictions for shoaling and breaking waves with a Boussinesq-type model

In this paper, the performance of a 1-D Boussinesq model is evaluated against laboratory data for its ability to predict surf zone velocity moments. Wave evolution over a plane beach and a complex bathymetry both extending into the surf-zone is examined for six cases. For the plane beach, these comprise two cases, a spilling and a plunging cnoidal wave. For the complex bathymetry, these comprise four cases of longer and short wavelengths (spilling and plunging breakers), with regular and irregular periodicity. The model evaluation places emphasis on parameters of the wave field that could be used for the prediction of sediment transport; orbital velocity, undertow, velocity skewness, kurtosis and asymmetry. It is found that, despite an overestimation of the depth-averaged horizontal velocity in the regular waves cases, the predicted higher order velocity moments and undertow are in good agreement with the laboratory data. A bispectral analysis demonstrates that the nonlinear transfers of energy amongst the low order harmonics are well reproduced, but energy exchanges with the higher harmonics are less well predicted. As a result, the model handles velocity moments better in the shorter wave tests than in the long wave cases where triad interactions are stronger. Of the four parameters describing wave breaking, the model behaviour is most sensitive to the critical wave front slope φ_B, especially with regard to velocity skewness and kurtosis predictions. It is also found that increasing the thickness of the surface roller for the case of plunging breakers improves the model's performance.     

Higher-order Boussinesq-type equations for interfacial waves in a two-fluid system

Interfacial waves propagating along the interface between a three-dimensional two-fluid system with a rigid upper boundary and an uneven bottom are considered. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. A set of higher-order Boussinesq-type equations in terms of the depth-averaged velocities accounting for stronger nonlinearity are derived. When the small parameter measuring frequency dispersion keeping up to lower-order and full nonlinearity are considered, the equations include the Choi and Camassa’s results (1999). The enhanced equations in terms of the depth-averaged velocities are obtained by applying the enhancement technique introduced by Madsen et al. (1991) and Schaffer and Madsen (1995a). It is noted that the equations derived from the present study include, as special cases, those obtained by Madsen and Schaffer (1998). By comparison with the dispersion relation of the linear Stokes waves, we found that the dispersion relation is more improved than Choi and Camassa’s (1999) results, and the applicable scope of water depth is deeper.     

完全非线性深水波的数值模拟

基于势流理论,并结合深水波质点运动从水面向下呈e指数衰减的特性,建立了完全非线性数值变深水槽模型,通过实时模拟活塞式造波机运动来产生波浪.采用时域高阶边界元法进行模拟,利用混合欧拉-拉格朗日方法和四阶Runge-Kutta方法追踪流体瞬时水面,应用镜像格林函数消除了水槽两个侧面的积分,在水槽末端布置人工阻尼层来消除反射...     

Two sets of higher-order Boussinesq-type equations for water waves

Based on the classical Boussinesq model by Peregrine [Peregrine, D.H., 1967. Long waves on a beach. J. Fluid Mech. 27 (4), 815–827], two parameters are introduced to improve dispersion and linear shoaling characteristics. The higher order non-linear terms are added to the modified Boussinesq equations. The non-linearity of the Boussinesq model is analyzed. A parameter related to h/L₀ is used to improve the quadratic transfer function in relatively deep water. Since the dispersion characteristic of the modified Boussinesq equations with two parameters is only equal to the second-order Padé expansion of the linear dispersion relation, further improvement is done by introducing a new velocity vector to replace the depth-averaged one in the modified Boussinesq equations. The dispersion characteristic of the further modified Boussinesq equations is accurate to the fourth-order Padé approximation of the linear dispersion relation. Compared to the modified Boussinesq equations, the accuracy of quadratic transfer functions is improved and the shoaling characteristic of the equations has higher accuracy from shallow water to deep water.     

A coupled-mode system with application to nonlinear water waves propagating in finite water depth and in variable bathymetry regions

A non-linear coupled-mode system of horizontal equations is presented, modelling the evolution of nonlinear water waves in finite depth over a general bottom topography. The vertical structure of the wave field is represented by means of a local-mode series expansion of the wave potential. This series contains the usual propagating and evanescent modes, plus two additional terms, the free-surface mode and the sloping-bottom mode, enabling to consistently treat the non-vertical end-conditions at the free-surface and the bottom boundaries. The present coupled-mode system fully accounts for the effects of non-linearity and dispersion, and the local-mode series exhibits fast convergence. Thus, a small number of modes (up to 5–6) are usually enough for precise numerical solution. In the present work, the coupled-mode system is applied to the numerical investigation of families of steady travelling wave solutions in constant depth, corresponding to a wide range of water depths, ranging from intermediate depth to shallow-water wave conditions, and its results are compared vs. Stokes and cnoidal wave theories, as well as with fully nonlinear Fourier methods. Furthermore, numerical results are presented for waves propagating over variable bathymetry regions and compared with nonlinear methods based on boundary integral formulation and experimental data, showing good agreement.     

Efficient computation method for two-dimensional nonlinear waves

I~IOWNUmerics wave tank is a ~ tOOl by which variouS nonlinear wave - ac interactionpwhleln can be treated in time domain. As the face stage, we f~ our efforts on the devious of an effita saution tee~ of fully nofor waves. aam numtried of water - wavesconnected PwhlemS, es~ in a wave tank, is alwayS comas Of PartS: wave generation, P~ and a~. The effi~ Of a nUm~ wave tankdepends not Only on the qwhty Of the wave abotion tecboaue, which allowS the lOng times~ and keePS finite tank leng…     

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

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

A theory for waves interacting with porous structures with multiple regions

This study presents an analytical solution for the problem of waves passing a submerged porous structure, using a multi-region method in the solution scheme considering the characteristics of geometry and composing materials of the porous structure. Using the flux and pressure conditions on horizontal boundaries and interfaces, the orthogonal property of wave motion within the porous layers through water depth is derived, and applied in the solution process. The flux and pressure conditions on vertical boundaries and interfaces are integrated to give a set of linear matrix equations, through which the unknown coefficients are solved. Comparisons of the present method with previous studies are preceded in verification, which suggests the validity and practicability of the present study, with a further expectation of extending our work to build a mild-slope equation over multiple-layer porous medium in the future.     

Deterministic and stochastic evolution equations for fully dispersive and weakly nonlinear waves

This paper presents a new and more accurate set of deterministic evolution equations for the propagation of fully dispersive, weakly nonlinear, irregular, multidirectional waves. The equations are derived directly from the Laplace equation with leading order nonlinearity in the surface boundary conditions. It is demonstrated that previous fully dispersive formulations from the literature have used an inconsistent linear relation between the velocity potential and the surface elevation. As a consequence these formulations are accurate only in shallow water, while nonlinear transfer of energy is significantly underestimated for larger wave numbers. In the present work we correct this inconsistency. In addition to the improved deterministic formulation, we present improved stochastic evolution equations in terms of the energy spectrum and the bispectrum for multidirectional waves. The deterministic and stochastic formulations are solved numerically for the case of cross shore motion of unidirectional waves and the results are verified against laboratory data for wave propagation over submerged bars and over a plane slope. Outside the surf zone the two model predictions are generally in good agreement with the measurements, and it is found that the accuracy of e.g., the energy spectrum and of the third-order statistics is considerably improved by the new formulations, particularly outside the shallow-water range.     

Wave propagation in a fully nonlinear numerical wave tank: A desingularized method

The problem of wave propagation in a fully nonlinear numerical wave tank is studied using desingularized boundary integral equation method coupled with mixed Eulerian–Lagrangian formulation. The present method is employed to solve the potential flow boundary value problem at each time step. The fourth-order predictor–corrector Adams–Bashforth–Moulton scheme is used for the time-stepping integration of the free surface boundary conditions. A damping layer near the end-wall of wave tank is added to absorb the outgoing waves with as little wave reflection back into the wave tank as possible. The saw-tooth instability is overcome via a five-point Chebyshev smoothing scheme. The model is applied to several wave propagations including solitary, irregular and random incident waves.     

Interactions between fully nonlinear water waves and cylinder arrays in a wave tank

Fully nonlinear interactions between water waves and vertical cylinder arrays in a numerical tank are studied based on a finite element method (FEM). The three-dimensional (3D) mesh is constructed through an extension of a 2D Delaunay surface grid along the vertical line. The velocity potential is obtained by solving a linear matrix system of FEM, and a difference scheme is then used to calculate the velocity on the free surface to track its movement. Waves and hydrodynamic forces are obtained for both bottom mounted and truncated cylinders. The simulations have provided many results to show the nature of mutual interference between cylinders in arrays and its effects on waves and forces at the nearly trapped mode frequency. The effect of the tank wall on waves and forces has been investigated, and the nonlinear features of waves and forces have also been discussed.     

完全非线性波浪与二维固定结构物作用的IGN-BEM耦合分析模型

为高效准确地对完全非线性波浪与二维固定结构物的相互作用进行模拟分析,建立了二维完全非线性时域耦合模型。耦合模型将计算域划分为靠近结构物的内域和远离结构物的外域,每个区域均采用满足完全非线性自由水面边界条件的波浪模型进行求解。在内域使用Laplace方程描述流体运动并采用高阶边界元法（BEM）对其进行求解;而在没有结构物的外域,波浪运动的控制方程为Irrotational Green-Naghdi（IGN）方程并采用有限元法（FEM）对其进行求解。内域和外域通过一段重叠区域进行耦合,从而实现模型间变量的传递。首先利用耦合模型分别对规则波的传播、直墙前立波的生成以及相关物理模型试验进行模拟,数值结果与精确解和试验结果的良好吻合验证了耦合模型耦合方式的合理性以及处理非线性问题的准确性;然后使用耦合模型模拟分析了波浪与固定结构物间的相互作用,并将结果与线性解析解以及完全非线性BEM模型的结果进行了对比分析,进一步证明了耦合模型的正确性与高效性。     

A parabolic equation extended to account for rapidly varying topography

In this paper, following the procedure outlined by Li (1994. An evolution equation for water waves. Coastal Engineering, 23, 227-242) and Hsu and Wen (2000. A study of using parabolic model to describe wave breaking and wide-angle wave incidence. Journal of the Chinese Institute of Engineers, 23(4), 515–527) and Hsu and Wen (2000) the extended refraction–diffraction equation is recasted into a time-dependent parabolic equation. This model, which includes higher-order bottom effect terms, is extended to account for a rapidly varying topography and wave energy dissipation in the surf zone. The importance of the higher-order bottom effect terms is examined in terms of the relative water depth. The present model was tested for wave reflection in a number of different environments, namely from a plane slope with different inclinations, from a patch of periodic ripples. The model was also tested for wave height distribution around a circular shoal and wave breaking on a barred beach. The comparison of predictions with other numerical models and experimental data show that the validity of the present model for describing wave propagation over a rapidly varying seabed is satisfactory.     

A coupled-mode technique for the transformation of ship-generated waves over variable bathymetry regions

In the present work, a coupled-mode technique is applied to the transformation of ship's waves over variable bathymetry regions, characterised by parallel depth-contours, without any mild-slope assumption. This method can be used, in conjunction with ship's near-field wave data in deep water or in constant-depth, as obtained by the application of modern (linearised or non-linear) ship computational fluid dynamic (CFD) codes, or experimental measurements, to support the study of wave wash generated by fast ships and its effects on the nearshore/coastal environment.

Under the assumption that the ship's track is straight and parallel to the depth-contours, and relatively far from the bottom irregularity, the problem of propagation–refraction–diffraction of ship-generated waves in a coastal environment is efficiently treated in the frequency domain, by applying the consistent coupled-mode model developed by Athanassoulis and Belibassakis [J. Fluid Mech. 1999;389] to the calculation of the transfer function enabling the pointwise transformation of ship-wave spectra over the variable bathymetry region.

Numerical results are presented for simplified ship-wave systems, obtained by the superposition of source–sink Havelock singularities simulating the basic features of the ship's wave pattern. The spatial evolution of the ship-wave system is examined over a smooth but steep shoal, resembling coastal environments, both in the subcritical and in the supercritical case. Since any ship free-wave system, either in deep water or in finite depth, can be adequately modelled by wavecut analysis and suitable distribution of Havelock singularities e.g. as presented by Scrags [21st Int. Conf. Offshore Mech. Arctic Eng., OMAE2002, Oslo, Norway, June 2002], the present method, in conjunction with ship CFD codes, supports the prediction of ship wash and its impact on coastal areas, including the effects of steep sloping-bed parts.     

适合复杂地形的高阶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.     

Transparent boundary condition for simulating nonlinear water waves by a particle method

Open boundaries are important when simulating water waves. In this study, a transparent boundary condition at an open boundary was developed for simulating nonlinear water waves propagating to a distant area using the Moving Particle Semi-implicit method. The novelty of this study is that the technique of wave analysis used in the experiment was introduced into the particle simulation to absorb incident waves; the simulation cost was reduced by employing inflow and outflow regions instead of a long dissipation region. Incident waves in front of the boundary were evaluated using Fourier analysis, and the particles on the transparent boundary were forced to move at the velocity of the analytical solution for Stokes waves in order to absorb the incident waves. The analysis was restricted to periodic waves. Wave propagation was simulated for two wave periods using the developed transparent boundary condition. The results showed that this transparent boundary transmitted the incident waves with small reflection and the simulation cost was lower than that for wave damping by a conventional highly viscous region.