Extreme wave crest distribution by direct numerical simulations of long-crested nonlinear wave fields

The wave crest height qualification checks are required during the wave calibration before the model test in wave basin. However, the reliable criteria of nonlinear wave crest probability distribution in 3-h duration (full-scale) has not been well established yet. We investigate wave crest-height statistics of long-crested nonlinear wave fields using high-order spectral (HOS) method, which can take the effects of both second-order bound waves and third-order free waves into account. The energy dissipation effects due to wave breaking were included by employing an eddy viscosity model. Sensitivity analyses to the wave breaking onset criterion have been performed. Validation is provided by comparing the obtained numerical results with the available calibration test data. Based on extensive and direct numerical simulations, semi-empirical single realization distributions for wave calibration have been developed through 3-parameter Weibull fitting and systematic regression analyses. Particular attention has been paid to the tail of upper bound of wave crest distributions. The effects of wave steepness and water depth on the maximum wave crest height in 3-h duration have been examined. It is found that with the increase of wave steepness, the extreme wave crest height increases until it reaches a critical value. In addition, for the scale water depth k_ph < 1.36, the maximum crest height decreases as the water depth increases, while in the opposite case the maximum crest height increases as the water depth increases. Moreover, it is confirmed that that the fourth-order nonlinearity does not have significant effects on the distribution of the wave crest height.     

Laboratory flume studies on monochromatic wave-fine sandy bed interactions Part 2. Sediment suspensions

As reported in preceding paper (Part 1. Soil Fluidization), the observed phenomena of sediment suspensions above a fluidized sandy bed of Sand II (d₅₀ = 0.092 mm) under monochromatic wave actions are quantitatively investigated. The suspended sediment concentration (SSC) at a single point within 5 cm above the bed was synchronously measured with water waves and bed soil's pore pressures with an intrusive optical sediment-concentration probe. The measurements show that SSC initiates several wave cycles after initiation of bed soil's fluidized response and grows to a peak value mainly in the post-fluidization phase. Under similar wave loadings in the same test series, SSC is usually higher over a resonantly fluidized (RF) bed than over a non-resonantly fluidized (NRF) bed. On the contrary, only relatively low SCC can be identified above an unfluidized bed. The analyses illustrate that to certain extent, peak values of SSC are directly proportional to the thickness of fluidized soil layer d_f. Values of d_f usually decrease with repeated fluidized response, longer consolidation periods, and in deeper water depths. Once the fluidized responses initiate, pore pressures are generally much significantly amplified in both shallow fluidized soil layers and near below the fluidized layer, especially during the resonance event. The resulting depth gradients of dynamic pore pressure amplitudes in shallow layers are likely to have caused higher initial rises of SSC in a RF bed than in the subsequent NRF bed. Those in deeper layer should have contributed to sustain the fluidization state for further SSC increments. Immediately after termination of wave loading, re-deposited suspended sediments always result in a typical flat bed form. For a pre-fluidized bed, wave-induced drastic sediment suspensions are still obtainable very near above the bed with even a rather thin fluidized surface soil layer.     

风暴过程中潮滩悬沙浓度和悬沙输运的变化及其动力机制——以长江三角洲南汇潮滩为例

悬沙浓度是淤泥质海岸重要的环境指标。为探讨潮滩悬沙浓度和悬沙输运对风暴事件的响应过程及其动力机制,于2014年9月"凤凰"台风过境前、中、后在长江三角洲南汇潮滩进行了现场观测,获得同步高分辨率的水深、波高、近底流速和浊度剖面时间序列(9个潮周期)。结果表明,风暴中平均和最大波高、波-流联合底床剪切应力、悬沙浓度和悬沙输运率可比平静天气高数倍;风暴期间高潮位低流速阶段悬沙沉降导致近底发育数十厘米厚的浮泥层(悬沙浓度大于10 g/L)。研究认为风暴事件中淤泥质海岸悬沙浓度和悬沙输运的剧烈变化其根本动力机制是风暴把巨大能量传递给近岸水体,进而显著增大波-流联合底床剪切应力,导致细颗粒泥沙再悬浮。     

Statistical studies on the wave-form and maximum height of large Tsunamis

The tide-gauge records of large tsunamis are classified into three types, A, B and C. The “A” type record is made up of one or a few large waves near the wave front. The “B” type record consists of one or a few wave groups. The “C” type is the combination of the “A” and “B” types. The data used are; the Kamchatka Tsunami of Nov. 4, 1952, the Aleutian Tsunami of March 9, 1957, the Chilean Tsunami of May 22, 1960 and the Alaska Tsunami of March 28, 1964. The A type occurs mostly at isolated islands in the Pacific Ocean and occasionally at continental coasts. The B type is mostly distributed on the continental coast and along the island-arc. The distribution of the C type differs from tsunami to tsunami. The relation between the delay time of the maximum wave and the the travel time of the wave front is as follows:

1. For the wave of the A type and the head wave of C type, the delay time (t _D ) is constant for all travel times.
2. For the first wave group of B and C types, the delay time (T ₁) is constant or slow decreases with travel time. For the second and third wave groups of B and C types, the definite decrease of delay times (T ₂ andT ₃) with travel time is observed.

The height (h) of the maximum wave of A and C types decreases generally with travel time. The maximum wave height along the path between Kamchatka and Chile, however, shows the increase. For all wave groups the wave heights (H ₁,H ₂ andH ₃) of B and C types increases with travel time. Some speculations on the causes of these features are also presented.     

Wave modelling in the vicinity of submerged breakwaters

This paper describes a simple method for modelling wave breaking over submerged structures, with the view of using such modelling approach in a coastal area morphodynamic modelling system.A dominant mechanism for dissipating wave energy over a submerged breakwater is depth-limited wave breaking. Available models for energy dissipation due to wave breaking are developed for beaches (gentle slopes) and require further modifications to model wave breaking over submerged breakwaters.In this paper, wave breaking is split into two parts, namely: 1) depth-limited breaking modelled using Battjes and Janssen's (1978) theory [Battjes, J.A. and Jannsen, J.P.F.M. (1978). Energy loss and setup due to breaking of random waves. Proceedings of the 16th Int. Conf. Coast. Eng., Hamburg, Germany, pp. 569-587.] and 2) steepness limited breaking modelled using an integrated form of the Hasselmann's whitecapping dissipation term, commonly used in fully spectral wind–wave models. The parameter γ₂, governing the maximum wave height at incipient breaking (H_max = γ₂d) is used as calibration factor to tune numerical model results to selected laboratory measurements. It is found that γ₂ varies mainly with the relative submergence depth (ratio of submergence depth at breakwater crest to significant wave height), and a simple relationship is proposed. It is shown that the transmission coefficients obtained using this approach compare favourably with those calculated using published empirical expressions.     

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.     

Interface depth used in a two-layer model of nonlinear internal waves

A two-layer model includes three parameters: interface depth h ₁, upper layer density \(\rho_{1}\) , and lower layer density \(\rho_{2}\) . Many theoretical and laboratorial studies of internal waves, as well as most numerical models, are based on the two-layer assumption. However, these three parameters cannot be directly measured because a pycnocline in the real ocean has finite thickness, and the densities in both the mixed layer and the deep ocean are not constant. In the present study, seven different methods are used to determine the interface depth of the two-layer model and compared with the depth of maximum vertical displacement: the depth of maximum buoyancy frequency (Ν _max), the depth where the first mode eigenfunction has its maximum (Φ_max), the depth where the lowest mode temperature empirical orthogonal function has its maximum, the depth where either the two-layer Korteweg–de Vries (KdV) or Benjamin–Ono equation has closest coefficients with their continuously stratified counterparts, and the same KdV approach with stratification replaced by two idealized distributions. The multi-ship measurement conducted near the Luzon Strait is used for deep ocean comparison, and two measurements conducted in the east of Dongsha Atoll are used for shallow water comparison. The results show that in the deep ocean, the KdV approach with idealized type I stratification gives the interface closest to the depth of maximum vertical displacement. In shallow waters, the KdV approach agrees with the measurement best.     

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.     

Experimental study of unidirectional irregular wave slamming on the three-dimensional structure in the splash zone

The experimental investigation of unidirectional random wave slamming on the three-dimensional structure in the splash zone is presented. The experiment is conducted in the marine environment channel in the State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology. The test wave is unidirectional irregular wave. The experiments are carried out with perpendicular random waves (β=0°) and oblique random waves (β=15°, 30°, 45°), the significant wave heights H_1/3 ranging from 7.5 to 20 cm with 2.5 cm increment, the peak wave periods T_p ranging from 0.75 to 2.0 s with 0.25 s increment, and the clearance of the model with respect to the significant wave height s/H_1/3 ranging from 0.0 to 0.5 with 0.1 increment. The statistical analysis results of different test cases are presented. The statistical distribution characteristics of the perpendicular irregular wave impact pressures are compared with that of the oblique irregular wave on the underside of the structure. The effect of the wave direction β on the wave impact forces on the underside of the structure is determined. The relation between the impact forces and the parameters such as the significant wave height, the relative structure width and the relative clearance of the structure is also discussed.     

Wave climate off northern Norway

The wave climate off northern Norway is considered and the investigation is based on wave measurements made at Tromsøflaket by means of a waverider buoy during the years 1977–1981. Data quality of waverider measurements is briefly commented upon; however, more emphasis is given to an evaluation of the long-term representativity of the actual measuring period and to a procedure accounting approximately for a lack of representativity. The wave climate is presented in terms of a smoothed joint probability density function of the significant wave height, H_s, and the spectral peak period, T_p. Based on this distribution a consistent design curve in the H_s, T_p space is established.     

Experimental study on the performance of twin plate breakwater

The wave transmission characteristics and wave induced pressures on twin plate breakwater are investigated experimentally in regular and random waves.A total of twenty pressure transducers are fixed on four surfaces of twin plate to measure the wave induced dynamic pressures.The spatial distribution of dynamic wave pressure is given along the surface of the twin plate.The uplift wave force obtained by integrating the hydrodynamic pressure along the structure is presented.Discussed are the influence of different incident wave parameters including the relative plate width B /L,relative wave height /i H a and relative submergence depth s /a on the non-dimensional dynamic wave pressures and total wave forces.From the investigation,it is found that the optimum transmission coefficient,t K occurs around B /L 0.41 ~ 0.43,and the twin plate breakwater is more effective in different water depths.The maximum of pressure ratio decreases from 1.8 to 1.1 when the relative submergence depth of top plate is increased from 0.8to +0.8.     

Prediction of wave pressures on smooth impermeable seawalls

Simple prediction methods are proposed to estimate the wave induced pressures on smooth impermeable seawalls. Based on the physics of the wave structure interaction, the sloped seawall is divided into a total of five zones (zones 1, 2 and 3 during run-up (corresponding pressures are called as positive pressures) and zones 4 and 5 during run-down (corresponding pressures are called negative pressures)) (Fig. 1). Zone 1 (0<z<d−H_i/2), where the wave pressure is governed by the partial reflection and phase shift; Zone 2 (d−H_i/2<z<d), where the effect of wave breaking and turbulence is significant; Zone 3 (d<z<Run-up height), where the pressure is induced by the run-up water; Zone 4 (Run-down<z<d), where the wave pressure is caused by the run-down effect and Zone 5 (0<z<d-Run down), where the negative wave pressures are due to partial reflection and phase shift effects. Here d is the water depth at the toe of the seawall, H_i is the incident wave height and z is the vertical elevation with toe of the seawall as origin and z is positive upward. For wave pressure prediction in zones 1 and 5, the empirical formula proposed by Ahrens et al. (1993) to estimate the wave reflection and Sutherland and Donoghue's recommendations (1998) for the estimation of phase shift of the waves caused by the sloped structures are used. Multiple regression analysis is carried out on the measured pressure data and empirical formulas are proposed for zones 2, 3 and 4. The recommendations of Van der Meer and Breteler (1990) and Schüttrumpf et al. (1994) for the prediction of wave run-down are used for pressure prediction at zone 4. Comparison of the proposed prediction formulas with the experimental results reveal that the prediction methods are good enough for practical purposes. The present study also shows a strong relation between wave reflection, wave run-up, wave run-down and phase shift of waves on wave pressures on the seawalls.     

Observation of wave-enhanced turbulence in the near-surface layer of the ocean during TOGA COARE

Dissipation rate statistics in the near-surface layer of the ocean were obtained during the month-long COARE Enhanced Monitoring cruise with a microstructure sensor system mounted on the bow of the research vessel. The vibration contamination was cancelled with the Wiener filter. The experimental technique provides an effective separation between surface waves and turbulence, using the difference in spatial scales of the energy-containing surface waves and small-scale turbulence. The data are interpreted in the coordinate system fixed to the ocean surface. Under moderate and high wind-speed conditions, we observed the average dissipation rate of the turbulent kinetic energy in the upper few meters of the ocean to be 3–20 times larger than the logarithmic layer prediction. The Craig and Banner (J. Phys. Oceanogr. 24 (1994) 2546) model of wave-enhanced turbulence with the surface roughness length from the water side z₀ parameterized according to the Terray et al. (J. Phys. Oceanogr. 26 (1996) 792) formula z₀=cH_s provides a reasonable fit to the experimental dissipation profile, where z is the depth (defined here as the distance to the ocean surface), c≈0.6, and H_s is the significant wave height. In the wave-stirred layer, however, the average dissipation profile deviates from the model (supposedly because of extensive removing of the bubble-disturbed areas close to the ocean surface). Though the scatter of individual experimental dissipation rates (10-min averages) is significant, their statistics are consistent with the Kolmogorov's concept of intermittent turbulence and with previous studies of turbulence in the upper ocean mixed layer.     

Stability of rubble-mound breakwater using H50 wave height parameter

The prediction of rubble mound breakwaters' stability is one of the most important issues in coastal and maritime engineering. The stability of breakwaters strongly depends on the wave height. Therefore, selection of an appropriate wave height parameter is very vital in the prediction of stability number. In this study, H₅₀, the average of the 50 highest waves that reach the breakwater in its useful life, was used to predict the stability of the armor layer. First, H₅₀ was used instead of the significant wave height in the most recent stability formulas. It was found that this modification yields more accurate results. Then, for further improvement of the results, two formulas were developed using model tree.To develop the new formulas, two experimental data sets of irregular waves were used. Results indicated that the proposed formulas are more accurate than the previous ones for the prediction of the stability parameter. Finally, the proposed formulas were applied to regular waves and a wide range of damage levels and it was seen that the developed formulas are applicable in these cases as well.     

Internal wave band kinetic energy production: flat vs. sloping bottom

Velocity and temperature measurements obtained with acoustic Doppler current profilers and thermistor strings are used to evaluate the production of internal wave band kinetic energy mainly in the frequency band σ>15 cpd. Results from a flat 19 m deep, vigorous tidal environment in a shelf sea are compared with energy production in a bottom boundary layer above a continental slope. In the tidal environment, maximum production occurs in the near-bottom and near-surface layers. A distinct mid-depth maximum in KE production occurs during a period when wind speeds exceed 10 m s⁻¹ and significant wave height ∼2 m. At the same time, no significant changes in the along-shore current speed take place but the cross-shore current, generated by strong stratification, is weakened. This suggests a direct energy input from the wind via surface waves into the water column turbulence. Maximum kinetic energy production in the frequency band σ>1.9 cpd, thus including the semidiurnal tide, occurs at mid-depth when strong stratification is present. The overall magnitude of internal wave band kinetic energy production agrees well with independent dissipation estimates obtained from microstructure profilers. Above the sloping bottom, KE production is somewhat larger than observed in the shallow tidal environment, despite rms currents being ∼50% smaller and wind effects being small. Above the sloping bottom KE shear production was comparable to buoyancy production. The latter was negligible at the shelf sea site.     

Focused wave evolution using linear and second order wavemaker theory

In this paper, the evolution of focused waves using different paddle displacements (piston type) under laboratory conditions is presented. It is well known that in intermediate water depths, linear paddle displacements will generate spurious, free, sub and super harmonics. Thus, a second order correction to suppress these spurious free sub and super harmonics was used to generate the focused waves. The focused waves were generated in the laboratory using a linear superimposition principle, in which the wave paddle displacement is derived based on the sum of a number of sinusoidal components at discrete frequencies, whose phases are accordingly set to focus at a particular location. For this method of generation, the second order wave maker theory proposed by Schäffer [24] can be easily adopted and was used in the present study. Two different centre frequencies (f_c = 0.68 Hz and 1.08 Hz) with three different bandwidth ratios (Δf/f_c = 0.5, 0.75 and 1.0) were tested in a constant water depth, to consider both narrow and broadband spectra. These test cases correspond to wave focusing packets propagating in intermediate and deep water regions. Further, for each wave packet, two different amplitudes were considered in order to analyze non-breaking and breaking cases. Thus, by systematically generating the wave packets using the linear and second order paddle displacements, the analysis was carried out for the spectral and temporal evolution of selected long waves. The temporal evolution of the selected harmonics was analyzed using the Inverse Fast Fourier Transform (IFFT), to show the propagation of the spurious, free, long waves. Further, the variations in energy for the lower, higher and primary frequency ranges are reported for different wave paddle displacements. The analysis revealed that for the broadband spectrum the differences are more pronounced when using linear paddle displacements. We have also noticed a shift in focusing/breaking location and time (i.e. premature) due to the increase in crest height using linear displacements. The experiment data used in this paper has been provided as a supplementary, which can be used to validate the numerical models.     

Laboratory observation of boundary layer flow under spilling breakers in surf zone using particle image velocimetry

A boundary layer flow under spilling breakers in a laboratory surf zone with a smooth bottom is investigated using a high resolution particle image velocimetry (PIV) technique. By cross-correlating the images, oscillatory velocity profiles within a viscous boundary layer of O(1) mm in thickness are resolved over ten points. Using PIV measurements taken for an earlier study and the present study, flow properties in the wave bottom boundary layer (WBBL) over the laboratory surf zone are obtained, including the mean velocities, turbulence intensity, Reynolds stresses, and intermittency of coherent events. The data are then used to estimate the boundary layer thickness, phase variation, and bottom shear stress. It is found that while the time averaged mass transport inside the WBBL is onshore in the outer surf zone, it changes to offshore in the inner surf zone. The zero Eulerian mass transport occurs at h/h_b ≈ 0.92 in the outer surf zone. The maximum overshoot of the streamwise velocity and boundary layer thickness are not constant across the surf zone. The bottom shear stress is mainly contributed by the viscous stress through mean velocity gradient while the Reynolds stress is small and negligible. The turbulence level is higher in the inner surf zone than that in the outer surf zone, although only a slight increase of turbulent intensity is observed inside the WBBL from the outer surf zone to the inner surf zone. The variation of phase inside and outside the WBBL was examined through the spatial velocity distribution. It is found the phase lead is not constant and its value is significantly smaller than previous thought. By analyzing instantaneous velocity and vorticity fields, a remarkable number of intermittent turbulent eddies are observed to penetrate into the WBBL in the inner surf zone. The size of the observed large eddies is about 0.11 to 0.16 times the local water depth. Its energy spectra follow the − 5/3 slope in the inertial subrange and decay exponentially in the dissipation subrange.     

Prediction of significant wave height using regressive support vector machines

Wave parameters prediction is an important issue in coastal and offshore engineering. In this literature, several models and methods are introduced. In the recent years, the well-known soft computing approaches, such as artificial neural networks, fuzzy and adaptive neuro-fuzzy inference systems and etc., have been known as novel methods to form intelligent systems, these approaches has also been used to predict wave parameters, as well. It is not a long time that support vector machine (SVM) is introduced as a strong machine learning and data mining tool. In this paper, it is used to predict significant wave height (H_s). The data set used in this study comprises wave wind data gathered from deep water locations in Lake Michigan. Current wind speed (u) and those belonging up to six previous hours are given as input variables, while the significant wave height is the output parameter. The SVM results are compared with those of artificial neural networks, multi-layer perceptron (MLP) and radial basis function (RBF) models. The results show that SVM can be successfully used for prediction of H_s. Furthermore, comparisons indicate that the error statistics of SVM model marginally outperforms ANN even with much less computational time required.     

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 γ.