2-D numerical analysis of near-field flow at low-crested permeable breakwaters

This paper describes the capability of a numerical model named COrnell BReaking waves And Structures (COBRAS) [Lin, P., Liu, P.L.-F., 1998. A numerical study of breaking waves in the surf zone. Journal of Fluid Mechanics 359, 239–264; Liu, P.L.-F., Lin, P., Chang, K.A., Sakakiyama, T., 1999. Numerical modeling of wave interaction with porous structures. Journal of Waterway, Port, Coastal and Ocean Engineering 125, 322–330, Liu, P.L.-F., Lin, P., Hsu, T., Chang, K., Losada, I.J., Vidal, C., Sakakiyama, T., 2000. A Reynolds averaged Navier–Stokes equation model for nonlinear water wave and structure interactions. Proc. Coastal Structures '99, 169–174] based on the Reynolds Averaged Navier–Stokes (RANS) equations to simulate the most relevant hydrodynamic near-field processes that take place in the interaction between waves and low-crested breakwaters. The model considers wave reflection, transmission, overtopping and breaking due to transient nonlinear waves including turbulence in the fluid domain and in the permeable regions for any kind of geometry and number of layers. Small-scale laboratory tests were conducted in order to validate the model, with different wave conditions and breakwater configurations. In the present study, regular waves of different heights and periods impinging on a wide-crested structure are considered. Three different water depths were tested in order to examine the influence of the structure freeboard. The experimental set-up includes a flow recirculation system aimed at preventing water piling-up at the lee of the breakwater due to overtopping. The applicability and validity of the model are examined by comparing the results of the numerical computations with experimental data. The model is proved to simulate with a high degree of agreement all the studied magnitudes, free surface displacement, pressure inside the porous structure and velocity field. The results obtained show that this model represents a substantial improvement in the numerical modelling of low-crested structures (LCS) since it includes many processes neglected previously by existing models. The information provided by the model can be useful to analyse structure functionality, structure stability, scour and many other hydrodynamic processes of interest.     

Numerical analysis of wave loads for coastal structure stability

The numerical model COBRAS-UC [Losada, I.J., Lara, J.L., Guanche,R., Gonzalez-Ondina, J.M. (2008). Numerical analysis of wave overtopping of rubble mound breakwaters. Coastal Engineering, Vol 55 (1), 47–62.] is used to carry out a two-dimensional analysis of wave induced loads on coastal structures. The model calculates pressure, forces and moments for two different cross-sections corresponding to a low-mound and a conventional rubble-mound breakwater with a crown-wall under regular and irregular incident wave conditions. Predicted results are compared with experimental information provided in Losada et al. [Losada, I.J., Lara, J.L., Guanche,R., Gonzalez-Ondina, J.M. (2008). Numerical analysis of wave overtopping of rubble mound breakwaters. Coastal Engineering, Vol 55 (1), 47–62.] and Lara et al. [Lara, J.L., Losada, I.J., Guanche, R. (2008). “Wave interaction with low mound breakwaters using a RANS model”. Ocean engineering (35), pp 1388–1400; doi:10.1016/j.oceaneng.2008.05.006.] on a 1:20 scale. Good agreement is found, and the differences between both typologies are explained in detail. Additionally, numerical results are also compared with several semi-empirical formulae recommended for design at both the 1:20 model scale and two prototype cross-sections. Results suggest that COBRAS-UC is able to provide realistic stability information that can be used to complete the approach based on currently existing methods and tools.     

Wave Reflection Coefficient Spectrum

1 .IntroductionWavereflectionfromstructuresisanimportantfactorforthedesignofthestructures .Therehavebeenmanyresearchesonwavereflectionwithregularwavesandmonochromaticirregularwaves .Miche( 1 951 )proposedanon dimensionalMichenumberMfornormallyincidentirregularwavesconsideringwavebreakingbecauseofthedeepeningofwavesteepnessontheslopeofbreakwaters .ThereflectioncoefficientisproportionaltoM ,i.e .,Kr ∝M =4g( 2π) 5/ 2tan5/ 2 α(Hsf2p) ( 1 )wheregisthegravitationalacceleration ,Hsthesignifican…     

Random wave-induced scour at the trunk section of a breakwater

This paper provides a stochastic method by which the random wave-induced scour depth at the trunk section of vertical-wall and rubble-mound breakwaters can be derived. Here the formulas for regular wave-induced scour depth provided by Xie [Xie, S.L., 1981. Scouring patterns in front of vertical breakwaters and their influence on the stability of the foundations of the breakwaters. Report. Department of Civil Engineering, Delft University of Technology, Delft, The Netherlands, September, 61 pp.] for vertical-wall breakwater and Sumer and Fredsøe [Sumer, B.M., Fredsøe, J., 2000. Experimental study of 2D scour and its protection at a rubble-mound breakwater. Coast. Eng. 40, 59–87] for rubble-mound breakwater are used. These formulas are combined with describing the waves as a stationary Gaussian narrow-band random process to derive the random wave-induced scour depth. Comparisons are made between the present method and the Sumer and Fredsøe [Sumer, B.M., Fredsøe, J., 2000. Experimental study of 2D scour and its protection at a rubble-mound breakwater. Coast. Eng. 40, 59–87.] random wave scour data for rubble-mound breakwater, as well as the Hughes and Fowler [Hughes, S.A., Fowler, J.A., 1991. Wave-induced scour predictions at vertical walls. ASCE Proc. Conf. Coastal Sediments vol. 91, 1886–1899] random wave scour data and formula for vertical-wall breakwater. A tentative approach to random wave-induced scour at a vertical impermeable submerged breakwater is also suggested.     

Examination of the influence of several factors on longshore current computation with random waves

Influence of various factors affecting the longshore currents induced by obliquely incident random waves is examined through numerical calculation. Seven numerical models for random wave breaking process are found to yield large differences in the wave heights in the surf zone and longshore current velocities. The turbulent eddy viscosity formulation by Larson and Kraus [Larson, M. and Kraus, N.C. (1991): Numerical model of longshore current for bar and trough beaches, J. Waterway, Port, Coastal, and Ocean Eng., ASCE, 117 (4), pp. 326-347.] functions almost equal to that by Battjes [Battjes, J.A. (1975): Modeling of turbulence in the surf zone, Proc. Symp. Modeling Techniques, pp. 1050–1061.], but the formulation by Longuet-Higgins [Longuet-Higgins, M.S. (1970): Longshore current generated by obliquely incident sea waves, 1 and 2, J. Geophys. Res., 75 (33), pp. 6779–6801.] produces excessive diffusion of longshore currents into the offshore zone. The generation and decay process of the surface roller is indispensable in the longshore current analysis. The random wave transformation model called PEGBIS (Parabolic Equation with Gradational Breaker Index for Spectral waves) by Goda [Goda, Y. (2004): A 2-D random wave transformation model with gradational breaker index, Coastal Engineering Journal, JSCE and World Scientific, 46 (1), pp. 1–38.] produced good agreement with several laboratory and field data of longshore currents.     

Cross-shore hydrodynamics within an unsaturated surf zone

This paper concerns the hydrodynamics induced by random waves incident on a steep beach. New experimental results are presented on surface elevation and kinematic probability density functions, cross-shore variation in wave heights, the fraction of broken waves and velocity moments. The surf zone is found to be unsaturated at incident wave frequencies, with a significant proportion of the incident wave energy remaining at the shoreline in the form of bores. Wave heights in both the outer and inner surf zones are best described by a full Rayleigh distribution [Thornton, E.B., Guza, R.T., 1983. Transformation of wave height distribution. J. Geophys. Res. 88, 5925–5938], rather than a truncated Rayleigh distribution as used by Battjes and Janssen (1978) [Battjes, J.A, Janssen, J.P., 1978. Energy loss and setup due to breaking of random waves. Proc. 16th Int. Conf. Coastal Eng. ASCE, New York, pp. 569–588]. A new parametric wave transformation model is outlined which provides explicit expressions for the fraction of broken waves and the energy dissipation rate within the surf zone. On steep beaches, the model appears to offer improved predictive capabilities over the original Battjes and Janssen model. Cross-shore variations in the velocity variance and velocity moments are best described using Linear Gaussian wave theory, with less than 20% of the velocity variance in the inner surf zone due to low frequency energy.     

The BCI criterion for the initiation of breaking process in Boussinesq-type equations wave models

The Breaking Celerity Index (BCI) is proposed as a new wave breaking criterion for Boussinesq-type equations wave propagation models (BTE).The BCI effectiveness in determining the breaking initiation location has been verified against data from different experimental investigations conducted with incident regular and irregular waves propagating along uniform slope [Utku, M. (1999). “The Relative Trough Froude Number. A New Criteria for Wave Breaking”. Ph.D. Dissertation, Dept. of Civil and Enviromental Engineering, Old Dominion University, Norfolk, VA; Gonsalves Veloso dos Reis, M.T.L. (1992). “Characteristics of waves in the surf zone”. MS Thesis, Department of Civil Engineering, University of Liverpool., Liverpool; Lara, J.L., Losada, I.J., and Liu, P.L.-F. (2006). “Breaking waves over a mild gravel slope: experimental and numerical analysis”. Journal of Geophysical Research, VOL 111, C11019] and barred beaches [Tomasicchio, G.R., and Sancho, F. (2002). “On wave induced undertow at a barred beach”. Proceedings of 28th International Conference on Coastal Engineering, ASCE, New York, 557–569]. The considered experiments were carried out in small-scale and large-scale facilities. In addition, one set of data has been obtained by the use of the COBRAS model based upon the Reynolds Averaged Navier Stokes (RANS) equations [Liu, P.L.-F., Lin, P., Hsu, T., Chang, K., Losada, I.J., Vidal, C., and Sakakiyama, T. (2000). “A Reynolds averaged Navier–Stokes equation model for nonlinear water wave and structure interactions”. Proceedings of Coastal Structures ‘99, Balkema, Rotterdam, 169–174; Losada, I.J., Lara, J.L., and Liu, P.L.-F. (2005). “Numerical simulation based on a RANS model of wave groups on an impermeable slope”. Proceedings of Fifth International Symposium WAVES 2005, Madrid].Numerical simulations have been performed with the 1D-FUNWAVE model [Kirby, J.T., Wei, G., Chen, Q., Kennedy, A.B., and Dalrymple, R.A. (1998). “FUNWAVE 1.0 Fully Nonlinear Boussinesq Wave Model Documentation and User's Manual”. Research Report No CACR-98-06, Center for Applied Coastal Research, University of Delaware, Newark]. With regard to the adopted experimental conditions, the breaking location has been calculated for different trigger mechanisms [Zelt, J.A. (1991). “The run-up of nonbreaking and breaking solitary waves”. Coastal Engineering, 15, 205–246; Kennedy, A.B., Chen, Q., Kirby, J.T., and Dalrymple, R.A. (2000). “Boussinesq modeling of wave transformation, breaking and run-up. I: 1D”. Journal of Waterway, Port, Coastal and Ocean Engineering, 126, 39–47; Utku, M., and Basco, D.R. (2002). “A new criteria for wave breaking based on the Relative Trough Froude Number”. Proceedings of 28th International Conference on Coastal Engineering, ASCE, New York, 258–268] including the proposed BCI.The calculations have shown that BCI gives a better agreement with the physical data with respect to the other trigger criteria, both for spilling and plunging breaking events, with a not negligible reduction of the calculation time.     

A numerical model for wave propagation in curvilinear coordinates

Using the perturbation method, a time dependent parabolic equation is developed based on the elliptic mild slope equation with dissipation term. With the time dependent parabolic equation employed as the governing equation, a numerical model for wave propagation including dissipation term in water of slowly varying topography is presented in curvilinear coordinates. In the model, the self-adaptive grid generation method is employed to generate a boundary-fitted and varying spacing mesh. The numerical tests show that the effects of dissipation term should be taken into account if the distance of wave propagation is large, and that the outgoing boundary conditions can be treated more effectively by introduction of the dissipation term into the numerical model. The numerical model is able to give good results of simulating wave propagation for waters of complicatedly boundaries and effectively predict physical processes of wave propagation. Moreover, the errors of the analytical solution deduced by Kirby et al. (1994) [Kirby, J.T., Dalrymple, R.A., Kabu, H., 1994. Parabolic approximation for water waves in conformal coordinate systems. Coastal Engineering 23, 185–213.] from the small-angle parabolic approximation of the mild-slope equation for the case of waves between diverging breakwaters in a polar coordinate system are corrected.     

RANS modelling applied to random wave interaction with submerged permeable structures

This paper is the second part of the work presented by Garcia et al. [Garcia, N., Lara, J.L., Losada, I.J., 2004. 2-D numerical analysis of near-field flow at low-crested breakwaters. Coastal Engineering 51 (10), 991–1020]. In the mentioned paper, flow conditions at low-crested rubble-mound breakwaters under regular wave attack were examined, using a combination of measured data of free surface, bottom pressure and fluid velocities from small-scale experiments and numerical results provided by a VOF-type model (COBRAS) based on the Reynolds-Averaged Navier–Stokes (RANS) equations. This paper demonstrates the capability of the COBRAS model to reproduce irregular wave interaction with submerged permeable breakwaters. Data provided by the numerical model are compared to experimental data of laboratory tests, and the main processes of wave–structure interaction are examined using both experimental and numerical results. The numerical model validation is carried out in two steps. First, the procedure of irregular wave generation is verified to work properly, comparing experimental and numerical data of different cases of irregular wave trains propagating over a flat bottom. Next, the validation of the numerical model for wave interaction with submerged rubble-mound breakwaters is performed through the simulation of small-scale laboratory tests on different incident wave spectra. Results show that the numerical model adequately reproduces the main aspects of the interaction of random waves with submerged porous breakwaters, especially the spectral energy decay at the structure and the spectrum broadening past the structure. The simulations give good results in terms of height envelopes, mean level, spectral shape, root-mean-square height for both free surface displacement and dynamic pressure inside the breakwater. Moreover, large-scale simulations have been conducted, on both regular and irregular incident wave conditions. The overall pattern of the wave interaction with a large-scale submerged breakwater is adequately reproduced by the numerical model. The processes of wave reflection, shoaling and breaking are correctly captured. The good results achieved at a near prototype scale are promising regarding the use of the numerical model for design purposes.     

A probability distribution for depth-limited extreme wave heights in a sea state

A model for the depth-limited distribution of the highest wave in a sea state is presented. The distribution for the extreme wave height is based on a probability density function (pdf) for depth-limited wave height distribution for individual waves [Méndez, F.J., Losada, I.J., Medina, R. 2004. Transformation model of wave height distribution. Coastal Eng, Vol. 50, 97:115.] and considers the correlation between consecutive waves. The model is validated using field data showing a good representation of the extreme wave heights in the surf zone. Some important statistical wave heights are parameterized obtaining useful expressions that can be used in further calculations.     

Lagrangian measurements and modelling of fluid advection in the inner surf and swash zones

New laboratory and field data are presented on fluid advection into the swash zone. The data illustrate the region of the inner surf zone from which sediment can be directly advected into the swash zone during a single uprush, which is termed the advection length. Experiments were conducted by particle tracking in a Lagrangian reference frame, and were performed for monochromatic breaking waves, solitary bores, non-breaking solitary waves and field conditions. The advection length is normalised by the run-up length to give an advection ratio, A, and different advection ratios are identified on the basis of the experimental data. The data show that fluid enters the swash zone from a region of the inner surf zone that can extend a distance seaward of the bore collapse location that is approximately equal to half of the run-up length. This region is about eight times wider than the region predicted by the classical swash solution of Shen and Meyer [Shen, M.C., Meyer, R.E., 1963. Climb of a bore on a beach. Part 3. Runup. Journal of Fluid Mechanics 16, 113–125], as illustrated by Pritchard and Hogg [Pritchard, D., Hogg, A.J., 2005. On the transport of suspended sediment by a swash event on a plane beach. Coastal Engineering 52, 1–23]. Measured advection ratios for periodic waves show no significant trend with Iribarren number, consistent with self-similarity in typical swash flows. The data are compared to recent characteristic solutions of the non-linear shallow water wave (NLSW) equations and both finite difference and finite volume solutions of the NLSW equations.     

Scour below pipelines and around vertical piles in random waves

An approach by which the scour depth and scour width below a fixed pipeline and scour depth around a circular vertical pile in random waves can be derived is presented. Here, the scour depth formulas by Sumer and Fredsøe [ASCE J. Waterw., Port, Coastal Ocean Eng. 116 (1990) 307] for pipelines and Sumer et al. [ASCE J. Waterw., Port, Coastal Ocean Eng. 114 (1992) 599] for vertical piles as well as the scour width formula by Sumer and Fredsøe [The Mechanics of Scour in the Marine Environment, World Scientific, Singapore, 2002] for pipelines combined with describing the waves as a stationary Gaussian narrow-band random process are used to derive the cumulative distribution functions of the scour depths and width. Comparisons are made between the present approach and random wave scour data. Tentative approaches to related random wave scour cases are also suggested.     

Scour around coastal structures: a summary of recent research

This paper summarizes the results of the European Union Marine Science and Technology (EU MAST) III project “Scour Around Coastal Structures” (SCARCOST). The summary is presented under three headings: (1) Introduction; (2) Flow and scour processes with the subheadings: flow and scour processes around vertical cylinders; flow and scour processes at detached breakwaters; flow and scour processes at submerged breakwaters; and the effect of turbulence on sediment transport; and (3) Sediment behaviour close to the structure with the subheadings: field measurement and analysis of wave-induced pore pressures and effective stresses around a bottom seated cylinder; non-linear soil modelling with respect to wave-induced pore pressures and gradients; wave-induced pressures on the bottom for non-linear coastal waves, including also wave kinematics; development of a numerical model (linear soil modelling) to calculate wave-induced pore pressures—the effect of liquefaction on sediment transport; penetration of blocks in non-consolidated fine soil; and cyclic stiffness of loose sand.The paper also includes a discussion of the role of scale effects in laboratory testing and the applicability of the results obtained in supporting engineering design.     

Book review

Geotechnical Aspects of Coastal and Offshore Structures. Edited by Yudhair and A. S. Balasubramaniam. Proceedings of the Symposium on Geotechnical Aspects of Coastal and Offshore Structures, Bangkok, December 14–18,1981. Rotterdam: A. A. Balkema, 1983. 280 pp. $40.00.

The Seaside Naturalist: A Guide to Nature Study at the Seashore. Written and illustrated by D. A. Coulombe. Englewood Cliffs, N.J.: Prentice‐Hall, Inc., 1984. 246 pp. $12.95.

Ichnology, Trace Fossils in Sedimentology and Stratigraphy. By A. A. Ekdale, R. G. Bromley, and S. G. Pemberton. Society of Economic Paleontologists and Mineralogists, Short Course No. 15. Tulsa, Okla.: S.E.P.M., 1984. 317 pp. S.E.P.M. members, $16.00; non‐members, $20.00.     

Field measurements of wave overtopping at the rubble mound breakwater of Rome–Ostia yacht harbour

This paper describes a new station for full-scale measurement of wave overtopping at the Rome yacht harbour rubble mound breakwater in Ostia (Italy) and the results of the successful first measurement campaign carried out during the winter season 2003–2004. The equipment and the research activities were supported by the EU project CLASH, focusing on scale effects for wave overtopping at coastal structures. The site is characterized by a very small tidal range, a long shallow foreshore and depth-limited breaking waves which interact with a shallow sloping porous rock structure. Overtopping water is collected by a steel tank installed on the crown slab behind the parapet wall. The measurement of water level variation inside the tank by means of two pressure transducers allows the calculation of individual overtopping volumes. Incident waves, sea levels and wind are also measured. During seven independent storms, more than 400 individual overtopping events were recorded and about 86 h of valid data are available. This extensive dataset is presented, discussed and then used for comparison with two commonly used overtopping prediction formulae based on small-scale model tests showing their tendency to underestimate the prototype results. A strong correlation between the hourly mean overtopping discharge and corresponding maximum volume is also presented. The paper generally confirms the validity of the approach used in Troch et al. (2004) [Troch, P., Geeraets, J., Van de Walle, B., De Rouck, J., Van Damme, L., Allsop, W., Franco, L., 2004. Full-scale wave overtopping measurements on the Zeebrugge rubble mound breakwater. Coastal Engineering 51, 609–628] for field measurement of wave overtopping.     

Internal generation of waves for extended Boussinesq equations

It is studied whether the mass transport or energy transport is the proper viewpoint for internally generating waves in the extended Boussinesq equations of Nwogu [J. Waterw., Port, Coastal Ocean Eng. 119 (1993) 618–638]. Numerical solutions of the Boussinesq equations with the internal generation of sinusoidal waves show that the energy transport approach yields the required wave amplitude properly while the mass transport approach yields wave amplitude different from the required one by the ratio of phase velocity to energy velocity. The waves which pass through the wave generation point do not cause any numerical distortion while the incident waves are generated. The technique of internal generation of waves shows its capability of generating nonlinear cnoidal waves as well as linear sinusoidal waves.     

Wave interaction with an upright breakwater structure

Breakwaters are often built in coastal waters to facilitate navigation and recreation, both inside and outside regions of the breakwater. This requires that the reflection and transmission characteristics of the structure be both minimized at the same time. This is achieved by a design that will allow dissipation of wave energy by multiple reflection. Such structures will need the knowledge of these characteristics in their design. Model tests were performed on a shallow water breakwater concept of this type to determine the reflection and transmission coefficients. The concept of the breakwater was to reduce both the reflection and transmission of waves. It was found that the breakwater design was effective at certain wave characteristics. Nondimensional loads and local pressures on the breakwater panels are also reported which will facilitate the structural design of such breakwaters.     

Discussion of “A function to determine wavelength from deep in to shallow water based on the length of the cnoidal wave at breaking” by J.P. Le Roux

An empirical modification to the Airy equation for wave celerity reduces to the expression for solitary waves in shallow water whilst retaining its usual form for deep water. The equation yields celerities in reasonable agreement with those for cnoidal waves in intermediate water depths. In this aspect, it is similar to the work described by Le Roux [Le Roux, J.P., 2007. A function to determine wavelength from deep into shallow water based on the length of the cnoidal wave at breaking. Coastal Engineering 54, 770–774]. The empirical modification has been widely applied in computer programs over the past 30 years.     

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.     

堆石防波堤不规则波浪反射系数试验研究

结合物理模型试验,分析斜坡坡度、波陡、相对水深、护面类型和破波参数等因素对堆石防波堤不规则波浪反射系数的影响规律。将常用的Van der Meer公式,Seelig公式,Postma公式和Davison公式计算值和实测值进行比较,并结合试验数据,基于有效波高和平均周期定义的Iribarren数,得出堆石防波堤不规则波浪反射系数经验公式。结果表明,该公式能较好地计算不规则波作用下块石和扭王块体护面堆石防波堤波浪反射系数。