基于桁架模型和多孔介质模型的柔性网衣和周围流场单向耦合方法研究

柔性网衣作为渔业养殖的重要组成部分,其水动力的计算分析对网箱的安全使用至关重要。基于桁架模型和多孔介质模型对柔性网衣和流场单向耦合计算的数值模拟方法进行了研究。通过code_aster软件基于桁架模型计算出网衣变形后节点坐标和网衣受力,然后将节点坐标导入OpenFOAM识别出网衣变形后的多孔介质区域,并基于多孔介质模型模拟水流作用下网衣周围的流场和计算网衣的受力。为了验证模型的正确性,将不同工况对圆形网衣的数值模拟结果与物理试验数据进行比较,包括code_aster中计算的网衣变形形状、OpenFOAM模拟的网衣周围流场以及在code_aster和OpenFOAM中计算的网衣受力。计算结果表明,圆形网衣在单纯来流下变形明显,与试验结果保持一致,证明了网衣变形和水动力单向耦合方法的准确性。     

Fluid—Structure Interaction Analysis of Flexible Plate with Partitioned Coupling Method

The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.     

A three-dimensional wake impingement model and applications on tandem oscillating foils

Potential flow based vortex numerical methods have been widely used in aerodynamics and hydrodynamics. In these methods, vortices shed from lifting bodies are traced by using vortex filaments or dipole panels. When the wake elements encounter a downstream body, such as a rudder behind a propeller or a stator behind a rotor, a treatment is necessary to divert the wake elements to pass by the body. This treatment is vital to make wake simulations realistic and to satisfy the non-penetration condition during wake body interaction. It also helps to avoid pure numerical disturbances such as when a vortex filament or an edge of a dipole panel passes through the collection point of a body element; this is a singularity for induced velocity and it will introduce a large numerical disturbance. This necessary treatment for three-dimensional problems with geometrical complexity has not been found to date. In this study, a wake impingement model was developed to divert wake elements to slip over the body surface, model the vortex/body interaction, and predict forces on fluctuating components. The model was also tested on configurations of oscillating foils in tandem with an existing panel method code. Simulation results with the wake impingement model are shown to be in closer agreement with limited published experimental data than those without the model. With the established wake impingement model, force fluctuations on the after body due to the wake vortex impingement were investigated based on a series of simulations. The series varied several parameters including distance between two foils, oscillating frequency, span, rear foil pitch angle, swap angle and vertical position.     

细长输流管内外流耦合振动特性研究

管内流动会影响输流管的振动响应,目前关于输流弹性管涡激振动方面的研究较少。基于计算流体力学(CFD)方法,开展内外流对细长输流弹性管振动特性影响的研究。首先在不考虑内流的情况下将弹性管涡激振动数值预报结果与模型试验数据进行对比,验证了数值方法的可靠性。再者考虑内外流耦合作用情况下,对不同内流流速下细长输流弹性管振动位移时—空分布、顺流向最大平均偏移、振动轨迹、内部横向涡的形成与分布等进行了对比分析。结果发现,与外流流速相比,内流流速的增加虽然难以改变弹性管的主振模态,但对沿管体的振动强度影响显著。顺流向最大偏移处管体运动轨迹发生明显的变形和跳跃。在剪切外流和均匀内流对弹性管的联合作用下,沿管跨方向模态间能量转换频繁,伴随着间歇性出现或消失的沿弹性管传播的行波组分,这主要归因于复杂的双重流固耦合系统(外流—管体,内流—管体)。在内流以附加质量力、离心力和科氏力形式的激励下,弹性管内二次流现象明显。在振动过程中,内部横向涡沿管壁生成、脱落并逐渐散布于整个横截面。     

海底管线柔性导流板的变形方程研究

海底输油管线是海洋油气田开发的生命线工程,在复杂的海洋环境中易发生破坏.为使管线实现自埋防护,提出用柔性材料代替刚性导流板安装在管线顶部,不仅可以增加阻水面积,加大冲刷深度和范围,而且可以减缓刚性导流板对管线上部水流的扰动强度和尾流涡旋的扰动强度,减小管线的振动.利用力学关系由变形方程推导出海底管线与海底间距离为零以及不为零情况下柔性导流板变形的二次曲面方程和自由端最大挠度.在单向流条件下将不同材料不同长度的柔性导流板安装在管线中轴上方开展试验研究,测量柔性导流板在水流作用下的变形曲面.对比分析利用公式计算柔性导流板的自由端最大挠度和曲面变形数据与试验结果,发现两者吻合较好.     

Numerical Prediction of Vortex Induced Vibrations on Top Tensioned Riser in Consideration of Internal Flow

In consideration of the effect of the internal flowing fluid and the external marine environmental condition on the vortex-induced vibration （VIV） of top tensioned riser （Till）, the differential equation is derived based on work-energy principles and the riser near wake dynamics is modeled by Facchinetti＇ s wake oscillator model. Then Galerkin＇ s finite element approximation is implemented to derive the nonlinear matrix equation of the coupled equations and file corresponding numerical programs are compiled which solve the coupled equations directly in the time domain. The comparison of the predicted results with the recent experimental results and the prediction of SHEAR7 is performed. The results show the validity of the proposed method on the prediction of VIV of deep water risers. The effect of internal flow on the dynamic characteristics and dynmnic response of the riser is analyzed and several valuable conelusions are drawn.     

柔性管卷管式铺管上卷有限元模拟

采用卷管法进行海底管道铺设过程中,管道首先通过牵引作用上卷于卷筒进行储存。管道与卷筒发生非线性接触,可能会产生复杂的塑性变形和局部屈曲。通过全尺寸柔性管力学性能试验获得柔性管轴力—应变以及弯曲—曲率等非线性力学性能关系,将试验所得的非线性材料性能参数导入建立的两种柔性管上卷ABAQUS有限元模型（梁—实体单元模型与壳和桁架—实体单元模型）,实现柔性管较大轴向抗拉刚度和较小抗弯刚度的同步模拟以及管道与卷筒的非线性接触响应特征。通过对比分析两种有限元模型数值模拟得到的管道弯矩、弯曲曲率、管道轴力、管道与卷筒的接触压强等数据,发现在管道上卷过程中管道沿副法线方向的SM3弯矩占据其弯曲变形主导地位;管道与卷筒之间的摩擦效应对于管道轴力的影响较为显著;管道与卷筒的最大接触压强主要发生在卷管过渡段区域。     

杭州湾海底管道冲刷自埋演化过程初步研究

针对杭州湾海底管道的实际工程状况,使用声学探测技术对海底管道冲刷自埋演化过程进行了现场检测,验证了海底管道从初始敷设状态依次发展为孔道冲刷、部分掩埋、尾流冲刷、尾流冲刷平衡、再次孔道冲刷和回淤自埋状态的演化过程.采用数值模拟方法分析了海底管道自埋演化过程中的典型状态,结果表明:管道下方的渗流和天然微孔隙是孔道冲刷起动的...     

适于模拟不规则水域波浪的缓坡方程两种数值模型比较

本文分析比较了适于不规则水域波浪模拟的椭圆型缓坡方程两种数值模型。两种数值模型均采用有限体积法离散,分别基于四叉树网格和非结构化三角形网格建立。首先结合近岸缓坡地形上波浪传播的经典物理模型实验对两种数值模型分别进行了验证,并结合计算结果对比分析了两种模型的计算精度和效率。计算结果表明,两种数值模型均可有效地模拟近岸波浪的传播变形;相对非结构化三角形网格下的模型,基于四叉树网格建立的数值模型在数值离散和求解过程中无需引入形函数、不产生复杂的交叉项,离散简单,易于程序实现,且节约计算存储空间,计算效率高。     

Mechanism for enhanced diffusivity in the deep-sea perpetual salt fountain

The mechanism of enhanced diffusivity occurring in the deep-sea perpetual salt fountain has been investigated experimentally and numerically. Some factors which possibly contribute to the enhanced diffusivity were found to be the pipe oscillation with ocean waves and its baffled wall surface. Field experiments in the ocean (Onagawa Bay of Miyagi, Japan) and numerical simulations were performed to study and confirm the dynamics of the flow and heat transport with enhanced diffusivity occurring in upwelling deep-sea water. The agreement between the field experimental data and the numerical solutions of an oscillating-wall boundary condition imposed on the baffled pipe is encouraging, and it indicates the baffled pipe surface subject to the oscillatory motion leads to the enhanced diffusivity. The buoyancy force and then upwelling velocity can be greatly increased by the enhanced diffusivity. The dominant mechanism is the occurrence of complicated vortices and vortex shedding leading to efficient mixing and enhanced diffusion.     

指数剪切流作用下海洋立管振动特性研究

细长海洋立管在复杂来流作用下的涡激振动（VIV）是海洋工程领域备受关注的热点问题。采用双向流—固耦合方法,对立管在指数剪切来流作用下的涡激振动实施了数值仿真研究。基于均方根振幅包络图的显性模态分析发现,与线性剪切流工况相比,当流速较小时,立管模型在指数剪切流作用下振幅最大值相对较小;流速较大时,指数剪切流的非线性分布致使立管的振动响应增强。通过对立管模型均方根振幅包络峰、谷处的振动响应频率分析发现,在波峰位置处振动频率单一且较为稳定,而在相邻波谷位置处多频共存现象显著,两者有显著差异。沿管的轴向波形主要表现为驻波主导和驻波—行波混合模式,行波一般间歇发生,流速越大发生频次愈高,其在横流向的传播方向通常由高流速区段传递至低流速区段。指数剪切流作用下,沿管体轴向各截面位置处的旋涡脱落模式差异显著,高流速段尾流区的旋涡发放具有较强的周期性,涡管完整且与管轴存在一定偏离角度,低流速段蜂窝状离散涡旋较多。     

波浪与抛石潜堤相互作用的MPS法数值模型研究

为了研究波浪与抛石潜堤相互作用过程中大自由表面变形和堤内渗流等强非线性紊流运动问题,利用改进的MPS法,建立了模拟波浪与抛石潜堤相互作用的MPS法数值计算模型。模型将抛石潜堤假定为均质多孔介质,采用Drew的二相流运动方程描述多孔介质内外的流体运动;通过在动量方程中增加非线性阻力项,并引入亚粒子尺度紊流模型,模拟波浪与可渗结构物相互作用过程中的紊流运动。选取“U”型管中多孔介质内渗流过程和孤立波与可渗潜堤相互作用两个典型的渗流问题,通过将数值计算结果与理论解和实测值的对比分析,对所提出的MPS法紊流渗流模型的模拟精度进行验证。结果表明:基于改进的MPS法构建的垂向二维紊流渗流模型可以很好地再现“U”型管中多孔介质内渗流以及波浪作用下可渗潜堤内外的复杂流场,显著缓解流-固界面处的压力震荡与粒子分布不均匀问题,实现了较高的模拟精度。     

Submarine debris flow impact on pipelines — Part II: Numerical analysis

Computational Fluid Dynamics (CFD) numerical analysis was employed to analyze the situations tested experimentally, as described in Part I. The methodology and results of the CFD analyses are discussed and compared with the observations made from the experiments. The numerical model performed satisfactorily with regard to obtaining the impact forces exerted on the model pipe as well as simulating the hydroplaning phenomenon and estimating slurry flow heights. The experimental results were combined with the results of the CFD analyses to develop a practical method to compute the drag force caused by a submarine debris flow impact on a pipeline. The CFD analyses provided some insight to the separated region characterization, but the attempt to analyze the vortex shedding phenomenon as observed in the experiments was unsuccessful. Additional studies are required for better understanding of both the separated region characteristics and vortex shedding.     

基于物理造波和主动消波的黏性数值波浪水槽开发

基于推板造波理论和摇板造波理论,在Open FOAM平台上采用重叠网格技术建立黏性数值波浪水槽,并使用一种结合SIMPLE算法和PISO算法的PIMPLE算法对数值模型进行求解。利用开发的数值模型通过数值收敛性测试和网格独立性测试分别重点研究了时间步长、库朗数和网格尺寸对数值精度和计算效率的影响。并对比研究了此数值模型分别嵌入层流模型和湍流模型的计算精度和计算效率。实现的规则波和二阶有限振幅波与理论结果和试验结果吻合,验证了此黏性数值波浪水槽的造波和主动消波功能。基于二维数值波浪水槽,进一步研究了三维数值造波,数值计算结果与理论结果吻合良好。研究结果不仅验证了重叠网格在二维和三维两相流体域中求解运动物体与流场交互的可靠性和正确性,而且为使用此黏性数值波浪水槽解决更复杂的海洋工程问题提供了依据。     

Laboratory, numerical, and theoretical modeling of the flow in a far wake in a stratified fluid

The far-wake flow past a sphere towed in a fluid with high Reynolds and Froude numbers and with a pycnocline-form salt-density stratification is studied in a laboratory experiment based on particle image velocimetry and in numerical and theoretical modeling. In the configuration under consideration, the axis of sphere towing is located under a pycnocline. Flow parameters, the profiles of density and average velocity, and the initial field of velocity fluctuation in numerical modeling are specified from the data of the laboratory experiment. The fields of fluid velocity at different times and the time dependences of integral parameters of wake flow, such as the average velocity at the axis and the transverse width of the flow, are obtained. The results of numerical modeling are in good qualitative and quantitative agreement with the data of the laboratory experiment. The results of the laboratory experiment and numerical modeling are compared to the predictions of a quasi-linear and quasi-two-dimensional theoretical model. The time evolution of both the average velocity at the axis and the transverse width of the wake is obtained with the model and is in good agreement with the experimental data. The results of numerical modeling also show that, under the effect of velocity fluctuation in the wake, internal waves whose spatial period is equal to the characteristic period of the wake’s vortex structure are excited efficiently in the pycnocline.     

Ships with ventilated cavitation in seaways and active flow control

Bottom ventilated cavitation has been proven as a very effective drag reduction technology for river ships and planning boats. The ability of this technology to withstand the sea wave impact usual for seagoing ships depends on the ship bottom shape and could be enhanced by some active flow control devices. Therefore, there is the need in numerical tools to estimate the effects of bottom changes and to design such devices. The fundamentals of active flow control for the ship bottom ventilated cavitation are considered here on the basis of a special model of cavitating flows. This model takes into account the air compressibility in the cavity, as well as the multi-frequency nature of the incoming flow in wavy seas and of the cavity response on perturbations by incoming flow. The numerical method corresponding to this model was developed and widely manifested with an example of a ship model tested in a towing tank at Froude numbers between 0.4 and 0.7.The impact of waves in head seas and following seas on cavities has been studied in the range of wavelengths from 0.45 to 1.2 of the model (or ship) length. An oscillating cavitator-spoiler was considered as the flow controlling devices in this study. The oscillation magnitude and the phase shift between cavitator oscillation and the incoming waves have been varied to determine the best flow control parameters. The main results of the provided computational analysis include oscillations of cavity surface, of the pressure in cavity and of the moment of hydrodynamic load on the cavitator. The major part of computations has been carried out for the flap oscillating at the frequency coinciding with the wave frequency, but the effect of a frequency shift is also analyzed.     

Prediction of hydrodynamic forces on submarine pipelines using an improved Wake II Model

The hydrodynamic force model for prediction of forces on submarine pipelines as described includes flow history effect (wake effects) and time dependence in the force coefficients. The wake velocity correction is derived by using a closed-form solution to the linearized Navier–Stokes equations for oscillatory flow. This is achieved by assuming that the eddy viscosity in the wake is only time dependent and of a harmonic sinusoidal form. The forces predicted by the new Wake (Wake II) Model have been compared to Exxon Production Research Company Wake Model in terms of time histories (force shape) and magnitudes of peak forces. Overall, the model predictions by the Wake II Model are satisfactory and represent a substantial improvement over the predictions of the conventional models. The conventional force models representing adaptations of Morison's equation with ambient velocity and constant coefficients give predictions that are in poor agreement with the measurements especially for the lift force component. The Wake II Force Model can be used for submarine pipeline on-bottom stability design calculations for regular waves with various pipe diameters.     

Prediction of hydrodynamic forces on submarine pipelines using an improved Wake II Model

The hydrodynamic force model for prediction of forces on submarine pipelines as described includes flow history effect (wake effects) and time dependence in the force coefficients. The wake velocity correction is derived by using a closed-form solution to the linearized Navier–Stokes equations for oscillatory flow. This is achieved by assuming that the eddy viscosity in the wake is only time dependent and of a harmonic sinusoidal form. The forces predicted by the new Wake (Wake II) Model have been compared to Exxon Production Research Company Wake Model in terms of time histories (force shape) and magnitudes of peak forces. Overall, the model predictions by the Wake II Model are satisfactory and represent a substantial improvement over the predictions of the conventional models. The conventional force models representing adaptations of Morison's equation with ambient velocity and constant coefficients give predictions that are in poor agreement with the measurements especially for the lift force component. The Wake II Force Model can be used for submarine pipeline on-bottom stability design calculations for regular waves with various pipe diameters.     

基于平行双板结构的圆管尾流抑制方法研究

通过数值模拟方法对基于平行双板的圆管尾流进行抑制研究。以二维不可压缩粘性流体方程为流动控制方程,并采用稳定化流体有限元方法求解。基于数值结果,详细探讨了平行双板的安装角度和板长等关键特征参数对尾流的抑制作用和影响。数值分析发现,平行双板结构对尾流的抑制作用明显优于传统单板控制。计算结果显示,基于平行双分隔板的尾流抑制最优位置为从圆管后驻点起向上游40°~50°处,并且与单板相比,在更短板长条件下,亦可完全抑制圆管的旋涡脱落,从而尾流从非定常涡旋流转入定常流动。当安装位置为θ≥80°时,被控圆管表面压力分布与无控圆管类似,说明此时平行双板控制无减阻效应。最后,基于数值模拟结果,对本方法的尾流高效控制效应,给出了相应的流动机制。     

低雷诺数下圆柱强迫振动特性光滑粒子流体动力学模拟

基于光滑粒子流体动力学(SPH)方法,开发了能够准确描述水流作用下圆柱强迫振动特性的数学模型。通过引入适合于无网格粒子法的开边界算法,来模拟出入流边界条件,建立了具有造流功能的SPH数值水槽。圆柱及计算域的上下边界均采用修正的动力边界条件进行模拟。借助于粒子位移矫正和压力修正算法,避免了圆柱周围流体粒子压力大幅震荡以及结构下游区域出现空腔等非物理性现象。使用典型的圆柱绕流数据,验证了所建SPH模型的计算性能,研究了固定圆柱在低雷诺数情况下的尾涡脱落模式和升阻力变化规律。明确了低雷诺数下强迫振动圆柱在频率锁定以及非锁定区间内的升力变化规律,量化了升力与外界激励频率之间的关系。