Stokes波的谱特性分析

应用理论推导及数值计算方法，对Stokes随机波的谱特性进行了分析。首先将波面方程，海水质点水平速度用一阶波面分量的非线性组合表示，应用平稳随机高阶短的降阶计算法则，得到了波面方程及海水质点水平速度与一阶波面分量的自相关函数之间的关系，从而确定了Stokes随机波浪的波浪谱密度及海水质点水平速度和加速度谱密度，进而求得有关波浪要素的均方根值。文章还应有数值计算方法，分析了波浪基本参数对均方根值的影响。     

三阶Stokes型随机波浪载荷谱研究

应用非线性谱分析理论，对三阶Stokes型随机波浪载荷谱进行了分析，将波面方程及海水质点的水平速度用一阶波面的非线性组表示，导出了随机波浪谱的表达式。为了便于求解随机波浪的载荷谱，将阻力项展开为幂级数式，并应用非线性谱分析理论，确定了幂级数的系数，进而将波浪载荷表示为一阶波面及其导数的非线性组合，最后得出波浪载荷谱密度的表达式，并应用数值分析方法，得出单位桩柱波浪力及总波浪力谱密度。     

Kinematics of extreme waves in deep water

The velocity profiles under crest of a total of 62 different steep wave events in deep water are measured in laboratory using particle image velocimetry. The waves take place in the leading unsteady part of a wave train, focusing wave fields and random wave series. Complementary fully nonlinear theoretical/numerical wave computations are performed. The experimental velocities have been put on a nondimensional form in the following way: from the wave record (at a fixed point) the (local) trough-to-trough period, T_TT and the maximal elevation above mean water level, η_m of an individual steep wave event are identified. The local wavenumber, k and an estimate of the wave slope, ε are evaluated from ω²/(gk)=1+ε², where ω=2π/T_TT and g denotes the acceleration of gravity. A reference fluid velocity, is then defined. Deep water waves with a fluid velocity up to 75% of the estimated wave speed are measured. The corresponding kη_m is 0.62. A strong collapse of the nondimensional experimental velocity profiles is found. This is also true with the fully nonlinear computations of transient waves. There is excellent agreement between the present measurements and previously published Laser Doppler Anemometry data. A surprising result, obtained by comparison, is that the nondimensional experimental velocities fit with the exponential profile, i.e. e^ky, y the vertical coordinate, with y=0 in the mean water level.     

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.     

Second-Order 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 plane sloping bottom is presented in this paper.For special case of slope angle β=π/2,this solution can 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 are also discussed.The unified analytic expressions with transfer functions for kinematic-dynamic elements of edge waves are also given.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 the second orders are derived.The boundary conditions for the 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.     

Derivation and application of new equations for radiation stress and volume flux

In this paper, new expressions of radiation stress and volume flux for long waves have been analytically derived by inclusion of higher-order surface elevations up to the sixth-order. To quantify these expressions, surface elevations along a beach are first simulated using the fully nonlinear Boussinesq-type model COULWAVE. Then, based on the large amount of numerical data, new equations for radiation stress and volume flux are statistically formulated. The research unveils the essential roles of the Ursell parameter, Irribarren number and wave steepness described by the local wave height, wave length and bottom slope. The study shows the importance of nonlinear wave properties in wave-induced currents and mean water levels (set-up/down). The higher-order formulations produce lower values for radiation stress and volume flux than calculated from the lower-order and linear waves. Case studies suggest that the new formulations produce an accurate estimation for mean water level. However, improvement on the computed current profiles is marginal for some cases. This implies that the accurate prediction of the current profile would require more than just the proposed improvement of the radiation stress and volume flux.     

基于二阶斯托克斯波理论的辐射应力垂向分布

基于二阶斯托克斯波理论推导了辐射应力的垂向分布表达式,通过算例讨论了辐射应力在深水和有限水深条件下的垂向分布规律,并与基于微幅波理论的辐射应力进行了比较.结果表明,在波浪非线性不强时,基于二阶斯托克斯波理论的辐射应力与基于微幅波理论的辐射应力表达式计算结果接近;而当水深较浅波浪非线性较强时,基于二阶斯托克斯波理论的辐射应力在近表面处明显大于基于微幅波理论的辐射应力.采用二阶斯托克斯波理论推导的波浪辐射应力更为合理地反映了波浪非线性效应.     

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.     

Validation of the three-wave quasi-kinetic approximation for the spectral evolution in shallow water

This paper aims at validating the three-wave quasi-kinetic approximation for the spectral evolution of weakly nonlinear gravity waves in shallow water. The problem is investigated using a one-dimensional numerical wave propagation model, formulated in the spectral representation. This model includes both a nonlinear triad interactions term and a wave breaking dissipation term. Some numerical tests were carried out in order to show the importance of using the triad nonlinear term in wave propagation spectral models, particularly to describe both behavior of the spectral integral parameters and of the spectral shape evolution in shallow water depth. Furthermore; a comparison against different set of experimental observations was carried out. Comparing the numerical results with the experimental observations made it possible to show the modeling efficiency of the three-wave quasi-kinetic approximation.     

Statistical characteristics of long waves nearshore

The random long wave runup on a beach of constant slope is studied in the framework of the rigorous solutions of the nonlinear shallow water theory. These solutions are used for calculation of the statistical characteristics of the vertical displacement of the moving shoreline and its horizontal velocity. It is shown that probability characteristics of the runup heights and extreme values of the shoreline velocity coincide in the linear and nonlinear theory. If the incident wave is represented by a narrow-band Gaussian process, the runup height is described by a Rayleigh distribution. The significant runup height can also be found within the linear theory of long wave shoaling and runup. Wave nonlinearity nearshore does not affect the Gaussian probability distribution of the velocity of the moving shoreline. However the vertical displacement of the moving shoreline becomes non-Gaussian due to the wave nonlinearity. Its statistical moments are calculated analytically. It is shown that the mean water level increases (setup), the skewness is always positive and kurtosis is positive for weak amplitude waves and negative for strongly nonlinear waves. The probability of the wave breaking is also calculated and conditions of validity of the analytical theory are discussed. The spectral and statistical characteristics of the moving shoreline are studied in detail. It is shown that the probability of coastal floods grows with an increase in the nonlinearity. Randomness of the wave field nearshore leads to an increase in the wave spectrum width.     

Numerical Modeling of Water Wave Interaction with A Soft Mud Bed

A vertical 2-D numerical model is presented for simulating the interaction between water waves and a soft mud bed. Taking into aceotmt nonlinear theology, a semi-empirical theological model is applied to this water-mud model, reflecting the combined visco-elasto-plastic properties of soft mud under such oscillatory external forces as water waves. In order to increase the resolution of the flow in the neighborhood of both sides of the inter-surface, a logarithmic grid in the vertical direction is employed for numerical treatment. Model verifications are given through comparisons between the calculated and the measured mud mass transport velocities as well as wave height changes.     

Intercomparison of wave-driven current models

Seven numerical models which simulate waves and currents in the surf-zone are tested for the case of a reduced-scale detached breakwater subjected to the action of regular waves with normal incidence. The computed wave heights, water levels and velocities are compared with measurements collected in an experimental wave basin. The wave height decay in the surf-zone is predicted reasonably well. Set-up and currents appear to be less well predicted. This intercomparison exercise shows that radiation stresses are systematically overestimated by formulations used in the models, mean bottom shear stresses are not always co-linear with the mean bottom velocity vector in shallow water, and turbulence modelling in the surf-zone requires a sophisticated     

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

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

An unsteady wave driver for narrowbanded waves: modeling nearshore circulation driven by wave groups

In this paper, we derive an unsteady refraction–diffraction model for narrowbanded water waves for use in computing coupled wave–current motion in the nearshore. The end result is a variable coefficient, nonlinear Schrödinger-type wave driver (describing the envelope of narrow-banded incident waves) coupled to forced nonlinear shallow water equations (describing steady or unsteady mean flows driven by the short-wave field). Comparisons with experimental data show that good accuracy can be obtained for cases of nonbreaking wave transformation. Numerical simulations show that the interaction of wave groups with longshore topographic nonuniformities generates strong edge wave resonances, providing a generating mechanism for low-order edge waves.     

Transparent boundary condition for simulating nonlinear water waves by a particle method

Open boundaries are important when simulating water waves. In this study, a transparent boundary condition at an open boundary was developed for simulating nonlinear water waves propagating to a distant area using the Moving Particle Semi-implicit method. The novelty of this study is that the technique of wave analysis used in the experiment was introduced into the particle simulation to absorb incident waves; the simulation cost was reduced by employing inflow and outflow regions instead of a long dissipation region. Incident waves in front of the boundary were evaluated using Fourier analysis, and the particles on the transparent boundary were forced to move at the velocity of the analytical solution for Stokes waves in order to absorb the incident waves. The analysis was restricted to periodic waves. Wave propagation was simulated for two wave periods using the developed transparent boundary condition. The results showed that this transparent boundary transmitted the incident waves with small reflection and the simulation cost was lower than that for wave damping by a conventional highly viscous region.     

Ephemeral sand waves in the hurricane surf zone

Airborne bathymetric LIDAR observations along the Florida panhandle after Hurricane Dennis (2005) show the first unequivocal observations of surf-zone sand wave trains.

These are found in depths of 5m along the trough of the hurricane bar, where hindcasts show strong longshore currents only during severe storms. The waves extend over tens of kilometers of coast after Dennis but are absent from the same area in four other datasets. Observed wavelength to water depth ratios are comparable to river dunes and tidal sand waves but height to depth ratios are smaller, with the largest wave heights around 0.1 times the water depth. The sand wave generation mechanism is hypothesized to be from wind-and-wave-induced longshore currents, which were hindcast to be large during Dennis, with destruction from water wave orbital velocities.     

Comparison of linear and nonlinear extreme wave statistics

An extremely large (“freak”) wave is a typical though rare phenomenon observed in the sea. Special theories (for example, the modulation instability theory) were developed to explain mechanics and appearance of freak waves as a result of nonlinear wave-wave interactions. In this paper, it is demonstrated that the freak wave appearance can be also explained by superposition of linear modes with the realistic spectrum. The integral probability of trough-to-crest waves is calculated by two methods: the first one is based on the results of the numerical simulation of a wave field evolution performed with one-dimensional and two-dimensional nonlinear models. The second method is based on calculation of the same probability over the ensembles of wave fields constructed as a superposition of linear waves with random phases and the spectrum similar to that used in the nonlinear simulations. It is shown that the integral probabilities for nonlinear and linear cases are of the same order of values     

Establishment of Numerical Wave Flume Based on the Second-Order Wave-Maker Theory

With growing computational power, the first-order wave-maker theory has become well established and is widely used for numerical wave flumes. However, existing numerical models based on the first-order wave-maker theory lose accuracy as nonlinear effects become prominent. Because spurious harmonic waves and primary waves have different propagation velocities, waves simulated by using the first-order wave-maker theory have an unstable wave profile. In this paper, a numerical wave flume with a piston-type wave-maker based on the second-order wave-maker theory has been established. Dynamic mesh technique was developed. The boundary treatment for irregular wave simulation was specially dealt with. Comparisons of the free-surface elevations using the first-order and second-order wave-maker theory prove that second-order wave-maker theory can generate stable wave profiles in both the spatial and time domains. Harmonic analysis and spectral analysis were used to prove the superiority of the second-order wave-maker theory from other two aspects. To simulate irregular waves, the numerical flume was improved to solve the problem of the water depth variation due to low-frequency motion of the wave board. In summary, the new numerical flume using the second-order wave-maker theory can guarantee the accuracy of waves by adding an extra motion of the wave board. The boundary treatment method can provide a reference for the improvement of nonlinear numerical flume.     

A Study of Bispectra in Relation to the Adequacy of Wave Spectral Model

In this paper.the characteristics of laboratory wind waves under various wind speeds and wat-er depths are studied.It is found that either the real or the imaginary part of the bispectrum can be relatedto the asymmetry of the wave profile,and the bicoherence is related to the ratio of nonlinear to linear wavecomponent.Occasionally,these two categories of nonlinear index lead to opposite inferences,because eachof them has its own significance and functions.The applicability of linear wave spectral model in oceanwaves becomes questionable only when strong nonlinearity is indicated by both of these two indexes.Thelinear spectral representation of wave fields does not necessarily become inadequate as water depth de-creases,and its appropriateness can be examined through the characteristics of the bispectrum.     

斯里兰卡南部海域次重力波特征研究

次重力波在近岸水动力运动过程中起到重要作用,深入地了解次重力波对预防其引发的近岸灾害有着重要意义。结合波浪现场观测方法和数值模拟方法研究了斯里兰卡南部海域次重力波特征。短波和自由次重力波通过海浪谱模型WAVEWATCH III模拟,而约束次重力波通过二阶非线性理论模拟,数值模拟结果与现场观测结果对比吻合良好。研究结果表明,斯里兰卡南部海域大部分时间以自由次重力波为主;在强涌浪海况下（短波波高大于2.5 m,周期大于15 s）约束次重力波逐渐逼近自由次重力波甚至占主导地位;斯里兰卡大陆架极其狭窄,对涌浪的能量损耗作用极为有限,使得近岸面对强劲的涌浪及其伴生的次重力波的侵袭,增大了次重力波引发近岸水动力灾害的风险。