The BCI criterion for the initiation of breaking process in Boussinesq-type equations wave models

The Breaking Celerity Index (BCI) is proposed as a new wave breaking criterion for Boussinesq-type equations wave propagation models (BTE).The BCI effectiveness in determining the breaking initiation location has been verified against data from different experimental investigations conducted with incident regular and irregular waves propagating along uniform slope [Utku, M. (1999). “The Relative Trough Froude Number. A New Criteria for Wave Breaking”. Ph.D. Dissertation, Dept. of Civil and Enviromental Engineering, Old Dominion University, Norfolk, VA; Gonsalves Veloso dos Reis, M.T.L. (1992). “Characteristics of waves in the surf zone”. MS Thesis, Department of Civil Engineering, University of Liverpool., Liverpool; Lara, J.L., Losada, I.J., and Liu, P.L.-F. (2006). “Breaking waves over a mild gravel slope: experimental and numerical analysis”. Journal of Geophysical Research, VOL 111, C11019] and barred beaches [Tomasicchio, G.R., and Sancho, F. (2002). “On wave induced undertow at a barred beach”. Proceedings of 28th International Conference on Coastal Engineering, ASCE, New York, 557–569]. The considered experiments were carried out in small-scale and large-scale facilities. In addition, one set of data has been obtained by the use of the COBRAS model based upon the Reynolds Averaged Navier Stokes (RANS) equations [Liu, P.L.-F., Lin, P., Hsu, T., Chang, K., Losada, I.J., Vidal, C., and Sakakiyama, T. (2000). “A Reynolds averaged Navier–Stokes equation model for nonlinear water wave and structure interactions”. Proceedings of Coastal Structures ‘99, Balkema, Rotterdam, 169–174; Losada, I.J., Lara, J.L., and Liu, P.L.-F. (2005). “Numerical simulation based on a RANS model of wave groups on an impermeable slope”. Proceedings of Fifth International Symposium WAVES 2005, Madrid].Numerical simulations have been performed with the 1D-FUNWAVE model [Kirby, J.T., Wei, G., Chen, Q., Kennedy, A.B., and Dalrymple, R.A. (1998). “FUNWAVE 1.0 Fully Nonlinear Boussinesq Wave Model Documentation and User's Manual”. Research Report No CACR-98-06, Center for Applied Coastal Research, University of Delaware, Newark]. With regard to the adopted experimental conditions, the breaking location has been calculated for different trigger mechanisms [Zelt, J.A. (1991). “The run-up of nonbreaking and breaking solitary waves”. Coastal Engineering, 15, 205–246; Kennedy, A.B., Chen, Q., Kirby, J.T., and Dalrymple, R.A. (2000). “Boussinesq modeling of wave transformation, breaking and run-up. I: 1D”. Journal of Waterway, Port, Coastal and Ocean Engineering, 126, 39–47; Utku, M., and Basco, D.R. (2002). “A new criteria for wave breaking based on the Relative Trough Froude Number”. Proceedings of 28th International Conference on Coastal Engineering, ASCE, New York, 258–268] including the proposed BCI.The calculations have shown that BCI gives a better agreement with the physical data with respect to the other trigger criteria, both for spilling and plunging breaking events, with a not negligible reduction of the calculation time.     

The effect of relative crest submergence on wave breaking over submerged slopes

Experiments were performed in a wave flume to measure the intensity, transmission and reflection of waves breaking over a submerged reef with an offshore gradient of 1:10. The results demonstrate that the relative water depth over the reef crest (h_c/H_o) is a dominant factor affecting the breaking characteristics. In particular it is found that as the relative crest submergence is reduced, there is a considerable increase in the intensity of wave breaking over the reef that can be quantified through measurements of the air cavity enclosed beneath the plunging jet. It is also shown that there is a corresponding decrease in wave transmission and reflection as the submergence is reduced.     

A note on wave energy dissipation over steep beaches

This paper revisits the derivation of the parametric surf zone model proposed by Baldock et al. [Baldock, T. E., Holmes, P., Bunker, S. & Van Weert, P. 1998 Cross-shore hydrodynamics within an unsaturated surf zone. Coast. Eng. 34, 173–196.]. We show that a consistent use of the proposed Rayleigh distribution for surf zone wave heights results in modification of the expressions for the bulk dissipation rate and enhanced dissipation levels on steep beaches and over-saturated surf zone conditions. As a consequence, the modification proposed herein renders the model robust even on steep beaches where it could otherwise develop a shoreline singularity.     

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.     

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.     

Simulation of wave breaking effects in two-dimensional elliptic harbor wave models

A technique is developed for including the effects of dissipation due to wave breaking in two-dimensional elliptic models based on the mild-slope wave equation. This involves exploration of convergence properties pertaining to iteration due to presence of the nonlinear wave breaking parameter in the governing equations as well as new boundary conditions that include wave-breaking effects. Five wave-breaking formulations are examined in conjunction with the resulting model, which is applied to tests involving a sloping beach, a bar-trough bottom configuration, shore-connected and shore-parallel breakwaters on a sloping beach, and two real-world cases. Model results show that three of the formulations, when used within the context of the modeling scheme presented here, provide excellent results compared to data.     

The wave transformation and breaking phenomena in shallow water

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

A new predictive formula for inception of regular wave breaking

Efforts are made to enhance the predictive formula for the inception of wave breaking. To achieve success, the existing formulas are extensively reviewed. They are categorized into four types, i.e., the McCowan type, the Miche type, the Goda type and the Munk type. The inherent relations among the different types are then exploited. The differences among each formula within a group are also discussed. Four representative formulas from the different types are chosen to compare with the measured data for a total number of 1193 cases reported in literatures. It is shown that Goda's and Ostendorf and Madsen's formulas are advantageous in general among the selected ones. Goda's formula, however, is found to be inaccurate as the beach slope becomes steeper than 1/10. Ostendorf and Madsen's formula is fairly good even for cases of very steep slopes, but its accuracy for the cases of ordinary slopes is not as good as Goda's. A new predictive formula for the inception of wave breaking is proposed. The unique index, defined by ψ′_b = (1.21 − 3.30λ_b)(1.48 − 0.54γ_b)ψ_b, where ψ_b = gH_b/C_b², H_b is the breaking wave height, C_b is the breaking wave celerity, λ_b is the breaking wave steepness, γ_b is the relative breaking wave height, and g is the gravity acceleration, is introduced. The incipient condition of wave breaking is then given by ψ′_b = 0.69. This formula is a significant improvement to the existing ones in terms of the accuracy. In addition, it is a local relation. Further verification shows that the proposed formula performs similarly well when applied to the field and to the waves over permeable bed.     

Characteristics of breaking irregular wave forces on a monopile

The substructures of offshore wind turbines are subjected to extreme breaking irregular wave forces. The present study is focused on investigating breaking irregular wave forces on a monopile using a computational fluid dynamics (CFD) based numerical model. The breaking irregular wave forces on a monopile mounted on a slope are investigated with a numerical wave tank. The experimental and numerical irregular free surface elevations are compared in the frequency-domain for the different locations in the vicinity of the cylinder. A numerical analysis is performed for different wave steepness cases to understand the influence of wave steepness on the breaking irregular wave loads. The wave height transformation and energy level evolution during the wave shoaling and wave breaking processes is investigated. The higher-frequency components generated during the wave breaking process are observed to play a significant role in initiating the secondary force peaks. The free surface elevation skewness and spectral bandwidth during the wave transformation process are analysed and an investigation is performed to establish a correlation of these parameters with the breaking irregular wave forces. The role of the horizontal wave-induced water particle velocity at the free surface and free surface pressure in determining the breaking wave loads is highlighted. The higher-frequency components in the velocity and pressure spectrum are observed to be significant in influencing the secondary peaks in the breaking wave force spectrum.     

Improved representation of breaking wave energy dissipation in parametric wave transformation models

An improved formulation to describe breaking wave energy dissipation is presented and incorporated into a previous parametric cross-shore wave transformation model [Baldock, T.E., Holmes, P., Bunker, S., Van Weert, P., 1998. Cross-shore hydrodynamics within an unsaturated surf zone. Coastal Engineering 34, 173–196]. The new formulation accounts for a term in the bore dissipation equation neglected in some previous modelling, but which is shown to be important in the inner surf zone. The only free model parameter remains the choice of γ, the ratio of wave height to water depth at initial breaking, and a well-established standard parameter is used for all model runs. The proposed model is compared to three sets of experimental data and a previous version of the model which was extensively calibrated against field and laboratory data. The model is also compared to the widely used model presented by Thornton and Guza (1983) [Thornton, E.B., Guza, R.T., 1983. Transformation of wave height distribution. Journal of Geophysical Research 88 (No.C10), 5925–5938].     

可渗透海床单桩孤立波破碎波浪力数值研究

当波浪传播至近岸浅水区时易发生破碎,波浪破碎后水体直接拍击单桩结构,其波浪作用力显著增大,可能导致结构失稳破坏。首先建立包括斜坡海床的数值波浪水槽,并与已有研究进行对比验证。进而开展考虑斜坡海床可渗透性的孤立波数值模拟,分析孤立波传播与浅水化破碎特征,着重研究竖直单桩上破碎波浪力的特性,及其与单桩位置、海床渗透率的关联性。数值研究发现:对于低渗透率海床,随着单桩位置由深水向岸线位置变动,其所受波浪力先增后减,所受波浪力最大的桩体位置随海床渗透率增加而从波浪破碎点前方移动至破碎点后1D处,而在高渗透率海床上不同位置处桩体所受波浪力均较小;随海床渗透率等梯度增加,海床消波作用逐渐增强,波浪破碎进程延缓,波浪破碎点向岸线方向加速移动,单桩上破碎波浪力呈整体下降趋势,但可能因波浪破碎点的位置变动导致部分位置桩体所受波浪力异常增大。     

Investigating the impact of surface wave breaking on modeling the trajectories of drifters in the northern Adriatic Sea during a wind-storm event

An accurate numerical prediction of the oceanic upper layer velocity is a demanding requirement for many applications at sea and is a function of several near-surface processes that need to be incorporated in a numerical model. Among them, we assess the effects of vertical resolution, different vertical mixing parameterization (the so-called Generic Length Scale –GLS– set of k–ε, k–ω, gen, and the Mellor–Yamada), and surface roughness values on turbulent kinetic energy (k) injection from breaking waves.First, we modified the GLS turbulence closure formulation in the Regional Ocean Modeling System (ROMS) to incorporate the surface flux of turbulent kinetic energy due to wave breaking. Then, we applied the model to idealized test cases, exploring the sensitivity to the above mentioned factors. Last, the model was applied to a realistic situation in the Adriatic Sea driven by numerical meteorological forcings and river discharges. In this case, numerical drifters were released during an intense episode of Bora winds that occurred in mid-February 2003, and their trajectories compared to the displacement of satellite-tracked drifters deployed during the ADRIA02-03 sea-truth campaign.Results indicted that the inclusion of the wave breaking process helps improve the accuracy of the numerical simulations, subject to an increase in the typical value of the surface roughness z₀. Specifically, the best performance was obtained using α_CH = 56,000 in the Charnok formula, the wave breaking parameterization activated, k–ε as the turbulence closure model. With these options, the relative error with respect to the average distance of the drifter was about 25% (5.5 km/day). The most sensitive factors in the model were found to be the value of α_CH enhanced with respect to a standard value, followed by the adoption of wave breaking parameterization and the particular turbulence closure model selected.     

A coupled numerical model for simulation of wave breaking and hydraulic performances of a composite seawall

A numerical model is developed by combining a porous flow model and a two-phase flow model to simulate wave transformation in porous structure and hydraulic performances of a composite type low-crest seawall. The structure consists of a wide submerged reef, a porous terrace at the top and an impermeable rear wall. The porous flow model is based on the extended Navier-Stokes equations for wave motion in porous media and k−ε turbulence equations. The two-phase flow model combines the water domain with the air zone of finite thickness above water surface. A unique solution domain is established by satisfying kinematic boundary condition at the interface of air and water. The free surface advection of water wave is modeled by the volume of fluid method with newly developed fluid advection algorithm. Comparison of computed and measured wave properties shows reasonably good agreement. The influence of terrace width and structure porosity is investigated based on numerical results. It is concluded that there exist optimum value of terrace width and porosity that can maximize hydraulic performances. The velocity distributions inside and in front of the structure are also investigated.     

Dynamic analysis of a vertical plate exposed to breaking wave impact

From the experimental studies in recent years, it has become known that when a wave breaks directly on a vertical faced coastal structure, high magnitude impact pressures are produced. The theoretical and experimental studies show that the dynamic response of such structures under wave impact loading is closely dependent on the magnitude and duration of the load history. The dynamic analysis and design of a coastal structure can be succeeded provided the design load history for the wave impact is available. Since these types of data are very scarce, it is much more convenient to follow a method which is based on static analysis for the dynamic design procedure. Therefore, to facilitate the dynamic design of a vertical plate that is exposed to breaking wave impact, a multiplication factor called “dynamic magnification factor” is herein presented which is defined as the ratio of the maximum value of the dynamic response to that found by static analysis. The computational results of the present study show that the dynamic magnification factor is a useful ratio to transfer the results of static analysis to the dynamic design of a coastal plate for the maximum impact pressure conditions of p_max/γH₀≤18.     

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.     

Verification of a VOF-based two-phase flow model for wave breaking and wave–structure interactions

The objective of the present study is to develop a volume of fluid (VOF)-based two-phase flow model and to discuss the applicability of the model to the simulation of wave–structure interactions. First, an overview of the development of VOF-type models for applications in the field of coastal engineering is presented. The numerical VOF-based two-phase flow model has been developed and applied to the simulations of wave interactions with a submerged breakwater as well as of wave breaking on a slope. Numerical results are then compared with laboratory experimental data in order to verify the applicability of the numerical model to the simulations of complex interactions of waves and permeable coastal structures, including the effects of wave breaking. It is concluded that the two-phase flow model with the aid of the advanced VOF technique can provide with acceptably accurate numerical results on the route to practical purposes.     

Experimental study of the working principal and efficiency of a superposed inclined planes wave absorber

The swell reflection leads to problems in full scale conditions, shipping near seawalls and inside harbors, and for experiments in basins or channels, as well. Consequently, many studies were carried out to define wave absorbers with best efficiency. In this paper, the three main working principals of the wave absorbers, breaking waves, viscous dissipation resonating mechanism, are detailed to see how it is possible to enhance the efficiency. Through this analysis, a wave absorber made up of several superposed inclined planes in front of a wall is defined and its working principal is explained. Through experiments, it is shown that in spite of its short size its efficiency is one of the best compared to the more classical wave absorbers. Moreover, this superposed inclined planes wave absorber presents two main advantages: its short size and the possibility to be adjusted to the swell caracteristics, which make it interesting to be used in real scale and in laboratory environment, as well.     

基于波形不对称性的波浪破碎指标

为了探寻波浪破碎与波形不对称性的关系,通过对1/200缓坡上波浪破碎实验研究结果的进一步分析,运用最小二乘法,拟合了波形不对称性参数与相对水深的关系,以及用波形不对称性参数表示的波浪破碎指标表达式。所得规则波的结果与Kjeldsen的深水波结果相同,而不规则波的结果比规则波的小。研究还表明,这一破碎指标与相对水深有关系,随着水深变浅,指标值增大。     

Testing and calibrating parametric wave transformation models on natural beaches

To provide coastal engineers and scientists with a detailed inter-comparison of widely used parametric wave transformation models, several models are tested and calibrated with extensive observations from six field experiments on barred and unbarred beaches. Using previously calibrated (“default”) values of a free parameter γ, all models predict the observations reasonably well (median root-mean-square wave height errors are between 10% and 20%) at all field sites. Model errors can be reduced by roughly 50% by tuning γ for each data record. No tuned or default model provides the best predictions for all data records or at all experiments. Tuned γ differ for the different models and experiments, but in all cases γ increases as the hyperbolic tangent of the deep-water wave height, H_o. Data from two experiments are used to estimate empirical, universal curves for γ based on H_o. Using the new parameterization, all models have similar accuracy, and usually show increased skill relative to using default γ.