Heat flow analysis of an FPSO topside model with wind effect taken into account: A wind-tunnel test and CFD simulation

The aim of this paper is to study the feasibility of a computational fluid dynamics (CFD) method to examine the effect of wind on the thermal-diffusion characteristics of floating production storage and offloading (FSPO) topside models subject to fire. It is motivated by the need to identify the fire loads on FPSO topsides, taking into account the effects of wind speed and direction as well as the effects of geometry of the FPSO topsides. The results of a wind-tunnel test and CFD simulation undertaken for these purposes on a 1/14-scale FPSO topside model of a VLCC class FPSO unit are reported here. In the wind-tunnel test, the locations of the heat source of the fire are varied, as are the speed and direction of the wind, and the temperature distribution is measured. CFD simulations, using the ANSYS CFX (2009) program, were performed on the test model, with the results compared with the experimental results. It is concluded that wind has a significant effect on the thermal-diffusion characteristics of the test model and that the CFD simulations are in good agreement with the experimental results. The insights developed in this study will be very useful for the fire engineering of FPSO topsides.     

多级可分离式AUV系统设计与仿真分析

多级可分离式自主水下航行器(MS-AUV)可搭载多种不同功能的载荷舱,并在目标海域完成载荷舱布放。其在海洋探测,海防等领域具有较高的应用价值。对MS-AUV进行了主体结构和分离机构的设计和研究,其中功能载荷舱与航行器本体采用真空连接方式,实现多级连接,降低了机构的复杂性,提高了载荷舱分离运动的稳定性和安全性。为了研究载荷舱与航行器分离的安全性,采用CFD和六自由度(6-DOF)刚体运动学的耦合方法来模拟其在水下的分离运动。在仿真过程中,采用弹簧光顺和局部重构混合网格生成方法的非结构化动态网格可以很好地适应大距离多体分离运动。仿真结果表明,在一定工况条件下,提高MS-AUV初始航行速度V0和分离机构弹簧刚度K有利于载荷与航行器本体分离。通过上述研究工作,验证了载荷分离方法的可行性,缩短研制周期,降低设计成本,对多级可分离式AUV样机的设计制造具有重要的参考价值。     

Frequency domain analysis of dynamic response of drag dominated offshore structures

The paper describes a method for stochastic representation of the hydrodynamic drag forces on offshore structures subjected to irregular waves. It is shown that, for the case of zero current, it is possible to construct a genuinely quadratic representation of the drag force which reproduces the statistical properties of the standard formulation of the drag force very closely, and which at the same time has sufficient flexibility to ensure a spectral density that accurately approximates the desired force spectrum. The distinct advantage of the new representation is that it brings dynamic analysis of extensive linear structures back into the frequency domain.     

基于CFD模拟的单点系泊船型网箱浮架结构分析与设计

针对单点系泊(SPM)网箱系泊力大于传统网箱,离岸环境承受风浪时风险更高的问题,研究设计了一种相较传统圆形网箱拥有更小系泊力的船型网箱,运用正交原理制定了三因素三水平的试验方案,对网箱浮架进行了计算机数值模拟,分析网箱在波流共同作用下的水动力特性变化。研究表明,在兼顾经济成本与网箱安全性的条件下,长24 m,宽9 m,头部角度为60°的网箱浮架具有较优良的水动力性能和相对较低的养殖成本。该网箱系统能随着潮流改变方向而在锚点周围移动,从而极大减小网箱底部的养殖残渣堆积并降低锚固安装成本,可为我国深水网箱的发展提供一定的理论参考依据。     

Numerical Study on the 3-D Complex Characteristics of Flow Around the Hull Structure of TLP

Vortex-induced motion is based on the complex characteristics of the flow around the tension leg platform (TLP) hull. By considering the flow field of the South China Sea and the configuration of the platform, three typical flow velocities and three flow directions are chosen to study the numerical simulation of the flow field characteristics around the TLP hull. Reynolds-averaged Navier–Stokes equations combined with the detached eddy simulation turbulence model are employed in the numerical study. The hydrodynamic coefficients of columns and pontoons, the total drag and lift coefficients of the TLP, the formation and development of the wake, and the vorticity iso-surfaces for different inlet velocities and current directions are discussed in this paper. The average value of the drag coefficient of the upstream columns is considerably larger than that of the downstream columns in the inlet direction of 0°. Although the time history of the lift coefficient demonstrates a “beating” behavior, the plot shows regularity in general. The Strouhal number decreases as the inlet velocity increases from the power spectral density plot at different flow velocities. The mean root values of the lift and drag coefficients of the front column decrease as the current direction increases. Under the symmetrical configuration of 45°, the streamwise force on C4 is the smallest, whereas the transverse force is the largest. The broken vortex conditions in current directions of 22.5° and 45° are more serious than that in the current direction of 0°. In addition, turbulence at the bottom of the TLP becomes stronger when the current direction changes from 0° to 45°. However, a high inlet velocity indicates a large region influenced by the broken vortex and shows the emergence of the wake behind the TLP under the same current angle.     

强迫摇荡运动三维球体水动力系数的CFD计算

采用计算流体力学(CFD)的雷诺平均法(RANS方法),并考虑了自由表面的影响(采用VOF方法模拟自由表面),来模拟三维球体的单自由度强迫摇荡运动,得到纵荡、升沉及横摇的附加质量与阻尼系数。该数值模拟结果与三维势流理论计算的结果进行比较,趋势基本一致,但CFD方法更能凸显粘性的作用。本文所应用方法能合理给出浮式结构物的水动力系数,更精确描述海洋浮式结构物周围的流场,可应用于船舶与海洋工程浮式结构物的水动力性能研究。     

潮汐汊口水流结构和动力特性研究

潮汐汊口作为河网的重要单元,关系到水沙以及可溶污染物、营养物质的输运和分配,对河口三角洲的水文过程、地貌沉积和生态环境,都有至关重要的影响。本文选取国内最大的珠江河网的分汊顶点-马口、三水分汊,基于一个大潮潮周期内走航ADCP测量,分析潮汐汊口各个支汊的流速过程及三维流速的时空分布特征;结合非稳态的调和分析算法和小波连续变换确定汊口处主要的潮汐组分,利用潮流分离技术,剥离出占主导作用的潮波信号并探讨潮波动力在潮汐汊口处的空间分布形态。研究发现,由于潮波在西江马口和北江三水汊道内传播的不对称性,两汊间的流速过程存在1.5 h左右的相位差。而连接西北江的岗根汊道,由于存在频繁的水体交换和弯道水流的特性,岗根断面的动力作用较为复杂,其最大垂向平均流速在涨落急时刻分别向左、右岸偏移,横向断面均出现了显著的顺时针二次流。     

Investigation of a longfin inshore squid's swimming characteristics and an underwater locomotion during acceleration

In the present study, locomotion of a real longfin inshore squid (Doryteuthis pealeii) was numerically investigated. Geometry of a real squid was obtained using computed tomography (CT) images. In addition to a two-dimensional axisymmetric squid model, a modified squid model with no cavities around her head and an ellipse shaped model were generated with a fineness ratio (the ratio of body length to maximum body diameter) of 7.56. These numerical models were exposed to an acceleration with two different velocity programs. Added mass and basset forces on bodies were calculated during acceleration of the squid models starting from rest. Pressure and viscous drag forces were also calculated due to pressure variation along the squid models and friction on the surface of the models. The effect of a nozzle diameter on jet velocities and propulsive efficiency at all bodies were evaluated when time dependent velocity profiles (from 0 to 10 m/s in 0.5 and 1 s time durations) were set for the inlet of computational domain. The modified squid model required least thrust force during acceleration phase of time dependent velocity profile compared to the other models while the 0.02 m nozzle diameter provided largest propulsive efficiency for all models.     

半潜式超大型浮式结构水动力系数研究

以半潜式超大型浮式结构在动力响应分析中的各水动力系数为研究对象,经理论推导得到D’Alembert动力学方程中的移动式海上基地(MOB)单模块运动的结构质量、结构附加质量、静恢复力系数的简易计算公式。以MOB的"三模块模型"为例,研究其在6级海况浪向角为0°~90°条件下,各模块的附加质量系数及静恢复力系数的历时规律,以实例MOB中的第1个模块为代表展示了计算结果,并统计其最大值与文献资料中的结果进行对比。结果表明:运用本理论公式计算的结果与文献中所得结果相似,可验证本理论公式的正确性、可行性与合理性,为求解半潜式超大型浮式结构模块动力响应位移及转角提供简便的方法。     

感潮河段潮流数值计算

应用圣维南方程组、普赖斯曼隐格式及双扫描解法等提出了感潮河段洪潮遭遇的非恒定流数值计算模型,并对感潮河网中所常见的边界条件进行了分析研究。计算与实测结果比较表明,所采用的方法有效可行。     

Mixed layer depth front and subduction of low potential vorticity water in an idealized ocean GCM

It is known that there is a front-like structure at the mixed layer depth (MLD) distribution in the subtropical gyre, which is called the MLD front, and is associated with the formation region of mode water. In the present article, the generation mechanism of the MLD front is studied using an idealized ocean general circulation model with no seasonal forcing. First, it is shown that the MLD front occurs along a curve where u _g ·∇T _s = 0 is satisfied (u _g is the upper ocean geostrophic velocity vector, T _s is the sea surface temperature and ∇ is the horizontal gradient operator). In other words, the front is the boundary between the subduction region (u _g ·∇T _s > 0) and the region where subduction does not occur (u _g ·∇T _s < 0). Second, we have investigated subduction of low potential vorticity water at the MLD front, which has been pointed out by past studies. Since u _g ·∇T _s = 0 at the MLD front, the water particles do not cross the outcrop at the MLD front. The water that is subducted at the MLD front has come from the deep mixed layer region where the sea surface temperature is higher than that at the MLD front. The temperature of the water in the deep mixed layer region decreases as it is advected eastward, attains its minimum at the MLD front where u _g ·∇T _s = 0, and then subducts under the warmer surface layer. Since the deep mixed layer water subducts beneath a thin stratified surface layer, maintaining its thickness, the mixed layer depth changes abruptly at the subduction location.     

受迫流动下的枝晶生长相场法模拟研究

基于耦合流场和热噪声的相场模型及合理高效的三维动态求解域加速算法,定量模拟了在受迫流动下枝晶的非对称生长及流速对迎流、背流两侧的温度分布和层流层分布的影响.计算结果表明,受迫流动使迎流、背流两侧温度的分布与层流层分布呈现不对称状态,导致迎流侧与背流侧的过冷度不同,而熔体施加于枝晶界面前沿迎流侧的力还不足以抑制过冷度的作用,结果造成枝晶迎流方向优先生长,从而产生倾向于散热方向的倾斜;同时,由于迎流侧的实际过冷度大于背流侧,有利于促进迎流一侧枝晶生长速度以及稳定侧向分枝生长,从而导致了侧向分枝的非对称生长.随     

Numerical simulation of viscous cavitating flow around a ship propeller

In the present study, cavitation and a ship propeller wake are reported by computed fluid dynamics based on viscous multiphase flow theory. Some recent validation results with a hybrid grid based on unsteady Navier-Stokes (N-S) and bubble dynamics equations are presented to predict velocity, pressure and vapor volume fraction in propeller wake in a uniform inflow. Numerical predictions of sheet cavitation, tip vortex cavitation and hub vortex cavitation are in agreement with the experimental data, same as numerical predictions of longitudinal and transversal evolution of the axial velocity. Blade and shaft rate frequency of propeller is well predicted by the computed results of pressure, and tip vortex is the most important to generate the pressure field within the near wake. The overall results indicate that the present approach is reliable for prediction of cavitation and propeller wake on the condition of uniform inflow.     

考虑诱导速度的螺旋桨敞水性能CFD模拟

螺旋桨工作时在其周围形成诱导速度场,诱导速度随到桨叶距离的增大而衰减。采用CFD方法模拟螺旋桨敞水性能时,只能截取有限尺度的流域进行计算,此时计算域边界上诱导速度并不为零,将进口速度设为进速是近似的。一般采用足够大的计算域,使螺旋桨前方及侧面边界尽量远离桨叶。为了在较小的计算域中实现螺旋桨敞水性能的准确预报,提出在设定进口速度时计入螺旋桨诱导速度的CFD模拟方法。应用升力面方法计算诱导速度,将进口速度设为进速与诱导速度之和。逐步减小计算域尺度,考察敞水性能及压力分布计算结果的变化情况及精度。算例比较表明:通过考虑诱导速度,可以大幅度减小进口与螺旋桨的距离而不降低计算精度。     

基于势流模型的双体船耐波性能研究

双体船由于其优良的快速性和横稳性等特性,在民用和军用中应用广泛。针对无限水深下迎浪、斜浪中自由运动的WigleyⅢ双体船,利用线性势流软件WAMIT,通过调整片体间距,对双体船的水动力系数及船体运动响应进行计算和比较,分析不同片体间距的双体船的耐波性。结果发现,片体间距较小的双体船附加质量系数较小,垂荡、纵摇阻尼系数显著增大,垂荡和纵摇运动响应幅值增大,且附加质量系数为负值时的系数绝对大小和频率范围随着片体间距的减小而变大。同时,片体间距越小的双体船间的流体干扰作用令自由面运动响应更加剧烈,迎浪下船体垂荡运动响应略有增大,纵摇运动响应略有减小;而斜浪下垂荡、纵摇、横摇运动响应均增大,船体的耐波性较差。当片体间距较大,在某一频率入射波下,双体船运动响应接近为零,且该频率随着片体间距的增大而减小。另外说明,相较于阻尼系数,附加质量系数对船体运动的影响更大。     

深水浊流形成过程的水动力学模拟:以南海北部琼东南盆地浊流形成过程为例

Turbidity currents represent a major agent for sediment transport in lakes, seas and oceans. In particu-lar, they formulate the most significant clastic accumulations in the deep sea, which become many of the world＇s most important hydrocarbon reservoirs. Several boreholes in the Qiongdongnan Basin, the north-western South China Sea, have recently revealed turbidity current deposits as significant hydrocarbon res-ervoirs. However, there are some arguments for the potential provenances. To solve this problem, it is es-sential to delineate their sedimentary processes as well as to evaluate their qualities as reservoir. Numerical simulations have been developed rapidly over the last several years, offering insights into turbidity current behaviors, as geologically significant turbidity currents are difficult to directly investigate due to their large scale and often destructive nature. Combined with the interpretation of the turbidity system based on high-resolution 3D seismic data, the paleotophography is acquired via a back-stripping seismic profile integrated with a borehole, i.e., Well A, in the western Qiongdongnan Basin; then a numerical model is built on the basis of this back-stripped profile. After defining the various turbidity current initial boundary conditions, includ-ing grain size, velocity and sediment concentration, the structures and behaviors of turbidity currents are investigated via numerical simulation software ANSYS FLUENT. Finally, the simulated turbidity deposits are compared with the interpreted sedimentary bodies based on 3D seismic data and the potential provenances of the revealed turbidites by Well A are discussed in details. The simulation results indicate that a sedimen-tary body develops far away from its source with an average grain size of 0.1 mm, i.e., sand-size sediment. Taking into account the location and orientation of the simulated seismic line, the consistence between normal forward simulation results and the revealed cores in Well A indicates that the turbidites should ha     

悬链线管体三维流场数值模拟研究

基于计算流体力学开源代码Open FOAM开展不同雷诺数(100≤Re≤3 900)下弯曲管体绕流流场三维数值模拟研究。对比分析不同雷诺数条件下弯曲管体的瞬时和时均流场特征,研究了弯曲管体不同位置处的时均压力系数和旋涡脱落频率。结果表明:当Re=100时,尾流中仅有展向涡存在;当300≤Re≤3 900时,流向涡出现且强度逐渐增加,展向涡减弱,随着局部曲率增大,沿弯管展向涡脱强度减小。随着雷诺数的增大,回流区沿流向的延伸尺度减小,但其沿展向的分布范围增大。沿弯管展向随着深度减小,前后驻点时均压力系数的绝对值递减,圆周时均压力系数变化幅度减小。在300≤Re≤3 900范围内,管体曲率会对展向涡脱频率产生一定影响。     

黄河口清水沟流路地貌特征研究

黄河入海流路自1976年5月改道至清水沟流路以来,其地貌形态一直处于剧烈的变化当中,以往的研究虽然关注了典型断面主槽的形态特征及过流能力的变化,但对于清水沟流路范围界限内河道断面整体形态及其沿程变化等地貌特征还缺乏系统研究.本文以2019年10月黄河口现行清水沟流路范围渔洼以下河段加测横断面为研究对象,辨识现状地貌分布...     

Euler-Euler two-phase flow simulation of tunnel erosion beneath marine pipelines

In this study an Euler-Euler two-phase model was developed to investigate the tunnel erosion beneath a submarine pipeline exposed to unidirectional flow. Both of the fluid and sediment phases were described via the Navier-Stokes equations, i.e. the model was implemented using time-averaged continuity and momentum equations for the fluid and sediment phases and a modified k−ε turbulence closure for the fluid phase. The fluid and sediment phases were coupled by considering the drag and lift interaction forces. The model was employed to simulate the tunnel erosion around the pipeline laid on an erodible bed. Comparison between the numerical result and experimental measurement confirms that the numerical model successfully predicts the bed profile and velocity field during the tunnel erosion. It is evident that the sediments are transported as the sheet-flow mode in the tunnel erosion stage. Also the transport rate under the pipe increases rapidly at the early stage and then reduces gradually at the end of the tunnel erosion beneath pipelines.