《海洋工程》2007（第25卷）,第1～4期,总目录

第25卷,第1期,2007年2月分层流体中内孤立波在潜浮式竖直薄板上透射和反射…………………………………………魏岗,尤云祥,缪国平,等(1)深水导管架下水扶正过程定量风险分析……………………………………………………………余建星,谭振东,田佳(9)导管架下水过程中浮力及浮心精确     

内孤立波对不同水深竖直圆柱体水平作用力分析

采用理论研究和试验研究相结合的方法,对两层流体中内孤立波对不同水深竖直圆柱体水平作用力进行了分析。采用重力塌陷法进行内孤立波造波试验,通过PIV(Particle Image Velocimetry)技术获取内孤立波波剖面的时间序列,利用测力系统对竖直圆柱体在内孤立波作用下的受力进行测量。通过内孤立波波剖面试验分析,选取合适的内孤立波理论并对竖直圆柱体所受内孤立波作用力进行计算,将计算结果与试验结果进行汇总分析。结果表明:上层流体中圆柱体所受内孤立波最大水平作用力明显大于下层流体中各段圆柱体,且方向相反;随着内孤立波振幅的增大,各段圆柱体所受内孤立波最大水平作用力也相应地增大;密度跃层处圆柱体所受水平作用力具有相较于其他位置处更为复杂的变化趋势;选用合适的内孤立波理论进行计算能够较好的模拟竖直圆柱体的实际水平受力。     

基于MCC理论的内孤立波数值模拟

基于MCC(Miyata-Choi-Camassa)理论,将内孤立波诱导上下层深度平均水平速度为入口条件,结合理想流体完全非线性欧拉方程,建立了两层流体中内孤立波生成与传播的CFD(Computational Fluid Dynamics)数值水槽方法。在系列数值模拟基础上,获得了内孤立波设计振幅与数值模拟振幅之间的相关关系,实现了在振幅可控条件下的内孤立波数值模拟。结果表明,在极限振幅范围内,数值模拟所得内孤立波波形均与MCC理论解吻合良好,而KdV(Korteweg-de Vries)和mKdV(modi-fied KdV)理论解只适用于小振幅的情况。同时,利用CFD数值模拟结果,对内孤立波诱导流速场特性进行了分析,结果表明内孤立波诱导水平速度在上层及波面下方流体层中垂向衰减很小,但在内界面与波面之间流体层中垂向衰减明显。     

内孤立波过山脊地形演化实验研究

通过模型实验,研究了下沉型内孤立波通过山脊地形演化特征。实验以三角形障碍物模拟海底山脊地形,采用两种密度的分层水,对上层流体和下层流体的高度比、密度等进行了量化处理。实验研究表明:KdV理论波形可较好模拟本次实验内孤立波波形,但随着内孤立波振幅的增大,误差增加;在内孤立波与障碍物微量作用、中度作用和破波作用三种程度的相互作用中,内孤立波过障碍物具有不同的波形变化和主波能量衰减率。     

二维浮体上孤立波透射与反射的解析解

本文应用渐近匹配法得到了浅水孤立波在二维薄浮体上透射与反射的解析解,给出了透射波波高和相位变化的规律以及反射波的特性,并探讨了浮体与水底间间隙对孤立波透射与反射的影响。在工程上对淡水域浮式防波堤设计有借鉴和指导意义。     

海表外强迫效应对连续层化流体中先驱内孤立波生成的影响

移动的海面或海底强迫在共振条件下会向强迫上游周期性地生成内孤立波,这种内孤立波被称为先驱内孤立波。本文利用三维全非线性、全频散的非静力模型NHWAVE,模拟连续分层流体受到海面外强迫作用而产生的先驱内孤立波,研究了正负外强迫效应对先驱内孤立波生成过程和形态特征的影响。结果显示,负效应强迫作用下先驱内孤立波具有更长的生成周期,这与两层流体中的结论相同;在强迫区内部,正效应强迫使水体的等密度面维持较稳定的结构,而当强迫效应为负时,强迫区内水体的等密度面随时间起伏较大。在正负两种效应强迫作用下,流场均维持了逆时针的回流结构,不过在负效应强迫作用下,流场更为复杂。     

两层流体中内孤立波与潜体相互作用数值模拟

提出基于双推板的内孤立波数值造波方法,对两层流体中内孤立波与潜体的相互作用进行了数值模拟,分析不同潜深下潜体所受的内孤立波载荷特性.结果表明,内孤立波对水下潜体的载荷作用是不可忽视的.特别地,当潜体位于内孤立波中时,其所受的垂向力要远比潜体位于内孤立波下时的大.     

内孤立波与Spar型浮式风机平台耦合作用研究

针对内孤立波与系泊浮式平台的耦合作用问题,本文以单立柱Spar型浮式风机平台为研究对象,研究了平台在内孤立波作用下的耦合运动特性。利用重叠网格技术及自编耦合动力程序,以计算流体动力学求解Navier-Stokes流体控制方程,建立了Spar型浮式风机平台与内孤立波(Internal solitary waves, ISWs)流固耦合作用数值模型。通过系列数值模拟,研究了Spar型浮式风机平台在内孤立波作用下的水平动态载荷和纵荡、横荡规律。研究表明,内孤立波对Spar平台会产生较显著的水平作用力,并产生大幅的水平位移;在水深一定时,同一流体分层比下,随着内孤立波振幅提高,纵荡运动幅值增大;当振幅一定时,随着流体分层比的增大,纵荡幅值减小;通过调整系泊系统,发现非线性系泊系统下平台纵荡幅值要大于线性系泊系统下的情况,即线性系泊系统会低估纵荡运动幅值。     

内孤立波过凹陷地形的实验探究

本文通过实验室水槽实验,讨论内孤立波经过凹陷地形的演化过程。实验在两层流体中进行,上下水层的厚度和密度被定量地控制,改变凹陷地形的宽度,并对每种地形分别进行了不同分层结构下具有不同初始振幅的内孤立波对照实验,研究内孤立波与凹陷地形的作用过程。研究发现,本文实验环境下凹型内孤立波实验波形与mKdV(Korteweg-deVries)理论波形更符合;而上凸型内孤立波实验,当非线性参数ε≤0.22时,KdV理论波形与实验波形符合较好,当非线性参数ε≥0.27时.mKdV理论与实验波形符合较好。下凹型内孤立波经凹陷地形过程中波形变化轻微,主波振幅有先减小后增大的趋势,且障碍比越大,变化趋势越明显。本实验结果可作为研究海洋中内孤立波与海底凹陷地形作用情况的参考。     

有限深两层流体内孤立波非线性数值模拟入口边界条件研究

基于有限深两层流体KdV(Korteweg-de Vries)、eKdV(extended KdV)和MCC(Miyata-Choi-Camassa)理论,以内孤立波诱导上下层深度平均水平速度为入口边界条件,采用理想流体完全非线性欧拉方程,建立了两层流体中内孤立波生成的CFD(Computational Fluid Dynamics)数值模拟方法。以系列数值模拟结果为依据,结合内孤立波非线性和色散参数的组合条件,给出了选择合适内孤立波理论解作为CFD数值模拟入口边界条件的方法,从而实现了振幅与波形可控的内孤立波完全非线性数值模拟。     

Experimental Study on Ocean Internal Wave Force on Vertical Cylinders in Different Depths

A series of experimental studies about the force of internal solitary wave and internal periodic wave on vertical cylinders have been carried out in a two-dimensional layered internal wave flume. The internal solitary waves are produced by means of gravitational collapse at the layer thickness ratio of 0.2, and the internal periodic waves are produced with rocker-flap wave maker at the layer thickness ratio of 0.93. The wave parameters are obtained through dyeing photography. The vertical cylinders of the same size are arranged in different depths. The horizontal force on each cylinder is measured and the vertical distribution rules are researched. The internal wave heights are changed to study the impact of wave heights on the force. The results show that the horizontal force of concave type internal solitary wave on vertical cylinder in the upper-layer fluid has the same direction as the wave propagating, while it has an opposite direction in the lower-layer. The horizontal force is not evenly distributed in the lower fluid. And the force at different depths increases along with wave height. Internal solitary wave can produce an impact load on the entire pile. The horizontal force of internal periodic waves on the vertical cylinders is periodically changed at the frequency of waves. The direction of the force is opposite in the upper and lower layers, and the value is close. In the upper layer except the depth close to the interface, the force is evenly distributed; but it tends to decrease with the deeper depth in the lower layer. A periodic shear load can be produced on the entire pile by internal periodic waves, and it may cause fatigue damage to structures.     

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

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

Two-Layer Non-Hydrostatic Model for Generation and Propagation of Interfacial Waves

When pycnocline thickness of ocean density is relatively small, density stratification can be well represented as a two-layer system. In this article, a depth integrated model of the two-layer fluid with constant density is considered,and a variant of the edge-based non-hydrostatic numerical scheme is formulated. The resulting scheme is very efficient since it resolves the vertical fluid depth only in two layers. Despite using just two layers, the numerical dispersion is shown to agree with the analytical dispersion curves over a wide range of kd, where k is the wave number and d the water depth. The scheme was tested by simulating an interfacial solitary wave propagating over a flat bottom, as well as over a bottom step. On a laboratory scale, the formation of an interfacial wave is simulated,which also shows the interaction of wave with a triangular bathymetry. Then, a case study using the Lombok Strait topography is discussed, and the results show the development of an interfacial wave due to a strong current passing through a sill.     

Green functions with pulsating sources in a two-layer fluid of finite depth

The derivation of Green function in a two-layer fluid model has been treated in different ways.In a two-layer fluid with the upper layer having a free surface,there exist two modes of waves propagating due to the free surface and the interface.This paper is concerned with the derivation of Green functions in the three dimensional case of a stationary source oscillating.The source point is located either in the upper or lower part of a two-layer fluid of finite depth.The derivation is carried out by the method of singularities.This method has an advantage in that it involves representing the potential as a sum of singularities or multipoles placed within any structures being present.Furthermore,experience shows that the systems of equations resulted from using a singularity method possess excellent convergence characteristics and only a few equations are needed to obtain accurate numerical results.Validation is done by showing that the derived two-layer Green function can be reduced to that of a single layer of finite depth or that the upper Green function coincides with that of the lower,for each case.The effect of the density on the internal waves is demonstrated.Also,it is shown how the surface and internal wave amplitudes are compared for both the wave modes.The fluid in this case is considered to be inviscid and incompressible and the flow is irrotational.     

张力腿平台内孤立波作用特性数值模拟

依据三类内孤立波理论KdV、eK dV和MCC的适用性条件,采用Navier-Stokes方程为流场控制方程,以内孤立波诱导上下层深度平均水平速度作为入口边界条件,建立了两层流体中内孤立波对张力腿平台强非线性作用的数值模拟方法。结果表明,数值模拟所得内孤立波波形及其振幅与相应理论和实验结果一致,并且在内孤立波作用下张力腿平台水平力、垂向力及力矩数值模拟结果与实验结果吻合。研究同时表明,张力腿平台内孤立波载荷由波浪压差力、粘性压差力和摩擦力构成,其中摩擦力很小,可以忽略;水平力的主要成分为波浪压差力和粘性压差力,粘性压差力与波浪压差力相比较小却不可忽略,流体粘性的影响较小;垂向力中粘性压差力很小,流体粘性影响可以忽略。     

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

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

Hydroelastic interaction between obliquely incident waves and a semi-infinite elastic plate on a two-layer fluid

The hydroelastic response of a semi-infinite thin elastic plate floating on a two-layer fluid of finite depth due to obliquely incident waves is investigated. The upper and lower fluids with different densities separated by a sharp and stable interface are assumed to be inviscid and incompressible and the motion to be irrotational. Simply time-harmonic incident waves of the surface and interfacial wave modes with a given angular frequency are considered within the framework of linear potential flow theory. With the aid of the methods of matched eigenfunction expansion and the inner product of the two-layer fluid, a closed system of simultaneous linear equations is derived for the reflection and transmission coefficients of the series solutions. Based on the dispersion relations for the gravity waves and the flexural–gravity waves in a two-layer fluid and Snell’s law for refraction, we obtain a critical angle for the incident waves of the surface wave mode and three critical angles for the incident waves of the interfacial wave mode, which are related to the existence of the propagating waves. Graphical representations of the series solutions show the interaction between the water waves and the plate. The effects of several physical parameters, including the density and depth ratios of the fluid and the thickness of the plate, on the wave scattering and the hydroelastic response of the plate are studied. It is found that the variation of the thickness of the plate may change the wave numbers and the critical angles. The density ratio is the main factor to influence the wave numbers of the interfacial wave modes. Finally, the stress state is considered.     

Laboratory modeling of the propagation of periodic internal waves over bottom slopes

The experimental investigation of the run-up of periodic internal waves in a two-layer fluid on the coastal slope is performed in an open hydrochannel at the Physical Department of the Lomonosov Moscow State University. The waves are produced by a wave generator. We study the transformation of waves, the vertical structure of the field of velocities of mass transfer, and the behavior of the parameters of internal waves propagating over the sloping bottom. It is shown that the run-up and breaking of internal waves are accompanied by periodic emissions of portions of the heavier fluid from the bottom layer upward along the slope. The Stokes drift velocity changes its sign as a function of depth. Moreover, both the wave length (the horizontal distance between the neighboring crests) and the height of waves over the sloping bottom (the elevation of the crest over the slope along the vertical) decrease as the wave approaches the coast.     

Numerical simulation of interaction between internal solitary waves and submerged ridges

Direct numerical simulations are performed to study the transformation of internal solitary waves (ISWs) of depression type propagating over an underwater ridge in a two-layer fluid system. Bottom ridges with relatively smooth vertex are employed to represent sills in natural lakes and oceans. Consistent with previous experiments, three interaction types (weak, moderate and strong) are observed to be based on the energy loss. In addition, the moderate interaction are found to be categorized into transmitted and reflected type according to their distinct transformation process. General flow characteristics for ISW–ridge interaction in the benthic boundary layer and in the pycnocline is monitored and analysed. A modified degree of blocking B_m considering both the nonlinear effect of incident ISWs and the blockage effect of the submerged ridge is proposed. Different ISW–ridge interactions are discovered to be linked with B_m. Maximum wave-induced velocities, wave energy losses, reflected and transmitted wave amplitudes are found to have a self-similar feature with B_m. The maximum energy loss is up to 35% and the maximum wave-induced velocity can reach 1.8 times of the phase speed of the incident ISW. Empirical equations are obtained based on the data fitting to predict some useful physical parameters during ISW–ridge interaction.