A potential based panel method for 2-D hydrofoils

A potential based panel method for the hydrodynamic analysis of 2-D hydrofoils moving beneath the free surface with constant speed without considering cavitation is described. By applying Green's theorem and the Green function method, an integral equation for the perturbation velocity potential is obtained under the potential flow theory. Dirichlet type boundary condition is used instead of Neumann type boundary condition. The 2-D hydrofoil is approximated by line panels which have constant source strength and constant doublet strength distributions. The free surface condition is linearized and the method of images is used for satisfying this free surface condition. All the terms in fundamental solution (Green function) of perturbation potential are integrated over a line panel. Pressure distribution, lift, residual drag and free surface deformations are calculated for NACA4412, symmetric Joukowski and van de Vooren profile types of hydrofoil. The results of this method show good agreement with both experimental and numerical methods in the literature for the NACA4412 and symmetric Joukowski profile types. The lift and residual drag values of the van de Vooren profile are also presented. The effect of free surface is examined by a parametric variation of Froude number and depth of submergence.     

A potential based panel method for 2-D hydrofoils

A potential based panel method for the hydrodynamic analysis of 2-D hydrofoils moving beneath the free surface with constant speed without considering cavitation is described. By applying Green's theorem and the Green function method, an integral equation for the perturbation velocity potential is obtained under the potential flow theory. Dirichlet type boundary condition is used instead of Neumann type boundary condition. The 2-D hydrofoil is approximated by line panels which have constant source strength and constant doublet strength distributions. The free surface condition is linearized and the method of images is used for satisfying this free surface condition. All the terms in fundamental solution (Green function) of perturbation potential are integrated over a line panel. Pressure distribution, lift, residual drag and free surface deformations are calculated for NACA4412, symmetric Joukowski and van de Vooren profile types of hydrofoil. The results of this method show good agreement with both experimental and numerical methods in the literature for the NACA4412 and symmetric Joukowski profile types. The lift and residual drag values of the van de Vooren profile are also presented. The effect of free surface is examined by a parametric variation of Froude number and depth of submergence.     

Performance analysis of 3D hydrofoil under free surface

The purpose of the present paper is to develop a potential-based panel method for determining the steady potential flow about three-dimensional hydrofoil under free surface. The method uses constant-strength doublets and source density distribution over the foil body surface and thereby Dirichlet-type boundary condition is used instead of Neumann-type condition. On the undisturbed free surface source density is used to meet the free surface condition that is linearised in terms of double-body model approach and is discretised by a one-side, upstream, four-point finite difference operator. After solving the doublets on the foil and sources on the free surface, the numerical results of pressure, lift and resistance coefficients and also wave profiles can then be calculated for different Froude number and depth of submergence to demonstrate the influence of free surface and aspect ratio effects on performance of the hydrofoil.     

基于大涡模拟的波状前缘水翼空化抑制研究

为了降低空化造成的水动力性能损失,基于仿生学原理,参考座头鲸鳍肢剖面形状,将前缘波浪构型引入到水翼设计中,研究波状前缘水翼的非定常空化特性,并探究前缘参数改变对空化控制的效果和规律。选用NACA63₄-021水翼为基准模型,进行前缘参数化重构,设计出3种不同的波状水翼进行对比研究。采用大涡模拟（LES）方法对空化流场进行精细化数值模拟,针对基准水翼和不同波幅与波长参数下的波状水翼开展了空化周期、升阻力系数、压力脉动以及流向涡结构的对比分析。结果发现,波状水翼在抑制空化和降低压力脉动方面都取得了显著效果。其中,3种不同的波状水翼空化抑制率分别为15.7%、18.6%和27.9%,压力脉动幅值分别降低了55.3%、67.3%和74.6%。分析表明,波浪前缘的引入使得空化的分区效应更加凸显,空化从波谷处初生,增大波幅或减小波长都可以加强对空化的抑制效果,并可以提高升力系数以及显著降低水翼表面的压力脉动。前缘波浪构型还将诱发向下游发展的对转涡结构,不同前缘参数的波状水翼涡结构的演化是相似的,空泡发展与溃灭的整个过程对涡结构的发展也具有显著影响。     

Isoparametric Boundary Element Method Applied to the Rectangular and Delta Hydrofoils near the Free Surface

In this paper, the hydrodynamic characteristics and flow field around rectangular and delta hydrofoils, moving with a constant speed beneath the free surface are numerically studied by means of isoparametric boundary element method (IBEM). The quantities (source and dipole strengths) and the geometry of the elements are represented by a linear distribution. Two types of three-dimensional hydrofoils (rectangular and delta) are selected with NACA4412 and symmetric Joukowski sections. Some numerical results of pressure distribution, lift, wave-making drag coefficients and velocity field around the hydrofoils are presented. Also, the wave pattern due to moving hydrofoil is predicted at different operational conditions. Comparisons are made between computational results obtained through this method and those from the experimental measurements and other numerical results which reveal good agreement.     

Deep ocean wave energy conversion using a cycloidal turbine

A lift based wave energy converter, namely, a cycloidal turbine, is investigated. This type of wave energy converter consists of a shaft with one or more hydrofoils attached eccentrically at a radius. The main shaft is aligned parallel to the wave crests and submerged at a fixed depth. In the two-dimensional limit, i.e. for large spans of the hydrofoil (or an array of these), the geometry of the converter is suitable for wave termination of straight crested Airy waves. Results from two-dimensional potential flow simulations, with thin hydrofoils modeled as either a point vortex or discrete vortex panel, are presented. The operation of the cycloidal turbine both as a wave generator as well as a wave-to-shaft energy converter interacting with a linear Airy wave is demonstrated. The impact on the performance of the converter for design parameters such as device size, submergence depth, and number of hydrofoils is shown. For optimal parameter choices, simulation results demonstrate inviscid energy conversion efficiencies of more than 99% of the incoming wave energy to shaft energy. This is achieved using feedback control to synchronize the rotational rate, blade pitch angle, and phase of the cycloidal wave energy converter to the incoming wave. While complete termination of the incoming wave is shown, the remainder of the energy is lost to harmonic waves traveling in the up-wave and down-wave directions.     

典型参量对全被动式振荡水翼获能规律的影响

全被动式振荡水翼是振荡水翼式潮流能装置的主要型式之一。对全被动式振荡水翼进行了数值模拟研究,构建了二维数值模型,研究了雷诺数及升沉刚度对全被动式振荡水翼水动力性能的影响。通过流场结构及水动力性能分析研究了典型参量对全被动式振荡水翼获能性能的影响机理,确定了维持其良好水动力性能的参数范围。研究发现,水翼的水动力性能对雷诺数及升沉刚度的变化较为敏感。雷诺数增大,水翼所需升沉刚度随之增加,且水翼可以在更大的参数范围下获得较优的水动力性能。另外,水翼可以在没有升沉刚度的情况下实现周期性运动,其获能甚至优于一些有升沉刚度的情况。最优工况下,平均功率系数和能量转换效率分别为1.07和27.48%。     

涡流发生器对潮流能水轮机翼型水动力学特性影响的数值模拟

为了改善潮流能水轮机叶片表面流动分离问题,提高其升阻比,本文通过在潮流能水轮机叶片表面加装涡流发生器,来研究涡流发生器对潮流能水轮机水动力学性能的影响。本文以NACA4418翼型为研究对象,分别建立了含VGs和不含VGs的三维模型,利用CFD方法研究了VGs的高度、长度以及相邻一对VGs之间的间距等多个方面对该翼型性能的影响。结果表明:VGs可以有效地提高翼型的最大升力系数;相邻VGs间距的增加对流动分离的抑制有积极影响。此外,通过对尾迹区流线和旋涡的分析,进一步揭示了尾迹区的流场特征。     

On the formulation of a finite element method for the stiffened multi-layered airfoil/hydrofoil structure: Post buckling analysis for the wings of underwater gliders

A finite element method is developed for the stiffened multi-layered airfoil/hydrofoil structure for the large deformation and finite strain problem. The kinematics of the airfoil/hydrofoil is set up. The Consistent Orthogonal Basis Function Space is applied for the airfoil/hydrofoil structure. Given the airfoil/hydrofoil configuration and boundary conditions, the basis function space can be uniquely determined, such that the diagonal mass matrix is obtained accurately and the basis functions are very identical with the mode shape functions of the structure. In order to satisfy the displacement compatibility condition between adjacent layers of the airfoil/hydrofoil, the traction degree-of-freedom is also induced.The post buckling analysis is presented for the wing (hydrofoil) structure of the underwater glider. The water pressure is applied on the outer surface and the critical buckling pressure is calculated. The post buckling equilibrium path is also given. The results are verified with ANSYS. The present study of the buckling analysis of the airfoil/hydrofoil under water pressure is helpful to the design of underwater glider.     

Analysis of a jet-controlled high-lift hydrofoil with a flap

A jet-controlled high-lift hydrofoil with a flap is investigated using both experimental and computational methods. Experiments were carried out in a cavitation tunnel to measure forces and moment acting on the hydrofoil, and surface pressure distribution. The measured data show the feasibility of such a device for marine applications. Computational studies have also been carried out in parallel with the measurements. The computational results are analyzed in terms of global and local quantities using available experimental data. The present computational results compare well with the well-known experimental data for circulation control flows. The results for flow around a hydrofoil with a blown flap further validate the concept behind the proposed device. The results of the study demonstrate the applicability of the technology to the design of practical control surfaces.     

Analysis of loads,motions and cavity dynamics during freefall wedges vertically entering the water surface

In this paper, theoretical models are developed and numerical methods are used to analyze the loads, motions and cavity dynamics for freefall wedges with different deadrise angles vertically entering the water surface at Froude numbers: 1 ≤ Fn < 9. The time evolutions of the penetration depth, the velocity and the acceleration are analyzed and expressed explicitly. The maximum and average accelerations are predicted. The theoretical results are compared with numerical data obtained through a single-fluid BEM model with globally satisfactory agreement. The evolution of the pressures on the impact side is investigated. Before flow separation, gravity and the acceleration of the wedge have negligible influence on the pressure on the impact side for large Froude numbers or small deadrise angles; with increasing the deadrise angle or decreasing Froude number, the effects of gravity and the acceleration of the wedge tend to become more important. Global loads, with the main emphasis on the drag coefficient, are also studied. It is found that for the light wedge, the transient drag coefficient has slow variation in the first half of the collapse stage and rapid variation in the last half of the collapse stage. For the heavy wedge, the transient drag coefficients vary slowly during the whole collapse stage and can be treated as constant. The characteristics of the transient cavity during its formation are investigated. The non-dimensional pinch-off time, pinch-off depth and submergence depth at pinch-off scale roughly linearly as the Froude number.     

Hydrodynamic performance improvement on small horizontal axis current turbine blade using different tube slots configurations

Vortex generators are used extensively as a passive flow control devices to delay or remove the boundary layer separation, which affects the hydrodynamic performance of the hydrofoil. In this paper, a new approach is introduced to overcome the boundary layer separation on the NACA S1210 hydrofoil. The outcome of tube slots combination in the S1210 hydrofoil on the boundary layer separation are numerically investigated. The performance is compared with respect to the force coefficients and glide ratio. The effects of tube slot inlet positions with different diameters on S1210 hydrofoil are presented here. The results show that the smaller diameter tube slots starting near the leading edge improves the hydrodynamics performance of the hydrofoil.     

Hydrodynamics of a submerged hydrofoil advancing in waves

The hydrodynamic problem of a hydrofoil travelling at constant speed in water waves has been investigated through velocity potential theory. The boundary conditions on the free surface have been linearized, and the effects are accounted for through the Green function. The overall problem is decomposed into the steady forward speed problem and periodic wave radiation and diffraction problems. Each of these problems is solved using the boundary integral equation over the hydrofoil surface together with a vortex sheet behind the trailing edge. The body surface boundary condition is imposed on its mean position. As a result the steady potential will contribute a well-known m_j term to the body surface boundary condition on the radiation problem. The numerical difficulty in dealing with this term is effectively resolved through a difference method. The effects of the thickness on the wave radiation and diffraction are investigated. The applicability of various reciprocity relationships in this problem is discussed.     

Numerical analysis of water exit for a sphere with constant velocity using the lattice Boltzmann method

In this paper, the three-dimensional water exit of a sphere with different vertical velocities is investigated numerically using the lattice Boltzmann method (LBM). In this method, the liquid-gas two-phase flow is simplified as a single-phase free surface flow. To capture the free surface, a mass tracking algorithm is incorporated into the LBM. The gravity as a body force is introduced in the form of calculating the equilibrium distribution with an altered velocity, while the surface tension is neglected. Besides, the employed bounce-back boundary conditions are used for a moving sphere. What’s more, the Wall-Adapting Local Eddy (WALE) viscosity model is employed to capture the turbulent structures of the flow and stabilize the simulation. The accuracy of the numerical results is demonstrated through comparisons with the previous numerical and experimental results in the literature. The results show that the spike height is significantly influenced under the Froude number (Fr) below 4.12 and slightly affected under the Fr varying from 4.12 to 8.24. After the sphere exits water totally, the evolution of the free surface waterfall can be described as two phases and becomes more intense with the Froude number increasing. The non-uniform distribution of velocity results in the breaking of the free surface after the sphere completely exits the water. Moreover, the Reynolds number greatly affects the wake dynamics and hydrodynamics acting on the sphere when it moves beneath the water surface.     

The effect of appendages on the hydrodynamic characteristics of an underwater vehicle near the free surface

An underwater vehicle typically has various appendages such as sail, rudders and hydroplanes. These appendages affect the hull hydrodynamic characteristics, including the resistance components and the form of the generated wave due to the motion of the vehicle near the free surface. The effect of the appendages on the hydrodynamic characteristics of an underwater vehicle near the free surface is studied. Initially the DARPPA SUBOFF submarine without the appendages is selected and hydrodynamic characteristics, including the friction resistance, viscous pressure resistance, wave resistance and shape of the created wave on the free surface are calculated for Froude numbers in the range of 0.128–0.84 and non-dimensional submergence depths 1.3, 2.2, 3.3 & 4.4. Then, by adding the appendages and comparing these two conditions, the effect of appendages is obtained. The results of computations indicate that the appendages cause a mean increase of about 16% in the total resistance. This increment is due to viscosity of fluid and also the interaction of the main hull with the appendages. There are no significant changes in the wave pattern and wave making resistance due to the presence of appendages.     

Hull hydrodynamic optimization of autonomous underwater vehicles operating at snorkeling depth

Traditionally autonomous underwater vehicles (AUVs) have been built with a torpedo-like shape. This common shaping is hydrodynamically suboptimal for those AUVs required to operate at snorkeling condition near the free surface. In this case, the wave resistance associated to the wavy deformation of the sea surface induced by the motion of the platform is an important component of the drag. This work has investigated the optimum hull shape of an underwater vehicle moving near the free surface. Specifically a first-order Rankine panel method has been implemented to compute the wave resistance on a body of revolution moving close to the free surface. A simulated annealing algorithm was then employed to search those set of parameters defining the hull shape that minimize the wave resistance. The optimization was constrained to keep constant the total volume of the vehicle. The total drag of scaled models of the torpedo-like and resulting optimum shapes was measured in the naval tank of the University of Trieste. Measurements showed a smaller resistance of the optimized shape in the range of the considered Froude numbers.     

URANS investigation of the interaction between the free surface and a shallowly submerged underwater vehicle at steady drift

An axisymmetric underwater vehicle (UV) at a steady drift angle experiences the complex three-dimensional crossflow separation. This separation arises from the unfavorable circumferential pressure gradient developed from the windward side toward the leeward side. As is well known, the separated flow in the leeward side gives rise to the formation of a pair of vortices, which affects considerably the forces and moments acting on the UV. In this regard, the main purpose of the present study is to evaluate the role of the leeward vortical flow structure in the hydrodynamic behavior of a shallowly submerged UV at a moderate drift angle traveling beneath the free surface. Accordingly, the static drift tests are performed on the SUBOFF UV model using URANS equations coupled with a Reynolds stress turbulence model. The simulations are carried out in the commercial code STARCCM+ at a constant advance velocity based on Froude number equal to F_n = 0.512 over submergence depths and drift angles ranging from h = 1.1D to h = ∞ and from β = 0 to β = 18.11°, respectively. The validation of the numerical model is partially conducted by using the existing experimental data of the forces and moment acting on the totally submerged bare hull model. Significant interaction between the low-pressure region created by the leeward vortical flow structure and the free surface is observed. As a result of this interaction, the leeward vortical flow structure appears to be largely responsible for the behavior of the forces and moments exerted on a shallowly submerged UV at steady drift.     

Implementation of PISO algorithm for simulating unsteady cavitating flows

The present study focuses on the implementation of PISO algorithm to simulate cavitating flows. For simulation of unsteady behaviors of cavitation which have practical applications, the development of unsteady PISO algorithm based on the non-conservative approach is investigated. The effects of mixture compressibility are considered to improve accuracy of simulations. For multi-phase simulation, single-fluid Navier–Stokes equations, along with the volume fraction transport equation, are employed. The bubble dynamics model is utilized to simulate phase change. To prove capabilities of the developed PISO algorithm to simulate cavitating flows, unsteady simulation of cavitation around NACA0015 hydrofoil, a two dimensional flat plat, and a three dimensional circular disk are performed. The frequency of flow, pressure distributions, cavitating vortex shedding, and cavity characteristics are analyzed to discern results accuracy. To investigate accuracy of results, comparisons with available published experimental data are made and good agreement is achieved.     

On the design of a longitudinal motion control system of a fully-submerged hydrofoil craft based on the optimal preview servo system

The current control system of a fully-submerged hydrofoil craft based on optimal feedback control has several problems: it has no good performance of contouring waves and it cannot reduce the effect of wave disturbance. In this paper, we apply the optimal preview servo system to the longitudinal control system of a fully-submerged hydrofoil craft in order to settle the problems of current control system. The control system based on optimal preview servo system is composed of feedback controller and feedforward controller. Also in order to design the control system, the necessary prediction range of future information, the weight function in performance index, the future reference input and wave disturbance are determined. Finally, the validity of the proposed control system in this paper is confirmed by simulation. The simulation results show that the control system has good performance of contouring waves in following seas and is effective to settle the problems of current control system of a fully-submerged hydrofoil craft.     

船舶螺旋桨尾流场的数值分析

利用基于速度势的低阶面元法计算船舶螺旋桨的尾流场。采用计算较为简捷的关于扰动速度势的基本积分微分方程，并采用双曲面形状的面凶以消除面元间的缝隙。Newton-Raphson迭代过程被用来在桨叶随边满足压力Kutta条件，使桨叶面上表面的压力在随边有良好的一致性。在计算面元的影响系数时，应用了Morino导出的解析计算公式，加快了数值计算的速度。从解面元法的基本积分方程得到的偶极强度和源汇强度，直接求得尾流场的速度分布。