基于大涡模拟的ROV导管推进器梢隙流动数值研究

导管推进器是一种普遍应用于无人遥控潜水器(ROV)等潜器中的特种推进器。在桨叶与导管之间的梢隙中存在非常复杂的流动,本研究基于大涡模拟(LES)对导管推进器的梢隙流动进行了数值模拟分析。通过对时间步长的收敛性研究,建立两套基于不同网格类型的计算模型。将计算与试验结果进行对比,比较两种不同类型网格模拟结果的差异,发现切割体网格能够更好地捕捉到泄涡的细节,并结合梢隙流场的原理分析了泄涡发展的过程,梢隙涡的驱动力是吸力面、压力面之间的压差。此外,随着进速系数增大,梢隙周向的涡管轴向分布范围减小,主泄涡发生位置延后,泄出涡的长度和数量都有所减少。     

升力面方法预报全方向推进器的非定常水动力性能

本文研究了全方向推进器非定常水动力性能的升力面预报方法,基于势流理论和格林定理建立了全方向推进器非定常水动力性能计算的数学模型,然后利用非定常涡格法对全方向推进器的非定常水动力性能进行了数值预报。本文的计算结果和日本水池模型试验结果、升力线方法计算结果对比表明,本文的计算值与试验值吻合良好。     

二维圆柱涡激振动数值模拟中湍流模式适用性的探讨

采用RANS结合4种湍流模型对低质量比单自由度涡激振动进行数值模拟,对比分析其对该问题的适用性。用四阶Runge-Kutta法离散运动方程,基于动网格技术处理圆柱振荡引起的网格运动,并对壁面条件的处理进行了细致分析。根据Williamson水槽实验,从振幅比、频率比响应、水动力系数及相位突变、尾涡模式等方面对比分析了4种模型的性能和表现,结果表明Standardκ-ω模型与实验差异较大;目前本类研究中运用较少的Realizableκ-ε模型也是可以适用于涡激振动计算的;κ-ω系列模型得到的最大振幅(0.55D)小于κ-ε系列模型的结果(0.87D);SSTκ-ω模型以及2种κ-ε模型都反映出锁定区振动频率fex与泄涡频率fst分离,其中SSTκ-ω模型较为接近实验结果;尾涡模式上,SSTκ-ω模型在各个分支与实验结果一致;总之,各种模型针对不同物理现象各有优势和缺陷。     

潮流能水轮机尾流场及涡特性DDES模拟

为研究潮流能水轮机尾流场流动特性及涡结构组成,基于DDES(Delayed Detached Eddy Simulation)模型对不同流速和转速的4种工况下水轮机尾流场进行数值模拟,并进一步探究其尾流场空间涡结构的变化特性。结果表明,水轮机的数值模拟结果与试验结果能够较好吻合。对比不同工况下的尾流场模拟结果可知:水轮机尾流区域流动复杂,延迟分离涡模拟方法能有效模拟水轮机旋转过程中产生的叶尖涡、叶尖脱落涡、轮毂涡等不同涡结构,并能完整观察到叶尖涡的产生、脱落、失稳、破碎过程。转速一定时,流速越大,叶尖脱落涡、轮毂涡的发展距离越远;流速一定时,转速越大,涡的发展距离越短。本文数值模拟计算结果可为实际海况中潮流能阵列水轮机的布局提供可靠依据。     

平板近壁面双孔通气气液两相流场的数值模拟研究

基于自编程序发展了一种平板近壁面通气气液两相湍流流动数值模拟方法,并对近壁面双孔通气气液两相流场进行了数值模拟分析。研究结果表明：建立的数值计算方法摆脱了对传统 RANS 湍流模型的依赖, 可以较好地捕捉气液界面的失稳现象,模拟结果与实验结果吻合较好。平板近壁面双孔通气气泡的掺混融合过程可以分为 3 个流动阶段：独立稳定发展阶段、初步掺混融合阶段以及充分融合阶段。流场中气泡的发展演化与旋涡结构相关,剪切层涡和反向旋转涡对在其掺混融合过程中起着主要作用。     

三维机翼分离流的计算方法

本文提出了一种用面元法计算绕三维机翼分离流的方法。将分离面与物面之间的流体考虑为不流动的死水区,应用面元法中的总速度势法（ＴｏｔａｌＰｏｔｅｎｔｉａｌＭｅｔｈｏｄ）求解布置于物面及分离面上的偶极子分布及周围流场和流体力。考虑到死水区内的流速为零,因此压力一定,自由分离面上的流速与涡强不变。分离面形状通过迭代方法确定。每次迭代中,分离点位置用二维边界层理论计算得出。计算结果与实验值进行了比较,吻合程度良好。     

改进的移动粒子半隐式法模拟楔形体入水砰击

移动粒子半隐式法(Moving-Particle Semi-Implicit Method,MPS)是一种新的基于拉格朗日(Lagrange)理念的无网格方法,适用于模拟自由液面的大变形和水流的喷射现象。用基于大涡模拟的改进MPS法首先模拟了矩形体的入水砰击,砰击压力的计算结果证明了这种方法的正确性,然后模拟了楔形体的匀速入水砰击,并与实验结果进行了对比,验证了大涡模拟改进MPS法在砰击问题中的适用性。     

带前缘结节叶片导管桨尾流不稳定性分析

导管桨的尾流不稳定性在其性能评价中非常重要,不但是其能否提供稳定推力的保证,而且也与螺旋桨的尾流噪声直接相关。为了改善导管桨的尾流,提高尾流稳定性,并优化导管桨的流场脉动,根据座头鲸鳍肢前缘结节的仿生原理,对导管桨叶片的导边进行改进,提出了两种仿生桨型,采用IDDES湍流模型对低进速系数下常规导管桨和仿生叶片导管桨进行数值模拟,探究叶片构型对导管桨性能和尾流不稳定性的影响。计算结果表明,前缘结节可以有效降低叶片受力波动的幅值和叶片所受合力的主频域峰值,具有较大结节的叶片对导管桨尾流有明显的优化作用,在尾流远场中扩大了流动稳定区,延后了尾流处涡破碎的发生,改善了能量谱密度的频域分布。进一步,大前缘结节叶片导管桨应用在低速工况下时,可以大量减少尾流泄涡区域的二次涡产生,这是由于前缘结节提升了相邻涡互感的强度,使得尾流更加稳定,而小结节叶片仿生桨型对导管桨尾流则无明显优化作用。研究方法和成果可为螺旋桨尤其是导管桨尾流不稳定性研究提供参考,不仅验证了前缘结节在导管桨叶片应用的合理性,而且揭示了其优化尾流稳定性的机理。     

Spar风电平台涡激运动特性研究

本文研究了无侧板单立柱Spar风电平台在不同约化速度(Ur)下的涡激运动特征.基于计算流体动力学进行数值模拟,采用剪应力运输(Shear Stress Transport,SST)k-ω湍流模型模拟漩涡脱落,通过自编动力积分程序和重叠网格实现涡激运动的实时流固耦合,对Ur为2～14时的平台涡激运动响应进行计算,分析涡激...     

中尺度暖涡对热带气旋强度变化的影响及作用机制

基于两组理想化数值试验，对比研究了分布于热带气旋不同位置处的海洋中尺度暖涡所引发的热带气旋强度变化的时空特征。研究发现，热带气旋中心附近的暖涡对热带气旋强度有增强作用，而位于热带气旋外围的暖涡则会抑制热带气旋的发展。本研究将暖涡增强（减弱）热带气旋强度的区域称为内（外）区。随着时间的推移，内（外）区暖涡对热带气旋强度的增强（减弱）幅度逐渐减小（增大），区域范围同步减小（增大）。内区暖涡增强了热带气旋的次级环流和结构对称性、增加了海气界面热通量，同时减弱了外围螺旋雨带，进而导致热带气旋强度增强；若暖涡在外区，其对热带气旋的作用相反，导致热带气旋强度减弱。由于理想化试验中热带气旋静止不动，因此研究结果可能只适用于传播速度较慢的热带气旋。本研究结果有助于更好地理解热带气旋和海洋中尺度暖涡之间的相互作用，并通过引入热带气旋外区暖涡的影响助力提高热带气旋强度预报工作。     

Computational and Empirical Investigation of Propeller Tip Vortex Cavitation Noise

In this study, non-cavitating and cavitating flow around the benchmark DTMB 4119 model propeller are solved using both viscous and potential based solvers. Cavitating and non-cavitating propeller radiated noises are then predicted by using a hybrid method in which RANS(Reynolds-averaged Navier-Stokes) and FWH(Ffowcs Williams Hawkings) equations are solved together in open water conditions. Sheet cavitation on the propeller blades is modelled by using a VOF(Volume of Fiuld) method equipped with Schnerr-Sauer cavitation model.Nevertheless, tip vortex cavitation noise is estimated by using two different semi-empirical techniques, namely Tip Vortex Index(TVI, based on potential flow theory) and Tip Vortex Contribution(TVC). As the reference distance between noise source and receiver is not defined in open water case for TVI technique, one of the outputs of this study is to propose a reference distance for TVI technique by coupling two semi-empirical techniques and ITTC distance normalization. At the defined distance, the starting point of the tip vortex cavitation is determined for different advance ratios and cavitation numbers using potential flow solver. Also, it is examined that whether the hybrid method and potential flow solver give the same noise results at the inception point of tip vortex cavitation.Results show that TVI method based on potential flow theory is reliable and can practically be used to replace the hybrid method(RANS with FWH approach) when tip vortex cavitation starts.     

Experimental analysis on the risk of vortex ventilation and the free surface ventilation of marine propellers

The paper presents a discussion of the ventilation inception and air drawing prediction of ships propellers, aiming to predict under what conditions ventilation will happen, and the actual physical mechanism of the ventilation.Three different types of ventilation inception mechanisms are included in our discussion: free surface vortex ventilation, ventilation by sucking down the free surface without forming a vortex as well as ventilation by propeller coming out of the water. Ventilation prediction is based on a series of model tests, where the propeller is tested in different levels of intermittent ventilation. The use of underwater video gives a visual understanding of the ventilation phenomena.Ventilation by vortex formation has analogies with other phenomena, such as the inlet vortex in pump sumps, ground vortex at the inlet of the aircraft engines and the Propeller Hull Vortex Cavitation (PHVC). The paper includes comparison between Propeller Hull Vortex Cavitation (PHVC) and Propeller Free Surface Vortex Ventilation (PFSVV) as well as comparison between PFSVV and vortex formations of aero engines during high power operation near a solid surface.Experimental data based on several different model tests shows the boundary between the vortex forming, non-vortex forming and free surface ventilation flow regimes. For comparison the following parameters, which determined the intensity of the hydrodynamic interaction between the propeller and free surface have been used: propeller load coefficient c_T, tip clearance ratio c/D, propeller submergence ratio h/R, ambient velocity V_i and flow cavitation/ventilation number σ_cav/σ_vent.     

Numerical simulation of viscous cavitating flow around a ship propeller

In the present study, cavitation and a ship propeller wake are reported by computed fluid dynamics based on viscous multiphase flow theory. Some recent validation results with a hybrid grid based on unsteady Navier-Stokes (N-S) and bubble dynamics equations are presented to predict velocity, pressure and vapor volume fraction in propeller wake in a uniform inflow. Numerical predictions of sheet cavitation, tip vortex cavitation and hub vortex cavitation are in agreement with the experimental data, same as numerical predictions of longitudinal and transversal evolution of the axial velocity. Blade and shaft rate frequency of propeller is well predicted by the computed results of pressure, and tip vortex is the most important to generate the pressure field within the near wake. The overall results indicate that the present approach is reliable for prediction of cavitation and propeller wake on the condition of uniform inflow.     

A reliable simulation for hydrodynamic performance prediction of surface-piercing propellers using URANS method

Surface Piercing Propellers (SPPs) are a particular kind of propellers which are partially submerged operating at the interface of air and water. They are more efficient than submerged propellers for the propulsion system of high-speed crafts because of larger propeller diameter, replacing cavitation with ventilation, decreasing the torque and higher efficiency. This study presents a reliable numerical simulation to predict SPP performance using Unsteady Reynolds-Averaged Navier–Stokes (URANS) method. A numerical study on 841-B SPP is performed in open water condition. The free surface is modeled by Volume of Fluid (VOF) approach and the sliding mesh technique is implemented to model the propeller rotational motion. The sliding mesh allows capturing the process of water entry and water exit of blades. The propeller hydrodynamic characteristics, the ventilation pattern and the time history of blade loads are validated through the comparison with available experimental data. For the studied case, it was found that the common grid independence study approach is not sufficient. The grid should be elaborately generated fine enough based on the flow pattern and turbulence modeling parameters in regions near the blade's tip, trailing and leading edges and over the suction side. Details of URANS simulations including optimal time-step size based on propeller revolution rate and the required number of propeller revolutions for periodical results are presented and discussed.     

Effect of various winglets on the performance of marine propeller

The tip vortex cavitation (TVC) is an issue of increasing interest, because the TVC plays an important role in propeller radiated noise and cavitation erosion. The marine propeller with winglets, which is inspired by the winglets of airfoil, is numerically investigated in the present paper. The blade tip of newly designed propeller tilts toward the pressure side. The difference between six propellers is the change of the rake angle at r/R = 1.0. The pressure coefficient, TVC, axial velocity field and helicity are analyzed. The numerical results show that the winglets of newly designed propeller scarcely affect the efficiency of propeller. The thrust coefficient gradually decreases with the increase in rake angle. As for the suction side, the pressure coefficient (C_p) of winglets propellers is higher than the conventional propeller in general. In addition, the winglets are beneficial to generate less cavitation behavior when the rake angle is small. However, as the rake angle is further increased, the cavitation behavior of winglets propeller is also increased, even larger than the conventional propeller. Therefore, it can be deduced that the winglets can be used to effectively improve the TVC characteristics to some extent.     

Relevance of transition turbulent model for hydrodynamic characteristics of low Reynolds number propeller

The propeller of an Autonomous Underwater Vehicle (AUV) operates at low Reynolds number in laminar to turbulent transition region. The performance of these propellers can be calculated accurately using RANSE solver with γ − Re_θ transition model. In this study, the global and local hydrodynamic characteristics of open and ducted propeller are investigated using the γ − Re_θ transition model. The capability of the γ − Re_θ transition model to capture laminar to turbulent transition on the surface of the open propeller is demonstrated by comparison with published experimental results. The application of transition model for the propeller K_a-4-70 inside the duct 19A shows that the centrifugal forces are dominant at low Reynolds number and the flow is mainly directed in the radial direction. The transition model is able to predict complex flow physics such as leading-edge separation, tip leakage vortex, and the separation bubble on outer surface of the duct. The accurate prediction of these flow phenomenon can lead to correct calculation of global hydrodynamic forces and moments acting on the propeller at low Reynolds number.     

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

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

Numerical investigation on the hydrodynamic characteristics of a marine propeller operating in oblique inflow

The hydrodynamic characteristics of a marine propeller operating in oblique inflow are investigated by using CFD method. Two propellers with different geometries are selected as the study subjects. RANS simulation is carried out for the propellers working at a wide range of advance coefficients and incidence angles. The effects of axial inflow and lateral inflow are demonstrated with the hydrodynamic force on the propeller under different working conditions. Based on the obtained flow field details, the hydrodynamic mechanism of propeller operating in oblique inflow is analyzed further. The trailing vortex wake of propeller is highly affected by the lateral inflow, resulting in the deflected development path and the circumferentially non-uniform structure, as well as the enhanced axial velocity in slipstream. Different flow patterns are observed on the propeller blade with the variation of circumferential position. Combined with the computed hydrodynamic forces and pressure distribution on propeller, the mechanism resulting in the increase of propulsive loads and the generation of propeller side force is explored. Finally, a systematic analysis is carried out for the propulsive loads and propeller side force as a function of axial and lateral advance coefficients. The major terms that play a dominant role in the modeling of propulsive loads and propeller side force are determined through the sensitivity analysis. This study provides a deeper insight into the hydrodynamic characteristics of propeller operating in oblique inflow, which is useful to the investigation of propeller performance during ship maneuvers.     

Analysis of wake behind a rotating propeller using PIV technique in a cavitation tunnel

A two-frame particle image velocimetry (PIV) technique is used to investigate the wake characteristics behind a marine propeller with 4 blades at high Reynolds number. For each of 9 different blade phases from 0° to 80°, 150 instantaneous velocity fields are measured. They are ensemble averaged to study the spatial evolution of the propeller wake in the region ranging from the trailing edge to one propeller diameter (D) downstream location. The phase-averaged mean velocity shows that the trailing vorticity is related to radial velocity jump, and the viscous wake is affected by boundary layers developed on the blade surfaces and centrifugal force. Both Galilean decomposition method and vortex identification method using swirling strength calculation are very useful for the study of vortex behaviors in the propeller wake region. The slipstream contraction occurs in the near-wake region up to about X/D=0.53 downstream. Thereafter, unstable oscillation occurs because of the reduction of interaction between the tip vortex and the wake sheet behind the maximum contraction point.     

基于重叠网格方法浮筒涡激运动特性数值分析

浮筒被广泛应用在海洋工程中,研究浮筒的涡激运动对于减少其对海洋平台构件的疲劳损坏具有指导意义。传统的动网格方法在处理浮筒转动运动时会因网格变形过大导致计算不收敛,采用了重叠网格方法以解决这一问题。数值试验采用了基于开源工具包Open FOAM自主开发的naoe-FOAM-SJTU求解器。分别进行了自由衰减数值试验和涡激运动数值试验。研究表明,随着折合速度的增加,浮筒的顺流向、横流向、垂荡和艏摇运动频率增加,且顺流向与垂荡频率相近,横流向与艏摇频率相近;其次,根据涡量场分布,表明浮筒前一时刻的泻涡会影响到浮筒之后时刻的周向涡量分布;最后,研究发现自由端对于浮筒尾流场泻涡有着显著影响,为将来探究减少浮筒涡激运动的方法提供指导。