内孤立波对小直径直立桩柱作用力的实验研究

内孤立波具有较大的振幅与较强的垂向剪切,能对海洋工程设施产生极大的破坏作用。本文设计实验研究了内孤立波与直立小直径桩柱的相互作用,采用粒子图像速度测量法(PIV)测量内孤立波的流速场,并采用自行设计的测力天平测量桩柱受力,测量分析了内孤立波对小直径直立桩柱产生作用力的实验值,与Morison公式计算的理论值比较。实验结果表明,第一模态内孤立波的流速方向以及作用力在桩柱的上下部分方向相反,产生很强的速度切变和扭力,对小直径直立桩柱造成破坏。通过与Morison公式计算的理论值比较,发现实验值与理论值的大小以及分布规律基本相同。     

沿岸流流速垂向分布的实验研究

进行了沿岸流模型实验,测量沿岸流流速的垂向分布.采用Faria等(1996)的方法对实验结果进行了分析,将水深分为上层和下层两部分:下层采用对数分布来表达沿岸流的垂向分布;上层考虑了波浪引起的自由表面的波动,使流速分布包括了修正的对数分布和质量输移速度两部分.对实验结果得到的相对粗糙度与sleath(1991)的经验公式结果进行了对比,也对摩阻流速结果与已有经验公式计算结果的进行了对比.     

两相流模型在南海内孤立波与沙波相互作用中的应用

为研究内孤立波与沙波的相互作用,本文对基于OpenFOAM的SedWaveFoam求解器进行改进,建立了内孤立波-泥沙运动欧拉两相流模型。在利用试验资料对模型进行验证的基础上,在南海北部典型代表性条件下,模拟分析了500 m水深位置沙波床面上内孤立波作用下的水动力变化和泥沙运动。结果表明,内孤立波逐渐离开沙波时,海底沙波背流面处出现与内孤立波背景流速反向的流速,在内孤立波导致的流场作用下,沙波床面上的泥沙悬起并运动到床面以上的水体中。振幅100 m的内孤立波可以导致床面以上14 m高的位置处出现约0.07 kg/m³的悬沙浓度。     

潮汐作用下渤海温跃层波动与起伏的数值研究

研究了潮汐(四大分潮 M₂,S₂,K₁,O₁同时输入)作用下渤海温跃层起伏与波动的三维数值模型(将海洋分为3层,即上混和层、跃层和下混和层),揭示了整个海区温跃层上界面处跃层起伏(在文中指每个时刻跃层波高的周期平均值)的地理分布及叠加在起伏之上的潮周期波动的时空变化,模拟出跃层波动与实测基本一致.结果发现大振幅的跃层波动均发生在海峡及近海地形突变之处.一般界面波动的波高大于甚至远远大于同一位置的表层潮波.从位相以及周期来看,潮波和跃层上下界面波动相互之间,有些海域一致,有些地方则相差甚远.跃层上界深度及厚度的梯度,对跃层起伏分布有一定的影响.跃层起伏还可能与海岸海底摩擦有关.     

Experimental study on internal waves in a stratified shear flow

This paper describes experiments on interfacial phenomena in a stratified shear flow having a sharp velocity shear at a density interface. The interface was visualized in vertical cross-section using dye, and the flow pattern was traced using aluminum powder. Two kinds of internal waves with different phase velocities and wave profiles were observed. They are here named p(positive)-waves and n(negative)-waves, respectively. By means of a two-dimensional visualization technique, the following facts have been confirmed regarding these waves. (1) The two kinds of waves propagate in the opposite direction relative to a system moving with the mean velocity at the interface, and their dispersion relations approximately agree with the two solutions of interfacial waves in a two-layer system of a linear basic shear flow. (2) The p-wave has sharp crests and flat troughs, and the n-wave has the reverse of this. This difference in wave profile is due to the finite amplitude effect. (3) Phase velocity of each wave lies within the range of the mean velocity profile, so that a critical layer exists and each wave has a “cat's eye” flow pattern in the vicinity of the critical layer, when observed in a system moving with the phase velocity. Consequently, these two waves are symmetrical with respect to the interface. The mechanisms of generation of these waves, and the entrainment process are discussed. It is inferred that when the “cat's eye” flow pattern is distorted and a stagnation point approaches the interface, entrainment in the form of a stretched wisp from the lower to the upper layer occurs for the p-wave, and from the upper to the lower layer for the n-wave.     

实验室造波条件对内孤立波发展影响的直接数值模拟

内孤立波具有振幅尺度大、能量集中的特点,其引起流场和密度场的迅速变化可能对海洋工程结构物以及水下潜体造成严重威胁。因此研究不同造波条件下生成的内孤立波运动的流场特征具有重要的学术意义和实际应用价值。采用直接数值模拟方法和给定的初始密度场密度跃迁函数,对重力塌陷激发内孤立波的运动过程进行研究,探讨了不同造波条件下,激发产生的内孤立波波型、涡度、振幅和水平速度等流场特征。结果表明:(1)直接模拟数值方法能够模拟内孤立波传播过程中的密度界面波型反转现象;(2)从定性和定量的角度,证实了不稳定内孤立波传播过程中存在能量的向后传递;(3)对于相同的台阶深度(水闸两侧初始密度界面的高度差),初始涡流保持相同,但是随着上下层水深比的减小,其强度下降显著;(4)台阶深度对初始涡流的垂直结构的影响要大于上下层水深比,且台阶深度对内孤立波的振幅、水平速度的影响显著。     

High-fidelity simulation of regular waves based on multi-moment finite volume formulation and THINC method

The performance of interFoam (a widely used solver within OpenFOAM package) in simulating the propagation of water waves has been reported to be sensitive to the temporal and spatial resolution. To facilitate more accurate simulations, a numerical wave tank is built based on a Navier–Stokes model, which employs the VPM (volume-average/point-value multi-moment) scheme as the fluid solver and the THINC/QQ method (THINC method with quadratic surface representation and Gaussian quadrature) for the free-surface capturing. Simulations of regular waves in an intermediate water depth are conducted and the results are assessed via comparing with the analytical solutions. The performance of the present model and interFoam solver in simulating the wave propagation is systematically compared in this work. The results clearly demonstrate that compared with interFoam solver, the present model significantly improves the dissipation properties of the propagating wave, where the waveforms as well as the velocity distribution can be substantially maintained while the waves propagating over long distances even with large time steps and coarse grids. It is also shown that the present model requires much less computation time to reach a given error level in comparison with interFoam solver.     

Flow under standing waves: Part 1. Shear stress distribution,energy flux and steady streaming

The conditions for energy flux, momentum flux and the resulting streaming velocity are analysed for standing waves formed in front of a fully reflecting wall. The exchange of energy between the outer wave motion and the near bed oscillatory boundary layer is considered, determining the horizontal energy flux inside and outside the boundary layer. The momentum balance, the mean shear stress and the resulting time averaged streaming velocities are determined. For a laminar bed boundary layer the analysis of the wave drift gives results similar to the original work of Longuet–Higgins from 1953. The work is extended to turbulent bed boundary layers by application of a numerical model. The similarities and differences between laminar and turbulent flow conditions are discussed, and quantitative results for the magnitude of the mean shear stress and drift velocity are presented. Full two-dimensional simulations of standing waves have also been made by application of a general purpose Navier–Stokes solver. The results agree well with those obtained by the boundary layer analysis. Wave reflection from a plane sloping wall is also investigated by using the same numerical model and by physical laboratory experiments. The phase shift of the reflected wave train is compared with theoretical and empirical models.     

相同尺度地形与背景流共振生成内孤立波机制

为研究内孤立波的地形和背景流共振机制,用地形和背景流共振机制计算了3个潜标观测的内孤立波(不同模态、不同波长)的流速和传播速度,并与观测到的内孤立波进行比较。潜标观测的第一模态内孤立波(波长分别为6.4和3.3km)都是下凹型内孤立波,2个内孤立波的传播速度约为1.4m/s、最大振幅约为48m,水平流向结构都是上层西北向、下层东南向,波长3.3km 的内孤立波波峰前后有更明显的下降流和上升流。用共振机制计算出的第一模态和第二模态纬向流速的垂向结构与观测相同,最大纬向流速出现的深度与观测一致,分别相差5和12m。用共振机制计算出的内孤立波传播速度与用 KdV 方程计算的传播速度相当,共振机制计算波速为0.66~1.21m/s,KdV 方程计算波速为0.79~1.40m/s。     

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

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

完全非线性深水波的数值模拟

基于势流理论,并结合深水波质点运动从水面向下呈e指数衰减的特性,建立了完全非线性数值变深水槽模型,通过实时模拟活塞式造波机运动来产生波浪.采用时域高阶边界元法进行模拟,利用混合欧拉-拉格朗日方法和四阶Runge-Kutta方法追踪流体瞬时水面,应用镜像格林函数消除了水槽两个侧面的积分,在水槽末端布置人工阻尼层来消除反射...     

Mud-Wave Interaction： A Viscoelastic Model

1 .IntroductionThe interaction between waves propagating onthe surface of a body of water andthe bed materialis a long standing coastal problem.The problemis practicallyimportant because the waves can be at-tenuated at a muchfaster rate whenthe bed materi…     

垂直分辨率对长江口海域M2分潮模拟的影响

基于EFDC(Environmental Fluid Dynamics Code)模式建立了长江口及其邻近海域的三维水动力学模型, 研究模型的垂直分辨率对该海域M₂分潮模拟的影响。结果表明:垂直分辨率的变化对M₂分潮传播方向的模拟结果影响较小, 但其可通过底摩擦和湍流耗散两个计算过程来影响潮能通量的模拟结果, 最终对长江口和杭州湾内的M₂分潮振幅产生显著的影响。最底层厚度较大时, 上层自由水体的高流速特征在最底层过于明显, 进而导致计算的底摩擦应力偏高, 此时提高底层的垂直分辨率会降低底摩擦对能量的耗散。另一方面, 垂直湍流混合作用会随垂直分辨率的增加而增强, 所以垂直分辨率增加到一定程度后, 上层自由水体的高流速会经由增强的湍流混合而更多的传入底层, 使计算的底摩擦应力随垂直分辨率的提高而有重新增加的趋势, 进而又增强底摩擦对潮能的耗散。     

Mass transport in a two-layer wave boundary layer

The transport of a chemical species under the pure action of surface progressive waves in the benthic boundary layer which is loaded with dense suspended sediments is studied theoretically. The flow structure of the boundary layer is approximated by that of a two-layer Stokes boundary layer with a sharp interface between clear water and a heavy fluid. The simplest model of constant eddy diffusivities is adopted and the exchange of matter with the bed is ignored. For a thin layer of heavy fluid, whose thickness is comparable to the surface wave amplitude and the Stokes boundary layer thickness, effective transport equations are deduced using an averaging technique based on the method of homogenization. The effective advection velocity is found to be equal to the depth-averaged mass transport velocity, while the dispersion coefficient can be shown to be positive definite. Explicit expressions for the transport coefficients are obtained as functions of fluid properties and flow kinematics. Physical discussions on their relations are also presented.     

Study of nonlinear currents induced by coincident surface waves in liquid

This is an experimental study of the mixing induced by coincident surface waves in a liquid. The main mechanism leading to the emergence of mixing was shown to be the middle currents generated by coincident waves. The regime of these currents strongly depends on the amplitude of surface waves. For waves of near-critical amplitudes, an intense turbulization of middle currents is observed. Patterns of the velocity field were obtained using the Particle Imaging Velocimetry (PIV) technique for different amplitudes of surface waves. The results obtained can be used to estimate mixing in the near-surface oceanic layer.     

Flow Separation and Vortex Dynamics in Waves Propagating over A Submerged Quartercircular Breakwater

The interactions of cnoidal waves with a submerged quartercircular breakwater are investigated by a Reynolds-Averaged Navier–Stokes (RANS) flow solver with a Volume of Fluid (VOF) surface capturing scheme (RANS-VOF) model. The vertical variation of the instantaneous velocity indicates that flow separation occurs at the boundary layer near the breakwater. The temporal evolution of the velocity and vorticity fields demonstrates vortex generation and shedding around the submerged quartercircular breakwater due to the flow separation. An empirical relationship between the vortex intensity and a few hydrodynamic parameters is proposed based on parametric analysis. In addition, the instantaneous and time-averaged vorticity fields reveal a pair of vortices of opposite signs at the breakwater which are expected to have significant effect on sediment entrainment, suspension, and transportation, therefore, scour on the leeside of the breakwater.     

浮式生产储油系统内孤立波载荷特性实验

将置于大尺度密度分层水槽上下层流体中的两块垂直板反方向平推,以基于MCC理论解的内孤立波诱导上下层流体中的层平均水平速度作为其运动速度,发展了一种振幅可控的双推板内孤立波实验室造波方法,并对内孤立波作用下浮式生产储油系统(FPSO)的载荷特性开展了系列实验。结果表明,无因次内孤立波水平力和力矩幅值均与先导内孤立子无因次振幅之间呈线性关系,其斜率与上下层流体深度比有关;在各种上下层流体深度比下,无因次内孤立波垂向力幅值与先导内孤立子振幅之间近似呈幂函数关系。     

浮冰在规则波和双色波下的运动特性研究

为研究浮冰在波浪下的运动特性,以单个聚氯乙烯（PVC）塑料圆板模拟浮冰,进行了一系列波浪试验研究。在规则波试验基础上,同时研究了双色波条件下的浮冰运动特性。结果显示：规则波条件下浮冰运动幅频响应算子首先随无量纲波长的增大而增大,当无量纲波长大于4.0后,运动幅频响应算子随无量纲波长的变化不明显,且趋于一定值;在双色波条件下,对应频率组成成分的浮冰运动幅频响应算子与规则波条件下随波长的变化规律一致。根据部分规则波试验结果提出预测浮冰慢漂速度的经验公式,并用已有的试验和余下的规则波与双色波试验结果进行验证。结果表明,经验公式对规则波和双色波条件下的浮冰慢漂速度的误差在20%以内,预测结果吻合较好。     

波浪与起伏水平板防波堤相互作用数值模拟

利用自主研发的基于紧致插值曲线CIP(constrained interpolation profile)方法的数学模型,开展规则波与起伏水平板防波堤相互作用的数值模拟研究。模型在笛卡尔直角坐标下建立,以CIP方法为流场基本求解器,分步求解Navier-Stokes方程,利用高精度的流体体积类型的THINC/SW (tangent of hyperbola for interface capturing with slope weighting)方法重构自由液面,采用浸入边界IBM(immersed boundary method)方法处理波浪与起伏板防波堤的耦合作用问题,通过动量源项造波方法模拟波浪的产生。重点关注波浪的浅水变形和板两端涡旋脱落的非线性现象,分析不同潜深、波要素下的板周围流场分布、板的运动响应和波浪的反透射系数。结果表明:起伏水平板主要通过能量反射、板上浅水变形和板两端的涡脱落消能,能有效减小板后波高,具有作为防波堤的可行性。