船体周围Stokes波浪流耦合VOF模型的数值模拟（英文）

The surface wave generated by flow around a ship hull moving near free surface of water is simulated numerically in this study. The three-dimensional implicit finite volume method(FVM) is applied to solve Reynolds averaged Navier-Stokes(RANS) equation. The realizable k-ε turbulence model has been implemented to capture turbulent flow around the ship hull in the free surface zone. The volume of fluid(VOF) method coupled with the Stokes wave theory has been used to determine the free surface effect of water. By using is a six degrees of freedom model, the ship hull's movement is numerically solved with the Stokes wave together. Under the action of Stokes waves on the sea, the interface between the air and water waves at the same regular pattern and so does the pressure and the vertical velocity. The ship hull moves in the same way as the wave. The amplitude of the ship hull's heave is less than the wave height because of the viscosity damping. This method could provide an important reference for the study of ships' movement, wave and hydrodynamics.     

Numerical Simulation of the Stokes Wave for the Flow around a Ship Hull Coupled with the VOF Model

The surface wave generated by flow around a ship hull moving near free surface of water is simulated numerically in this study. The three-dimensional implicit finite volume method(FVM) is applied to solve Reynolds averaged Navier-Stokes(RANS) equation. The realizable k-ε turbulence model has been implemented to capture turbulent flow around the ship hull in the free surface zone. The volume of fluid(VOF) method coupled with the Stokes wave theory has been used to determine the free surface effect of water. By using is a six degrees of freedom model, the ship hull’s movement is numerically solved with the Stokes wave together. Under the action of Stokes waves on the sea, the interface between the air and water waves at the same regular pattern and so does the pressure and the vertical velocity. The ship hull moves in the same way as the wave. The amplitude of the ship hull’s heave is less than the wave height because of the viscosity damping. This method could provide an important reference for the study of ships’ movement, wave and hydrodynamics.     

Improving longitudinal motion prediction of hybrid monohulls with the viscous effect

A new method improves prediction of the motion of a hybrid monohull in regular waves. Stem section hydrodynamic coefficients of a hybrid monohull with harmonic oscillation were computed using the Reynolds Averaged Navier-Stokes Equations (RANSE). The governing equations were solved using the finite volume method. The VOF method was used for free surface treatment, and RNGК-ε turbulence model was employed in viscous flow calculation. The whole computational domain was divided into many blocks each with structured grids, and the dynamic process was treated with moving grids. Using a 2-D strip method and 2.5D theory with the correction hydrodynamic coefficients allows consideration of the viscous effect when predicting longitudinal motion of a hybrid monohull in regular waves. The method is effective at predicting motion of a hybrid monohull, showing that the viscous effect on a semi-submerged body cannot be ignored.     

Interaction of two-dimensional impulsively started airfoils

Continuous vortieity panels were used to model general unsteady inviscid, incompressible, two-dimensional flows. The geometry of the airfoil was approximated by series of short straight segments having endpoints that lie on the actual surface. A piecewise linear, continuous distribution of vorticity over the airfoil surface was used to generate disturbance flow. The no-penetration condition was imposed at the midpoint of each segment and at discrete times. The wake was simulated by a system of point vortices, which moved at local fluid velocity. At each time step, a new wake panel with uniform vortieity distribution was attached to the trailing edge, and the condition of constant circulation around the airfoil and wake was imposed. A new expression for Kutta condition was developed to study the interference effect between two impulsively started NACA0012 airfoils. The tandem arrangement was found to he the most effective to enhance the lift of the rear airfoil. The interference effect between tidal turbine blades was shown clearly.     

粘流方法的螺旋桨尾流场数值分析(英文)

The computational fluid dynamics(CFD) method is used to numerically simulate a propeller wake flow field in open water.A sub-domain hybrid mesh method was adopted in this paper.The computation domain was separated into two sub-domains,in which tetrahedral elements were used in the inner domain to match the complicated geometry of the propeller,while hexahedral elements were used in the outer domain.The mesh was locally refined on the propeller surface and near the wake flow field,and a size function was used to control the growth rate of the grid.Sections at different axial location were used to study the spatial evolution of the propeller wake in the region ranging from the disc to one propeller diameter(D) downstream.The numerical results show that the axial velocity fluctuates along the wake flow;radial velocity,which is closely related to vortices,attenuates strongly.The trailing vortices interact with the tip vortex at the blades’ trailing edge and then separate.The strength of the vortex shrinks rapidly,and the radius decreases 20% at one diameter downstream.     

海岸桩结构布置对波浪参数影响的数值模拟研究（英文）

In this study, the passage of waves through pile groups with different arrangements is investigated using a three-dimensional(3D)numerical model. For the simulations, waves of three different heights of 36, 58, and 81 mm, a fixed period of 0.88s, and a fixed wave length of 1.128 m were used. To simulate the waves and flow pattern through the piles, Reynolds-averaged Navier–Stokes(RANS) equations of fluid motion were solved based on the finite volume method(FVM). Piles were defined as obstacles in the rectangular domain using the fractional area/volume obstacle representation(FAVOR) method. The volume-of-fluid(VOF) and re-normalization group(RNG) methods were used to simulate the free surface and turbulence phenomenon, respectively. By performing different numerical simulations, the effect of coastal pile arrangements on wave pattern was studied and was compared with existing experimental data, and an acceptable agreement was achieved.     

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.     

Numerical solutions for a two-dimensional airfoil undergoing unsteady motion

Continuous vorticity panels are used to model general unsteady inviscid, incompressible, and two-dimensional flows. The geometry of the airfoil is approximated by series of short straight segments having endpoints that lie on the actual surface. A piecewise linear, continuous distribution of vorticity over the airfoil surface is used to generate disturbance flow. The no-penetration condition is imposed at the midpoint of each segment and at discrete times. The wake is simulated by a system of point vortices, which move at local fluid velocity. At each time step, a new wake panel with uniform vorticity distribution is attached to the trailing edge, and the condition of constant circulation around the airfoil and wake is imposed. A new expression for Kutta condition is developed to study (i) the effect of thickness on the lift build-up of an impulsively started airfoil, (ii) the effects of reduced frequency and heave amplitude on the thrust production of flapping airfoils, and (iii) the vortex-airfoil interaction. This work presents some hydrodynamic results for tidalstream turbine.     

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

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.     

Self-propulsion computations using a speed controller and a discretized propeller with dynamic overset grids

A method that can be used to perform self-propulsion computations of surface ships is presented. The propeller is gridded as an overset object with a rotational velocity that is imposed by a speed controller, which finds the self-propulsion point when the ship reaches the target Froude number in a single transient computation. Dynamic overset grids are used to allow different dynamic groups to move independently, including the hull and appendages, the propeller, and the background (where the far-field boundary conditions are imposed). Predicted integral quantities include propeller rotational speed, propeller forces, and ship’s attitude, along with the complete flow field. The fluid flow is solved by employing a single-phase level set approach to model the free surface, along with a blended k−ω/k−ɛ based DES model for turbulence. Three ship hulls are evaluated: the single-propeller KVLCC1 tanker appended with a rudder, the twin propeller fully appended surface combatant model DTMB 5613, and the KCS container ship without a rudder, and the results are compared with experimental data obtained at the model scale. In the case of KCS, a more complete comparison with propulsion data is performed. It is shown that direct computation of self-propelled ships is feasible, and though very resource intensive, it provides a tool for obtaining vast flow detail.     

船桨干扰定常空化性能数值模拟

文章采用CFD技术计算了船桨干扰定常空化性能。首先计算了不同空化数下的NACA66(MOD)水翼和DTRC4381螺旋桨的定常空化性能,同时预报了KCS船的伴流场、阻力性能和船体表面压力分布。为了验证计算船后螺旋桨空化流的可行性,文中对KCS船和KP505螺旋桨进行了整体考虑,并计算了不同空化数下的船后螺旋桨的空化覆盖面积。计算表明,文中采用的计算方法合理,计算结果与试验结果吻合。     

基于滑移网格的带桨水面船自航性能预报研究

应用FINE/Marine软件对KCS船、KP505桨以及考虑自由液面的船桨组合体进行数值计算,并计算其自航性能。利用滑移网格技术和随体网格来实现船桨之间的相互耦合。考虑到原有的自航性能数据处理方法并不代表实际情况,文中借鉴强制自航法的概念提出了一种新的船舶自航点求解方法预报船舶自航性能,并与模型试验结果进行比较,吻合良好,其中推力减额系数、伴流分数以及船舶推进效率的计算误差分别为0.5%、2.18%、6.76%。本文研究为预报船舶自航性能提供了一种新的研究手段。     

重叠网格在船舶CFD中的应用研究

文章采用RANS方法和重叠网格计算了带自由液面的船舶绕流问题。计算网格采用重叠网格的型式,自由液面的模拟采用单相Level-Set方法,Reynolds应力采用k-ω模型,采用体单元有限差分方法和PISO算法求解RANS方程。文中简要描述了重叠网格和单相Level Set自由液面模拟方法的数学模型及求解。通过对S60单体船型约束模和自由模型采用重叠网格的数值求解及与试验结果的比较表明该重叠可较好地模拟带自由液面船舶自由态绕流问题。另采用重叠网格对一简单的双体船和三体船进行了数值求解,计算结果也表明重叠网格和单相Level-Set方法可较好地模拟带自由液面的船舶绕流问题。     

基于RANS法的船舶阻力及自由运动预报

文章介绍了基于RANS法求解船舶自由运动水动力性能的方法。通过对流域设计、网格划分方法、运动求解法及湍流模型的选择等进行完整的讨论与分析,提出了一套求解船舶自由运动的RANS方法。以低速船(KCS型船模)及高速船(Fridsma滑行艇)为计算对象,计算结果验证了RANS法在预报船舶自由航行时水动力性能的实用性。     

KCS迎浪中耐波性的CFD数值模拟

采用计算流体动力学(CFD)软件STAR-CCM,对国际集装箱标模KCS在迎浪规则波中的运动响应和波浪增阻等主要耐波性能进行数值计算。采用动态重叠网格处理船体运动,同时利用VOF方法模拟自由面波形演化情况。计算结果显示:动态重叠网格具有较强的适应性,能较好地处理船体的大幅运动,计算得到的船体运动幅值能与试验结果相一致,静水和波浪中的阻力计算结果与试验结果能较好地吻合;从自由面波形方面看,KCS船体兴波与来波存在较强的相互干扰,当波峰撞击船首时,波面起伏较大,数值计算能较好地模拟这些流动现象。     

一种基于体积力螺旋桨模型的自航计算方法

文章以集装箱船模型KCS(KRISO Container Ship)为研究对象,基于通用CFD软件FLUENT 12.0.16,采用螺旋桨体积力模型,实现了模型尺度下实船自航点的全粘带自由面计算。该文的计算包括静水拖航计算、自航计算,以及扭矩对自航计算结果的影响分析。通过与试验值的分析比较,验证了该方法的有效性。由于采用了完全结构化的六面体网格,带对称面的计算域单元数仅为33万。研究结果还表明:考虑扭矩的全流场计算对于改善流动细节,获得更精确的计算结果是有利的,但计算量有所增加。     

用于船舶伴流场预报的一种非几何相似的模型（英文）

The scale effect leads to large discrepancies between the wake fields of model-scale and actual ships, and causes differences in cavitation performance and exciting forces tests in predicting the performance of actual ships. Therefore, when test data from ship models are directly applied to predict the performance of actual ships, test results must be subjected to empirical corrections. This study proposes a method for the reverse design of the hull model. Compared to a geometrically similar hull model, the wake field generated by the modified model is closer to that of an actual ship. A non- geometrically similar model of a Korean Research Institute of Ship and Ocean Engineering(KRISO)'s container ship(KCS) was designed. Numerical simulations were performed using this model, and its results were compared with full-scale calculation results. The deformation method of getting the wake field of full-scale ships by the non-geometrically similar model is applied to the KCS successfully.     

A non-geometrically similar model for predicting the wake field of full-scale ships

The scale effect leads to large discrepancies between the wake fields of model-scale and actual ships, and causes differences in cavitation performance and exciting forces tests in predicting the performance of actual ships. Therefore, when test data from ship models are directly applied to predict the performance of actual ships, test results must be subjected to empirical corrections. This study proposes a method for the reverse design of the hull model. Compared to a geometrically similar hull model, the wake field generated by the modified model is closer to that of an actual ship. A non- geometrically similar model of a Korean Research Institute of Ship and Ocean Engineering (KRISO)’s container ship (KCS) was designed. Numerical simulations were performed using this model, and its results were compared with full-scale calculation results. The deformation method of getting the wake field of full-scale ships by the non-geometrically similar model is applied to the KCS successfully.     

绕船体粘性自由面流动的数值计算

考虑自由面的船舶粘性流动计算由于其复杂性和对计算机硬伯的较高要求，只是到了最近几年才得到较大发展，而在国内尚属起步阶段。本文使用VOF方法计算了绕Wigley船型的粘性自由面流动，流模式选择了SSTκ-ω模型。计算获得了船体上波形以及船波等值线图，同时考察了三个不同横截面上的流动状况。计算得到的阻力还和发表的试验结果作了对比，结果令人鼓舞。