The vertical structure of combined wave–current flow

A simple numerical model, based on the Reynolds stress equations and k–ε turbulence closure scheme, is developed for the coastal wave and current bottom boundary layer. The current friction velocity is introduced to account for the effect of currents on waves. The implicit Crank–Nicolson finite difference method discretizes the governing equations. Vertical changing step grids with the constant ratio for two adjacent spatial steps are used together with the equal time steps in the modeling. Vertical profiles of mean current velocity and wave velocity amplitude are obtained. These modeled results are compared with the laboratory experimental data of Van Doorn [1981. Experimental investigation of near bottom velocities in water waves with and without a current. Report M1423, Delft Hydraulics Laboratory, Delft, The Netherlands; 1982. Experimenteel onderzoek naar het snelheidsveld in de turbulente bodemgrenslaag in een oscillerende stroming in een golftunnel. Report M1562, Delft Hydraulics Laboratory, Delft, The Netherlands]. It has been shown that modeled and observed (Van Doorn, T., 1981. Experimental investigation of near bottom velocities in water waves with and without a current. Report M1423, Delft Hydraulics Laboratory, Delft, The Netherlands; 1982. Experimenteel onderzoek naar het snelheidsveld in de turbulente bodemgrenslaag in een oscillerende stroming in een golftunnel. Report M1562, Delft Hydraulics Laboratory, Delft, The Netherlands) mean velocity profiles within the wave and current bottom boundary layer are in better agreement than outside. Modeled and observed (Van Doorn, T., 1981. Experimental investigation of near bottom velocities in water waves with and without a current. Report M1423, Delft Hydraulics Laboratory, Delft, The Netherlands) wave velocity amplitude profiles within the wave and current bottom boundary layer are in better agreement than outside. Modeled wave velocity amplitudes are in good agreement with the laboratory experimental data of Van Doorn [1982. Experimenteel onderzoek naar het snelheidsveld in de turbulente bodemgrenslaag in een oscillerende stroming in een golftunnel. Report M1562, Delft Hydraulics Laboratory, Delft, The Netherlands].     

Field investigations of the drift of artificial thin films on the sea surface

The results of field investigations of the drift of surface films are presented. The experiments were performed at the southern coast of the Crimea in the area of an oceanographic platform. It is shown that the drift of a thin film is described well by a linear combination of the vectors of the coastal current and wind, and the coefficient of the current velocity is close to unity. The total velocity of wind and wave drifts is 0.013 of the wind velocity magnitude or 0.38 of the air friction velocity. Analysis of experimental data demonstrates an increase in the film spreading rate with increasing wind velocity.     

波浪作用下异重流运动特性研究

本文通过水槽试验和理论分析手段对波浪作用下异重流运动特性进行了研究，应用流体运动基本方程推导了异重流在水平底面运动时的运动速度、厚度和流量表达式。结果表明，计算结果与水槽试验较一致     

3D scour below pipelines under waves and combined waves and currents

This paper presents the results of an experimental investigation on three-dimensional local scour below a rigid pipeline subjected to wave only and combined wave and current conditions. The tests were conducted in a conventional wave flume. The major emphasis of the investigation was on the scour propagation speed (free span expansion rate) along the pipeline after local scour was initiated at a controlled location. The effects of flow ratio (steady current velocity vs. combined waves/current velocity), flow incidence angle and pipeline initial embedment depth on free span expansion rate were investigated. It was observed that the scour along the pipeline propagated at a constant rate under wave only conditions. The scour propagation rate decreased with increasing embedment depth, however, increased with the increasing Keuglegan–Carpenter (KC) number. Under combined wave and current conditions, the effect of velocity ratio on scour propagation velocity along the pipeline was quantified. Empirical relationships between the scour propagation rate (V_h) and key parameters such as the KC number and embedment depth (e/D) were established based on the testing results.     

浅化波浪层流边界层流速分布特性的数值分析

建立了同时考虑波致雷诺应力和时均水平压强梯度影响的二阶波浪边界层数学模型,模型计算得到的浅化波浪层流边界层内瞬时流速剖面、振荡速度幅值和时均流速剖面均与水槽实验数据吻合较好,在此基础上探讨了浅化波浪边界层流速分布特性及其影响机制。随着波浪的浅化变形,边界层内时均流速剖面"底部向岸、上部离岸"的变化特征越来越明显。这是二阶对流项引起的波致雷诺应力和离岸回流引起的时均水平压强梯度共同作用的结果,在床面附近由波致雷诺应力占主导作用并趋于引起向岸流动,在上部区域由时均水平压强梯度占主导作用并趋于引起离岸流动。     

Numerical simulation of wave-induced laminar boundary layers

The three-dimensional numerical model with σ-coordinate transformation in the vertical direction is applied to the simulation of surface water waves and wave-induced laminar boundary layers. Unlike most of the previous investigations that solved the simplified one-dimensional boundary layer equation of motion and neglected the interaction between boundary layer and outside flow, the present model solves the full Navier–Stokes equations (NSE) in the entire domain from bottom to free surface. A non-uniform mesh system is used in the vertical direction to resolve the thin boundary layer. Linear wave, Stokes wave, cnoidal wave and solitary wave are considered. The numerical results are compared to analytical solutions and available experimental data. The numerical results agree favorably to all of the experimental data. It is found that the analytical solutions are accurate for both linear wave and Stokes wave but inadequate for cnoidal wave or solitary wave. The possible reason is that the existing analytical solutions for cnoidal and solitary waves adopt the first-order approximation for free stream velocity and thus overestimate the near bottom velocity. Besides velocity, the present model also provides accurate results for wave-induced bed shear stress.     

珊瑚礁海岸波流运动特性整体物理模型实验研究

根据现场地形在港池中建立三维珊瑚礁?潟湖?裂口海岸定床整体物理模型,采用波高传感器、流速仪和表面流速测量系统分别测量了规则波作用下珊瑚礁海岸不同位置的波浪和流场特征。结果表明:礁坪上,波高在向岸方向逐渐减小,总减小幅度为86.7%,增水先增大后减小,沿礁坪下降幅度为65.9%,水流以向岸流为主,存在着先增大后减小的趋势;潟湖中,波高靠近裂口处较大,中部最大值约为两侧最小值的2.8倍,增水则靠近裂口处最小,相比两侧最大值下降了25.5%,水流主要为对称地指向裂口的沿岸流,流速从两侧到裂口先增大后减小;裂口中波高变化不大,增水在靠近潟湖处最大,为礁坪上增水的47.6%,水流主要为离岸流,流速同样是先增大后减小。量化分析了环流驱动力、辐射应力与波面压力梯度的空间变化规律,发现礁坪上向岸流变化是平均水位梯度和辐射应力相互作用的结果,在裂口中的离岸流驱动力主要为辐射应力,而潟湖中的沿岸流变化由平均水位梯度决定的。     

Numerical study of pollutant movement in waves and wave-induced long-shore currents in surf zone

Water waves, wave-induced long-shore currents and movement of pollutants in waves and currents have been numerically studied based on the hyperbolic mild-slope equation, the shallow water equation , as well as the pollutant movement equation, and the numerical results have also been validated by experimental data. It is shown that the long-shore current velocity and wave set-up increase with the increasing incident wave amplitude and slope steepness of the shore plane ; the wave set-up increases with the in- creasing incident wave period;and the pollutant morement proceeds more quiekly with the increasing incident wave amplitude and slope steepness of the shore palane. In surf zones, the long-shore currents induced by the inclined incident waves have effectively affected the pollutant movement.     

基于TerraSAR-X卫星数据的内孤立波参数和海表流速信息提取的探索研究

本文利用TerraSAR-X(TSX)卫星于2010年4月22日在南海东沙岛附近海域获取的数据进行海洋内孤立波动力要素和海表流速信息的提取研究。基于TSX数据的后向散射强度信息,利用经验模态分解法得到内孤立波半波宽度,再利用两层模型法和参数化法计算得到内孤立波振幅和相速度。反演结果显示,利用参数化方法得到的振幅(约21~39 m)和两层模型法得到的相速度(约1.07 m/s)与历史实测资料较为一致。进而利用TSX的顺轨干涉数据获取研究海域内的多普勒速度,再分别采用M4S模型法和直接分离法处理,进而提取海表流速。结果显示,两种方法得到的海表流速的全场平均值较为一致,均为1.10 m/s左右。M4S模型法对流速最大值的改变量较大而直接分离法对流速最小值的改变量较大。M4S模型对内孤立波波峰线区域海表流速的修正大于无内孤立波的海域。最后,基于KdV方程计算得到内孤立波引起的表面流的流速约为0.28 m/s,对反演出的海表流速贡献占比23%。     

Laboratory measurements of uniform longshore currents

Experiments on uniform longshore currents in a wave basin are described. The measurements were performed in a basin with a pumped recirculation through openings in the wave guides. Minimal return flows in the offshore region of the basin are found to be accompanied with longshore currents which are virtually uniform alongshore. Three-dimensional longshore current velocity distributions were measured with much attention to quality control. Detailed experimental results are presented for different wave fields, two beach slopes and two beach roughnesses.     

Numerical Simulation of Local Scour Around A Large Circular Cylinder Under Wave Action

A horizontal two- dimensional numerical model is developed for estimation of sediment transport and sea bed change around a large circular cylinder under wave action. The wave model is based on an elliptic mild slope equation. The wave-induced current by the gradient of radiation stress is considered and a depth integrated shallow water equation is applied to the calculation of the current. The mass transport velocity and the bed shear stress due to streaming are considered, which are important factors affecting the sediment transport around a structure due to waves, especially in reflective areas. Wave-current interaction is taken into account in the model for computing the bed shear stress. The model is implemented by a finite element method. The results of this model are compared with those from other methods and agree well with experimental data.     

The effect of wave-induced stress on current profiles

Current velocity profiles in the presence of non-breaking waves on a horizontal bottom are studied. Particular consideration is given to the derivations of measured current profiles from the standard logarithmic profiles near the mean water surface. The deviations are found to be due mainly to the wave-induced second-order stress which was generally neglected in the former models. The available experimental data indicate that the wave-induced second-order stress is a linear function of elevation and depends on the wave parameters, the current strength and the angle between the waves and the current. A semi-empirical model is developed and gives good agreements with experimental measurements of current profiles near the mean water surface.     

荣喜近岸区波流耦合作用下的流场模拟

一方面将波浪对底部剪切应力、表面拖曳力系数,辐射应力以及表面混合长度的影响引入至COHERENS.另一方面又将水动力模型COHERENS和第三代波浪模型SWAN耦合,使两模型能够随时互相交换水流、水位以及波浪信息,最终获得波流耦合模型COHERENS-SWAN并将其应用于荣喜近岸区波流共同存在情况下的波流作用模拟研究.计算所得的流速、流向和水位与实测数据吻合较好.     

波流共同作用下海底人工边坡动态响应分析

为了研究波流共同荷载作用下开挖基槽附近海床动态响应和液化破坏情况，提出一个二维耦合计算模型，采用雷诺时均纳维-斯托克斯（RANS）方程描述波浪运动情况，通过设定侧边界条件实现稳定流场。海床部分通过求解Biot固结方程，得到波流荷载下海床中的应力和位移情况。将模型计算结果与水槽试验数据和解析解进行比较，验证了波流模型和海床模型的有效性。在此模型基础上，分析得到了开挖之后海床新的应力和固结状态。同时，通过参数分析得到了波流耦合情况下波浪形态的变化，以及海流对海床液化情况和孔压情况的影响。最后，通过线性回归计算得到最大液化深度与流速的拟合关系曲线。计算结果可用于判断基槽开挖后不规则海床的液化情况，对相关研究和实际工程具有一定参考意义。     

水动力时空变化对近岸风浪演化的影响——以渤海湾西南岸为例

采用SWAN模型和ADCIRC模型建立了风浪、潮汐和水流联合作用的耦合数值模式,并通过渤海湾西南岸实测资料对该模式进行了验证。利用该模式分析了近岸区水位和流场时空变化对风浪模拟结果的影响,计算结果表明水位变化对近岸区风浪模拟结果有显著影响,特别是中等大风过程高潮位时波高受水位影响的变化幅值可达0.5m以上,且水深越浅影响越大。但在岸滩平缓的近岸海域由于流速、流向的时空变化不太剧烈,流场作用和波浪辐射应力作用对波浪场的影响都基本可以忽略。在模拟近岸风浪过程时,应选用耦合模式。     

Numerical Modeling of Wave and Current Interaction in the Tidal Inlet Area

This paper investigates the phenomena of wave refraction and diffraction in the slowly varying topography, as well as the current deflection due to wave actions. A numerical model is developed based on depth integrated mean continuity equation and momentum equations, and a 3rd-order wave equation which governs the wave diffraction, refraction and interaction with current. Examples to examine the above model are given comparing with the laboratory data or the numerical results of other researchers. An example simulating the inlet area shows the interesting velocity field which may be used as a pioneer to further study on the nearshore hydrodynamics and sedimentation.     

非线性波传播的新型数值模拟模型及其实验验证 引入变换速度变量

以一种新型的含变换速度变量的Boussinesq型方程为控制方程组,采用五阶Runge-Kutta-England格式离散时间积分,采用七点差分格式离散空间导数,并采用恰当的出流边界条件,从而建立了非线性波传播的新型数值模拟模型.对均匀水深水域内波浪传播的数值模拟,说明在引入变换速度后进一步增大了模型的水深适用范围.对潜堤地形上波浪传播的数值模拟说明,在引入变换速度后进一步提高了模型的数值模拟精度.     

近岸沿岸流及污染物运动的数值模拟

基于双曲型缓坡方程和近岸浅水方程对近岸波浪斜向入射破碎所生成的沿岸流及污染物在沿岸波流作用下的运动进行了数值模拟,并对数值模拟结果进行了验证分析。数值模拟结果表明,在相近工况参数下,随着入射波高的增大,沿岸流流速和平均水面升高值均明显增大;随着岸坡坡度的增加,沿岸流流速和平均水面升高值明显增大;随着入射波浪周期的增大,平均水面升高值明显增大。在沿岸缓坡区域,由斜向入射波浪破碎所产生的沿岸流对污染物的运动起着不可忽略的影响。     

Modeling of propagation and transformation of transient nonlinear waves on a current

A novel theoretical approach is applied to predict the propagation and transformation of transient nonlinear waves on a current. The problem was solved by applying an eigenfunction expansion method and the derived semi-analytical solution was employed to study the transformation of wave profile and the evolution of wave spectrum arising from the nonlinear interactions of wave components in a wave train which may lead to the formation of very large waves. The results show that the propagation of wave trains is significantly affected by a current. A relatively small current may substantially affect wave train components and the wave train shape. This is observed for both opposing and following current. The results demonstrate that the application of the nonlinear model has a substantial effect on the shape of a wave spectrum. A train of originally linear and very narrow-banded waves changes its one-peak spectrum to a multi-peak one in a fairly short distance from an initial position. The discrepancies between the wave trains predicted by applying the linear and nonlinear models increase with the increasing wavelength and become significant in shallow water even for waves with low steepness. Laboratory experiments were conducted in a wave flume to verify theoretical results. The free-surface elevations recorded by a system of wave gauges are compared with the results provided by the nonlinear model. Additional verification was achieved by applying a Fourier analysis and comparing wave amplitude spectra obtained from theoretical results with experimental data. A reasonable agreement between theoretical results and experimental data is observed for both amplitudes and phases. The model predicts fairly well multi-peak spectra, including wave spectra with significant nonlinear wave components.     

Distribution of eddy scales for wave current combined flow

This paper illustrates the modulation of the eddy scale distribution due to superimposition of surface wave on only current flow. Time series data of three-dimensional velocity components were measured in a laboratory flume by a three-dimensional (3D) 16-MHz micro-acoustic Doppler velocimeter (Micro-ADV). The velocity time series of only current case and waves following the current were analysed to obtain the phase-averaged mean velocities, turbulent intensities, and Reynolds stress. The probability density function of phase-averaged stream-wise and vertical velocity fluctuations showed bimodal oscillations towards the free surface for higher frequency surface waves. It was revealed that surface waves along the current effectively decrease the intermittency of turbulence of the only current flow. Surface wave changed the intermittent structure of only current flow by modulation of the energy cascade mechanism of the only current flow by introduction of wave induced length scales. Also the scale of the finer dissipative eddies were prominently enhanced by the increase in surface wave frequency. Wavelet analysis of time series of velocity signals provided information on the eddy scale and their frequency of occurrence. It was found that the large eddies are carried by the crest regions of the progressive wave while the small scale eddies are carried by the trough regions.