Large Eddy Simulation for Wave Breaking in the Surf Zone

In this paper, (he large eddy simulation method is used combined with the marker and cell method to study the wave propagation or shoaling and breaking process. As wave propagates into shallow water, the shoaling leads lo the increase of wave height, and then at a certain position, the wave will be breaking. The breaking wave is a powerful agent for generating turbulence, which plays an important role in most of the fluid dynamic processes throughout the surf zone, such as transformation of wave energy, generation of near-shore current and diffusion of materials. So a proper numerical model for describing the turbulence effect is needed. In this paper, a revised Smagorinsky subgrid-scale mode! is used to describe the turbulence effect. The present study reveals that the coefficient of the Smagorinsky model for wave propagation or breaking simulation may be taken as a varying function of the water depth and distance away from the wave breaking point. The large eddy simulation model presented in this pape     

Large Eddy Simulation for Plunge Breaker and Sediment Suspension

Breaking waves are a powerful agent for generating turbulence that plays an important role in many fluid dynamical processes, particularly in the mixing of materials. Breaking waves can dislodge sediment and throw it into suspension, which will then be carried by wave-induced steady current and tidal flow. In order to investigate sediment suspension by breaking waves, a numerical model based on large-eddy-simulation (LES) is developed. This numerical model can be used to simulate wave breaking and sediment suspension. The model consists of a free-surface model using the surface marker method combined with a two-dimensional model that solves the flow equations. The turbulence and the turbulent diffusion are described by a large-eddy-simulation (LES) method where the large turbulence features are simulated by solving the flow equations, and a subgrid model represents the small-scale turbulence that is not resolved by the flow model. A dynamic eddy viscosity subgrid scale stress model has been used for the     

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.     

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.     

Bottom friction and its effects on periodic long wave propagation

A new set of Boussinesq-type equations describing the free surface evolution and the corresponding depth-integrated horizontal velocity is derived with the bottom boundary layer effects included. Inside the boundary layer the eddy viscosity gradient model is employed to characterize Reynolds stresses and the eddy viscosity is further approximated as a linear function of the distance measured from the seafloor. Boundary-layer velocities are coupled with the irrotational velocity in the core region through boundary conditions. The leading order boundary layer effects on wave propagation appear in the depth-integrated continuity equation to account for the velocity deficit inside the boundary layer. This formulation is different from the conventional approach in which a bottom stress term is inserted in the momentum equation. An iterative scheme is developed to solve the new model equations for the free surface elevation, depth-integrated velocity, the bottom stress, the boundary layer thickness and the magnitude of the turbulent eddy viscosity. A numerical example for the evolution of periodic waves propagating in one-dimensional channel is discussed to illustrate the numerical procedure and physics involved. The differences between the conventional approach and the present formulation are discussed in terms of the bottom frictional stress and the free surface profiles.     

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.     

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).     

Downward transfer of momentum by wind-driven waves

A model for the downward transfer of wind momentum is derived for growing waves. It is shown that waves, which grow due to an uneven pressure distribution on the water surface or a wave-coherent surface shear stress have horizontal velocities out of phase with the surface elevation. Further, if the waves grow in the x-direction, while the motion is perhaps time-periodic at any fixed point, the Reynolds stresses associated with the organized motion are positive. This is in agreement with several field and laboratory measurements which were previously unexplained, and the new theory successfully links measured wave growth rates and measured sub-surface Reynolds stresses. Wave coherent air pressure (and/or surface shear stress) is shown to change the speed of wave propagation as well as inducing growth or decay. From air pressure variations that are in phase with the surface elevation, the influence on the waves is simply a phase speed increase. For pressure variations out of phase with surface elevation, both growth (or decay) and phase speed changes occur. The theory is initially developed for long waves, after which the velocity potential and dispersion relation for linear waves in arbitrary depth are given. The model enables a sounder model for the transfer to storm surges or currents of momentum from breaking waves in that it does not rely entirely on ad-hoc turbulent diffusion. Future models of atmosphere-ocean exchanges should also acknowledge that momentum is transferred partly by the organized wave motion, while other species, like heat and gasses, may rely totally on turbulent diffusion. The fact that growing wind waves do in fact not generally obey the dispersion relation for free waves may need to be considered in future wind wave development models.     

Efficient computation of surf zone waves using the nonlinear shallow water equations with non-hydrostatic pressure

A numerical method for non-hydrostatic, free-surface, irrotational flow governed by the nonlinear shallow water equations including the effects of vertical acceleration is presented at the aim of studying surf zone phenomena. A vertical boundary-fitted grid is used with the water depth divided into a number of layers. A compact finite difference scheme is employed for accurate computation of frequency dispersion requiring a limited vertical resolution and hence, capable of predicting the onset of wave breaking. A novel wet–dry algorithm is applied for a proper handling of moving shoreline. Mass and momentum are strictly conserved at discrete level while the method only dissipates energy in the case of wave breaking. The numerical results are verified with a number of tests and show that the proposed model using two layers without ad-hoc assumptions enables to resolve propagating nonlinear shoaling, breaking waves and wave run-up within the surf and swash zones in an efficient manner.     

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.     

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.     

Large eddy simulation of breaking waves

A numerical model is used to simulate wave breaking, the large scale water motions and turbulence induced by the breaking process. The model consists of a free surface model using the surface markers method combined with a three-dimensional model that solves the flow equations. The turbulence is described by large eddy simulation where the larger turbulent features are simulated by solving the flow equations, and the small scale turbulence that is not resolved by the flow model is represented by a sub-grid model. A simple Smagorinsky sub-grid model has been used for the present simulations. The incoming waves are specified by a flux boundary condition. The waves are approaching in the shore-normal direction and are breaking on a plane, constant slope beach. The first few wave periods are simulated by a two-dimensional model in the vertical plane normal to the beach line. The model describes the steepening and the overturning of the wave. At a given instant, the model domain is extended to three dimensions, and the two-dimensional flow field develops spontaneously three-dimensional flow features with turbulent eddies. After a few wave periods, stationary (periodic) conditions are achieved. The surface is still specified to be uniform in the transverse (alongshore) direction, and it is only the flow field that is three-dimensional.The turbulent structures are investigated under different breaker types, spilling, weak plungers and strong plungers. The model is able to reproduce complicated flow phenomena such as obliquely descending eddies. The turbulent kinetic energy is found by averaging over the transverse direction. In spilling breakers, the turbulence is generated in a series of eddies in the shear layer under the surface roller. After the passage of the roller the turbulence spreads downwards. In the strong plunging breaker, the turbulence originates to a large degree from the topologically generated vorticity. The turbulence generated at the plunge point is almost immediately distributed over the entire water depth by large organised vortices. Away from the bed, the length scale of the turbulence (the characteristic size of the eddies resolved by the model) is similar in the horizontal and the vertical direction. It is found to be of the order one half of the water depth.     

Numerical simulation of breaking waves by a multi-scale turbulence model

In this paper, a two-dimensional multi-scale turbulence model is proposed to study breaking waves. The purpose of developing this model is to produce a relatively accurate model with moderate computer requirements. The free surface is tracked by the VOF technique, the log-law profile for the mean velocity is applied at the bottom. Comparing with the Reynolds-Averaged Navier-Stokes models (RANS), the present model shows improving agreement with experimental measurements in terms of surface elevations, particle velocities, wave height distributions and undertow profiles. The subgrid scale (SGS) turbulent transport mechanism is also discussed in the paper. It is found that turbulent production and dissipation are of the same order, but turbulent production is primarily located at the wavefront and above the wave trough, whereas turbulent dissipation is primarily located at the back face of a wave, indicating that in these regions, the assumption of equilibrium is not correct. Below the trough level, the local equilibrium assumption is reasonable. Turbulent convection and diffusion are of the same order at the trough level. Above the trough level, turbulent convection dominates. Under the spilling breaking wave, turbulent kinetic energy is continue to dissipate in the bore region, whereas under the plunging breaking wave, the turbulent kinetic energy is dissipated very rapidly within one wave period.     

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.     

Numerical Simulation of Spatial Lag Between Wave Breaking Point and Location of Maximum Wave-Induced Current

A quasi three-dimensional numerical model of wave-driven coastal currents with the effects of surface rollers is developed for the study of the spatial lag between the location of the maximum wave-induced current and the wave breaking point.The governing equations are derived from Navier-Stokes equations and solved by the hybrid method combining the fractional step finite different method in the horizontal plane with a Galerkin finite element method in the vertical direction.The surface rollers effects are considered through incorporating the creation and evolution of the roller area into the free surface shear stress.An energy equation facilitates the computation process which transfers the wave breaking energy dissipation to the surface roller energy.The wave driver model is a phase-averaged wave model based on the wave action balance equation.Two sets of laboratory experiments producing breaking waves that generated longshore currents on a planar beach are used to evaluate the model's performance.The present wave-driven coastal current model with the roller effect in the surface shear stress term can produce satisfactory results by increasing the wave-induced nearshore current velocity inside the surf zone and shifting the location of the maximum longshore current velocity landward.     

An integral contravariant formulation of the fully non-linear Boussinesq equations

In this paper we propose an integral form of the fully non-linear Boussinesq equations in contravariant formulation, in which Christoffel symbols are avoided, in order to simulate wave transformation phenomena, wave breaking and nearshore currents in computational domains representing the complex morphology of real coastal regions. Following the approach proposed by Chen (2006), the motion equations retain the term related to the approximation to the second order of the vertical vorticity. A new Upwind Weighted Essentially Non-Oscillatory scheme for the solution of the fully non-linear Boussinesq equations on generalised curvilinear coordinate systems is proposed. The equations are rearranged in order to solve them by a high resolution hybrid finite volume–finite difference scheme. The conservative part of the above-mentioned equations, consisting of the convective terms and the terms related to the free surface elevation, is discretised by a high-order shock-capturing finite volume scheme in which an exact Riemann solver is involved; dispersive terms and the term related to the approximation to the second order of the vertical vorticity are discretised by a cell-centred finite difference scheme. The shock-capturing method makes it possible to intrinsically model the wave breaking, therefore no additional terms are needed to take into account the breaking related energy dissipation in the surf zone. The model is verified against several benchmark tests, and the results are compared with experimental, theoretical and alternative numerical solutions.     

One-dimensional modelling of individual waves and wave-induced longshore currents in the surf zone

A probabilistic model ( -model) was developed to describe the propagation and transformation of individual waves (wave by wave approach). The individual waves shoal until an empirical criterion for breaking is satisfied. Wave height decay after breaking is modelled by using an energy dissipation method. Wave-induced set-up and set-down and breaking-associated longshore currents are also modelled. Laboratory and field data were used to calibrate and verify the model. The model was calibrated by adjusting the wave breaking coefficient (as a function of local wave steepness and bottom slope) to obtain optimum agreement between measured and computed wave height. Four tests carried out in the large Delta flume of Delft Hydraulics were considered. Generally, the measured H_1/3-wave heights are reasonably well represented by the model in all zones from deep water to the shallow surf zone. The fraction of breaking waves was reasonably well represented by the model in the upsloping zones of the bottom profile. Verification of the model results with respect to wave-induced longshore current velocities was not extensive, because of a lack of data. In case of a barred profile the measured longshore velocities showed a relatively uniform distribution in the (trough) zone between the bar crest and the shoreline, which could to some extent be modelled by including space-averaging of the radiation force gradient, horizontal mixing and longshore water surface gradients related to variations in set-up. In case of a monotonically upsloping profile the cross-shore distribution of the longshore current velocities is reasonably well represented.     

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.     

Application of a Boussinesq model for the computation of breaking waves: Part 2: Wave-induced setdown and setup on a submerged coral reef

Based on a set of Boussinesq-type equations with improved linear frequency dispersion characteristics in deeper water, the present paper incorporates the simplified effect of spilling wave breaking into the equations. The analysis is restricted to a single horizontal dimension but the method can be extended to include the second horizontal dimension. Inside the surf zone the vertical variation of the horizontal velocity profile is assumed to be composed of an (initially unknown) organised velocity component below the roller and a surface roller travelling with the wave celerity. This leads to a new set of equations which is capable of simulating the transformation of waves before, during and after wave breaking. The model is calibrated and verified by comparison with several wave flume measurements. The results show that the model produces sound physical results.