FPSO甲板上浪研究现状

浮式生产储油轮（FPSO）是当今海洋石油开发的主流设施。FPSO通常采用单点系泊系统,船首暴露在波浪作用下,定位于特定海域进行长达数年的长期作业,上浪的风险很大。甲板上浪及其产生的载荷已是当前国际FPSO工程和研究领域的热点之一。介绍了国内外在试验研究和理论研究两方面对FPSO甲板上浪研究的进展情况,并建议我国开展这方面的研究。     

垂荡式波浪能装置运动模型数值分析

垂荡式波浪能装置在海洋可再生能源开发中被广泛应用,通过浮子与摇杆在垂荡方向的相对运动吸收波浪能。在以往相关的运动预报数值分析中,通常基于微幅波假设,仅考虑浮子与摇杆在垂荡方向的运动,忽略摇杆其他自由度运动。建立并求解了垂荡式波浪能装置的非线性联合运动方程组,分析垂荡式波浪能装置的波浪载荷、浮子与平台连接处的受力情况。数值计算系统的运动响应,并将计算结果与已有的试验数据进行比较验证,结果表明数值模拟的垂荡式波能装置的运动响应与试验结果相符合。最后,应用本计算方法分析PTO(power take-off)参数对波能装置发电性能的影响。     

缓坡地形拍岸浪数值模拟研究

采用基于完全非线性Boussinesq方程建立的精细化海浪数值模型——Funwave-TVD模式及水槽实验数据进行拍岸浪数值模拟研究,并对不同入射情况拍岸浪进行数值模拟得出:拍岸浪缓坡传播过程中,均符合波高先增大直至破碎后迅速减小这一现象;在典型坡度1.91°(1∶30)情况下,进行了大量数值模拟,当仅入射波高增大时,破碎波高增大,破碎位置向深水区移动;当仅入射周期增大时,除上述特征外,当周期增大到一定程度时,这种趋势显著减弱;同时,当入射周期和波高较大时,非线性效应增强,模拟结果出现不稳定。下一步将利用实验与数值模拟相结合的优势,对前人拍岸浪预报方法和理论进行深入分析研究,为实现拍岸浪精细化预报提供优化解决方案。     

斜坡上波浪破碎与越浪非静压数值模拟

为合理确定防浪建筑物的越浪量,基于含非静水压力梯度项的非线性浅水方程建立了近岸波浪越浪数值模型。通过采用域内造波、消波并结合波前静压假定的破碎模型,模拟了规则波和不规则波在斜坡上的波浪传播变形,并在此基础上进行了越浪量数值计算。数值计算结果与物理模型实验结果表明,非静压模型可合理地描述波浪破碎点位置、破碎后的波高、增减水以及斜坡上的堤后越浪量。数值模型具有较高的计算精度和计算效率,可为实际工程防浪建筑物越浪以及堤顶高程的设计提供一种新的数值研究手段。     

高桩挡板透空式防波堤消浪性能数值研究

以RANS方程为控制方程,基于有限体积法,在动量方程中添加源项,建立了具有造波‐消波功能的数值波浪水槽。利用建立的源项造波数值波浪水槽,模拟了高桩挡板透空式防波堤在规则波作用下的消浪效果,完整地再现了堤前堤后的流态,分析了挡板相对入水深度对透浪系数的影响。在与试验值及拉帕公式对比后发现,数模计算结果与试验值较接近,拉帕公式偏大。研究了堤顶相对宽度、相对水深、相对波高对透浪系数的影响并提出修正公式,修正公式与计算值和试验值吻合较好。     

碟形越浪式波能发电装置的三维数值模拟

基于计算流体力学软件的三维数值模拟技术,分析了碟形越浪式波能发电装置的越浪性能,通过构建基于水气两相VOF(Volume of Fluid)模型的三维数值波浪水槽对该装置进行三维数值模拟研究,数值计算结果与物理试验结果相互对比验证较为吻合,验证了所构建的三维数值波浪水槽的可靠性,通过考察装置的坡度、导流叶片个数、干舷高度对越浪性能的影响确定装置的最优结构参数。结果表明,在装置的斜坡面边缘增加回流板可减少波浪的反射,提高装置的越浪性能。在数值模拟中将装置的斜坡面边缘处安装回流板对碟形越浪式波能发电装置参数进行优化,通过分析回流板的长度对装置越浪性能的影响来探索最优回流板长度。     

多消浪室开孔沉箱消浪特性的数值模拟研究

本文基于雷诺平均的Navier-Stokes方程和k-ε模型求解湍流流动,采用流体体积法(Volume of Fluid,VOF)追踪自由表面运动,建立无反射波浪数值水槽,对多消浪室开孔沉箱的消浪特性进行数值模拟研究.将单消浪室和多消浪室开孔沉箱反射系数和结构前波面分布的数值分析结果与物理模型试验结果进行对比验证,两者...     

新型直立式透空堤消浪性能数值研究

为使防波堤同时具有良好的掩护效果和水体交换能力,提出了两种带有透浪通道的新型直立式防波堤。基于Fluent求解器建立了三维数值波浪水槽,通过与试验结果对比,验证了该数值水槽求解波浪与透空堤作用具有较高的精度。对两种防波堤在规则波作用下的透浪特性进行了研究,结果表明:透射系数K_t与透空率呈正线性相关,且可通过调整透浪通道间距,使相同透空率下K_t降低20%～30%。对同一结构,K_t随相对波长的增大而显著增大,但受相对波高的影响较小。在透空率大于0.16后,异型沉箱防波堤的消浪性能明显优于错位沉箱。基于数值计算结果,给出了以上两种透空堤波浪透射系数的经验公式。     

运动坐标系下的海浪数值模拟

以运动坐标系中的海浪谱为靶谱应用静止坐标系中随机波面的模拟方法─—线性叠加法和线性过滤法模拟匀速运动坐标系中的随机波面,给出了两种方法的原理及模拟手续。通过模拟实验发现,在高速运动的坐标系中线性叠加法优于线性过滤法,在低速情形,线性过滤法好于线性叠加法,并在模拟技术上给出了各参数的选取方法。     

基于格子Boltzmann方法的斜坡堤越浪数值模拟研究

准确确定越浪量对于斜坡堤设计有重要意义。利用格子Boltzmann方法(LBM),并采用主动吸收式速度入口造波、出流边界消波、VOF方法追踪自由表面以及静态Smagorinsky模型模拟紊流运动,建立二维数值波浪水槽,对光滑斜坡堤上规则波与不规则波越浪进行数值模拟。模拟结果与试验值及其他数值模型结果比较表明,二维LBM数值波浪水槽具有模拟斜坡堤越浪的能力,但对于破碎较为剧烈的越浪过程模拟,该模型还存在一定的不足,未来可从提高自由表面模型精度等方面进一步改善其性能。     

Numerical prediction of green water incidents

Green water loads on moored or sailing ships occur when an incoming wave significantly exceeds the freeboard and water runs onto the deck. In this paper, a Navier–Stokes solver with a free surface capturing scheme (i.e. the VOF model; [Hirt and Nichols, 1981]) is used to numerically model green water loads on a moored FPSO exposed to head sea waves. Two cases are investigated: first, green water on a fixed vessel has been analysed, where resulting waterheight on deck, and impact pressure on a deck mounted structure have been computed. These results have been compared to experimental data obtained by [Greco, 2001] and show very favourable agreement. Second, a full green water incident, including vessel motions has been modelled. In these computations, the vertical motion has been modelled by the use of transfer functions for heave and pitch, but the rotational contribution from the pitch motion has been neglected. The computed water height on deck has been compared to the experimental data obtained by [Buchner, 1995a] and it also shows very good agreement. The modelling in the second case was performed in both 2-D and 3-D with very similar results, which indicates that 3-D effects are not dominant.     

Numerical study of the behavior of heated water discharged into the ocean

The effect of the earth's rotation on the behavior of heated water discharged from gigantic nuclear power plants was investigated numerically. It is shown that the effect is significant even for the existing scale of discharge. Horizontal spreading of the heated water is suppressed in offshoreward and northward directions (when heated water is discharged in an eastward direction from a westfacing coast in the northern hemisphere), but accelerated in the southward direction, and the shape of the formed warm water region becomes asymmetric. In this paper, several calculations were done using exaggerated Coriolis parameters (f) in order to demonstrate the effect of the earth's rotation. An increase off enhances above effect, but the situation of the inner region of the produced warm water mass is much different from that in the outer region near its margin. The southward expansion of surface isothermal contours in the inner region increases with an increase inf, but those in the outer region do not. This results in an increase of the sharpness of the front which is generated near the southern edge of the warm water region. We have not found simple parameters which describe the temperature and velocity distributions in the warm water mass produced.     

Numerical study of the behavior of heated water discharged into the ocean

Matsuno and Nagata (1987) showed numerically that the spreading characteristics of the discharged heated water from power plants of existing scale is significantly influenced by the earth's rotation effect. Although the effect of the enhancement of the Coriolis parameter on the spreading characteristics of the formed warm water mass was discussed in order to demonstrate the rotation effect, other parameters such as the density difference between the discharged heated water and ambient water and the vertical eddy viscosity and diffusivity were fixed. In this paper, the dependence of the spreading characteristics on these parameters is examined. Then, it is shown that the overall shape of the formed warm water mass and density and velocity structure strongly depend on these parameters. Also, it is indicated that the behavior of the warm water mass under the rotation effects is too complicated to be described with a few parameters. For example, the internal radius of deformation seems to be one of the determinative parameters, and the increase of the density difference between discharged water and ambient water has a similar effect on the shape of the formed warm water mass as a decrease of the Coriolis parameter. However, a change of the two values has different effects on the detailed density structure and current structure of the warm water mass. The Prandtl number seems to determine some aspects of the veloczty field such as velocity magnitude and width of the southward flowing current zone. However, other features such as the thickness of the warm water mass are not determined by the Prandtl number.     

楔形体在波浪中自由入水的数值模拟

物体入水时波浪的影响不可忽略,基于流体力学模型采用VOF法,并利用自定义函数,模拟了楔形体的自由入水过程;同时结合推波板原理及海绵层消波理论实现了数值水槽的造消波,完成了波浪中楔形体自由入水的模拟,计算了楔形体入水时所受的水作用力、自由液面变化及物面压强分布等,研究了不同波高、周期以及在波浪不同位置入水时对楔形体的影响。结果表明:本文建立的数值模型可很好地模拟楔形体入水造成的射流及空泡的形成发展过程,波浪对楔形体入水的影响主要由波浪内部流场变化及表面波形决定,在波浪不同位置处入水对楔形体受力及入水形态均有较大影响。     

大洋采矿输送软管动力特性的数值研究

在大量勘探数据的基础上，对中国东太平洋7．5万km^2矿区采矿环境进行研究，得出集矿机是否安全行驶，直接关系着大洋采矿的成败。影响集矿机行驶的因素很多，其中输送软管是重要因素之一。通过理论分析计算和小型模拟实验相结合的研究方法，进行了软管与集矿机受力分析与仿真计算研究，开发了相应的应用软件。总结了大洋采矿系统中输送软管动态特性数值研究的最新研究成果，对我国大洋多金属结核矿产资源的开发，深海采矿系统的设计和回采策略的制定有一定的参考价值和实用意义。     

半封闭海湾的水交换数值模拟研究

基于采用不规则三角网格和有限体积方法的FVCOM模式,建立半封闭海湾-湛江湾附近海域的三维潮汐潮流数值模型,通过验证,结果与观测数据符合良好,并在此模型基础上对湛江湾的水交换状况进行了数值模拟,将湛江湾划分成3个区域,针对各区域进行了水交换能力研究,研究结果表明:由于湛江湾内不同区域的地形和地理位置变化较大,使得湛江湾内不同区域交换能力相差较大,其中,靠近湛江湾出口处交换能力最强,特呈岛以北海域交换能力最弱,交换时间与物质的初始浓度无关,与投放时刻和外源强迫密切相关,在治理湛江湾环境时,应分区进行,注意选择污染物排放时间和位置。     

Numerical study on water waves and wave-induced longshore currents in Obaky coastal water

In this paper, the water waves and wave-induced longshore currents in Obaky coastal water which is located at the Mediterranean coast of Turkey were numerically studied. The numerical model is based on the parabolic mild-slope equation for coastal water waves and the nonlinear shallow water equation for the wave-induced currents. The wave transformation under the effects of shoaling, refraction, diffraction and breaking is considered, and the wave provides radiation stresses for driving currents in the model. The numerical results for the water wave-induced longshore currents were validated by the measured data to demonstrate the efficiency of the numerical model. Then the water waves and longshore currents induced by the waves from main directions were numerically simulated and analyzed based on the numerical results. The numerical results show that the movement of the longshore currents was different while the wave propagated to a coastal zone from different directions.     

Numerical and experimental study on seakeeping performance of ship in finite water depth

Vessels operating in shallow waters require careful observation of the finite-depth effect. In present study, a Rankine source method that includes the shallow water effect and double body steady flow effect is developed in frequency domain. In order to verify present numerical methods, two experiments were carried out respectively to measure the wave loads and free motions for ship advancing with forward speed in head regular waves. Numerical results are systematically compared with experiments and other solutions using the double body basis flow approach, the Neumann-Kelvin approach with simplified m-terms, and linearized free surface boundary conditions with double-body m-terms. Furthermore, the influence of water depths on added mass and damping coefficients, wave excitation forces, motions and unsteady wave patterns are deeply investigated. It is found that finite-depth effect is important and unsteady wave pattern in shallow water is dependent on both of the Brard number τ and depth Froude number F_h.     

人工岛对海湾水环境影响数值研究

针对大连湾人工岛的建设预案,在人工岛多种长度及宽度的组合条件下,运用MIKE21 flow model模块进行人工岛建设前、后的海湾流场及污染物扩散运动模拟,探讨研究嵌入人工岛后海湾水体交换、污染物扩散及污染物存留率的变化规律,进而评价人工岛对海湾水环境的影响.在不同的人工岛功能定位前提下,综合考虑人工岛和水环境的诸多因素,分别给出了可供规划参考的人工岛推荐方案.     

四叉树网格下的椭圆型缓坡方程数值模型研究

波浪是近岸海域关键的水动力因素之一。考虑到近岸地形复杂、波浪演化显著的特点,建立了四叉树网格体系下的椭圆型缓坡方程数值模型,采用有限体积法对模型进行数值离散,应用GPBiCG(m, n)算法求解离散后的控制方程。模型中根据波浪波长布局计算网格,生成多层次四叉树网格,对复杂计算域有较好的适应性,并且在离散和方程求解中无需引入形函数、不产生复杂的交叉项,节约了存储空间和计算时间。将模型成功应用于物理模型实验及Acapulco海湾的波浪场数值模拟,结果表明该模型能够准确、高效地模拟近岸波浪场,可为近岸波浪场的模拟提供一定的理论和技术支持。