Numerical simulation of random wave slamming on structures in the splash zone

The numerical investigation of random wave slamming on superstructures of marine structures in the splash zone is presented in this paper. The impact pressures on the underside of the structure are computed based on the improved volume of fluid method (VOF). The governing equations are Reynolds time-averaged equations and the two equation k– model. The third order upwind difference scheme is applied to the convection term to reduce the effect of numerical viscosity. The numerical wave flume with random wave-maker suitable for VOF is established. Appropriate moving contact-line boundary conditions are introduced to the model wave in contact with and separated from the underside of structure. Parametric studies have been carried out for different incident waves, structure dimensions and structure clearance. The numerical results are verified by the experimental results.     

A ten-year data set for fetch- and depth-limited wave growth

This paper presents the key results from a ten-year data set for Lake IJssel and Lake Sloten in The Netherlands, containing information on wind, storm surges and waves, supplemented with SWAN 40.51 wave model results. The wind speeds U₁₀, effective fetches x and water depths d for the data set ranged from 0–24 m s^− 1, 0.8–25 km and 1.2–6 m respectively. For locations with non-sloping bottoms, the range in non-dimensional fetch x? ( = gxU₁₀^− 2) was about 25–80,000, while the range in dimensionless depth d? ( = g d U₁₀^− 2) was about 0.03–1.7. Land–water wind speed differences were much smaller than the roughness differences would suggest. Part of this seems due to thermal stability effects, which even play a role during near-gale force winds. For storm surges, a spectral response analysis showed that Lake IJssel has several resonant peaks at time scales of order 1 h. As for the waves, wave steepnesses and dimensionless wave heights H? ( = gH_m0U₁₀^− 2) agreed reasonably well with parametric growth curves, although there is no single curve to which the present data fit best for all cases. For strongly depth-limited waves, the extreme values of d? (0.03) and H_m0 / d (0.44) at the 1.7 m deep Lake Sloten were very close to the extremes found in Lake George, Australia. For the 5 m deep Lake IJssel, values of H_m0 / d were higher than the depth-limited asymptotes of parametric wave growth curves. The wave model test cases of this study demonstrated that SWAN underestimates H_m0 for depth-limited waves and that spectral details (enhanced peak, secondary humps) were not well reproduced from H_m0 / d = 0.2–0.3 on. SWAN also underestimated the quick wave response (within 0.3–1 h) to sudden wind increases. For the remaining cases, the new [Van der Westhuysen, A.J., Zijlema, M., and Battjes, J.A., 2007. Nonlinear saturation-based whitecapping dissipation in SWAN for deep and shallow water, Coast. Eng., 54, 151–170] SWAN physics yielded better results than the standard physics of Komen, G.J., Hasselmann, S., Hasselmann, K., 1984. On the existence of a fully developed wind-sea spectrum. J. Phys. Oceanogr. 14, 1271–1285, except for persistent overestimations that were found for short fetches. The present data set contains many interesting cases for detailed model validation and for further studies into the evolution of wind waves in shallow lakes.     

Experimental study of irregular wave impact on structures in the splash zone

The experimental investigation of irregular wave slamming on structure members with large dimension in the splash zone is presented. The time-domain and frequency-domain analyses results of the irregular wave impact pressure on the subface of the structure under various case studies are presented. The influence of different incident wave parameters and relative clearance s/H_1/3 on the impact pressure is discussed. The empirical formula of the characteristic impact pressure P_c, P_t and the spectral moment of the impact pressure m₀ on the structure are proposed. The empirical expression of the characteristic peak impact pressure P_c and the spectral moment of the impact pressure m₀ on the structure is given.     

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.     

Profiles of fully developed (Airy) waves in different water depths

Airy waves have a sinusoidal profile in deep water that can be modeled by a time series at any point x and time t, given by η(x,t) = (H_o/2) cos[2πx/L_o − 2πt/T_w], where H_o is the deepwater height, L_o is the deepwater wavelength, and T_w is the wave period. However, as these waves approach the shore they change in form and dimension so that this equation becomes invalid. A method is presented to reconstruct the wave profile showing the correct wavelength, wave height, wave shape, and displacement of the water surface with respect to the still water level for any water depth.     

Characteristic friction coefficient and scale effects in oscillatory porous flow

This work presents a simple method to evaluate the performance of a porous breakwater when it is impinged with normal incidence by a non-breaking monochromatic wave train. It is based on: 1) a potential flow model for wave interaction with permeable structures and 2) a set of experimental tests on a rectangular porous structure with uniform granular distribution. A characteristic friction diagram is obtained considering wave energy balance in a control volume, minimising the error between the numerical model and the experimental results for the wave transmission coefficient. Results show that, for large breakwater widths, the reflection process reaches a saturation regime before the waves exit the structure at a distance from the seaside between the interval 0.2 < x/L < 0.45. For larger breakwater widths, the reflection coefficient is almost constant (except for “resonant” conditions) and wave transmission decreases exponentially. Under such conditions, the wave propagation through the porous medium depends on the relative diameter D/L and the porosity of the material; the dependence on the relative breakwater width B/L and the ratio diameter wave height D/H is weak. This diagram intends to be useful for preliminary engineering studies of breakwater's efficiency and performance and as an adequate selection criteria of the experimental stone diameter to minimize scale effects in laboratory studies.     

The effect of temporal wave averaging on the performance of an empirical shoreline evolution model

The effect of using time-averaged wave statistics in a simple empirical model for shoreline change is investigated. The model was first calibrated with a six-year time series of hourly wave conditions and weekly shoreline position at the Gold Coast, Australia. The model was then recalibrated with the hourly waves averaged over intervals up to 1 year. With wave averaging up to 2 days, model performance was approximately constant (squared correlation r² ~ 0.61–0.62), with only small changes in the values of empirical model parameters (e.g. the beach response coefficient c varied by less than 4%). With between 2 and 40 day averaging, individual storms are not resolved; model skill decreased only modestly (r² ~ 0.55), but c varied erratically by up to 40% of the original value. That is, optimal model coefficients depend on wave averaging, an undesirable result. With increased averaging (> 40 days) seasonal variability in the wave field is not resolved well and model skill declined markedly. Thus, temporal averaging of wave conditions increases numerical efficiency, but over-averaging degrades model performance and distorts best-fit values of model free parameters.     

Field and laboratory measurements of mean overtopping discharges and spatial distributions at vertical seawalls

Within the CLASH project, wave overtopping at the vertical seawall at Samphire Hoe was measured by HR Wallingford (HRW), and compared laboratory tests in 2 & 3 dimensions carried out at the University of Edinburgh and HRW. At Samphire Hoe, overtopping volumes were captured in three volumetric tanks capable of measuring wave-by-wave and total overtopping volumes. The three tanks were placed progressively farther back from the seawall edge so that the spatial distribution of the overtopping discharges could be determined. The field measurement equipment was successfully deployed on three occasions, and measured overtopping discharges ranged from that barely considered to be hazardous to the public to over q = 3.0 l/s/m. The 2d testing at Edinburgh was modelled at a scale of 1:40, and the 3d model at HRW was modelled at 1:20. For both sets of laboratory tests, a range of conditions, representative of the storm wave conditions and water levels, was reproduced in addition to a set of parametric conditions. The storm conditions allowed a direct comparison between the field and laboratory measurements, and the parametric conditions were used to test the generic overtopping behaviour of the structure. For both sets of laboratory tests, mean overtopping discharges and the spatial distribution were measured separately. Analysis of the distribution data relates the proportion of the discharge that has landed as a function of (x / L_o); where x is the distance behind the crest, and L_o is the offshore wavelength. Analysis of the field, 2d & 3d laboratory data, and empirical prediction methods have not identified any scale effects for overtopping discharges at vertical and near-vertical seawalls.     

Wave boundary layer over a stone-covered bed

This paper summarizes the results of an experimental investigation on wave boundary layers over a bed with large roughness, simulating stone/rock/armour block cover on the sea bottom. The roughness elements used in the experiments were stones the size of 1.4cm and 3.85cm in one group of experiments and regular ping-pong balls the size 3.6cm in the other. The orbital-motion-amplitude-to-roughness ratio at the bed was rather small, in the range a/k_s = 0.6–3. The mean and turbulence properties of the boundary-layer flow were measured. Various configurations of the roughness elements were used in the ping-pong ball experiments to study the influence of packing pattern, packing density, number of layers and surface roughness of the roughness elements. The results show that the friction factor seems to be not extremely sensitive to these factors. The results also show that the friction factor for small values of the parameter a/k_s does not seem to tend to a constant value as a/k_s → 0 (contrary to the suggestion made by some previous investigators). The present friction-factor data indicates that the friction factor constantly increases with decreasing a/k_s. An empirical expression is given for the friction factor for small values of a/k_s. The results further show that the phase lead of the bed friction velocity over the surface elevation does not seem to change radically with a/k_s, and found to be in the range 12°–23°. Furthermore the results show that the boundary-layer turbulence also is not extremely sensitive to the packing pattern, the packing density, the number of layers and the surface roughness of the roughness elements. There exists a steady streaming near the bed in the direction of wave propagation, in agreement with the existing work. The present data indicate that the steady streaming is markedly smaller in the case where the ping-pong balls are aligned at 45° to the wave direction than in the case with 90° alignment. Likewise, it is found that the steady streaming is relatively smaller in the case of the one-layer ping-pong-ball roughness than in the case of the two-layer situation.     

Effect of bed roughness on turbulent boundary layer and net sediment transport under asymmetric waves

This paper presents an investigation of the roughness effects in the turbulent boundary layer for asymmetric waves by using the baseline (BSL) k–ω model. This model is validated by a set of the experimental data with different wave non-linearity index, N_i (namely, N_i = 0.67, N_i = 0.60 and N_i = 0.58). It is further used to simulate asymmetric wave velocity flows with several values of the roughness parameter (a_m/k_s) which increase gradually, namely from a_m/k_s = 35 to a_m/k_s = 963. The effect of the roughness tends to increase the turbulent kinetic energy and to decrease the mean velocity distribution in the inner boundary layer, whereas in the outer boundary layer, the roughness alters the turbulent kinetic energy and the mean velocity distribution is relatively unaffected. A new simple calculation method of bottom shear stress based on incorporating velocity and acceleration terms is proposed and applied into the calculation of the rate of bed-load transport induced by asymmetric waves. And further, the effect of bed roughness on the bottom shear stress and bed-load sediment transport under asymmetric waves is examined with the turbulent model, the newly proposed method, and the existing calculation method. It is found that the higher roughness elements increase the magnitude of bottom shear stress along a wave cycle and consequently, the potential net sediment transport rate. Moreover, the wave non-linearity also shows a big impact on the bottom shear stress and the net sediment transport.     

Effective fetch and non-linear four-wave interactions during wave growth in slanting fetch conditions

Wave growth in slanting fetch (with wind blowing obliquely off a coast) is investigated with 7 years worth of routine wave measurements in Lake IJssel in The Netherlands and with the SWAN wave model. Two aspects are considered in particular for this case: the validity of the concept of effective fetch and the role of the non-linear four-wave interactions. For slanting and parallel fetch conditions, we found some significant deviations from the effective fetch assumption, leading to 20–35% mismatch in either the peak period T_p or the significant wave height H_m0 respectively. However, the effect of discrepancies between various widely accepted wave growth formulas turned out to be even more important. The wave directions during slanting fetch are significantly ‘steered’ by the coastline, especially in the first kilometre(s) off the coast. The role of the non-linear four-wave interactions is investigated by running the SWAN (version 40.41) wave model with three different quadruplet formulations. Exact quadruplet methods (Xnl) yielded relatively strong wave steering, despite the four-wave interactions being relatively weak. Application of Xnl did not lead to better overall agreement with measurements — improvements for the mean wave period T_m01 were offset by some deterioration for the wave height H_m0.     

Wave interaction with a submarine pipeline in clayey soil due to random waves

The wave pressure and uplift force due to random waves on a submarine pipeline (resting on bed, partially buried and fully buried) in clayey soil are measured. The influence of various parameters viz., wave period, wave height, water depth, burial depth and consistency index of the soil on wave pressures around and uplift force on the submarine pipeline was investigated. The wave pressures were measured at three locations around the submarine pipeline (each at 120° to the adjacent one). It is found that the wave pressure and uplift force spectrum at high consistency index of the soil is smaller compared to that of low consistency index. Just burying the pipeline (e/D=1.0) in clayey soil reduces the uplift force to less than 60% of the force experienced by a pipeline resting on the seabed (e/D=0.0) for I_c=0.33.     

Observations of wave pump efficiency

Currently available data on wave pump efficiency is reviewed. The obtainable efficiency is an important consideration in the design of practical devices for the extraction of wave energy and the analysis of natural systems (e.g., coral flats and rip currents). We find that the peak efficiency is 0.5 for very steep (∼ 40–45°) ramps where the waves break over the top of the ramp. For flatter (< 30°) ramps, the breaking process is more gradual and the peak efficiency is less than 0.1. We have identified natural atoll lagoon systems where the flushing is wave driven and successfully modeled it as driven by a wave pump. The same is the case for rip currents. For both of these natural systems, the pump efficiency is around 0.035. In addition a numerical swash model is used to estimate wave pump efficiency and is seen to match the experimental results for natural systems or breaking wave scenario.     

Combined classifier–quantifier model: A 2-phases neural model for prediction of wave overtopping at coastal structures

A 2-phases neural prediction method for wave overtopping is developed. The ‘classifier’ predicts whether overtopping occurs or not, i.e. q = 0 or q > 0. If the classifier predicts overtopping q > 0, then the ‘quantifier’ is used to determine the mean overtopping discharge. The overtopping database set up within the EC project CLASH (De Rouck, J., Geeraerts, J., 2005. CLASH Final Report, Full Scientific and Technical Report, Ghent University, Belgium) is used to train the networks of the prediction method.     

Discussion of “Analytic solution of long wave propagation over a submerged hump” by Niu and Yu (2011)

Recently, Niu and Yu (2011) presented an analytical solution of the long wave refraction by a submerged circular hump. The geometry of the hump was assumed to be axi-symmetric and the water depth over the hump region was described by a positive constant plus a power function of the radial distance with an arbitrary value of the power exponent, i.e., h = h₁ + βr^s, where h₁ is the water depth at the crest of the hump. Their general hump is an extension of the paraboloidal hump (i.e., s = 2) studied by Zhang and Zhu (1994) and Zhu and Harun (2009). Because of this extension in the topography of the hump, the problem to seek a general analytical solution to the long-wave equation becomes much more complicated and the solution technique need to be more skillful, especially for the case with the exponent s being a rational, see Eq. (17) in Niu and Yu (2011).     

Alternative forms of the higher-order Boussinesq equations: Derivations and validations

An alternative form of the Boussinesq equations is developed, creating a model which is fully nonlinear up to O(μ⁴) (μ is the ratio of water depth to wavelength) and has dispersion accurate to the Padé [4,4] approximation. No limitation is imposed on the bottom slope; the variable distance between free surface and sea bottom is accounted for by a σ-transformation. Two reduced forms of the model are also presented, which simplify O(μ⁴) terms using the assumption ε = O(μ^2/3) (ε is the ratio of wave height to water depth). These can be seen as extensions of Serre's equations, with dispersions given by the Padé [2,2] and Padé [4,4] approximations. The third-order nonlinear characteristics of these three models are discussed using Fourier analysis, and compared to other high-order formulations of the Boussinesq equations. The models are validated against experimental measurements of wave propagation over a submerged breakwater. Finally, the nonlinear evolution of wave groups along a horizontal flume is simulated and compared to experimental data in order to investigate the effects of the amplitude dispersion and the four-wave resonant interaction.     

Sub aerial random wave pressure and layer thickness on impermeable sea walls

Investigations on sub aerial wave pressures and layer thickness on plane impermeable and non-overtopping seawallns were carried out by using physical model studies. Seawalls with slopes of 1:3, 1:4 and 1:6 were used. JONSWAP spectrum with significant wave height, H_s from 0.08 to 0.2 m and peak periods, T_p from 1.5 to 6.0 s and a constant water depth of 0.7 m is used. Based on extensive measurements, empirical formulas for practical applications are proposed to predict the maximum, significant and mean sub aerial random wave pressure and layer thickness (thickness of water layer perpendicular to the still water level on the run-up zone) by using the surf similarity parameter, significant wave height and elevation on the sub aerial region as inputs. It is found that the maximum layer thickness is 1.11 times the significant layer thickness and maximum sub Arial wave pressure is 1.06 times the significant wave pressures. The predictive equations based on extensive measurements can be used for the design of non-overtopping seawalls.     

斜向波浪作用下双层水平板式防波堤波浪荷载试验研究

In this study, systematic physical model tests were performed to investigate the wave forces on the twin-plate breakwater under irregular waves. Based on the experimental results, the effects of the relative plate width B/L,wave height Hs/D and incident angle θ0 on the wave forces were analyzed and discussed. The results showed that:(1) The envelopes of the total wave pressure were generally symmetrical along the direction of plate width under the incident angles(θ0) being 0°, 15°, 30°, 45° and 60°. In particular, the envelopes of wave pressure atθ0=30° were larger than all other cases.(2) The synchronous pressure distribution of the breakwater under oblique wave action was more complicated comparing to the normal incident waves.(3) Based on data analysis, an empirical formula was obtained to estimate the total vertical force of the twin-plate breakwater.This empirical formula can be a good reference for the design basis of engineering applications under specified wave conditions.     

Variability of wave-induced ripple migration in wave-flume experiments and its implications for sediment transport

A thorough discussion of results from laboratory experiments with regular waves sheds light on the gap that lies between the sediment transport associated with ripple migration and the performance of a standard bedload transport formula in terms of bed shear concept. It is found that the extent of deviations of the bedload transport formula by Ribberink (1998) from the measured rate of sediment transport associated with ripple migration becomes systematically apparent under conditions of increasing settling time factor Ω_s (= η/(w₀T); η is the ripple height, w₀ the settling velocity and T the wave period). Re-examination of previous two field studies demonstrates a further reinforcement for phase-lag argument addressed in this paper.     

Improvement of sand activation depth prediction under conditions of oblique wave breaking

This study presents sand activation depth (SAD) measurements recently obtained on two contrasting beaches located along the Atlantic coast of France: the gently sloping, high-energy St Trojan beach where wave incidence is usually weak, and the steep, low-energy Arçay Sandspit beach where waves break at highly oblique angles. Comparisons between field measurements and predictions from existing formulae show good agreement for St Trojan beach but underestimate the SAD on the Arçay Sandspit beach by 40–60%. Such differences suggest a strong influence of wave obliquity on SAD. To verify this hypothesis, the relative influence of wave parameters was investigated by means of numerical modelling. A quasi-linear increase of SAD with wave height was confirmed for shore-normal and slightly oblique wave conditions, and a quasi-linear increase in SAD with wave obliquity was also revealed. Combining the numerical results with previously published relations, both a new semi-empirical and an empirical formula for the prediction of SAD were developed which showed good SAD predictions under conditions of oblique wave breaking. The new empirical formula for the prediction of SAD (Z ₀) takes into account the significant wave height (H _s), the beach face slope (β) and the wave angle at breaking (α), and is of the form $ Z_{0} = 1.6\tan {\left( \beta \right)}H^{{0.5}}_{{\text{s}}} {\sqrt {1 + \sin {\left( {2\alpha } \right)}} } This study presents sand activation depth (SAD) measurements recently obtained on two contrasting beaches located along the Atlantic coast of France: the gently sloping, high-energy St Trojan beach where wave incidence is usually weak, and the steep, low-energy Ar?ay Sandspit beach where waves break at highly oblique angles. Comparisons between field measurements and predictions from existing formulae show good agreement for St Trojan beach but underestimate the SAD on the Ar?ay Sandspit beach by 40–60%. Such differences suggest a strong influence of wave obliquity on SAD. To verify this hypothesis, the relative influence of wave parameters was investigated by means of numerical modelling. A quasi-linear increase of SAD with wave height was confirmed for shore-normal and slightly oblique wave conditions, and a quasi-linear increase in SAD with wave obliquity was also revealed. Combining the numerical results with previously published relations, both a new semi-empirical and an empirical formula for the prediction of SAD were developed which showed good SAD predictions under conditions of oblique wave breaking. The new empirical formula for the prediction of SAD (Z ₀) takes into account the significant wave height (H _s), the beach face slope (β) and the wave angle at breaking (α), and is of the form . The use of a dataset from the literature demonstrates the predictive skill of these new formulae for a wide range of wave heights, wave incidence and beach gradients.