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.     

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.     

Surf zone dynamics simulated by a Boussinesq type model. Part I. Model description and cross-shore motion of regular waves

This is the first of three papers on the modelling of various types of surf zone phenomena. In this first paper, part I, the model is presented and its basic features are studied for the case of regular waves. The model is based on two-dimensional equations of the Boussinesq type and it features improved linear dispersion characteristics, possibility of wave breaking, and a moving boundary at the shoreline. The moving shoreline is treated numerically by replacing the solid beach by a permeable beach characterized by an extremely small porosity. Run-up of nonbreaking waves is verified against the analytical solution for nonlinear shallow water waves. The inclusion of wave breaking is based on the surface roller concept for spilling breakers using a geometrical determination of the instantaneous roller thickness at each point and modelling the effect of wave breaking by an additional convective momentum term. This is a function of the local wave celerity, which is determined interactively. The model is applied to cross-shore motions of regular waves including various types of breaking on plane sloping beaches and over submerged bars. Model results comprise time series of surface elevations and the spatial variation of phase-averaged quantities such as the wave height, the crest and trough elevations, the mean water level, and the depth-averaged undertow. Comparisons with physical experiments are presented. The phaseaveraged balance of the individual terms in the momentum and energy equation is determined by time-integration and quantities such as the cross-sectional roller area, the radiation stress, the energy flux and the energy dissipation are studied and discussed with reference to conventional phase-averaged wave models. The companion papers present cross-shore motions of breaking irregular waves, swash oscillations and surf beats (part II) and nearshore circulations induced by breaking of unidirectional and multidirectional waves (part III).     

A wave-resolving model for nearshore suspended sediment transport

This paper presents a wave-resolving sediment transport model, which is capable of simulating sediment suspension in the field-scale surf zone. The surf zone hydrodynamics is modeled by the non-hydrostatic model NHWAVE (Ma et al., 2012). The turbulent flow and suspended sediment are simulated in a coupled manner. Three effects of suspended sediment on turbulent flow field are considered: (1) baroclinic forcing effect; (2) turbulence damping effect and (3) bottom boundary layer effect. Through the validation with the laboratory measurements of suspended sediment under nonbreaking skewed waves and surfzone breaking waves, we demonstrate that the model can reasonably predict wave-averaged sediment profiles. The model is then utilized to simulate a rip current field experiment (RCEX) and nearshore suspended sediment transport. The offshore sediment transport by rip currents is captured by the model. The effects of suspended sediment on self-suspension are also investigated. The turbulence damping and bottom boundary layer effects are significant on sediment suspension. The suspended sediment creates a stably stratified water column, damping fluid turbulence and reducing turbulent diffusivity. The suspension of sediment also produces a stably stratified bottom boundary layer. Thus, the drag coefficient and bottom shear stress are reduced, causing less sediment pickup from the bottom. The cross-shore suspended sediment flux is analyzed as well. The mean Eulerian suspended sediment flux is shoreward outside the surf zone, while it is seaward in the surf zone.     

Quasi-3D modelling of nearshore currents

Existing concepts of wave-induced nearshore current models, in the cross-shore vertical plane (2DV) and depth-integrated (2DH), are combined to a quasi-3D mathematical model. This combination is tested for reproducing correct results in 2DV and 2DH situations. The importance of the various contributions to the wave-induced secondary circulation in the vertical plane is investigated for realistic parameter ranges, which leads to the conclusion that both the non-breaking and the breaking fraction of a random wave field in the surf zone generate important secondary currents.Additional computations show the relevance of a 3D-approach of nearshore currents, even in seemingly simple situations like a plane sloping beach with obliquely incident waves.     

沿岸流中混合系数的实验研究

进行了在规则波作用下破波带内水平混合系数测量的物理模型实验。通过在沿岸流流场中投放墨水点源和采用CCD摄像机摄像,测量了点源扩散过程。利用水深平均二维扩散方程近似解析解得到了由实验结果计算混合系数的方法。实验结果表明:扩散系数仅依赖于当地水深,与波浪参数(周期和波高)无关。沿岸流沿水深变化产生的离散作用导致顺流方向(纵向)混合系数远大于横流方向(横向)混合系数。横流方向混合系数中由波浪产生的扩散系数占总扩散系数约40%,其余为波浪破碎引起的湍流产生的扩散系数。     

Laboratory study of wave and turbulence velocities in a broad-banded irregular wave surf zone

The characteristics of wave and turbulence velocities created by a broad-banded irregular wave train breaking on a 1:35 slope were studied in a laboratory wave flume. Water particle velocities were measured simultaneously with wave elevations at three cross-shore locations inside the surf zone. The measured data were separated into low-frequency and high-frequency time series using a Fourier filter. The measured velocities were further separated into organized wave-induced velocities and turbulent velocity fluctuations by ensemble averaging. The broad-banded irregular waves created a wide surf zone that was dominated by spilling type breakers. A wave-by-wave analysis was carried out to obtain the probability distributions of individual wave heights, wave periods, peak wave velocities, and wave-averaged turbulent kinetic energies and Reynolds stresses. The results showed that there was a consistent increase in the kurtosis of the vertical velocity distribution from the surface to the bottom. The abnormally large downward velocities were produced by plunging breakers that occurred from time to time. It was found that the mean of the highest one-third wave-averaged turbulent kinetic energy values in the irregular waves was about the same as the time-averaged turbulent kinetic energy in a regular wave with similar deep-water wave height to wavelength ratio. It was also found that the correlation coefficient of the Reynolds stress varied strongly with turbulence intensity. Good correlation between u′ and w′ was obtained when the turbulence intensity was high; the correlation coefficient was about 0.3–0.5. The Reynolds stress correlation coefficient decreased over a wave cycle, and with distance from the water surface. Under the irregular breaking waves, turbulent kinetic energy was transported downward and landward by turbulent velocity fluctuations and wave velocities, and upward and seaward by the undertow. The undertow in the irregular waves was similar in vertical structure but lower in magnitude than in regular waves, and the horizontal velocity profiles under the low-frequency waves were approximately uniform.     

Numerical study on vertical structures of undertow inside and outside the surf zone

Nearshore shoaling and breaking waves can drive a complex circulation system of wave-induced currents. In the cross-shore direction, the local vertical imbalance between the gradient of radiation stress and that of pressure due to the setup drives an offshore flow near the bottom, called ‘undertow’, which plays a significant role in the beach profile evolution and the structure stability in coastal regions. A 1DV undertow model was developed based on the relationship between the turbulent shear stress and t...     

Experimental study of the instabilities of alongshore currents on plane beaches

A laboratory experiment on alongshore currents was conducted for two plane beaches, with gradients 1:40 and 1:100, to investigate the instability of alongshore currents. Complicated and strongly unstable alongshore current motions were observed. In order to clearly examine the spatial and temporal variations of the shear instability of the currents, digital images from a charge-coupled device (CCD) recorded the deformations of dye batches released in the surf zone. Some essential characteristics of the shear instability were obtained from analyses of images showing the temporal variation of the dye patches.A high-resolution spectral analysis technique (the maximum entropy method, or MEM) was used to analyze the dominant frequency of the observed oscillation, along with the trigonometric regression method for determining the variations of the oscillation strength in the cross-shore direction. The propagation speed of the dye patch was obtained by tracking the movement over time of fixed locations in the dye patch, such as its peak, in the longshore direction. This data was then fitted linearly.Alongshore and cross-shore velocity time series acquired from sensors showed clearly that large-amplitude, long-period (about 50 s or 100 s) oscillations were present for all sensors deployed in the cross-shore direction under regular and irregular wave conditions. The analysis found that the maximum shear wave amplitude was approximately one-sixth of the maximum for the mean alongshore current, and occurred approximately at the position of the maximum of the mean alongshore current for irregular waves. The spatial structure of the shear waves was studied by analyzing collected images of the dye patches. The phase velocity of the meandering movements was obtained by measuring the magnitude of the oscillations of the dye patches in the alongshore direction with respect to time. The results suggest that the propagation speed of the shear instability was approximately one-half to three-quarters of the maximum mean longshore current for both regular and irregular waves.Linear instability analysis theory was applied to the characteristics of alongshore current instability, which suggested that there are two instability modes related to the observed oscillations: the frontshear mode observed for the 1:100 slope, and the backshear mode observed for the 1:40 slope. Theoretical analyses agreed with the experimental results in both cases. The velocity profile of the mean longshore current was found to affect the instability mode significantly, leading to further investigations on the influence of the velocity profiles and to provide support for the above conclusions.     

Magnitude and cross-shore distribution of bed return flow measured on natural beaches

Field measurements of cross-shore currents 0.25 m from the bed were made on two natural beaches under a range of incident wave conditions. The results indicated the presence of a relatively strong, offshore-directed mean current, both within and seaward of the surf zone. Typical velocities within the surf zone were of the order of 0.2–0.3 m/s. This bed return flow, or “undertow”, represents a mass conservation response, returning water seaward that was initially transported onshore in the upper water column, primarily above the trough of the incident waves. The measurements demonstrated that the bed return flow velocity increases with the incident wave height. In addition, the crossshore distribution of the bed return flow is characterised by a mid-surf zone maximum, which exhibits a strong decrease in velocity towards the shoreline and a more gradual decay in the offshore direction. Several bed return flow models based on mass continuity were formulated to predict the cross-shore distribution of the bed return flow under an irregular wave field and were compared with the field data. Best agreement was obtained using shallow water linear wave theory, after including the mass transport associated with unbroken waves. The contribution of the unbroken waves enables net offshore-directed bottom currents to persist outside the region of breaking waves, providing a mechanism, other than rip currents, to transport sediment offshore beyond the surf zone.     

Breakpoint generated surf beat induced by bichromatic wave groups

This paper presents new experimental data on 2-D surf beat generation by a time-varying breakpoint induced by bichromatic wave groups. The experimental investigation covers a broad range of wave amplitudes, short wave frequencies, group frequencies and modulation rates. The data include measurements of incident and outgoing wave amplitudes, breakpoint position, shoreline run-up and the cross-shore structure of both the short and long wave motion. Surf beat generation is shown to be in good agreement with theory [Symonds, G., Huntley, D.A., Bowen, A.J., 1982. Two dimensional surf beat: long wave generation by a time-varying breakpoint. J. Geophys. Res. 87, 492–498]. In particular, surf beat generation is dependent on the normalised surf zone width, which is a measure of the phase relationship between the seaward and shoreward breakpoint forced long waves, and linearly dependent on the short wave amplitude. The cross-shore structure of the long wave motion is also consistent with theory; at maximum and minimum surf beat generation, the mean breakpoint coincides with the nodal and anti-nodal points, respectively, for a free long wave standing at the shoreline. A numerical solution, using measured data as input, additionally shows that the phase relationship between the incident bound long wave and the outgoing breakpoint forced wave is consistent with the time-varying breakpoint mechanism.     

Hyperbolic Mild Slope Equations with Inclusion of Amplitude Dispersion Effect： Regular Waves

A new form of hyperbolic mild slope equations is derived with the inclusion of the amphtude dispersion of nonlinear waves. The effects of including the amplitude dispersion effect on the wave propagation are discussed. Wave breaking mechanism is incorporated into the present model to apply the new equations to surf zone. The equations are solved nu- merically for regular wave propagation over a shoal and in surf zone, and a comparison is made against measurements. It is found that the inclusion of the amplitude dispersion can also improve model＇ s performance on prediction of wave heights around breaking point for the wave motions in surf zone.     

Process-based model for nearshore hydrodynamics, sediment transport and morphological evolution in the surf and swash zones

Hydrodynamics and sediment transport in the nearshore zone were modeled numerically taking into account turbulent unsteady flow. The flow field was computed using the Reynolds Averaged Navier–Stokes equations with a k–ε turbulence closure model, while the free surface was tracked using the Volume-Of-Fluid technique. This hydrodynamical model was supplemented with a cross-shore sediment transport formula to calculate profile changes and sediment transport in the surf and swash zones. Based on the numerical solutions, flow characteristics and the effects of breaking waves on sediment transport were studied. The main characteristic of breaking waves, i.e. the instantaneous sediment transport rate, was investigated numerically, as was the spatial distribution of time-averaged sediment transport rates for different grain sizes. The analysis included an evaluation of different values of the wave friction factor and an empirical constant characterizing the uprush and backwash. It was found that the uprush induces a larger instantaneous transport rate than the backwash, indicating that the uprush is more important for sediment transport than the backwash. The results of the present model are in reasonable agreement with other numerical and physical models of nearshore hydrodynamics. The model was found to predict well cross-shore sediment transport and thus it provides a tool for predicting beach morphology change.     

C.C. Mei, The applied dynamics of ocean surface waves, Wiley, New York (1983) xvii + 740 pp., £ stg. 66.75.

A theoretical model is developed for wave heights and set-up in a surf zone. In the time-averaged equations of energy and momentum the energy flux, radiation stress and energy dissipation are determined by simple approximations which include the effect of the surface roller in the breaker and the actual shape of the waves. To the first approximation the roller represents a volume of water moving with the wave speed. It significantly changes both energy flux and radiation stress of the surf zone waves. The equations of energy and momentum are solved simultaneously to give the wave height variation and the set-up. Comparison with measurements shows good agreement. Also the transitions immediately after breaking are analyzed and shown to be in accordance with the above-mentioned ideas and results.     

Wave heights and set-up in a surf zone

A theoretical model is developed for wave heights and set-up in a surf zone. In the time-averaged equations of energy and momentum the energy flux, radiation stress and energy dissipation are determined by simple approximations which include the effect of the surface roller in the breaker and the actual shape of the waves. To the first approximation the roller represents a volume of water moving with the wave speed. It significantly changes both energy flux and radiation stress of the surf zone waves. The equations of energy and momentum are solved simultaneously to give the wave height variation and the set-up. Comparison with measurements shows good agreement. Also the transitions immediately after breaking are analyzed and shown to be in accordance with the above-mentioned ideas and results.     

Laboratory study of wave and turbulence characteristics in narrow-band irregular breaking waves

Wave elevations and water particle velocities were measured in a laboratory surf zone created by the breaking of a narrow-band irregular wave train on a 1/35 plane slope. The incident waves form wave groups that are strongly modulated. It is found that the waves that break close to the shoreline generally have larger wave-height-to-water-depth ratios before breaking than the waves that break farther offshore. After breaking, the wave-height-to-water-depth ratio for the individual waves approaches a constant value in the inner surf zone, while the standard deviation of the wave period increases as the still water depth decreases. In the outer surf zone, the distribution of the period-averaged turbulent kinetic energy is closely correlated to the initial wave heights, and has a wider variation for narrow-band waves than for broad-band waves. In the inner surf zone, the distribution of the period-averaged turbulent kinetic energy is similar for narrow-band waves and broad-band waves. It is found that the wave elevation and turbulent kinetic energy time histories for the individual waves in a wave group are qualitatively similar to those found in a spilling regular wave. The time-averaged transport of turbulent kinetic energy by the ensemble-averaged velocity and turbulence velocity under the irregular breaking waves are also consistent with the measurements obtained in regular breaking waves. The experimental results indicate that the shape of the incident wave spectrum has a significant effect on the temporal and spatial variability of wave breaking and the distribution of turbulent kinetic energy in the outer surf zone. In the inner surf zone, however, the distribution of turbulent kinetic energy is relatively insensitive to the shape of the incident wave spectrum, and the important parameters are the significant wave height and period of the incident waves, and the beach slope.     

Modeling wave�Ccurrent bottom boundary layers beneath shoaling and breaking waves

The boundary layer characteristics beneath waves transforming on a natural beach are affected by both waves and wave-induced currents, and their predictability is more difficult and challenging than for those observed over a seabed of uniform depth. In this research, a first-order boundary layer model is developed to investigate the characteristics of bottom boundary layers in a wave–current coexisting environment beneath shoaling and breaking waves. The main difference between the present modeling approach and previous methods is in the mathematical formulation for the mean horizontal pressure gradient term in the governing equations for the cross-shore wave-induced currents. This term is obtained from the wave-averaged momentum equation, and its magnitude depends on the balance between the wave excess momentum flux gradient and the hydrostatic pressure gradient due to spatial variations in the wave field of propagating waves and mean water level fluctuations. A turbulence closure scheme is used with a modified low Reynolds number k-ε model. The model was validated with two published experimental datasets for normally incident shoaling and breaking waves over a sloping seabed. For shoaling waves, model results agree well with data for the instantaneous velocity profiles, oscillatory wave amplitudes, and mean velocity profiles. For breaking waves, a good agreement is obtained between model and data for the vertical distribution of mean shear stress. In particular, the model reproduced the local onshore mean flow near the bottom beneath shoaling waves, and the vertically decreasing pattern of mean shear stress beneath breaking waves. These successful demonstrations for wave–current bottom boundary layers are attributed to a novel formulation of the mean pressure gradient incorporated in the present model. The proposed new formulation plays an important role in modeling the boundary layer characteristics beneath shoaling and breaking waves, and ensuring that the present model is applicable to nearshore sediment transport and morphology evolution.     

Hyperbolic Mild Slope Equations with Inclusion of Amplitude Dispersion Effect: Random Waves

New hyperbolic mild slope equations for random waves are developed with the inclusion of amplitude dispersion. The frequency perturbation around the peak frequency of random waves is adopted to extend the equations for regular waves to random waves. The nonlinear effect of amplitude dispersion is incorporated approximately into the model by only considering the nonlinear effect on the carrier waves of random waves, which is done by introducing a representative wave amplitude for the carrier waves. The computation time is greatly saved by the introduction of the representative wave amplitude. The extension of the present model to breaking waves is also considered in order to apply the new equations to surf zone. The model is validated for random waves propagate over a shoal and in surf zone against measurements.