Numerical modelling of deformation of multi-directional random waves over a varying topography

Two numerical formulations of the breaking phenomenon were implemented in a numerical model for random wave propagation based on the elliptic formulation of the mild-slope equation. The randomness of the wave field was simulated based on a spectral component method, in which the 3-D spectrum is discretised in components of equal energy. One of the breaking process formulations is based on the concept of breaking each independent spectral component. The other is based on the distribution of the local amount of energy dissipated through the independent spectral components. The model based on the concept of breaking each independent spectral component produces the best estimates of the wave field, when the numerical results are compared with laboratory data.     

Generation of random waves in time-dependent extended mild-slope equations using a source function method

We develop techniques of numerical wave generation in the time-dependent extended mild-slope equations of Suh et al. [1997. Time-dependent equations for wave propagation on rapidly varying topography. Coastal Engineering 32, 91–117] and Lee et al. [2003. Extended mild-slope equation for random waves. Coastal Engineering 48, 277–287] for random waves using a source function method. Numerical results for both regular and irregular waves in one and two horizontal dimensions show that the wave heights and the frequency spectra are properly reproduced. The waves that pass through the wave generation region do not cause any numerical disturbances, showing usefulness of the source function method in avoiding re-reflection problems at the offshore boundary.     

A finite difference method for effective treatment of mild-slope wave equation subject to non-reflecting boundary conditions

A numerical model for coastal water wave motion that includes an effective method for treatment of non-reflecting boundaries is presented. The second-order one-way wave equation to approximate the non-reflecting boundary condition is found to be excellent and it ensures a very low level of reflection for waves approaching the boundary at a fairly wide range of the incidence angle. If the Newman approximation is adopted, the resulting boundary condition has a unique property to allow the free propagation of wave components along the boundary. The study is also based on a newly derived mild-slope wave equation system that can be easily made compatible to the one-way wave equation. The equation system is theoretically more accurate than the previous equations in terms of the mild-slope assumption. The finite difference method defined on a staggered grid is employed to solve the basic equations and to implement the non-reflecting boundary condition. For verification, the numerical model is then applied to three coastal water wave problems including the classical problem of plane wave diffraction by a vertical circular cylinder, the problem of combined wave diffraction and refraction over a submerged hump in the open sea, and the wave deformation around a detached breakwater. In all cases, the numerical results are demonstrated to agree very well with the relevant analytical solutions or with experimental data. It is thus concluded that the numerical model proposed in this study is effective and advantageous.     

A new kind of nonlinear mild-slope equation for combined refraction-diffraction of multifrequency waves

This paper presents the numerical solution of a new nonlinear mild-slope equation governing waves with different frequency components propagating in a region of varying water depth. There are two new nonlinear equations. The linear part of the equations is the mild-slope equation, and one of the models has the same non-linearity as the Boussinesq equations. The new equations are directly applicable to the problems of nonlinear wave-wave interactions over variable depth. The equations are first simplified with the parabolic approximation, and then solved numerically with a finite difference method. The Crank-Nicolson method is used to discretize the models. The numerical models are applied to a set of published experimental cases, which are nonlinear combined refraction-diffraction with generation of higher harmonic waves. Comparison of the results shows that the present models generally predict the measurements better than other nonlinear numerical models which have been applied to the data set.     

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.     

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 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.     

Assessment of open boundary conditions on the elliptic formulation of the mild-slope equation

The performance of open boundaries in a finite differences scheme of the elliptic mild-slope equation is assessed. The wave propagation results show that lowest order parabolic radiation boundary conditions, unlike sponge layers combined with first order radiation boundary conditions, are an efficient alternative to first order radiation boundary conditions in order to improve the accuracy of the numerical solution of the problem.     

一种推广的缓坡方程

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

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

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

Refraction and diffraction of random waves through breakwater

A physical model study of combined refraction and diffraction of waves through a breakwater gap at different incident angles was conducted. Both regular and random waves with narrow and broad frequency and direction spreading were studied. Besides the presence of a mild bottom slope in the lee of the breakwater, the distribution of wave heights across the width of a navigation channel inside the model harbor was also simulated. In addition to contributing to an understanding of the phenomenon of refraction and diffraction of random waves, the relatively complete set of data obtained can serve as a benchmark for testing of numerical models.     

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

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

Improvement on open boundaries on a time dependent numerical model of wave propagation in the nearshore region

An improvement on the simulation of outgoing waves on a time dependent numerical model for water wave propagation in the nearshore region is presented. The governing equations consist of a system of first order partial differential equations (PDEs), the equation of continuity and the equation of motion. A comparative study of first order radiation boundary conditions (BCs) and first order radiation BCs combined with sponge layers is presented for cases where outgoing waves leave the numerical domain of calculation through the open boundary. A reduction of spurious reflections from the numerical open boundaries can be obtained with an irrelevant increase in terms of computational cost.     

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

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

Transient response of harbours to long waves under resonance conditions

This paper presents the results of a study aimed at quantifying the time–response of harbour basins to long waves under resonance conditions. On the basis of numerical simulations reproducing long waves in the yacht harbour of Rome (Ostia, Italy), it shows that the results valid for periodic forcing waves, acting for an infinitely long time, as those provided by models based on elliptic equations like the Helmoltz and the mild-slope equations, can be misleading with respect to the more realistic ones that can be obtained using time-varying wave equations. Taking advantage of the similarity between the processes studied here and a simple one-dimensional resonator, a method is also proposed to roughly estimate a time–response parameter of each mode of the harbour, using results from steady-state numerical model results, commonly applied for studying harbour resonance in engineering practice. On the basis of further numerical simulations, aimed at reproducing schematic harbour layouts, the effect on resonance of the position of the entrance and of an outer harbour is studied. The results indicate that the effects of design solutions to reduce resonance, by placing the entrance at the middle of the harbour, or using the outer harbour as a resonator, can be correctly evaluated only when considering the time needed for the oscillations to fully develop.     

Mild-slope equation for water waves propagating over non-uniform currents and uneven bottoms

I~IOXThe interaction between surface waves and ambient currents and nearshore topography lies atone of the heat of morphological medelling. Accurate predictions of how wave propagates overcurrents and topography, and of the consequent erosion and dePOSition of sand on a beach or tidalflat are vital when assessing how a coastline may be affected by changing conditions.The mild-slope equation was introduced by Berkhoff (1972) as a way of approximating therefraction-diffraction of linearized s…     

Probability distribution of random wave forces in weakly nonlinear random waves

Based on the second-order random wave theory, the joint statistical distribution of the horizontal velocity and acceleration is derived using the characteristic function expansion method. From the joint distribution and the Morison equation, the theoretical distributions of drag forces, inertia forces and total random wave forces are determined. The distribution of inertia forces is Gaussian as that derived using the linear wave model, whereas the distributions of drag forces and total random forces deviate slightly from those derived utilizing the linear wave model. It is found that the distribution of wave forces depends solely on the frequency spectrum of sea waves associated with the first order approximation and the second order wave–wave interaction.     

Modelling linear wave transformation induced by dissipative structures—Random waves

The relevant theory is presented and numerical results are compared with the analytical solution for the interaction of non-breaking waves with an array of vertical porous circular cylinders on a horizontal bed. The extension to the cases of unidirectional and multidirectional waves is obtained by means of a transfer function. The influence of the mechanical properties of porous structures and wave irregularity on wave transformation is analysed. Results for unidirectional and multidirectional wave spectra are compared to those obtained for regular waves. The model presented reproduces well the analytical results and provides a tool for analysing several engineering problems.     

Long-term distribution of seabed shear stresses beneath random waves

The long-term distribution of seabed shear stresses under random waves is presented. The approach combines short-term distribution of maximum bottom shear stresses with a joint frequency table of significant wave height and peak period. An example of application is given where the long-term probability of exceeding a given level of the maximum bottom shear stress in the central North Sea is presented. The example includes estimation of the return period of the critical shear stress for sheet flow conditions, as well as the bottom shear stresses associated with the 1, 10 and 100 years return periods.