动波浪壁圆柱绕流三维数值模拟

利用计算流体力学软件Fluent开展了三维动波浪壁圆柱绕流的数值模拟,建立了三维运动波浪壁圆柱模型,通过C语言自编程序实现波浪壁面的运动控制,并保证壁面变形时网格的高质量。在来流速度u=0.125 m/s、雷诺数Re=12 500的情况下,开展了动波浪壁波动速度w=0、0.062 5、0.125、0.187 5 m/s四个工况的计算分析,并比较了不同波动速度对流场结构、升力、阻力特性的影响。结果表明:动波浪壁圆柱能有效抑制流动的分离,消除交替脱落的尾涡,从而消除周期振荡的升力;在消除卡门涡街的同时,圆柱后驻点处的涡量值随波动速度增加而增加,其原因在于波形移动加大了壁面流体的速度,从而减小了圆柱前后的压力差,减小了阻力;随着波动速度的增大,平均阻力系数呈明显下降趋势,当波动速度为来流速度的1.5倍时,平均阻力系数相对于光滑圆柱下降了53.76%。     

Generation of internal waves by moving non-stationary disturbances in a real stratified ocean

The plane problem on the generation of linear internal waves by a moving area of time-harmonic surface pressures in a continuously-stratified ocean of constant depth is considered. An analytical relation has been derived for forced internal waves off the site of their generation in the form of an internal wave field superposition corresponding to individual vertical modes. The possible wave regimes are determined. For the Brunt-Väisälä frequency distribution in the North Atlantic, the generation conditions and amplitudes of diverse radiated waves are numerically determined.Translated by Vladimir A. Puchkin.     

Transmission and reflection of the Rossby waves in a stratified ocean with bottom topography

Transmission and reflection coefficients are calculated for Rossby waves incident on a bottom topography with constant slope in a continuously stratified ocean. The characteristics of the coefficients are interpreted in terms of the quasigeostrophic waves on the slope. In the parameter range where only the barotropic Rossby waves can propagate in the region outside the slope, the bottom trapped wave plays the same role as the topographic Rossby wave in a homogeneous ocean, and hence the transmission is weak unless phase matching takes place. When both of the barotropic and baroclinic Rossby waves can propagate outside the slope, the total transmission can be strong. The bottom trapped wave affects the transmission and reflection, and it leads to the possibility that the Rossby wave is transmitted as a mode different from the incident mode. When the number of the wavy modes on the slope is smaller than that of the Rossby wave modes outside the slope, strong reflection occurs.The results for an ocean with linear distribution of the squared Brunt-Väisälä frequency are compared to those in a uniformly stratified ocean. The weakening of the stratification near the bottom is almost equivalent to reducing the effect of the slope.     

A spectral model for unidirectional nonlinear wave propagation over arbitrary depths

A weakly-nonlinear and dispersive wave equation recently developed by the authors is used for formulating a spectral-type unidirectional wave propagation model describing spectral transformations of narrow-band waves travelling over arbitrary depths. The essential characteristics of the model equation are recapitulated first and then the spectral domain representation in terms of spatially varying harmonic amplitudes is presented. The resulting evolution equations are used to simulate the experiments concerning harmonic generation in shallow water and nonlinear random wave transformations over a submerged bar. Furthermore, the spectral model predictions are compared with the field measurements in nearshore with satisfactory results.     

Shelf waves with diurnal tidal frequency at the Greenland shelf edge

Tidal analysis was carried out on current measurements at a “cross-shelf” transect off Greenland at 71°N. The diurnal tides manifest themselves mainly as a barotropic continental shelf wave, travelling southward along the shelf slope. This follows from the amplitude distribution of the diurnal tidal components and from the rotation sense of the tidal ellipses at different cross-slope locations, as calculated with simple two-dimensional models. The well organized cross-slope pattern of the velocity amplitudes is absent in observations further north near 75°N. These observations suggest that the local vanishing of the group velocity, which is caused by topography, is of importance for the existence and local amplification of these continental shelf waves with diurnal tidal frequency.     

Second-order wavemaker theory for multidirectional waves

The present paper develops the complete second-order wavemaker theory for the generation of multidirectional waves in a semi-infinite basin. The theory includes superharmonics and subharmonics and is valid for a rotational as well as a translatory serpent-type wave-board motion. The primary goal is to obtain the second-order motion of the wave paddles required to get a prescribed multidirectional irregular wave field correct to second order, i.e. to suppress spurious free-wave generation. The wavemaker theory is a 3D extension of the full second-order wavemaker theory for wave flumes by Schäffer (1996).     

利用WAVEWATCH和SWAN嵌套计算珠江口附近海域的风浪场*

基于CCMP(Cross Calibrated Multi-platform)卫星遥感海面风场数据,通过将WAVEWATCH和SWAN (Simulating WAves Nearshore)模型嵌套的方法,数值模拟了珠江口附近海域的风浪场。将总计10个月的数值模拟的有效波高、波周期和波向分别与相应的观测值进行了定量比较。结果说明,有效波高的平均绝对误差为15.4cm,分散系数SI为0.240,相关系数为0.925;波周期的平均绝对误差为1.9s,分散系数SI为0.433,相关系数为0.636;波向的平均绝对误差为23.9°。计算的波高和波向与观测结果的变化趋势相吻合。由于第三代海浪模式本身的缺陷,导致所计算的波周期偏小。总体说来,本文所采用的数值模式能较好地模拟珠江口附近海域的风浪场。另外,还设计了6个算例以探讨采用不同的计算方法和风场对计算结果精度的影响。结果表明使用本文的数值方法和高精度的CCMP风场确实可以提高计算结果的精度。     

河口、海岸环境下床沙输运的初步实验研究

量测了在水深0.5m的情况下波浪、水流通过水平沙床所产生的床沙输移。实验波高为0.15m，周期分别为1.4和2.0s。首先量测纯波浪下的床沙输移，然后量测波流共同作用下，水流与波浪行进方向一致其稳定速度分别为0.02,0.04,0.06m/s时的床沙输移，结果表明，2个沙槽所得的总输沙率在波周期2.0s时最大，净输沙率在波周期1.4s时最大，将波叠加在速度为0.02m/s的水流上时，2种波型的净输沙率都增加约1倍，水流为0.04m/s和0.06m/s时2种波型的净输沙率分别都减少。     

层结流体中的非线性惯性重力内波

根据地球流体力学基本方程组,在密度垂直层结的情况下,引进行波坐标,研究非线性定形波在相平面上的几何拓扑结构。严格论证了不存在定形孤立波,并通过Ｈａｍｉｌｔｏｎ函数及其角作用变把行波系统化成最简形式,由此而得到非线性惯性重力内波的解析解。     

Generation of spatial internal waves in a rotating ocean by moving non-stationary disturbances

In a general linear statement, the kinematic structure of spatial internal waves generated by a uniformly moving area of oscillating surface pressures in a continuously stratified ocean of constant depth is studied. The earth's rotation effects are considered. Possible types of ocean wave fields with a constant Brunt-Väisälä frequency are examined. The wave regimes for individual modes of internal and gyroscopic waves are classified on the basis of estimating the integrals asymptotically.Translated by V. Puchkin.     

A coupled-mode system with application to nonlinear water waves propagating in finite water depth and in variable bathymetry regions

A non-linear coupled-mode system of horizontal equations is presented, modelling the evolution of nonlinear water waves in finite depth over a general bottom topography. The vertical structure of the wave field is represented by means of a local-mode series expansion of the wave potential. This series contains the usual propagating and evanescent modes, plus two additional terms, the free-surface mode and the sloping-bottom mode, enabling to consistently treat the non-vertical end-conditions at the free-surface and the bottom boundaries. The present coupled-mode system fully accounts for the effects of non-linearity and dispersion, and the local-mode series exhibits fast convergence. Thus, a small number of modes (up to 5–6) are usually enough for precise numerical solution. In the present work, the coupled-mode system is applied to the numerical investigation of families of steady travelling wave solutions in constant depth, corresponding to a wide range of water depths, ranging from intermediate depth to shallow-water wave conditions, and its results are compared vs. Stokes and cnoidal wave theories, as well as with fully nonlinear Fourier methods. Furthermore, numerical results are presented for waves propagating over variable bathymetry regions and compared with nonlinear methods based on boundary integral formulation and experimental data, showing good agreement.     

Diffraction of a long Kelvin wave at the apex angle formed by intersecting shores

This paper discusses the diffraction of a low-frequency Kelvin wave (having a frequency comparable to that of a baroclinic Rossby wave) at the apex angle. The problem is handled, with the Coriolis parameter being variable, by expanding into series over the minor parameter ε characteristic of the ratio between the deformation radius and the length of an incident wave travelling along a latitudinal shore. An analytical expression has been derived for a wave translating along meridional shore. It is demonstrated that when the latitudinal shore is located above the critical latitude, the energy flux through the horizontal cross-section transported by the wave coincides with the corresponding flux in the incident wave. If, on the other hand, it is situated below the critical latitude, then part of the energy is consumed by offshore Rossby waves. Translated by Vladimir A. Puchkin.     

不同结构形式潜堤上的随机波浪运动

潜堤在海岸保护和促潍保淤中具有重要的意义。工程效益和经济投入受到潜堤结构形式和高度的影响。基于RANS方程和流体体积（VOF）方法的模型,计算不同潜堤上的波浪运动。潜堤形式包括矩形结构、梯形结构、半圆结构和槽形结构。模拟了不规则波越堤时波浪的破碎特征、破波类型以及波高的变化。通过分析流场和波高分布的变化,得到潜堤上波能散耗的规律,并且得出最优潜堤结构形式。讨论了相对堤顶高度、破碎点位置、相对堤顸宽度、波陡和潜堤结构形式等对波能耗散的影响。     

Vortex-Induced Vibrations of A Long Flexible Cylinder in Linear and Exponential Shear Flows

A numerical study based on a wake oscillator model was conducted to determine the response performance of vortex-induced vibration(VIV) on a long flexible cylinder with pinned-pinned boundary conditions subjected to linear and exponential shear flows. The coupling equations of a structural vibration model and wake oscillator model were solved using a standard central finite difference method of the second order. The VIV response characteristics including the structural displacement, structural frequency, structural wavenumber, standing wave behavior,travelling wave behavior, structural velocity, lift force coefficient and transferred energy from the fluid to the structure with different flow profiles were compared. The numerical results show that the VIV displacement is a combination of standing waves and travelling waves. For linear shear flow, standing waves and travelling waves dominate the VIV response within the low-velocity and high-velocity zones, respectively. The negative values of the transferred energy only occur within the low-velocity zone. However, for exponential shear flow, travelling waves dominate the VIV response and the negative energy occurs along the entire length of the cylinder.     

基于Holland台风模型及三重嵌套海浪模式的台风浪数值模拟研究

现有的风场资料存在台风中心附近风速偏低的问题。为改进台风期间风场数据, 使用Holland经验台风模型结合多平台交叉校准数据(cross-calibrated multi-platform, CCMP)及欧洲中期天气预报中心的再分析数据(European Centre for Medium-range Weather Forecasts Reanalysis data, ERA5)风场资料, 研究了不同台风最大风速半径(maximum wind radius of the typhoon, RMW)、Holland B参数对模拟效果的影响, 确定了最优模拟参数, 并以改进后的风场驱动三重嵌套海浪模型对台风“威马逊”发生期间的台风浪进行模拟。模拟结果与实测数据对比表明, (1)改进的风场资料与实测结果更为接近, 作为海浪模式驱动项可更好地模拟台风期间波浪状况; (2)三重嵌套海浪模型的波浪模拟效果优于单独的海浪模型。     

The statistical distribution of nearbed wave orbital velocity in intermediate coastal water depth

The statistical distribution of wave orbital velocity in intermediate coastal water depth has been quantitatively determined from the comprehensive field velocity data collected near the seabed in this study. Two ocean ADV current meters, which were mounted at 0.5 m above the seabed on two separate stainless steel tripods sitting on the seabed, were used to measure instantaneous water particle velocities at a 2 Hz sampling rate for 17.07 min every hour in two coastal water depths of 11 m and 23 m in nine field deployments over a period of 2 years. The zero-crossing method is applied to analyse the field velocity data collected in each field deployment to obtain a large sample of wave orbital velocity amplitudes of individual waves. Based on the collected field velocity data, it is found that the histogram of instantaneous wave orbital velocities perfectly follows the Gaussian distribution as commonly assumed, while the histogram of wave orbital velocity amplitudes is less accurately described by the Rayleigh distribution than the modified Rayleigh and the Weibull distribution. It is also found that large orbital velocity amplitudes are generally overestimated by the Rayleigh distribution, but well predicted by the modified Rayleigh and the Weibull distribution. The expected value of maximum orbital velocity in a velocity record of finite size is also derived from the three distributions and found to agree well with the present field data.     

Far fields of internal gravity waves from oscillating sources of disturbances

We considered the problem dealing with the far field of internal gravity waves generated by the source oscillating in the vertical. It is shown that, at a fixed observation point, we first record the arrival of the main wave front of a separate mode followed by the arrival of two fronts of an additional field of this wave mode, which is caused by the oscillations of the source, the asymptotics near the additional wave front being expressed via Airy function.     

Generation of a double Kelvin wave resulting from the reflection of the Rossby wave by the jump layer

This paper considers the evolution of a spatially-localized divergent Rossby wave field near the depth jump. If the jump magnitude is comparable to the depth, Rossby waves are fully reflected and a double Kelvin wave is then generated. The Rossby waves and the double Kelvin wave are described by the first- and zero-approximation fields of the asymptotic expansion, respectively. Over the characteristic Rossby wave period, the level elevation produced by the double Kelvin wave spreads over an extensive area, theraby making up for the change in the total fluid mass of the Rossby waves.Translated by Vladimir A. Puchkin.     

A method of submersion in the problem on the determination of dispersion curves of internal waves in the presence of a current with vertical velocity shear

The previously proposed method [1] of submerging aimed at determining internal wave dispersion characteristics is generalized to consider the case of a background horizontal flow with a vertical velocity shear. The results of calculating, dispersion curves of discrete spectrum modes for the profiles of the Brunt-Väisälä frequency and the mean current velocity are given; they agree well with the results of studies carried out in ref. [5].Translated by Mikhail M. Trufanov.     

Internal flow structure of short wind waves

The internal flow structure of wind waves in a wind-wave tunnel was investigated on the bases of the measured vorticity distributions, streamline patterns, internal pressure fields, and stress distributions at the water surface for some waves in the field. In part I the experimental method and the internal vorticity structure relative to the individual wave crests are described. The measured vorticity distributions of distinct waves (waves with waveheight comparable with or larger than that of significant wavesH _1/3) in the field indicate that the surface vorticity layer is extraordinarily thickened near the crest, and the vorticity near the water surface shows a particularly large value below the crest. The flow near the crest of distinct waves is found to be in excess of the phase speed in a very thin surface layer, and the tangential stress distribution has a dominant peak near the crest. It is argued that the occurrence of the region of high vorticity in distinct waves is associated with the local generation of vorticity near the crest by tangential stress which attains a peak, under the presence of excess flow.