流体与结构砰击的全非线性模拟:综述（英文）

This paper presents a review of the work on fluid/structure impact based on inviscid and imcompressible liquid and irrotational flow. The focus is on the velocity potential theory together with boundary element method(BEM). Fully nonlinear boundary conditions are imposed on the unknown free surface and the wetted surface of the moving body. The review includes(1) vertical and oblique water entry of a body at constant or a prescribed varying speed, as well as free fall motion,(2) liquid droplets or column impact as well as wave impact on a body,(3) similarity solution of an expanding body. It covers two dimensional(2D), axisymmetric and three dimensional(3D) cases. Key techniques used in the numerical simulation are outlined, including mesh generation on the multivalued free surface, the stretched coordinate system for expanding domain, the auxiliary function method for decoupling the mutual dependence of the pressure and the body motion, and treatment for the jet or the thin liquid film developed during impact.     

水面舰船附近气泡运动数值模拟

当气泡在自由面和水面舰船附近运动时,两者的存在均会对气泡脉动产生影响。基于不可压缩势流理论,采用边界元方法对边界积分方程进行求解。针对建立自由面需要大量的网格,且处理自由面与水面结构交界面时数值不稳定,采用考虑自由面效应的格林函数取代基本格林函数,通过与自由面附近气泡轴对称模型的计算结果进行对比,验证了该方法的有效性。通过对舰船结构附近的气泡运动射流特性进行模拟,发现当药包在舷侧和自由面附近爆炸时,气泡射流可能不会完全作用在舰船上,自由面效应的存在削弱了气泡的打击能力。     

A CIP-based method for numerical simulations of violent free-surface flows

A CFD model is proposed for numerical simulations of extremely nonlinear free-surface flows such as wave impact phenomena and violent wave–body interactions. The constrained interpolation profile (CIP) method is adopted as the base scheme for the model. The wave–body interaction is treated as a multiphase problem, which has liquid (water), gas (air), and solid (wave-maker and floating body) phases. The flow is represented by one set of governing equations, which are solved numerically on a nonuniform, staggered Cartesian grid by a finite-difference method. The free surface as well as the body boundary are immersed in the computation domain and captured by different methods. In this article, the proposed numerical model is first described. Then to validate the accuracy and demonstrate the capability, several two-dimensional numerical simulations are presented, and compared with experiments and with computations by other numerical methods. The numerical results show that the present computation model is both robust and accurate for violent free-surface flows.     

基于无奇异边界元法模拟三维全非线性液舱晃荡

文章基于全非线性势流理论对三维液舱晃荡进行了数值模拟,其控制方程由无奇异边界积分方程法(Desingular-ized Boundary Integral Equation Method,DBIEM)进行离散求解,在求解全非线性的自由面微分方程时,文中采用混合欧拉—拉格朗日法(Mixed Eulerian-Lagrangian,MEL)和四阶Adams-Bashforth-Moulton(ABM4)预报—修正方法,为了避免结果发散即增强数值稳定性,文中采用B样条法来光顺自由面.在微幅水平激励下,该文中得到的结果与解析解吻合较好.     

全非线性方法对于沉浸双圆柱定常前行运动引起自由表面流的研究（英文）

The free surface flow generated by twin-cylinders in forced motion submerged beneath the free surface is studied based on the boundary element method. Two relative locations, namely, horizontal and vertical, are examined for the twin cylinders. In both cases, the twin cylinders are starting from rest and ultimately moving with the same constant speed through an accelerating process. Assuming that the fluid is inviscid and incompressible and the flow to be irrotational, the integral Laplace equation can be discretized based on the boundary element method. Fully-nonlinear boundary conditions are satisfied on the unknown free surface and the moving body surface. The free surface is traced by a Lagrangian technique. Regriding and remeshing are applied, which is crucial to quality of the numerical results. Single circular cylinder and elliptical cylinder are calculated by linear method and fully nonlinear method for accuracy checking and then fully nonlinear method is conducted on the twin cylinder cases, respectively. The generated wave elevation and the resultant force are analysed to discuss the influence of the gap between the two cylinders as well as the water depth. It is found that no matter the kind of distribution, when the moving cylinders are close to each other, they suffer hydrodynamic force with large absolute value in the direction of motion. The trend of force varying with the increase of gap can be clearly seen from numerical analysis. The vertically distributed twin cylinders seem to attract with each other while the horizontally distributed twin cylinders are opposite when they are close to each other.     

带有月池的非对称结构物水波辐射和激励力（英文）

A highly efficient "hybrid integral-equation method" for computing hydrodynamic added-mass, wave-damping, and wave-exciting force of general body geometries with a vertical axis of symmetry is presented. The hybrid method utilizes a numerical inner domain and a semi-infinite analytical outer domain separated by a vertical cylindrical matching boundary.Eigenfunction representation of velocity potential is used in the outer domain; the three-dimensional potential in the inner domain is solved using a "two-dimensional" boundary element method with ring sources and ring dipoles to exploit the body symmetry for efficiency. With proper solution matching at the common boundary, both radiation and diffraction potentials can be solved efficiently while satisfying the far-field radiation condition exactly. This method is applied to compute the hydrodynamic properties of two different body geometries: a vertical-walled moonpool with a bottom plate that restricts the opening and a spar-like structure with a diverging bottom opening inspired by designs of floating Oscillating Water Columns. The effects of the size of the bottom opening on the hydrodynamic properties of the body are investigated for both geometries. The heave motion of the floater as well as the motion of the internal free surface under incident wave excitation are computed and studied for the spar-like structure.     

纵摇容器中液体晃荡的非线性数值模拟

在海洋工程领域,液体晃荡是一种普遍存在的物理现象。对于船舶而言,转动比平动有着更重要的影响。该文针对纵摇容器中的液体晃荡问题,采用高阶边界元方法建立自由水面满足完全非线性边界条件的时域数学模型。通过大地坐标系和随体坐标系之间的坐标变换,使得计算域仅控制在随体坐标系内。求解中采用半混合欧拉—拉格朗日方法追踪流体瞬时水面,运用四阶龙格库塔方法更新下一时间步的波面和速度势。通过与已发表试验和数值结果的对比,验证了建立模型的准确性。进而开展大量数值试验研究容器纵摇运动频率、纵摇转动中心和容器中布置一竖直隔板对晃动波面与荷载的影响。     

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

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

基于全时域势流理论的船舶与液舱晃荡耦合运动的数值计算

基于三维线性有航速时域势流理论计算船体时域运动外域波浪力,同时采用三维全非线性时域势流理论来计算舱内液体的非线性晃荡所诱导力与力矩,进而建立了波浪中载液船舶耦合运动方程。该方法能够完整地考虑波浪、船体和液舱晃荡之间的实时耦合作用。研究结果表明:通过模型实验和数值模拟计算的对比,数值模拟计算能够清晰显现出液舱晃荡对船体全局运动影响,无航速船体运动RAO与模型实验结果吻合良好,有航速运动计算合乎预期。     

基于镜像SPH的三维容器液体晃荡研究

基于SPH方法对液体在三维矩形容器中的非线性晃荡问题进行了数值模拟研究。采用镜像粒子边界建立三维模型,成功解决传统SPH方法在容器壁产生缝隙的问题。对横荡和横摇两种容器运动形式进行数值模拟,在液面时间历程、不同时刻液体形态、点压强以及液体晃荡作用力等方面与其他学者的实验研究结果和数值计算结果进行对比验证,表明此三维模型的正确性。     

浮体升沉和纵摇运动对甲板上浪载荷影响数值研究（英文）

In this paper, the influence of heave and pitch motions on green water impact on the deck is numerically investigated. The vessel motions are determined using a potential theory based method and provided as input to finite volume based CFD computations of green water phenomenon. A dynamic mesh approach is adopted to determine instantaneous body positioning in the fluid domain.Detailed validation studies with published experimental results for 2D and 3D fixed vessel cases are initially performed to validate the present numerical approach before studying the moving vessel problem. The results show that inclusion of heave and pitch motion changes the disturbed wave field near the bow which influences the free surface as well as the impact loading due to green water. The effect of wave steepness on green water impact is also investigated and it is seen that the present numerical method is capable of capturing green water load. It is observed that the effects of vessel motions on green water load are not negligible and one should consider this effect too. The incorporation of vessel motions in the vertical plane affects the green water loading on the deck.     

带有挡板的注液舱室地震响应分析（英文）

In thispaper, the effects of a rigid baffle on the seismic response of liquid in a rigid cylindrical tank are evaluated. A baffle is an annular plate which supplies a kind of passive control on the effects of ground excitation. The contained liquid is assumed incompressible, inviscid and has irrotational motion. To estimate the seismic response, the method of superposition of modes has been applied. To analyze the rigid tank response, Laplace's equation is considered as the governing equation of the fluid domain, in both time and frequency domains. The boundary element method(BEM) is employed to evaluate the natural modes of liquid in a cylindrical tank. To gain this goal, the fluid domain is divided into two upper and lower parts partitioned by the baffle. Linearized kinematic and dynamic boundary conditions of the free surface of the contained liquid have been considered.     

Seismic response of liquid-filled tank with baffles

In thispaper, the effects of a rigid baffle on the seismic response of liquid in a rigid cylindrical tank are evaluated. A baffle is an annular plate which supplies a kind of passive control on the effects of ground excitation. The contained liquid is assumed incompressible, inviscid and has irrotational motion. To estimate the seismic response, the method of superposition of modes has been applied. To analyze the rigid tank response, Laplace＆#39;s equation is considered as the governing equation of the fluid domain, in both time and frequency domains. The boundary element method （BEM） is employed to evaluate the natural modes of liquid in a cylindrical tank. To gain this goal, the fluid domain is divided into two upper and lower parts partitioned by the baffle. Linearized kinematic and dynamic boundary conditions of the free surface of the contained liquid have been considered.     

Numerical investigation of the gate motion effect on a dam break flow

The effects of a moving gate on dam break flow experiments are numerically considered in this paper. The dam break flow process is modelled by a sudden gate release. A Constrained Interpolation Profile (CIP)-based model is employed to solve the Navier-Stokes equations with the free surface boundary condition to deal with the water-air-gate interactions. In addition, the VOF-type THINC/SW method is used to capture the free surface. The movable gate is simulated by an immersed boundary, and a gate motion formula is proposed from existing experimental data and applied in the dam break flow computations. The difference between the water profile and pressure domain with/without gate motion is examined. The dam breaks with different gate motions are discussed as well. The numerical results show that the proposed gate motion formula correlates well with the numerical modelling, and the gate motion has a significant influence in the water collapse process and cannot be neglected in the study of dam break.     

一种计算非线性兴波阻力的基于NURBS的高阶面元法(英文)

文章开发了一种用于求解非线性兴波阻力问题的基于非均匀有理B样条(NURBS)的高阶面元法.该方法采用NURBS基函数表达物面几何和边界面上的源强分布.为了满足非线性自由面边界条件,开发了一个迭代算法;辐射条件通过采用网格错位的数值技术得到满足.以Wigley数学船型为例,在一系列的船速下计算了非线性兴波阻力.数值结果表明所提出的方法对于计算兴波阻力是精确和有效的.     

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

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.     

能量耗散效应的多域边界元法(英文)

The wave diffraction and radiation around a floating body is considered within the framework of the linear potential theory in a fairly perfect fluid.The fluid domain extended infinitely in the horizontal directions but is limited by the sea bed,the body hull,and the part of the free surface excluding the body waterplane,and is subdivided into two subdomains according to the body geometry.The two subdomains are connected by a control surface in fluid.In each subdomain,the velocity potential is described by using the usual boundary integral representation involving Green functions.The boundary integral equations are then established by satisfying the boundary conditions and the continuous condition of the potential and the normal derivation across the control surface.This multi-domain boundary element method(MDBEM) is particularly interesting for bodies with a hull form including moonpools to which the usual BEM presents singularities and slow convergence of numerical results.The application of the MDBEM to study the resonant motion of a water column in moonpools shows that the MDBEM provides an efficient and reliable prediction method.     

简单Green函数法求解三维水动力系数

本文提出用简单Green函数1/r求解无航速的三维物体在水面、水下或在无界域中作各向振荡时水动力系数的数值方法。对于无穷远处扩散条件转换成辐射面条件的问题作了专门的考虑,即用圆柱将物体包围在其内部,该圆柱将流场分成内、外两部分,外部域中速度势的特征函数展开式与内部域中用Green公式求得的解相匹配。 根据本计算方法编制的程序对一些典型情况作了试算,其中包括在水面振荡的半球、方盒和圆盘的水动力系数曲线;水面及水底对圆球水动力系数的影响等。为考核本方法,与其它计算方法得到的结果作了比较。对方盒的三维结果与切片法作了比较。     

Time-domain TEBEM method for mean drift force and moment of ships with forward speed under the oblique seas

A numerical method for solving 3D unsteady potential flow problem of ship advancing in waves is put forward. The flow field is divided into an inner and an outer domain by introducing an artificial matching surface. The inner domain is surrounded by ship wetted surface and matching surface as well as part of the free surface. The free surface condition for the inner domain is formulated by perturbation about the double-body flow or uniform incoming flow assumption. The outer domain is surrounded by matching surface and the rest free surface as well as infinite far-field radiation boundary. The free surface condition for the outer domain is formulated by perturbation about uniform incoming flow. The simple Green function and transient free surface Green function are used to form the boundary integral equation (BIE) for the inner and outer domains, respectively. Taylor Expansion Boundary Element Method (TEBEM) is utilized to solve the double-body flow and inner domain and outer domain unsteady flow BIE. Matching conditions for the inner domain flow and outer domain flow are enforced by the continuity of velocity potential and normal velocity on the matching surface. Direct pressure integration on ship wetted surface is used to obtain the first-order and second-order wave forces (moments). The numerical predictions on the displacement, added resistance, sway mean drift force and yaw mean drift moment of the modified KVLCC2 ship at different forward speeds are investigated by the proposed TEBEM method. It is also compared with the other numerical results. The physical tank experiment results are also developed to validate the accuracy of numerical tank results. Compared with the experiment solutions, a good agreement can be obtained by TEBEM method.

     

桁架式Spar平台与系泊/立管系统的全时域非线性耦合动态分析

开发了对浮式平台系统进行耦合动态分析的全时域程序。采用二阶时域方法计算水动力荷载,在此方法中,对物面边界条件和自由水面边界条件进行泰勒级数展开,利用Stokes摄动展开分别建立相应的一阶、二阶边值问题,而且此边值问题的计算域不随时间变化。采用高阶边界元方法计算每一时刻流场中的速度势,利用四阶预报校正法对二阶自由水面边界条件进行数值积分。在自由表面加入一个人工阻尼层来避免波浪的反射。对于系泊缆索/立管/张力腿的动力分析,在一个总体坐标系中对控制方程进行描述,采用基于细长杆理论的有限元方法进行求解。在耦合动态分析中,采用Newmark方法对平台和系泊缆索/立管/张力腿的运动方程同时进行求解。利用开发的耦合分析程序对一个桁架式Spar平台的运动响应进行了数值模拟,给出了平台的位移和系泊缆索/立管上端点的张力,并得到了一些重要结论。