1. 引言
减小在流体中运动物体的粘性阻力一直是人们关注的热点。粘性阻力分为摩擦阻力和压差阻力。摩擦阻力对外形不敏感,早期研究工作主要是通过减小表面粗糙度 [1] [2] 来减小摩擦阻力,随着湍流理论的发展以及新材料的不断涌现,出现了种类繁多的减阻方法与思路,如消除或减小壁面黏附作用 [3] [4] 、自适应壁 [4] [5] 、优化或改变近壁面流动特性、改变近壁层流体物性等 [6] 。其中,在水中引入气体,在物体表面形成气体隔层/润滑的思想很早就有,与其它添加物质法相比,在物体表面引入空气显然具有无污染的优点。在水中引入气体有两种方式,即空化法和通气法。通气有顺流向和逆流向两种方式。顺流向通气是空气润滑意识下的自然选择,逆流向通气是在形成超空腔的指导思想下提出。逆流向喷气也可以在运动体表面形成空气润滑层,但是,需要极大的喷流速度才能获得足够的动量使得物体表面具有较高的气层覆盖率及较理想的减阻效果。前苏联著名学者Sedov提出,在物体前缘逆流喷水可以获得助推力 [7] ,其原理在于逆流喷水可使流动驻点脱离前缘,改善驻点区域高压、减小压差阻力。
基于Sedov的理论和逆流喷气减阻的思想,本文提出了逆流气液混合喷方式的构想。在物体前缘逆流喷水同时携带气体,一方面逆流喷水使流动驻点脱离前缘从而降低压差阻力,获得助推力;另一方面,与顺流喷气相比,逆流喷气既容易保持物体外形又可以在摩擦阻力最大的物体前缘形成空气润滑,且在喷射的液体的携带下,能得到更好的气层覆盖率,优势明显。此外必须要指出的是在物体前缘设置喷射装置必然要消耗能量,会某种程度上削弱逆流喷水/气的减阻节能效果。因此如何选择喷水/气口面积,喷水/气速度达到较佳的减阻节能效果是本文思考的主要方向。
据此本文针对SUBOFF潜艇模型,利用FLUENT软件进行了迎流喷水/气的减阻技术的数值模拟研究。研究不同喷水/气口面积、喷水/气速度对减阻效果的影响。
2. 数值计算模型的建立
2.1. 模型尺寸及网格划分
本文针对David Taylor研究中心与美国国防预研规划署(NSWCCD)共同发布的SUBOFF [8] 潜艇模型展开研究。潜艇模型艇体总长L为4.356 m (缩尺比为1:24),最大直径D为0.508 m。喷气/水口设置在潜艇艏部。
计算域为圆柱体,如图1所示。沿潜艇纵向,来流入口距艇首取一倍艇长L使得湍流充分发展,潜艇尾部至流域出口取3 L以使得尾流场充分发展;沿潜艇径向取8 D以避免边界的干扰。采用六面体结构化网格,在距离艇体表面径向1.5 D的范围内对网格进行加密,且由外向内逐步加密;特别在艇体首尾线型变化较大的部位进行局部加密,以期能准确捕捉潜艇周围的流动特性。总网格数约90万。
Figure 1. The sketch map computational domain and cells
图1. 潜艇计算域与网格示意图
2.2. 边界条件及数值计算方法
流域入口为速度入口,来流速度为5.144 m/s,即雷诺数为2.32 × 107;出口为压力出口边界条件;潜艇表面为无滑移边界;顶部和底部边界为对称边界。
湍流模型采用RNG两方程模型,压力与速度的耦合使用SIMPLEC算法;采用有限体积法离散控制方程;水汽交界面采用VOF方法跟踪。
2.3. 减阻率与节能率的计算方法
在潜艇艏部逆流喷水可减小压差阻力,使潜艇运动消耗功率减小;但另一方面,逆流喷射的设置必然要消耗掉能量,使潜艇运动消耗功率增加。因此逆流喷射最终是否能够达到减阻节能的效果,要综合考虑这两方面的影响。
如果不喷射流体,则原潜艇受到的粘性总阻力Rt0为:
(1)
式中Rf0为潜艇摩擦阻力、RP0为压差阻力;分别由潜艇表面的切应力与压力积分求和得出。
在潜艇艏部设置喷口,逆流喷射流体时,考虑到喷射的流体会对喷射管道等装置产生力的作用,再通过力的传导,对潜艇产生一定影响。而喷射的流体会对喷射管道等装置产生的力主要与喷口处流体压力以及流体冲量有关 [1] ,故将喷口压力以及流体冲量计入到潜艇的压差阻力中,所以设置喷射后潜艇的总压差阻力RtPj,其计算式为:
(2)
式中RP根据潜艇壁面压力积分计算;Ri 表示喷口处的压力;RC为喷口喷出的流体对潜艇的冲量,其计算式为:
(3)
式中ρj为喷射流体的密度,vj为喷射流体的速度,Sj为喷口面积。
则有喷射潜艇总阻力Rtj为:
(4)
式中Rfj为有喷射潜艇表面的摩擦阻力。
故减阻率μ为:
(5)
不喷射流体时潜艇航行消耗的功率Ps0为:
(6)
式中V为潜艇航行速度。
有喷射流体潜艇运动所消耗的总功率Psj为:
(7)
式中E为单位时间通过喷口的流体能量,其计算式为
(8)
其中Q表示通过喷口的流体流量。
因此节能率η的计算式为:
(9)
3. 水下航行体迎流单相喷水数值模拟
根据前述计算方法,为了便于对比分析,首先对不喷射时原潜艇模型流场进行了数值模拟,得到了其阻力值。其次对潜艇前缘迎流单相喷水进行了数值模拟,确定了较佳的喷水口面积以及喷水速度。在此基础上结合气液两相流的VOF模型,对潜艇前缘迎流气液同喷进行了数值模拟研究,得出喷气口面积、喷气流量对潜艇流场特性以及阻力的影响规律。
3.1. 喷水口设置
为了讨论喷口面积对阻力性能的影响,在艇体前缘分别设置面积不同的圆形喷口,如图2所示。喷口面积见表1。
Figure 2. The sketch map of submarine at bow
图2. 潜艇艏部喷口局部放大示意图
分别计算不同喷口面积,不同喷水量时的流场,具体计算条件为:首先固定喷水口面积(Sw = 0.0004 m2),改变喷水速度(流量);然后固定喷水速度(3 m/s),改变喷水口面积。
3.2. 水下航行体迎流单相喷水计算结果分析
不喷射流体时,原型潜艇模型的摩擦阻力Rf0为217.65 N,压差阻力RP0为28.02 N,粘性总阻力Rt0为245.67 N,其运动消耗功率PS0为1263.73 W。
同一喷水口面积(0.0004 m2)不同喷水速度以及同一喷水速度(3 m/s)不同喷口面积时潜艇表面压力沿艇身分布如图3,图4所示。
由图3可知,沿整个艇身,喷水与不喷水两种情况下压力值基本重合,迎流喷水对潜艇艇身的压力分布影响不大;迎流喷水使得艏部压力明显降低,且喷水速度越大,压力降低越多。不喷水时,潜艇艏部前端为驻点,该处压力值最大;在潜艇前缘迎流喷水导致艏部驻点脱离潜艇表面且前移,潜艇前端不再是驻点,故其艏部压力有明显降低。
而图4则表明喷水速度不变,仅改变喷水口面积时,对潜艇最前缘的压力影响不大,最前缘压力的大小与喷水速度的关系较大,与喷水口面积的大小基本无关。
此外,从图3、图4中均可看出,喷水区域过后,压力有一个明显下降再回升的变化,反应出漩涡产生低压现象,这些都对减阻产生有利的影响。
不同喷水速度,不同喷水面积潜艇前缘迎流单相喷水数值模拟结果如表2和表3所示。
Figure 3. Pressure with longitudinal direction x of submarine at different injection velocity
图3. 不同喷水速度潜艇表面压力沿潜艇纵向分布
Figure 4. Pressure with longitudinal direction x of submarine at different spouts area
图4. 不同喷水口面积潜艇表面压力沿潜艇纵向分布
Table 2. The result of simulation at different water injection velocity
表2. 不同喷水速度时数值模拟结果(Sw = 0.0004 m2)
Table 3. The result of simulation at different spouts area (vw = 3 m/s)
表3. 不同喷口面积的数值模拟结果(喷水速度为3 m/s)
从表2可知,潜艇在5.144 m/s来流速度下,固定喷口面积(0.0004 m2),随着喷流速度的增加,潜艇总阻力呈现先减后增的趋势;迎流喷水主要改善的是潜艇艏部压力但摩擦阻力也略有降低,潜艇减阻率μ、节能率η随喷流速度的增加呈现先增后减的规律,在喷水速度为3 m/s时有最大值2.61%与2.19%。
表3表明,潜艇在5.144 m/s来流速度下,固定喷水速度(3 m/s),总阻力随着喷口面积的增加而逐步降低,减阻率μ随着喷口面积的增加显著提高,在喷水口面积为13.83 × 10−3 m2时,减阻率可达60%。
但随着喷口面积的增加,喷射的流体对潜艇的反冲力RC以及单位时间内通过喷口的流体能量E也随之增加。为了更清楚的揭示不同喷流面积下减阻率与喷射该流量的水所需付出的能量的关系,引入减阻效能比这一参数,它表示付出单位喷射流体需消耗的能量后所获得的减阻率,计算结果表明减阻效能比(减阻率/E)随喷口面积的增加逐渐降低。
因此可以得出,尽管喷口面积越大,减阻率越高减阻效果越好,但所耗费的能量较多,故并非喷口面积越大越好。此外,也应考虑到在实际情况下,潜艇前端开口对结构强度、流噪声的影响等等,喷口面积也不宜过大。故选取喷水口面积SW = 1.57 × 10−3 m2 ,喷水速度vW = 3 m/s,来进行气液混喷的数值研究。
4. 水下航行体迎流气液混喷数值模拟
迎流喷水可以减小潜艇的压差阻力,但其对潜艇的摩擦阻力影响较小,而在SUBOFF裸壳潜艇的总阻力中摩擦阻力占总阻力的88.6%,是其总阻力的主要组成部分。为减小潜艇的摩擦阻力,考虑在潜艇艏部设置圆形/圆环形喷口,内部圆形喷口喷水,外部圆环形喷口喷气,这样所喷射的液体将气体包裹在潜艇表面,以期减小潜艇表面的摩擦阻力。
艏部中心圆形喷口喷水,根据前文结果,喷口面积SW选取为1.57 × 10−3 m2,分别计算4种不同喷气口面积(SA1 = 1.94 ×10−3 m2, SA2 = 4.64 ×10−3 m2, SA3 = 8.10 × 10−3 m2, SA3 = 12.16 × 10−3 m2),每种面积计算7种不同喷气量共计28种的两相流流场。
4.1. 水下航行体迎流气液混喷压力分布
喷气口面积一定(SA3 = 8.10 × 10−3 m2),不同喷气量以及喷气流量一定(5.83 × 10−3 m3/s)不同喷气口面积下潜艇表面压力分布分别如图5和图6所示。
由图5和图6可以看出,不论是固定喷气口面积改变喷气流量,还是固定喷气流量改变喷气口面积,对潜艇艇身的压力分布的影响都不太大,仅对潜艇艏部喷气口附近的压力有所影响。
Figure 5. Pressure with longitudinal direction x of submarine at different air flow rate (SA3 = 8.10 × 10−3 m2)
图5. 不同喷水口面积潜艇表面压力沿潜艇纵向分布(SA3 = 8.10 × 10−3 m2)
Figure 6. Pressure with longitudinal direction x of submarine at different spouts area (Q = 5.83 × 10−3 m3/s)
图6. 不同喷口面积潜艇表面压力分布(Q = 5.83 × 10−3 m3/s)
图5表明喷气面积固定的情况下,改变喷气流量对潜艇表面压力影响很小,不同喷气流量下潜艇表面的压力曲线基本是重合的。并且都出现了压力下降再回升的变化,这主要是因为迎流喷水在艇艏处产生了漩涡低压区。
由图6可见喷气流量固定的情况下,改变喷气口面积同样主要影响喷口附近的压力。喷气口面积较小时,压力有更加明显的下降再回升的变化。
图7给出了喷气流量为5.83 × 10−3 m3/s,不同喷气面积时潜艇艏部的流速矢量图。
Figure 7. Velocity of submarine at different area of jet exhaust (Q = 5.83 × 10−3 m3/s)
图7. 不同喷气口面积时潜艇艏部速度矢量图(Q = 5.83 × 10−3 m3/s)
艏部迎流喷射的流体在来流的作用下,会顺流向后流动,同时由于流体的粘性会带走部分艏部的液体,周围的流体要补充这部分被带走的流体所占据的空间,会在潜艇前缘产生漩涡。而迎流喷气不断增加喷气口面积时,迎流喷出的气体不断补充了部分被带走的流体所占据的空间,从而使得艏部因迎流喷水产生的漩涡区的到改善,流场不那么紊乱,且喷气口面积越大时这种现象越明显。
4.2. 水下航行体迎流气液混喷阻力计算结果与分析
表4给出了喷气口面积为SA4 = 12.26 × 10−3 m2不同喷气量的计算结果,其余三种喷气口面积,不同喷气量的计算结果与此类似,在此不再赘述。
Table 4. The result of simulation at different air flow rate (SA4 = 12.26 × 10−3 m2)
表4. 不同喷气量计算结果(SA4 = 12.26 × 10−3 m2)
同一喷气流量下(Q = 13.24 × 10−3 m3/s)不同喷口面积时的计算结果如表5所示。
Table 5. The result of simulation at different the area of air injection (Q = 13.24 × 10−3 m−3/s)
表5. 不同喷气口面积计算结果(Q = 13.24 × 10−3 m−3/s)
从表4和表5可以看出,不论喷气流量与喷气口面积的大小,潜艇前缘气液同喷时减阻效果都优于单相喷水情况,减阻率、节能率都比单相喷水时大。
喷气口面积不变时,随着喷气流量的增加,摩擦阻力基本保持相同的梯度显著降低。相比不喷气情况,随着喷气流量的增加,摩擦阻力最大可降低约24%;总阻力随着喷气流量的增加而减小,相比仅迎流喷水情况,气液混喷时总阻力最大可降低约20%,而相比不喷水不喷气时,气液混喷时总阻力最大可降低约30%。
喷气面积不变时,减阻率与节能率都随着喷气流量的增加而增加。除了喷口面积为SA1 的情况减阻率与节能率在喷气流量较大时变化稍平缓外,其余喷口面积的流量增大约7倍时,减阻率与节能率可增加约20%左右。整体而言,迎流气液同喷的情况下,不同喷气口面积下减阻率与节能率可达到10%~30%。
喷气流量不变时,增加喷气口的面积,摩擦阻力同样有所减小,但没有随流量的增加减小的显著,喷口面积增大约6倍时,摩擦阻力最大可降低约11%。总的来说,迎流气液同喷对减小摩擦阻力效果显著;总阻力也随着喷气口面积的增大而降低,但没有随流量的增加减小的显著,喷口面积增大约6倍时,总阻力最大可降低约12%。喷气流量不变时,减阻率与节能率同样随着喷口面积的增加而增加。
综上所述,不论是增加喷气口面积还是增加喷气流量,对减小潜艇的总阻力都是有利的,且都明显优于迎流单相喷水时的阻力性能。在实践中,可以结合实际情况,综合考虑喷口对潜艇结构强度、总布置以及流噪声等方面的影响,选取合适的喷口面积与喷气流量的组合来达到降低潜艇总阻力的目的。
5. 结论
本文基于Sedov理论,提出了在物体前缘迎流喷水同时携带气体的构想。并针对SUBOFF潜艇模型,在来流速度为5.144 m/s,雷诺数为2.32 × 107 的条件下,利用FLUENT软件,对潜艇前缘单相喷水以及气液同喷进行了数值模拟研究。研究结果表明:
1) 迎流喷射流体可以使得潜艇艏部驻点脱离潜艇表面前移,从而减小其压差阻力。迎流气液混合喷既可使流动驻点脱离前缘从而降低压差阻力,又可在潜艇表面形成空气润滑,达到减少摩擦阻力的目的。
2) 单相喷水且喷水面积固定时,减阻率以及节能率随喷水速度的增加呈先增后减的变化特点。在喷水速度为3 m/s时减阻率与节能率有最大值,分别为2.61%与2.19%。喷水速度固定时,单相喷水潜艇阻力随喷水面积的增加而减小,减阻率随着喷水面积的增加显著增加,在喷口面积为13.83 × 10−3 m2时减阻率可达60.62%。但减阻效能比随着喷水面积的增加,呈现递减的变化特点。
3) 气液同喷且喷气口面积不变时,摩擦阻力随着喷气流量的增加显著降低。相比仅单相喷水情况,摩擦阻力最大可降低约24%,总阻力最大可降低约20%;喷气流量不变时,摩擦阻力最大可降低约11%,总阻力最大可降低约12%。整体而言,迎流气液同喷的情况下,减阻率与节能率可达到10%~30%。
数值模拟研究表明潜艇前缘迎流喷水/气减阻性能良好,但想要将这种减阻技术成功的应用于实践中,其节能效果尚需实验验证;此外在潜艇前缘开口势必对潜艇结构强度以及潜艇噪声等产生影响,这些都需要进一步的研究。
基金项目
国家自然科学基金(51309182)。