围阱内流体振荡载荷试验研究

很多水面舰船为吊放声呐、潜水系统等特殊设备,设置了围阱结构。舰船航行过程中,围阱随边的壁面将受到底部剪切层的冲击,且围阱内的流体在外部来流和波浪的激励下将发生剧烈晃荡,流体的晃荡运动会对围阱内壁面产生较大的砰击载荷。针对该问题,开展了模型试验,对围阱内剪切层冲击及流体晃荡的载荷进行研究,分析剪切层冲击及流体振荡载荷随流速、围阱开口形式等参数的变化规律,通过对比围阱内不同位置测点的测试结果,确定流体晃荡载荷较大的区域,指导围阱结构设计。     

Effects of excitation angle and coupled heave–surge–sway motion on fluid sloshing in a three-dimensional tank

Sloshing waves in moving tanks is an important engineering problem, and most studies of this phenomenon have focused on tanks that are excited by forcing motion in a limited number of directions and with fixed excitation frequencies throughout the forcing. In practice, the excitation comprises multiple degree of freedom motion that potentially couples surge, sway, heave, pitch, roll, and yaw motions. In the present study, a time-independent finite difference method is used to simulate fluid sloshing in three-dimensional tanks filled to an arbitrary depth for various excitation frequencies and multiple degree of freedom motion. The numerical scheme developed here was verified by rigorous benchmark tests. The coupled motions of surge and sway are simulated for various excitation angles, frequencies and water depths. Five kinds of sloshing waves found under coupled surge–sway motions: diagonal, single-directional, square-like, swirling, and irregular waves. The effect of excitation angle on the frequency responses of different sloshing waves is analyzed and discussed in the present study. Further, the components of horizontal force of various sloshing waves are also presented. The coupled effect of surge, sway and heave motions is also discussed, and the results show that unstable sloshing occurs when the excitation frequency of the heave motion is twice the fundamental natural frequency. Moreover, the effects of heave motion on the different types of sloshing waves are explored. It is found that heave motion causes all of the sloshing waves to change type.     

基于MPS-FEM耦合方法研究波浪与结构的相互作用（英文）

Nowadays, an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea. As a result, some phenomena related to the violent fluid-elastic structure interactions(e.g., hydrodynamic slamming on marine vessels, tsunami impact on onshore structures, and sloshing in liquid containers) have aroused huge challenges to ocean engineering fields. In this paper, the moving particle semi-implicit(MPS) method and finite element method(FEM) coupled method is proposed for use in numerical investigations of the interaction between a regular wave and a horizontal suspended structure. The fluid domain calculated by the MPS method is dispersed into fluid particles, and the structure domain solved by the FEM method is dispersed into beam elements. The generation of the 2D regular wave is firstly conducted, and convergence verification is performed to determine appropriate particle spacing for the simulation. Next, the regular wave interacting with a rigid structure is initially performed and verified through the comparison with the laboratory experiments. By verification, the MPS-FEM coupled method can be applied to fluid-structure interaction(FSI) problems with waves. On this basis, taking the flexibility of structure into consideration, the elastic dynamic response of the structure subjected to the wave slamming is investigated, including the evolutions of the free surface, the variation of the wave impact pressures, the velocity distribution,and the structural deformation response. By comparison with the rigid case, the effects of the structural flexibility on wave-elastic structure interaction can be obtained.     

VLCC液舱晃荡仿真及结构响应

当船体在波浪中运动的频率与液舱内液体振动的固有频率相近时,舱内液体将会发生剧烈的运动,即晃荡产生了.在船舶结构设计中,必须要考虑这种晃荡载荷及结构的响应.该文针对一超大型油轮(VLCC)液舱晃荡问题进行了三维仿真模拟.在DNV基于压力的晃荡评估方法的基础上,提出了应用有限元方法和任意拉格朗日-欧拉耦合技术进行该船全生命期内液舱晃荡仿真及强度评估的工程实用的步骤方法.     

Ship motion-sloshing interaction using a field method

Even in relatively calm waters, low amplitude wave-induced motions of an LNG carrier may induce large amplitude liquid sloshing inside the ship's partially filled tanks, and the interaction between ship motions and sloshing may affect the ship's seekeeping properties. A computational procedure, here referred to as the RANS-RANS method, was developed to account for this interaction, and this method was then employed to predict the free surface flow inside the tanks and the corresponding motions of the ship in regular head and beam waves. This method coupled a compressible VoF technique with a generic wave generation and absorption scheme to obtain wave-induced ship motions with and without considering the effects of sloshing in the ship's tanks. Systematic grid studies were performed to obtain a sufficiently fine grid needed to yield converging predictions. The resulting wave patterns, ship motions, and internal sloshing elevations were compared with results obtained from a computational method, here referred to as the RANS-BEM method, that relied on a boundary element method to obtain ship motions. This RANS-RANS method was validated against model test measurements.     

Structural safety assessment procedure for membrane-type LNG CCS considering hydroelasticity effect

As environmental regulations have become more onerous, the demand for LNG and LNG carriers has increased. The LNG cargo containment system (LNG CCS) is one of the most important facilities in LNG carriers, and many membrane types of LNG CCS have been developed so far. Traditionally, sloshing model tests are performed and a series of statistical approaches are used to obtain design sloshing pressures. Then, these design loads are utilized to analyze the structural safety of LNG CCS. In the sloshing model test, the flat and almost rigid plate is used as the LNG CCS instead of a real model, and thus, the hydroelasticity effect cannot be considered. In the present research, the effect of hydroelasticity on sloshing pressure was investigated in a fluid-structure interaction simulation. A wet drop simulation was conducted, and its results subsequently were compared with wet drop experimental test results to ensure their validity. Then, two types of structure model, namely a flat-rigid plate model and a flat-flexible CCS model, were generated to investigate the effect of hydroelasticity. Also, a fluid hitting analysis model was devised to realize the sloshing phenomenon, and the two generated structure models were applied as the structure domain in the analysis. In the fluid hitting simulations, it was found that the hydroelasticity effect on sloshing pressure is significant. Thus, correction factors, which are quantitative values of the effect of hydroelasticity, were derived (and are proposed herein). Finally, a structural safety assessment procedure for consideration of the hydroelasticity effect was derived (and is suggested herein).     

基于ALE有限元法二维液体晃荡的数值模拟

随着各类液货船的不断增加,晃荡现象在船舶与海洋工程领域引起了广泛的关注.文章从不可压粘性流体的N-S方程组出发,运用任意拉格朗日(ALE)有限元法对水平激励下的二维液体晃荡现象进行了数值模拟.首先通过计算值与解析解的比较,验证了数值方法的有效性,进而分析了激励振幅、激励频率对晃荡现象的影响,最后计算了作用于壁面上的力和力矩.     

粘性流体晃荡与弹性结构的相互耦合作用

提出了处理粘性流体晃荡与弹性结构相互耦合作用的理论及相应数值计算方法.流体运动采用N-S方程描述,控制方程采用有限差分法离散,并由超松弛迭代法求解.液体自由表面通过流体体积法进行重构.为了考虑液舱结构变形对液体晃荡的影响,应用了FAVOR技术.由此建立了描述液体晃荡与结构相互耦合作用的水弹性理论,并由相应的计算方法进行了系列运算,以考察二维液舱不同结构刚度对液体晃荡的影响.计算结果有助于进一步理解此类问题的物理现象.     

考虑晃荡耦合船舶运动的晃荡压力算法（英文）

The effect of coupling between sloshing and ship motions in the evaluation of slosh-induced interior pressures is studied. The coupling between sloshing loads and ship motions is modelled through a hybrid algorithm which combines a potential flow solution based on transient Green function for the external ship hydrodynamics with a viscous flow solution based on a multiphase interface capturing volume of fluid(VOF) technique for the interior sloshing motion. The coupled algorithm accounts for full nonlinear slosh forces while the external forces on the hull are determined through a blended scheme of linear radiationdiffraction with nonlinear Froude-Krylov and restoring forces. Consideration of this level of nonlinearities in ship motions is found to have non-negligible effects on the slosh-coupled responses and slosh-induced loads. A scheme is devised to evaluate the statistical measure of the pressures through long-duration simulation studies in extreme irregular waves. It is found that coupling significantly influences the tank interior pressures, and the differences in the pressures between coupled and uncoupled cases can be as much as 100% or more. To determine the RAO over the frequency range needed for the simulation studies in irregular waves, two alternative schemes are proposed, both of which require far less computational time compared to the conventional method of finding RAO at each frequency, and the merits of these are discussed.     

A numerical study of the effects of bottom and sidewall stiffeners on sloshing behavior considering roll resonant motion

In the present study, aiming to provide practical guidance for the structural design, the effects of the small bottom and sidewall stiffeners on sloshing loads in shallow and intermediate liquid conditions are investigated numerically. The interaction between the highly nonlinear free surface flow and a large number of small stiffeners is modeled based on the Moving Particle Semi-implicit method. 2D models of a rectangular clean tank and tanks with small sidewalls or bottom stiffeners are considered. For filling ratios ranging from 5% to 20%, harmonic angular motions with periods around the analytical sloshing resonant ones are applied. Simulation snapshots and computed pressure time series, pressure variation, wave run-up, maximum fluid velocity, and kinetic energy were provided to clarify the complex interaction between the small baffles and the shallow free surface flow. As a result, except for the tanks with bottom stiffeners in shallow liquid conditions, in which the response under actual resonant frequency should be considered, the response computed using the clean tank excited by its resonant period can be adopted as a conservative estimation for the sloshing loads of actual tanks with small stiffeners.     

A fully coupled ship motion and sloshing analysis in various container geometries

In the present study a novel modeling approach is presented to solve the combined internal sloshing and sea-keeping problem. The model deals with interesting effects arising due to the coupled interaction between the sloshing in partially filled containers of several geometries and the ship motion. The study is very important for the liquid cargo carrier operating in rough sea or under different environmental conditions. The resulting slosh characteristics that include transient pressure variation, free surface profiles and hydrodynamic pressure over the container walls have been reported in this study. In addition, the effects of coupled ship response and sloshing on ship motion parameters have also been investigated. The equations of motion of fluid, considered inviscid, irrotational, and partially compressible, are expressed in terms of the pressure variable alone. A finite difference-based iterative time-stepping technique is employed to advance the coupled solution in the time domain. Several parameters of interest, including the container parameters, level of liquid, thrusters modeling and some important environmental factors are investigated.     

液舱内三维液体非线性晃荡的数值模拟

晃荡是一种非常复杂的非线性液体流动现象,载液船舶的晃荡问题颇受关注,因为在外界的激励下液舱内会产生剧烈的晃荡现象,巨大的晃荡冲击力会造成结构的破坏.文中建立了三维晃荡数学模型及数值计算模型,借助处理自由表面的VOF(volume of fluid)方法对液舱内液体晃荡的自由表面进行追踪,编制程序实现了液舱内三维液体非线性晃荡的数值模拟,并就三维刚性液舱内粘性液体的自由晃荡和强迫晃荡做了分析,讨论了液体不同粘性系数对晃荡的影响.模拟结果证明了三维晃荡理论的可行性.     

多种激励耦合作用下的二维液舱晃荡(英文)

In this research,liquid sloshing behavior in a 2-D rectangular tank was simulated using ANSYS-FLUENT software subject to single or multiple-coupled external excitations(such as sway coupled with roll,and sway and roll coupled with heave).The volume of fluid(VOF) method was used to track the free surface of sloshing.External excitation was imposed through the motion of the tank by using the dynamic mesh technique.The study shows that if the tank is subjected to multiple coupled excitations and resonant excitation frequencies,liquid sloshing will become violent and sloshing loads,including impact on the top wall,will be intensified.     

升船机提升系统中船厢、水体和钢缆相互耦合作用分析

本文在作者[1]先前采用水弹性理论建立的对整个提升系统的动力特性进行数值研究的理论模型的基础上,首先对船厢变形状态下的缆力表达式作了更为详细的推导,然后以中国船舶科学研究中心的1∶30的缩比模型为对象,从整体角度,对其流固耦合问题,主要是提升过程中船厢内水体的晃荡、提升缆的受力和船厢动态响应,作了数值研究.在提升系统受到船厢内一非对称水波扰动的假定下,分别对船厢简化为刚性和考虑实际弹性的两种情况进行了计算、分析和对比,得出了一些有益的结论.计算结果表明,对升船机这样一个复杂系统,文献[1]中所建立的理论模型及本文的数值计算方法是可行的.     

预应力和流固耦合效应对水下壳结构振动特性影响研究

对于水下结构,流体对结构模态频率的影响主要体现在预应力效应与流固耦合效应两方面.为研究水下壳结构固有频率对这2种效应的敏感度,本文首先计算某水下壳体在空气中的模态,然后研究壳结构在不同水深工作产生的静水压对其固有频率的影响,最后使用声固耦合的方法计算了流固耦合效应下的壳体湿模态.计算结果表明:流固耦合效应占主导因素,但随着壳体工作水深(预应力)的增大,预应力效应对水下壳体频率的影响也逐渐显著.因而,在开展水下结构动力学设计时,不仅要重点考虑流固耦合效应对结构固有频率的影响,还应关注结构实际工作环境下的应力状态,分析其对结构固有频率的影响.     

液舱内大幅晃荡引起的压强预报

本文基于流体体积(VOF)法就部分装载液体的液舱内的晃荡压强进行了二维数值计算.为了模拟大幅晃荡引起的冲击压强,重点对数值计算较为敏感的速度边界条件进行了数值处理.通过对某一矩形液舱内的液体晃荡的计算,得到了自由表面位置和压强的时间历经曲线.计算结果同实验值的比较显示:本方法可以用于计算预报大幅晃荡引起的载荷.     

Numerical simulation of violent sloshing by a CIP-based method

In this article, a new computational fluid dynamics simulation approach based on the constraint interpolation profile (CIP) method is applied to tackle the violent sloshing problem. The present study considers two-dimensional sloshing phenomena in a rectangular tank. By the proposed method, the sloshing problem is viewed as a multiphase problem that includes water and air flows. A stationary Cartesian grid is used and the free surface is solved by an interface capturing method. The CIP combined unified procedure (CCUP) scheme was adopted for the flow solver, and both the CIP scheme and the CIP conservative semi-Lagrangian with cubic interpolation polynomial (CIP-CSL3) scheme were used for interface capturing. For validation of the numerical method, a physical experiment was conducted with a rectangular tank for several frequencies and filling heights. A convergence check was first performed for the method. The numerical simulation results on violent sloshing show that the use of the CIP-CSL3 scheme as an interface capturing procedure gives much better results for the pressures and free-surface profiles than the conventional CIP scheme.     

具有高阶导数自由表面边界条件的两层流体无波势和多级的构建（英文）

There is a large class of problems in the field of fluid structure interaction where higher-order boundary conditions arise for a second-order partial differential equation. Various methods are being used to tackle these kind of mixed boundary-value problems associated with the Laplace's equation(or Helmholtz equation) arising in the study of waves propagating through solids or fluids. One of the widely used methods in wave structure interaction is the multipole expansion method. This expansion involves a general combination of a regular wave, a wave source, a wave dipole and a regular wave-free part. The wave-free part can be further expanded in terms of wave-free multipoles which are termed as wave-free potentials. These are singular solutions of Laplace's equation or two-dimensional Helmholz equation. Construction of these wave-free potentials and multipoles are presented here in a systematic manner for a number of situations such as two-dimensional non-oblique and oblique waves, three dimensional waves in two-layer fluid with free surface condition with higher order partial derivative are considered. In particular, these are obtained taking into account of the effect of the presence of surface tension at the free surface and also in the presence of an ice-cover modelled as a thin elastic plate. Also for limiting case, it can be shown that the multipoles and wave-free potential functions go over to the single layer multipoles and wave-free potential.     

Construction of wave-free potentials and multipoles in a two-layer fluid having free-surface boundary condition with higher-order derivatives

There is a large class of problems in the field of fluid structure interaction where higher-order boundary conditions arise for a second-order partial differential equation. Various methods are being used to tackle these kind of mixed boundary-value problems associated with the Laplace’s equation (or Helmholtz equation) arising in the study of waves propagating through solids or fluids. One of the widely used methods in wave structure interaction is the multipole expansion method. This expansion involves a general combination of a regular wave, a wave source, a wave dipole and a regular wave-free part. The wave-free part can be further expanded in terms of wave-free multipoles which are termed as wave-free potentials. These are singular solutions of Laplace’s equation or two-dimensional Helmholz equation. Construction of these wave-free potentials and multipoles are presented here in a systematic manner for a number of situations such as two-dimensional non-oblique and oblique waves, three dimensional waves in two-layer fluid with free surface condition with higher order partial derivative are considered. In particular, these are obtained taking into account of the effect of the presence of surface tension at the free surface and also in the presence of an ice-cover modelled as a thin elastic plate. Also for limiting case, it can be shown that the multipoles and wave-free potential functions go over to the single layer multipoles and wave-free potential.     

三维完全非线性波浪水槽的数值模拟

用有限元求解拉普拉斯方程,采用满足完全非线性自由表面条件的半拉格朗日法追踪流体自由表面,离散单元采用20节点的六面体二次等参数单元,建立了三维完全非线性数值波浪水槽.把数值计算结果与水面初始升高产生箱体内流体运动解析解和二阶斯托克波理论解进行了对比,结果表明该模型是稳定的、守恒的,能精确模拟非线性波浪的产生和传播,为今后研究非线性波浪对船舶等非规则物体的作用提供了参考.