Implementation and evaluation of alternative wave breaking formulas in a coastal spectral wave model

This paper describes methods and results of research for incorporating four different parameterized wave breaking and dissipation formulas in a coastal wave prediction model. Two formulations assume the breaking energy dissipation to be limited by the Rayleigh distribution, whereas the other two represent the breaking wave energy by a bore model. These four formulations have been implemented in WABED, a directional spectral wave model based on the wave action balance equation with diffraction, reflection, and wave–current interaction capabilities. Four parameterized wave breaking formulations are evaluated in the present study using two high-quality laboratory data sets. The first data set is from a wave transformation experiment at an idealized inlet entrance, representing four incident irregular waves in a slack tide and two steady-state ebb current conditions. The second data set is from a laboratory study of wave propagation over a complex bathymetry with strong wave-induced currents. Numerical simulation results show that with a proper breaking formulation the wave model can reproduce laboratory data for waves propagating over idealized or complicated bathymetries with ambient currents. The extended Goda wave breaking formulation with a truncated Rayleigh distribution, and the Battjes and Janssen formulation with a bore model produced the best agreement between model and data.     

Prediction of impact pressure induced by breaking waves on vertical cylinders in random seas

This paper presents a method to statistically predict the magnitude of impact pressure (including extreme values) produced by deep water waves breaking on a circular cylinder representing a column of an ocean structure. Breaking waves defined here are not those whose tops are blown off by the wind but those whose breaking is associated with steepness. The probability density function of wave period associated with breaking waves is derived for a specified wave spectrum, and then converted to the probability density function of impact pressure. Impacts caused by two different breaking conditions are considered; one is the impact associated with waves breaking in close proximity to the column, the other is an impact caused by waves approaching the column after they have broken. As an example of the application of the present method, numerical computations are carried out for a wave spectrum obtained from measured data in the North Atlantic.     

Hyperbolic Mild Slope Equations with Inclusion of Amplitude Dispersion Effect： Regular Waves

A new form of hyperbolic mild slope equations is derived with the inclusion of the amphtude dispersion of nonlinear waves. The effects of including the amplitude dispersion effect on the wave propagation are discussed. Wave breaking mechanism is incorporated into the present model to apply the new equations to surf zone. The equations are solved nu- merically for regular wave propagation over a shoal and in surf zone, and a comparison is made against measurements. It is found that the inclusion of the amplitude dispersion can also improve model＇ s performance on prediction of wave heights around breaking point for the wave motions in surf zone.     

Wave Breaking Phenomena of Irregular Waves Combined with Opposing Current

- Experimental study and theoretical analysis show that the critical value of relative wave height (H/ d)b given by Goda and the critical wave steepness (H/ L)b given by Michell and Miche can be adopted as the spilling breaking indices of regular waves. According to the same principle, a systematic theoretical analysis and experiment of irregular waves have been done by the authors in order to solve the breaking problem of irregular waves. It is indicated that the authors' method for determining wave breaking of regular waves can also be used for irregular waves.     

Numerical models for evolution of extreme wave groups

The paper considers the application of two numerical models to simulate the evolution of steep breaking waves. The first one is a Lagrangian wave model based on equations of motion of an inviscid fluid in Lagrangian coordinates. A method for treating spilling breaking is introduced and includes dissipative suppression of the breaker and correction of crest shape to improve the post breaking behaviour. The model is used to create a Lagrangian numerical wave tank, to reproduce experimental results of wave group evolution. The same set of experiments is modelled using a novel VoF numerical wave tank created using OpenFOAM. Lagrangian numerical results are validated against experiments and VoF computations and good agreement is demonstrated. Differences are observed only for a small region around the breaking crest.     

Approach to Calculating Methods of Wave Breaking Depth and Height

- Engineers in coastal engineering have been paying much attention to the research subject on wave breaking. In this paper, previous research results on the calculating methods of wave breaking depth and height are enumerated, the laws of wave transformation before and after wave breaking are investigated, an adequate supposition is made, and the effect of beach slope, bed friction and breaking turbulence on wave breaking is considered. By applying the theory of wave energy dissipation rate and combining with proper formulas of solitary waves, a new calculating formula of wave breaking depth and height on the movable bed is derived and examined with the data from experimental pools of different sizes, and it is proved to be of practicability.     

Effects of wave breaking on wave statistics for deep-water random wave train

The experimental studies of the breaking effects on wave statistics for deep-water random waves are presented. It is especially focused on the behavior of kurtosis of surface elevations due to wave breaking. Wave breaking suppresses the maximum limit of kurtosis of the surface elevation, although skewness depends on characteristic wave steepness. The mean instantaneous wave steepness of breaking waves defined using the zero-down-crossing method was much lower than expected from the Stokes waves.     

波浪作用下方块石护面稳定性研究

护面是海堤和护岸的重要结构,直接抵御波浪作用,可采用人工块体、块石等,种类繁多。采用紧密排列方块石作为护面结构是一种景观性较好的型式,依据方块石厚度不同能抵御不同大小的波浪作用。干砌条石及干砌块石护面曾有一些规范给出过计算方法,但现行规范没有相关内容可供设计参考,已有计算方法的理论分析还存在不足。当波浪与斜坡堤相互作用时,方块石护面出现位移或脱落可能发生在波浪回落最低阶段、波浪破碎打击阶段及破后爬高水流作用阶段,通过研究得到了不同阶段波浪对方块石护面作用力的计算方法。在波浪回落最低阶段,考虑了护面及其下方垫层渗透性影响,通过理论分析建立了低渗透护面浮托压强计算模型,采用物模试验将计算结果与试验测量值进行了对比分析,结果表明总体趋势符合,量值接近;在波浪破碎冲击阶段,基于射流冲击作用原理,提出了波浪在斜坡面破碎冲击压强计算方法,通过试验分析了波浪破碎水深波高比与破波相似参数的关系,利用浅水波理论计算了波浪破碎冲击水流流速;在爬高水流作用阶段,提出了水流引起的方块石护面垂直浮托力及水平拖曳力的计算方法,通过试验结果拟合了浮托力系数和拖曳力系数,验证了水流作用下护面的受力特征。最后,针对方块...     

波群内单个波的波陡分布与波破碎

对波群内单个波的波陡分布和波破碎进行了实验研究。研究结果是,波群中波动的最大振幅出现在波群前部而不是出现在波群中央,这种不对称性导致波群前部单个波出现大波陡的概率大于后部单个波出现大波陡的概率;进一步的波破碎统计发现波群前部单个波破碎的频率是后部单个波破碎频率的4倍。因此认为,波群结构的不对称性能够导致单个波发生破碎的...     

近岸波浪折射-绕射-破波耗散联合模式的有限元数值研究

建立了近岸波浪折射-绕射-破波耗散的有限元数值模式。采用的有限元方法为改进的混合元法,其中外域开边界条件得到改进,内域有限元为伽廖金有限元。用理论解检验了所建立的数值模式,并将该模式应用到一个模型港湾。     

基于激波捕捉类Boussinesq模型对岛礁地形上波浪传播的数值模拟

为了探究激波捕捉类Boussinesq模型在模拟岛礁地形上规则波和不规则波传播的可行性,采用基于完全非线性Boussinesq方程并具有激波捕捉能力的数值模型Funwave-TVD对规则波和不规则波在岛礁地形上的传播进行了数值模拟,通过与试验数据对比,分析模型中空间步长的影响,验证模型在模拟波高、平均水位分布以及波谱空间演变的能力,结果表明:采用合适的空间步长,模型能较好地模拟规则波和不规则波在岛礁地形上的传播和演化过程。对于规则波,较小的空间步长可改善破碎点处波高峰值的预测,并能更好地预测波浪破碎后波高的空间分布,相比结合经验破碎的Boussinesq模型,Funwave-TVD能更好地模拟规则波在岛礁地形上的破碎,以及破碎以后行进涌波的再生成过程;对于不规则波,Funwave-TVD总体而言能较好地模拟涌浪有效波高、次重力波的生成及空间演化和平均水位,但会低估礁坪上次重力波波高,较粗的空间步长也会低估礁坪上涌浪有效波高。     

Large-wave simulation of three-dimensional,cross-shore and oblique,spilling breaking on constant slope beach

In this work, the large-wave simulation (LWS) method is adapted for application in spilling wave breaking over a constant slope beach. According to LWS, large scales of velocities, pressure and free-surface elevation are numerically resolved, while the corresponding unresolved scale effects are taken into consideration by a subgrid scale (SGS) model for wave and eddy stresses. The model may be not fully applicable in very shallow water, close to the shoreline, where the unresolved, turbulent, free-surface oscillation is of the same order with the water depth. Time integration of the Euler equations is achieved by a two-stage fractional scheme, combined with a hybrid scheme for spatial discretization, consisting of finite difference and pseudospectral approximation methods. Model parameters are calibrated by comparison to available experimental data of free-surface elevation and velocities in the surf zone for cross-shore incoming waves. The action of the wave SGS stresses in the outer coastal and surf zones initiates breaking and generates appropriate vorticity, in the form of an eddy structure (surface roller), at the breaking wavefront. At incipient breaking, both advection and gravity contribute to the vorticity flux at the free surface, while only after the full development of the surface roller, the effect of advection becomes stronger. The SGS model is also utilized to simulate propagation, refraction and breaking of oblique incoming waves. The gradual breaking and dissipation of wave crestlines and the surface roller structure along the breaking wavefront are automatically captured without any empirical input, such as data for the roller shape or the wave propagation angle at breaking.     

海浪特性及其传播变形研究新进展

综述第27届国家海岸工程会议（ICCE2000）上有关海浪特性及其传播变形研究的新进展，内容包括海浪的时域特性，海浪的特征，海浪的传播变形和波浪的破碎，反射等四部分。     

Breaking of Solitary Waves on Uniform Slopes

The propagation,shoaling and breaking of solitary waves on mild slopes are simulated byboundary element method.In this paper,the criterion of breaking solitary waves on mild slopes is discussed.The criterion is that the ratio of horizontal velocity of water particles on the wave crest to wave celerity equalsone.However,the case that the ratio of horizontal velocity of water particles on the wave crest to wave ce-lerity is below one but the front face of wave profile becomes vertical is also considered as a breaking criteri-on.According to the above criteria,the breaking index for slopes 1:10 to 1:25 is studied.The result is com-pared to other researchers'.The deformation of solitary waves on slopes is discussed and the distribution offluid velocities at breaking is shown.     

Experimental investigation of breaking criteria of deepwater wind waves under strong wind action

The present study describes an experimental investigation of breaking criteria of deepwater wind waves under strong wind action. In a wind wave flume, waves were generated using different wind speeds and measured at different locations to obtain wave trains of no, intermittent, or frequent breaking. Water particle movement and free surface elevation were measured simultaneously using a PIV system and a wave gauge, respectively. For wind waves, not all the waves measured at a fixed location are breaking waves, and the breaking of a larger wave is not guaranteed. However, the larger the wave height, the larger the probability of breaking. In order to take as many breaking waves as possible for the cases of frequent breaking, we used the waves whose heights were close to the highest one-tenth wave height. The experimental results showed that the geometric or kinematic breaking criteria could not explain the occurrence of breaking of wind waves. On the other hand, the vertical acceleration beneath the wave crest was close to the previously suggested limit value, −0.5g, when frequent breaking of large waves occurred, indicating that the dynamic breaking criterion would be good for discriminating breaking waves under a strong wind action.     

The wave transformation and breaking phenomena in shallow water

Ｔｈｅｗａｖｅｔｒａｎｓｆｏｒｍａｔｉｏｎａｎｄｂｒｅａｋｉｎｇｐｈｅｎｏｍｅｎａｉｎｓｈａｌｌｏｗｗａｔｅｒ￥ＬｉＹｕｃｈｅｎｇ（１．ＤａｌｉａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ，Ｄａｌｉａｎ１１６０２３，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：...     

Laboratory study of wave and turbulence characteristics in narrow-band irregular breaking waves

Wave elevations and water particle velocities were measured in a laboratory surf zone created by the breaking of a narrow-band irregular wave train on a 1/35 plane slope. The incident waves form wave groups that are strongly modulated. It is found that the waves that break close to the shoreline generally have larger wave-height-to-water-depth ratios before breaking than the waves that break farther offshore. After breaking, the wave-height-to-water-depth ratio for the individual waves approaches a constant value in the inner surf zone, while the standard deviation of the wave period increases as the still water depth decreases. In the outer surf zone, the distribution of the period-averaged turbulent kinetic energy is closely correlated to the initial wave heights, and has a wider variation for narrow-band waves than for broad-band waves. In the inner surf zone, the distribution of the period-averaged turbulent kinetic energy is similar for narrow-band waves and broad-band waves. It is found that the wave elevation and turbulent kinetic energy time histories for the individual waves in a wave group are qualitatively similar to those found in a spilling regular wave. The time-averaged transport of turbulent kinetic energy by the ensemble-averaged velocity and turbulence velocity under the irregular breaking waves are also consistent with the measurements obtained in regular breaking waves. The experimental results indicate that the shape of the incident wave spectrum has a significant effect on the temporal and spatial variability of wave breaking and the distribution of turbulent kinetic energy in the outer surf zone. In the inner surf zone, however, the distribution of turbulent kinetic energy is relatively insensitive to the shape of the incident wave spectrum, and the important parameters are the significant wave height and period of the incident waves, and the beach slope.     

Probabilistic Models for the Probability of Wave Breaking and Whitecap Coverage Based on Kinematic Breaking Criterion

More and more researches show that neither the critical downward acceleration nor the critical slope of water waves is a universal constant. On the contrary, they vary with particular wave conditions. This fact moders the models either for the probability of wave breaking B or for the whitecap coverage W based on these criteria difficult to apply. In this paper and the one which follows we seek to develop models for the prediction of both B and W based on the kinematical criterion. First, several joint probabihstic distribution functions (PDFs) of wave characteristics are derived, based on which the breaking properties B and W are estimated. The estimation is made on the assumption that a wave breaks ff the horizontal velocity of water particles at its crest exceeds the local wave celerity, and whitecapping occurs in regions of fluid where water particles travel faster than the waves. The consequent B and W depend on wave spectral moments of orders 0 to 4.Then the JONSWAP spectrum is used to represent the fetch-limited sea waves in deep water, so as to relate the probahility of wave breaking and the whitecap coverage with wind parameters. To this end, the time-averaging technique proposed by Glazman (1986) is applied to the estimation of the spectral moments involved, and furthermore, the theoretical models are compared with available observations collected from published literature. From the comparison, the averaging time scale is determined. The final models show that the probability of wave breaking as well as the whitecap coverage depends on the dimensionless fetch. The agreement between these models and the database is reasonable.     

Experimental study on suspended sediment concentration and its vertical distribution under spilling breaking wave actions in silty coast

In this paper,flume experiments are focused on sediment transport inside and outside the surf zone.According to the energy dissipation balance principle of sediment-laden flow and the similarity between energy dissipation of spilling breaking wave and hydraulic jump,formulas are proposed to predict time averaged suspended sediment concentration under both non-breaking and breaking waves.Assuming that the sediment diffusion coefficient,which is related with energy dissipation,is proportional to water depth,formulas are proposed to predict close-to-bed suspended sediment concentration and vertical distribution of suspended sediment under spilling breaking waves,and the prediction shows a good agreement with the measurement.     

Modeling waves and longshore transport potential in Half Moon Bay, Grays Harbor, Washington, USA

A local-scale phase-resolving wave transformation model with CGWAVE is established in connection with a regional-scale coupled STWAVE-ADCIRC wave-current model for its application in the Half Moon Bay, Grays Harbor. Wave transformation from offshore to the harbor entrance is simulated by the STWAVE model which includes wave-current interaction. The STWAVE results provide incident wave conditions for the local-scale CGWAVE model at its outer boundary. A simple method is developed to take into ac- count the lateral variation of wave height in constructing the model's wave boundary conditions. The model was validated for three wave condition cases which yielded good agreement with field data. The validated model was applied to predicting nearshore waves in the Half Moon Bay and longshore transport parameters along the wave breaking line for the existing condition and three engi- neering alternatives. A comparative analysis indicated that storm waves that have a combination of long period and large height are the most destructive to the crenulate shoreline in the Half Moon Bay; both 152 m jetty extension (Alt. 2) and diffraction mound enlargement ( Alt. 3) would significantly reduce breaking wave height and longshore transport potential in the southwest comer of Half Moon Bay.