WAVECALC: an Excel-VBA spreadsheet to model the characteristics of fully developed waves and their influence on bottom sediments in different water depths

The generation and growth of waves in deep water is controlled by winds blowing over the sea surface. In fully developed sea states, where winds and waves are in equilibrium, wave parameters may be calculated directly from the wind velocity. We provide an Excel spreadsheet to compute the wave period, length, height and celerity, as well as horizontal and vertical particle velocities for any water depth, bottom slope, and distance below the reference water level. The wave profile and propagation can also be visualized for any water depth, modeling the sea surface change from sinusoidal to trochoidal and finally cnoidal profiles into shallow water. Bedload entrainment is estimated under both the wave crest and the trough, using the horizontal water particle velocity at the top of the boundary layer. The calculations are programmed in an Excel file called WAVECALC, which is available online to authorized users. Although many of the recently published formulas are based on theoretical arguments, the values agree well with several existing theories and limited field and laboratory observations. WAVECALC is a user-friendly program intended for sedimentologists, coastal engineers and oceanographers, as well as marine ecologists and biologists. It provides a rapid means to calculate many wave characteristics required in coastal and shallow marine studies, and can also serve as an educational tool.     

A function to determine wavelength from deep into shallow water based on the length of the cnoidal wave at breaking—Reply to discussion by T.S. Hedges

The equations of Hedges [Hedges, T.S., 2009. Discussion of “A function to determine wavelength from deep into shallow water based on the length of the cnoidal wave at breaking” by J.P. Le Roux, Coastal Eng.], although yielding similar wavelengths, are not consistent with the fact that the horizontal water particle velocity in the wave crest should equal the wave celerity at breaking over a nearly horizontal bottom.     

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.     

Accuracy of asymptotic series for the kinematics of high solitary waves

The accuracy of several asymptotic series expansions for wave speed and particle velocity under the crest of a solitary wave (on a fluid at rest) up to maximum height is investigated. The very accurate numerical results of Williams (1985) are the measure for our comparisons. The results are based on a scaling of calculated properties of long periodic waves to the case of solitary waves.For wave speeds the classical Boussinesq–Rayleigh expression gives good agreement up to a relative wave height of, say, 0.3. An asymptotic fourth-order expression based on Fenton (1990) can be used up to a relative wave height of 0.7, whereas the corresponding fifth-order expression is slightly less accurate.The Eulerian particle velocity profile under the wave crest is examined using a cnoidal wave expression from Fenton (1990) in the limit of the solitary wave. For low waves a `consistent' (i.e. properly truncated) fifth-order expression and an `inconsistent' ditto both coincide with Williams' results. Beginning at medium high waves, the consistent expression surprisingly exhibits oscillations in the velocity profile, and the oscillations become stronger as the wave gets higher. The inconsistent expression, however, yields the same shape as Williams' profile, but is displaced parallel to this, resulting in slightly larger velocities. For high waves also the inconsistent expression begins to differ in shape from Williams' profile, and asymptotic theory fails. Only for low waves `lowest order theory' gives acceptable results. We show analytically that for the highest wave the particle velocity profile has a horizontal tangent at the water surface; this is corroborated by Williams' numerical results.We also study the particle velocity at the wave crest as a function of wave height. It is shown that the variation has a vertical tangent for the highest wave. Two fifth-order asymptotic series for this velocity, based on the wave speed through the Bernoulli equation, show very good agreement with Williams up to a relative wave height of about 0.6.It is finally shown that it is possible to produce very accurate rational-function approximations to Williams' results for the wave speed as well as for the particle velocity at the wave crest.     

Kinematics of extreme waves in deep water

The velocity profiles under crest of a total of 62 different steep wave events in deep water are measured in laboratory using particle image velocimetry. The waves take place in the leading unsteady part of a wave train, focusing wave fields and random wave series. Complementary fully nonlinear theoretical/numerical wave computations are performed. The experimental velocities have been put on a nondimensional form in the following way: from the wave record (at a fixed point) the (local) trough-to-trough period, T_TT and the maximal elevation above mean water level, η_m of an individual steep wave event are identified. The local wavenumber, k and an estimate of the wave slope, ε are evaluated from ω²/(gk)=1+ε², where ω=2π/T_TT and g denotes the acceleration of gravity. A reference fluid velocity, is then defined. Deep water waves with a fluid velocity up to 75% of the estimated wave speed are measured. The corresponding kη_m is 0.62. A strong collapse of the nondimensional experimental velocity profiles is found. This is also true with the fully nonlinear computations of transient waves. There is excellent agreement between the present measurements and previously published Laser Doppler Anemometry data. A surprising result, obtained by comparison, is that the nondimensional experimental velocities fit with the exponential profile, i.e. e^ky, y the vertical coordinate, with y=0 in the mean water level.     

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.     

Wave force on vertically submerged circular thin plate in shallow water

The paper presented is a solution of shallow water wave force, using small amplitude linear wave theory on two-dimensional vertically submerged circular thin plates under three different configurations: (1) a surface-piercing circular thin plate, (2) a submerged circular thin plate, and (3) a bottom-standing circular thin plate. Finally Morison's equation is used for the determination of wave force which is based on the linear wave theory. The plate is submerged in water near the shore on uniformly sloping bottom. The solution method is confined in a finite domain, which contains both the region of different depth of water and the plate. Laplace's equation and boundary value problems are solved in a finite domain, by the method of separation of variables and the small amplitude linear wave theory. The variation of horizontal force by single particle, total horizontal force and moment with respect to the wave amplitude are obtained at different depth of water and at different wave period. It is observed that the force and moment are converging with the increase of wave period and the gradients of force and moment with respect to the wave amplitude are extremely high for lower wave period.     

High order Hamiltonian water wave models with wave-breaking mechanism

Based on the Hamiltonian formulation of water waves, using Hamiltonian consistent modelling methods, we derive higher order Hamiltonian equations by Taylor expansions of the potential and the vertical velocity around the still water level. The polynomial expansion in wave height is mixed with pseudo-differential operators that preserve the exact dispersion relation. The consistent approximate equations have inherited the Hamiltonian structure and give exact conservation of the approximate energy. In order to deal with breaking waves, we extend the eddy-viscosity model of Kennedy et al. (2000) to be applicable for fully dispersive equations. As breaking trigger mechanism we use a kinematic criterion based on the quotient of horizontal fluid velocity at the crest and the crest speed. The performance is illustrated by comparing simulations with experimental data for an irregular breaking wave with a peak period of 12 s above deep water and for a bathymetry induced periodic wave plunging breaker over a trapezoidal bar. The comparisons show that the higher order models perform quite well; the extension with the breaking wave mechanism improves the simulations significantly.     

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.     

The continuity equation for longshore current velocity with breaker angle adjusted for a wave-current interaction

The continuity equation for mean longshore current velocity, V = gmT sin 2θ_b, agrees with selected field and laboratory data covering a wide range of conditions. Agreement between continuity equation and data is improved by eliminating those laboratory data which imply deep-water wave crests at angles near or greater than 90 degrees to the shoreline. Agreement between continuity equation and data is further improved by adjusting breaker angles to account for convection of the breaker point by the longshore current. Breaker point convection increases breaker angle by an amount predictable from the analysis developed here. This increase in angle is significant in those laboratory experiments with breaking wave crests at high angles to the shoreline.In the continuity equation, m is bottom slope, T is wave period, and θ_b is breaker angle, but breaker height does not appear. According to radiation stress theory, mean velocity does depend on breaker height, but only weakly. Consistency between the two approaches would require a dimensionless velocity, C_b/gT, to be relatively constant, which it is. (The same dimensionless velocity appears in the analyses of breaker point convection.) The continuity equation is functionally independent of friction and mixing, in keeping with its derivation from simple conservation of mass considerations. The equation has no adjustable coefficients. The degress of agreement with data and the internal consistency of the analysis suggests that it is a good predictor of mean velocity in ordinary longshore currents.     

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.     

One-dimensional modelling of individual waves and wave-induced longshore currents in the surf zone

A probabilistic model ( -model) was developed to describe the propagation and transformation of individual waves (wave by wave approach). The individual waves shoal until an empirical criterion for breaking is satisfied. Wave height decay after breaking is modelled by using an energy dissipation method. Wave-induced set-up and set-down and breaking-associated longshore currents are also modelled. Laboratory and field data were used to calibrate and verify the model. The model was calibrated by adjusting the wave breaking coefficient (as a function of local wave steepness and bottom slope) to obtain optimum agreement between measured and computed wave height. Four tests carried out in the large Delta flume of Delft Hydraulics were considered. Generally, the measured H_1/3-wave heights are reasonably well represented by the model in all zones from deep water to the shallow surf zone. The fraction of breaking waves was reasonably well represented by the model in the upsloping zones of the bottom profile. Verification of the model results with respect to wave-induced longshore current velocities was not extensive, because of a lack of data. In case of a barred profile the measured longshore velocities showed a relatively uniform distribution in the (trough) zone between the bar crest and the shoreline, which could to some extent be modelled by including space-averaging of the radiation force gradient, horizontal mixing and longshore water surface gradients related to variations in set-up. In case of a monotonically upsloping profile the cross-shore distribution of the longshore current velocities is reasonably well represented.     

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.     

Transformation of a strongly nonlinear wave in a shallow-water basin

The transformation of a nonlinear wave in shallow water is investigated analytically and numerically within the framework of long-wave theory. It is shown that the nonlinearity parameter (the Mach number), which is defined as the ratio of the particle velocity in the wave to the propagation velocity, can be well above unity in a deep trough and that a jump appears initially in the trough. It is demonstrated that shockwave amplitudes at large times change in accordance with the prediction of weakly nonlinear theory. The shock front generates a reflected wave, which, in turn, transforms into a shock wave if the initial amplitude is large enough. The amplitude of the reflected wave is proportional to the cube of the initial amplitude (as predicted by weakly nonlinear theory) over a wide range of amplitudes except for the case of anomalously strong nonlinearity. When there is a sign-variable sufficiently intense initial perturbation, the basic wave transforms into a positive shock pulse (crest) and the reflected wave turns into a negative pulse (trough).     

斜向波浪对半圆型防波堤作用的实验研究

通过正向与斜向波浪对半圆型防波堤(不开孔出水堤)的实验研究,给出了其水力特性及单位堤长所受的无因次水平波浪力、竖向波浪力和波浪浮托力随各主要影响因素的变化规律,特别指出了斜向波浪力可大于正向波浪力及水平波浪力中波谷作用大于波峰作用发生的条件。     

An experimental investigation of plunging breaker boundary layers in the transformation zone

Experimental investigation is made on the boundary layers of the transformation zone (i.e. the region between the last symmetrical wave profile depth and the breaking point) of plunging breakers propagating on a smooth beach with 1/12 uniform slope. Using a laser anemometer, the particle velocities are measured at four verticals along the transformation zone for three different steepnesses of waves within the plunging breaker range. The boundary layer flow in the transformation zone is found mostly of turbulent character and vertical distribution of particle velocities does not seem to conform to the classical law of the wall distribution given for steady-flow boundary layers. The results show that free-stream particle velocities, in the boundary layer of the breaker under the crest phase, increase considerably as the wave progresses towards the breaking point. The boundary layer thickness, defined as the velocity-affected region, remains constant throughout the transformation zone but it decreases with increasing deep-water wave steepness for the particular beach slope tested.     

Wave crest and trough distributions in a broad-banded directional wave field

It is well established that the modulational instability enhances the probability of occurrence for extreme events in long crested wave fields. Recent studies, however, have shown that the coexistence of directional wave components can reduce the effects related to the modulational instability. Here, numerical simulations of the Euler equations are used to investigate whether the modulational instability may produce significant deviations from second-order statistical properties of surface gravity waves when short crestness (i.e., directionality) is accounted for. The case of a broad-banded directional wave field (i.e. wind sea) is investigated. The analysis is concentrated on the wave crest and trough distribution. For completeness a comparison with a unidirectional wave field is presented also. Results will show that the distributions based on second-order theory provide a good estimate for the simulated crest and trough height also at low probability levels.     

直立堤浮托力的概率特征

本文基于1987和1988年夏季在古镇口港同步连续观测的波浪、波压力和浮托力资料,分析研究了浮托力的概率特征、浮托力沿堤底面的变化及其谱特征。结果表明,浮托力幅度和周期的累积分布实际上不随测点在堤底面上的位置变化,可用公式表示。浮托力沿堤底面宽度近似呈线性衰减,前趾浮托力最大,后趾不为零,呈梯形状。文中还给出了计算最大峰、谷总浮托力的方法。     

Wave forces and moments on an intake well

A sea water intake well of size 20 m diameter and 15.5 m height in a water depth of 9.8 m is proposed north of the Visakhapatnam Port for a project to extract magnesia from sea water. A 1:25 scale model of the intake well was tested in the wave basin of the Ocean Engineering Centre, Indian Institute of Technology, Madras to measure the wave forces and moments on the intake well and the variation of water levels inside and outside the well. Accordingly, an intake well model of 0.8 m diameter and 0.62 m height was fabricated and fixed over a false bottom in a wave basin. The well model was subjected to the action of both regular waves for two test conditions, intake well inlet closed during installation and intake well inlet open. The experimental results on wave forces and moments were compared with the results of the Linear Diffraction Theory. The water level inside the well was measured to determine the submergence of suction pipes of pumps and location of the inlet opening of the intake well. The wave crest elevation in front of the well was also measured in order to fix the deck level of the well so as to avoid water overspill onto the deck. The salient results of the present study are presented and discussed in this paper.