Exceedance probability for wave overtopping on a fixed deck

Detailed laboratory measurements were made of the instantaneous free surface elevation in front of a fixed deck and the instantaneous free surface elevation, velocity, and overtopping rate at the leading edge of the deck. The study showed that the exceedance probabilities for the normalized maximum instantaneous overtopping rate and the normalized overtopping volume were predicted by a simple exponential curve. The measured exceedance probability seaward of the deck compared well with the nonlinear theory of Kriebel and Dawson (Kriebel D.L., Dawson T.H., 1993. Nonlinearity in wave crest statistics. In: Proceedings Ocean Wave Measurement and Analysis. American Society of Civil Engineers, pp. 61–75). Conditional sampling of the crest heights seaward of the deck gave a normalized probability distribution similar to that of the maximum water level measured on the deck for each overtopping event. However, the values used to normalize each distribution were not the same.     

Dynamic properties of green water event in the overtopping of extreme waves on a fixed dock

A statistical model is developed to predict wave overtopping volume and rate of extreme waves on a fixed deck. The probability density function for the volume and rate of overtopping water are formulated based on the truncated Weibull distribution with the assumption of local sinusoidal profile for small amplitude waves. Sensitivity to the wave nonlinearity parameter and deck clearance is discussed. The statistical model is compared to laboratory data of the instantaneous free surface elevation measured in front of a fixed deck, and overtopping volume and overtopping rate measured at the leading edge of the deck. The statistical theory compared well with the measured exceedance probability seaward of the deck. The model prediction of the exceedance probability of deck overtopping gave qualitatively good agreement for large overtopping values.     

Numerical study of nonlinear freak wave impact underneath a fixed horizontal deck in 2-D space

Freak waves are extreme and unexpected surface waves with huge wave heights that may lead to severe damage to ships and offshore structures. However, few researches have been conducted to investigate the impact underneath fixed horizontal decks caused by freak waves. To study these phenomena, a 2-D numerical wave tank is built in which nonlinear freak waves based on the Peregrine breather solution are generated. As a validation, a regular-wave-induced underneath impact is simulated and compared to the existing experimental measurements. Then the nonlinear freak-wave-induced impact is investigate with different values of deck clearance above the mean free surface. In addition, a comparative simulation of a “large” regular wave based on the 2nd-order Stokes wave theory with the same crest height and wave length of the nonlinear freak wave is carried out to reveal the unique features of the nonlinear freak-wave-induced impact. By applying a fluid–structure interaction (FSI) algorithm in which the bottom deck and front side wall are simplified as Euler beams in 2-D and discretized by the finite element method (FEM), the hydroelastic effects are considered during the impact event. The vertical force acting underneath the bottom deck, the transversal force acting on the front side wall, the structural displacements of the elastic deck and wall are analyzed and discussed respectively, from which meaningful conclusions are drawn.     

基于Boussinesq水波模型的聚焦波模拟

基于最高导数为3阶的单层Boussinesq方程,建立了聚焦波的时域波浪计算模型。数值模型求解采用了预报?校正的有限差分法。对于时间差分格式,预报和校正分别采用3阶Adams-Bashforth格式和4阶Adams-Moulton格式。首先,针对不同水深条件下水槽中传播的强非线性波进行模拟,并将数值结果与流函数的数值解析解进行了比较,结果表明无论是波面位移、波面处的水平速度和垂向速度均与解析解符合较好,最大波峰面的速度分布伴随水深的增加与解析解吻合程度变差,非线性速度分布的适用范围与线性解析解适应范围kh<3.5基本一致。其次,对深水聚焦波演化进行了模拟研究,研究中聚焦波的生成采用在边界点累加不同频率线性规则波的方法。应用聚焦波物理模型实验结果验证模型,计算聚焦位置处的波面位移和沿水深的速度分布与实验结果的对比表明,波面位移吻合程度较好,垂向的水平速度分布基本吻合。最后,保持中心频率(周期)不变,数值模拟了周期范围变化下最大聚焦波峰面以及波峰面水平速度的变化趋势,结果表明波峰面值和波峰面水平速度随着周期范围缩小而增大。     

An extension of the Airy theory for linear waves into shallow water

As a fully developed (Airy) wave propagates from deep into shallow water, its crest becomes more peaked while the trough flattens out. The median crest diameter MCD, defined as the distance between the wave flanks under the crest at a level halfway between the crest and trough, therefore decreases relative to the similarly defined median trough diameter MTD, which remains constant up to the breaking point. The MCD is directly related to other wave characteristics, which enables water particle velocities to be calculated for any water depth without having to recur to more complex, higher-order Stokes, cnoidal or Fenton theories. Over a nearly horizontal bottom, most fully developed wave characteristics can be expressed as functions of the wave period T_w. It is shown that the horizontal particle velocity at the bottom under the breaker crest is at least 9 times faster than under the breaker trough, which explains why sediment is transported landward under fair weather conditions. The proposed equations also shed new light on the formation of spilling, plunging and surging/collapsing breakers.     

Numerical Prediction of Regular Wave Breaking on Very Gentle Slopes

Regular wave deformation and breaking on very gende slopes is calculated by Mixed-Eulerian-Lagrangian procedure. The velocity potentials and their normal derivatives on the boundary are calculated through the mixed 0-1 boundary element method. The wave elevation and the potentials of time-stepping integration are detertnined by the 2nd-order Taylor expansion at the nodes of free surface boundary elements. During calculation the x-coordinates of the free surface element nodes are supposed to remain unchanged, i.e. the partial derivatives of wave elevation and potentials with respect to x are considered as zero. The numerical results of asymmetric parameters of breaking waves are verified by experimental study. It is shown that when the wave asymmetry is weak, the maximum horizontal velocity of water particales occurs at the wave peak and, the average ratio of this maximum velocity to wave celerity is 0.96. However, when the wave asymmetry is strong, the maximum horizontal velocity of water particles occu     

Effects of following and opposing vertical current shear on nonlinear wave interactions

A Navier-Stokes solver in OpenFOAM^® is combined with the Volume of Fluid (VOF) surface capturing method to investigate the wave interaction with depth-varying currents in intermediate and shallow waters. A special attention is paid to the separate effect of vertical current shear on near resonant triad wave interactions. It was found that in the presence of following vertical current shear, the wave exhibits a sharper crest and flatter trough, and the opposite is true in the presence of opposing vertical current shear. Our model results indicate that the wave steepness at which the current shear starts to affect the crest elevation is greater in deeper water than in shallower water. We found that adding vertical current shear to the uniform current further enhances the relative harmonic wave energy and the extent of triad interaction in the following current while weakens them in the opposing current. As a result, following and opposing current shear may cause wave to break at a lower and higher sea state respectively. Due to the increased wave nonlinearity in the presence of a following current shear, a linear superposition of the individual wave and current velocities is no longer adequate to represent the total horizontal velocity close to the free surface.     

Evolution of wave shape over a low-crested structure

This paper investigates the evolution of wave shape over a low-crested structure (LCS) using a 2-D RANS-VOF model. The model predictions of surface elevation and wave skewness and asymmetry are in good agreement with the recent measurements collected in a small scale wave channel at the University of Cantabria (UCA). The empirical formulae relating wave skewness and asymmetry to local Ursell number by Peng et al. (2009) have been extended to include the effect of wave reflection and the ramp in front of LCS and a wider range of Ursell number in the present study. In the presence of LCS, wave skewness decreases slightly above the seaward slope, then increases rapidly up to a maximum value above the structure crest, and decreases drastically above the leeward slope. Wave asymmetry decreases sharply above the seaward slope to a negative minimum value at the structure crest, and then increases rapidly to a positive value above the leeward slope. Our bispectral analysis indicates that sum interactions increase skewness and decrease asymmetry while difference interactions have opposite effects and that the former dominate above the seaward slope and on the structure crest but the latter dominate above the leeward slope of LCS. The observed wave shape evolution over a LCS can be attributed to the changes in the interplay of sum and difference interactions. We found that incident wave height and wave period, relative structure freeboard, structure crest width and structure porosity are the controlling factors for wave shape evolution over LCS. This study provides new insights on the role of wave skewness and asymmetry in the breakwaters stability and sediment transport around the structure and on the beaches behind it.     

Large-scale turbulence under a solitary wave: Part 2: Forms and evolution of coherent structures

The instantaneous turbulent velocity field produced by a broken solitary wave propagating on a 1 in 50 plane slope was measured in the longitudinal transverse plane in the middle part of the water column and near the bottom using a stereoscopic particle image velocimetry system. These measurements showed that large-scale turbulence first arrived in the form of a downburst of turbulent fluid. In the middle of the water column, the downbursts arrived shortly after the wave crest had passed. Each downburst was accompanied by two counter-rotating vortices. The latter grew rapidly in size to become a prominent feature of the flow field. Each vortex had a typical length scale of 1/2 to 1 water depth, and carried most of the turbulent kinetic energy in the region between the vortices. Near the bottom, the counter-rotating vortices were not as well defined and covered only a small plane area compared to the entire flow structure. The turbulent fluid descending from above diverged at the bed and the resulting flow structure developed an elongated shape as the source of down-flow travelled onshore with the broken wave. It was found that the transverse spacing between adjacent downbursts ranged from 2 to 5 times the local still water depth. Since vortices cannot end in the interior of the fluid, the counter-rotating vortices must extend to the free surface in the form of a vortex loop. It was suggested that these vortex loops were produced by bending and stretching of primary vorticity generated in the wave breaking process, possibly as a result of three-dimensional water surface deformation. The vortex loops were then carried downward by the falling water from the broken wave.     

An experimental investigation of the velocity field under cnoidal waves over the asymmetric rippled bed

1 IntroductionIn coastal areas a ubiquitous phenomenon is theformation of ripples in the seabed. It is now widelyaccepted that the flow and sediment transport overseabed are vital in relation to erosion, surface wavedissipation and pollution dispersion et…     

底部透空不透水建筑物波压力研究

底部透空不透水建筑物是近海工程中的一种新型结构,与传统水工建筑物相比,其能灵活适应基床不发生变形、阻止波浪进入后方造成破坏。本文通过构建物理模型,在波浪水槽中进行试验,测试了底部透空不透水直立建筑物不规则波作用,得到了不同水深、相同波浪要素时,该型直立建筑物迎浪面和底面的水平总力最大时各测点正向波压力和负向波压力数据。对比分析了不同水深下底部透空不透水直立建筑物的水平总力及垂直力,总结了波峰和波谷作用下波浪力沿建筑物表面的分布规律,可为海岸工程结构设计提供依据。     

Numerical predictions of water–air wave slam using incompressible–compressible smoothed particle hydrodynamics

The high-speed impact between a body and water is an important practical problem, whether due to wave impact on a structural deck or wall, or due to a moving body such as a ship or aircraft hitting water. The very high pressures exerted are difficult to predict and the role of air may be significant. In this paper, numerical simulations are undertaken to investigate the impact of a rigid horizontal plate onto a wave crest and, in the limit, onto a flat water surface. A two-phase incompressible–compressible smoothed particle hydrodynamics (SPH) method for water and air, respectively, is applied where the water phase imposes kinematics on the air phase at the air–water interface and the air phase imposes pressures on the water at the interface. Results are compared with experimental measurements undertaken using a drop rig positioned over a wave flume so that a horizontal plate impacts the water surface in free flight. Numerical predictions of impact pressure are quite accurate; air is shown to have a significant cushioning effect for impact on to flat water and this reduces for waves as the ratio of wave height to wavelength increases.     

Modelling of the spatial evolution of extreme laboratory wave crest and trough heights with the NLS-type equations

The statistical properties of long-crested nonlinear wave time series measured in an offshore basin have been analyzed in different aspects such as the distributions of surface elevation, wave crest, wave trough, and wave period. Comparison with linear, second-order and third-order theoretical models indicates that although bound wave effects also contribute to the deviation from a Gaussian process, it is the modulational instability that primarily determines the discrepancy in the evolution process in the presence of strong nonlinearity. Interestingly enough, wave crest is more sensitive to the quasi-resonant four-wave interaction effect than wave trough and the scaled maximal wave crest presents a linear regression model with the coefficient of kurtosis. Meanwhile, the estimation of the observed statistical properties is reconstructed on the basis of an ensemble of 100 wave series simulated by the NLS-type equations and compared favourably with the experimental results in most cases. Moreover, with the increased third-order nonlinear effect the difference between NLS and Dysthe simulations is enlarged and mainly reflected on the distribution of wave crest.     

Breaking of Solitary Waves on Uniform Slopes

The propagation,shoaling and breaking of solitary waves on mild slopes are simulated byboundary element method.In this paper,the criterion of breaking solitary waves on mild slopes is discussed.The criterion is that the ratio of horizontal velocity of water particles on the wave crest to wave celerity equalsone.However,the case that the ratio of horizontal velocity of water particles on the wave crest to wave ce-lerity is below one but the front face of wave profile becomes vertical is also considered as a breaking criteri-on.According to the above criteria,the breaking index for slopes 1:10 to 1:25 is studied.The result is com-pared to other researchers'.The deformation of solitary waves on slopes is discussed and the distribution offluid velocities at breaking is shown.     

Comparative study on the accuracy of solitary wave generations in an ISPH-based numerical wave flume

In the present study, six solitary wave generations by different mathematical approximations are investigated using a piston type wave maker at dimensionless amplitudes ranging from 0.1 to 0.6 and two water depths. Incompressible smoothed particle hydrodynamics is used to simulate solitary wave propagation along the fixed depth channel. The present numerical results are compared with analytical results and experimental data in terms of free surface displacements, fluid particle velocity, phase speed, paddle motion, etc. The present mesh-free numerical results of wave profile variations over time proved that “Rayleigh” has the lowest relative wave height variation. However, its solitary wave has notable phase lead, while “Third order” and “Ninth order” have the least wave lags. Furthermore, the record of present numerical free surface elevation at different distances and the loss of amplitude of the main pulse showed that regarding both of them, “Ninth order” has supremacy over five others. Considering the numerical velocity components of generated solitary wave, “Third order” and “Ninth order” trace analytical results more accurately than other four ones, whereas “Rayleigh” is the most accurate one in predicting the maximum runup. Finally, the paddle motion, its velocity, and displacement, as well as phase speed and outskirts decay coefficient are also compared and discussed intensely.     

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.     

Wave impacts on structures with rectangular geometries: Part 2 decks,baffles and seawalls with impermeable or porous surfaces

This paper considers wave impacts on baffles, on baffles or decks adjacent to a vertical wall, and on porous seawalls and/or sea beds. For seawalls and vertical baffles, impacts can occur in steep waves, whilst a deck can be struck from below by a rising wave crest either in open sea or in a tank with standing waves (sloshing). A simple analytical model for the pressure impulse, P, due to a wave of idealized geometry and dynamics is developed and applied to the following geometries with impermeable surfaces:

• •horizontal wave impact onto a vertical wall with a deck at the waterline,
• •vertical wave impact under a deck in the same configuration (equivalent to vertical water impact of a horizontal plate),
• •horizontal wave impact onto a surface-piercing vertical baffle in open sea,
• •as for 3. but with the baffle in front of a wall,
• •as for 4. but with a deck extending from the vertical wall to the baffle,
• •bottom-mounted baffle in front of a wall with impact occurring on the wall.

We also consider cases that complement part 1 of this paper to include the effect on impacts on a seawall with a porous sea bed and/or sea wall with/without a berm. Finally we reconsider case 3) above but with a porous baffle.The method uses eigenfunction expansions in each of the rectangular regions that satisfy some of the impermeable or porous surface conditions, and a simplified free-surface condition. Their unknown coefficients are determined from the impact boundary condition, impermeable or porous boundary conditions and by matching the solutions, in any two neighbouring rectangles, along their common boundary. Although the fluid motion is treated rather crudely, the method yields the pressure impulse throughout the entire region. Impulses, I, and moment impulses, M, on all or parts of the structure are also presented.     

A Numerical Investigation of the Reduction of Solitary Wave Runup by A Row of Vertical Slotted Piles

To improve the current understanding of the reduction of tsunami-like solitary wave runup by the pile breakwater on a sloping beach, we developed a 3D numerical wave tank based on the CFD tool OpenFOAM in this study. The Navier Stokes equations were applied to solve the two-phase incompressible flow, combined with an LES model to solve the turbulence and a VOF method to capture the free surface. The adopted model was firstly validated with existing empirical formulas for solitary wave runup on the slope without the pile structure. It is then validated using our new laboratory observations of the free surface elevation, the velocity and the pressure around a row of vertical slotted piles subjected to solitary waves, as well as the wave runup on the slope behind the piles. Subsequently, a set of numerical simulations were implemented to analyze the wave reflection, the wave transmission, and the shoreline runup with various offshore wave heights, offshore water depths, adjacent pile spaces and beach slopes. Finally, an improved empirical equation accounting for the maximum wave runup on the slope was proposed by taking the presence of the pile breakwater into consideration.     

Laboratory study of combined wave overtopping and storm surge overflow of a levee

Combined wave overtopping and storm surge overflow of a levee with a trapezoidal cross section was studied in a two-dimensional laboratory wave/flow flume at a nominal prototype-to-model length scale of 25-to-1. The goal of this study was to develop design guidance in the aftermath of Hurricane Katrina. Time series of water depth at two locations on the levee crown and flow thickness at five locations on the landward-side slope were measured along with horizontal velocity near the landward edge of the crown. New equations are presented for average overtopping discharge, distribution of instantaneous discharge, and distribution of individual wave volumes. Equations are also given for mean flow thickness, RMS wave height, mean velocity, and velocity of the wave front down the landward-side slope.     

Study on wave-induced setup over fringing reefs in the presence of a reef crest

It has been well observed that a reef crest (ridge) may be present at the reef edge, but so far very few published studies focusing on the effects of such reef-crest on the wave dynamics over fringing reefs. To understand the role of a reef-crest configuration in determining breaking-wave induced setup over the reef flat, a series of experiments were carried out in a wave flume using an idealized fringing reef model with a reef crest. Experimental results were reported for a trapezoidal reef crest with five reef-crest widths under a series of monochromatic waves. Also examined was the reef without a reef crest. Data analysis shows that larger energy dissipation associated with smaller surfzone width around the reef edge occurred with a wider reef crest. The maximum wave-induced setup on the reef flat in the presence of the reef crest was significantly larger than that without, and it also increased with increasing reef-crest width. The reef-crest submergence was found to be a primary parameter controlling the magnitude of wave setup on the reef flat provided that the reef crest was sufficient wide. An alternative semi-analytical 1DH model based on the balance of cross-shore momentum was proposed. The model was validated by present laboratory data as well as three existing 1DH laboratory studies. Comparing with other two representative semi-analytical models in the literature showed that the proposed model was capable of better reproducing the maximum wave-induced setup on the reef flat for a variety of reef profiles with/without a reef crest, different reef-crest water levels, as well as both monochromatic and spectral waves. The model parameter was physically related to the two characteristic lengths in the surf zone and its value was dependent on the fore-reef slope as well as the presence of a reef crest. The 1DH model was also satisfactorily applied to a fringing reef in field conditions where the effects of fore-reef friction and back-reef lagoon were not important.