论文部分内容阅读
空化是发生在液体中复杂的非定常流动物理现象,涉及到多相流、可压缩以及相间交换等多个方面,对水力机械、水利工程、船舶工程、水下兵器和核工业等众多领域均有显著的影响。为了更加准确的对空化性能进行预测和评估,很有必要针对空化初生时的瞬态机理特性进行研究,分析空化初生时的瞬态特性。本文分别以液氩、水作为微观研究对象,研究液氩、水内部单个空化核生长过程中的热力学参数变化。以NACA4412翼形为宏观研究对象,将空化机理微观研究结论应用到宏观二维水翼的空化瞬态特性变化中,对其进行数值模拟计算,并运用试验研究方法验证瞬态二维水翼空化的模拟计算结果的准确性。主要内容如下:1.查阅相关资料,介绍了空化、气泡成核等相关的理论基础知识及国内外最新研究进展。并介绍了分子动力学(MD)研究方法,以及相变和气体状态方程的理论基础知识,叙述了分子动力学(MD)分子力场、边界条件以及分子动力学(MD)模拟的求解过程,阐述了在计算过程中的各种系综的定义以及适用类型,并详细介绍了分子动力学软件的选择以及应用,同时采用LAMMPS软件完成本文第三、四章的模拟计算工作,在模拟计算中采用Lennard-Jones(12-6)力场,边界条件采用周期性边界条件。2.在NVT系综下,对存在不同初始尺寸空化核的Lennard-Jones流体进行了研究,得到NVT系综下空化核发展的演变过程,并对其分子势能、系统密度、分子径向分布函数、系统压力以及系统总能量等相关的热力学参数进行分析,结果表明:在空化初生阶段,由于流体的局部压力降低而造成液相分子不断的进入空化核,促使了空化核的生长。通过分子势能研究发现,空化核位置处的分子势能较高,而液相分子势能在空化核生长初期变化较大。在空化核成长后期,液相分子势能相对较为平稳。液相密度和界面密度变化比较明显,而空泡中心区域的密度变化甚微。液相区域径向分布函数的峰强度就越大,峰宽度越窄,该区域分子间距rbin越小。随着空化核初始尺寸的减小,系统压力和总能量达到平衡的时间越长,平衡后的压力值和能量值也就越小。而且在空化核生长初期,系统的总能量和压力的变化波动较大。3.在NPT系综下,对存在空化核的水分子计算域进行了研究,得到空化核初始尺寸与系统初始压力对空化初生的影响关系,并对空化核发展的演变过程、系统压力以及系统能量、分子径向分布函数等相关的热力学参数进行分析,结果表明:在系统初始压力一定时,空化核尺寸存在临界值,当空化核初始尺寸小于该临界值时,即使液体受到负压的影响,由于分子之间的作用力和氢键力束缚空化核不能生长,因此空化不会发生。空化在微观层面上不仅是压力和空化核参数的变化,更重要的原因是系统初始压力造成系统能量增大,结构失稳,从而造成空化核的变化。在内部空化核尺寸一定时,外部压力存在一临界值,当外部压力值大于该临界值时,计算域内部不会产生相应的负压,或者产生负压值较小,不足以使系统失稳,在分子之间的作用力和氢键力束缚下空化核不能生长,因此空化不会发生。伴随着空化核的生长,系统能量和分子势能总有增大的趋势,而该增大趋势影响了计算域内部结构的稳定性,促使空化核迅速增大。4.通过第四章中纳米尺度空化核生长变化过程以及目前研究的微米尺度空化核生长过程对Z-G-B空化模型进行了优化,并在此基础上选用RNG k-ε湍流模型以及优化后的空化模型对二维水翼空化初生过程进行模拟计算。在空化初生时,由于不同攻角产生的压力场略有不同,因此空化核周围液体失稳的可能性不同,造成水翼上表面空化核生长发育的几率不同。攻角越大,产生的低压极值越小,空化核周围液体所受到的负压逐渐增大,从而势能增大,空化核快速生长。通过微观空化核研究发现,在不同压力场作用下,空化核生长几率不同,因此在不同攻角影响下,空泡发展规模和强度不同。在空化发展过程中,因尾部高压产生的回射流对空泡的形状以及脱落有着重要影响。空化形成后期,低压和高压的相互影响直接造成脱落空泡破裂与发展,因而影响到水流的流态分布以及旋涡的产生。空化初生时只是受到低压作用的影响,而空化再生则受到低压作用和紊乱流场的双重作用。