基于改进的Logvinovich法对带横摇运动非对称楔形体入水砰击力的计算(英文)

The water entry problem of an asymmetric wedge with roll motion was analyzed by the method of a modified Logvinovich model (MLM). The MLM is a kind of analytical model based on the Wagner method, which linearizes the free surface condition and body boundary condition. The difference is that the MLM applies a nonlinear Bernoulli equation to obtain pressure distribution, which has been proven to be helpful to enhance the accuracy of hydrodynamic loads. The Wagner condition in this paper was generalized to solve the problem of the water entry of a wedge body with rotational velocity. The comparison of wet width between the MLM and a fully nonlinear numerical approach was given, and they agree well with each other. The effect of angular velocity on the hydrodynamic loads of a wedge body was investigated.     

Free Surface Flow Generated by Submerged Twin-cylinders in Forward Motion Using a Fully Nonlinear Method

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.     

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

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.     

液舱晃荡与浮体非线性耦合运动分析(英文)

Nonlinear interactions among incident wave, tank-sloshing and floating body coupling motion are investigated. The fully nonlinear sloshing and body-surface nonlinear free surface hydrodynamics is simulated using a Non-Uniform Rational B-Spline (NURBS) higher-order panel method in time domain based on the potential theory. A robust and stable improved iterative procedure (Yan and Ma, 2007) for floating bodies is used for calculating the time derivative of velocity potential and floating body motion. An energy dissipation condition based on linear theory adopted by Huang (2011) is developed to consider flow viscosity effects of sloshing flow in nonlinear model. A two-dimensional tank model test was performed to identify its validity. The present nonlinear coupling sway motion results are subsequently compared with the corresponding Rognebakke and Faltinsen (2003)’s experimental results, showing fair agreement. Thus, the numerical approach presented in this paper is expected to be very efficient and realistic in evaluating the coupling effects of nonlinear sloshing and body motion.     

Time-domain simulation for water wave radiation by floating structures （Part A）

Direct time-domain simulation of floating structures has advantages： it can calculate wave pressure fields and forces directly; and it is useful for coupled analysis of floating structures with a mooring system. A time-domain boundary integral equation method is presented to simulate three-dimensional water wave radiation problems. A stable form of the integration free-surface boundary condition （IFBC） is used to update velocity potentials on the free surface. A multi-transmitting formula （MTF） method with an artificial speed is introduced to the artificial radiation boundary （ARB）. The method was applied to simulate a semi-spherical liquefied natural gas （LNG） carrier and a semi-submersible undergoing specified harmonic motion. Numerical parameters such as the form of the ARB, and the time and space discretization related to this method are discussed. It was found that a good agreement can be obtained when artificial speed is between 0.6 and 1.6 times the phase velocity of water waves in the MTF method. A simulation can be done for a long period of time by this method without problems of instability, and the method is also accurate and computationally efficient.     

船体周围Stokes波浪流耦合VOF模型的数值模拟（英文）

The surface wave generated by flow around a ship hull moving near free surface of water is simulated numerically in this study. The three-dimensional implicit finite volume method(FVM) is applied to solve Reynolds averaged Navier-Stokes(RANS) equation. The realizable k-ε turbulence model has been implemented to capture turbulent flow around the ship hull in the free surface zone. The volume of fluid(VOF) method coupled with the Stokes wave theory has been used to determine the free surface effect of water. By using is a six degrees of freedom model, the ship hull's movement is numerically solved with the Stokes wave together. Under the action of Stokes waves on the sea, the interface between the air and water waves at the same regular pattern and so does the pressure and the vertical velocity. The ship hull moves in the same way as the wave. The amplitude of the ship hull's heave is less than the wave height because of the viscosity damping. This method could provide an important reference for the study of ships' movement, wave and hydrodynamics.     

Fully nonlinear simulation for fluid/structure impact: A review

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.     

Numerical Simulation of the Stokes Wave for the Flow around a Ship Hull Coupled with the VOF Model

The surface wave generated by flow around a ship hull moving near free surface of water is simulated numerically in this study. The three-dimensional implicit finite volume method(FVM) is applied to solve Reynolds averaged Navier-Stokes(RANS) equation. The realizable k-ε turbulence model has been implemented to capture turbulent flow around the ship hull in the free surface zone. The volume of fluid(VOF) method coupled with the Stokes wave theory has been used to determine the free surface effect of water. By using is a six degrees of freedom model, the ship hull’s movement is numerically solved with the Stokes wave together. Under the action of Stokes waves on the sea, the interface between the air and water waves at the same regular pattern and so does the pressure and the vertical velocity. The ship hull moves in the same way as the wave. The amplitude of the ship hull’s heave is less than the wave height because of the viscosity damping. This method could provide an important reference for the study of ships’ movement, wave and hydrodynamics.     

三维水翼升力线理论(英文)

Prandtl’s lifting line theory was generalized to the lifting problem of a three-dimensional hydrofoil in the presence of a free surface. Similar to the classical lifting theory, the singularity distribution method was utilized to solve two-dimensional lifting problems for the hydrofoil beneath the free surface at the air-water interface, and a lifting line theory was developed to correct three-dimensional effects of the hydrofoil with a large aspect ratio. Differing from the classical lifting theory, the main focus was on finding the three-dimensional Green function of the free surface induced by the steady motion of a system of horseshoe vortices under the free surface. Finally, numerical examples were given to show the relationship between the lift coefficient and submergence Froude numbers for 2-D and 3-D hydrofoils. If the submergence Froude number is small free surface effect will be significant registered as the increase of lift coefficient. The validity of these approaches was examined in comparison with the results calculated by other methods.     

A parametric study on water-entry of a twin wedge by boundary element method

High speed catamarans are used for pleasure, racing as well as passenger transportations. Optimal design of these crafts requires knowledge of sea loads exerted on their structures. The total load may be estimated by integration of loads exerted on a series of two-dimensional sections along the hull. In order to access the cross-sectional loads, the problem may be simplified to solve the water-entry problem of a twin hull. In this paper, water-entry problem of a twin wedge at constant vertical water-entry speed is studied. The problem is solved in the framework of potential theory using boundary element method where gravity effect on the flow is neglected. A simplified model based on Wagner theory is employed. Free surface elevation and pressure distribution on the body in different deadrise angles have been evaluated. A parametric study has been done to investigate effects of deadrise angle, distance between demi-hulls and free surface elevation on maximum pressure coefficient. Finally, a regression formula for maximum pressure coefficient has been proposed. Results of parametric study reveal that as time advances the interaction between two demi-hull gets more severe, besides the interaction effect on pressure coefficient is nonlinear.     

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.

     

采用THOBEM方法研究波流与结构物相互作用的数值模拟(英文)

To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method (HOBEM). By assuming small flow velocities, the velocity potential could be expressed for linear and higher order components by perturbation expansion. A 4th-order Runge-Kutta method was applied for time marching. An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves. Validation of the numerical method was carried out on run-up, wave exciting forces, and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder. The results were in close agreement with the results of a frequency-domain method and a published time-domain method. The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.     

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

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

带自由面流体运动的拉格朗日有限元分析

介绍了拉格朗日描述下的非定常不可压缩带自由面流体运动的有限元分析.在时间积分中采用速度修正法,采用速度修正分步法不仅可以让速度和压力均采用同阶插值函数,而且可以使算法变得更简单.用四边形单元对所求区域加以划分,并借助于Galerkin加权余量法导出相应的有限元方程组.通过实例计算,阐明了用该方法来分析带自由面流体运动的有效性及其实用性.     

Two-dimensional numerical simulation of an elastic wedge water entry by a coupled FDM-FEM method

Hydroelastic behavior of an elastic wedge impacting on calm water surface was investigated. A partitioned approach by coupling finite difference method (FDM) and finite element method (FEM) was developed to analyze the fluid structure interaction (FSI) problem. The FDM, in which the Constraint Interpolation Profile (CIP) method was applied, was used for solving the flow field in a fixed regular Cartesian grid system. Free surface was captured by the Tangent of Hyperbola for Interface Capturing with Slope Weighting (THINC/SW) scheme. The FEM was applied for calculating the structural deformation. A volume weighted method, which was based on the immersed boundary (IB) method, was adopted for coupling the FDM and the FEM together. An elastic wedge water entry problem was calculated by the coupled FDM-FEM method. Also a comparison between the current numerical results and the published results indicate that the coupled FDM-FEM method has reasonably good accuracy in predicting the impact force.     

江海直达船船首入水砰击研究

基于有限体积法的数值仿真方法对3种不同江海直达船船首的入水砰击问题进行了研究。建立包含气、水流场的三维有限元模型,研究了在自由液面的变化、空气层的作用下、不同垂向速度时3种船首入水砰击问题。结果表明,数值计算仿真的方法可以有效地模拟江海直达船船首入水过程中自由液面和空气层的变化情况。在此方法下,相同速度时西瓜首的砰击峰值优于其他船首形式。     

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

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

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

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