Numerical simulation of the flow around two-dimensional partially cavitating hydrofoils

In the present study, a new approach is applied to the cavity prediction for two-dimensional （2D） hydrofoils by the potential based boundary element method （BEM）. The boundary element method is treated with the source and doublet distributions on the panel surface and cavity surface by usethe of the Dirichlet type boundary conditions. An iterative solution approach is used to determine the cavity shape on partially cavitating hydrofoils. In the case of a specified cavitation number and cavity length, the iterative solution method proceeds by addition or subtraction of a displacement thickness on the cavity surface of the hydrofoil. The appropriate cavity shape is obtained by the dynamic boundary condition of the cavity surface and the kinematic boundary condition of the whole foil surface including the cavity. For a given cavitation number the cavity length of the 2D hydrofoil is determined according to the minimum error criterion among different cavity lengths, which satisfies the dynamic boundary condition on the cavity surface. The NACA 16006, NACA 16012 and NACA 16015 hydrofoil sections are investigated for two angles of attack. The results are compared with other potential based boundary element codes, the PCPAN and a commercial CFD code （FLUENT）. Consequently, it has been shown that the results obtained from the two dimensional approach are consistent with those obtained from the others.     

二维水翼的设计与优化

通过指定沿叶剖面的压力分布或速度分布来设计翼剖面是水翼设计优化的有效途径.在给定二维水翼表面压力分布的条件下,通过已知物体表面的扰动势分布,利用格林定理,建立求解物体表面法向速度的Dirichlet型边值方程.假定一个初始剖面,通过积分得到扰动势分布,利用面元法求解方程得到法向速度.物面边界条件用于剖面的形状修正,利用迭代计算方法可达到设计目标.通过调整压力分布,可对翼型进行优化.文中对对称翼及不对称翼进行设计,并优化设计了剖面A, B和C.     

有限水深三维脉动源格林函数数值算法研究

由于柯西主值积分的奇异性和贝塞尔函数的振荡性,有限水深情况下复杂格林函数及其导数的精确数值求解一直是浅水中波浪水动力计算的难点,因此寻找格林函数的精确数值解显得非常重要.通过对格林函数奇异项进行变形推导,文中给出了一种去掉了奇点的新积分形式.另外通过改进前人推导的Gauss-Laguerre积分方法,给出了一种改进的新Gauss-Laguerre积分方法.格林函数及其导数的数值结果显示文中给出的两种新方法可以有效地计算复杂格林函数及其导数值.最后对这两种方法、级数解以及传统的Gauss-Laguerre积分方法的计算精度和效率进行了比较研究,结果显示文中给出的两种方法的计算精度高于传统的Gauss-Laguerre积分方法,但其计算效率低于级数解.但在接近于0的近场附近级数解失真,而文中提出的改进的新Gauss-Laguerre积分方法可以获得正确结果.同时当频率和水深均较大时级数解失真,而文中提出的方法也可以获得正确结果.最后针对这些方法的优缺点,该文提出了建议的策略用于计算有限水深格林函数.     

水深测量数据的组织管理与三维可视化

文章应用GIS理论思想,对水深测量数据格式进行深入分析,提出应用空间数据库来组织管理水深数据的新方法,给出了数据结构。在此基础上,对水深测量数据的应用价值进行深入挖掘,提出利用水深离散点数据进行水下地形建模,从而实现水深数据的三维可视化直观表达,并就此给出了地形建模过程和渲染步骤以及三维场景效果优化方法。     

The applicability of hydrofoils as a ship control device

Centrifugal forces are commonly created when ships turn, which may cause a ship to capsize in a critical situation. A mathematical model has been developed to optimize the stability coefficients for ship, with the aim to prevent capsizing and to increase ship maneuverability in high-speed water craft. This model can be used to develop algorithms for control system improvement. The mathematical model presented in this paper optimized the use of multipurpose hydrofoils to reduce heeling and the trimming moment, maintaining an upright ship’s position and lessening the resistance via transverse force. Conventionally, the trimming and heeling of a ship are controlled using ballast water; however, under variable sea conditions it is sometimes difficult to control a ship’s motion using ballast water. In this case, a hydrofoil would be more stable and maneuverable than a ballast tank controlled vessel. A movable hydrofoil could theoretically be adapted from moveable aerofoil technology. This study proves the merit of further investigation into this possibility.     

热力学效应对非定常空化流动结构影响的实验研究

为了研究热力学效应对空化流动结构的影响,采用实验的方法研究了水温为6℃、25℃和45℃时绕水翼的非定常空化流场。利用高速全流场显示技术,观察了不同水温下绕水翼的空化流动形态,并利用数字粒子图像测速仪(DPIV)测量了其速度和涡量分布。结果表明:随着水温的升高,空穴区域减小,空穴长度减短。在相同的空化阶段,热力学效应对空泡的脱落周期影响较小,脱落周期基本不变。热力学效应对速度和涡量影响较大的区域主要集中在空化区域及其与主流区的交界面处。随着水温的升高,低速高脉动区域逐渐减小,且对应的速度值略有升高,导致速度梯度减小,湍流脉动强度降低。同时发现,涡量区域对应于流场中具有速度梯度的区域,上下涡量区域随着水温的升高减小,涡量值降低。     

二维水翼型空化流的数值计算

In order to predict the effects of cavitation on a hydrofoil, the state equations of the cavitation model were combined with a linear viscous turbulent method for mixed fluids in the computational fluid dynamics (CFD) software FLUENT to simulate steady cavitating flow. At a fixed attack angle, pressure distributions and volume fractions of vapor at different cavitation numbers were simulated, and the results on foil sections agreed well with experimental data. In addition, at the various cavitation numbers, the vapor fractions at different attack angles were also predicted. The vapor region moved towards the front of the airfoil and the length of the cavity grew with increased attack angle. The results show that this method of applying FLUENT to simulate cavitation is reliable.     

柴油机主轴承盖应力和应变的3D有限元分析

本文应用AutoCAD软件对柴油机悬挂式主轴承盖进行三维实体造型,用Auto3Dmesh程序对三维实体进行自动网格划分,并进行适当网格调整,用SAPP93软件及自编程序进行了载荷、约束处理和三维有限元计算,并对计算结果进行了可视化后处理,取得了令人满意的结果。     

天津港二突堤转角区域码头结构与土体三维有限元分析

近年码头调查发现,天津港多个突堤转角区域的码头结构变形和横向错位情况十分严重,这与以前发现和认识的规律有所不同。文中从这些问题入手,以天津港二突堤为对象,系统地研究了该突堤转角区域码头结构和岸坡相互作用机理和变形特点,并指出突堤转角区域的码头结构受力比较复杂,其变形具有空间变化特点,传统的二维分析和设计方法是不够全面的。     

三维波对船舶运动响应影响研究

通常情况,实际海况的波浪更接近三维波,而在船舶运动预报时,目前主要考虑二维波(长峰波)中的响应,不能真正反映船舶在实际海浪中的运动.文中利用不同方法(频域方法、时域方法以及模型)研究了三维波的方向扩散性对船舶运动响应的影响.三种方法得到的三维波对船舶运动的影响是一致的.结果显示,在研究实际海况中的船舶运动时,考虑三维波的方向扩散性是非常重要的.     

时域有限水深格林函数及其导数的数值计算

精确快速地求取有限水深格林函数及它的导数是应用时域法求解有限水深水动力学问题的关键.通过对有限水深格林函数在不同区域渐近特性的分析,给出了使用多维Chebyshev多项式和渐近展开式快速近似计算有限水深格林函数及它的导数的方法.近似方法与直接积分方法的比较表明近似方法的计算工作量大幅减少,可以给出足够精度的计算结果.     

New method for ship finite element method preprocessing based on a 3D parametric technique

A new method for ship finite element method (FEM) preprocessing is presented as well as its program development. The method is applicable for all kinds of ships at different levels, such as a whole ship, cargo hold parts or detailed structures. The 3D parametric technique is used when creating ship structures, which improves the modeling efficiency greatly and makes the model easy to modify. A 3D geometric constraint solver is developed to solve the constraint system of the parametric model. A meshing procedure is presented to automatically convert the parametric structure model into a finite element model, by which high quality mesh is generated in the stress concentrated area. It also becomes possible to create finite element models for different levels from the same structure model. Using this method, the engineers avoid much of the complex and laborious work of FEM preprocessing, which consumes a very significant amount of time in finite element analysis, and can pay more attention to post-processing. This method has proved to be practical and highly efficient by several engineering trials. This work is sponsored by “Liaoning BaiQianWan Talents Program”.     

On prediction of longitudinal attitude of planing craft based on controllable hydrofoils

The purpose of this research study was to examine the attitude response of a planing craft under the controllable hydrofoils.Firstly,a non-linear longitudinal attitude model was established.In the mathematical model,effects of wind loads were considered.Both the wetted length and windward area varied in different navigation conditions.Secondly,control strategies for hydrofoils were specified.Using the above strategies,the heave and trim of the planing craft was adjusted by controllable hydrofoils.Finally,a simulation program was developed to predict the longitudinal attitudes of the planing craft with wind loads.A series of simulations were performed and effects of control strategies on longitudinal attitudes were analyzed.The results show that under effects of wind loads,heave of fixed hydrofoils planing craft decreased by 6.3%,and pitch increased by 8.6% when the main engine power was constant.Heave decreased by less than 1% and trim angle decreased by 1.7% as a result of using variable attack angle hydrofoils;however,amplitude changes of heave and pitch were less than 1% under the control of changeable attack angle hydrofoils and longitudinal attitude.     

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

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

The effect of sloshing on the sway motions of 2D rectangular cylinders in regular waves

In this paper, we investigated the effect of sloshing on the sway motions of two-dimensional rectangular cylinders in regular waves, bearing in mind possible applications for LNG-FPSO and LNG-FSRU. First, we carried out experiments for two models with different drafts, or the same draft but different filling ratios, in which the models were firmly connected to each other. The sway motion was measured with a noncontact video camera. This is an extension of Rognebakke and Faltinsen’s work for a single model (J Ship Res 47(3):208–221, 2003). It was found that the sway motion became small when the incident wave frequency was close to the lowest natural frequency of each model. The sway motion greatly increased when the wave frequency was higher than this frequency. The measured data were compared with numerical results obtained by a single-dominant multi-modal method; relatively good agreement was noted. However, the numerical results deviated from the experimental results near the lowest natural frequency of the smaller model, which was believed to be due to overturning waves, as observed during the experiment. Since this is out of the valid range for the single-dominant multi-modal method, other, more appropriate, methods such as the multi-dominant modal method must be applied instead.     

A 3D, finite element model for baroclinic circulation on the Vancouver Island continental shelf

This paper describes the development and application of a 3-dimensional model of the barotropic and baroclinic circulation on the continental shelf west of Vancouver Island, Canada. A previous study with a 2D barotropic model and field data revealed that several tidal constituents have a significant baroclinic component (the K1 in particular). Thus we embarked on another study with a 3D model to study the baroclinic effects on the residual and several selected tidal constituents.The 3D model uses a harmonic expansion in time and a finite element discretization in space. All nonlinear terms are retained, including quadratic bottom stress, advection and wave transport (continuity nonlinearity). The equations are solved as a global and a local problem, where the global problem is the solution of the wave equation formulation of the shallow water equations, and the local problem is the solution of the momentum equation for the vertical velocity profile. These equations are coupled to the advection-diffusion equation for density so that density gradient forcing is included in the momentum equations. However, the study presented here describes diagnostic calculations for the baroclinic residual circulation only.The model is sufficiently efficient that it encourages sensitivity testing with a large number of model runs. In this sense, the model is akin to an extension of analytical solutions to the domain of irregular geometry and bottom topography where this parameter space can be explored in some detail.In particular, the consequences of the sigma coordinate system used by the model are explored. Test cases using an idealized representation of the continental shelf, shelf break and shelf slope, lead to an estimation of the velocity errors caused by interpolation errors inherent in the sigma coordinate system. On the basis of these estimates, the computational grid used in the 2D model is found to have inadequate resolution. Thus a new grid is generated with increased accuracy in the region of the shelf break. However, even with increased resolution, spurious baroclinic circulation seaward of the shelf break and in the vicinity of Juan de Fuca canyon remained a significant problem when the pressure gradient terms were evaluated using the σ coordinate system and using a realistic density profile.With the new grid, diagnostic calculations of the barotropic and baroclinic residual circulation are performed using forcing from the observed σ_t (density) field and from the gradient of this field.     

Scattering of surface water waves involving semi-infinite floating elastic plates on water of finite depth

Two problems of scattering of surface water waves involving a semi-infinite elastic plate and a pair of semi-infinite elastic plates,separated by a gap of finite width,floating horizontally on water of finite depth,are investigated in the present work for a two-dimensional time-harmonic case.Within the frame of linear water wave theory,the solutions of the two boundary value problems under consideration have been represented in the forms of eigenfunction expansions.Approximate values of the reflection and transmission coefficients are obtained by solving an over-determined system of linear algebraic equations in each problem.In both the problems,the method of least squares as well as the singular value decomposition have been employed and tables of numerical values of the reflection and transmission coefficients are presented for specific choices of the parameters for modelling the elastic plates.Our main aim is to check the energy balance relation in each problem which plays a very important role in the present approach of solutions of mixed boundary value problems involving Laplace equations.The main advantage of the present approach of solutions is that the results for the values of reflection and transmission coefficients obtained by using both the methods are found to satisfy the energy-balance relations associated with the respective scattering problems under consideration.The absolute values of the reflection and transmission coefficients are presented graphically against different values of the wave numbers.     

水深对沙波紊流特性影响的数值分析

运用FLUENT软件雷诺应力紊流模型(RSM)及SIMPLEC算法,对不同水深下沙波紊流特性进行数值模拟分析,模拟结果与实测结果吻合良好。在此基础上对4种不同相对水深下纵、垂向流速、压强及雷诺切应力分布规律进行研究。研究结果表明:相对水深较大时紊流特性受沙波地形影响较小,只有近底水体受沙波作用;相对水深较小时整个水体都会受到影响。在沙波尺度不变的条件下,紊流强度一般随水深的增加而减小,不同位置处的沙波紊流特性受水深影响程度及变化规律均不同。最大压强水平位置受水深影响的变化规律与再附点基本一致,可将最大压强作为再附点位置的判断依据。     

一种二维态势显示与三维视景一体化系统设计与实现

在介绍当前二维态势编辑软件与三维视景仿真软件优缺点的基础上,提出一种二维态势编辑显示与三维视景显示一体化的思想来解决独立的二维态势与独立的三维视景显示的不足之处,有机地结合了二维态势的宏观性、整体性、抽象性和三维视景的现实性、直观性的优点,克服了二维态势的抽象多义性和三维虚拟场景漫游的迷失感,并且编程实现该设想.     

VB与TrueGRID软件在船闸闸首三维有限元网格自动剖分中的应用

利用VB语言设计的友好的人机交互界面,获取船闸闸首模型的参数化数据并对数据进行分析处理,生成TrueGRID软件能识别的命令流;编制VB与TmeGRID的接口程序,导入命令流至TrueGRID软件,生成可视化的船闸闸首三维有限元网格.提高了有限元计算前处理的工作效率.