随机波浪下斯托克司漂移速度是否引起水面漂移物离散的研究

漂浮于自由水面的污染物的的迁移、扩散会受到天然随机海浪的影响。之前的研究（以Herterich和Hasselmann（1982）为代表）普遍认为,随机波浪作用下的斯托克司漂移速度会引起水面污染物的离散,这个离散甚至有可能跟风和海流引起的离散同一量级。本研究就随机波浪作用下的斯托克司漂移速度是否会引起水面漂移物的离散进行理论和试验探讨。从理论推导可知,随机波浪下的质量输移速度是个定常分量,因此它不会随时间变化而引起水面漂移物的离散。随后我们在实验室水槽中进行了漂移物在随机波浪（P-M谱）作用下的漂移过程的测量。试验结果也印证了随机波浪作用下的斯托克司漂移速度不会引起水面漂移物离散的结论。     

Numerical simulation of tsunami waves generated by deformable submarine landslides

This paper presents a new submarine landslide model based on the non-hydrostatic wave model NHWAVE of Ma et al. (2012). The landslide is modeled as a water–sediment mixture. The dense plume is driven by baroclinic pressure forcing introduced by spatial density variations. The model is validated using laboratory measurements of turbidity currents and of water wave generation by a granular landslide. The model is then utilized to study the dependence of landslide motion and associated tsunami wave generation on parameters including sediment settling velocity, initial depth of the landslide and slide density. Model results show that the slide motion and water waves which it generates are both sensitive to these parameters. The relative tsunamigenic response to rigid and deformable landslides of equal initial geometry and density is also examined. It is found that the wave energy is mostly concentrated on a narrow band of the dominant slide direction for the waves generated by rigid landslides, while directional spreading is more significant for waves generated by deformable landslides. The deformable landslide has larger speed and acceleration at the early stage of landslide, resulting in larger surface waves. The numerical results indicate that the model is capable of reasonably simulating tsunami wave generation by submarine landslides.     

On the evolution and run-up of breaking solitary waves on a mild sloping beach

This paper presents new laboratory experiments carried out in a supertank (300 m × 5 m × 5.2 m) of breaking solitary waves evolution on a 1:60 plane beach. The measured data are employed to re-examine existing formulae that include breaking criterion, amplitude evolution and run-up height. The properties of shoreline motion, underwater particle velocity and scale effect on run-up height are briefly discussed. Based on our analyses, it is evidently found that there exist five zones during a wave amplitude evolution course on the present mild slope. A simple formula which is capable of predicting maximum run-up height for a breaking solitary wave on a uniform beach with a wide range of beach slope (1:15–1:60) is also proposed. The calculated results from the present model agree favorably with available laboratory data, indicating that our method is compatible with other predictive models.     

适用于曲边界的浅水非线性波浪数值模型

实际工程中存在大量的曲边界,因此在曲边界上的计算准确性可以考察出一个数值模型的实用价值。利用Beji的改进型Boussinesq方程建立了一个有限元方法的数值波浪模型。造波方面采用Fenton提出的非线性规则波浪解;在墙边界处,以求解法线方向和切线方向的速度和导数代替求解x、y方向的速度和导数,从而使边界条件直接适用、严格满足,保证了对曲边界计算的准确性。"重开始广义极小残量法"的使用保证了求解方程组的效率和精度,使造波和边界处理方法的有效性和准确性得到了合理地诠释。通过与试验数据、他人数值结果、解析解的比对,显示出该模型计算稳定、结果准确,真正体现出了有限元方法对曲边界适用的优势。     

随机波浪下泰勒离散系数的时域解

利用Wolk提出的粒子追踪方程,通过等分频率法划分不规则波谱,利用MATLAB做粒子运动模拟计算,得到无因次化泰勒离散系数K/D随时间t变化的曲线;通过与Huang等得到的P-M谱的泰勒离散系数K/D计算结果比较证明了本计算方法的可靠性。采用该方法研究了不规则波条件下,波序列(同一谱型不同波面序列)和谱型(谱峰周期、有效波高、谱峰升高因子)对波浪离散系数的影响;计算结果表明:同一谱型不同波序列对泰勒纵向离散系数稳定值和稳定时间无影响;不规则波谱峰周期越大,纵向离散系数K/D越小,稳定时间越短;有效波高越大,纵向离散系数K/D越大,稳定时间越长;谱峰升高因子越大,泰勒离散系数K/D越大,稳定时间越长;与规则波相比,不规则波的泰勒离散系数K/D的值略小10%~30%。     

Propagation of shallow water waves in an open parabolic channel using the WKB perturbation technique

A singular perturbation analysis based on the WKB technique to study the hydrodynamic performance of periodic ocean waves that are incident on an open parabolic channel of constant depth is proposed. We derive a linear model to predict the propagation of the long ocean waves into the channel. In this manner, the spatial distribution for the surface elevation of the ocean waves inside the channel as a function of two dimensionless parameters, namely, a kinematical parameter, κ and a geometrical parameter ε, is governed by a second-order ordinary differential equation. The kinematical parameter κ denotes the ratio of the potential head, due to gravity, to the kinetic head of the ocean waves along the longitudinal axis of the parabolic channel. Meanwhile, ε is a dimensionless geometrical parameter that represents a characteristic ratio of the parabolic channel. Using matching conditions, simple expressions for the reflection and transmission coefficients are obtained.     

短峰畸形波生成、演化过程的外部特征研究

在实际海况中,由于波浪组分方向不同,使海浪常表现为三维短峰波,与二维长峰波相比,有诸多异同特征。为了研究三维畸形波生成和演化过程中异常大波的外部特征以及大波间的关联,本文采用色散和方向聚焦方法在物理水池中模拟了短峰畸形波的整个生成和演化过程。研究结果表明,生成和演化过程中出现的异常大波沿着方向分布函数峰值方向演化,具有一定程度“对称”特征;根据异常大波的特征参数可将该过程分成3～4个典型的阶段,这意味着可以根据异常大波的特征参数判断其所处阶段及预测其演化趋势;三维畸形波生成和演化过程可能会跳过“波群”或者“深谷”阶段,与二维畸形波相比,时−空范围较小,这意味着在三维波浪场中,与畸形波相关的深谷和波群的发生概率较低。     

Numerical modelling of deformation of multi-directional random waves over a varying topography

Two numerical formulations of the breaking phenomenon were implemented in a numerical model for random wave propagation based on the elliptic formulation of the mild-slope equation. The randomness of the wave field was simulated based on a spectral component method, in which the 3-D spectrum is discretised in components of equal energy. One of the breaking process formulations is based on the concept of breaking each independent spectral component. The other is based on the distribution of the local amount of energy dissipated through the independent spectral components. The model based on the concept of breaking each independent spectral component produces the best estimates of the wave field, when the numerical results are compared with laboratory data.     

福建长乐人工震源波的观测与分析

使用DCS-302数字磁带记录地震仪在近场(震中距△≤5km)观测人工爆破源产生的地震波,在不同方位、不同震中距和土层条件下得到了9条近场地面运动加速度记录。通过回放及计算机处理研究了近场爆破地震波记录波形特征、频谱特性、振幅特征、衰减情况以及加速度最大幅值和相应的卓越周期。     

Tsunami-like solitary waves impinging and overtopping an impermeable seawall: Experiment and RANS modeling

This study investigates tsunami-like solitary waves impinging and overtopping an impermeable trapezoidal seawall on a 1:20 sloping beach. New laboratory experiments are performed for describing three typical cases: a turbulent bore rushes inland and subsequently impacts and overtops the seawall (Type 1); a wave directly collapses on the seawall and then generates overtopping flow (Type 2); and, a wave straightforwardly overtops the seawall crown and collapses behind the seawall (Type 3). A two-dimensional volume of fluid (VOF) type model called the COBRAS (COrnell BReaking And Structure) model, which is based on the Reynolds-Averaged Navier–Stokes (RANS) equations and the k–ε turbulence closure solver, is validated by experimental data and then applied to investigate wave dynamics for which laboratory data are unavailable. Additionally, a set of numerical experiments is conducted to examine the dynamic wave acting force due to waves impacting the seawall. Effects of wave nonlinearity and freeboard are elucidated. Special attention is given to a distinct vortex evolutionary behavior behind the seawall, in which the dynamic properties of entrapped air-bubbles are briefly addressed experimentally and numerically.     

深水波浪破碎时波浪演化特征实验研究

采用能量聚焦的方式产生深水破碎波,并通过增加输入波陡使发生不同强度的波浪破碎现象。实验中,沿水槽中心位置布置22个浪高仪,分析波浪传播过程中的波面演化特征。对水槽不同位置处波面数据进行波能谱与小波能谱分析,发现在聚焦波传播过程中,低频能量部分保持相对稳定,而一次谐波高频部分先逐渐拓宽,经过破碎区域后又逐渐恢复。能量在高频部分有所损失,这种现象在破碎时更加明显,且破碎强度越大,越显著。波浪未破碎时,由于波浪传播过程中高频部分拓宽,导致聚焦前后特征频率略有增加,特征群速和特征周期略有减小;当波浪破碎时,由于破碎导致的能量损失比较明显,且卷破时更加明显,导致破碎后特征频率减小,特征群速和特征周期增大。     

Measurement and modeling of bed shear stress under solitary waves

Direct measurements of bed shear stresses (using a shear cell apparatus) generated by non-breaking solitary waves are presented. The measurements were carried out over a smooth bed in laminar and transitional flow regimes (~ 10⁴ < R_e < ~ 10⁵). Measurements were carried out where the wave height to water depth (h/d) ratio varied between 0.12 and 0.68; maximum near bed velocity varied between 0.16 m/s and 0.51 m/s and the maximum total shear stress (sum of skin shear stress and Froude–Krylov force) varied between 0.386 Pa and 2.06 Pa. The total stress is important in determining the stability of submarine sediment and in sheet flow regimes. Analytical modeling was carried out to predict total and skin shear stresses using convolution integration methods forced with the free stream velocity and incorporating a range of eddy viscosity models. Wave friction factors were estimated from skin shear stress at different instances over the wave (viz., time of maximum positive total shear stress, maximum skin shear stress and at the time of maximum velocity) using both the maximum velocity and the instantaneous velocity at that phase of the wave cycle. Similarly, force coefficients obtained from total stress were estimated at time of maximum positive and negative total stress and at maximum velocity. Maximum positive total shear stress was approximately 1.5 times larger than minimum negative total stress. Modeled and measured positive bed shear stresses are well correlated using the best convolution model, but the model underestimates the data by about 4%. Friction factors are dependent on the choice of normalizing using the maximum velocity, as is conventional, or the instantaneous velocity. These differ because the stress is not in phase with the velocity in general. Friction factors are consistent with previous data for monochromatic waves, and vary inversely with the square-root of the Reynolds number. The total shear stress leads the free stream fluid velocity by approximately 50°, whereas the skin friction shear stress leads by about 30°, which is similar to that reported by earlier researchers.     

On low-frequency waves in the surf and swash

Low-frequency waves in the surf and swash zones on various beach slopes are discussed using numerical simulations. Simulated surface elevations of both primary waves and low-frequency waves across the surf zone were first compared with experimental data and good agreement found. Low-frequency wave characteristics are then discussed in terms of their physical nature and their relationship to the primary wave field on a series of sea bottom slopes. Unlike primary waves, low-frequency wave energy increases towards the shoreline. Low-frequency waves in the surf and swash are a function of incident waves and the sea bottom slope and hence the saturation level of the surf zone. Wave energy on a gently sloping beach is dominated by low-frequency waves while primary waves play a significant role on a steep beach. Low-frequency wave radiation from the surf zone on a given beach depends on primary wave frequency and beach slope. However, a very poor correlation was found between surf similarity parameter and low-frequency wave radiation.     

Effect of bed roughness on turbulent boundary layer and net sediment transport under asymmetric waves

This paper presents an investigation of the roughness effects in the turbulent boundary layer for asymmetric waves by using the baseline (BSL) k–ω model. This model is validated by a set of the experimental data with different wave non-linearity index, N_i (namely, N_i = 0.67, N_i = 0.60 and N_i = 0.58). It is further used to simulate asymmetric wave velocity flows with several values of the roughness parameter (a_m/k_s) which increase gradually, namely from a_m/k_s = 35 to a_m/k_s = 963. The effect of the roughness tends to increase the turbulent kinetic energy and to decrease the mean velocity distribution in the inner boundary layer, whereas in the outer boundary layer, the roughness alters the turbulent kinetic energy and the mean velocity distribution is relatively unaffected. A new simple calculation method of bottom shear stress based on incorporating velocity and acceleration terms is proposed and applied into the calculation of the rate of bed-load transport induced by asymmetric waves. And further, the effect of bed roughness on the bottom shear stress and bed-load sediment transport under asymmetric waves is examined with the turbulent model, the newly proposed method, and the existing calculation method. It is found that the higher roughness elements increase the magnitude of bottom shear stress along a wave cycle and consequently, the potential net sediment transport rate. Moreover, the wave non-linearity also shows a big impact on the bottom shear stress and the net sediment transport.     

An integrated numerical method for simulation of drifted objects trajectory under real-world tsunami waves

The present study focused on tracing tsunami-drifted objects under a real tsunami based on an integrated numerical method. Instead of a solitary wave that is much shorter and steeper than real-world tsunami waves, an extra-long tsunami wave is represented here in a nearshore region using a new approach. To this end, propagation of a seismic tsunami from the source to the nearshore region was simulated using two-dimensional depth-averaged equations. When the waves reached the target coastal area, the time series of the free surface of the tsunami was approximated by a theoretical relation based on a combination of several solitons, which were then used to solve the linearized trajectory equation of the wave-maker to generate the intended time series of the tsunami wave. Finally, in a nearshore model, the movement of drifted bodies under the generated tsunami wave was simulated based on the smoothed-particle hydrodynamics (SPH) method. In order to verify the accuracy of the proposed method in tracing the drifted bodies under a real tsunami, the giant fish-oil tank, which was transported about 300 m during the 2011 Tohoku tsunami of Japan, was selected as the benchmark. The results demonstrate that the time series of the long tsunami wave was successfully generated by the piston wave-maker in the GPU-based SPH model, and the proposed approach can be regarded as a suitable alternative for reproduction of a real tsunami. The results also showed that the simulated fish-oil tank properly followed the estimated trajectory in Ishinomaki but it was transported more than the reported distance, which was expected due to absence of a holding connection between the tank and the ground in the SPH model. It should be emphasized that this study is one of the first studies on three-dimensional tracing of a tsunami-drifted body during a real event, and the tracing can be more accurate in further simulations by applying higher-resolution topography data and faster computation systems that help include more details in the nearshore model.     

Extended mild-slope equation for random waves

A time-dependent extended mild-slope equation is derived from the elliptic equation of Chamberlain and Porter [J. Fluid Mech. 291 (1995) 393] using the Taylor series technique. Numerical tests are made on a horizontally one-dimensional case for regular waves over sloping beds and for both regular and irregular waves over a ripple patch. Numerical results prove that the proposed model gives accurate results for both regular and irregular waves over rapidly varying topography.     

Wave types of landslide generated impulse waves

Subaerial landslide generated impulse waves were investigated in a prismatic wave channel. Seven governing parameters, namely the still water depth, slide impact velocity, slide thickness, bulk slide volume, bulk slide density, slide impact angle, and grain diameter, were systematically varied. The generated impulse waves are nonlinear, intermediate- to shallow-water waves involving a small to considerable fluid mass transport. The Stokes wave, cnoidal wave, solitary wave, and bore theories were applied to describe the observed maximum waves. The theoretical and observed features of these four wave types are highlighted. A diagram allows to predict the wave type directly as a function of the slide parameters, the slide impact angle, and the still water depth.     

可控深水破碎波的实验室生成方法研究

实验室一般采用波浪聚焦方法生成深水破碎波,通过各组分波浪的波幅叠加生成一个波高显著增大的大波,使其波陡超过极限波陡发生破碎。利用该方法生成深水破碎波浪的破碎次数通常并不唯一,导致波浪破碎后的流场特征不明显;造波参数不易于选取导致研究工况的设置难度大,直接影响深水破碎精细化实验的效果和效率。本文采用聚焦波理论计算波面,并利用上跨零点法定义的波高和波长计算理论波陡,结合物理模型实验统计波浪沿程破碎次数与剧烈程度,研究以JONSWAP谱为造波输入谱型时,聚焦波幅、谱峰频率、频宽等造波输入参数对于波浪破碎情况的影响,从而建立深水波浪破碎次数与造波输入参数之间的近似定量关系,为实验造波参数的选取提供参考,提高实验效率。     

不同破碎波对沙质海床作用的实验研究

破碎波对近海海岸地形以及海岸建筑物影响强烈,通过物理模型实验对孤立波、规则波作用下破碎带的床面形态以及孔隙水压力进行分析。破碎波冲击海床,破碎处床面上形成沙坝和沙坑,与规则波相比,孤立波破碎时对床面的冲刷更加剧烈,床面形成的沙坝和沙坑尺度更大,且土体内孔隙水压力幅值也较大。同时研究了波面变化对孔隙水压力的影响,发现波面变化历时曲线与孔隙水压力历时曲线相似,与孔隙水压力梯度历时曲线更为相似,说明波面变化更能反映海床内部孔隙水压力梯度的变化。通过探讨波浪与海床之间相互耦合作用,发现破碎带地形变化使得波浪出现不同破碎类型,分析得出卷破波比崩破波作用下孔隙水压力幅值大。     

近岸波浪破碎区不规则波浪的数值模拟

基于近岸不规则波浪传播的抛物型缓坡方程和两类波浪破碎能量损耗因子,对近岸波浪破碎区不规则波浪的波高分布进行了数值模拟,并结合实验结果对数值模拟结果进行了验证分析,结果表明采用两类波浪破碎能量损耗因子所模拟的破碎区波高与实测值均吻合良好,波浪破碎能量损耗因子及波浪破碎指标对破碎区波浪波高分布影响较明显。