黄海龙,王 聪,余德磊,张 鹤,卢佳兴
(1.北京宇航系统工程研究所,北京 100076; 2.哈尔滨工业大学 航天学院,哈尔滨 150001)
射弹并联入水是指两个或多个射弹沿空间平行线同时或在极短时间间隔内穿越自由水面进入水中的过程.射弹在入水过程中,当流体低压区压力降低至水的饱和蒸汽压时,发生空化形成空泡,空泡向后发展覆盖整个弹体,形成超空泡[1].处于超空泡包裹之中的射弹流体阻力急剧下降,因此超空泡射弹得到国内外学者的广泛关注[2].
国内外对超空泡射弹的研究手段主要集中于理论、数值和试验研究.早期的理论研究大多基于Logvinovich等[3]的空泡截面独立膨胀原理对空泡发展过程进行描述.Lee等[4]利用能量守恒原理对射弹高速入水空泡生成、发展过程和闭合特性进行研究.Al’ev[5]基于细长体理论研究了细长锥形体在考虑压缩性影响时的空泡特征.文献[6-7]在计算中应用了均值多相流模拟多相流场,给出了稳态和瞬态情况下轴对称流动、通气和自然超空泡形态的计算.国内在理论方面研究较早,早期陈九锡等[8]利用有限差分法对N-S方程直接求解,获得了平头射弹垂直入水时的空泡生成、发展、闭合过程及流场变化过程.马庆鹏等[9]对射弹高速入水时的弹道特性、空泡形态等进行了数值模拟研究.在试验方面,Mostafa等[10]通过水平发射入水试验发现不同的空化数下超空泡射弹形成的空泡形态不同.Hrubes[11]通过高速摄影技术获得超空泡射弹的水下运动情况,研究了射弹在运动过程中的尾拍现象.文献[12-14]系统开展大量超空泡射弹试验研究,对实验中获取的空泡形态进行测量分析,总结出由一定适用范围的超空泡形态计算的经验公式.曹伟等[15]对空化器、弹径不同的多种弹丸进行射弹试验研究.金大桥等[16]利用轻气炮射弹实验装置开展了通气超空泡水下射弹试验研究.赵成功等[17]开展细长体模型水下高速运动的试验研究.
国内外的研究大多基于单个射弹,随着 “饱和攻击”[18]的提出,多发及队列射弹入水逐渐被世界各国所重视.目前尚未从国外公开文献中查到多射弹入水的研究.国内的相关研究也极少.王志东等[19]和孙亚亚[20]对连发射弹进行了数值模拟研究. 本文将通过高速摄像方法和高压动力并联发射装置进行双射弹的上、下方向并联入水试验,研究了双射弹并联入水过程中的双空泡演化特征及其机理.
试验系统如图1所示.试验系统主要由动力加载系统、发射控制系统、光学测试系统和试验水槽及防护系统组成.其中发射系统主要包括高压气瓶、高压气室、双炮管、高压管路系统;发射控制系统主要由电磁阀与控制器组成,控制器能够提供电信号触发高速摄像和电磁阀,从而保证气室触发装置和摄像触发装置同步,能够更有效的拍摄射弹水下运动过程;光学测试系统包括高速摄像系统、照明系统及炮口通断靶测速系统;开放式水槽尺寸为1.2 m×1.2 m×10.0 m,顶部开口,前侧面采用较高强度的钢化玻璃,底部和其余侧壁由钢板和相应支架组成;保护及回收系统包括多层钢丝网和回收池等,保护及回收系统布置于水槽末端,当射弹达到水槽末端时,会撞击上钢丝网减速并进入回收池.
图1 试验系统总体方案
试验中采用两种截锥型射弹和一种圆柱型射弹. 其中截锥射弹模型如图2所示.3种射弹模型的具体尺寸见表1.试验中通过调整发射管间距改变射弹间距;通过加载不同的气室压强改变射弹入水初速度.
图2 截锥型射弹模型示意
表1 射弹模型尺寸
射弹空泡形态主要受空化数σ和弗劳德数Fr影响,对于本实验,由于射弹速度较高,弗劳德数Fr较大,即重力影响较小,在本文研究中予以忽略,故本实验空泡形态主要由空化数σ决定,定义如下:
式中:p为流场的参考压力;pv为当前温度下的饱和蒸汽压,对于水的空化,即为水的饱和蒸汽压;u为流场的来流速度,本文中即为射弹速度. 双发射弹入水时其速度存在较小的差异,因此定义双射弹初始空化数时采用其平均速度计算空化数表2列出了部分试验工况.
表2 并联入水试验部分工况
图3展示了工况1的双空泡演化过程.不同于单独入水时空泡的自由发展,并联入水时射弹空泡发展受到耦合影响,另一发射弹对空泡的附近流场形成扰动,甚至直接破坏空泡界面. 图3对双空泡的附近流场和空泡界面扰动现象进行了细致研究.在2t0时刻,下方射弹率先入水形成较大的入水空泡并且将上方射弹包含在内,此时上方射弹在下方射弹空泡内受到的阻力作用很小,这将导致上方射弹追上下方射弹,在穿过下方射弹的空泡壁面后产生新的入水空泡.在3t0~4t0时间内,上方空泡迅速扩张,下方射弹空泡壁面被上方射弹破坏后再次受到空泡扩张扰动,导致界面严重失稳且伴随着大量的汽水混合物. 另外,在3t0时刻下方射弹空泡下侧出现凸起的流场结构,分析其产生原因可能是上方射弹空泡扩张时对下方射弹空泡界面流场造成强扰动.在7t0时下方射弹头部空泡变得光滑,说明其受到的扰动减弱.上方射弹空泡后半部分尺寸“过度膨胀”,这种现象是由于空泡下方存在大量的脱落气团,使得空泡扩张时受到的流场压力较小.运动至12t0时射弹偏移较大,双空泡发展时的相互影响减弱,此时下方射弹空泡变得透明且在尾部闭合,上方射弹空泡混合了下方射弹部分泄露气体,其体积进一步增大.至20t0时,下方射弹远离上方射弹空泡发展区域,上方射弹空泡不断被拉长后由于泡内压力降低,在环境压力作用下收缩溃灭.
图3 双空泡演化过程中的完整穿越现象
图4为工况2的双空泡演化过程. 与图3不同,先入水射弹产生的入水空泡直径尺寸较小,在相同间距下未能完全包裹后发射弹.在运动过程中上方射弹空泡较短时间内维持部分在空泡内,部分在水域中,如3t0时刻,上方射弹部分与下方射弹的空泡壁面接触,部分穿越时上方射弹持续冲击着下方射弹空泡壁面,导致下方射弹空泡表面产生不透明的汽水混合物. 运动至6t0时刻,上方射弹空泡扩张且空泡表面较透明,在运动时其空泡扩张处于主导地位,趋向于挤压下方射弹空泡后得到扩张. 随着射弹继续运动,在8t0时刻,上方射弹空泡尾部尺寸较前方略大,分析形成原因是由于下方射弹空泡内部从空气中得到的补充气体迅速减小,当下方射弹空泡体积增大时泡内压力降低,所以上方射弹尾部区域射弹空泡处于低压区,空泡尺寸较自由发展时大,而在射弹前部双空泡耦合剧烈(图4中椭圆标识),其膨胀受到抑制后空泡尺寸较单空泡发展小.在这种发展模式下形成了上方射弹前部分尺寸较小,尾部区域尺寸较大的波浪状轮廓. 运动至13t0时刻,下方空泡内部压力降低,上方射弹空泡在扩张过程中占据优势,空泡轮廓偏向椭圆,下方射弹空泡不断以气团形式向两侧脱落.至18t0时刻,双发射弹的耦合程度减弱,双空泡相对完整.此时先发的下方射弹的脱落空泡呈发散状发展,后发的上方射弹脱落空泡基本单股发展.
图4 双空泡演化过程中的部分穿越现象
综上所述,并联入水过程中双空泡的穿越现象产生条件为先入水空泡产生的空泡尺寸较大,造成后发射弹入水后全部或部分处于空泡内部.在本文中保持射弹间距不变,改变了射弹空化器直径,发现双空泡由完全穿越过渡到了部分穿越.由试验图片还观察到了穿越过程中出现的“空泡扩张冲刷”和“过度膨胀”现象.
图5为B型射弹在不同间距下的流场结构.从图5中可以观察到双空泡之间的吸引现象,即同单独入水时的脱落空泡运动相比,并联入水时射弹脱落空泡朝着相互靠近的方向运动,图5中先入水的下方射弹脱落空泡的偏移现象较明显.图5中的虚线表示单独入水时脱落空泡的发展趋势,可以直观的观察到双脱落空泡的吸引现象.另外,从图5中可以观察到上、下脱落空泡的径向尺寸有较大差异,并且后发的上方射弹脱落空泡尺寸较大,说明在形成入水空泡过程中后发的上方射弹空泡含气量更多.双脱落空泡在发展一段时间后逐渐混合.
图5 B型射弹在不同间距下的双空泡吸引现象
图6为A型射弹在不同间距下的双空泡流场结构.同B型射弹相比,A型射弹入水后产生的空泡尺寸较小,从图6中同样看到双空泡之间的吸引现象,具体特征和B型射弹稍有差异,在较小间距下下方射弹的脱落空泡部分向上发展,部分沿下游发展,此分流现象的形成机理将在下文说明.当间距增大至2.0D0时下方射弹的脱落空泡呈单股发展,脱落空泡朝着互相靠近的方向偏移.同图5相同先入水射弹的脱落空泡偏移量更大,从图6中可以发现,双空泡闭合位置附近的脱落空泡处于不同的流场环境中.其中空泡耦合段(矩形区域标识)受后发射弹空泡演化影响具有一定速度,其压力较低,脱落空泡亦对流场产生一定影响.另外,射弹受到的静水压力存在差距(上方射弹上侧环境压力较小),进一步揭示了空泡耦合段的压力更低,在较大的压差驱动下先发射弹的脱落空泡出现较大偏移.
图6 A型射弹在不同间距下的双空泡吸引现象
为研究空泡吸引现象随射弹间距的变化规律,定义量纲一的距离d*如下:
式中:D为射弹轴间距;d为后发射弹闭合点与另一股脱落空泡的距离(如图6中所示);D0为射弹直径.通过提取不同间距工况下的相同时刻的距离d,得到以下两种射弹的空泡靠拢幅度随着射弹间距的变化特点.从图7中发现A型射弹的曲线呈现先增大后减小的特点,而B型射弹随着间距增大,脱落空泡偏移量逐渐增大. 从而可以得知,对于A型射弹,其脱落空泡相互靠拢,并且靠拢的幅值随着间距增大而增大(在D=2.0D0~3.0D0内),当到达某个临界后,随着间距增大,脱落空泡之间的靠拢趋势减弱;对于B型射弹,其入水空泡尺寸较大,在相同间距下,空泡耦合范围较大,双射弹之间的流体区域较小,如图7中该间距下某时刻的空泡演化状态所示.因此当A型射弹达到临界间距时,B型射弹尚未达到临界距离,表现为曲线随着射弹间距增大而继续增大.
图7 双空泡吸引吸现象随射弹间距变化规律
由上述分析可知双空泡之间的吸引现象不仅和射弹间距有关,还与空泡尺寸有关. 通过试验结果对双空泡的吸引现象进行初步分析:当双射弹间距较大时,射弹入水瞬间其头部产生完整空泡.在射弹之间存在有限的流体区域,射弹运动时排开周围液体,同时空泡的演化对流场造成扰动,使得中间流体产生一定的扰动速度.根据伯努利原理知此区域压力降低.在压差力驱动下射弹的脱落气团向低压区发展,形成图7中的吸引现象,并且脱落空泡的靠拢幅度随着间距的增大而增大,当达到临界间距时,靠拢幅度随着间距增大而减小.在后续研究中,将对此吸引机理进行深入分析.
图8为工况7中的并联入水过程空泡演化过程.从图8中可以看到下方射弹率先入水,与上方射弹前后间距l=0.614L0.此工况下出现了射弹空泡被截断的现象.在2t0时刻射弹空泡发展尚未受到上方射弹的影响,其形态基本和单发射弹相同.当运动至6t0时,从图8中看到双发射弹的空泡形态存在很大的差异:上方的后发射弹空泡基本呈椭圆形并包裹着射弹(图中实线标识),下方先发射弹空泡未包裹射弹,并且其上、下方的轮廓不对称(图中虚线标识).造成空泡不对称的原因解释如下:在相同位置截面处,上方射弹空泡扩张速度远大于下方射弹空泡扩张速度,这将对下方射弹空泡发展造成抑制作用.由于双射弹间距过近,这种抑制效应将下方射弹空泡截断.如图8中6t0所示,下方射弹的截断位置在上方射弹头部正下方附近,说明了下方射弹空泡被截断缘于上方射弹空泡扩张的“挤压”.当下方射弹空泡上侧空泡被截断后,下侧空泡受到其影响后亦闭合于弹身某处.同上侧空泡相比下侧空泡未受到上方射弹作用,其闭合位置偏向尾部,形成了不对称的空泡轮廓.运动至10t0时,下方射弹空泡由于尾部闭合处产生的回射流向前运动,整个射弹空泡变得不透明,空泡界面产生波动.此时上方射弹的空泡较透明,与之前时刻相比,其尺寸进一步增大.射弹进一步运动,在16t0时可以清楚观察到双射弹尾迹相互吸引的趋势.其中上方射弹脱落气团在运动时逐渐偏下,下方射弹的脱落气团偏向上方,双脱落空泡在发展过程中逐渐混合(图中椭圆所示).
图9为工况8的并联入水过程空泡演化过程.从图9中可看到下方射弹率先入水,上、下射弹的前后间距l=0.25L0.该工况下发现弱化的截断现象,表现为凹型空泡轮廓.在3t0时刻,上方射弹空泡基本对称而下方射弹空泡出现较大变形(图中红线标识):其上侧空泡在上方射弹对应位置处的空泡尺寸普遍较小,其原因是下方射弹空泡扩张时受到上方射弹的抑制作用.与图8相比射弹的间距较大,并且前后距离较近时相同截面位置处的空泡膨胀速率差距较小,因此形成空泡上侧的凹型结构,未产生图8中下方射弹空泡被截断的现象.在5t0时刻,上方射弹空泡基本透明,从图9中可以明显观察到其下侧尺寸大于上侧,最大尺寸位置距离下方射弹尾端较远处.同时下方射弹空泡在空泡尾部出现分流现象:部分向上方发展,小部分基本沿水平向下游发展.在8t0时刻观察到下方射弹闭合点处的双泄气通道.射弹继续运动,双射弹空泡均闭合后相对独立发展.
图8 双空泡演化过程中的截断现象
图9 双空泡演化过程中的弱截断现象
在工况9(图6中的Δ=1.5D0)中双空泡演化规律与图9类似,均出现了弱截断和分流现象.分析此现象的形成机理:图10中的Ⅰ表示双空泡发展时的强耦合区域,入水后在Ⅰ区域内射弹存在前后间距且速度差距很小,因此在相同位置截面处上方射弹空泡膨胀速度大于下方射弹,根据独立膨胀原理可知3t0时刻在Ⅰ中下方射弹空泡部分扩张部分处于收缩阶段,上方射弹处于膨胀阶段.当双空泡均处于扩张阶段时,对应的区域的空泡均受到抑制,空泡尺寸较正常演化时小;当下方射弹收缩而上方射弹扩张时,下方射弹的收缩过程加快.这种耦合状态导致下方射弹空泡上侧出现未充分发展后的凹状结构.相对于截断现象,双射弹前后间距l较小,上方射弹对下方射弹空泡的抑制作用较弱,不足以将下方射弹空泡在弹身处截断,但是下方射弹空泡收缩加快,造成其空泡壁面触碰到射弹表面.在下方射弹发生尾拍后,其尾部产生大量分散的脱落气团.在Ⅱ区域,下方射弹空泡基本闭合,由伯努利原理可知此区域流场压力较低,因此上方射弹空泡收缩过程相对减缓,表现为此区域内上方射弹空泡下侧尺寸较大.同时较快的流速导致压力较低,加快了脱落气团向上方发展的趋势,形成图中出现的分流现象.
图10 双空泡耦合发展时的流场结构示意
1)建立了双射弹并列、连续发射试验系统并通过试验验证了该系统的有效性,为进一步开展并联入水多空泡融合与演化特性,以及多射弹运动干扰等奠定基础.
2)高速射弹并联入水过程中,其空泡演化过程存在穿越、吸引及截断3种典型现象.空化器直径和射弹间距是影响双空泡耦合演化的重要因素.相同射弹轴向间距时,随着射弹空化器直径减小,双空泡之间的耦合程度减弱,完全穿越转变成部分穿越现象;相同的空化器直径时,随着双射弹的轴间距增大,射弹的脱落空泡之间的吸引现象逐渐减弱,但是脱落空泡的偏移幅值随着轴间距的增大呈现先增大后减小的变化趋势.
3)当双射弹存在前后距离时,先发射弹易出现空泡的截断现象.随着轴向间距增大,截断现象将减弱.根据独立膨胀原理可知,相同条件下射弹前后间距减小时截断现象亦将减弱.