A combined method for the hydrodynamic characteristics of planing crafts

The prediction of the total resistance of planing crafts at high speeds is very important. In this paper, a combined method is investigated for determining the hydrodynamic characteristics of planing crafts in the calm water. The study consists of a potential-based boundary element method (BEM) for the induced pressure resistance, the boundary layer theory for the frictional resistance and practical method for the spray resistance. The planing surface is represented by a number of elements with constant velocity potential at each element. The unknown-induced pressure is obtained by using the free surface elevation condition and the Kutta condition at the transom stern. Hydrodynamic-induced resistance and lift are determined by the calculated dynamic pressure distributions. The boundary layer analysis method is based on calculations of the momentum integral equation applied to obtain the frictional resistance. A particular practical approach is introduced to present the region of the upwash geometry for the spray. A numerical program has been developed for the present research and applied to the hull form of the craft. Four different hull forms of Series 62 model 4666 planing craft are presented. It is shown that the present combined method is efficient and the results are in good agreement with the experimental measurements over a wide range of volumetric Froude numbers.     

低雷诺数下圆柱强迫振动特性光滑粒子流体动力学模拟

基于光滑粒子流体动力学(SPH)方法,开发了能够准确描述水流作用下圆柱强迫振动特性的数学模型。通过引入适合于无网格粒子法的开边界算法,来模拟出入流边界条件,建立了具有造流功能的SPH数值水槽。圆柱及计算域的上下边界均采用修正的动力边界条件进行模拟。借助于粒子位移矫正和压力修正算法,避免了圆柱周围流体粒子压力大幅震荡以及结构下游区域出现空腔等非物理性现象。使用典型的圆柱绕流数据,验证了所建SPH模型的计算性能,研究了固定圆柱在低雷诺数情况下的尾涡脱落模式和升阻力变化规律。明确了低雷诺数下强迫振动圆柱在频率锁定以及非锁定区间内的升力变化规律,量化了升力与外界激励频率之间的关系。     

The spray volume shed by an uncambered planing hull in steady planing

Earlier papers (Payne, 1981 a,b,c) have developed what might be called a virtual mass theory which in principle permits the forces on any planing hull form to be calculated. In the present paper, this methodology is extended to calculate the thickness and momentum of the jet or spray sheet thrown off by the planing surface. For a two-dimensional flat planing plate—the only case where comparison is possible—the theory gives essentially the same result as that of Pierson and Leshnover (1948). For a three-dimensional flat plate and prismatic hulls, the results seem physically reasonable.For the small trim angles associated with efficient planing, on a weightless inviscid fluid the total pressure drag of any hull can be reduced to close to zero by deflecting the jet rearwards and parallel to the undisturbed surface, the residual resistance being due to the cross-flow force which varies as (trim angle)².     

Numerical investigation of wave elevation and bottom pressure generated by a planing hull in finite-depth water

A numerical investigation of the bottom pressure and wave elevation generated by a planing hull in finite-depth water is presented. While the existing literature addresses the free-surface deformation and pressure field at the seafloor independently, this work proposes a direct comparison between the two hydrodynamic quantities. The dependence of the pressure disturbances at the ocean floor from the waves generated at the free-surface by a planing hull is studied for several values of both the depth and hull Froude numbers. The methodology employed is Smoothed Particle Hydrodynamics (SPH), a numerical technique based on the discretization of the continuum fields of hydrodynamics through mesh-less particles. The SPH code herein chosen is initially validated against experimental data for transom-stern flow. Subsequently, numerical simulations are presented for a planing hull in high-speed regimes. The results show a direct correlation between surface wave dynamics and hydrodynamic pressure disturbances at the seafloor as the value of the Froude number is varied. This is assessed by studying the inverse dependence of the low-pressure wake angle with the Froude number and by comparison of SPH results with similar works in the cited literature.     

聚焦波作用下平面网衣结构的水动力特性研究

由于沿海区域的限制以及愈加严重的环境污染,渔业养殖正从近海走向深远海。深远海海域的海况更加恶劣,给养殖装备的设计与性能评估带来新的挑战。为解决该问题,对极端波浪与养殖装备网衣结构的相互作用开展研究。基于waves2Foam建立数值波浪水池,极端波浪模拟采用基于NewWave理论的聚焦波模型,网衣结构模拟采用多孔介质模型,并通过与Morison模型计算的网衣受力等效分析,获得多孔介质模拟网衣结构阻力系数的直接估计方法。然后将多孔介质模型嵌入waves2Foam中,开展聚焦波与网衣结构相互作用的数值模拟,同时开展水槽试验,验证数值模拟的准确性。基于数值模拟结果,系统地分析了不同网衣密实度及不同波浪参数下网衣结构的升阻力特性以及网衣结构对波浪场的扰动规律。研究表明：聚焦波波峰幅值和网衣密实度对网衣结构的升阻力影响较大,且升力峰值出现在阻力为0的时刻;网衣结构对聚焦波的时空演化特性有影响,改变了聚焦波波形。     

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.     

Interaction between submarine landslides and suspended pipelines with a streamlined contour

Recently, the security and stability of submarine pipelines have attracted much attention in ocean engineering. In this paper, pipelines with a streamlined contour (wedge, airfoil, double-ellipse, and arc-angle hexagon) are designed in hopes of defending against the impact of submarine landslides, and the computational fluid dynamics (CFD) approach is used to investigate the interaction between submarine landslides and streamlined pipelines. The results show that the peak interactional force is more representative of the hazard level of pipelines imposed by submarine landslides. It is also found that the streamlined pipelines possess a significant advantage in reducing the drag force and lift force of landslide–pipeline interaction with a maximum lessening percentage of 66.32 and 40.17%, compared with a conventional circular pipeline. In addition, the influence of applying streamlined pipelines to engineering is briefly discussed, and the empirical equation for estimating the drag force and lift force of streamlined pipelines induced by landslides is recommended based on the numerical test results.     

The interceptor hydrodynamic analysis for controlling the porpoising instability in high speed crafts

A well-known instability of the high-speed planing crafts is the porpoising instability. This instability involves periodic, coupled heave/pitch oscillations possibly experienced in a planing vessel at high speeds. The porpoising can be controlled by using external devices. Interceptors are vertical blades installed symmetrically at the aft of the craft and have been introduced as a trim control appendage. Here, based on numerical methods and Savitsky porpoising theory, the effects of hydrodynamic interceptors on the porpoising control are investigated. Using computational fluid dynamics, the pressure distribution created by interceptor and its effects on porpoising are computed and then discussed. To model the flow around the vessel model, the Reynolds Average Navier Stokes (RANS) equations are applied. The work deals with craft with and without an interceptor at different heights. A dynamic grid mode involving two degrees of freedom is used. The results show that the interceptor causes an intense pressure at the stern bottom. It also decreases the trim and resistance of the vessel and increases the lift force coefficient which directly affects the porpoising instabilities. Based on the results, the interceptor can completely control the porpoising phenomenon.     

Computation of Viscous Uniform and Shear Flow over A Circular Cylinder by A Finite Element Method

The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from lO0 to lO00. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem.     

An automated computational method for planing hull form definition in concept design

This paper describes the development of a computer-based method for producing chined planing boat hull forms adequate to be applied in concept design. The method is based on a principle where the designer specifies a small set of critical parameters he/she wishes obtain or keep preserved and generates a complete hull form, without the traditional skilled recourse of giving stations point by point. From this set of parameters a detailed and faired drawing with offsets is generated very quickly. The method allows, in its execution mode, the flexibility to modify, adjust and enlarge the input set of parameters. The method was created to allow both (1) automated hull form definitions when integrated to an existing computer system and (2) quick but detailed preliminary calculations of stability, lift and drag, volumes and internal space allocations, sea-keeping estimates, etc., all with very reasonable precision. As application examples some planing boat hull forms are generated. Some are typical and others less usual. The later ones are defined to show the method's limits, in order to validate it.     

Large eddy simulations of a circular cylinder at Reynolds numbers surrounding the drag crisis

Large eddy simulations of the flow around a circular cylinder at high Reynolds numbers are reported. Five Reynolds numbers were chosen, such that the drag crisis was captured. A total of 18 cases were computed to investigate the effect of gridding strategy, turbulence modelling, numerical schemes and domain width on the results. It was found that unstructured grids provide better resolution of key flow features, when a ‘reasonable’ grid size is to be maintained.When using coarse grids for large eddy simulation, the effect of turbulence models and numerical schemes becomes more pronounced. The dynamic mixed Smagorinsky model was found to be superior to the Smagorinsky model, since the model coefficient is allowed to dynamically adjust based on the local flow and grid size. A blended upwind-central convection scheme was also found to provide the best accuracy, since a fully central scheme exhibits artificial wiggles, due to dispersion errors, which pollute the solution.Mean drag, fluctuating lift Strouhal number and base pressure are compared to experiments and empirical estimates for Reynolds numbers ranging from 6.31 × 10⁴ to 5.06 × 10⁵. In terms of the drag coefficient, the drag crisis is well captured by the present simulations, although the other integral quantities (rms lift and Strouhal number) show larger discrepancies. For the lowest Reynolds number, the drag is seen to be more sensitive to the domain width than the spanwise grid spacing, while at the higher Reynolds numbers the grid resolution plays a more important role, due to the larger extent of the turbulent boundary layer.     

两层粘性流体中圆柱体受迫振荡数值模拟

研究两层粘性流体中无限长水平圆柱体的受迫振荡问题。在湍流模式下,采用VOF方法追踪两层流体的内界面,基于动网格技术模拟圆柱体的运动边界,对均匀流中横向振荡圆柱体的绕流场进行了数值模拟。计算受迫振荡圆柱体的升力系数、阻力系数随时间的演化曲线和圆柱体的尾涡分布,以及圆柱体的受迫振荡激发两层流体内界面的扰动,并与均匀流体的情况进行了比较分析。研究表明,流体的两层分层效应对受迫振荡圆柱体的升阻力系数和尾涡分布特性都有显著影响,在水下输油气管道涡激振动特性的工程评估中,应考虑流体的密度分层效应。     

A numerical modeling of hydrodynamic characteristics of various planing hull forms

A numerical algorithm based on the boundary element method (BEM) is presented for predicting the hydrodynamic characteristics of the various planing hull forms. The boundary integral equation is derived using Green's theorem on the wetted body surface and the free surface. The ventilation function at the transom is estimated with Doctor's empirical formula. This function is defined as the transom zone free surface boundary condition. The combined boundary integral equation and modified free surface boundary condition are simultaneously solved to determine the dipole on the wetted hull surface and the source on the free surface. The method is applied to investigate three examples of planing hulls, which include flat-plates, as well as wedge-shaped and variable deadrise planing hulls. Their hydrodynamic characteristics are calculated for different speeds. Computational results are presented and compared with existing theories and experiments. On the whole, the agreement between the present method and the selected experimental and numerical data is satisfactory.     

动波浪壁圆柱绕流三维数值模拟

利用计算流体力学软件Fluent开展了三维动波浪壁圆柱绕流的数值模拟,建立了三维运动波浪壁圆柱模型,通过C语言自编程序实现波浪壁面的运动控制,并保证壁面变形时网格的高质量。在来流速度u=0.125 m/s、雷诺数Re=12 500的情况下,开展了动波浪壁波动速度w=0、0.062 5、0.125、0.187 5 m/s四个工况的计算分析,并比较了不同波动速度对流场结构、升力、阻力特性的影响。结果表明:动波浪壁圆柱能有效抑制流动的分离,消除交替脱落的尾涡,从而消除周期振荡的升力;在消除卡门涡街的同时,圆柱后驻点处的涡量值随波动速度增加而增加,其原因在于波形移动加大了壁面流体的速度,从而减小了圆柱前后的压力差,减小了阻力;随着波动速度的增大,平均阻力系数呈明显下降趋势,当波动速度为来流速度的1.5倍时,平均阻力系数相对于光滑圆柱下降了53.76%。     

Studies on wave-subsea pipeline interaction

Estimation of wave induced loadings on subsea pipelines involves the evaluation of the hydrodynamic coefficients of drag, inertia, vertical inertia and lift. These coefficients vary with the methodology adopted and the order of wave theory used to compute the water particle kinematics. This paper presents a new and simple method to determine the hydrodynamic coefficients for subsea pipelines in the drag and inertial flow regime. The results are reported for three relative clearances of a pipeline model from the simulated ocean floor.     

An accurate four-quadrant nonlinear dynamical model for marinethrusters: theory and experimental validation

This paper reports two specific improvements in the finite-dimensional nonlinear dynamical modeling of marine thrusters. Previously reported four-quadrant models have employed thin airfoil theory considering only axial fluid flow and using sinusoidal lift/drag curves. First, we present a thruster model incorporating the effects of rotational fluid velocity and inertia on thruster response. Second, we report a novel method for experimentally determining nonsinusoidal lift/drag curves. The model parameters are identified using experimental thruster data (force, torque, and fluid velocity). The models are evaluated by comparing experimental performance data with numerical model simulations. The data indicates that thruster models incorporating both reported enhancements provide superior accuracy in both transient and steady-state responses     

On the shape of the wake behind a planing boat

Recent developments in slender planing theory permit rather precise estimates of forces and moments to be made. So using a methodology first suggested by Epshtein for flat plate wakes, the modern theory is employed to find the longitudinal shape of a flat plate's wake. The result is in better agreement with experiment than previous formulations. The same approach is then applied (for the first time) to a prismatic wake, and the agreement with experiment is found to be excellent. Finally, the equations developed are used to deduce the wave energy developed by the passage of the hull, and from this, its wave drag. Although there are no experimental data available for comparison, the ratio (wave drag)/(total pressure drag) follows the same trend with Froude number as earlier two-dimensional solutions by Sedov and Squire, and tends to zero as Froude number tends to infinity.     

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.     

The differences between a wing and a planing plate in two-dimensional flow

The normal force coefficient on a flat planing surface having arbitrary heave and pitch motion in two-dimensional flow is compared with the lift coefficient of a thin wing in an infinite fluid. Despite the totally different derivations, they are found to be identical (at large Froude numbers and low trim angles and allowing for the wing's interaction with twice as much fluid) at low reduced frequencies. For higher frequency motions, the wing's angle of attack induced lift and its pitch and heave damping are less than those of a planing surface, but the acceleration terms remain identical. The differences at the higher reduced frequencies are due to the fact that, in invisad irrotational flow, the planning plate cannot leave a vortex wake, whereas a wing does.It seems to follow that the “virtual mass” planing hull analysis can be applied to “quasi-static” problems involving wings and bodies in an infinite fluid without the slenderness restriction originally imposed by Jones (1946). Certainly, it is remarkable that the so called “quasi-steady” forces on a two-dimensional wing can be obtained in a few lines of elementary analysis. On the other hand, the method fails entirely when used to compute the pitching moment on a two-dimensional plate, even though it has been found to give good results for the three-dimensional case (Payne, 1981c).This work is offered as a very incomplete study of an intriguing relationship between two very different bodies of analysis. Much more work will need to be done before the relationship between the two approaches will be fully understood.     

Experimental and numerical investigations on hydrodynamic and aerodynamic characteristics of the tunnel of planing trimaran

The planing trimaran possesses distinctive hybrid hydrodynamic and aerodynamic performance due to the presence of tunnel. The research described in this paper was carried out based on the observation of wave characteristics of a planing trimaran model in towing tests, in which the resistance drops as soon as the wave surface separates from tunnel roof. In order to gain a deeper insight into the relationship between wave flow and forces in tunnel region, a comprehensive series of viscous CFD simulations considering free-surface and 2-DOF motion of the hull (heave and pitch) have been performed for the tested model at the volume based Froude numbers ranging from 3.16 to 5.87. The calculated results were validated by comparison with experimental data and showed good agreement. Numerical results of wave contours, longitudinal wave cuts and lifting force distributions at the calculated speeds were presented for the analysis of ventilation process in tunnel region and the corresponding variation of tunnel forces. It is found that, for the speeds higher than Froude number of 4.52, the aerodynamic forces provide major tunnel lift and mainly act on the straight section of the tunnel. And, therefore, numerical simulations of two modified models have also been performed for the analysis of influence of straight section length on the hydrodynamic and aerodynamic performance of planing trimaran.