近自由面的超空泡流动的数值模拟,本文是关于模拟方面毕业论文提纲范文和空泡和数值模拟和数值模拟研究有关参考文献格式范文.
模拟论文参考文献：
摘 要:为研究物体在近自由面的运动规律,采用Fluent软件,选择VOF多相流模型和标准kＧε 湍流模型,针对物体在２６、１０６、２０６、４０６mm 四个不同水深下分别以６０、１００、２００m/s的初速度运动的超空泡流场进行数值模拟,研究水深对空泡形态的影响,计算出在自由液面上出现的先导波浪.结果表明:随着水深的增加,空泡的对称性越来越好,自由面上的波浪高度逐渐减小,空泡长度和空泡直径的变化与Loinovich空泡截面独立原理所导出的半经验公式拟合曲线更加吻合;随着物体速度的增加,自由面上的波浪高度也在增加.文中分析了先导波的起始距离以及水深对它的影响,可为高速射弹超空泡的形态特性和流体动力学特性研究提供参考.
关键词:空泡形态;先导波浪;自由面;波浪高度
中图分类号:O３５２;O３５９ 文献标志码:A 文章编号:１６７３Ｇ３８５１(２０１８)０５Ｇ０３０４Ｇ０８ ０ 引 言
超空泡现象在水中高速鱼雷、潜水艇发射战略战术导弹和反鱼雷鱼雷武器系统中有着重要的作用,水中高速运动物体在超空泡状态下运动时,可以减少阻力达９０％.高速鱼雷作为对舰攻击的杀手锏武器,暴露出其在空气和水中连续机动隐身的性能,这也必然涉及超空泡流动.另外,空泡的溃灭又会给机器带来损伤,如产生水泵和水轮机的气蚀问题.因此,详细研究超空泡的流动机理,具有重要的实际意义,且一直是国内外的热门研究课题之一[１Ｇ２].Chen等[３]通过求解RANS方程,模拟了超空泡在粘性浅水区域中的各种行为,包括靠近自由面处空泡和自由面的变形,并与实验数据进行了对比.超空泡产生时的空化数一般小于０．１,实现这个条件有三种途径[５]:a)保证航行体的速度大于４５m/s;b)降低空泡周围流场压力;c)采用人工通气的方法增加空泡内压力.水下航行体在高速运行时周围液体汽化,产生包裹航行体的空化泡,即自然超空泡,此时运动体在空泡内运行.在深水位时,由于空泡形态受空腔内外压差及自由面的干扰较小,所以空泡形态才可以看成一个近似的椭球体,但是在接近水面时,自由面的影响将不会被忽略.自由面的波动将会引起空泡表面压力的分布,进而影响航行体的运动轨迹.
近年来,在国防科技技术和海洋领域以及自然科学领域方面,对超空泡的研究一直在不断进行.施红辉等[４Ｇ８]设计了最高速度可达１００~１２０m/s的模拟高速射弹实验装置,开展高速弹体水平入水实验,同时也进行了相关数值模拟,获得水平超空泡流动中自然超空泡的形态特性和一些水动力学特性.施红辉等[６]采用高速照相机记录了在２８０、３４、１８mm三种不同水位下的弹体高速运动的过程,分析了水深对速度衰减的影响.施红辉等[７]通过大量实验对浅水区域中高速射弹进行研究,获得了近自由面处空泡形态的变化,并与半经验公式进行了对比.陈波等[８]通过采用Fluent软件模拟研究了浅水区时自由面对超空泡形状的影响,以及深水域时空化数对空泡形态的影响.马庆鹏等[９]对不同头型运动体高速入水的空泡形态进行了数值模拟研究,得到不同头型条件下高速入水运动参数及空泡形态发展规 律和流场的速度分布规律. 朱棒棒等[１０] 采用Fluent软件对不同头型物体高速入水过程进行了数值模拟,分析了入水空泡夹带空气的原因.黄海龙等[１１]通过采用结构化网格对圆盘空化器进行了三维数值模拟,解决了带攻角圆盘空化器形成超空泡的问题,并分析了带攻角圆盘空化器生成超空泡的一些特征.Truscott等[１２]对入水问题研究现状进行了综述,主要包括实验、理论和数值等方面的研究.
虽然目前国内外超空泡现象的研究已经有了许多研究成果,但是很少有文献考虑在浅水区自由面对超空泡流动及其诱导的波浪的影响.本文利用Fluent流体仿真软件,对物体在自由面附近的运动过程进行了数值模拟,将绘制的不同水深下的空泡直径的变化曲线与Loinovich和Schenko等[１３Ｇ１４]推导出的经验公式拟合的曲线作相关验证,进一步给出无量纲波浪高度随时间变化情况,通过对比不同工况下空泡形态及自由面的波动情况,分析自由面对超空泡流动的影响,给出先导波浪随时间的变化过程.模拟结果可为高速射弹超空泡的形态特性和流体动力学特性研究提供基础和参考.
１ 数值方法
１．１ 控制方程
本文为了解气、液、汽三相之间的相互作用,对近自由面超空泡流场进行了数值模拟.其中基本控制方程包括[８]:连续性方程和动量方程,选择VOF多相流模型和标准kＧε 湍流模型;设置气、液、汽三相;采用SchnerrＧSauer空化模型.
连续性方程为:
１．２ 数值模拟方法
对于二维流场来说,空泡形态可以用弹体轴线截面处空泡轮廓来描述,陈波等[８]采用的二维数值模拟方法,其模拟结果验证了二维流场的合理性,因此,本文也采用二维数值模拟方法对超空泡流场进行计算.本文给出了数值计算的物体模型,如图１所示.数值模拟中所使用的物体头型统一采用平头型,弹身为圆柱体[８].计算区域及其网格划分,如图２所示.计算区域所取的是１９００mm×９５０mm 的矩形区域,采用“H”型进行网格加密处理.如图２所示,计算区域的左端设置成空气和水速度入口,右端设置成空气和水压力出口,上端设置为symmetry,下端和弹体周边均为壁面.压力耦合和速度耦合均采用SIMPLE 算法,而梯度的求解则采用Leastsquarescellbased(最小二乘法),压力的求解则采用PRESTO！ 算法,至于动量离散则采用Secondorderupwind(二阶迎风格式).
２ 网格无关性验证
超空泡的形态特征的描述主要是通过空泡直径、空泡长度和空泡外形轮廓三个参数来呈现.图３给出了数值模拟结果中波浪高度H 和水深h(h 为射弹中心线到水面的距离)以及超空泡的长度Lc 和直径Dc 的定义.通过对参数的无量纲化定义,可以得到无量纲的空泡长度Lc为空泡长度Lc 与弹体直径Dn 的比值,无量纲的空泡直径(Dc为空泡直径Dc 与弹体直径Dn 的比值,无量纲的波浪高度(H为波浪高度H 与弹体直径Dn 的比值.
为了验证网格对数值模拟的影响,针对水深２６mm的工况,采用对自由面及弹体周围进行网格细化,得到相应的四个工况case１、case２、case３、case４,网格数分别是５２５８３２、５５０７９２、６００７１２、７００５５２.通过Fluent软件及后处理得到这四种工况的水相图,如图４(a)－(d)所示.
将数值模拟结果所得的水相图导入到AutoCAD软件中,经过多次测量获得波浪高度的数据,将这些数据导入Origin软件中,输出无量纲的波浪高度的曲线图,如图５所示.通过对不同网格数的数值模拟,从图５中可以得出,在相同的水深等条件下,不同网格数对波浪高度的影响;对比四种工况,case１和case２有相同的走势;case３和case４有相同的走势,从而得证网格无关性.但是考虑到实际计算周期等其他条件的限制,又要为了使计算精度较高,同时又能更好地捕捉自由液面波动,对周围和自由面附近做了网格加密处理,选取网格数为６００７１２万,时间步长为２×１０－６s.
３ 数值模拟结果与分析
３．１ 水深h 的影响
文中计算了四个水深(即h＝２６、１０６、２０６、４０６mm)下的超空泡形状及自由面随时间的变化过程.图６是四个水深下的水相图.如图６所示,深水处的超空泡的外形轮廓是一个上下基本对称的椭球形(见图６(c)－(d)).然而,近自由面的工况下,自然超空泡为上表面形状出现上凸状态,下表面略显平坦(见图６(a)).
图７(a)－(d)给出了在四个水深下的,空泡上下壁面的流线;图８给出了了在四个水深下的空泡外形轮廓.比较图７和图８可知,是空泡形状引起了自由面的抬升,Wang等[１６]也提出了类似的观点.当然,本文计算也指出,当水深大于２０６mm 后,空泡轮廓对自由面的抬起作用大大减小.
从图６(a)可知,数值模拟算出了超空泡在自由面上拉出了一个倾斜的先导波浪,波浪的前端位置超过了物体或超空泡的前端位置,这是由于快速运动的物体表面形成的超空泡对浅水层挤压造成的.施红辉等[７]已经通过实验证实这一结论.另外,从图６(a)中的计算结果可知,空泡右上部分的水层已经很薄,预示在实际情况中空泡很可能破裂,而空气进入超空泡中并扩张其尺寸.这一点在实验中也得到了证实[４].在图３中给出了自由面上的起波位置的定义,即波前距离S.图９(a)－(c)给出了S 的计算结果,它们大致呈现这样的趋势:随着水深的增加,自由面上的起波位置距离物体头部越远;随着速度的增加,自由面上的起波位置距离物体头部越近.
自由面的变化过程用自由面波动引起的最大高度,即波浪高度H (见图３)来描述.图１０给出了波浪高度随时间的变化,由图中可知:在深水位时,自由面的波动不是太明显,而在浅水区时,自由面的波动很大.针对h＝２６、１０６、２０６、４０６mm 四个水深,分别模拟了初速度为１００m/s和２００m/s的情况,进一步分析速度对自由面的影响,给出自由面引起的波浪高度的曲线图,如图１１(a)－(d)所示.在图１１中可以得出:在相同水深的条件下,随着速度的增加,自由面的波动较大.所以,通过模拟可以得出:自由面的变化不仅要受到水深的影响,而且还要受到物体速度的影响.由于水的声速为１５００m/s,当物体运动速度为２００m/s时,马赫数只有０．１３,属于不可压流动.又因为本文采用的理想流体,所以这三种介质的压缩性不用考虑.
由图１２可以得出,在浅水位h＝２６mm 时,计算出最大空泡直径要明显大于Loinovich 和Schenko的半经验公式的值.这是由于当水平运动超空泡经过浅水区时,自由面将向上抬起,空泡形状上下不对称;而且上部薄水层受空泡拉动,不再满足Loinovich的空泡独立膨胀原理.然而从数值量级上,计算结果与半经验公式的是一致的.
４ 结 论
本文选择平头射弹模型,计算了水深对超空泡形状的影响.弹体在h＝２６、１０６、２０６、４０６mm 四种不同水深下高速运行时,随着水深的增加,空泡对称性越来越好.超空泡运动诱发自由面上产生波浪,其最大高度随着水深的增加而减小;随着物体速度的增加而增加.数值模拟还计算出了在自由面上的波浪,是一个先导波浪,即波浪的前端位置超过了物体或超空泡的前端位置,这与施红辉等[７]的实验结果相吻合.
本文首次分析了先导波的起始距离以及水深对它的影响.在浅水区时,自由面的起波位置靠近空泡前端,速度越大越靠近;在深水时,自由面的起波位置远离空泡前端,趋于无穷远处,且空泡形状更接近一个椭球体,可为高速射弹超空泡的形态特性和流体动力学特性研究提供参考.
本文点评：此文是一篇适合空泡和数值模拟和数值模拟研究论文写作的大学硕士及关于模拟本科毕业论文，相关模拟开题报告范文和学术职称论文参考文献.
周恩来和邓颖超：那个年代的爱情
作为新中国的开路人,周恩来的一生精彩壮阔 他也有柔情的一面,他和邓颖超半个多世纪的爱情传奇,至今还为人们称道 现今社会,无论50后、60后,还是70后、80后,爱情和婚姻生活常面临一些普遍问题缺少话题.
艾叶泡脚的养生功效
艾叶泡脚使用方便,成为很多人养生的最佳选择 据百科名医网介绍,用艾叶水泡脚能有效地祛虚火、实火,可以治疗口腔溃疡、咽喉肿痛、牙周炎、牙龈炎、中耳炎……头面部反复发作的与虚火、实火有关的疾病,除此之外,.
植物弦歌（25）自由行走的花
文／管弦和蒲公英的相知相守,是在一寸一寸飞翔的时光中 记忆中的蒲公英是安静的,常常在绿草丛中静静地开 她的根牢牢地钉在地里,碧绿的长圆形的边缘类似羽状的叶子舒展着 冬末春初抽花茎,她的顶端会生出花序,.
老军医吴孟超：妙手仁心铸医魂
他是中国肝脏外科之父,主刀完成了我国第一例成功的肝脏手术；他翻译了第一部中文版的肝脏外科入门专著；他制作了中国第一具肝脏血管的铸型标本；他创造了间歇性肝门阻断切肝法和常温下无血切肝法；他完成了世界上第.