Boussinesq modelling of solitary wave and N-wave runup on coast

A total variation diminishing Lax–Wendroff scheme has been applied to numerically solve the Boussinesq-type equations. The runup processes on a vertical wall and on a uniform slope by various waves, including solitary waves, leading-depression N-waves and leading-elevation N-waves, have been investigated using the developed numerical model. The results agree well with the runup laws derived analytically by other researchers for non-breaking waves. The predictions with respect to breaking solitary waves generally follow the empirical runup relationship established from laboratory experiments, although some degree of over-prediction on the runup heights has been manifested. Such an over-prediction can be attributed to the exaggeration of the short waves in the front of the breaking waves. The study revealed that the leading-depression N-wave produced a higher runup than the solitary wave of the same amplitude, whereas the leading-elevation N-wave produced a slightly lower runup than the solitary wave of the same amplitude. For the runup on a vertical wall, this trend becomes prominent when the wave height-to-depth ratio exceeds 0.01. For the runup on a slope, this trend is prominent before the strong wave breaking occurs.     

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.     

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     

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

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

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.     

Large eddy simulation of spilling and plunging breakers

A Navier–Stokes solver with a free surface model is used for simulating wave breaking, undertow, and turbulence in breaking waves. The free surface model is based on the Volume of Fluid concept. Turbulence scales larger than the grid scale are simulated directly while turbulence scales smaller than the grid scale are represented by a sub-grid scale model. Two different approaches for the sub-grid scale model have been applied, which are the Smagorinsky model and a model based on a k-equation for the sub-grid scale turbulence. The waves approach the shore in shore-normal direction and break on a plane constant sloping beach. Periodic spilling and plunging breakers are simulated for 20 and 16 wave periods, respectively. The set-up, undertow, and turbulence levels are compared to experimental results. Despite the rather coarse resolution of the computational domain, satisfactory results for the wave height decay and undertow have been obtained. However, the turbulence levels are over-predicted when using the standard values of the model parameters and a complete answer to this problem has not been found. Furthermore, the evolution of vorticity over the wave period has been studied. It shows that at the initial breaking point vorticity is generated around the vertical as well as around the transverse axis. Later vorticity around the longitudinal axis (offshore–onshore direction) is generated, probably through deformation of vorticity around the other axis.     

Wave characteristic and morphologic effects on the onshore hydrodynamic response of tsunamis

While the destruction caused by a tsunami can vary significantly owing to near- and onshore controls, we have only a limited quantitative understanding of how different local parameters influence the onshore response of tsunamis. Here, a numerical model based on the non-linear shallow water equations is first shown to agree well with analytical expressions developed for periodic long waves inundating over planar slopes. More than 13,000 simulations are then conducted to examine the effects variations in the wave characteristics, bed slopes, and bottom roughness have on maximum tsunami run-up and water velocity at the still water shoreline. While deviations from periodic waves and planar slopes affect the onshore dynamics, the details of these effects depend on a combination of factors. In general, the effects differ for breaking and non-breaking waves, and are related to the relative shift of the waves along the breaking–non-breaking wave continuum. Variations that shift waves toward increased breaking, such as steeper wave fronts, tend to increase the onshore impact of non-breaking waves, but decrease the impact of already breaking waves. The onshore impact of a tsunami composed of multiple waves can be different from that of a single wave tsunami, with the largest difference occurring on long, shallow onshore topographies. These results demonstrate that the onshore response of a tsunami is complex, and that using analytical expressions derived from simplified conditions may not always be appropriate.     

Experimental investigation of breaking criteria of deepwater wind waves under strong wind action

The present study describes an experimental investigation of breaking criteria of deepwater wind waves under strong wind action. In a wind wave flume, waves were generated using different wind speeds and measured at different locations to obtain wave trains of no, intermittent, or frequent breaking. Water particle movement and free surface elevation were measured simultaneously using a PIV system and a wave gauge, respectively. For wind waves, not all the waves measured at a fixed location are breaking waves, and the breaking of a larger wave is not guaranteed. However, the larger the wave height, the larger the probability of breaking. In order to take as many breaking waves as possible for the cases of frequent breaking, we used the waves whose heights were close to the highest one-tenth wave height. The experimental results showed that the geometric or kinematic breaking criteria could not explain the occurrence of breaking of wind waves. On the other hand, the vertical acceleration beneath the wave crest was close to the previously suggested limit value, −0.5g, when frequent breaking of large waves occurred, indicating that the dynamic breaking criterion would be good for discriminating breaking waves under a strong wind action.     

Breaking wave prediction with boundary elements and finite volumes for use with small boat capsize studies: Convergence and resource requirements

The prediction of extreme breaking waves forms the foundation of many fields of research. The authors have recently completed a study in the capsize and re-righting of sailing yachts using breaking wave prediction to enhance experimental results. As breaking wave prediction is only the beginning of any research program a prediction method is required to be both accurate and computationally inexpensive. This paper describes the investigation of two methods varying in computational demand. It has been concluded that a non-linear free surface boundary element method is immediately realisable with application to a research program requiring a large number of predicted waves. A finite volume approach is realisable, but its engineering application across numerous waves is difficult.     

Experimental study of the breaking limit of multi-directional random waves passing over an impermeable submerged breakwater

This study investigates experimentally the breaking wave height of multi-directional random waves passing over an impermeable submerged breakwater. Experiments have been conducted in a three-dimensional wave basin equipped with a multi-directional random wave generator. A special type of wave gauge has been newly devised to record the water surface elevations in the breaker zone as accurately as possible. The records are analyzed to estimate the location and limit of wave breaking. Comparisons have also been made with the results of regular waves. The influence of the incident wave conditions on the breaking wave height normalized by the breakwater dimensions has been investigated. Empirical formulae have been presented to estimate the breaking limit of multi-directional random waves based on the experimental records. The formulae have been tested and found to work well not only for multi-directional random waves, but for regular waves as well.     

Surface tension effects on energy dissipation by small scale,experimental breaking waves

Much of the existing knowledge about breaking waves comes from physical model experiments scaled using Froude's law. A widely held assumption is that surface tension effects are not significant at typical laboratory scales and specifically for waves longer than 2 m. Since, however, smaller wavelengths are not untypical in small to medium scale laboratory facilities, a consideration of surface tension effects is indeed important. Although some emphasis has been given in the past, little is known regarding the importance of surface tension following impingement of the breaking-wave crest and especially on the overall energy dissipation by laboratory breaking waves.     

An integrated three-dimensional model of wave-induced pore pressure and effective stresses in a porous seabed: II. Breaking waves

The evaluation of the wave-induced liquefaction potential is particularly important for coastal engineers involved in the design of marine structures. Most previous investigations of the wave-induced liquefaction have been limited to two-dimensional non-breaking waves. In this paper, the integrated three-dimensional poro-elastic model for the wave-seabed interaction proposed by [Zhang, H., Jeng, D.-S., 2005. An integrated three-dimensional model of wave-induced pore pressure and effective stresses in a porous seabed: I. A sloping seabed. Ocean Engineering 32(5/6), 701–729.] is further extended to simulate the seabed liquefaction potential with breaking wave loading. Based on the parametric study, we conclude: (1) the liquefaction depth due to breaking waves is smaller than that of due to non-breaking waves; (2) the degree of saturation significantly affects the wave-induced liquefaction depth, and no liquefaction occurs in full saturated seabed, and (3) soil permeability does not only significantly affect the pore pressure, but also the shear stresses distribution.     

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.     

On the application of a Boussinesq model to river flows including shocks

SERR-1D is a 4th-order finite volume 1D Boussinesq model including wave breaking energy dissipation through extra diffusive-like terms. This model has been primarily conceived to compute wave propagation in coastal areas and has been validated for breaking and non-breaking waves propagating over uneven bathymetries (Cienfuegos et al., 2005, 2006a, b, 2007). The present paper aims at investigating the ability of SERR-1D to simulate challenging fluvial hydraulic applications such as sudden gate operation in open channels generating short waves, dam-break flows and a steady hydraulic jump over a bump. The performance of the absorbing-generating boundary condition implemented in SERR-1D is first analysed in the context of fluvial applications where relatively short waves must be evacuated from the computational domain without producing spurious reflection. Next, by comparing numerical results to analytical and experimental dam-break test cases we show that the model is able to reproduce the overall features of these flows, but that additional care should be paid to the representation of energy dissipation and front speed in order to accurately represent bore dynamics.     

结合椭圆型缓坡方程模拟近岸波流场

波浪向近岸传播的过程中,由波浪破碎效应所产生的近岸波流场是近岸海域关键的水动力学因素之一.结合近岸波浪场的椭圆型缓坡方程和近岸波流场数学模型对近岸波浪场及由斜向入射波浪破碎后所形成的近岸波流场进行了数值模拟.计算中考虑到波浪向近岸传播中由于波浪的折射、绕射、反射等效应使局部复杂区域波向不易确定,采用结合椭圆型缓坡方程所给出的波浪辐射应力公式来计算波浪产生的辐射应力,在此基础上耦合椭圆型缓坡方程和近岸波流场数学模型对近岸波流场进行数值模拟,从而使模型综合考虑了波浪的折射、绕射、反射等效应且避免了对波向角的直接求解,可以应用于相对较复杂区域的近岸波流场模拟.     

Application of VOF in Calculating Wave Energy Loss due to Spilling Breaker

- A large amount of experimental analysis and systematical theoretical calculation has been done by the authors to solve the problem of wave transformation and breaking, considering the effect of both current and topography, but only the wave energy loss due to spilling breaker in the surf zone has been discussed in this paper. Based on test result analysis and calculation with the Stream Function Wave Theory, the wave velocity field at breaking points has been obtained, and it is used to calculate the wave heights after breaking by the VOF (Volume of Fluid) method, in which the governing equations are continuity equation and Navier-Stokes Equation for imcompressible fluids. In the present paper, the improved VOF technique is used to calculate the wave heights of stable regular waves after breaking. Results fit the test data well, which shows that the VOF method is suitable to numerical simulation of regular waves after breaking. Besides, the breaker coefficient B of regular waves in the bore model is a     

Numerical simulation of wave overtopping of coastal structures using the non-linear shallow water equations

A one-dimensional high-resolution finite volume model capable of simulating storm waves propagating in the coastal surf zone and overtopping a sea wall is presented. The model (AMAZON) is based on solving the non-linear shallow water (NLSW) equations. A modern upwind scheme of the Godunov-type using an HLL approximate Riemann solver is described which captures bore waves in both transcritical and supercritical flows. By employing a finite volume formulation, the method can be implemented on an irregular, structured, boundary-fitted computational mesh. The use of the NLSW equations to model wave overtopping is computationally efficient and practically flexible, though the detailed structure of wave breaking is of course ignored. It is shown that wave overtopping at a vertical wall may also be approximately modelled by representing the wall as a steep bed slope. The AMAZON model solutions have been compared with analytical solutions and laboratory data for wave overtopping at sloping and vertical seawalls and good agreement has been found. The model requires more verification tests for irregular waves before its application as a generic design tool.     

波群内单个波的波陡分布与波破碎

对波群内单个波的波陡分布和波破碎进行了实验研究。研究结果是,波群中波动的最大振幅出现在波群前部而不是出现在波群中央,这种不对称性导致波群前部单个波出现大波陡的概率大于后部单个波出现大波陡的概率;进一步的波破碎统计发现波群前部单个波破碎的频率是后部单个波破碎频率的4倍。因此认为,波群结构的不对称性能够导致单个波发生破碎的...     

Mathematical analysis of long-wave breaking on open channels with bottom friction

In this paper we study the breaking of long waves propagating along an open channel with linear friction on the bottom. The equations governing the wave propagation consist of a pair of first-order nonlinear hyperbolic partial differential equations (PDEs). We first transformed the PDEs into a pair of ordinary differential equations (ODEs) along the characteristic directions by means of a pair of Riemann invariants. By analyzing the ODEs, we found that the breaking of waves can be identified by the singularity of the derivative of the Riemann invariants. Thus, we derived an analytical solution for the derivative of the Riemann invariants. Then, a breaking criterion and an analytical formula for the estimation of breaking time were developed and validated through numerical experiments. It is also shown in the paper that the present model includes the previous model neglecting bottom friction as a special case.     

白冠破碎空间分布模式的统计推断

白冠空间分布模式分析是海浪破碎统计研究的前提。本文提出利用空间点过程统计分析工具研究白冠空间分布模式,并结合实际白冠破碎观测录像资料,计算观测数据的L-函数和K-函数,与Markov chain Monte Carlo(MCMC)方法生成的模拟包迹进行比较,推断得出其白冠空间分布模式类型为空间齐次Poisson过程。研究表明空间点过程统计分析工具适用于白冠破碎研究。