基于激波捕捉类Boussinesq模型对岛礁地形上波浪传播的数值模拟

为了探究激波捕捉类Boussinesq模型在模拟岛礁地形上规则波和不规则波传播的可行性,采用基于完全非线性Boussinesq方程并具有激波捕捉能力的数值模型Funwave-TVD对规则波和不规则波在岛礁地形上的传播进行了数值模拟,通过与试验数据对比,分析模型中空间步长的影响,验证模型在模拟波高、平均水位分布以及波谱空间演变的能力,结果表明:采用合适的空间步长,模型能较好地模拟规则波和不规则波在岛礁地形上的传播和演化过程。对于规则波,较小的空间步长可改善破碎点处波高峰值的预测,并能更好地预测波浪破碎后波高的空间分布,相比结合经验破碎的Boussinesq模型,Funwave-TVD能更好地模拟规则波在岛礁地形上的破碎,以及破碎以后行进涌波的再生成过程;对于不规则波,Funwave-TVD总体而言能较好地模拟涌浪有效波高、次重力波的生成及空间演化和平均水位,但会低估礁坪上次重力波波高,较粗的空间步长也会低估礁坪上涌浪有效波高。     

Experimental Study of Wave Breaking on Gentle Slope

—An experimental study of regular wave and irregular wave breaking is performed on a gentleslope of 1:200.In the experiment,asymmetry of wave profile is analyzed to determine its effect on wavebreaker indices and to explain the difference between Goda and Nelson about the breaker indices of regu-lar waves on very mild slopes.The study shows that the breaker index of irregular waves is under less influ-ence of bottom slope i,relative water depth d/L_0 and the asymmetry of wave profile than that of regularwaves.The breaker index of regular waves from Goda may be used in the case of irregular waves, whilethe coefficient A should be 0.15.The ratio of irregular wavelength to the length calculated by linear wavetheory is 0.74.Analysis is also made on the waveheight damping coefficient of regular waves after break-ing and on the breaking probability of large irregular waves.     

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.     

Breaking Criteria for Regular Waves on Gentle Slopes

Nelson pointed out that the wave breaking criterion(H/d)_b for gentle slopes(i<1/100),inwhich H is the wave height and d is the water depth at the breaking point,is smaller than that for beachslopes of i>1/100),i.e.,the value of(H/d)_b for gentle slopes may be smaller than 0.6.Goda indicatedthat the wave breaking criterion given by himself is a result based on theoretical study,so it should be alsocorrect and can be used for gentle beaches,i.e.,the value of(H/d)_b for gentle slopes may be still largerthan 0.7.By use of high order nonlinear wave theory,this problem is analyzed in this study and the reasonwhy there is a large difference between different studies is explained.Moreover,the bottom energy loss dur-ing wave propagation is considered and the critical slope for the absence of wave breaking is also analyzed.     

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.     

Violent breaking wave impacts. Part 1: Results from large-scale regular wave tests on vertical and sloping walls

As part of an investigation into the detailed characteristics of wave impacts, experimental data are presented for the impact pressures and forces generated by waves up to 1.7 m high breaking onto a vertical wall and a wall inclined at 27° to the vertical. Particular attention is given to the influence of entrained and entrapped air and, by selecting regular wave conditions that produce impacts, trends are identified for highly variable phenomena that could easily be missed when masked by the even greater variability associated with irregular waves.     

Determination of wave energy dissipation factor and numerical simulation of wave height in the surf zone

Based on the time-dependent mild slope equation including the effect of wave energy dissipation, an expression for the energy dissipation factor is derived in conjunction with the wave energy balance equation. The wave height of regular and irregular waves is numerically simulated by use of the parabolic mild slope equation considering the energy dissipation due to wave breaking. Comparison of numerical results with experimental data shows that the expression for the energy dissipation factor is reasonable. The effects of the wave breaking coefficient on the breaking point and the distribution of wave height after breaking are discussed through the study of a specific experimental topography.     

Wave Forces Acting on Vertical Walls

Regular and irregular wave forces acting on vertical walls are studied by a previously developed numerical model. The computed wave forces are compared with the available experimental data to verify the numerical model, and satisfactory agreements are obtained. The variation of wave forces with incident angles and the shape of simultaneous pressure distribution are investigated, and the comparisons between numerical results and Goda＇ s predictions are also carried out. It is concluded that the maximum wave forces acting on the unit length of vertical wall is often induced by the obliquely incident waves instead of normally incident waves, while Goda＇ s formula may be inapplicable for oblique wave incidence. The shape of simultaneous pressure distribution is not significantly influenced by incident angles, and it can be favorably predicted by Goda＇ s formula. When regular wave heights are taken as the same as irregular wave height H1%, the irregular wave forces Ph. 1% are slightly larger than regular wave forces in most cases.     

Theoretical Model and Random Simulation for Nonlinear Interaction of Irregular Waves with Vertical Wall

According to the theoretical solutions for the nonlinear three-dimensional gravity surface waves and their interactions with vertical wall previously proposed by the lead author, in this paper an exact second-order random model of the unified wave motion process for nonlinear irregular waves and their interactions with vertical wall in uniform current is formulated, the corresponding theoretical nonlinear spectrum is derived, and the digital simulation model suitable to the use of the FFT (Fast Fourier Tansform) algorithm is also given. Simulations of wave surface, wave pressure, total wave pressure and its moment are performed. The probability properties and statistical characteristics of these realizations are tested, which include the verifications of normality for linear process and of non-normality for nonlinear process; the consistances of the theoretical spectra with simulated ones; the probability properties of apparent characterstics, such as amplitudes, periods, and extremes (maximum and minimum, positive and negative extremes). The statistical analysis and comparisons demonstrate that the proposed theoretical and computing models are realistic and effective, the estimated spectra are in good agreement with the theoretical ones, and the probability properties of the simulated waves are similar to those of the sea waves. At the same time, the simulating computation can be completed rapidly and easily.     

Laboratory study of the nonlinear transformation of irregular waves over a mild slope

This paper considers the nonlinear transformation of irregular waves propagating over a mild slope （1：40）. Two cases of irregular waves, which are mechanically generated based on JONSWAP spectra, are used for this purpose. The results indicate that the wave heights obey the Rayleigh distribution at the offshore location; however, in the shoaling region, the heights of the largest waves are underestimated by the theoretical distributions. In the surf zone, the wave heights can be approximated by the composite Weibull distribution. In addition, the nonlinear phase coupling within the irregular waves is investigated by the wavelet-based bicoherence. The bicoherence spectra reflect that the number of frequency modes participating in the phase coupling increases with the decreasing water depth, as does the degree of phase coupling. After the incipient breaking, even though the degree of phase coupling decreases, a great number of higher harmonic wave modes are also involved in nonlinear interactions. Moreover, the summed bicoherence indicates that the frequency mode related to the strongest local nonlinear interactions shifts to higher harmonics with the decreasing water depth.     

Laboratory Study of the Nonlinear Transformation of Irregular Waves over A Mild Slope

This paper considers the nonlinear transformation of irregular waves propagating over a mild slope (1?40). Two cases of irregular waves, which are mechanically generated based on JONSWAP spectra, are used for this purpose. The results indicate that the wave heights obey the Rayleigh distribution at the offshore location; however, in the shoaling region, the heights of the largest waves are underestimated by the theoretical distributions. In the surf zone, the wave heights can be approximated by the composite Weibull distribution. In addition, the nonlinear phase coupling within the irregular waves is investigated by the wavelet-based bicoherence. The bicoherence spectra reflect that the number of frequency modes participating in the phase coupling increases with the decreasing water depth, as does the degree of phase coupling. After the incipient breaking, even though the degree of phase coupling decreases, a great number of higher harmonic wave modes are also involved in nonlinear interactions. Moreover, the summed bicoherence indicates that the frequency mode related to the strongest local nonlinear interactions shifts to higher harmonics with the decreasing water depth.     

海岸沙坝剖面和滩肩剖面特征研究

通过物理模型实验研究了海岸沙坝剖面和滩肩剖面的形成和演化过程,给出了稳定的沙坝剖面和滩肩剖面的几何特征。实验中考虑了两种初始坡度(1∶20和1∶10)和不同波高的规则波和不规则波,讨论了不同初始坡度海岸上破碎波空间分布特征。结果表明,沙坝产生后存在向岸和离岸两种运动形态,但最终将停留在稳定位置。稳定的沙坝剖面对应不同初始坡度和波浪存在不同的大沙坝和小沙坝分布。沙坝剖面由长时间小波高波浪序列作用后可转化为稳定滩肩剖面,该剖面不依赖于波浪和初始坡度。实验也给出了平衡剖面与理论曲线的对比以及剖面上泥沙粒径的分布。     

波浪作用下沙坝不稳定性实验研究

对波群、规则波和不规则波作用下沙质海岸上沙坝运动进行了实验研究。研究沙坝产生的机理和不稳定性,以及不同波浪下所产生的海岸平衡剖面形状。实验中发现,波浪作用下沙坝处于向岸和离岸运动的不稳定状态,其原因可由地形不稳定性分析。探讨了波群调制系数对沙坝运动的影响,讨论了不同波况对海岸剖面形态的影响。     

Second-Oorder Analytic Solutions of Nonlinear Interactions of Edge Waves on A Plane Sloping Bottom

Based on the full water-wave equation, a second-order analytic solution for nonlinear interaction of short edge waves on a constant plane sloping bottom is presented in this paper. For special case of slope angle b=p/2, this solution can be reduced to the same order solution of deep water gravity surface waves traveling along parallel coastline. Interactions between two edge waves including progressive, standing and partially reflected standing waves were also discussed. The unified analytic expressions with transfer functions for kinematic-dynamic elements of edge waves were also discussed. The random model of the unified wave motion processes for linear and nonlinear irregular edge waves is formulated, and the corresponding theoretical autocorrelation and spectral density functions of the first and second orders are derived. The boundary conditions for the determining determination of the parameters of short edge wave are suggested, that may be seen as one special simple edge wave excitation mechanism and an extension to the sea wave refraction theory. Finally some computation results are demonstrated.     

多向不规则波浪的确定性模拟

波浪波动时间过程及波列的模拟,对于开展实际波浪对于工程建筑物的作用具有重要的意义。本文采用线性叠加的单叠加模型,建立了多向不规则波浪的确定性模拟方法。基于理论模拟的规则波、单向不规则波和多向不规则波,验证了波浪确定性模拟方法的有效性。定性地对比分析了模拟波列和已知波列的一致性;定量地研究了模拟波浪在空间范围rr/L_s的误差分布情况(rr表示指定位置与给定位置的空间距离,L_s为有效波长)。并且建议,采用本文方法进行波浪确定性模拟时,最佳的浪高仪间距应小于0.12L_s。     

Transformation and Breaking of Irregular Waves on Very Gentle Slopes

—Based on theoretical analysis.numerical calculation.and experimental study,this paper dis-cusses breaker indices of irregular waves.transformation of wave spectrum.characteristics and computa-tion of breaking waves.as well as the critical beach slope under which waves will not break.Computed re-sults are in good agreement with laboratory physical model test data and ocean wave field measurements.     

Theoretical and Experimental Study of Breaking Wave on Sloping Bottoms

This paper studies the continuous evolution of breaking wave for the surface water waves propagating on a sloping beach. A Lagrangian asymptotic solution is derived. According to the solution coupled with the wave breaking criteria and the equations of water particles motion, the wave deformation and the continuous wave breaking processes for the progressive water waves propagating on a sloping bottom can be derived. A series of experiments are also conducted to compare with the theoretical solution. The results show that the present solution can reasonably describe the plunging or spilling wave breaking phenomenon.     

Wave-Current Forces on Small Square Cylinders (Part II)

Based on model tests, the lift and resultant forces on small square cylinders caused by waves (regular and irregular) and currents are analyzed in this paper. The lift and resultant force coefficients CL and Cf related to KC number and the effect of direction of wave propagation are also given, which may be useful for practical engineering application.     

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     

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.