A Finite-Difference Approach to the Time-Dependent Mild-Slope Equation

A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity.The Euler predictor-corrector method and the three-point finite-difference method with varying spatial steps are adopted to discretize the time derivatives and the two-dimensional horizontal ones,respectively,thus leading both the time and spatial derivatives to the second-order accuracy.The boundary conditions for the present model are treated on the basis of the general conditions for open and fixed boundaries with an arbitrary reflection coefficient and phase shift.Both the linear and nonlinear versions of the numerical model are applied to the wave propagation and transformation over an elliptic shoal on a sloping beach,respectively,and the linear version is applied to the simulation of wave propagation in a fully open rectangular harbor.From comparison of numerical results with theoretical or experimental ones,it is found that they are in reasonable agreement.     

一般曲线坐标系下波浪传播的数值模拟

在曲线坐标系下,建立了缓变水深水域波浪传播的数值模拟模型.模型适宜于复杂变化的边界形状,克服了各种代数坐标变换的局限性.在建立模型时,将原始的椭圆型缓坡方程的近似型式——依赖时间变化的抛物型方程,作为控制方程,既克服了一般抛物近似方法的缺点,又便利了方程的求解;从开边界条件、不同反射特性的固壁边界条件相统一的表达式出发,对边界条件进行处理;用ADI法数值求解控制方程.对模型的验证表明,数值解与物模实验值吻合良好,模型对于具有复杂边界的工程实际有较强的适应性.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

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.     

Wave Numerical Model for Shallow Water

The history of forecasting wind waves by wave energy conservation equation is briefly des-cribed.Several currently used wave numerical models for shallow water based on different wave theoriesare discussed.Wave energy conservation models for the simulation of shallow water waves are introduced,with emphasis placed on the SWAN model,which takes use of the most advanced wave research achieve-ments and has been applied to several theoretical and field conditions.The characteristics and applicabilityof the model,the finite difference numerical scheme of the action balance equation and its source termscomputing methods are described in detail.The model has been verified with the propagation refractionnumerical experiments for waves propagating in following and opposing currents;finally.the model is ap-plied to the Haian Gulf area to simulate the wave height and wave period field there,and the results arecompared with observed data.     

Practical radiation stress calculations connected with equations of wave propagation

Numerical models of combined surface gravity wave refraction, diffraction and reflection can be solved conveniently in terms of the water surface displacement, η, and a vertically integrated, wave-induced, water particle velocity, Q. However, the normal formulation for the radiation stress components, expressed in terms of the wave energy, wave number and water depth, is correct only for linear progressive waves.This paper describes a method of calculating the radiation stress for a linear progressive wave plus an arbitrary reflected or back-scattered wave in terms of variables η and Q. The calculations are related to a finite-difference scheme. Correction factors are given which compensate for the errors introduced by the use of finite differentials in the calculation of certain elements of radiation stress.The theory upon which the analytical method is based is only exact for water of uniform depth. However, results are presented which show that the errors are not significant (typical error less than ± 2%) for bed slopes of less than 1:3.     

任意曲线边界条件下缓变水深水域波浪传播的数值模拟

缓坡方程被广泛地应用于描述波浪的传播变形计算,目前一般采用矩形网格求解.将计算域剖分为任意四边形网格,以格林公式为基础,在变量沿单元边界线性变化的假定下,对双曲型的波能守恒方程、波数矢无旋性方程进行离散,同时通过等参单元变换推求节点偏导数值以离散椭圆型光程函数方程,从而建立了任意曲线边界条件下缓变水深水域波浪传播的数值模拟模型.将模型应用于平行直线型等深线地形,并将计算域剖分为不规则四边形网格,对不同入射角、底坡、波高等多种组合情况比较了数值解与解析解,结果表明两者一致.应用于复杂边界的实例,数值模拟结果与物模实验值基本吻合.     

Numerical Simulation of Non-Linear Wave Propagation in Waters of Mildly Varying Topography with Complicated Boundary

In this paper, the characteristics of different forms of mild slope equations for non-linear wave are analyzed, and new non-linear theoretic models for wave propagation are presented, with non-linear terms added to the mild slope equations for non-stationary linear waves and dissipative effects considered. Numerical simulation models are developed of non-linear wave propagation for waters of mildly varying topography with complicated boundary, and the effects are studied of different non-linear corrections on calculation results of extended mild slope equations. Systematical numerical simulation tests show that the present models can effectively reflect non-linear effects.     

Generalized boussinesq model for periodic non-linear shallow-water waves

This paper presents the development of a generalized Boussinesq (gB) model for the periodic non-linear shallow-water waves. An incident cnoidal wave solution for the gB model is derived and applied to the wave simulation. A set of radiation boundary conditions is also established to transmit effectively the cnoidal waves out of the computational domain. The classical solutions of the second-order cnoidal waves are discussed within the content of the KdV equation and the generalized Boussinesq equations. An Euler's predictor-corrector finite-difference algorithm is used for numerical computation. The propagation of normally incident cnoidal waves in a channel is studied. The simulated wave profiles agree well with the analytical results. The temporal and spatial evolution of an obliquely incident cnoidal wave is also modelled. The phenomenon of Mach reflection is discussed.     

Experimental and Numerical Study on the Hydrodynamic Characteristics of Solitary Waves Passing Over A Submerged Breakwater

In this study, solitary waves passing over a submerged breakwater are investigated both experimentally and numerically. A total of 9 experimental conditions are carried out, including different incident wave heights and water depths. Numerical simulations are performed using a high-order finite-difference model solving Navier–Stokes (N–S) equations. The predicted water wave elevation, velocity and pressure show good agreement with experimental data, verifying the accuracy and capacity of the numerical model. Furthermore, parametric studies are conducted by numerical modelling to examine the effects of the geometrical features of submerged dike on hydrodynamic characteristics around the breakwater.     

Numerical Modeling of Wave Diffraction-Refraction in Water of Varying Current and Topography

—A numerical model for wave diffraction-refraction in water of varying current and topogra-phy is proposed,and time-dependent wave mild-slope equation with a dissipation term and correspondingequivalent governing equations are presented.Two different expressions of parabolic approximations forthe case of the absence of current are also given and analyzed.The influence of current on the results ofsimulation of waves is discussed.Some examples show that the present model is better than others in simu-lating wave transformation in large water areas.And they also show that the influence of current shouldbe taken into account,on numerical modeling of wave propagation in water of strong current and coastalareas,otherwise the modeling results will be largely distorted.     

Wave transformation by axisymmetric three-dimensional bathymetric anomalies with gradual transitions in depth

The development of an analytic model (Axisymmetric 3-D Step Model) for the propagation of linear water waves over an axisymmetric bathymetric anomaly in arbitrary water depth is presented. The Axisymmetric 3-D Step Model is valid in a region of uniform depth containing an axisymmetric bathymetric anomaly with gradual transitions in depth allowed as a series of steps approximating arbitrary slopes. The velocity potential is calculated by applying matching conditions at the interface between regions of constant depth. The velocity potential obtained determines the wave field in the domain for monochromatic incident waves of linear form. A second analytic model (3-D Shallow Water Exact Model) is developed for comparison within the shallow water limit.The Axisymmetric 3-D Step Model determines the wave transformation caused by the processes of wave refraction, diffraction and reflection. Wave transformation is demonstrated in plots of the relative amplitude for bathymetric anomalies in the form of pit or a shoal, highlighting areas of wave sheltering and wave focusing. Anomalies of constant volume, but variable cross-section are employed to isolate the effect of the transition slope on the wave transformation.Comparisons to a shallow water model, numerical models, and experimental data verify the results of the Axisymmetric 3-D Step Model for several bathymetries including both pits and shoals. Also included are estimates of the energy reflection induced by an axisymmetric depth anomaly. The 3-D Axisymmetric Step Model has been applied previously to account for nearshore transformation (sloping bathymetry) and associated shoreline changes [C.J. Bender, R.G. Dean, Coastal Engineering 51 (2004) 1143].     

An explicit finite difference model for simulating weakly nonlinear and weakly dispersive waves over slowly varying water depth

In this paper, a modified leap-frog finite difference (FD) scheme is developed to solve Non linear Shallow Water Equations (NSWE). By adjusting the FD mesh system and modifying the leap-frog algorithm, numerical dispersion is manipulated to mimic physical frequency dispersion for water wave propagation. The resulting numerical scheme is suitable for weakly nonlinear and weakly dispersive waves propagating over a slowly varying water depth. Numerical studies demonstrate that the results of the new numerical scheme agree well with those obtained by directly solving Boussinesq-type models for both long distance propagation, shoaling and re-fraction over a slowly varying bathymetry. Most importantly, the new algorithm is much more computationally efficient than existing Boussinesq-type models, making it an excellent alternative tool for simulating tsunami waves when the frequency dispersion needs to be considered.     

Boundary element modeling of multidirectional random waves in a harbor with partially reflecting boundaries

A numerical model is presented for the prediction of the wave field due to the diffraction of directional random waves in a harbor of arbitrary shape with partially reflecting boundaries. The water depth is assumed uniform and the method is based upon the superposition of diffraction solutions for monochromatic waves obtained by a two-dimensional boundary integral equation approach. The incident wave conditions are specified using a discrete form of the Mitsuyasu directional spectrum. The present numerical model has been validated through comparisons with previous experimental data and theoretical results for both regular and random wave diffraction by offshore breakwaters and in harbors. Good agreement was obtained in all cases. Based on these comparisons it is concluded that the present numerical model is an accurate and efficient tool to predict the wave field inside a harbor or around a breakwater in many practical applications.     

Numerical modeling of the propagation of nonlinear acoustic-gravity waves in the middle and upper atmosphere

A numerical algorithm for modeling the vertical propagation and breaking of nonlinear acoustic-gravity waves (AGWs) from the Earth’s surface to the upper atmosphere is described in brief. Monochromatic variations in the vertical velocity at the Earth’s surface are used as an AGW source in the model. The algorithm for solving atmospheric hydrodynamic equations is based on three-dimensional finite-difference analogues of fundamental conservation laws. This approach selects physically correct generalized solutions to hydrodynamic equations. A numerical simulation is carried out in an altitude region from the Earth’s surface to 500 km. Vertical profiles of the background temperature, density, and coefficients of molecular viscosity and heat conduction are taken from the standard atmosphere models. Calculations are made for different amplitudes of lower-boundary wave forcing. The AGW amplitudes increase with altitude, and waves may break in the middle and upper atmosphere.     

A comparison of several wave spectra for the random wave diffraction by a semi-infinite breakwater

A comparison of the diffraction of multidirectional random waves using several selected wave spectrum models is presented in this paper. Six wave spectrum models, Bretschneider, Pierson–Moskowitz, ISSC, ITTC, Mitsuyasu, and JONSWAP spectrum, are considered. A discrete form for each of the given spectrum models is used to specify the incident wave conditions. Analytical solutions based on both the Fresnel integrals and polynomial approximations of the Fresnel integrals and numerical solutions of a boundary integral approach have been used to obtain the two-dimensional wave diffraction by a semi-infinite breakwater at uniform water depth. The diffraction of random waves is based on the cumulative superposition of linear diffraction solution. The results of predicted random wave diffraction for each of the given spectrum models are compared with those of the published physical model presented by Briggs et al. [1995. Wave diffraction around breakwater. Journal of Waterway, Port, Coastal and Ocean Engineering—ASCE 121(1), 23–35]. Reasonable agreement is obtained in all cases. The effect of the directional spreading function is also examined from the results of the random wave diffraction. Based on these comparisons, the present model for the analysis of various wave spectra is found to be an accurate and efficient tool for predicting the random wave field around a semi-infinite breakwater or inside a harbor of arbitrary geometry in practical applications.     

Correlation of extreme waves and water levels using a third-generation wave model and a 3D flow model

A procedure to correlate extreme wave heights and extreme water levels in coastal waters using numerical models together with joint probability analysis has been proposed. A third-generation wave model for wave simulation and a three-dimensional flow model for water level simulation are coupled through the surface atmospheric boundary layer. The model has been calibrated and validated against wind, wave and water level data collected in the coastal waters of Hong Kong. The annual maximum wave height and the concomitant water level have been obtained by simulating the annual extreme typhoon event for 50 consecutive years. The results from bivariate extreme value analysis of the simulated data show that the commonly used empirical method may lead to underestimation of the design water level.     

非线性波浪波面追踪的一种新模式

基于Laplace方程的Green积分表达式和波面BemouUi方程所建立的非线性波动数学模型,是一个时域上具有初始值的边值问题,而精确地追踪自由表面的波动位置,给出波面运动瞬时的波面高度和波面势函数,是建立时域内非线性波浪数值模式的基础。本文采用0-1混合型边界元剖分计算域边界并离散Laplace方程的Green积分表达式,采用有限元剖分自由水面并推导满足自由表面非线性边界条件的波面有限元方程,联立计算域内以节点波势函数和波面位置高度的时间增量为未知量的线性方程组,通过时步内的循环迭代,给出每个时步上的波面位置和波面势函数,从而建立了一种新的非线性波浪波面追踪模式。数值造波水槽内的波浪试验表明,其数值模拟结果具有良好的计算精度。     

Internal generation of waves on an arc in a rectangular grid system

This paper presents a technique to generate waves at oblique angles in finite difference numerical models in a rectangular grid system by using internal generation technique [Lee, C., Suh, K.D., 1998. Internal generation of waves for time-dependent mild-slope equations. Coast. Eng. 34, 35–57.] along an arc-shaped line source. Tests were made for four different types of wave generation layouts. Quantitative experiments were conducted under the following conditions: the propagation of waves on a flat bottom, the refraction and shoaling of waves on a planar slope, and the diffraction of waves to a semi-infinite breakwater. Numerical experiments were conducted using the extended mild-slope equations of Suh et al. [Suh, K.D., Lee, C., Park, W.S., 1997. Time-dependent equations for wave propagation on rapidly varying topography. Coast. Eng. 32, 91–117.]. The fourth layout type consisting of two parallel lines connected to a semicircle showed the best solutions, especially for a small grid size. This technique is useful for the numerical simulation of irregular waves with broad-banded directional spectrum using conventional spectral wave models for the reasonable estimation of bottom friction and wave-breaking.