Discussion of “A simple method to determine breaker height and depth for different deepwater wave height/length ratios and sea floor slopes”, by J.P. Le Roux [Coastal Engineering 54 (2007) 271–277]

In the recent paper by J.P. Le Roux [Coastal Engineering 54 (2007) 271–277], the author provides a simplified approach to calculating the depth, length, and height of waves at the onset of depth-induced breaking (i.e. at the breaker line). However, the proposed methodology and the comparisons to other methods suffer from a large number of inconsistencies and basic calculation errors. In addition, there are a number of erroneous physical interpretations and many of the conclusions are based on erroneous data. The remaining conclusions are either not new or based on circular logic, such as to render them moot. In the following, we will not attempt to point out all the errors or inconsistencies that we found, instead we focus on major points of contention.     

Shear and Richardson number in a mode-water eddy

Measurements of stratification and shear were carried out as part of the EDDIES tracer release experiment in mode-water eddy A4 during the summer of 2005. These measurements were accomplished using both shipboard instrumentation and a drifting mooring. A strong relationship between shear intensity and distance from the center of the eddy A4 was observed with the shipboard ADCP. Diapycnal diffusivity at the SF₆ tracer isopycnal prior to and during the release was estimated from the drifting mooring to be 2.9×10⁻⁶ m² s⁻¹. Diffusivity increased by an order of magnitude to 3.2×10⁻⁵ m² s⁻¹ during the period of the final tracer survey in early September, which was similar to the value estimated from the tracer analysis for the whole experiment (3.5×10⁻⁵ m² s⁻¹, [Ledwell, J.R., McGillicuddy Jr., D.J., Anderson, L.A., 2008. Nutrient flux into an intense deep chlorophyll layer in a mode-water eddy. Deep-Sea Research II, this issue [doi:10.1016/j.dsr2.2008.02.005]].     

The statistical distribution of nearbed wave orbital velocity in intermediate coastal water depth

The statistical distribution of wave orbital velocity in intermediate coastal water depth has been quantitatively determined from the comprehensive field velocity data collected near the seabed in this study. Two ocean ADV current meters, which were mounted at 0.5 m above the seabed on two separate stainless steel tripods sitting on the seabed, were used to measure instantaneous water particle velocities at a 2 Hz sampling rate for 17.07 min every hour in two coastal water depths of 11 m and 23 m in nine field deployments over a period of 2 years. The zero-crossing method is applied to analyse the field velocity data collected in each field deployment to obtain a large sample of wave orbital velocity amplitudes of individual waves. Based on the collected field velocity data, it is found that the histogram of instantaneous wave orbital velocities perfectly follows the Gaussian distribution as commonly assumed, while the histogram of wave orbital velocity amplitudes is less accurately described by the Rayleigh distribution than the modified Rayleigh and the Weibull distribution. It is also found that large orbital velocity amplitudes are generally overestimated by the Rayleigh distribution, but well predicted by the modified Rayleigh and the Weibull distribution. The expected value of maximum orbital velocity in a velocity record of finite size is also derived from the three distributions and found to agree well with the present field data.     

Modified Moving Particle Semi-implicit methods for the prediction of 2D wave impact pressure

As a gridless particle method, the MPS (Moving Particle Semi-implicit) method has proven useful in a wide variety of engineering applications including free-surface hydrodynamic flows. Despite its wide range of applicability, the MPS method suffers from some shortcomings such as non-conservation of momentum and spurious pressure fluctuation. By introducing new formulations for the pressure gradient and a new formulation of the source term of the Poisson Pressure Equation (PPE), and by allowing a slight compressibility, we have proposed modified MPS methods for the prediction of wave impact pressure on a coastal structure. The improved performance of the modified methods is shown through the simulation of numerous wave impact problems (including the impacts by a dam break flow, a flip-through and two cases of slightly-breaking waves) in comparison with the experimental data.     

Characteristics of developing waves as a function of atmospheric conditions,water properties,fetch and duration

For any specific wind speed, waves grow in period, height and length as a function of the wind duration and fetch until maximum values are reached, at which point the waves are considered to be fully developed. Although equations and nomograms exist to predict the parameters of developing waves for shorter fetch or duration conditions at different wind speeds, these either do not incorporate important variables such as the air and water temperature, or do not consider the combined effect of fetch and duration. Here, the wind conditions required for a fully developed sea are calculated from maximum wave heights as determined from the wind speed, together with a published growth law based on the friction velocity. This allows the parameters of developing waves to be estimated for any combination of wind velocity, fetch and duration, while also taking account of atmospheric conditions and water properties.     

Higher harmonics induced by a submerged horizontal plate and a submerged rectangular step in a wave flume

The decomposition of a monochromatic wave over a submerged plate is investigated experimentally in a wave flume. Bound and free higher harmonic modes propagating upstream and downstream the structure are discriminated by means of moving resistive probes. The first-order analysis shows a resonant behaviour linked to the ratio of the plate's width and the fundamental mode wavelength over the plate. The second-order analysis shows an energy transfer from the fundamental mode towards free harmonics propagating downstream the structure. This transfer is linked to the ratio of the width of the plate and the bound harmonic wavelength over the plate. We also performed experiments with a submerged step to compare the efficiency of both structures. The submerged plate is shown to be a more efficient breakwater than the step, at the first as well as the second-order.     

Large-scale laboratory measurements of solitary wave inundation on a 1:20 slope

The run-up flow and related pressure of solitary waves breaking on a 1:20 plane beach were investigated experimentally in a super tank (300 m × 5 m × 5.2 m). Swash flow measurements of flow velocities are compared with an existing analytical solution. By incorporating an analytical solution, the hydrodynamic pressure for a quasi-steady flow state is determined and compared with laboratory data. Concerning the evident extra pressure exerted by the impact of swash flow, an empirical drag coefficient for a circular plate is also suggested in the present study.     

Prototype measurements and small-scale model tests of wave overtopping at shallow rubble-mound breakwaters: the Ostia-Rome yacht harbour case

The paper presents the comparison between the results of small-scale model tests and prototype measurements of wave overtopping at a rubble-mound breakwater. The specific structure investigated is the west breakwater of the yacht harbour of Rome at Ostia (Italy) and is characterized by a gentle seaward slope (1/4) and by a long, shallow foreshore. The laboratory tests firstly aimed at carefully reproducing two measured storms in which overtopping occurred and was measured. The tests have been carried out in two independent laboratories, in a wave flume and in a wave basin, hence using a two-dimensional (2-D) and a three-dimensional (3-D) setup. In the 2-D laboratory tests no overtopping occurred during the storm reproductions; in the 3-D case discharges five to ten times smaller than those observed in prototype have been measured. This indicates the existence of model and scale effects. These effects have been discussed on the basis of the results of several parametric tests, which have been carried out in both laboratories, in addition to the storm reproductions, varying wave and water level characteristics. Final comparison of all the performed tests with 86 prototype measurements still suggests the existence of scale and model effects that induce strong underestimation of overtopping discharge at small scale. The scale reproduction of wave breaking on the foreshore, together with the 3-D features of the prototype conditions and the absence of wind stress in the laboratory measurements, have been individuated as the main sources of scale and model effects. The paper also provides a comparison between the data and a largely used formula for wave overtopping discharges in the presence of structures similar to the one at hand. The suitable value of a roughness factor that appears in that formula is investigated and good agreement is found with other recent researches on rubble-mound breakwaters.     

白浪覆盖率模式中参量的研究

在波面位移为正态过程的假定下,推导出一种以平均周期和风速为参量的白浪覆盖率公式W=1-Φ[5.11094[-T/U10]0.7576].依据摩擦风速和U10的表达式,进一步推导出白浪覆盖率依赖于摩擦风速的形式W=1-Φ[0.5227[-T/U]0.7576]].考虑到在实际应用中,经常需要用波龄描述波浪的状态,将白浪覆盖率表示成一种形式简单的波龄的函数W=1-Φ(3.6496ξ0.7576),与Monahan等的海上测量数据符合良好.     

Numerical simulation of nearshore circulation on field topography under random wave environment

In order to investigate surf zone hydrodynamics through two-dimensional numerical simulations of nearshore circulation under random wave environment, a nearshore circulation model, SHORECIRC, and a random wave model, SWAN, were combined and utilized. Using this combined model, a numerical simulation of the October 2, 1997 SandyDuck field experiment was performed. For this simulation, field topography and an input offshore spectrum were constructed using observed data sets synchronized with the experiment. The wave-breaking model in SWAN was modified by using breaker parameters varied according to bottom slope. The simulation results were compared with the experimental data, which revealed a well-developed longshore current, as well as with results using other combinations which were SHORECIRC and its original monochromatic wave-driver, and SHORECIRC and the default of SWAN. The results from the novel combined model agreed well with the experimental data. The results of the present simulation also indicate that alongshore field topography influences shear fluctuation of longshore currents.