毕业论文网
关于免费研究论文范文在这里免费下载与阅读,为您的研究相关论文写作提供资料。
绕过带凹槽的圆盘空化器的超空泡流
摘 要:采用VOF的方法,对带有不同凹槽数空化器的航行体诱导产生的三维超空泡,进行了数值模拟.分析了圆盘形空化器凹槽数对所形成的三维超空泡形状的影响;将数值模拟的结果与Loinovich超空泡截面独立扩张原理的结果进行对比,发现数值计算结果与Loinovich半经验公式计算结果具有较好的一致性;得出航行体的阻力系数随时间变化曲线.研究发现:随着空化器凹槽数的增加,航行体阻力系数以及超空泡的无量纲直径和长度均呈现出减小的趋势;带凹槽空化器产生的超空泡是典型的三维复合型空泡,其特征是:在来流方向上,不同位置处超空泡横截面的形状是不同的;但是超空泡具有很好的自修复特性,能够在较短时间内恢复成光滑椭球形.该研究的计算数据以及结论可为航行体空化器的设计提供参考.
关键词:不同凹槽圆盘空化器;超空泡截面独立扩张原理;VOF方法;阻力系数;自修复特性
中图分类号:O３５２;O３５９ 文献标志码:A 文章编号:１６７３Ｇ３８５１(２０１８)０５Ｇ０３１２Ｇ０７
０ 引 言
航行体在水下高速运动时,头部空化器能够诱导产生超空泡,包裹整个航行体,大幅度降低航行体的阻力.但水下航行体超空泡受到扰动流场时,容易溃灭,极其不稳定.Schenko等[１]提出一种带有浅凹槽的、多边形的空化器,该空化器能够产生稳定力矩,且具有抗扰动的作用.在此基础上,施红辉等[２]研究了水下航行体表面开槽数量对超空泡流场影响,发现随着航行体表面的凹槽数的增加,超空泡的最大无量纲直径和长度会减小,但无法确定航行体表面开槽对无量纲尺寸的影响,也未对开槽航行体超空泡水动力学特性进行分析等.所以针对超空泡航行体的主要沾湿部位———头部空化器,进行了深入的研究,能够达到超空泡的最大减阻效果.超空泡减阻增速的优势,已经被美俄等国用于高速水中兵器的设计中,但我国水下超空泡流动的研究仍亟待提高[３Ｇ６],只有夯实超空泡流动技术,才能进入空(气)中到水下连续机动的超空泡兵器开发阶段.
超空泡流动是通过空化器的脱体绕流发生的,因此空化器的几何形状对超空泡流动有着十分重要的影响[７].先前研究主要针对圆盘、圆锥、矩形板、卵巢型等空化器形状的空化器,其中圆盘型空化器是应用最为广泛[８].针对以上空化器形状,国内外许多专家学者也进行了许多实验、理论和数值模拟方面的研究.如罗格维诺维奇等[９]在２０世纪６０年代通过理论和实验研究给出了圆盘空化器阻力系数、空化器直径和空化数与自然超空泡形态的函数关系.Schenko[１０]在实验基础上研究了锥形空化器,通过外推法给出了空化数趋近于零时的阻力系数,并与近似公式进行了比较.栗夫园等[１１]在水洞试验中研究了锥形空化器,分析了锥形空化器的锥角对航行体流体动力学特性的影响,对圆盘空化器超空泡流体力学特性有待进一步分析.徐清沐等[１２]基于雷诺平均NＧS方程和均质平衡流理论对带凹槽和凸起的圆盘空化器进行了数值模拟,并给出了空化数对空泡凹凸形态和阻力系数的影响.金大桥等[１３]分析了水下航行体结构参数对自然超空泡减阻特性的影响,获得了超空泡航行体的阻力系数变化趋势.牟斌等[１４]利用空化流动隐式求解方法对三维尖锥空化器形成的自然空泡的流场特征进行研究.
为了进一步研究超空泡航行体最大限度的减阻机理,本文选取不同凹槽数的圆盘空化器航行体为研究对象,应用商业软件Fluent对其所形成的三维超空泡进行数值模拟,验证三维超空泡的自恢复特性,还分析了空化器凹槽数对航行体流体动力特性的影响.
１ 数值模拟方法
１．１ 控制方程
对超空泡射弹的数值模拟是基于连续性方程和动量方程,以及湍流模型采用标准kＧε 模型,多相流模型采用VOF模型,空化模型采用SchnerrＧSauer空化模型.
空化数表征超空泡的最为重要的无量纲参数之一,其定义为:
１．２ 计算域、网格划分及边界条件
数值模拟采用复杂表面形状带有圆盘形空化器的航行体,以Yang等[１７]的航行体为原型,又根据试验中高压的口径为参考设计,总长为４８mm,最大直径为６mm,空化器是直径５mm、厚度１mm的圆盘;在空化器边缘处开槽数分别为０、２、４、６,凹槽为１mm×１mm×１mm 的正方体,如图１所示;航行体模型的几何尺寸参数如图２所示.
计算域是一个４００mm×２００mm×２００mm 的长方形水箱,如图３所示.航行体头部空化器到前方来流的速度入口的距离为８０mm,位于计算域中间处.航行体的速度为１００m/s,压力值为一个标准大气压,四周壁面设为无滑移壁面边界条件,近壁面处采用标准壁面函数,见图３.通过Fluent中的Patch功能,可以设置计算域的初始相.
网格的划分如图４所示,网格全部采用六面体结构网格划分,为了能更清晰地显示出超空泡壁面处的汽液交界面,在航行体模型的四周进行了局部网格加密,最终得到的０、２、４、６道槽航行体计算域网格的节点数分别为７２７８１２、７０６３０７、７７４５７６、８９６７９２个.在数值计算中采用VOF的多相流方法捕获汽液交界面,并且采用标准的kＧε 湍流模型,与连续性方程、动量方程等形成封闭的方程组,进行数值模拟的求解.此外使用压力基求解器,流动为非定常流动,流动中设置空气、水、水蒸气三相,并考虑重力的影响.
通过改变航行体头部圆盘空化器上的凹槽数量,逐个进行数值模拟的计算,分析不同凹槽数空化器的航行体对超空泡流场的影响,并进行了水动力学分析.不同凹槽数空化器航行体网格划分情况示意图如图５所示.
１．３ 网格无关性验证
选取０道槽航行体的计算域进行网格划分,得到节点数分别为５０万、７０万和１００万左右三种密度的网格,分别进行数值模拟计算.将计算结果进行后处理分析分析,绘制不同网格密度下航行体阻力系数Cd 随时间变化的曲线图,如图６所示.
对图６分析可知,中等程度的网格质量(７０万左右)与较高程度的网格质量(１００万左右)的计算精度相差不大,最大偏差为４．９％,最小偏差为１．２％,故选取节点数为７０万左右的网格能够满足计算的要求.
２ 计算结果分析
２．１ 计算结果
头部空化器开０、２、４、６道凹槽的航行体在１００mm水深处诱导产生的超空泡数值模拟的三维和截面超空泡的水相图如图７—图１０所示.每幅图由两部分构成,分别是超空泡三维形态图和超空泡纵向截面图,每一幅表示的时刻T 依次为０．２、０．４、０．８、１．２、１．６、２．０、２．４ms.
从图７—图１０可知,当T＝０．２ms时,在航行体的空化器处、肩部以及尾部同时形成局部空泡;当T＝０．４ms时,空化器、肩部以及尾部空泡逐渐融为一体;当T＝０．８ms时,形成一个完全包裹航行体的超空泡.综合超空泡的截面水相图来看:在速度相同的情况下,不同凹槽数空化器航行体的超空泡的初生、发展过程比较接近.但是从空泡的三维形态图,可以看出,不同的凹槽数空化器航行体形成的超空泡会根据相应的凹槽数形成凸凹形状.然而,随着空泡长度的增加,凸凹形状就会很快被修复成光滑椭球形,超空泡所具有的这种特性被称作超空泡的自修复特性.
分别对０、２、４、６道槽空化器航行体形成超空泡进行定量分析,分别测量T＝２．４ms时刻的超空泡的最大直径及长度,并进行无量纲转化,即得最大无量纲直径Dc 分别为３．１２０、２．８２９、２．８１９、２．７０４;最大无量纲长度Lc 分别为３２．３７２、３０．２７８、３０．０３４、３０．１１８.通过对比可知,随着航行体圆盘空化器凹槽数的增加,其诱导产生的超空泡的最大无量纲直径和长度会呈减小的趋势.超空泡的最大无量纲直径Dc＝Dc/Dn,最大无量纲长度Lc＝Lc/Dn,其中Dc、Lc 是空泡的最大直径和长度,Dn 为航行体空化器的直径.
图１１给出不同凹槽航行体的距离圆盘空化器的距离为１mm 超空泡的横截面图;图１２给出了４、６道槽航行体的距离圆盘空化器的距离为１、３、５、７、９、１１mm 处超空泡的横截面图.为了节约篇幅,在最大空泡直径Dc 处(约在８０mm 处)的结果没有给出.
将图１１和图１２对比分析后可知:a)圆盘空化器有诱导产生超空泡的作用,且在超空泡的初生阶段,其截面形状会根据凹槽数的不同,形成不同凸凹形状的超空泡,但最终都变成圆形.b)图１２给出了４、６道槽航行体超空泡横截面从四角形、六角星形变化为圆形的过程.可知,在较短距离范围内,超空泡的横截面发生了较为显著的变化,体现出超空泡有很强的自修复能力,超空泡所具有的这种自动修复自身形态,我们称之为超空泡的自修复特性.这与超空泡所具有的自动缝合空洞的能力相似,均是由于超空泡处于高度非平衡状态,具有很大的熵差才驱动了这种自修复特性[１８].
图１３给出了超空泡横截面自修复成圆形时所处的无量纲位置x/Dn,其中横坐标代表不同的凹槽数.图中对比分析了文献[１,１９]中航行体超空泡的自修复特性结果,文献中航行体速度为４０m/s,在航行体表面开槽,开槽数同上.由图１３可知,尽管计算工况的速度不同,但都出现了超空泡的自修复现象.
２．２ 与半经验公式的比较
Loinovich通过动量定理结合大量的实验数据,对超空泡的外形提出如下半经验公式[６]:
分析图１４可知:上述两个时刻的超空泡前沿外形轮廓与Loinovich公式得出的轮廓曲线具有较好的一致性,这表明了数值模拟的准确性;而后沿是由于计算域限制,数值模拟并未达到超空泡的最大直径和长度,但从选取的两个时刻可以看出,数值模拟结果与经验公式结果是一个逐渐趋向的过程,在足够大的计算域中的某一时刻,数值模拟结果将会与半经验公式完全重合.
３ 航行体阻力系数分析
航行体超空泡在水下航行时,所受到的阻力F主要包括压差阻力Fp 和粘性阻力Ff 两部分,压差阻力主要取决于航行体头部与航行体尾部的流场压力,而粘性阻力与航行体表面积及其沾湿介质的密度有关[２０].综上所述,本文中考虑了压差阻力系数、粘性阻力系数以及总阻力系数.如下式:
其中:S 为特征面积,这里取圆柱段的横截面积.式(１４)—(１６)中的压差阻力、粘性阻力和总阻力可以在通过对数值模拟的后处理计算求得.
图１５—图１７ 是超空泡航行体的压差阻力系数、粘性阻力系数以及总阻力系数随时间的变化曲线.
根据图１５—图１７,并结合超空泡的形成、发展的水相图(图７—图１０)可知:a)在T＝０~０．６ms时间内,航行体的阻力系数急剧下降;这是由于超空泡逐渐形成,航行体的头部与尾部压差以及沾湿面积都急剧下降,致使航行体所受到的阻力急剧降低的缘故.b)在T＝０．６~２．４ms时间内,航行体的阻力系数趋于一条平滑的直线;这是因为在T＝０．６ms时间以后,超空泡完全包裹航行体,航行体的头部与尾部压差以及沾湿面积都保持不变,故其阻力系数均是直线.c)对单幅图来分析可知,对于不同的凹槽数航行体而言,随着凹槽数的增加,航行体的压差阻力、粘性阻力和总阻力系数均呈减小趋势.
４ 结 论
通过对不同凹槽数的圆盘空化器航行超空泡的三维流场的数值研究,将计算得到的数据进行后处理分析、与Loinovich半经验公式进行对比以及对其流体动力学进行分析,得出如下结论:
a)随着航行体空化器的凹槽数的增加,超空泡无量纲直径和长度呈现一定的减小趋势.
b)尽管不同凹槽数的空化器造成了非轴对称的三维超空泡形状,但在沿流向不长的距离处,超空泡能够很快自动修复形成光滑椭球形超空泡,这即是超空泡的自修复特性.
c)随着空化器凹槽数的增加,其受到的压差阻力和粘性阻力均呈减小趋势.
研究论文范文结:
关于对不知道怎么写研究论文范文课题研究的大学硕士、相关本科毕业论文研究论文开题报告范文和文献综述及职称论文的作为参考文献资料下载。