两层流体中振荡水平圆柱潜体水动力特性

研究了有限深两层流体中水平圆柱潜体的振荡辐射问题。在线性势流理论框架内,建立求解辐射势的多极展开理论,提出附加质量和阻尼系数的计算方法,进行了数值计算分析,并且与均匀流中的情况进行比较。结果表明,在圆柱潜体的某个振荡频率范围内,流体的分层效应对其附加质量和阻尼系数是有重要影响的。同时,对水平圆柱潜体振荡产生的自由面和内界面波动特性进行了数值分析。     

两层流体中多个振荡水平圆柱潜体水动力特性

采用解析方法研究了无限水深两层流体中多个振荡水平圆柱潜体的水动力特性。在线性势流理论框架内,基于多极子方法建立了辐射势的解析表达式,在此基础上导出了附加质量和阻尼系数的计算公式,采用边界元方法对这种解析方法进行验证,同时研究了两种不同工况下多个振荡水平圆柱潜体水动力的特性,结果表明,两层流体的密度比、圆柱的淹没深度以及圆柱的排列方式和间距等参数的变化对水平圆柱群附加质量和阻尼系数有很大的影响。     

两层流体中多个圆柱浮体波浪力的解析方法

采用解析方法研究了线性入射波作用下两层流体中多个圆柱形淹没浮体的渡浪力特性.首先基于多极子展开方法,建立了散射势函数的解析表达式,并进一步得到浮体散射渡浪力的计算公式,然后利用边界元方法验证了本文的解析解,最后分析了不同参数的变化对双圆柱形浮体结构波浪力的特有影响.     

两层流体中内波作用下Spar平台运动响应

研究两层流体中Spar平台在内波作用下的运动响应问题。在线性势流理论框架,提出在内波作用下Spar平台运动响应及分段式悬链线系泊张力特性的计算方法。数值分析两层流体内界面位置、入射内波的波长以及系泊索初始预张力对Spar平台运动响应及其系泊索张力特性的影响规律,结果表明内波对Spar平台纵摇运动响应的影响是小的,但对Spar平台纵荡与垂荡运动响应及其系泊索张力的影响是不可忽视的。因此,在Spar平台的设计中,考虑内波的影响是重要的。     

航行体近水面航行附加质量与阻尼的数值研究

航行体近水面航行时会引起周围流场中流体的惯性运动,动能会随波浪扩散而耗散,航行体近水面运动时固有周期和幅值衰减率受水的影响可以用附加质量和阻尼系数描述。运用计算流体力学软件STAR-CCM+开展航行体近水面自由横摇与纵摇衰减运动的数值仿真,根据不同初始角度条件下仿真得到的自由衰减时历曲线,分别计算横摇与纵摇的附加质量与阻尼系数,并结合实际运动情况,对不同初始角度条件下附加质量与阻尼系数的变化成因分别进行分析。     

大深度分层流体中二维淹没浮体的波浪力分析

研究了大深度分层流体中二维任意形状淹没浮体的波浪力特性。首先基于一种合适的格林函数,采用边界积分方程法研究了流体中浮体对水波散射问题,然后通过单个淹没圆柱体的透射能和反射能与解析方法结果的比较,对所提出的方法进行了验证,最后分析了在不同的几何和物理条件下几种形状的浮体对波浪力的特有影响,得到了一些有意义的结果,这对分层海洋中淹没浮体的设计具有重要的参考价值。     

连续层化流体中内波破碎的动力学机制的数值研究

对单频内波在白噪声的干扰下引起的不稳定现象进行了数值研究,分析了不同振幅和不同Pr数(Prandtl number)对静力不稳定、剪切不稳定、对流不稳定出现的时间、地点及持续时间的影响.结果显示,静力不稳定总是比其它2种不稳定先出现,其持续时间也是最长的;振幅越大(Pr＝0.72)时,内波越不稳定,导致3种不稳定出现的时间明显提前;Pr数增大(A＝0.256)时,3种不稳定出现的时间也略有提前.     

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

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

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

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

两层粘性流体中圆柱体受迫振荡数值模拟

研究两层粘性流体中无限长水平圆柱体的受迫振荡问题。在湍流模式下,采用VOF方法追踪两层流体的内界面,基于动网格技术模拟圆柱体的运动边界,对均匀流中横向振荡圆柱体的绕流场进行了数值模拟。计算受迫振荡圆柱体的升力系数、阻力系数随时间的演化曲线和圆柱体的尾涡分布,以及圆柱体的受迫振荡激发两层流体内界面的扰动,并与均匀流体的情况进行了比较分析。研究表明,流体的两层分层效应对受迫振荡圆柱体的升阻力系数和尾涡分布特性都有显著影响,在水下输油气管道涡激振动特性的工程评估中,应考虑流体的密度分层效应。     

Second-order random interfacial wave solutions for two-layer fluid with a free surface

1 Introduction Interfacial waves travelling along the interface between two fluids of different densities can be often observed in subsurface layers of the ocean since the upper subsurface layer is warmer over much of the o- cean (Umeyama, 2002). They are…     

Added mass and damping of a sphere section in heave

The Norwegian wave-power buoy¹ consists of a half-immersed floating sphere which is open to the sea at the bottom end. It is a two degree-of-freedom device involving the independent motion of the outer rigid sphere and the pressure across the internal free surface. A simpler model of the device is to represent the oscillatory flow through the bottom opening as another rigid body motion being that of the curved surface which would complete the sphere.The wave-induced forces on this surface and also the outer spherical surface due to independent oscillations of either are determined semi-analytically using a simple extension of the method of Havelock,² recently simplified and generalised by Hulme,³ for the complete half-immersed sphere. In particular the 2 × 2 added mass and damping matrices are determined as a function of frequency and relative size of bottom opening to sphere radius. These quantities are essential in any theoretical analysis of a multi-degree of freedom wave-energy device.     

Influence of the earth rotation on the interfacial wave solutions and wave-induced tangential stress in a two-layer fluid system

Interfacial waves and wave-induced tangential stress are studied for geostrophic small amplitude waves of two-layer .uid with a top free surface and a .at bottom. The solutions were deduced from the general form of linear .uid dynamic equations of two-layer .uid under the f -plane approximation, and wave-induced tangential stress were estimated based on the solutions obtained. As expected, the solutions derived from the present work include as special cases those obtained by Sun et al. (2004. Science in China, Ser. D, 47(12): 1147–1154) for geostrophic small amplitude surface wave solutions and wave-induced tangential stress if the density of the upper layer is much smaller than that of the lower layer. The results show that the interface and the surface will oscillate synchronously, and the in.uence of the earth’s rotation both on the surface wave solutions and the interfacial wave solutions should be considered.     

Hydroelastic response of a circular plate in waves on a two-layer fluid of finite depth

The hydroelastic response of a circular, very large floating structure （VLFS）, idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interfacial wave modes, of a given wave frequency, on a two-layer fluid of finite and constant depth. In linear potential-flow theory, with the aid of angular eigenfunction expansions, the diffraction potentials can be expressed by the Bessel functions. A system of simultaneous equations is derived by matching the velocity and the pressure between the open-water and the plate-covered regions, while incorporating the edge conditions of the plate. Then the complex nested series are simplified by utilizing the orthogonality of the vertical eigenfunctions in the open-water region. Numerical computations are presentedto investigate the effects of different physical quantities, such as the thickness of the plate, Young＇s modulus, the ratios ofthe densities and of the layer depths, on the dispersion relations of the flexural-gravity waves for the two-layer fluid.Rapid convergence of the method is observed, but is slower at higher wave frequency. At high frequency, it is found that there is some energy transferred from the interfacial mode to the surface mode.     

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.