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.     

规则波在浅水中运动及破碎的数值计算

利用一维Boussinesq方程描述了在浅水中的波浪运动以及破碎情况。在方程中引入了表面翻滚的概念，认为翻滚的水体是以波速运动的，翻滚的作用表现在水平速度的垂直分布上，产生了附加迁移动量项。通过对Airy波在浅水中运动以及破碎情况的研究，得出的一些结论与Schaeffer和合田的成果吻合良好。     

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

Characteristics of breaking irregular wave forces on a monopile

The substructures of offshore wind turbines are subjected to extreme breaking irregular wave forces. The present study is focused on investigating breaking irregular wave forces on a monopile using a computational fluid dynamics (CFD) based numerical model. The breaking irregular wave forces on a monopile mounted on a slope are investigated with a numerical wave tank. The experimental and numerical irregular free surface elevations are compared in the frequency-domain for the different locations in the vicinity of the cylinder. A numerical analysis is performed for different wave steepness cases to understand the influence of wave steepness on the breaking irregular wave loads. The wave height transformation and energy level evolution during the wave shoaling and wave breaking processes is investigated. The higher-frequency components generated during the wave breaking process are observed to play a significant role in initiating the secondary force peaks. The free surface elevation skewness and spectral bandwidth during the wave transformation process are analysed and an investigation is performed to establish a correlation of these parameters with the breaking irregular wave forces. The role of the horizontal wave-induced water particle velocity at the free surface and free surface pressure in determining the breaking wave loads is highlighted. The higher-frequency components in the velocity and pressure spectrum are observed to be significant in influencing the secondary peaks in the breaking wave force spectrum.     

Large-wave simulation of three-dimensional,cross-shore and oblique,spilling breaking on constant slope beach

In this work, the large-wave simulation (LWS) method is adapted for application in spilling wave breaking over a constant slope beach. According to LWS, large scales of velocities, pressure and free-surface elevation are numerically resolved, while the corresponding unresolved scale effects are taken into consideration by a subgrid scale (SGS) model for wave and eddy stresses. The model may be not fully applicable in very shallow water, close to the shoreline, where the unresolved, turbulent, free-surface oscillation is of the same order with the water depth. Time integration of the Euler equations is achieved by a two-stage fractional scheme, combined with a hybrid scheme for spatial discretization, consisting of finite difference and pseudospectral approximation methods. Model parameters are calibrated by comparison to available experimental data of free-surface elevation and velocities in the surf zone for cross-shore incoming waves. The action of the wave SGS stresses in the outer coastal and surf zones initiates breaking and generates appropriate vorticity, in the form of an eddy structure (surface roller), at the breaking wavefront. At incipient breaking, both advection and gravity contribute to the vorticity flux at the free surface, while only after the full development of the surface roller, the effect of advection becomes stronger. The SGS model is also utilized to simulate propagation, refraction and breaking of oblique incoming waves. The gradual breaking and dissipation of wave crestlines and the surface roller structure along the breaking wavefront are automatically captured without any empirical input, such as data for the roller shape or the wave propagation angle at breaking.     

Boussinesq-type modelling using an unstructured finite element technique

A model for solving the two-dimensional enhanced Boussinesq equations is presented. The model equations are discretised in space using an unstructured finite element technique. The standard Galerkin method with mixed interpolation is applied. The time discretisation is performed using an explicit three-step Taylor–Galerkin method. The model is extended to the surf and swash zone by inclusion of wave breaking and a moving boundary at the shoreline. Breaking is treated by an existing surface roller model, but a new procedure for the detection of the roller thickness is devised. The model is verified using four test cases and the results are compared with experimental data and results from an existing finite difference Boussinesq model.     

A numerical model of wave overtopping on seadikes

This paper describes the development of a numerical model for wave overtopping on seadikes. The model is based on the flux-conservative form of the nonlinear shallow water equations (NLSW) solved with a high order total variation diminishing (TVD), Roe-type scheme. The goal is to reliably predict the hydrodynamics of wave overtopping on the dike crest and along the inner slope, necessary for the breach modelling of seadikes. Besides the mean overtopping rate, the capability of simulating individual overtopping events is also required. It is shown theoretically that the effect of wave breaking through the drastic motion of surface rollers in the surfzone is not sufficiently described by the conventional nonlinear shallow water equations, neglecting wave setup from the mean water level and thus markedly reducing the model predictive capacity for wave overtopping. This is significantly improved by including an additional source term associated with the roller energy dissipation in the depth-averaged momentum equation. The developed model has been validated against four existing laboratory datasets of wave overtopping on dikes. The first two sets are to validate the roller term performance in improving the model prediction of wave overtopping of breaking waves. The last two sets are to test the model performance under more complex but realistic hydraulic and slope geometric conditions. The results confirm the merit of the supplemented roller term and also demonstrate that the model is robust and reliable for the prediction of wave overtopping on seadikes.     

A function to determine wavelength from deep into shallow water based on the length of the cnoidal wave at breaking—Reply to discussion by T.S. Hedges

The equations of Hedges [Hedges, T.S., 2009. Discussion of “A function to determine wavelength from deep into shallow water based on the length of the cnoidal wave at breaking” by J.P. Le Roux, Coastal Eng.], although yielding similar wavelengths, are not consistent with the fact that the horizontal water particle velocity in the wave crest should equal the wave celerity at breaking over a nearly horizontal bottom.     

Numerical models for evolution of extreme wave groups

The paper considers the application of two numerical models to simulate the evolution of steep breaking waves. The first one is a Lagrangian wave model based on equations of motion of an inviscid fluid in Lagrangian coordinates. A method for treating spilling breaking is introduced and includes dissipative suppression of the breaker and correction of crest shape to improve the post breaking behaviour. The model is used to create a Lagrangian numerical wave tank, to reproduce experimental results of wave group evolution. The same set of experiments is modelled using a novel VoF numerical wave tank created using OpenFOAM. Lagrangian numerical results are validated against experiments and VoF computations and good agreement is demonstrated. Differences are observed only for a small region around the breaking crest.     

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.     

An extension of the Airy theory for linear waves into shallow water

As a fully developed (Airy) wave propagates from deep into shallow water, its crest becomes more peaked while the trough flattens out. The median crest diameter MCD, defined as the distance between the wave flanks under the crest at a level halfway between the crest and trough, therefore decreases relative to the similarly defined median trough diameter MTD, which remains constant up to the breaking point. The MCD is directly related to other wave characteristics, which enables water particle velocities to be calculated for any water depth without having to recur to more complex, higher-order Stokes, cnoidal or Fenton theories. Over a nearly horizontal bottom, most fully developed wave characteristics can be expressed as functions of the wave period T_w. It is shown that the horizontal particle velocity at the bottom under the breaker crest is at least 9 times faster than under the breaker trough, which explains why sediment is transported landward under fair weather conditions. The proposed equations also shed new light on the formation of spilling, plunging and surging/collapsing breakers.     

Implementation and evaluation of alternative wave breaking formulas in a coastal spectral wave model

This paper describes methods and results of research for incorporating four different parameterized wave breaking and dissipation formulas in a coastal wave prediction model. Two formulations assume the breaking energy dissipation to be limited by the Rayleigh distribution, whereas the other two represent the breaking wave energy by a bore model. These four formulations have been implemented in WABED, a directional spectral wave model based on the wave action balance equation with diffraction, reflection, and wave–current interaction capabilities. Four parameterized wave breaking formulations are evaluated in the present study using two high-quality laboratory data sets. The first data set is from a wave transformation experiment at an idealized inlet entrance, representing four incident irregular waves in a slack tide and two steady-state ebb current conditions. The second data set is from a laboratory study of wave propagation over a complex bathymetry with strong wave-induced currents. Numerical simulation results show that with a proper breaking formulation the wave model can reproduce laboratory data for waves propagating over idealized or complicated bathymetries with ambient currents. The extended Goda wave breaking formulation with a truncated Rayleigh distribution, and the Battjes and Janssen formulation with a bore model produced the best agreement between model and data.     

Establishment,validation and discussions of a three dimensional wave-induced current model

A three-dimensional numerical model was established to simulate the wave-induced currents. The depth-varying residual momentum, surface roller, wave horizontal and vertical turbulent mixing effects were incorporated as major driving forces. A surface roller evolution model considering the energy transfer, roller density and bottom slope dissipation was developed. The expression of the wave-induced horizontal turbulent mixing coefficient proposed by Larson and Kraus (1991) was extended to three-dimensional form. Plenty of experimental cases were used to validate the established model covering the wave setup, undertow, longshore currents and rip currents. Validation results showed the model could reasonably describe the main characteristics of different wave-induced current phenomena. The incorporation of surface roller for breaking waves should not be neglected in the modeling of surfzone hydrodynamics. The wave-induced turbulent mixing affects the structures of wave-induced current either in horizontal or in vertical directions. Sensitivity analysis of the major calibration parameters in the established model was made and their ranges were evaluated.     

Examination of empirical formulas for wave shoaling and breaking on steep slopes

Estimation of the wave height transformation of shoaling and breaking is essential for the nearshore hydrodynamics and the design of coastal structures. Many empirical formulas have been well recognized to the wave height transformation, but most of them were only applicable for gentle slopes. This paper reports the experimental results of wave shoaling and breaking over the steep slopes to examine the applicability of the previous empirical formulas. Two steep bottom slopes of 1/3 and 1/5, and one gentle slope of 1/10 were conducted in the present experiments. Experimental results show that the shoaling distance of steep slopes become short and the surface waves may be partially reflected from the steep bottom, thus the estimation of wave shoaling using the well-known previous formula did not conform completely to the experimental results. The previous empirical formulas for the wave breaking criteria were also examined, and the modified equations to the steep beaches were proposed in this work. A numerical model was finally adopted to calculate the wave height transformation in the surf zone by introducing the modified breaking index.     

High order Hamiltonian water wave models with wave-breaking mechanism

Based on the Hamiltonian formulation of water waves, using Hamiltonian consistent modelling methods, we derive higher order Hamiltonian equations by Taylor expansions of the potential and the vertical velocity around the still water level. The polynomial expansion in wave height is mixed with pseudo-differential operators that preserve the exact dispersion relation. The consistent approximate equations have inherited the Hamiltonian structure and give exact conservation of the approximate energy. In order to deal with breaking waves, we extend the eddy-viscosity model of Kennedy et al. (2000) to be applicable for fully dispersive equations. As breaking trigger mechanism we use a kinematic criterion based on the quotient of horizontal fluid velocity at the crest and the crest speed. The performance is illustrated by comparing simulations with experimental data for an irregular breaking wave with a peak period of 12 s above deep water and for a bathymetry induced periodic wave plunging breaker over a trapezoidal bar. The comparisons show that the higher order models perform quite well; the extension with the breaking wave mechanism improves the simulations significantly.     

The FUNWAVE model application and its validation using laboratory data

In this paper, the performance of a well known and widespread Boussinesq model is evaluated for its ability to predict deep waters and shoaling zone velocity and elevation, comparing them with laboratory data. The model evaluation places emphasis on the parameters of the wave field that could be used for the prediction of coastal phenomena. It was found that the model proved its ability to predict regular wave velocity and elevation both in deep waters and shoaling zones right up to the breaking region. In contrast to previous works found in literature, the comparisons in this case are not limited to modelled and measured wave elevations, but also involve horizontal velocities. As a result, experimental and computed time series of wave elevations and horizontal velocities and their corresponding phase-averaged values are investigated along the wave channel, showing a good performance of the model.     

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

海浪破碎对海洋上混合层中湍能量收支的影响

海浪破碎产生一向下输入的湍动能通量,在近海表处形成一湍流生成明显增加的次层,加强了海洋上混合层中的湍流垂向混合。为了研究海浪破碎对混合层中湍能量收支的影响,文中分析了海浪破碎对海洋上混合层中湍流生成的影响机制,采用垂向一维湍封闭混合模式,通过改变湍动能方程的上边界条件,引入了海浪破碎产生的湍动能通量,并分别对不同风速下海浪破碎的影响进行了数值研究,分析了混合层中湍能量收支的变化。当考虑海浪破碎影响时,近海表次层中的垂直扩散项和耗散项都有显著的增加,该次层中被耗散的湍动能占整个混合层中耗散的总的湍能量的92.0%,比无海浪破碎影响的结果增加了近1倍;由于平均流场切变减小,混合层中的湍流剪切生成减小了3.5%,形成一种存在于湍动能的耗散和垂直扩散之间的局部平衡关系。在该次层以下,局部平衡关系与壁层定律的结论一致,即湍动能的剪切生成与耗散相平衡。研究结果表明,海浪破碎在海表产生的湍动能通量影响了海洋上混合层中的各项湍能量收支间的局部平衡关系。     

Analysis of some key parametrizations in a beach profile morphodynamical model

A numerical process-based model to forecast beach profile morphodynamics has been developed. In the present paper, an analysis of various modelling approaches and key parametrizations involved in the estimation of the wave-driven current and the suspended sediment concentration is carried out.Several resolution techniques for the 1DV horizontal (i.e., in the x-direction perpendicular to coastline) momentum equation governing the Mean Horizontal Velocity (MHV) are analysed. In the first kind of techniques, the mean horizontal velocity is computed from the momentum equation, whereas the Mean Water Level (MWL) is computed using a parametrization of the depth-averaged momentum equation. Two boundary or integral conditions are thus needed. In the second kind, both mean horizontal velocity and mean water level gradient in the x-direction are the unknowns of the momentum equation, thus, three boundary or integral conditions are needed. Various additional conditions are discussed. We show that using a technique of the first kind is equivalent to imposing the difference between the surface and the bottom shear stresses in the 1D vertical equation. Both techniques lead to results that are in good agreement with the Delta Flume experimental data, provided the Stokes drift flow discharge is imposed as an additional condition. The influence of the breaking roller model and of the turbulent viscosity parametrization are also analysed.Suspended sediment transport by the mean current and wave-induced bedload transport are taken into account in the sediment flux. Three turbulent diffusivity parametrizations are compared for suspended sediment concentration estimations. A linear profile for the turbulent diffusivity taking into account the wave bottom shear stress and the surface wave breaking turbulence production is shown to give the best results. Using experimental data, we put forward the poor estimation of the bottom sediment concentration given by the three implemented parametrizations. We thus propose a new parametrization relying on a Shields parameter based on the breaking roller induced surface shear stress. Using this new parametrization, the bottom profile used in the tests keeps its two bars which disappear otherwise. However, the morphodynamical model still overestimates the bars offshore motion, a bias already observed in other models.