采用高阶边界元法的楔形体变速度入水数值模拟(英文)

A high order boundary element method was developed for the complex velocity potential problem. The method ensures not only the continuity of the potential at the nodes of each element but also the velocity. It can be applied to a variety of velocity potential problems. The present paper, however, focused on its application to the problem of water entry of a wedge with varying speed. The continuity of the velocity achieved herein is particularly important for this kind of nonlinear free surface flow problem, because when the time stepping method is used, the free surface is updated through the velocity obtained at each node and the accuracy of the velocity is therefore crucial. Calculation was made for a case when the distance S that the wedge has travelled and time t follow the relationship s=Dtα, where D and α are constants, which is found to lead to a self similar flow field when the effect due to gravity is ignored.     

Numerical simulation of the water entry of a structure in free fall motion

To solve the problems concerning water entry of a structure, the RANS equations and volume of fluid （VOF） method are used. Combining the user-defined function （UDF） procedure with dynamic grids, the water impact on a structure in free fall is simulated, and the velocity, displacement and the pressure distribution on the structure are investigated. The results of the numerical simulation were compared with the experimental data, and solidly consistent results have been achieved, which validates the numerical model. Therefore, this method can be used to study the water impact problems of a structure.     

Wave load computation in direct strength analysis of semi-submersible platform structures

A wave load computation approach in direct strength analysis of semi-submersible platform structures was presented in this paper. Considering the differences in shape of pontoon, column and beam, the combination of accumulative chord length cubic parameter spline theory and analytic method was adopted for generating the wet surface mesh of platform. The hydrodynamic coefficients of platform were calculated by the three-dimensional potential flow theory of the linear hydrodynamic problem for platform with low forward speed. The equation of platform motions was established and solved in frequency domain, and the responses of wave-induced loads on the platform can be obtained. With the interpolation method being utilized, the pressure loads on shell elements for finite element analysis(FEA) were converted from those on the hydrodynamic computation mesh, which pave the basis for FEA with commercial software. A computer program based on this method has been developed ,and a calculation example of semi-submersible platform was illustrated. Analysis results show that this method is a satisfying approach of wave loads computation for this kind of platform.     

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

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.     

一种水翼空化的被动控制方法（英文）

This paper introduces a new idea of controlling cavitation around a hydrofoil through a passive cavitation controller called artificial cavitation bubble generator(ACG). Cyclic processes, namely, growth and implosion of bubbles around an immersed body, are the main reasons for the destruction and erosion of the said body. This paper aims to create a condition in which the cavitation bubbles reach a steady-state situation and prevent the occurrence of the cyclic processes. For this purpose, the ACG is placed on the surface of an immersed body, in particular, the suction surface of a 2 D hydrofoil. A simulation was performed with an implicit finite volume scheme based on a SIMPLE algorithm associated with the multiphase and cavitation model. The modified k-ε RNG turbulence model equipped with a modification of the turbulent viscosity was applied to overcome the turbulence closure problem. Numerical simulation of water flow over the hydrofoil equipped with the ACG shows that a low-pressure recirculation area is produced behind the ACG and artificially generates stationary cavitation bubbles. The location, shape, and size of this ACG are the crucial parameters in creating a proper control. Results show that the cavitation bubble is controlled well with a well-designed ACG.     

3-D computational method of wave loads on turret moored FPSO tankers

A three-dimensional method of calculating wave loads of turret moored FPSO (Floating Production Storage and Offloading) tankers is presented. The linearized restoring forces acting on the ship hull by the mooring system are calculated according to the catenary theory, which are expressed as the function of linear stiffness coefficients and the displacements of the upper ends of mooring chains. The hydrodynamic coefficients of the ship are calculated by the three-dimensional potential flow theory of the linear hydrodynamic problem for ships with a low forward speed. The equations of ship motions are established with the effect of the restoring forces from the mooring system included as linear stiffness coefficients. The equations of motions are solved in frequency domain, and the responses of wave-induced motions and loads on the ship can be obtained. A computer program based on this method has been developed, and some calculation examples are illustrated. Analysis results show that the method can give satisfying prediction of wave loads.     

Numerical computation of motions and structural loads for large containership using 3D Rankine panel method

In this paper, we present the results of our numerical seakeeping analyses of a 6750-TEU containership, which were subjected to the benchmark test of the 2 nd ITTC–ISSC Joint Workshop held in 2014. We performed the seakeeping analyses using three different methods based on a 3D Rankine panel method, including 1) a rigid-body solver, 2) a flexible-body solver using a beam model, and 3) a flexible-body solver using the eigenvectors of a 3D Finite Element Model(FEM). The flexible-body solvers adopt a fully coupled approach between the fluid and structure. We consider the nonlinear Froude–Krylov and restoring forces using a weakly nonlinear approach. In addition, we calculate the slamming loads on the bow flare and stern using a 2D generalized Wagner model. We compare the numerical and experimental results in terms of the linear response, the time series of the nonlinear response, and the longitudinal distribution of the sagging and hogging moments. The flexible-body solvers show good agreement with the experimental model with respect to both the linear and nonlinear results, including the high-frequency oscillations due to springing and whipping vibrations. The rigid-body solver gives similar results except for the springing and whipping.     

Numerical simulation of nonlinear wave force on a quasi-ellipse caisson

A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume of fluid (VOF) method was employed to trace the free surface. The partial cell method was used to deal with the irregular boundary typical of this type of problem during first-time wave interaction with the structure, and a satisfactory result was obtained. The numerical model was verified and used to investigate the effects of the relative wave height H/d, relative caisson width kD, and relative length-width ratio B/D on the wave forces of the quasi-ellipse caisson. It was shown that the relative wave height H/d has a significant effect on the wave forces of the caisson. Compared with the non-dimensional inline wave force, the relative length-width ratio B/D was shown to have significant influence on the non-dimensional transverse wave force.     

舰船气泡尾流数值模拟(英文)

Based on a volume of fluid two-phase model imbedded in the general computational fluid dynamics code FLUENT6.3.26, the viscous flow with free surface around a model-scaled KRISO container ship (KCS) was first numerically simulated. Then with a rigid-lid-free-surface method, the underwater flow field was computed based on the mixture multiphase model to simulate the bubbly wake around the KCS hull. The realizable k-ε two-equation turbulence model and Reynolds stress model were used to analyze the effects of turbulence model on the ship bubbly wake. The air entrainment model, which is relative to the normal velocity gradient of the free surface, and the solving method were verified by the qualitatively reasonable computed results.     

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.     

A comparison of hydrodynamic impacts prediction methods with two dimensional drop test data

The effect of hydrodynamic impacts on a ship's bow structure is an important design consideration. In addition, to possible failure to local structure, such impacts can generate slam-induced whipping loads, the magnitude of which when combined with ordinary wave induced loads, could lead to first passage failure of a ship's hull girder. This paper presents the results of a recent study that compares select methods of predicting hydrodynamic impacts with available test data of 2-D wedge shaped section shapes. Results based on peak pressures at different impact velocities are presented. The results are compared with traditional prediction theories of Wagner and Chuang.     

Wave impact underneath horizontal decks

The problem of water impact on a fixed horizontal platform deck from regular incident waves was studied. Two-dimensional potential flow was assumed, and the resulting boundary-value problem was solved by three alternative numerical methods, a method based on a generalization of the impact theory by Wagner, and two different nonlinear boundary-element methods. The Wagner-based method used a von Karman approach during the water exit phase, i.e., when the wetted surface decreases. Experiments of the impact on an idealized platform deck have been performed to validate the theory. Comparisons show that the Wagner-based method yields good results for the water entry phase, when the wetted deck area increases, but poor results when the wetted area diminishes. The boundary-element methods compare well with experiments for the entire impact process. A Kutta condition is necessary at the aft body–free surface intersection during water exit and when the flow separates from the aft edge. Gravity effects matter for the water exit phase.     

二维物体入水砰击问题的理论方法研究

入水问题的研究在水动力学和海洋工程等领域有着重要的现实意义。文章采用Wagner模型的思想,将其中的相当平板理论改为椭圆拟合,得到了不同斜升角的二维楔形体匀速入水时的湿表面无量纲压力分布和砰击力,并与其它理论方法和数值计算结果相比较。椭圆拟合方法得到的压力分布在楔形体斜升角较小时与数值计算结果吻合良好,砰击力与数值计算结果在较大斜升角的情况下仍然很接近,比其它理论方法的适用范围更广,可以作为工程上估算砰击力大小的一种新方法。     

基于显式有限元方法的二维楔形体砰击压力系数预报

利用LS-DYNA软件对二维楔形体入水问题进行研究。基于显式有限元方法,选用任意拉格朗日-欧拉算法,建立了包含空气、水和楔形体的完全耦合的二维有限元模型,研究流场的射流现象与压力变化情况,预报了二维楔形体砰击压力系数,并与已公开发表的模型试验结果吻合较好,从而为后续砰击载荷计算提供可靠的方法。     

波流作用下圆柱体入水特性的三维数值模拟研究

结构在实际海况中入水受到多种载荷的共同作用,同时还伴随波浪作用的影响,因此该过程是一个强非线性的过程。针对结构在波、流中入水过程的特点,将入射波(波、流)引入非线性双渐进法,研究三维刚体圆柱体在波、流及波流联合作用下入水过程运动响应及姿态的变化,计算结果与试验结果符合得较好,非线性双渐进法适用于分析三维刚体波、流中入水问题。结果表明在近波面附近,结构受波浪作用明显,入水相位、浪级、流速及波流速度矢量差异对结构入水运动速度及轨迹影响明显。     

考虑定常兴波影响的辐射问题去奇Rankine源方法数值计算

随着中高速船舶在海洋资源开发中的广泛应用,航速对船舶在波浪中摇荡运动的影响受到重视。本文基于势流理论,采用非均匀有理B样条(NURBS)表达物面,应用去奇Rankine源方法求解流场速度势;推导了兴波问题的二阶近似非线性自由面边界条件,以及考虑二阶近似非线性兴波影响的辐射速度势的线性自由面和物面边界条件。本文编程首先数值计算了Wigley船型的近似非线性兴波问题,将计算的线性和非线性兴波阻力系数与文献结果做了比较;然后数值计算了计及线性定常兴波影响的下潜圆球有航速辐射问题,将计算的附加质量和阻尼系数与文献结果比较。本文计算结果与文献结果吻合程度较好,验证了本文数值计算方法的有效性。文中还给出了兴波波形和辐射波形。     

Hydroelastic analysis of a rectangular plate subjected to slamming loads

A hydroelastic analysis of a rectangular plate subjected to slamming loads is presented. An analytical model based on Wagner theory is used for calculations of transient slamming load on the ship plate. A thin isotropic plate theory is considered for determining the vibration of a rectangular plate excited by an external slamming force. The forced vibration of the plate is calculated by the modal expansion method. Analytical results of the transient response of a rectangular plate induced by slamming loads are compared with numerical calculations from finite element method. The theoretical slamming pressure based on Wagner model is applied on the finite element model of a plate. Good agreement is obtained between the analytical and numerical results for the structural deflection of a rectangular plate due to slamming pressure. The effects of plate dimension and wave profile on the structural vibration are discussed as well. The results show that a low impact velocity and a small wetted radial length of wave yield negligible effects of hydroelasticity.     

楔形与船首外飘剖面的砰击研究

文章通过分析MLM解析砰击模型,提出了一种更适用于大底升角楔形体砰击的半解析砰击模型。通过采用数值仿真技术对外飘剖面入水问题的研究,指出了现有规范方法存在的不足,并扩展了半解析方法来实现外飘剖面砰击压力峰值系数的预报。在楔形剖面与船首外飘剖面的入水砰击研究中,对比了半解析砰击模型和数值仿真技术的计算结果,验证了半解析方法的可行性。