分层海洋中运动物体生成的内波

该文提供分层海洋中运动物体生成内波的理论、数值和实验研究的背景资料的综述 ,并给出理论研究的 Green函数、正交模态模式和极点回避等方法的评述。数值研究给出数值求解Wronsky行列式的源求解方法及介绍直接差分计算含源内波方程和积分方程的方法。有关源致内波的实验研究 ,作者介绍物体在分层流水槽运动时生成内波的两个典型实验及其主要实验结果。最后评述船舶在海洋表面处生成表面波的研究方法 ,以此指导水下 Rankine卵形体生成内波的实验研究。     

基于快速行进算法的南海内波传播模型

南海内波水平二维空间位置的估计和预测具有应用价值。一个简单而精确的估计南海非线性内孤立波地理位置的经验模型已由Jackson(2009)提出,该模型将快速行进算法应用于海洋内波传播并给出内波相速度与深度关系的经验模式函数。本文尝试用KdV方程计算非线性内波相速度,以探讨海水层化对相速度和内波传播的影响。初步结果表明,海水层化的时空差异在内波传播过程中不可忽略。本文进一步将考虑海水层化的非线性内波相速度代替文献中的经验模式函数,建立不依赖于现场实测数据的基于快速行进算法的南海内波传播模型,并给出计算个例。本文模型的计算结果从个例来看与文献一致。     

单一海脊地形对海洋内波生成与传播影响的分析

在海洋内波多发区,海底地形变化是影响海洋内波生成、传播和演化的重要因素。本文基于不可压缩原始N-S方程,在非静压近似条件下,通过建立适应于非线性海洋内波研究的非静压海洋动力学模型,并将其应用于正压潮驱动下的内孤立波生成和传播的数值模拟研究。根据模拟结果,研究了一类单海脊地形拓扑结构变化对内孤立波生成和传播的影响;分析并讨论了地形拓扑结构参数变化与生成的内孤立波传播至特定位置的抵达时间、强度等特征参数之间的变化关系;提出内孤立波生成之前在海脊一侧形成"L-下陷"结构的观点,并揭示了与该观点相合的能量"积聚"和"释放"机制。     

海洋内孤立波预警监测识别技术及其在流花16-2油田群开发中的应用

我国南海内孤立波发生频繁,由于其振幅大、流速强、能量集中等特点,已成为深海油气开发海上安装工程必须考虑的重大风险之一。依托现场守护船,提出一种内波流结构单体式监测、识别技术,采用KdV方程实现内波流预警,并成功应用于流花16-2油田群开发项目海管管线终端(PLET)安装、浮式生产储卸油装置(FPSO)锚系及水下管缆回接、空气潜水等高风险作业的内波流预警安全保障中,同时基于此次内波流预警监测结果对该海域的内波生成源和生成机制进行了探讨。应用期间采用该技术共监测、识别并发布内波预警450次。应用结果表明,约93%的内波预警时刻误差在10 min以内,平均误差为±3.90 min。统计分析结果表明,该海域的内孤立波主要来源于巴士海峡而非局地的潮地相互作用,大致可以分为两类:一类是直接产生于巴士海峡,西传至该海域;另一类是局地产生,主要是由巴士海峡产生的内波西传至陆坡浅化分裂生成或是由巴士海峡产生的内潮西传至陆坡由于非线性变陡机制产生。该内波流监测技术在流花16-2油田群开发项目的成功应用,可为南海其他类似项目提供直接的借鉴。     

失稳背景流下对称型海洋内波的生成演变

本文采用无海底地形但考虑海洋跃层和剪切背景流的二维非静力准不可压缩方程组的数值模式,开展失稳垂向剪切背景流下线性和非线性对称型海洋内波生成演变的数值实验,并对结果进行分析、比较和讨论。研究结果表明,线性内波强度随积分时间始终呈指数增长,并有内波的对称不稳定;而非线性内波强度则在发展期呈准线性增长,最终进入稳定期。线性增长比非线性增长要快得多,非线性效应具有维稳作用。对该线性和非线性对称型内波,在跃层附近位密度扰动均有大值中心,即其为跃层所俘获,这与实际观测相一致;流函数与位密度扰动两者均有很好配合,位密度扰动的正、负中心分别相应于流函数的上升、下沉运动,表明有从海底向上的斜对流发生,且以跃层为顶盖。对线性内波来说,随积分时间增加,其波形大体不变,其正、负振幅也大体相同,并有符号相反原地增长的两个倾斜环流圈,而在它们之间则有较强倾斜上升流。非线性内波波形随积分时间改变,倾斜环流圈数目也在增加,最终形成负环流强于正环流的结果,并导致流函数、位密度扰动水平梯度剧增,其可视为间断。     

A Possible Mechanism of Ocean Finestructure Part Ⅰ: Energy and Coherence

依据在作者的前文（１９９８）中对于考虑旋转向量水平分量（ｆ～）的大洋内波方程求出的渐近解,文中对Ｓ场（考虑ｆ～条件下的总的内波场）,ＳＩ场（不计ｆ～影响的内波场）和ＳＦ场（Ｓ场中与ｆ～相关联的部分）导出了能量谱和相干谱。此理论谱与Ｍｕｌｅｒ等（１９７８）的ＩＷＥＸ实验的观测结果相比较表明,ＳＦ场的水平动能谱和铅直位移谱基本上与水平流细结构和铅直位移细结构的观测谱一致。作者推测ＳＦ场对应于中纬度大洋可逆细结构。但是对于相干谱的观测结果,在不考虑内波非线性相互作用和破碎过程的条件下,仅考虑ＳＦ场的影响还不能给予满意的解释。     

海南岛南部海区非线性内波特征分析

2004年6月在南海海南岛南部海域进行了内波观测实验,实验中使用了3条温度链组成链阵,来估计内波在该海区传播的方向和速度。实验中发现了非线性内波群的活动。对潮汐的测量表明,非线性内波群出现在海区的涨潮时刻。计算得到实验海区内波(非线性波群)速度约为0.54m/s。测量有效数据从29日10:00~30日14:00,只发现了一组非线性内波群,说明该海区的非线性内波群不同于别的海区所观测到具有半日潮周期的非线性内波。对非线性内波波形分析可以看出,KdV(Korteweg-de Vries)方程的dn2(x,t)解能比较好的描述得到的内波波形。数据表明非线性内波群在爬坡传播过程中随波包范围的不断扩张,波峰间距离有减小的趋势。对高频线性内波(非线性内波出现前后)的分析表明,高频线性内波中周期为15min~12min的波动成份的能量比较突出,这与非线性内波群中单个波包周期基本吻合,说明非线性内波群在传播扩张过程中其衰落的尾迹分散到了两个非线性波包之间的水体中。     

基于遥感与现场观测数据的南海北部内波传播速度

南海北部是全球海洋中内波最为活跃、生成和演变机制较为复杂的海域,本文利用多源卫星遥感数据(MODIS、GF-1、ENVISAT ASAR、RADARSAT-2)和现场观测数据开展了南海北部内波传播速度的研究。通过匹配捕获同一条内波的相邻两幅遥感图像,由内波的空间位移和时间间隔反演传播速度,并以0.5°×0.5°网格给出了南海北部内波传播速度的分布图。研究结果表明,内波传播速度受背景流场、水体层结和底地形变化等多因素影响,特别是水深。在南海北部由东至西、由南至北方向,内波传播速度逐渐递减。深海区内波传播速度最大,可达3m/s以上;内波在向西大陆架传播过程中,随着水深变浅速度逐渐减慢,传播速度为1—2m/s;大陆架浅海的内波传播速度较小,仅为零点几米每秒。同时,利用Kd V方程反演了内波传播速度理论值,对遥感数据提取的内波传播速度进行了精度验证,结果较为一致。     

基于多源遥感数据的日本海内波特征研究

日本海特殊的地理位置和复杂的地形使得该海域内波表征极为复杂,遥感是大范围观测内波的有效手段,已被广泛应用于内波的探测研究。本文利用MODIS、GF-1和ENVISAT ASAR遥感影像,开展了日本海内波特征研究。通过提取内波波峰线,生成了日本海内波空间分布图;获取了内波的波峰线长度和传播速度,并基于非线性薛定谔方程反演了内波振幅。研究结果表明,日本海内波分布范围宽广,不仅大陆架沿海区内波分布密集,深海盆地也探测到了大量内波;日本海北部45°N附近海域有少量内波出现,利用高分影像探测到朝鲜陆架浅海区有大量小尺度内波,大和海盆、大和隆起的西南部海域没有发现内波。日本海内波波峰线长达100多千米,深海区的传播速度大于1 m/s;浅海区内波振幅约10 m左右,深海区可达60 m以上。     

分层流水槽箱内垂荡板式内波造波模拟研究

提出了一种新型内波造波方法——箱内垂荡板式造波.应用FLUENT软件的有限体积法,采用k-ε两方程湍流模型封闭RANS方程作为求解控制方程,建立二维数值内波水槽.通过模仿物理造波机运动的数值模拟结果比较表明:箱内垂荡板式造波机生成的内波与理论解十分吻合,与双板联动式造波和平板拍击式造波相比,控制简单,对自由水面的扰动小...     

一种新型二阶全非线性Boussinesq方程及其数值验证

通过改进二阶全非线性 Boussinesq 波浪方程中的色散项,得到了一组没有改变原方程的数学形式但适用于更大变化水深的新方程,其色散性能和变浅性能都比原方程有了很大改进,所适用的水深范围更大,能更好地描述从深水到近岸浅水处的波浪传播;并基于新方程建立了波浪数值模型,通过模拟波浪从浅水到深水的传播变形来验证新方程的有效性.     

波浪沿斜面传播过程的数值模拟

精确模拟非线性波沿斜面传播过程非常困难,为此论文从势函数的边界积分方程出发,建立了一种时域内二维波浪模拟的数值模型,主要用来模拟完全非线性波浪的传播变形过程。论文的数值模型使用高阶二维边界元方法,采用可调节时间步长的基于二阶显式泰勒展开的混合欧拉-拉格郎日时间步进来求解带自由表面的线性或完全非线性波浪传播问题。在计算区域一端造出线性或非线性的周期性波浪,另一端采用消除反射波的人工粘性吸收边界。通过与现有理论比较证明了论文数值方法所得结果是准确可靠的。     

Gerris在海洋大振幅内孤立波数值模拟中的应用

本文运用基于自适应网格的流体动力学开源软件Gerris,来建立基于Boussinesq近似下的二维不可压缩Euler方程组的数值模型,以模拟不同层化条件下稳定状态的完全非线性大振幅内孤立波。文中比较了完全非线性的用Gerris实现的Euler模型与弱非线性的KdV理论模型在刻画大振幅内孤立波结构及特征参数上的差异,说明在模拟大振幅内孤立波时,高阶非线性不应忽略。Euler模型模拟结果表明,完全非线性大振幅内孤立波的等密度面半宽度随深度变化,这使得基于KdV方程解析解、利用卫星SAR(Synthetic Aperture Radar)图像提取内孤立波极值间距来反演内波振幅的可行性存疑,需要重新评估。此外,本文用两组实测数据验证了用Gerris实现的Euler模型模拟大振幅内波的有效性。     

THE NONLINEAR INTERNAL GRAVITY WAVES IN STRATIFIED FLUID

In this paper, starting from the equations of the nonlinear internal gravity waves in stratified fluid, using the method of the Taylor expansion nearby the equilibrium point for the nonlinear terms, we find the analytical solutions for nonlinear internal gravity waves. The linear internal gravity waves and solitary waves are its special cases. The nonlinear internal gravity waves satisfy the well-known KdV (Karteweg-de Vries) equation. The nonlinear internal gravity waves are different from linear waves in character. The former dispersive relation contains the amplitude, but the latter does not. The larger the amplitude and the wavelength the faster are waves for the nonlinear internal gravity waves. The smaller the stability of the stratification, the larger is the wavelength (or the width). Some phenomena such as squall line, cumulus, turbulent mass structure in atmosphere, and thermocline in ocean have these natures.     

三阶Stokes型随机波浪载荷谱研究

应用非线性谱分析理论，对三阶Stokes型随机波浪载荷谱进行了分析，将波面方程及海水质点的水平速度用一阶波面的非线性组表示，导出了随机波浪谱的表达式。为了便于求解随机波浪的载荷谱，将阻力项展开为幂级数式，并应用非线性谱分析理论，确定了幂级数的系数，进而将波浪载荷表示为一阶波面及其导数的非线性组合，最后得出波浪载荷谱密度的表达式，并应用数值分析方法，得出单位桩柱波浪力及总波浪力谱密度。     

Instability of Two Interacting Quasi-Monochromatic Waves in Shallow Water

The nonlinear interactions of waves with a double-peaked power spectrum have been studied in shallow water.The starting point is the prototypical equation for nonlinear unidirectional waves in shallow water,i.e.the Korteweg de Vries equation.By means of a multiple-scale technique two defocusing coupled Nonlinear Schrdinger equations are derived.It is found analytically that plane wave solutions of such a system are unstable for small perturbations,showing that the existence of a new energy exchange mechanism which can influence the behavior of ocean waves in shallow water.     

A Modified Form of Mild-Slope Equation with Weakly Nonlinear Effect

Nonlinear effect is of importance to waves propagating from deep water to shallow water.Thenon-linearity of waves is widely discussed due to its high precision in application.But there are still someproblems in dealing with the nonlinear waves in practice.In this paper,a modified form of mild-slope equa-tion with weakly nonlinear effect is derived by use of the nonlinear dispersion relation and the steady mild-slope equation containing energy dissipation.The modified form of mild-slope equation is convenient to solvenonlinear effect of waves.The model is tested against the laboratory measurement for the case of a submergedelliptical shoal on a slope beach given by Berkhoff et al,The present numerical results are also comparedwith those obtained through linear wave theory.Better agreement is obtained as the modified mild-slope e-quation is employed.And the modified mild-slope equation can reasonably simulate the weakly nonlinear ef-fect of wave propagation from deep water to coast.     

Resonant interaction of waves of continuous and discrete spectra in the simplest model of a stratified shear flow

The equations of dynamics of eddy—wave disturbances of two-dimensional stratified flows in an ideal incompressible fluid that are written in a Hamiltonian form are used to study the resonant interaction of waves of discrete and continuous spectra. A gravity—shear wave generated at a jump of the density and vorticity of the undisturbed flow and a wave generated at a weak vorticity jump, which is similar to a wave of a continuous spectrum, participate in the interaction. The equations are written in terms of normal variables to obtain the system of evolution equations for the amplitudes of the interacting waves. The stability condition for eddy—wave disturbances is derived within the framework of the linear theory. It is shown that a cubic nonlinearity may lead to the stabilization of unstable disturbances if the coefficient of the nonlinear term is positive.     

Deterministic and stochastic evolution equations for fully dispersive and weakly nonlinear waves

This paper presents a new and more accurate set of deterministic evolution equations for the propagation of fully dispersive, weakly nonlinear, irregular, multidirectional waves. The equations are derived directly from the Laplace equation with leading order nonlinearity in the surface boundary conditions. It is demonstrated that previous fully dispersive formulations from the literature have used an inconsistent linear relation between the velocity potential and the surface elevation. As a consequence these formulations are accurate only in shallow water, while nonlinear transfer of energy is significantly underestimated for larger wave numbers. In the present work we correct this inconsistency. In addition to the improved deterministic formulation, we present improved stochastic evolution equations in terms of the energy spectrum and the bispectrum for multidirectional waves. The deterministic and stochastic formulations are solved numerically for the case of cross shore motion of unidirectional waves and the results are verified against laboratory data for wave propagation over submerged bars and over a plane slope. Outside the surf zone the two model predictions are generally in good agreement with the measurements, and it is found that the accuracy of e.g., the energy spectrum and of the third-order statistics is considerably improved by the new formulations, particularly outside the shallow-water range.     

An Irregular Wave Generating Approach Based on naoe-FOAM-SJTU Solver

In this paper, a wave generating approach for long-crest irregular waves in a numerical tank by our in-house solver naoe-FOAM-SJTU is presented. The naoe-FOAM-SJTU solver is developed using an open source tool kit, OpenFOAM. Reynolds-averaged Navier?Stokes (RANS) equations are chosen as governing equations and the volume of fluid (VOF) is employed to capture the two phases interface. Incoming wave group is generated by imposing the boundary conditions of the tank inlet. A spectrum based correction procedure is developed to make the measured spectrum approaching to the target spectrum. This procedure can automatically adjust the wave generation signal based on the measured wave elevation by wave height probe in numerical wave tank. After 3 to 4 iterations, the measured spectrum agrees well with the target one. In order to validate this method, several wave spectra are chosen and validated in the numerical wave tank, with comparison between the final measured and target spectra. In order to investigate a practical situation, a modified Wigley hull is placed in the wave tank with incoming irregular waves. The wave-induced heave and pitch motions are treated by Fourier analysis to obtain motion responses, showing good agreements with the measurements.