摘要：潜射导弹/射弹超空泡出水时超空泡与界面的作用机理目前还不完全清楚，而非定常流重力作用下的超空泡形状的描述也需要进一步探讨。为此，进行了水下高速射弹超空泡出水的实验研究，用高速摄影仪拍摄了航行体出水过程中空泡形状变化的图像，测量分析了各种不同头型航行体诱导的空泡尺寸、轮廓和生成溃灭形态的变化规律，并将研究结果与已有的理论公式和半经验公式进行了对比。研究表明，对所研究的工况，Savchenko半经验公式的适用范围可以扩展到空化数小于0.12；此外，实验测得的超空泡归一化长度亦高于考虑弗劳德数影响的Vasin的理论值，文中采用一系数对该公式进行了修正。实验研究还发现，航行体头型对超空泡的尺寸影响较大，超空泡与自由液面接触后，空泡直径有一定程度的增大。
关键词：超空泡；出水；高速射弹；超空泡形状
中图分类号：O359
文献标志码：A
文章编号：0253-987X（2015）03-0067-07
当物体在水下高速运动时，物体周围的液体压力降低，当物体速度达到一定程度时，液体会发生气化，在物体周围形成空泡，当空泡包裹住整个物体表面时，这种现象被称为超空泡。水中运动物体诱导的超空泡流动是一个典型的气一液两相流问题。利用超空泡技术可以大大降低水下运动物体的黏性阻力，实现减阻提速的效果，在水中兵器、水下运输、海上急救等方面有着广泛的应用。
水下航行体超空泡出水问题是超空泡流动的一个重要分支，航行体的运动方向为垂直亦或是倾斜向上，主要应用在水下潜艇或者水下潜射导弹上。国外的Waugh、Xing-Kaeding等对航行体带空泡出水过程中自由液面和空泡的溃灭以及重力的影响规律进行了实验研究；国内的鲁传敬教授利用CFD方法以及动网格技术对细长体出水过程中航行体受力波动、液面的升高、航行体尾部拖水等问题进行了研究，王一伟等对典型航行体模型水下发射的全过程进行数值模拟，得到了航行体周围空泡演化过程、表面压力分布变化、航行体关键截面的全时程载荷响应。
空泡的形态及其形成机理和影响因素研究一直是该领域的研究重点，因为空泡的几何形态直接关系到水下航行体的受力情况，对航行体的稳定性有着重要的影响。国内外学者开展了大量的研究工作，Munzer-Reichardt、Logvinovich、Savchenko等提出了描述空泡形状的理论公式和半经验公式，Parishev和Vasin提出了考虑重力效应时的空泡形状理论公式；曹伟、周素云等利用实验手段对水平自然射弹空泡的形态进行了研究，张学伟、刘海曼等应用数值仿真方法研究了空泡的形态。自然空化空泡在出水过程中，不可避免地会受到重力和自由面的影响，因此其形态变化更为复杂，目前国内外相关方面的研究还比较少。本文利用高速射弹模型实验和高速摄影方法，研究了不同模型和参数下自然空化超空泡出水形态尺寸的变化规律。
1实验装置
实验的研究工作是在一套高速航行体出水实验系统中进行的，实验装置如图1所示。该系统主要由发射装置、观察装置和测量装置3部分组成。其中，测量装置包括工控计算机、高速摄影仪和照明灯，发射装置由高压气源、高压气缸、电磁阀和发射管等构成。发射装置为一台一级轻气炮，在实验时，通过调节高压气缸内气体压力驱动发射管中航行体高速运动。本实验采用高压氮气驱动，压力范围为0～2.5MPa，可以使航行体加速到100m/s以上。观察装置包括航行体捕获器、观察窗、法兰、球阀、防水隔膜和法兰盖板。航行体从发射管加速射出后，垂直向上地冲破防水隔膜进入水箱，在水中形成超空泡；此时，采用高速摄影仪（美国Cooke公司的pco.1200s）记录整个过程，摄影所需的光源由3盏1kW的照明灯提供。实验装置的具体构造和实验步骤，这里不再赘述。
本研究对多种模型的航行体进行了射弹实验，此处只给出几种。实验采用的航行体模型的外形及其尺寸如图2所示，表1给出了相关参数，模型的材料为镁铝合金。在本文中将采用1#～4#来分别代表所研究的不同航行体模型。
进行后期的实验数据处理时，需要确定拍摄图片尺寸与实际尺寸之间的缩放比例，本文采取的方法是在实验前拍摄直尺图像（见图3），根据直尺上的刻度取得图像的缩放比例。
由于空泡、水、有机玻璃观察窗的折射会导致图像各部分的比例不均匀，相关研究表明这种误差的影响较小，一般在0.3%～0.6%，可不予考虑。
空化数是空泡流动中一个很重要的参数，其表达式如下
2实验结果和分析
2.1自然超空泡从生成到溃灭的形态变化
图4示出了4#射弹模型（空化器直径为4.5mm的圆台型航行体）的出水全过程。本文实验的拍摄速度为4000帧/s，即相邻两张照片的时间间隔为0.25ms。航行体经水箱下方打开的球阀进入水箱，从下向上运动。当航行体进入可视观察窗时，航行体的速度为63m/s。因为航行体的速度较高，在航行体周围产生了较大尺寸的超空泡。在图4（16）中可看到完整的超空泡。航行体在阻力和自身重力的作用下（其中阻力起主导作用）速度逐渐减小，对应的超空泡的尺寸也有所减小，这在后续给出的结果中也可看出。
该工况因为航行体入水速度较大，因而一直到出水时刻（见图4（23））超空泡形状都较稳定，并未出现自尾部至头部的空泡溃灭和再生成，当航行体出水时才发生空泡的溃灭。
当航行体从下向上运动时，包裹航行体的超空泡因为航行体的运动，超空泡所在的水位也在不断地减小，即使在同一时刻，超空泡自身在长度方向上受到的水压也是不同的。因此，为了研究出水超空泡的形状变化特性，选取了3个水深进行研究。以水面自由面为O点，向上为正方向，水深h分别取-265.43、-153.77、-81.89mm。
定义航行体头部到达固定水位的时刻为0，用航行体特征直径（Dn=6mm，见图2）对空泡直径Dc进行归一化。航行体以不同速度穿过选取的不同水位时，归一化直径随时间的变化关系如图5所示。由图可知：航行体穿过固定水位时，空泡经历了快速扩张、缓慢变化、快速收缩的过程；空泡生长到溃灭的时间为3.5～4.25ms。-265.43、-153.77mm水位处的空泡归一化直径都近似地呈抛物线形式，在t=2ms时达到峰值；h=-81.89mm水位处，空泡直径并没有像其他水位那样在2ms时达到峰值，而是接着扩张到2.5ms时才达到峰值。其原因可能是：在其他两个水位处，空泡从形成到溃灭，航行体都在水下，尚未越过自由水面，而对于h=-81.89mm的水位，图5中横坐标的2ms时刻正好对应图4（23）所示的航行体出水以及超空泡开始接触空气的时刻，在自由液面和浮力等的影响下，空泡头部的纵向体积增长被抑制，空泡尾部横向地扩张，从而导致在此水位空泡直径的继续膨胀增长，延缓了溃灭。 2.2出水自然空化超空泡的尺寸变化
超空泡的形状会影响航行体在水中所受到的阻力。超空泡的长度和直径采用航行体特征直径Dn进行归一化。
图6给出了4种航行体的归一化长度随空化数的变化与理论和半经验公式的对比。据式（1）计算空化数σ时，式中的p∞取航行体头部所在位置处的水静压强和当地大气压之和。须注意的是，此处并未考虑Fr的影响（Fr的影响随后讨论）。图中带符号的4条实线为本文的实际测量数据；实线为根据Logvinovich所提出的空泡截面独立膨胀原理计算得到的结果；长虚线为根据Savchenko提出的半经验公式计算得到的结果，它适用的空化数范围在1×10
-4～1.2×10
-2内；短虚线（几乎与实线相重合）是根据Munzer-Reichardt模型计算得到的结果。该模型是赖卡特基于低阶势流理论方程导出的，他所给出的空泡形状如下
从图6中可看到，所有的实验数据都要高于Logvinovich和Reichardt的理论公式值，而低于Savchenko的半经验值。
2#和4#航行体与Savchenko所给出的数值非常接近。1#航行体的实验数据介于Logvinovich和Savchenko的值之间，变化趋势与理论和半经验公式相似。
从图6中还可看到，对于上述3种工况（即1#、2#和4#），随空化数的增大，空泡的归一化尺寸都逐渐减小，且变化趋势基本一致。对比1#和2#的结果发现，在相同的空化数时，2#航行体的归一化空泡长度要大于1#的值（这两种工况进行归一化的空化器直径是相同的）。3#航行体的超空泡归一化长度随空化数的变化趋势较为平缓，这可能是因为：一方面，理论公式只是用来描述圆盘空化器航行体产生的空泡，而3#为60°锥角航行体，它的头部绕流情况和平头空化器头部绕流相差很大，产生的空泡要小一些；另一方面，60°锥头时超空泡的尾部溃灭影响范围较大，这使得超空泡的尺寸难于测量。
图7为实验测得的几种工况超空泡归一化直径与理论和半经验数据的对比。
从图7中可看到，4#航行体的实验数据与Savchenko给出的半经验公式的数值符合很好，但对于其他3种工况，则或高于半经验公式给出的较大值，或低于理论公式给出的较低值。其原因一方面可能是因为本文所实验的航行体的最大直径不大，只有6mm，而在实验时存在着各种测量上的误差，如实验时空泡、水、有机玻璃观察窗的折射，以及实验后超空泡尺寸测量时的误差。除此之外，当超空泡尺寸相对较小时，超空泡的尾部溃灭也会影响超空泡的测量结果。
对于重力场中的超空泡流动，当考虑重力效应时，需引入Fr。Vasin给出了考虑Fr时超空泡长度的计算公式
式中：a=2；k为本文引入的一修正系数，当k=1时，即为Vasin所给出的公式。
图8给出了实验数据和Vasin公式计算结果的对比。从图中可看到，对所实验的4种航行体，超空泡的归一化长度都要高于Vasin所给出的计算值。
对比图8与图6可以发现，图8中的实线与Logvinovich和Reichardt的理论值接近。因此，本文在式（6）中引入一系数k=2对Vasin表达式进行修正。图8中的虚线即为本文得到的k=2时的经验表达式，用上限公式标注。对于垂直出水超空泡，当考虑重力效应时，超空泡归一化长度可能介于k=1和k=2之间。
2.3出水自然空化超空泡的轮廓分析
如图9所示，以航行体头部为坐标原点，空泡轴线为x轴，x处的空泡直径为D（x），分别测量不同x位置的空泡直径，即可得到在某空化数下的空泡轮廓。
为方便研究比较，对空泡轮廓参数同样采用航行体特征直径Dn进行归一化处理，具体定义如下
x=x/Dn；D=D（x）/Dn
此处，只给出了一个工况下（平头长径比为12时）3个空化数实验得到的结果，3个空化数分别为0.147、0.155、0.163。图10给出了实际测量的3个空化数下的空泡图像，从该图中并不能直观地看出不同空化数时空泡轮廓曲线的不同。
图11给出了3个空化数的空泡外形轮廓曲线的对比图。从图中可以看出：空泡的相对直径随轴向相对坐标x先增大后减小，呈抛物线变化；相同x下空化数越小，相对直径D越大；x8时，各空化数下D的差距变大。相同模型的航行体在不同空化数时，空泡前部分轮廓曲线大致相同，后部分轮廓随空化数增大而减小，空化数的变化对空泡尾部影响较大。
图12给出了空化数为0.155时实验测得的空泡外形轮廓曲线与Reichardt和Logvinovich公式的理论值对比。Reichardt和Logvinovich的理论曲线为对称的椭圆弧形，两种理论曲线峰值相同，但Reichardt曲线的其余部分低于Logvinovich的理论曲线。
在该实验工况下，实验数据和理论公式值吻合较好，各空化数下空泡轮廓的前半部分实验数值比Logvinovich的理论数值小，但比Reichardt的理论数值大，空泡轮廓的后半部分实验数值更接近Reichardt的理论数值。实验测得的空泡曲线不对称，空泡后半部分的直径要小于前半部分对应位置处的直径，其原因可能是垂直出水的空泡受到重力的作用，在空泡表面上的纵向压力梯度使得尾部空泡收缩所致。
3结论
（1）在高速摄影照片中设定欧拉控制面，让空泡通过控制面，得知了空泡直径经历了快速扩张、缓慢变化、快速收缩的过程；自由液面对临出水空泡的径向收缩有抑制作用，并促使空泡直径有一定程度的增大，空泡溃灭时间增加。
（2）对本文所给出的工况，Savchenko所给出的半经验公式的适用范围可以扩展到空化数为0.12的情况，此时实验数据与半经验值符合较好。实验测得的超空泡归一化长度值都高于Logvinovich和Reichardt的理论计算值。相同空化数下，航行体头型对空泡的归一化尺寸影响较大。减小航行体头部锥角会大幅度地降低空泡尺寸。
（3）对所实验的几种工况，当引入考虑重力效应的弗劳德数时，实验测得的超空泡归一化长度基本上在Vasin的理论计算值的1倍到2倍之间。
（4）出水超空泡的轮廓在最大空泡直径处的两侧不完全对称，后半部分尺寸略小。这可能是因为在重力的作用下，垂直出水空泡表面的纵向压力梯度使得尾部空泡收缩所致。
