Application of a fuzzy inference system for the prediction of longshore sediment transport

A fuzzy inference system (FIS) and a hybrid adaptive network-based fuzzy inference system (ANFIS), which combines a fuzzy inference system and a neural network, are used to predict and model longshore sediment transport (LST). The measurement data (field and experimental data) obtained from Kamphuis [1] and Smith et al. [2] were used to develop the model. The FIS and ANFIS models employ five inputs (breaking wave height, breaking wave angle, slope at the breaking point, peak wave period and median grain size) and one output (longshore sediment transport rate). The criteria used to measure the performances of the models include the bias, the root mean square error, the scatter index and the coefficients of determination and correlation. The results indicate that the ANFIS model is superior to the FIS model for predicting LST rates. To verify the ANFIS model, the model was applied to the Karaburun coastal region, which is located along the southwestern coast of the Black Sea. The LST rates obtained from the ANFIS model were compared with the field measurements, the CERC [3] formula, the Kamphuis [1] formula and the numerical model (LITPACK). The percentages of error between the measured rates and the calculated LST rates based on the ANFIS method, the CERC formula (K_sig = 0.39), the calibrated CERC formula (K_sig = 0.08), the Kamphuis [1] formula and the numerical model (LITPACK) are 6.5%, 413.9%, 6.9%, 15.3% and 18.1%, respectively. The comparison of the results suggests that the ANFIS model is superior to the FIS model for predicting LST rates and performs significantly better than the tested empirical formulas and the numerical model.     

Coefficient of horizontal eddy diffusion of a tracer on the water surface as a function of the wave age

A relationship between the statistical parameters of horizontal diffusion and the parameters of the energy-containing part of the frequency spectrum of sea-surface elevations is found depending on the wave age Ω and the ratio between the wind speed at 10 m and the phase velocity of the peak of a wave. It has been observed in [1–7] that the diffusion coefficient K(r) of a patch of size r increases as r ^β, where 1.15 < β < 4/3, and the patch area S(t) increases with time as t ^γ, where 2 < γ < 3. As was calculated in [15], in the energy-containing part of the elevation frequency spectrum, S(ω) ～ ω^?n , where n = 13/3 for young waves with Ω > 2, n = 4 for waves with 1.2 < ψ < 2, and n = 11/3 for developed waves with 0.83 < Ω < 1.2. It is found that β = (n + 1)/4 and γ = 8/(7 ? n). These relations explain the entire set of observed exponents: β = 4/3 and γ = 3 for young waves and β = 1.15 and γ = 2.34 for large sizes (up to 1000 km) and times (up to a month) when it is found here that β = 7/6 and γ = 2.4.     

Verification of a parabolic equation for propagation of weakly-nonlinear waves

Verification tests of linear wave propagation models for cases of waves propagating through a refractive focus typically overpredict peak amplitudes in the vicinity of the focus. Using a parabolic equation model for the combined refraction-diffraction of weakly-nonlinear waves [valid for Ursell number U_r < O(1)], we show that, in a particular experiment, much of the discrepancy between laboratory measurements and linear wave model predictions can be accounted for by including the effect of nonlinearity, which is important in the region of the focus.     

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.     

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

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

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.     

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.     

Phase-lag effects in sheet flow transport

The inception of the sheet flow regime as well as the effects of the phase lag when the sheet flow regime is established were investigated for oscillatory flows and combined steady and oscillatory flows. A new criterion for the inception of sheet flow is proposed based on around 300 oscillatory flow cases from experiments. This criterion was introduced in the Camenen and Larson [Camenen, B., Larson, M., 2005. A bedload sediment transport formula for the nearshore. Estuarine, Coastal and Shelf Science 63, 249–260.] bed load formula in order to take into account phase-lag effects in the sheet flow regime. The modification of the Camenen and Larson formula significantly improves the overall agreement with data and yields a correct behavior in relation to some of the main governing parameters, which are the median grain size d₅₀, the orbital wave velocity U_w, and the wave period T_w. The calibration of the new formula was based on more than 200 experimental data values on the net sediment transport rate for a full wave cycle. A conceptual model was also proposed to estimate the ratio between sediment transport rate with and without phase lag, (r_pl = q_s,net / q_s,net,ϕ=0). This simple model provides accurate results and may be used together with any quasi-steady model for bed load transport.     

Characteristic friction coefficient and scale effects in oscillatory porous flow

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

Study on far field wave patterns and their characteristics of Havelock form green function

A new mathematical integral representation including five integrals about the far field wave shape function of Havelock form translating-pulsating source is obtained by performing variable substitution. Constant-phase curves and propagation wave patterns are investigated by applying stationary phase analysis method to the new representation. Some findings are summarized as follows: (1) when 0< <0.25 (where is the Strouhal number), three types of stationary phase curves corresponding to three propagation wave patterns such as fan wave pattern, inner V and outer V wave patterns, are found in the integral representation. (2) When >0.25, besides three types of wave patterns such as a ring-faning wave pattern, a fan wave pattern and an inner V wave pattern, a new one called parallel wave pattern is also found which not only exists in the integrals about the ring-fan wave and fan wave, but also in the integrals whose interval is 0, . In addition, Characteristics about these parallel waves such as mathematical expressions, existence conditions, propagation directions and wave lengths are obtained, and cancellation relationships between these parallel waves are stated, which certificates the fact that there are no parallel waves existing in the far field.     

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.     

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.     

A coupled numerical model of wave interaction with porous medium

A new coupling model of wave interaction with porous medium is established in which the wave field solver is based on the two dimensional Reynolds Averaged Navier-Stokes (RANS) equations with a kε closure. Incident waves, which could be linear waves, cnoidal waves or solitary waves, are produced by a piston-type wave maker in the computational domain and the free surface is traced through the Piecewise Linear Interface Construction-Volume of Fluid (PLIC-VOF) method. Nonlinear Forchheimer equations are adopted to calculate the flow field within the porous media. By introducing a velocity–pressure correction equation, the wave field and the porous flow field are highly and efficiently coupled. The two fields are solved simultaneously and no boundary condition is needed at the interface of the internal porous flow and the external wave. The newly developed numerical model is used to simulate wave interaction with porous seabed and the numerical results agree well with the experimental data. The additional numerical tests are also conducted to study the effects of seabed thickness, porosity and permeability coefficient on wave damping and the pore water pressure responses.     

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.     

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

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

Formation of beach cusps in a wave tank

Beach cusps with a longshore spacing of 20 to 150 cm have been built by the continuous action of incident waves on a steep laboratory beach floor covered uniformly with a thin bed of glass beads. Breaking of incident waves was observed to induce vortices on the bed by interacting with swash motion along the beach face. Beach cusps formed when the value of a dimensionless parameter H_b/sgT_i² became smaller than 0.042; H_b is the breaking height of the incident waves, T_i their period, s the beach slope and g the acceleration due to gravity. This critical value occurred at a nearly central part of the generation region 0.003 < H_b/sgT_i² < 0.068 for plunging breakers presented by Galvin (1968). Breaking-wave-induced vortices rather than breaker types controlled the movement of bed material in the nearshore zone. Most of the measured spacings of beach cusps, including previous observations, were in good agreement with half a wavelength of the zero-mode subharmonic edge wave, which is generated on the beach by the refraction of incident waves and has twice the period of the waves. The role of edge waves at each stage of cusp formation still remains as an important problem to be clarified.     

Wave height parameter for damage description of rubble-mound breakwaters

In this paper it will be shown that the wave height parameter H₅₀, defined as the average wave height of the 50 highest waves reaching a rubble-mound breakwater in its useful life, can describe the effect of the wave height on the history of the armor damage caused by the wave climate during the structure's usable life.Using Thompson and Shuttler (Thompson, D.M., Shuttler, R.M., 1975. Riprap design for wind wave attack: A laboratory study on random waves. HRS Wallingford, Report 61, UK) data it will be shown that H₅₀ is the wave parameter that best represents the damage evolution with the number of waves in a sea state. Using this H₅₀ parameter, formulae as van der Meer (van der Meer, J.W., 1988. Rock slopes and gravel beaches under wave attack. PhD Thesis. Technical University of Delft) and Losada and Giménez-Curto (Losada, M.A., Gimenez–Curto, L.A., 1979. The joint effect of the wave height and period on the stability of rubble mound breakwaters using Iribarren's number. Coastal Engineering, 3, 77–96) are transformed into sea-state damage evolution formulae. Using these H₅₀-transformed formulae for regular and irregular sea states it will be shown how damage predictions are independent of the sea state wave height distribution.To check the capability of these H₅₀-formulae to predict damage evolution of succession of sea states with different wave height distributions, some stability tests with regular and irregular waves have been carried out. After analysing the experimental results, it will be shown how H₅₀-formulae can predict the observed damage independently of the sea state wave height distribution or the succession of sea states.     

Modelling of sand transport under wave-generated sheet flows with a RANS diffusion model

A 1DV-RANS diffusion model is used to study sand transport processes in oscillatory flat-bed/sheet flow conditions. The central aim is the verification of the model with laboratory data and to identify processes controlling the magnitude and direction (‘onshore’/‘offshore’) of the net time-averaged sand transport. The model is verified with a large series of measured net sand transport rates, as collected in different wave tunnels for a range of wave-current conditions and grain sizes. Although not all sheet flow details are represented in the 1DV-model, it is shown that the model is able to give a correct representation of the observed trends in the data with respect to the influence of the velocity, wave period and grain diameter. Also detailed mean sediment flux profiles in the sheet flow layer are well reproduced by the model, including the direction change from ‘onshore’ to ‘offshore’ due to a difference in grain size from 0.34 mm (medium sand) to 0.13 mm (fine sand). A model sensitivity study with a selected series of net transport data shows that the stirring height of the suspended sediment ε_s/w_s strongly controls the magnitude and direction of the net sediment transport. Inclusion of both hindered settling and density stratification appears to be necessary to correctly represent the sand fluxes for waves alone and for waves + a superimposed current. The best agreement with a large dataset of net transport measurements is obtained with the 1DV-RANS model in its original settings using a Prandtl–Schmidt number σ_ρ = 0.5.