Narrow banded wave propagation from very deep waters to the shore

A fully nonlinear Boussinessq-type model with several free coefficients is considered as a departure point. The model is monolayer and low order so as to simplify numerical solvability. The coefficients of the model are here considered functions of the local water depth. In doing so, we allow to improve the dispersive and shoaling properties for narrow banded wave trains in very deep waters. In particular, for monochromatic waves the dispersion and shoaling errors are bounded by ~ 2.8% up to kh = 100, being k the wave number and h the water depth. The proposed model is fully nonlinear in weakly dispersive conditions, so that nonlinear wave decomposition in shallower waters is well reproduced. The model equations are numerically solved using a fourth order scheme and tested against analytical solutions and experimental data.     

Boussinesq–Green–Naghdi rotational water wave theory

Using Boussinesq scaling for water waves while imposing no constraints on rotationality, we derive and test model equations for nonlinear water wave transformation over varying depth. These use polynomial basis functions to create velocity profiles which are inserted into the basic equations of motion keeping terms up to the desired Boussinesq scaling order, and solved in a weighted residual sense. The models show rapid convergence to exact solutions for linear dispersion, shoaling, and orbital velocities; however, properties may be substantially improved for a given order of approximation using asymptotic rearrangements. This improvement is accomplished using the large numbers of degrees of freedom inherent in the definitions of the polynomial basis functions either to match additional terms in a Taylor series, or to minimize errors over a range. Explicit coefficients are given at O(μ²) and O(μ⁴), while more generalized basis functions are given at higher order. Nonlinear performance is somewhat more limited as, for reasons of complexity, we only provide explicitly lower order nonlinear terms. Still, second order harmonics may remain good to kh ≈ 10 for O(μ⁴) equations. Numerical tests for wave transformation over a shoal show good agreement with experiments. Future work will harness the full rotational performance of these systems by incorporating turbulent and viscous stresses into the equations, making them into surf zone models.     

Parameterization of nonlinear shallow water waves over sloping bottoms

Investigation of the bottom slope effects on the nonlinear transformation of irregular waves, which are generated based on JONSWAP spectra, is carried out in a physical wave flume with three slopes (β = 1/15, 1/30, 1/45). The slope effects on the estimation of representative wave height are examined first. To obtain a better estimation of wave height, the slope effect should be considered when slope is larger than 1/30. The nonlinear parameters (bicoherence, skewness and asymmetry) are estimated by using the wavelet-based bispectrum, and the empirical formulae regarding these nonlinear parameters as a function of the local Ursell number are derived based on the present data measured on each slope. The results indicate that the slopes have a negligible effect on the variations of the skewness. The fitted coefficients of the formulae for the other parameters on slope β = 1/15 are clearly different from the results on the slopes β = 1/30 and 1/45, indicating that slope influence on the parameterization cannot be ignored when β > 1/30. Hence, new formulae considering the slope effect are presented. Furthermore, the empirical formulae for the data in surf zone are recommended.     

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.     

A numerical study of nonlinear wave run-up on a vertical plate

A finite difference model based on a recently derived highly-accurate Boussinesq-type formulation is presented. Up to the third-order space derivatives in terms of the velocity variables are retained, and the horizontal velocity variables are re-formulated in terms of a velocity potential. This decreases the total number of unknowns in two horizontal dimensions from seven to five, simplifying the implementation, and leading to increased computational efficiency. Analysis of the embedded properties demonstrates that the resulting model has applications with errors of 2 to 3% for (wavenumber times depth) kh ≤ 10 in terms of dispersion and kh ≤ 4 in terms of internal kinematics. The stability and accuracy of the discrete linearised systems are also analysed for both potential and velocity formulations and the advantages and disadvantages of each are discussed. The velocity potential model is then used to study physically demanding problems involving highly nonlinear wave run-up on a bottom-mounted (surface-piercing) plate. New cases involving oblique incidence are considered. In all cases, comparisons with recent physical experiments demonstrate good quantitative accuracy, even in the most demanding cases, where the local wave steepness can exceed (waveheight divided by wavelength) H / L = 0.20. The velocity potential model is additionally shown to have numerical advantages when dealing with wave–structure interactions, requiring less smoothing around exterior structural corners.     

How much data is enough? The importance of morphological sampling interval and duration for calibration of empirical shoreline models

The ability to robustly predict future shoreline position under the influence of changing waves and sea-level rise is a key challenge to scientists and engineers alike. While extrapolating a linear trend out in time is a common baseline approach, the recent development of a number of empirical shoreline models allows the prediction of storm and annual-scale variability as well. The largest constraint in applying these models is the availability of high quality, adequate duration data sets in order to calibrate model free parameters. This contribution outlines several such models and discusses the monitoring programs required to calibrate and hindcast shoreline change from 1 to 10 years at two distinct beach types: a storm-dominated site and the second exhibiting a large seasonal variability. The seasonally-dominated site required longer data sets but was less sensitive to sampling interval, while the storm-dominated site converged on shorter, more frequently sampled data sets. In general, calibration based on a single year of observed shorelines resulted in a large range of model skill and was not considered robust. Monitoring programs of at least two years, with shorelines sampled at dt ≤ 30 days were sufficient to determine initial estimates of calibration coefficients and hindcast short-term (1–5 years) shoreline variability. In the presence of unresolved model processes and noise, hindcasting longer (5 + years) data sets required longer (5 + years) calibration data sets, particularly when sampling intervals exceeded 60 days.     

Comparison of time-dependent extended mild-slope equations for random waves

The extended mild-slope equations of Suh et al. [Suh, K.D., Lee, C., Park, W.S., 1997. Time-dependent equations for wave propagation on rapidly varying topography. Coastal Eng., 32, 91–117] and Lee et al. [Lee, C., Kim, G., Suh, K.D., 2003. Extended mild-slope equation for random waves. Coastal Eng., 48, 277–287] are compared analytically and numerically to determine their applicability to random wave transformation. The geometric optics approach is used to compare the two models analytically. In the model of Suh et al., the wave number of the component wave with a local angular frequency ω is approximated with an accuracy of O(ω − ω¯) at a constant water depth, where ω¯ is the carrier frequency of random waves. In the model of Suh et al., however, the diffraction effects and higher-order bottom effects are considered only for monochromatic waves, and the shoaling coefficient of random waves is not accurately approximated. This inaccuracy arises because the model of Suh et al. was derived for regular waves. In the model of Lee et al., all the parameters of random waves such as wave number, shoaling coefficient, diffraction effects, and higher-order bottom effects are approximated with an accuracy of O(ω − ω¯). This approximation is because the model of Lee et al. was developed using the Taylor series expansion technique for random waves. The result of dispersion relation analysis suggests the use of the peak and weighted-average frequencies as a carrier frequency for Suh et al. and Lee et al. models, respectively. All the analytical results are verified by numerical experiments of shoaling of random waves over a slightly inclined bed and diffraction of random waves through a breakwater gap on a flat bottom.     

Hydrodynamic performance of double rows of piles suspending horizontal c shaped bars

The wave transmission, reflection, and energy dissipation of the double rows of vertical piles suspending horizontal steel C shaped bars are experimentally and theoretically studied under normal regular waves. Different wave and structural parameters are investigated e.g. the wave length, the C shaped bars draft and spacing, the supporting piles diameter and spacing, and the space between the double rows. Also, the theoretical model based on an eigenfunction expansion method is developed to study the hydrodynamic breakwater performance. In order to examine the validity of the theoretical model, the theoretical results are compared with the experimental and theoretical results obtained by different authors. Comparison between experiments and predictions showed that theoretical model provides a good estimate to the different hydrodynamic coefficients when the friction factors of the upper and the lower parts are f_U = 1.5 and f_L = 0.75. The present breakwater physical model gives efficiency near other similar systems of different shapes.     

Numerical study of turbulence and wave damping induced by vegetation canopies

Vegetation canopies control mean and turbulent flow structure as well as surface wave processes in coastal regions. A non-hydrostatic RANS model based on NHWAVE (Ma et al., 2012) is developed to study turbulent mixing, surface wave attenuation and nearshore circulation induced by vegetation. A nonlinear k − ϵ model accounting for vegetation-induced turbulence production is implemented to study turbulent flow within the vegetation field. The model is calibrated and validated using experimental data from vegetated open channel flow, as well as nonbreaking and breaking random wave propagation in vegetation fields. It is found that the drag-related coefficients in the k − ϵ model C_fk and C_fϵ can greatly affect turbulent flow structure, but seldom change the wave attenuation rate. The bulk drag coefficient C_D is the major parameter controlling surface wave damping by vegetation canopies. Using the empirical formula of Mendez and Losada (2004), the present model provides accurate predictions of vegetation-induced wave energy dissipation. Wave propagation through a finite patch of vegetation in the surf zone is investigated as well. It is found that the presence of a finite patch of vegetation may generate strong pressure-driven nearshore currents, with an onshore mean flow in the unvegetated zone and an offshore return flow in the vegetated zone.     

The influence of sea state on formation speed of alongshore variability in surf zone sand bars

The formation time of alongshore morphological variability in surf zone sand bars has long been known to differ from one beach to the other and from one post-storm period to another. Here we investigate whether the type of sea state, i.e. distant swell waves or locally generated short period wind sea, affects the formation time of the emerging alongshore topographic variability.A numerical modeling approach is used to examine the emergence of alongshore variability under different shore-normal wave forcing. A research version of Delft3D, operating on the time-scale of wave groups, is applied to a schematised bathymetry with a single bar. The model is then used to investigate several wave scenarios, examining the impact of peak period, frequency spread and directional spread on the formation time of alongshore variability.Results show that an increase in wave period has a large effect, changing the formation time up to O (250%) in case the wave period is changed from a representative value for the Dutch coast (T_p ~ 5–6 s) to an Australian South East coast value (T_p ~ 10–12 s). In contrast, modifications in the directional and frequency spread of the wave field result only in a minor change in the formation time.Examination of hydrodynamics and potential sediment transport shows that the variations in formation time are primarily related to changes in the magnitude of the time-averaged flow conditions. Variations in the magnitude of very low frequency (f < 0.004 Hz) or infragravity (0.004 < f < 0.04 Hz) surf zone flow velocities do not affect the mean sediment transport capacity. Consequently the formation speed of patterns is primarily governed by positive feedback between mean flow and morphology, and low frequency flow fluctuations are of minor importance.These findings indicate that the development of alongshore topographic variability may be faster at swell dominated open coasts, primarily due to the occurrence of longer period swell. Also, at a given site, the arrival of a long wave period swell after a storm can accelerate the emergence of variability.     

Surface wave predictions in weakly nonlinear directional seas

We have employed laboratory and numerical experiments in order to investigate propagation of waves in both long and short-crested wave fields in deep water. For long-crested waves with steepness, ϵ = k_ca_c = 0.1 (a fairly extreme case), reliable prediction can be performed with the modified nonlinear Schrödinger equation up to about 40 characteristic wavelengths. For short-crested waves the accuracy of prediction is strongly reduced with increasing directional spread.     

Numerical simulations of wave–current flow in an ocean basin

Wave–current flow is a phenomenon that is present in many practical engineering situations. Over the past several decades, this type of flow has been increasingly investigated under controlled laboratory conditions. This paper presents a numerical study of wave–current flow in the ocean basin of the LabOceano (COPPE/UFRJ). A homogeneous multiphase model based on the RANS equations and the k–ɛ turbulence model implemented in ANSYS-CFX code were used. A cross section of the ocean basin was represented. A regular wave with a height of 0.08 m and a period of 1.80 s (i.e., a wave steepness of H/L = 0.016), propagating on favourable currents, was simulated. The behaviour of the free surface elevation over time and the streamlines along the basin for wave and wave–current flows were presented. The numerical results were compared to the non-viscous theory given by the Rayleigh equation applied to the problem of wave–current interaction. Good agreement was found between the wave length estimated by the numerical results and the analytical solutions, with a deviation of less than 2%.     

Analytical investigation of hydrodynamic performance of a dual pontoon WEC-type breakwater

Based on the linear potential flow theory and matching eigen-function expansion technique, an analytical model is developed to investigate the hydrodynamics of two-dimensional dual-pontoon floating breakwaters that also work as oscillating buoy wave energy converters (referred to as the integrated system hereafter). The pontoons are constrained to heave motion independently and the linear power take-off damping is used to calculate the absorbed power. The proposed model is verified by using the energy conservation principle. The effects of the geometrical parameters on the hydrodynamic properties of the integrated system, including the reflection and transmission coefficients and CWR (capture width ratio, which is defined as the ratio of absorbed wave power to the incident wave power in the device width). It is found that the natural frequency of the heave motion and the spacing of the two pontoons are the critical factors affecting the performance of the integrated system. The comparison between the results of the dual-pontoon breakwater and those of the single-pontoon breakwater shows that the effective frequency range (for condition of transmission coefficient K_T < 0.5 and the total capture width ratio η_total > 20%) of the dual-pontoon system is broader than that of the single-pontoon system with the same total volume. For the two-pontoon system, the effective frequency range can be broadened by decreasing the draft of the front pontoon within certain range.     

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.     

Wave attenuation in mangroves: A quantitative approach to field observations

Coastal mangroves, dwelling at the interface between land and sea, provide an important contribution to reducing risk from coastal hazards by attenuating incident waves and by trapping and stabilizing sediments. This paper focusses on relations between vegetation densities, wave attenuation rates, sediment characteristics and sedimentation rates in mangroves. These processes were studied along two cross-shore transects through mangroves fringing estuaries in the southern Andaman region of Thailand. Volumetric vegetation densities in these mangroves were ranging up to 32‰, depending on the water depth. Generalized total wave attenuation rates increased from 0.002 m^− 1 in the sparsely vegetated forest fringes with Avicennia and Sonneratia species, up to 0.012 m^− 1 in the dense Rhizophora vegetation in the back of the forests. The total wave attenuation rates integrate effects of shoaling and energy losses due to various bio-physical interactions within the mangrove ecosystem. Wave attenuation in the mangroves is presumably dominated by energy losses due to vegetation drag, since wave attenuation due to bottom friction and viscous dissipation on the bare mudflats is significantly lower than those inside the mangrove vegetation.Additionally, wave attenuation in the mangroves was found to facilitate enhanced net sediment deposition and a gradual fining of the bed material. These findings corroborate the coastal defence function of mangroves by quantifying their contribution to wave attenuation and sediment trapping. The explicit linking of these properties to vegetation composition and structure facilitates modelling studies investigating the mechanisms determining the coastal defence capacities of mangroves.     

Wave climate scatter performance of a cycloidal wave energy converter

A lift based cycloidal wave energy converter (WEC) was investigated using potential flow numerical simulations in combination with viscous loss estimates based on published hydrofoil data. This type of wave energy converter consists of a shaft with one or more hydrofoils attached eccentrically at a radius. The main shaft is aligned parallel to the wave crests and submerged at a fixed depth. The operation of the WEC as a wave-to-shaft energy converter interacting with straight crested waves was estimated for an actual ocean wave climate. The climate chosen was the climate recorded by a buoy off the north-east shore of Oahu/Hawaii, which was a typical moderate wave climate featuring an average annual wave power P_W = 17 kWh/m of wave crest. The impact of the design variables radius, chord, span and maximum generator power on the average annual shaft energy yield, capacity factor and power production time fraction were explored. In the selected wave climate, a radius R = 5 m, chord C = 5 m and span of S = 60 m along with a maximum generator power of P_G = 1.25 MW were found to be optimal in terms of annual shaft energy yield. At the design point, the CycWEC achieved a wave-to-shaft power efficiency of 70%. In the annual average, 40% of the incoming wave energy was converted to shaft energy, and a capacity factor of 42% was achieved. These numbers exceeded the typical performance of competing renewables like wind power, and demonstrated that the WEC was able to convert wave energy to shaft energy efficiently for a range of wave periods and wave heights as encountered in a typical wave climate.     

Geometry and migration characteristics of bedforms under waves and currents: Part 2: Ripples superimposed on sandwaves

This paper presents results of a series of detailed measurements of geometric and migrating characteristics of ripples superimposed upon sandwaves under the action of combined waves and currents. Velocity measurements within the fluid, surface wave characteristics and 3D mapping of the bottom were recorded with an Acoustic Doppler Velocimeter (ADV), an acoustic water level sensor and a 32 composite element array of sub-aquatic acoustic sensors, respectively. Bottom records were statistically analyzed to obtain height, length and migration rates of ripples. Experiments examined ripple heights and wavelengths for the mobility factors (as defined in Eq. (4)) and the Reynolds wave number within the ranges 10 < ψ < 88 and 16 × 10³ < Re_w < 5 × 10⁵, respectively. Measured values were compared with laboratory and field data together with semi-empirical and analytical formulae from the literature. Good correlation was obtained when plotting measured ripple length and length in dimensionless form as a function of the Reynolds wave number Re_w. Under a given hydraulic condition, it was observed that ripples with different geometric characteristics may coexist at different locations over the sandwave. Ripple steepness is presented as a function of the Shields parameter although characterized with rather large scatter. Finally, average ripple migration speed is presented as a function of the Shields parameter and the mobility number.     

Modeling wave-mud interaction on the central chenier-plain coast,western Louisiana Shelf,USA

The strong coupling between hydrodynamics and seafloors on shallow muddy shelves, and resulting bed reworking, have been extensively documented. On these shelves, spectral wave transformation is driven by a complex combination of forcing mechanisms that include nonlinear wave interactions and wave energy dissipation induced by fluid-mud at a range of frequencies. Wave-mud interaction is investigated herein by using a previously validated nonlinear spectral wave model and observations of waves and near-bed conditions on a mildly-sloping seafloor off the muddy central chenier-plain coast, western Louisiana Shelf, United States. Measurements were made along a cross-shelf transect spanning 1 km between 4 and 3 m water depths. The high-resolution observations of waves and near-bed conditions suggest presence of a fluid mud layer with thickness sometimes exceeding 10 cm under strong long wave action (1 meter wave height with 7 s peak period at 4 meter depth). Spectral wave transformation is modeled using the stochastic formulation of the nonlinear Mild Slope Equation, modified to account for wave-breaking and mud-induced dissipation. The model is used in an inverse manner in order to estimate the viscosity of the fluid mud layer, which is a key parameter controlling mud-induced wave dissipation but complicated to measure in the field during major wave events. Estimated kinematic viscosities vary between 10⁻⁴-10⁻³ m²/s. Combining these results of the wave model simulations with in-depth analysis of near-bed conditions and boundary layer modeling allows for a detailed investigation of the interaction of nonlinear wave propagation and mud characteristics. The results indicate that mud-induced dissipation is most efficient when the wave-induced resuspensions of concentrations  > 10 g/L settle due to relatively small bottom stresses to form a fluid mud layer that is not as thin and viscous as a consolidated seafloor in absence of wave action but also not as thick and soft as a near-bed high concentration layer that forms during strong wave action.