Interaction of flexural gravity waves with shear current in shallow water

In the present study, the effect of shear current on the propagation of flexural gravity waves is analyzed under the assumptions of linearized shallow-water theory. Explicit expressions for the reflection and transmission coefficients associated with flexural gravity wave scattering by a step discontinuity in both water depth and current speed are derived. Further, trapping and scattering of flexural gravity waves by a jet-like shear current with a top-hat profile are examined and certain limiting conditions for the waves to exist are derived. The effects of change in water depth, current speed, incident wavelength and the angle of incidence on the group and phase velocities as well as on the reflection and transmission characteristics are analyzed through different numerical results.     

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.     

Velocity potential formulations of highly accurate Boussinesq-type models

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.     

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.     

Formula for predicting bedload transport rate in oscillatory sheet flows

The method of Wang [Wang, Y.-H., 2007. Formula for predicting bedload transport rate in oscillatory sheet flow. Coastal Engineering 54, 594–601] to predict the wave friction factor is discussed. At the threshold of sediment entrainment the proposed equation produces a wide scatter of data points when plotted against an equation based on the Shields parameter. It is shown that a better correlation coefficient can be obtained by calculating the critical wave friction factor directly from the wave period, as well as the sediment and water properties.     

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.     

A temporally second-order accurate Godunov-type scheme for solving the extended Boussinesq equations

A numerical scheme for solving the class of extended Boussinesq equations is presented. Unlike previous schemes, where the governing equations are integrated through time using a fourth-order method, a second-order Godunov-type scheme is used thus saving storage and computational resources. The spatial derivatives are discretised using a combination of finite-volume and finite-difference methods. A fourth-order MUSCL reconstruction technique is used to compute the values at the cell interfaces for use in the local Riemann problems, whilst the bed source and dispersion terms are discretised using centred finite-differences of up to fourth-order accuracy. Numerical results show that the class of extended Boussinesq equations can be accurately solved without the need for a fourth-order time discretisation, thus improving the computational speed of Boussinesq-type numerical models. The numerical scheme has been applied to model a number of standard test cases for the extended Boussinesq equations and comparisons made to physical wave flume experiments.     

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.