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

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

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

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

Numerical analysis of the influence of fixed hydrofoil installation position on seakeeping of the planing craft

This paper analyzes the hydrodynamic performance of a planing craft with a fixed hydrofoil in regular waves. Numerical simulations are carried out based on a RANS-VOF solver to study the hydrodynamic performance of the planing craft and the influence of the fixed hydrofoil on its seakeeping. To validate the numerical method, a series of hydrodynamic experiments of a bare planing craft without the hydrofoil were carried out, from which the seakeeping performance of the planing craft was recorded, the numerical method based on overset grid was compared with the experiment and verified reliable. Eight hydrofoil design cases were then studied, whereby, their seakeeping performance in regular wave conditions were predicted through the numerical method which has been verified reliable and compared with each other. Effects of hydrofoil parameters, such as angle of attack and installation height, on the seakeeping performance were investigated. Finally, the suitable installation parameters which can optimize the performance of hydrofoil and reduce the negative influence are verified. The influence of the speed on the effect of the hydrofoil and the flow field around the planing craft are also investigated.     

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.     

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.     

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.     

Calculation of hydrofoil craft take-off speed including influence of foil size, foil angle and propeller shaft angle

The sensitivity of the hydrofoil take-off speed to changes in foil size and angle has not been well documented. This paper first presents the verification of the UNO-HYF1 with published performance data. Following the satisfactory validation, the results are presented in hydrofoil take-off to show the change with speed when the foil area, foil angle and propeller shaft angle are changed. The calculations indicate that varying the hydrofoil area has the largest impact on hydrofoil take-off speed.     

Wind tunnel experiments on a fully appended laminar flow submersible for oceanographic survey

This paper reports experimental measurements of the hydrodynamics forces acting on a low drag submersible. The model was manufactured using standard, cost-effective methods without undue attention paid to the roughness or waviness of the model's surface. It was seen, using flow visualisation techniques, that the laminar boundary layer extended over 70% of the vehicle length; furthermore, no separation of the flow was seen over the remaining 30% of the body. Hydrodynamic forces were measured, and a significant reduction in drag was found compared with standard turbulent submersible shapes. The low drag performance was quite resilient to changes in incidence, showing a large drag “bucket” for angles of incidence between ± 4°. With the addition of a set of cruciform fins at the rear of the body, incorporating control surfaces, it was found that for a centre of gravity position forward of 41.3% of the body length, the vehicle would be statically stable both stick-free and stick-fixed. Finally, using a pitot tube wake rake the velocity profile of the boundary layer at the propeller location was found to be that of an attached, turbulent boundary layer closely matching a power law. The importance of full interference measurement tests to allow for the support strut effects is clearly demonstrated.     

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

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

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.     

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.     

Effects of bottom topography and density stratification on variations of western boundary currents caused by periodic change in wind stress

Variations of the western boundary currents induced by a periodic change in wind stress are studied in a two-layer model with a continental slope along the western boundary. The variation of the total transport of the western boundary current over the continental slope shows a considerable phase lag with the wind stress and a decrease in amplitude compared with for the flat bottom ocean, though the interior barotropic response is to adjust almost instantaneously to the wind stress. The total transport variation of the western boundary current is well approximated by the upper layer transport variation. That is, almost complete separation of the upper- and lower-layer flows takes place over the slope, and only the upper layer flow contributes to the change in total transport of the western boundary current. Contributions of the interior barotropic and baroclinic responses to the upper layer transport variation depend on the forcing period. With decrease in the forcing period, the barotropic response becomes relatively important for determining the upper layer transport variation although the amplitude of the variation is smaller.     

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.     

High-speed submerged and surface piercing cavitating hydrofoils, including tandem case

The iterative method that is originally developed before both for two- and three-dimensional single cavitating hydrofoils moving with a constant speed under free surface is applied to the case of high-speed (Froude number up to 6.5) and some figures are given. The method is also extended to include the surface piercing hydrofoils (vertical struts) and the case of tandem hydrofoils into the calculations. The iterative nonlinear method based on the Green's theorem allows separating the cavitating hydrofoil problem(s) and the free surface problem. These two (or three in the case of tandem hydrofoil) problems are solved separately, with the effects of one on the other being accounted for in an iterative manner. The cavitating hydrofoil surface(s) and the free surface are modeled with constant strength dipole and constant strength source panels. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream and transverse boundaries. The cavitation number is expressed in terms of Froude number and the submergence depth of the hydrofoil from the free surface. An algebraic grid on the free surface has been described to get a smooth transition between the panels along the direction of uniform inflow and to have a long distance in the downstream direction depending on the wave-length (or Froude number) while keeping the number of panels fixed. First, the method is validated in the case of surface piercing hydrofoil. Then, the effects of high Froude number and the submergence depth of the hydrofoil from free surface on the results are discussed and some figures are given for interested engineers and designers. The method is later applied to the case of tandem hydrofoils and the effects of one hydrofoil on the other are discussed.     

Two-range scaling for tip vortex cavitation inception

The suggested scaling concept for cavitation inception index takes into account boundary layer separation caused by tip vortices, as well as difference in dependencies of boundary layer thickness on Reynolds number for laminar and turbulent boundary layer. Provided asymptotic analysis gives two scaling laws with different values of power for high and low Reynolds numbers. Deduced values correspond well to known experimental data.     

Effects of bottom topography and density stratification on the formation of western boundary currents

A wind-driven, general circulation for a two-layer ocean with continental shelf-slope along the western boundary is studied numerically. Special attention is focused on the formation process of the western boundary current in the subtropical gyre. The western boundary current develops in the upper layer along the western boundary on the shelf-slope with a bottom trapped poleward flow in the lower layer. The poleward undercurrent is concentrated approximately along the contour lines of the potential vorticity,f/D, wheref is the Coriolis parameter andD the depth of the ocean. The separation of upper- and lower-layer flows on the shelf-slope represents a typical transient response. As the response approaches a steady state, the poleward undercurrent decreases in amplitude, and the motion tends to be confined to the upper layer. The flow pattern becomes similar to that found in a flat bottom ocean. A steady-state response is expected to be isostatic (no motion in the lower layer), even on the shelf-slope, as conservation of potential vorticity would suggest.The remarkable increase in transport of the western boundary current produced by the formation of an anticyclonic vortex on the shelf-slope extending throughout the hemisphere (Holland, 1973) does not occur in the wind-driven general circulation.     

Dynamic and diffusive growth of microbubbles near a two-dimensional hydrofoil

A technique for predicting the bubble growth along a two-dimensional hydrofoil with traveling bubble cavitation is presented. The method is based on the dynamic response of ambient microbubbles to the flow field and the subsequent diffusion of dissolved air into the flow field cavities. The bubble growth model is divided into three components, including the prediction of 1) the hydrofoil surface pressure distribution, 2) the ambient microbubble response to the pressure distribution, and 3) the diffusive mass flow rate. The hydrofoil velocity and pressure field is determined by two-dimensional thin airfoil theory. The microbubble response to the pressure field is given by the Rayleigh-Plesset equation with the addition of a mass diffusion term. The diffusion of dissolved gasses into the cavitation bubbles is determined by a solution to the steady-state diffusion equation under spherically symmetric convective flow. Results are given for the bubble wake of a NACA 66-006 (a = 0.8meanline) hydrofoil with traveling bubble cavitation. The effect of the relative velocity of the cavitation bubbles with respect to the surrounding water is investigated as well as the significance of the mass diffusion term in the Rayleigh-Plesset equation.     

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.     

Upwelling into the thermocline of the Pacific Ocean

The established “island rule” and the recently introduced “separation formula” are combined to yield an analytical expression for the total upwelling into the thermocline in the Pacific. The combination of the two is achieved with the use of a hybrid model containing a stratified upper layer, a thick (slowly moving) homogenous intermediate layer and an inert lower layer. Both the upper and the intermediate layers are subject to diabatic cooling and heating (which need not be specified) and there is an exchange of mass between the two active layers. An attempt is made to examine the above analytical (hybrid) model numerically. Ideally, this should be done with a complete two-and-a-half layer model (with upwelling and downwelling), but such a model is much too complex for process-oriented studies (due to the required parameterization of vertical mixing). Consequently, we focus our attention on verifying that the separation formula and the island rule are consistent with each other in a much simpler, layer-and-a-half model (without upwelling). We first verified that the new “separation formula” provides a reasonable estimate of the wind-induced transport in an island-free basin. We then compare the wind-induced transport predicted by the separation formula and the island rule in an idealized basin containing an island. We show that in these idealized situations the two methods give results that are consistent with each other and the numerics. We then turned to an application of the (hybrid) two-and-a-half layer model to the Pacific where, in contrast to the idealized layer-and-a-half models (where the two methods address the same water mass), the two methods address two different water masses. While the separation formula addresses only thermocline water (σ_θ<26.20), the island rule addresses all the water down to 27.5σ_θ (i.e., both the upper and intermediate layer). This is why the application of the two methods to the Pacific gives two different results — an application of the formula gives zero warm water transport whereas an application of the island rule gives 16 Sv. Namely, the difference between the amount predicted by the island rule (16 Sv) and the amount predicted by the separation formula (zero) enters the Pacific as intermediate water and is then somehow upwelled into the thermocline. The upwelling should take place north of the southern western boundary currents separation (40°S).     

理想流体层中悬浮隧道管体动水荷载研究-P波

针对理想流体层中的悬浮隧道管体及锚索体系,借助波动方程理论对平面P-波作用下的悬浮隧道,通过引入其管体只在竖直方面上运动和产生刚体位移而不发生变形以及只考虑海水和管体质量产生的动力效应等假定,同时考虑管体上下部海水和锚索刚度及间距的影响,导出了在理想流体层中悬浮隧道管体受到的动水荷载的计算方程组和边界条件。结合一座典型的悬浮隧道,通过参数研究,分析了上部海水厚度h,P-波频率值ω,入射角度α,锚索刚度K_c和锚索间距L对动水荷载值的影响,结果表明:上部海水的存在有利于减小动水荷载值,P波频率值与海水频率一致时将产生共振,入射角度的增大会大幅减小动水荷载值,较小的锚索刚度和较大的锚索间距也会减小动水荷载值。