Wave shoaling and diffraction in current over a mild-slope

The wave relative frequency in the coordinate system moving with current and the angle between the direction of wave propagation and that of current are computed based on the wave dispersion relation. The current field is computed by solving the depth averaged shallow water equations. The wave field is computed by solving the mildslope equation which has taken the current‘s effect into account. A numerical model is established using a finite element method for simulating the wave shoaling and diffraction in current over a mild-slope, and the numerical results are reasonable to compare with the experimental data.     

An Analytical Solution for Wave Transformation Over Axi-Symmetric Topography

The present study considers wave scattering phenomena around a cylindrical island mounted on a general axisymmetric topography or a general submerged truncated axi-symmetric shoal based on the mild-slope equation. The method of separation of variables and Taylor series expansion are invoked to find the approximate solution to the variable water depth region which varies proportionally to an arbitrary power of radial distance. Validations against the solutions for the combined wave refraction and diffraction around a cylindrical island mounted on a paraboloidal shoal of Liu et al. in 2004 and the scattering and trapping of wave energy by a submerged truncated paraboloidal shoal of Lin and Liu in 2007 show excellent agreements as the power of radial distance being equal to two. For the solutions of wave refraction and diffraction around a cylindrical island mounted on a shoal with depth proportionally to an arbitrary power of radial distance, good agreements with Zhai et al.'s(2013) solutions are demonstrated. Since the robustness of the assumption of a general axi-symmetric geometry based on an arbitrary power variability of the radial distance, the present solution can be very conveniently employed to investigate the effects of bottom topography on wave scattering and trapping patterns.     

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.     

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.     

抛物型缓坡方程的变分及数值模拟

对线性水波的折射一绕射问题应用变分原理,对非等深、具有缓坡和不连续的底被导出了一种修改的抛物型缓坡方程近似模型,可预测三维地形上波浪的折射一绕射。同抛物型缓坡方程的线性方程进行了对比。通过数值模拟方法进行数值求解,表明本方法可用于地形条件下的波浪折射一绕射问题。     

A Modified Form of Mild-Slope Equation with Weakly Nonlinear Effect

Nonlinear effect is of importance to waves propagating from deep water to shallow water.Thenon-linearity of waves is widely discussed due to its high precision in application.But there are still someproblems in dealing with the nonlinear waves in practice.In this paper,a modified form of mild-slope equa-tion with weakly nonlinear effect is derived by use of the nonlinear dispersion relation and the steady mild-slope equation containing energy dissipation.The modified form of mild-slope equation is convenient to solvenonlinear effect of waves.The model is tested against the laboratory measurement for the case of a submergedelliptical shoal on a slope beach given by Berkhoff et al,The present numerical results are also comparedwith those obtained through linear wave theory.Better agreement is obtained as the modified mild-slope e-quation is employed.And the modified mild-slope equation can reasonably simulate the weakly nonlinear ef-fect of wave propagation from deep water to coast.     

Unsteady ship waves in shallow water of varying depth based on Green–Naghdi equation

Based on Green–Naghdi equation this work studies unsteady ship waves in shallow water of varying depth. A moving ship is regarded as a moving pressure disturbance on free surface. The moving pressure is incorporated into the Green–Naghdi equation to formulate forcing of ship waves in shallow water. The frequency dispersion term of the Green–Naghdi equation accounts for the effects of finite water depth on ship waves. A wave equation model and the finite element method (WE/FEM) are adopted to solve the Green–Naghdi equation. The numerical examples of a Series 60 (C_B=0.6) ship moving in shallow water are presented. Three-dimensional ship wave profiles and wave resistance are given when the ship moves in shallow water with a bed bump (or a trench). The numerical results indicate that the wave resistance increases first, then decreases, and finally returns to normal value as the ship passes a bed bump. A comparison between the numerical results predicted by the Green–Naghdi equation and the shallow water equations is made. It is found that the wave resistance predicted by the Green–Naghdi equation is larger than that predicted by the shallow water equations in subcritical flow , and the Green–Naghdi equation and the shallow water equations predict almost the same wave resistance when , the frequency dispersion can be neglected in supercritical flows.     

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.     

A finite element model for wave refraction and diffraction

A two-dimensional hybrid finite element method is developed to study the scattering of water waves by an island and to calculate wave forces and moments on offshore structures. The offshore structure, which could be either semi-submerged or fully extended in the water, is assumed to be stationary. The numerical model is based on the mild-slope equation. It can be applied to both long-wave and short-wave problems. A special treatment for the problem with the semi-submerged structure is introduced. Comparisons are given with existing analytical solutions and other numerical results. The present model is shown to be an efficient and accurate method for the solution of wave refraction and diffraction problems.     

Wave reflection from vertical breakwater with porous structure

This paper presents numerical solutions for the wave reflection from submerged porous structures in front of the impermeable vertical breakwater. A new time-dependent mild-slope equation involves the parameters of the porous medium including the porosity, the friction factor and the inertia coefficient, etc. is derived for solving the boundary value problem. A comprehensive comparison between the present model and the existing analytical solution for the case of simple rectangular geometries of the submerged structure is performed first. Then, more complicated cases such as the inclined and trapezoidal submerged porous structures in front of the vertical breakwater with sloping bottom are considered. This study also examines the effects of the permeable properties and the geometric configurations of the porous structure to the wave reflection. It is found that the submerged porous structure with trapezoidal shape has more efficiency to reduce the wave reflection than that of triangular shape. The numerical results show that the minimum wave reflection is occurred when the breakwater is located at the intermediate water depth.     

An analytic solution to the mild slope equation for waves propagating over an axi-symmetric pit

An analytic solution to the mild slope equation is derived for waves propagating over an axi-symmetric pit located in an otherwise constant depth region. The water depth inside the pit decreases in proportion to an integer power of radial distance from the pit center. The mild slope equation in cylindrical coordinates is transformed into ordinary differential equations by using the method of separation of variables, and the coefficients of the equation in radial direction are transformed into explicit forms by using the direct solution for the wave dispersion equation by Hunt (Hunt, J.N., 1979. Direct solution of wave dispersion equation. J. Waterw., Port, Coast., Ocean Div., Proc. ASCE, 105, 457–459). Finally, the Frobenius series is used to obtain the analytic solution. Due to the feature of the Hunt's solution, the present analytic solution is accurate in shallow and deep waters, while it is less accurate in intermediate depth waters. The validity of the analytic solution is demonstrated by comparison with numerical solutions of the hyperbolic mild slope equations. The analytic solution is also used to examine the effects of the pit geometry and relative depth on wave transformation. Finally, wave attenuation in the region over the pit is discussed.     

New Numerical Scheme for Simulation of Hyperbolic Mild-Slope Equation

The original hyperbolic mild-slope equation can effectively take into account the combined effects of wave shoaling, refraction, diffraction and reflection, but does not consider the nonlinear effect of waves, and the existing numerical schemes for it show some deficiencies. Based on the original hyperbolic mild-slope equation, a nonlinear dispersion relation is introduced in present paper to effectively take the nonlinear effect of waves into account and a new numerical scheme is proposed. The weakly nonlinear dispersion relation and the improved numerical scheme are applied to the simulation of wave transformation over an elliptic shoal. Numerical tests show that the improvement of the numerical scheme makes efficient the solution to the hyperbolic mild-slope equation. A comparison of numerical results with experimental data indicates that the results obtained by use of the new scheme are satisfactory.     

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.     

Numerical implementation of the harmonic modified mild-slope equation

A numerical solver is presented of the modified time-independent mild-slope equation, which incorporates energy dissipation. Using a second-order parabolic approximation, the following external boundary conditions are modelled: open and fully transmitting to both incoming and outgoing waves; partially reflecting, and; fully absorbing. Discretisation of the governing equation and boundary conditions is by means of a second-order accurate central difference scheme. The resulting sparse-banded matrix is solved using an inexpensive banded solver with Gaussian elimination. The numerical predictions are in excellent agreement with the analytical solution for the interaction of non-breaking waves with an array of vertical surface-piercing circular cylinders on a horizontal bed. Results are compared with those for the same array on various seabed topographies. The model is robust and can be used for wave propagation in complex geometries. It has fewer restrictions associated with wave obliqueness at boundaries than traditional models based on the mild-slope equation.     

A theoretical formulation for wave propagations over uneven bottom

A first-order solution is obtained with a mild-slope parameter within the framework of the linear wave theory using the classical multiple scale method. This solution satisfies (1) the governing equations and free surface condition; (2) the coastal bottom condition in terms of multiple scale with its form analytical. As to two-dimensional wave propagating towards a slope, its solution reduces to the famous Biesel [1952. Study of wave progression in water of gradually varying depth. Gravity Waves, vol. 521. US National Bureau of Standards Circular, pp. 243–253] solution.     

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.     

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.     

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.     

Phase-decoupled refraction–diffraction for spectral wave models

Conventional spectral wave models, which are used to determine wave conditions in coastal regions, can account for all relevant processes of generation, dissipation and propagation, except diffraction. To accommodate diffraction in such models, a phase-decoupled refraction–diffraction approximation is suggested. It is expressed in terms of the directional turning rate of the individual wave components in the two-dimensional wave spectrum. The approximation is based on the mild-slope equation for refraction–diffraction, omitting phase information. It does therefore not permit coherent wave fields in the computational domain (harbours with standing-wave patterns are excluded). The third-generation wave model SWAN (Simulating WAves Nearshore) was used for the numerical implementation based on a straightforward finite-difference scheme. Computational results in extreme diffraction-prone cases agree reasonably well with observations, analytical solutions and solutions of conventional refraction–diffraction models. It is shown that the agreement would improve further if singularities in the wave field (e.g., at the tips of breakwaters) could be properly accounted for. The implementation of this phase-decoupled refraction–diffraction approximation in SWAN shows that diffraction of random, short-crested waves, based on the mild-slope equation can be combined with the processes of refraction, shoaling, generation, dissipation and wave–wave interactions in spectral wave models.     

缓坡方程与Boussinesq方程特征的分析比较

在对缓坡方程和Boussinesq方程研究的基础上，从方程的基本形式和特征以及频散关系等方面对二者进行了分析和比较，明确了线性缓坡方程在频散性上要好于非线性Boussinesq方程。此外还对Boussinesq型模型与抛物型缓坡方程模型在Berkhoff椭圆地形的计算结果及其精度也进行比较，计算结果与实测数据吻合很好，说明这两种模型都可以用于模拟近岸波浪传播过程所发生的各种变形。但由于各自控制方程对各物理过程的处理不同，因此各有特征。