2. 华中科技大学 船舶与海洋工程学院,湖北 武汉 430074
2. School of Naval Architecture and Ocean Engeering , Huazhong University of Science and Technology, Wuhan 430074, China
螺旋桨作为潜艇的主要推进装置,深入了解潜艇螺旋桨下游尾流场的流动特性,对于评估潜艇在水下的推进性能至关重要。国内外围绕螺旋桨的流动研究集中在包括螺旋桨空化[1-2]、螺旋桨入流[3-4]、螺旋桨自由表面-波浪相互作用[5]、和螺旋桨-方向舵相互作用[6]在内的广泛方面。
近年来,随着CFD技术的发展,使得通过数值模拟方法研究螺旋桨尾流精细流场特性成为可能。Kumar等[7]采用大涡模拟方法研究了尾流不稳定机理,讨论了螺旋桨尾流从近场到远场的轴向演化。李鹏程等[8-9]则采用RANS 和LES方法分别对螺旋桨轮毂涡进行了数值模拟,分析了毂涡强度与桨毂环量的关系,基于Oseen涡的切向速度分布形式建立了毂涡模型。王恋舟等[10]对比了不同湍流模型下螺旋桨尾流动力学特性,揭示了螺旋桨尾流不稳定性的触发机理,提出了螺旋桨梢涡演化模型。Posa等[11]采用LES模拟对比了不同载荷下螺旋桨尾流涡旋结构和湍动能轴向演变,研究了涡旋结构对湍动能变化的影响。
本文采用LES方法计算不同载荷条件下螺旋桨尾流场水动力特性,验证了计算方法的准确性,分析了指定载荷条件下螺旋桨尾流梢涡和随边涡的发展演变,以及梢涡和轮毂涡强度关系,揭示了近场梢涡和轮毂涡湍动能产生的原因,以及湍动能的变化规律。
1 数值计算方法 1.1 数学模型大涡模拟理论核心思想是通过过滤函数将大尺度涡和小尺度涡分离开,大尺度涡进行直接模拟,小尺度涡用模型来封闭。滤波后的连续性方程和N-S方程表达式为:
$ \frac{{\partial \rho }}{{\partial {\text{t}}}} + \frac{\partial }{{\partial {x_i}}}\left( {\rho \overline {{u_i}} } \right) = 0,$ | (1) |
$ \frac{\partial }{{\partial t}}\left( {\rho \overline {{u_i}} } \right) + \frac{\partial }{{\partial t}}\left( {\rho \overline {{u_i}} \overline {{u_j}} } \right) = \frac{{\partial {\sigma _{ij}}}}{{\partial {x_j}}} - \frac{{\partial \overline p }}{{\partial {x_i}}} + \frac{{\partial {\tau _{ij}}}}{{\partial {x_j}}} 。$ | (2) |
式中:
$ {\tau _{ij}} = \rho \overline {{u_i}{u_j}} - \rho \overline {{u_i}} \overline {{u_j}}。$ | (3) |
需用亚格子应力模型进行模拟,本文采用壁面自适应的局部涡粘性 (WALE)模型[12] 。
1.2 计算模型本文研究对象E1619桨是意大利INSEAN水池为潜艇设计的7叶大侧斜螺旋桨[13],直径D为0.26 m,毂径比dh/ D为0.226。计算域与边界条件如图1所示。
![]() |
图 1 E1619螺旋桨计算域 Fig. 1 Calculation domain of propeller E1619 |
其中螺旋桨前方6D处为速度入口,后方12D处为压力出口,螺旋桨表面为无滑移壁面,距螺旋桨6D处为对称平面。
1.3 计算参数采用STAR-CCM+软件结合2 046万非结构网格,首先采用RANS中的
采取固定转速、改变来流流速的方法,取进速系数J的值为0.56、0.7、0.84,将仿真结果与已有的实验数据进行了比较,计算结果如图1所示(为避免计算结果不收敛,取螺旋桨旋转15圈流场稳定后的时均统计值)。其中进速系数J、推力系数KT、扭矩系数KQ定义如下:
$ J = \frac{{{U_0}}}{{nD}} ,$ | (4) |
$ {K_T} = \frac{T}{{\rho {n^2}{D^4}}} ,$ | (5) |
$ {K_Q} = \frac{Q}{{\rho {n^2}{D^4}}} 。$ | (6) |
可以看出大涡模拟很好捕捉了螺旋桨推力和扭矩系数,计算误差都在5%以内,结果比较理想。而实验与计算结果存在差异,可能是由计算建模螺旋桨与实际模型存在差异、离散格式、网格精细、求解方法等原因造成的。
![]() |
图 2 敞水螺旋桨水动力系数 Fig. 2 Hydrodynamic coefficient of propeller in open water |
对于进速系数J=0.7工况,图3给出了螺旋桨尾流在z=0截面上的轴向速度分布云图(无量纲化ux/U0)。在图3中,尾流轴向流动的加速是明显的,由于轮毂的几何形状,在边界层上产生了逆压力梯度,导致在内半径处产生了很大的速度缺陷区域和一个相当强的轮毂涡。除了轮毂涡区,螺旋桨尾迹的轴向速度比自由流高,这使得它不同于其他转子尾迹,如风力涡轮机等。图4给出的压力场(无量纲化Cp=P−P0/0.5 ρU02)等值线轮廓上,可以清楚地看到涡核是一个低压区域,尤其是轮毂涡内部的区域压力最低,因此更容易发生空化。
![]() |
图 3 轴向速度分布云图 Fig. 3 Axial velocity distribution nephogram |
![]() |
图 4 压力等值线分布云图 Fig. 4 Pressure contour distribution nephogram |
图5为螺旋桨尾流在z=0截面上的涡量幅值分布云图(无量纲化|ω|D/U0)。可以看出,涡量场在近场中显示出明显的梢涡和随边涡,随着向下游发展,这些旋涡结构变得不稳定,最终破裂形成远尾流。而在梢涡向下游发展的过程中,在(Ⅰ)处存在明显的梢涡互感与合并现象。为了突出尾流梢涡梢涡互感不稳定性触发与合并过程,图6采用Q准则(无量纲等值面QD2/U02=20),显示了由无量纲轴向速度着色的三维涡结构。为了便于分析对梢涡进行了编号,可以看到7号梢涡与1号梢涡、2号梢涡与3号梢涡、4号梢涡与5+6号合并梢涡在(Ⅰ)处的合并过程。由无量纲轴向速度大小分析可知,相邻梢涡发生互感后,其合并的本质是相邻的下游梢涡将上游梢涡推向尾流轴向速度较高的内半径区域,最终包裹上游梢涡使其靠近下游梢涡完成合并,王恋舟等[9]对E1658螺旋桨使用延迟分离涡模拟(DDES)计算结果中观察到类似的梢涡合并现象。
![]() |
图 5 涡量幅值分布云图 Fig. 5 Vorticity magnitude distribution nephogram |
![]() |
图 6 无量纲轴向速度着色的涡量等值面,QD2/U02=20 Fig. 6 Non-dimensional Isosurface QD2/U02=20,countour of non-dimensional axial velocity |
图7为螺旋桨下游x/D=−0.1,−0.25时涡量场的轴向分量云图(无量纲化 ωxD/U0)。不难发现,轮毂涡流比梢涡更强。这是由于轮毂涡流是由桨毂涡流和叶片根部涡流合并产生的,因此其强度是梢涡的N倍。在螺旋桨近场尾流中,轮毂涡流速度和压力更低,导致其耗散速度明显较梢涡更快,这也意味着N在不断减小,从x/D=−0.1截面上轮毂涡强度最高是梢涡的5倍,减小到x/D= −0.25截面上的3倍。同时,轴向涡流有很明显的螺旋形状,并且梢涡与轮毂涡螺旋方向相反,即梢涡螺旋方向为逆时针,轮毂涡螺旋方向为顺时针。这在图8给出的涡量在z=0截面上的分量(无量纲化ωzD/U0)也可以看到,梢涡与轮毂涡的符号相反。在图8中,除了能看到梢涡的互感和合并现象之外,还发现梢涡与相邻上游梢涡的随边涡符号相同,这导致其相互吸引。随着尾流向下游的发展,由于梢涡互感而合并产生的梢涡强度更强,会吸引随边涡更加靠近其轴向平面。
![]() |
图 7 涡量场不同位置的轴向分量云图 Fig. 7 Axial component of vorticity at different positions |
![]() |
图 8 涡量的z向分量分布云图 Fig. 8 z-component of vorticity distribution nephogram |
图9给出了螺旋桨尾流湍动能轴向演变(无量纲化k/U02)。图中螺旋桨近场尾流梢涡湍动能值较大,这是由于桨叶尾流剪切层所造成的。外半径尾流处湍流动能的流向演变受梢涡的影响很大,大约在x/D=−1范围内,梢涡非常连贯,与螺旋桨叶片的旋转同步。因此,湍动能的局部最大值发生在梢涡的核心,主要与其中的速度波动有关。而在x/D=−1附近,明显缺失了一个梢涡湍动能极大值。这是因为该处梢涡与下游梢涡发生互感,与下游梢涡合并,导致该处梢涡断裂,速度脉动减小,缺少了极大值。同理在x/D=−1.7附近,梢涡湍动能极大值明显较上游梢涡极大值更大。这是因为该处梢涡是与上游梢涡合并所产生的,导致速度脉动变大,湍动能增大。而在x/D=−3后,尤其是x/D=−3.4后,湍动能发生明显增大,造成这个现象的原因有2个:一是梢涡完全分解,产生了较大的湍动能,这从图10给出的x/D=−3处截面涡量幅值和湍动能可以看出;二是x/D=−3.4后网格分辨率不足,导致涡旋耗散过多,湍动能增大。在轮毂涡流中,随着向下游移动,湍动能先减小,然后增大。螺旋桨附近轮毂涡流的湍动能值较高是轮毂上尾流分离产生的不稳定性大尺度波动造成的。轮毂涡流向下游移动时,湍动能急剧下降,随后在梢涡完全失稳后,湍动能增加。梢涡失稳导致轮毂涡流振荡,从而产生湍流。这一现象支持了Felli等[12]的假设,即螺旋桨尾流中的梢涡和轮毂涡流不稳定性之间存在因果关系。
![]() |
图 9 湍动能分布云图 Fig. 9 Turbulent kinetic energy distribution nephogram |
![]() |
图 10 x/D=−3处涡量幅值和湍动能分布云图 Fig. 10 Vorticity magnitude and turbulent kinetic energy at x/D=−3 |
本文基于LES湍流模型,计算了不同工况下螺旋桨敞水水动力特性,通过与实验数据对比,验证了计算方法的准确性,对指定工况下七叶大侧斜螺旋桨尾流场、涡旋特性、湍动能进行了分析。结果表明:
1)尾流轴向流动加速是明显的,由于轮毂的存在,产生了一个速度缺陷区和强烈的轮毂涡,并且在压力场可见涡核都处于低压区域。
2)在向下游发展的过程中,梢涡会发生互感,其合并的本质是下游梢涡将上游梢涡推向尾流轴向速度较高的内半径区域。轮毂涡与梢涡螺旋方向相反,且轮毂涡强度明显更高,这是由于其包含了叶根涡。梢涡与上游相邻随边涡符号相同,而相互吸引。
3)近场梢涡湍动能较高是桨叶尾流剪切层所造成的,轮毂涡的湍动能值较高是轮毂上尾流分离产生的不稳定性大尺度波动造成的。在向下游发展的过程中,梢涡失稳导致轮毂涡流振荡,湍动能增大,从而产生湍流。
[1] |
杨琼方, 王永生, 张志宏. 非均匀进流对螺旋桨空化水动力性能的影响[J]. 水动力学研究与进展A辑, 2011, 26(5): 538-550. YANG Q F, WANG Y S, ZHANG Z H. Effects of non-uniform inflow on propeller cavitation hydrodynamics[J]. Hydrodynamics Research and Progress A, 2011, 26(5): 538-550. DOI:10.3969/j.issn1000-4874.2011.05.004 |
[2] |
ZHU Z. Numerical study on characteristic correlation between cavitating flow and skew of ship propellers[J]. Ocean Engineering, 2015, 99: 63-71. DOI:10.1016/j.oceaneng.2014.12.023 |
[3] |
HIGGINS R J, ZAREV A, BARAKOS G N, et al. Numerical investigation of a two-bladed propeller inflow at yaw[J]. Journal of Aircraft, 2020, 57(2): 292-304. DOI:10.2514/1.C035647 |
[4] |
RAD R G, SHAFAGHAT R, YOUSEFI R. Numerical investigation of the immersion ratio effects on ventilation phenomenon and also the performance of a surface piercing propeller[J]. Applied Ocean Research, 2019, 89: 251-260. DOI:10.1016/j.apor.2019.05.024 |
[5] |
CARRICA P M, SADAT-Hosseini H, STERN F. CFD analysis of broaching for a model surface combatant with explicit simulation of moving rudders and rotating propellers[J]. Computers & Fluids, 2012, 53: 117-132. |
[6] |
WANG L, GUO C, XU P, et al. Analysis of the wake dynamics of a propeller operating before a rudder[J]. Ocean Engineering, 2019, 188: 106250.
|
[7] |
KUMAR P, MAHESH K. Large eddy simulation of propeller wake instabilities[J]. Journal of Fluid Mechanics, 2017, 814: 361-396. DOI:10.1017/jfm.2017.20 |
[8] |
李鹏程, 刘登成, 董郑庆, 等. 基于定常RANS方法的螺旋桨毂涡结构分析[J]. 舰船科学技术, 2020, 42(9): 3439-39. LI P C, LIU D C, DENG Z Q, et al. Analysis of hub vortex structure based on steady RANS method[J]. Ship science and technology, 2020, 42(9): 3439-39. DOI:10.3404/j.issn.1672-7649.2020.05.007 |
[9] |
李鹏程, 周伟新, 刘登成, 等. 基于大涡模拟的螺旋桨毂涡结构分析与建模[J]. 中国造船, 2020, 61(3): 131-139. LI P C, LIU D C, DENG Z Q, et al. Modeling of propeller hub vortex structure based on large eddy simulation[J]. Shipbuilding of China, 2020, 61(3): 131-139. DOI:10.3969/j.issn.1000-4882.2020.03.015 |
[10] |
王恋舟, 吴铁成, 郭春雨. 螺旋桨梢涡不稳定性机理与演化模型研究[J]. 力学学报, 2021, 53(8): 2267-2278. WNG L Z, WU T C, GUO C Y. Study on instability mechanism and evolution modelof propeller tip vortices[J]. Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(8): 2267-2278. DOI:10.6052/0459-1879-21-151 |
[11] |
POSA A, BROGLIA R, FELLI M, et al. Characterization of the wake of a submarine propeller via large-eddy simulation[J]. Computers & fluids, 2019, 184: 138-152. |
[12] |
NICOUD F, DUCROS F. Subgrid-scale stress modelling based on the square of the velocity gradient tensor[J]. Flow Turbulence and Combustion, 1999, 62(3): 183-200. DOI:10.1023/A:1009995426001 |
[13] |
CHASE N, CARRICA P M. Submarine propeller computations and application to self-propulsion of DARPA suboff[J]. Ocean Engineering, 2013, 60: 68-80. DOI:10.1016/j.oceaneng.2012.12.029 |
[14] |
FELLI M, CAMUSSI R, FELICE D F. Mechanisms of evolution of the propeller wake in the transition and far fields[J]. Journal of Fluid Mechanics, 2011, 682: 5-53. DOI:10.1017/jfm.2011.150 |