The influence of seaward boundary conditions on swash zone hydrodynamics

The influence of the seaward boundary condition on the internal swash hydrodynamics is investigated. New numerical solutions of the characteristics form of the nonlinear shallow-water equations are presented and applied to describe the swash hydrodynamics forced by breaking wave run-up on a plane beach. The solutions depend on the specification of characteristic variables on the seaward boundary of the swash zone, equivalent to prescribing the flow depth or the flow velocity. It is shown that the analytical solution of Shen and Meyer [Shen, M.C., Meyer, R.E., 1963. Climb of a bore on a beach. Part 3. Runup. J. Fluid Mech. 16, 113–125] is a special case of the many possible solutions that can describe the swash flow, but one that does not appear appropriate for practical application for real waves. The physical significance of the boundary conditions is shown by writing the volume and momentum fluxes in terms of the characteristic variables. Results are presented that illustrate the dependence of internal flow depth and velocity on the boundary condition. This implies that the internal swash hydrodynamics depend on the shape and wavelength of the incident bore, which differs from the hydrodynamic similarity inherent in the analytical solution. A solution appropriate for long bores is compared to laboratory data to illustrate the difference from the analytical solution. The results are important in terms of determining overwash flows, flow forces and sediment dynamics in the run-up zone.     

Alongshore fluid motions in the swash zone of a sandy and gravel beach

Field experiments were conducted on a low-gradient, high-energy sandy beach (Truc Vert, France) and a steep, low-energy gravel beach (Slapton, UK) to examine alongshore-directed currents within the swash zone. At Truc Vert, data were collected over 33 tidal cycles with offshore significant wave heights of 1–4 m and periods of 5–12 s. At Slapton data were collected during 12 tides with wave heights of 0.3–1 m and periods of 4–9 s. The swash motion was predominantly at infragravity frequencies at Truc Vert and incident frequencies at Slapton.     

A Lagrangian model for boundary layer growth and bed shear stress in the swash zone

A new model for the boundary layer development and associated skin friction coefficients and shear stress within the swash zone is presented. The model is developed within a Lagrangian reference frame, following fluid trajectories, and can be applied to both laminar flow and smooth turbulent flow. The model is based on the momentum integral approach for steady, flat-plate boundary layers, with appropriate modifications to account for the unsteady flow regime and flow history. The model results are consistent with previous measurements of bed shear stress and skin friction coefficients within the swash zone. These indicate strong temporal and spatial variation throughout the swash cycle, and a clear distinction between the uprush and backwash phase. This variation has been previously attributed the unsteady flow regime and flow history effects, both of which are accounted for in the new model. Fluid particle trajectories and velocity are computed using the non-linear shallow water wave equations and the boundary layer growth across the entire swash zone is estimated. Predictions of the bed shear stress and skin friction coefficients agree reasonably well with direct bed shear stress measurements reported by Barnes et al. (Barnes, M.P., O’Donaghue, T., Alsina, J.M., Baldock, T.E., 2009. Direct bed shear stress measurements in bore-driven swash. Coastal Engineering 56 (8), 853–867) and, for a given flow velocity, give stresses which are consistent with the bias toward uprush sediment transport which has consistently been observed in measurements. The data and modelling suggest that the backwash boundary layer is initially laminar, which results in the late development of significant bed shear during the backwash, with a transition to a turbulent boundary layer later in the backwash. A new conceptual model for the boundary layer structure at the leading edge of the swash is proposed, which accounts for both the no-slip condition at the bed and the moving wet–dry interface. However, further development of the Lagrangian Boundary Layer Model is required in order to include bore-generated turbulence and to account for variable roughness and mobile beds.     

Sediment suspension events in the inner surf and swash zone. Measurements in large-scale and high-energy wave conditions

The present study presents a database of hydrodynamic properties and suspended sediment concentration collected within the inner surf and swash zones aiming to improve the current understanding of the sediment dynamics occurring within the beach area closest to the shoreline.     

Numerical study on wave-induced filtration flow across the beach face and its effects on swash zone sediment transport

A numerical model, coupling an analysis of beach groundwater flow with an analysis of swash wave motion over a uniform slope, is presented. Model calculations are performed to investigate the variations of swash-induced filtration flows across the beach face for different input parameters. Swash zone sediment transport under the influence of such filtration flow across the beach face is investigated through modification of effective weight of sediment particle and modification of swash boundary layer thickness. These effects are quantified based on a bed load transport model with a modified Shields parameter.     

Measurement of wave-by-wave bed-levels in the swash zone

A technique is described to observe and quantify wave-by-wave bed-level changes in the swash zone. The ultrasonic instrument system is non-contact with the beach face surface being measured and the sensors remain outside of the fluid flows causing sediment movement. Sensor resolution combined with the electronic noise inherent within a digital network data-logging system results in a (conservative) measurement accuracy of ± 1 mm, equating to a couple of sand grain diameters in height. Illustrative field results demonstrate the practical use of the instrumentation, and a simple data pre-processing method to separate swashes and intervening bed-level ‘events’ is discussed. These example data reveal rather complex fluctuations of the bed observed over time periods of minutes to hours. Rather strikingly, gross bed-level changes per wave are revealed to be up to many orders of magnitude larger than the observed net rate of beach face evolution. It is outlined how observations of successive bed-level changes at multiple locations within a dense grid, combined with a consideration of sediment continuity, will now enable the total net sediment transported per uprush–backwash to be quantified.     

Large-scale experiments on beach profile evolution and surf and swash zone sediment transport induced by long waves,wave groups and random waves

New large-scale laboratory data are presented on the influence of long waves, bichromatic wave groups and random waves on sediment transport in the surf and swash zones. Physical model testing was performed in the large-scale CIEM wave flume at UPC, Barcelona, as part of the SUSCO (swash zone response under grouping storm conditions) experiment in the Hydralab III program (Vicinanza et al., 2010). Fourteen different wave conditions were used, encompassing monochromatic waves, bichromatic wave groups and random waves. The experiments were designed specifically to compare variations in beach profile evolution between monochromatic waves and unsteady waves with the same mean energy flux. Each test commenced with approximately the same initial profile. The monochromatic conditions were perturbed with free long waves, and then subsequently substituted with bichromatic wave groups with different bandwidth and with random waves with varying groupiness. Beach profile measurements were made at half-hourly and hourly intervals, from which net cross-shore transport rates were calculated for the different wave conditions. Pairs of experiments with slightly different bandwidth or wave grouping show very similar net cross-shore sediment transport patterns, giving high confidence to the data set. Consistent with recent small-scale experiments, the data clearly show that in comparison to monochromatic conditions the bichromatic wave groups reduce onshore transport during accretive conditions and increase offshore transport during erosive conditions. The random waves have a similar influence to the bichromatic wave groups, promoting offshore transport, in comparison to the monochromatic conditions. The data also indicate that the free long waves promote onshore transport, but the conclusions are more tentative as a result of a few errors in the test schedule and modifications to the setup which reduced testing time. The experiments suggest that the inclusion of long wave and wave group sediment transport is important for improved near-shore morphological modeling of cross-shore beach profile evolution, and they provide a very comprehensive and controlled series of tests for evaluating numerical models. It is suggested that the large change in the beach response between monochromatic conditions and wave group conditions is a result of the increased significant and maximum wave heights in the wave groups, as much as the presence of the forced and free long waves induced by the groupiness. The equilibrium state model concept can provide a heuristic explanation of the influence of the wave groups on the bulk beach profile response if their effective relative fall velocity is larger than that of monochromatic waves with the same incident energy flux.     

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.     

On low-frequency waves in the surf and swash

Low-frequency waves in the surf and swash zones on various beach slopes are discussed using numerical simulations. Simulated surface elevations of both primary waves and low-frequency waves across the surf zone were first compared with experimental data and good agreement found. Low-frequency wave characteristics are then discussed in terms of their physical nature and their relationship to the primary wave field on a series of sea bottom slopes. Unlike primary waves, low-frequency wave energy increases towards the shoreline. Low-frequency waves in the surf and swash are a function of incident waves and the sea bottom slope and hence the saturation level of the surf zone. Wave energy on a gently sloping beach is dominated by low-frequency waves while primary waves play a significant role on a steep beach. Low-frequency wave radiation from the surf zone on a given beach depends on primary wave frequency and beach slope. However, a very poor correlation was found between surf similarity parameter and low-frequency wave radiation.     

海滩冲流带高频振动地形动力过程分析

在粤东汕尾后江湾海滩冲流带布设2条观测剖面,共计6个观测点,对滩面冲流带在约一个潮周期内的高频振动进行了观测,取得采样频率分别为1min/次和6min/次的滩面数据各1组.结合同步观测的碎波带波浪潮汐数据,分析探讨了海滩高频振动特征.分析认为在涌浪条件下,滩面高频振动的日内变化主要受到潮位变化过程的控制,涨潮堆积,落潮侵蚀.利用交叉谱分析的结果表明滩面高程变化滞后于潮位变化.滩面下部比上部振动幅度大,变化复杂,滩角脊部比凹部活跃.波群对滩面高频振动有显著影响,特别是波高大于有效波高的波群.滩面高频振动没有表现出明显的泥沙逐渐向陆地堆积过程,有一定的振动周期.     

Wave run-up on slender piles in design conditions — Model tests and design rules for offshore wind

Wave run-up on foundations is a very important factor in the design of entrance platforms for offshore wind turbines. When the Horns Reef 1 wind turbine park in Denmark was designed the vertical wave run-up phenomenon was not well known in the industry, hence not sufficiently considered in the design of Horns Reef 1. As a consequence damage was observed on the platforms. This has been the situation for several sites and design tools for platform loads are lacking. As a consequence a physical model test study was initiated at Aalborg University to clarify wave run-up on cylindrical piles for different values of diameter to water depth ratios (D/h) and different wave heights to water depth ratios (H/h) for both regular and irregular waves. A calculation model is calibrated based on stream function theory for crest kinematics and velocity head stagnation theory. Due to increased velocities close to the pile an empirical factor is included on the velocity head. The evaluation of the calculation model shows that an accurate design rule can be established even in breaking wave conditions. However, calibration of a load model showed that it was necessary to increase the run-up factor on the velocity head by 40% to take into account the underestimation of run-up for breaking or nearly breaking waves given that they produce thin run-up wedges and air entrainment, two factors not coped with by the measurement system.     

Observation and measurement of the bottom boundary layer flow in the prebreaking zone of shoaling waves

In this paper, the characteristics of the bottom boundary layer flow induced by nonlinear, asymmetric shoaling waves, propagating over a smooth bed of 1/15 uniform slope, is experimentally investigated. Flow visualization technique with thin-layered fluorescent dye was first used to observe the variation of the flow structure, and a laser Doppler velocimeter was then employed to measure the horizontal velocity, U.The bottom boundary layer flow is found to be laminar except within a small region near the breaking point. The vertical distribution of the phase-averaged velocity U at each phase is non-uniform, which is directly affected by the mean velocity, . The magnitude of increases from zero at the bottom to a local positive maximum at about z/δ2.02.5 (where z is the height above the sloping bottom and δ is the Stokes layer thickness), then decreases gradually to zero at z/δ6.07.0 approximately, and finally becomes negative as z/δ increases further. Moreover, as waves propagate towards shallower water, the rate of increase in the maximum onshore oscillating velocity component is greater than that of the offshore counterpart except near the breaking point. The free stream velocities in the profiles of the maximum onshore and offshore oscillating velocity components, and are found to appear at z/δ≥6.0. This implies that, if the Stokes layer thickness is used as a length scale, the non-dimensionalized boundary layer thickness remains constant in the pre-breaking zone. Although is greater than and the asymmetry of the maximum free stream velocities (i.e. ) increases with decrease of water depth, a universal similar profile can be established by plotting z/δ versus ( ) or ( ). The final non-dimensional profile is symmetric and unique for the distributions of the maximum onshore and offshore oscillating velocity components within the bottom boundary layer, which are induced by nonlinear, asymmetric shoaling waves crossing the pre-breaking zone.     

珠江三角洲快速客船船行波模拟及其爬高研究

采用不规则波生波技术模拟珠江三角洲快速客船实测船行波,分别在波浪水槽和波浪水池中进行了船行波模拟试验研究。确定了不同波要素、不同入射角以及不同护岸坡比、平台宽窄、护面层结构等对船行波在斜坡护岸上的爬高影响,取得了船行波在斜坡式护岸上爬高的经验关系式。为珠江三角洲内河航道快速客船船行波作用下护坡的设计提供参考依据。     

Dynamics of cobbles in the shoaling region of a surf zone

The motion of large bottom particles (cobbles/mines) was studied in the laboratory under simulated surf conditions. A series of experiments was conducted in a large wave tank, 32×0.9×1.8 m, equipped with a computer-controlled wave maker and a sloping beach. As a first step, a solid impermeable beach with artificial roughness was used in the experiments. Cobbles of different size were placed along the floor and their evolution with time was studied and compared with the model predictions. Onshore and offshore mean motions of cobbles, as well as steady oscillations with zero mean displacement, were observed for different conditions. To explain the results of observations a theoretical model was advanced. The model takes into account all main governing parameters (size and density of cobbles, bottom slope, dynamic and static friction at the bottom, background flow characteristics, etc.). Standard parameterizations were used for a pressure accelerating term, drag, lift and other nonlinear forces. For the range of parameters used in the experiments, satisfactory agreement between the measured and calculated values of the cobble displacements as a function of time was obtained. The model is practically insensitive to the vertical accelerating pressure term but sensitive to the dynamic and static friction. One of the most important variables in the model, which is known with the least accuracy, is the virtual mass coefficient for disk-shaped cobbles moving with variable velocity along a solid boundary.     

On the evolution and run-up of breaking solitary waves on a mild sloping beach

This paper presents new laboratory experiments carried out in a supertank (300 m × 5 m × 5.2 m) of breaking solitary waves evolution on a 1:60 plane beach. The measured data are employed to re-examine existing formulae that include breaking criterion, amplitude evolution and run-up height. The properties of shoreline motion, underwater particle velocity and scale effect on run-up height are briefly discussed. Based on our analyses, it is evidently found that there exist five zones during a wave amplitude evolution course on the present mild slope. A simple formula which is capable of predicting maximum run-up height for a breaking solitary wave on a uniform beach with a wide range of beach slope (1:15–1:60) is also proposed. The calculated results from the present model agree favorably with available laboratory data, indicating that our method is compatible with other predictive models.     

Surf zone dynamics simulated by a Boussinesq type model. Part II: surf beat and swash oscillations for wave groups and irregular waves

This is the second of three papers on the modelling of various types of surf zone phenomena. In the first paper the general model was described and it was applied to study cross-shore motion of regular waves in the surf zone. In this paper, part II, we consider the cross-shore motion of wave groups and irregular waves with emphasis on shoaling, breaking and runup as well as the generation of surf beats. These phenomena are investigated numerically by using a time-domain Boussinesq type model, which resolves the primary wave motion as well as the long waves. As compared with the classical Boussinesq equations, the equations adopted here allow for improved linear dispersion characteristics and wave breaking is modelled by using a roller concept for spilling breakers. The swash zone is included by incorporating a moving shoreline boundary condition and radiation of short and long period waves from the offshore boundary is allowed by the use of absorbing sponge layers. Mutual interaction between short waves and long waves is inherent in the model. This allows, for example, for a general exchange of energy between triads rather than a simple one-way forcing of bound waves and for a substantial modification of bore celerities in the swash zone due to the presence of long waves. The model study is based mainly on incident bichromatic wave groups considering a range of mean frequencies, group frequencies, modulation rates, sea bed slopes and surf similarity parameters. Additionally, two cases of incident irregular waves are studied. The model results presented include transformation of surface elevations during shoaling, breaking and runup and the resulting shoreline oscillations. The low frequency motion induced by the primary-wave groups is determined at the shoreline and outside the surf zone by low-pass filtering and subsequent division into incident bound and free components and reflected free components. The model results are compared with laboratory experiments from the literature and the agreement is generally found to be very good. Finally the paper includes special details from the breaker model: time and space trajectories of surface rollers revealing the breakpoint oscillation and the speed of bores; envelopes of low-pass filtered radiation stress and surface elevation; sensitivity of surf beat to group frequency, modulation rate and bottom slope is investigated. Part III of this work (Sørensen et al., 1998) presents nearshore circulations induced by the breaking of unidirectional and multi-directional waves.     

Experimental investigations of the bed evolution in wave flumes: Performance of 2D and 3D optical systems

Optical systems can provide very accurate measurements of bottom morphology in wave flumes. However, it is often necessary, e.g. when laser scanners are used, to stop the experiments and disturbing significantly the sandy bed itself, by emptying the flume. In the present work measurement strategies based on computer vision techniques which permit measurements also in the presence of water are applied in wave flumes at small and at large scale. Such techniques, based on the use of structured light, are demonstrated to be able to perform measurements of 2D and 3D bed evolution also in a very active area, such as the swash zone, where the alternating presence and absence of water makes it difficult to recover the bed morphology in a dynamic way.     

Modelling of sediment movement in the surf and swash zones

Under the action of marine currents, non-cohesive sediments evolve by bed-load, by saltation or suspension depending on their granulometry. Several authors have considered that the movement of sediment...     

Morphodynamics and cobbles behavior in and near the surf zone

The evolution of an initially flat sandy slope and the dynamics of large objects (cobbles/mines) emplaced on it are studied in a laboratory wave tank under simulated surf conditions. Upon initiation of wave forcing, the initially flat beach undergoes bedform changes before reaching a quasi-steady morphology characterized by a system of sand ripples along the slope and a large bar near the break point. Although the incoming wave characteristics are held fixed, the bottom morphology never reaches a strict steady state, but rather slowly changes due to the migration of ripples and bar transformation. When the wave characteristics are changed, the bedform adjusts to a new quasi-steady state after a suitable adjustment time. Studies conducted by placing model cobbles/mines on the evolving sandy bottom subjected to wave forcing show four distinct scenarios: (i) periodic cobble oscillations with zero mean displacement and small scour around the cobbles, (ii) mean onshore motion of relatively light cobbles, (iii) periodic burial of relatively heavy cobbles when their sizes are comparable to those of sand ripples, and (iv) the burial of relatively large cobbles under the bar, when the bar migrates due to changes of incoming waves. Quantitative data on the characteristics and dynamics of the bedform, including ripple-formation front propagating down the slope, ripple growth and drift, and flow around ripples, are presented. Physical explanations are provided for the observations.