空泡流现象的非平衡分子动力学模拟和机理研究

论文摘要

流体中物体超高速运动的微观机理是现代科技最前沿的研究领域之一,无论在军事还是在民用技术中都具有极高的应用价值,特别是超空泡现象的研究,对水下航行体的减阻效应和超高速鱼雷的设计至关重要,有着重要的科学价值。因此本工作得到了国家自然科学基金项目《超高速流体的汽液相变与界面现象研究》（20473073）的资助。本论文运用非平衡分子动力学模拟方法,研究了微观体系条件下流体分子与固体壁面高速相对运动过程中流体流粒子的分布规律、状态性质、界面现象和摩擦情况等,特别是对开放体系超空泡现象的计算机模拟的实现,对超空泡流动中分子运动规律的观察和微观机理的探讨,为真实环境中物体超高速运动的研究提供了重要的基础数据。为了逐步深入研究不同体系不同条件下的流动,我们模拟了从简单窄孔到开放体系,从简单分子到甲烷和水等真实流体,从单系统到复合系统等模型,应用NEMD方法计算了流体压力、密度、速度、界面张力、摩擦系数等物理化学性质,通过对流场细节的观察,提出了一系列有关高速流动现象的形成机理,并与宏观流体力学理论进行了关联。为了验证分子动力学模拟程序和所建分子模型,本工作首先利用所建立的NEMD模型,代入平衡流动条件,模拟了宏量条件下甲烷的PVT关系,模拟结果与实验测量值或文献模拟值都能较好地吻合,说明本工作的模拟方法和设置基本正确。随后,对微孔无阻碍流进行了非平衡模拟,获得了符合宏观水动力学的模拟结果,验证了模拟程序与模型在非平衡流动模拟中的可行性。在上述研究的基础上,通过设定阻碍物模型,模拟计算了受限于窄缝孔中甲烷流体在160K温度、不同力场作用以及不同阻碍物条件下的阻碍流流场细节。模拟结果发现在高速条件下阻碍流会在流经阻碍物时形成高密区和低密区,同时引发相应的压力变化。流速对低密区形成有重要影响;阻碍物前端形态会影响流体两相形成和前端压力梯度,因此钝头阻碍物与其他流线形阻碍物相比将受到更大的压差阻力,但对两相界面表面张力影响不大;增加阻碍物长度或外力场强度会增大两相间表面张力。同时流动中观察到的摆动波现象符合Rayleigh-Taylor理论,表明该模型可以作为数值模拟的补充,研究宏观现象的微观机理。为了将流体模型应用到更有实际意义的超空泡流动模拟中,本工作引入了新的开放体系模型。通过将EMD模拟体系与NEMD模拟体系结合,模拟系统由硅酸盐狭缝孔扩展到开放体系,壁面限制作用被消除,外加力场改为恒流速流动,使得流体更为贴近真实流动状念。通过引入局部空化数σl概念,模拟体系被划分为多个独立方格,以观察超空泡的形成过程以及不同因素对超空泡现象的影响。模拟结果表明,微观条件下,空泡流的传统判据空化数σ＜0.1仍然成立。产生低局部空化数的σl区域与实际上的空泡发生区域在空间上并不重叠;低σl区域可以引起微小气泡的形成,虽然并不足以成为气化核,但也能给流体带来周期性的密度降;当低σl区域在运动速度方向上的尺度达到2流体分子直径时,可以在其后方形成稳定的超空泡;水下物体运动速度是引起超空泡的决定因素。在以上模拟工作成功的基础上,进行了水分子的超空泡模拟,试图进一步挖掘该模型的实用价值。应用SPC/E水分子模型,模拟并研究了不同空化器设置和不同运动速度下,超空泡流动的密度分布、局部空化数分布、空泡含气百分比和表面摩擦系数,并于数值模拟结果进行了对比。统计结果表明,平头空化器在相同流速条件下较为容易产生优良超空泡,空泡内含气百分比较高,能进一步降低壁面摩擦力。空化器外形对减阻效应有很大影响;超空泡现象可以降低航行体50%-90%的表面摩擦力。同时通过与数值模拟结果的对比,验证了模型的可行性。本论文完成了既定的研究目标,为日后进一步研究更加复杂和更加实际的空泡流体系,以及为研究流体复杂流动状态的微观机理提供了有力的研究工具。

论文目录

• ABSTRACT
• 摘要
• 第一章 引言
• 1.1 空泡流及超空泡现象
• 1.1.1 空泡流概念
• 1.1.2 空泡流及水下超空泡的国外研究
• 1.1.3 空泡流和水下超空泡现象的国内研究
• 1.2 计算机模拟方法及应用
• 1.2.1 Monto-Carlo模拟
• 1.2.2 平衡分子动力学模拟
• 1.2.3 非平衡分子动力学模拟
• 1.3 模拟体系
• 1.3.1 微孔体系
• 1.3.2 阻碍流和空泡流等复杂流动体系
• 1.3.3 水和甲烷
• 1.4 选题依据及主要内容
• 1.4.1 本工作选题依据
• 1.4.2 论文主要内容
• 本章符号列表
• 本章参考文献
• 第二章 分子动力学模拟
• 2.1 分子动力学的理论基础
• 2.2 分子动力学的基本算法
• 2.2.1 初始化设置
• 2.2.2 Lennard-Jones相互作用势能
• 2.2.3 Mishin镶嵌原子势
• 2.2.4 不同种分子间相互作用势能
• 2.3 NEMD统计热力学实现
• 2.3.1 解牛顿运动方程
• 2.3.2 NEMD中热力学量的统计
• 2.3.3 NEMD温度校正
• 2.3.4 欧拉角、转动矩阵和四元数
• 2.3.5 力矩
• 2.3.6 长程作用力处理
• 2.4 分子动力学模拟模型
• 2.4.1 周期边界条件
• 2.4.2 最小映像法
• 2.4.3 壁面模型
• 2.4.4 水分子模型
• 2.5 其他模拟技巧
• 2.5.1 约化单位
• 2.5.2 时间步长
• 2.5.3 加速算法
• 2.6 本章小结
• 本章符号列表
• 本章参考文献
• 第三章 微孔阻碍流的分子动力学模拟
• 3.1 前言
• 3.2 模拟细节
• 3.2.1 NEMD压力计算
• 3.2.2 NEMD界面张力计算
• 3.2.3 外加力场
• 3.2.4 NEMD模拟模型
• 3.2.5 其他参数
• 3.3 结果与讨论
• 3.3.1 EMD宏量模拟和流体模型验证
• 3.3.2 圆柱阻碍物微孔NEMD模拟
• 3.3.3 三棱柱阻碍物微孔NEMD模拟
• 3.3.4 射弹型阻碍物微孔NEMD模拟
• 3.3.5 表面张力计算
• 3.4 本章小结
• 本章符号列表
• 本章参考文献
• 第四章 开放体系下甲烷流体自然空泡流的NEMD研究
• 4.1 前言
• 4.2 模拟细节
• 4.2.1 开放体系NEMD模型
• 4.2.2 空化数及局部空化数计算
• 4.2.3 其他模拟参数
• 4.2.4 NEMD超空泡流动模拟流程图
• 4.3 结果与讨论
• 4.3.1 阻碍流模型模拟与验证
• 4.3.2 三角形空化器不同速度下的空泡流模拟
• 4.3.3 不同空化器临界速度下的空泡流模拟
• 4.3.4 不同流体密度下的空泡流模拟
• 4.3.5 周期震荡现象
• 4.4 本章小结
• 本章符号列表
• 本章参考文献
• 第五章 开放体系水中空泡流的NEMD模拟及自然超空泡减阻机理研究
• 5.1 前言
• 5.2 模拟细节
• 5.2.1 水分子模型
• 5.2.2 刚体分子模型转动惯量计算
• 5.2.3 其他模拟参数
• 5.3 结果与讨论
• 5.3.1 三角形空化器超空泡流模拟
• 5.3.2 不同空化器形态对超空泡的影响
• 5.3.3 减阻效果比较
• 5.4 本章小结
• 本章符号列表
• 本章参考文献
• 第六章 结论与展望
• 6.1 结论
• 6.2 展望
• 致谢
• 附录 已发表论文