水下高速航行体超空泡减阻的大涡模拟与实验研究
水下兵器发展过程中，提高水下射弹、鱼雷的航行速度一直是摆在设计研究人员面前的最大障碍。过去研究人员采用优化航行体结构、添加航行体表面涂层、增加和改进动力系统推进方式等常规方法来提高水下航行体速度，研究结果表明由于水下阻力极大，这些方法很难大幅度提高水下航行体的运动速度。随着俄罗斯超空泡高速鱼雷的问世，打破了以往常规方法的限制，成功实现了水下航行体速度的提升，并且获得了非常喜人的效果。目前随着世界各国对于海上军事主动权的争夺日趋紧张，超空泡技术极具军事应用发展前景的特点，深深地吸引了各大军事强国的关注，开展了大量相关研究，并且取得了一定的研究成果。但是空泡流动是一种考虑相交、粘性、湍流运动、界面可压缩性的极其复杂多相流动现象，使得空泡流动的研究发展在很大程度上受到制约，其相关实验研究、数值预测和空泡机理研究仍然是目前多相流领域内的研究难点与热点。
针对这种情况，本文利用自行开发设计的大涡模拟算法程序，进行了水下高速航行体超空泡减阻问题的空泡流动数值模拟;利用空泡数值模拟理论结果设计航行体模型，研制水下测速仪器，开展水下超空泡航行体实验研究，具体研究工作内容如下:
首先，基于国内外关于超空泡技术的研究现状和发展动态，结合本文主要研究工作的实际工程应用背景，提出本文研究所要解决问题和解决这些问题所采用的研究方法与技术途径。
其次，通过对国内外关于空泡流动的数值模拟方法的分析研究，结合空泡流动特点，在均相流动理论的框架下，建立基于三维Navier-Stokes方程且考虑气液间相变的大涡模拟求解计算方法。详细分析了大涡模拟操作中常用滤波函数的特点和不同形式亚网格应力模型的优缺点，并且对本文采用的亚网格湍动能k-△模型进行了推导。讨论了可以处理气液间相变问题的不同空化模型。在此基础上，详细介绍了本文采用的数值离散方法，其中包括大涡模拟控制方程离散、压力-速度-密度的耦合算法以及PISO算法中压力修正方程的推导过程。另外详细分析了计算域和数值计算边界条件，讨论了不同网格类型的优缺点，并且针对复杂多相流场计算中所需计算量很大的问题，设计了基于MPI方法的并行算法。
第三，为验证本文计算方法的适用性和探究低速航行体空泡流动特性，开展了绕低速NACA0015水翼云空化流动的数值模拟研究。计算了空泡形成过程与周期性脱落频率，将计算结果与实验结果对比，两者吻合较好，证明了本文计算方法的适用性。在此基础上，分析了包括初始涡、回射涡、空泡分离以及空泡再生的非定常周期性脱落过程，研究了Kunz空化模型和Sauer空化模型对空泡形态的影响以及分析了不同空化数对水动力特性的影响。
第四，为验证本文计算方法的正确性和探究高速航行体空泡流动特性，分别从空泡形态特性与阻力特性两方面，开展了绕不同外形结构三维航行体模型超空泡流动的数值模拟研究。通过计算得到超空泡由航行体头部初生至完全包裹航行体的非定常形成过程。将计算的超空泡无量纲特征长度、厚度与Logvinovich实验结果对比，两者吻合较好，证明了计算结果的正确性。在此基础上，详细分析了不同空化数、空化器直径、空化器形状、航行体攻角以及航行体尾翼形状对超空泡形态特性和阻力特性的影响。
最后，为进一步验证本文数值方法的正确性，基于高速航行体空泡形态特性与阻力特性的数值分析结果，结合圆盘空化器具有较高减阻效率和良好空化能力的特点，设计了带尾翼结构的圆盘空化器多种外形结构航行体模型，开展了高速航行体从空中斜射入水中形成超空泡的实验研究。通过水下高速摄像获得了小空化数下自然超空泡形成图像。讨论了航行体空中运动稳定性与航行体能否形成超空泡之间的关系。利用水下测速靶上的航行体穿孔位置，分析超空泡航行体的水下弹道稳定性。获得了不同空化数和空化器直径下超空泡航行体的平均阻力系数，并且分析了空化数和空化器直径对平均阻力系数的影响。将实验结果与数值计算结果对比，两者吻合较好，验证了数值计算方法的正确性。在此基础上，计算了水下超空泡航行体的减阻率，分析了超空泡减阻效果。
综合上述关于高速航行体超空泡减阻问题的数值模拟与实验研究可知，本文实现的数值计算方法能够较好地模拟水下高速航行体超空泡流动;在空化数σ=5.9×10-4、空化器相对直径(D)n=2.3×10-1条件下，水下超空泡航行体可获得最高约为94.67％的减阻率。
大涡模拟;PISO算法;超空泡减阻;水下高速航行体;数值模拟
南京理工大学
博士
流体力学
谭俊杰;易文俊
2013
中文
U661.1
144
2013-12-31（万方平台首次上网日期，不代表论文的发表时间）