基于改进缓坡方程的波浪传播数值模拟

用变分原理导出考虑底坡一阶导数平方项和二阶曲率项影响的缓坡方程,对传统缓坡方程作了改进,提高波浪在海底地形变化剧烈、水深较浅时数值模拟精度。数值计算与已有实验室试验资料比较表明,该模型可以较好地模拟有剧烈变化的海底地形的波浪传播,比传统缓坡方程模型计算结果在精度上有明显提高。     

Experimental verification of horizontal two-dimensional modified mild-slope equation model

In order to verify modified mild-slope equation models in a horizontal two-dimensional space, a hydraulic experiment is made for surface wave propagation over a circular shoal on which water depth varies substantially. A horizontal two-dimensional numerical model is also constructed based on the hyperbolic equations that have been developed from the modified mild-slope equation to account for the substantial depth variation. Comparison between experimental measurements and numerical results shows that the modified mild-slope equation model is capable of producing accurate results for wave propagation in a region where water depth varies substantially, while the conventional mild-slope equation model gives large errors as the mild-slope assumption is violated.     

非结构化网格下椭圆型缓坡方程的数值求解

椭圆型缓坡方程是一种用线性波浪理论研究近岸波浪传播变形的有效波浪数学模型。非结构化网格下的有限容积法不仅对复杂边界的适应性好,还能保证迭代求解过程的守恒性。建立了非结构化网格下的椭圆型缓坡方程数值模型。在模型中采用非结构化网格下的有限容积法对椭圆型缓坡方程进行了数值离散,结合GPBiCG(m,n)算法求解离散方程。数值计算结果表明,该数值模型可有效地用于模拟近岸缓坡区域复杂边界下波浪的传播。     

Numerical implementation and sensitivity analysis of a wave energy converter in a time-dependent mild-slope equation model

Several Wave Energy Converters (abbreviated as WECs) have intensively been studied and developed during the last decade and currently small farms of WECs are getting installed. WECs in a farm are partly absorbing, partly redistributing the incident wave power. Consequently, the power absorption of each individual WEC in a farm is affected by its neighbouring WECs. The knowledge of the wave climate around the WEC is needed to predict its performance in the farm. In this paper a technique is developed to implement a single and multiple WECs based on the overtopping principle in a time-dependent mild-slope equation model. So far, the mild-slope equations have been widely used to study wave transformations around coastal and offshore structures, such as breakwaters, piles of windmills and offshore platforms. First the limitations of the WEC implementation are discussed through a sensitivity analysis. Next the developed approach is applied to study the wave height reduction behind a single WEC and a farm. The wake behind an isolated WEC is investigated for uni- and multidirectional waves; it is observed that an increase of the directional spread leads to a faster wave redistribution behind the WEC. Further the wake in the lee of multiple WECs is calculated for two different farm lay-outs, i.e. an aligned grid and a staggered grid, by adapting the performance of each WEC to its incident wave power. The evolved technique is a fast tool to find the optimal lay-out of WECs in a farm and to study the possible influence on surrounding activities in the sea.     

一种推广的缓坡方程

从流体力学基本方程出发，假定水流的涡量和垂向流速分量小量，推导出考虑非均匀水流的推广的缓坡方程，该方程中包含了△h^2h项和（△hh)^2项。在方程中引入底摩擦项、风能输入项和非线性项，其中风能输入项的推导考虑了风浪与涌浪的区别，风流情况依据青岛海洋大学的风浪成长经验关系，涌浪情况依据Snyder等人的观测结果。经过上述推广后，得到综合考虑折射、绕射、反射、非均匀水流、底摩擦损耗、风能输入及波浪非线性的推广的缓坡方程。     

四叉树网格下的椭圆型缓坡方程数值模型研究

波浪是近岸海域关键的水动力因素之一。考虑到近岸地形复杂、波浪演化显著的特点,建立了四叉树网格体系下的椭圆型缓坡方程数值模型,采用有限体积法对模型进行数值离散,应用GPBiCG(m, n)算法求解离散后的控制方程。模型中根据波浪波长布局计算网格,生成多层次四叉树网格,对复杂计算域有较好的适应性,并且在离散和方程求解中无需引入形函数、不产生复杂的交叉项,节约了存储空间和计算时间。将模型成功应用于物理模型实验及Acapulco海湾的波浪场数值模拟,结果表明该模型能够准确、高效地模拟近岸波浪场,可为近岸波浪场的模拟提供一定的理论和技术支持。     

Three dimensional application of the complementary mild-slope equation

The complementary mild-slope equation (CMSE) is a depth-integrated equation, which models refraction and diffraction of linear time-harmonic water waves. For 2D problems, it was shown to give better agreements with exact linear theory compared to other mild-slope (MS) type equations. However, no reference was given to 3D problems. In contrast to other MS-type models, the CMSE is derived in terms of a stream function vector rather than in terms of a velocity potential. For the 3D case, this complicates the governing equation and creates difficulties in formulating an adequate number of boundary conditions. In this paper, the CMSE is re-derived using Hamilton's principle from the Irrotational Green–Naghdi equations with a correction for the 3D case. A parabolic version of it is presented as well. The additional boundary conditions needed for 3D problems are constructed using the irrotationality condition. The CMSE is compared with an analytical solution and wave tank experiments for 3D problems. The results show very good agreement.     

Numerical simulation of wave-induced currents combined with parabolic mild-slope equation in curvilinear coordinates

Researches on breaking-induced currents by waves are summarized firstly in this paper. Then, a combined numerical model in orthogonal curvilinear coordinates is presented to simulate wave-induced current in areas with curved boundary or irregular coastline. The proposed wave-induced current model includes a nearshore current module established through orthogonal curvilinear transformation form of shallow water equations and a wave module based on the curvilinear parabolic approximation wave equation. The wave module actually serves as the driving force to provide the current module with required radiation stresses. The Crank-Nicolson finite difference scheme and the alternating directions implicit method are used to solve the wave and current module, respectively. The established surf zone currents model is validated by two numerical experiments about longshore currents and rip currents in basins with rip channel and breakwater. The numerical results are compared with the measured data and published numerical results.     

A time-dependent numerical model of the mild-slope equation

On the basis of the previous studies, the simplest hyperbolic mild-slope equation has been gained and the linear time-dependent numerical model for the water wave propagation has been established combined with different boundary conditions. Through computing the effective surface displacement and transforming into the real transient wave motion, related wave factors will be calculated. Compared with Lin’s model, analysis shows that calculation stability of the present model is enhanced efficiently, because the truncation errors of this model are only contributed by the dissipation terms, but those of Lin’s model are induced by the convection terms, dissipation terms and source terms. The tests show that the present model succeeds the merit in Lin’s model and the computational program is simpler, the computational time is shorter, and the computational stability is enhanced efficiently. The present model has the capability of simulating transient wave motion by correctly predicting at the speed of wave propagation, which is important for the real-time forecast of the arrival time of surface waves generated in the deep sea. The model is validated against analytical solution for wave diffraction and experimental data for combined wave refraction and diffraction over a submerged elliptic shoal on a slope. Good agreements are obtained. The model can be applied to the theory research an d engineering applications about the wave propagation in a biggish area.     

An efficient finite-difference approach to the mild-slope equation

A system of differential equations, the stationary part of which can be reduced to the elliptic mild-slope equation, is derived. The transient terms make the system of equations hyperbolic and similar to the system of equations governing nearly horizontal flow. The highly efficient ADI algorithm for the latter is used iteratively to find the stationary solution. By extracting the time-harmonic part and using a varying time step in the iterations the computational time is reduced greatly as compared with previous techniques.     

Extended mild-slope equation for random waves

A time-dependent extended mild-slope equation is derived from the elliptic equation of Chamberlain and Porter [J. Fluid Mech. 291 (1995) 393] using the Taylor series technique. Numerical tests are made on a horizontally one-dimensional case for regular waves over sloping beds and for both regular and irregular waves over a ripple patch. Numerical results prove that the proposed model gives accurate results for both regular and irregular waves over rapidly varying topography.     

A note on the accuracy of the mild-slope equation

The mild-slope equation is a vertically integrated refraction-diffraction equation, used to predict wave propagation in a region with uneven bottom. As its name indicates, it is based on the assumption of a mild bottom slope. The purpose of this paper is to examine the accuracy of this equation as a function of the bottom slope. To this end a number of numerical experiments is carried out comparing solutions of the three-dimensional wave equation with solutions of the mild-slope equation.For waves propagating parallel to the depth contours it turns out that the mild-slope equation produces accurate results even if the bottom slope is of order 1. For waves propagating normal to the depth contours the mild-slope equation is less accurate. The equation can be used for a bottom inclination up to 1:3.     

Analytic study to wave scattering by a general Homma island using the explicit modified mild-slope equation

In this paper, an exact analytic solution in terms of Taylor series to the explicit modified mild-slope equation (EMMSE) for wave scattering by a general Homma island is constructed and the convergence of the series solution is analyzed. To validate the new analytic solution, comparisons are made against the existing solutions including analytic solutions to both the long-wave equation and Helmholtz equation, approximate analytic solutions to the modified mild-slope equation, numerical solutions to the mild-slope equation and experimental solutions. Because of the use of the governing equation EMMSE together with mass-conserving matching conditions along the toe of the shoal, the present model is valid for not only waves in the whole spectrum from long waves to short waves but also bathymetries with the maximal seabed slope being as high as 4.27:1. Since the general Homma island is an extension of the original Homma island, the present solution can be very conveniently used to study the effects of bottom topography on combined refraction and diffraction. It is found that the larger the shoal size is, the more significant the wave amplification against the cylinder is.     

利用有限元方法求解双曲型缓坡方程

本文提出了一种双曲型缓坡方程的有限元计算方法 ,在建立有限元积分方程时通过在造波线处加入脉动源项来实现内部造波 ,并在开边界处利用阻尼层吸波 ,减少了在边界处由于数值处理引起的误差。数值计算结果与实测值吻合良好。本方法可用于大区域波浪场的计算中     

考虑风能输入的抛物型缓坡方程

在Radder和Kirby发展的波浪折射绕射缓坡方程抛物型模型基础上,对这种模型进行了改进,改进后的模型除可以考虑波浪传播过程中的底摩阻损耗、非线性作用外,加入了风能输入对波浪传播的影响。基于风能输入项的波浪模型数值计算结果表明,在纯风浪情况下的计算结果与传统的风浪计算方法结果一致,在波浪传播过程中由于风的作用,将导致波高比无风作用下计算的波高大。     

结合抛物型缓坡方程计算波浪辐射应力

将波浪辐射应力与抛物型缓坡方程中的待求变量联系起来,提出了一种计算辐射应力的新方法,并用有限差分法对控制方程进行了数值求解。数值结果表明这种方法精度高、编程简单、求解快速,可用于实际大区域波浪辐射应力的计算。     

椭圆型缓坡方程的一个有效的有限元解

将海绵层消波的方法用于有限元方法中,提出了缓坡方程的一种有效的有限元求解方法.在应用有限元法求解椭圆型缓坡方程时,通过在方程中加摩阻项,并在入射边界(波浪由此边界进入计算域)处使用不连续单元,将绕射势从总势中分离出来,在边界上利用海绵层进行消波处理,有效地消除了由于引用放射边界条件引起的误差和数值反射现象.     

A practical alternative to the “mild-slope” wave equation

The “mild-slope” equation which describes wave propagation in shoaling water is normally expressed in an elliptic form. The resulting computational effort involved in the solution of the boundary value problem renders the method suitable only for small sea areas. The parabolic approximation to this equation considerably reduces the computation involved but must omit the reflected wave. Hence this method is not suited to the modelling of harbour systems or areas near to sea walls where reflections are considerable. This paper expresses the “mild-slope” equation in the form of a pair of first-order equations, which constitute a hyperbolic system, without the loss of the reflected wave. A finite-difference numerical scheme is described for the efficient solution of the equations which includes boundaries of arbitrary reflecting power.     

Tide simulation using the mild-slope equation with Coriolis force and bottom friction

Since the mild-slope equation was derived by Berkhoff (1972),the researchers considered various mechanism to simplify and improve the equation,which has been widely used for coastal wave field calculation.Recently,some scholars applied the mild-slope equation in simulating the tidal motion,which proves that the equation is capable to calculate the tide in actual terrain.But in their studies,they made a lot of simplifications,and did not consider the effects of Coriolis force and bottom friction on tidal wave.In this paper,the first-order linear mild-slope equations are deduced from Kirby mild-slope equation including wave and current interaction.Then,referring to the method of wave equations’ modification,the Coriolis force and bottom friction term are considered,and the effects of which have been performed with the radial sand ridges topography.Finally,the results show that the modified mild-slope equation can be used to simulate tidal motion,and the calculations agree well with the measurements,thus the applicability and validity of the mild-slope equation on tidal simulation are further proved.     

Simulation of small-amplitude frequency-dispersive transient waves by means of the mild-slope equation

A numerical model based on the mild-slope equation is applied to reproduce the propagation of small-amplitude transient waves. The model makes use of the Fourier Transform to convert the time-dependent hyperbolic equation into a set of elliptic equations in the frequency domain. The results of two available experimental studies on tsunamis generated by landslides are used to validate the model, which appears to be able of carefully reproducing the effects of the frequency dispersion. An example application of tsunamis propagating around the Stromboli island is also presented to show the applicability of the present approach to real life scenarios. It is finally discussed how this model could be applied as support to a tsunami early warning system.