论文部分内容阅读
原子束源出口的流场决定了原子束源的工作效率,微管道内部流场影响着MEMS元器件的性能,动态气体锁内部的流场直接影响了其在真空腔室内对污染气体的抑制效果,所以对微结构内的流场和真空环境下流场的研究变得越来越重要。由于稀薄气体效应,传统流体仿真软件已经无法求解该类问题,并且直接求解流场的Boltzmann方程非常困难。本文使用了直接仿真蒙特卡洛方法(DSMC,Direct Simulation Monte Carlo method)求解微管道在稀薄气体环境下的流场问题,并提出了DSMC亚音速边界条件新的处理方法。然后使用DSMC方法及其亚音速边界条件研究了圆柱形毛细管出口流场的角度分布,并提出了在圆柱形毛细管出口增加圆锥形末端以减小管道出口分子束发散的方法。为了使用DSMC方法求解微管道内的亚音速流动问题,首先研究了亚音速边界条件,其中提出了恒流量入口和泵出口边界这两种亚音速边界新的处理方法。恒流量入口边界中提出了插入分子总数和分子在入口平面上分布的新处理方法。泵出口边界则由多孔出口边界实现,并提出了边界上分子移除概率新的计算方法。最后改进的恒流量入口、泵出口、常用的压力入口/出口和放气壁面等亚音速边界条件都被集成到开源流体仿真软件OpenFOAM的求解器dsmcFoam中。使用DSMC方法和集成在dsmcFoam中的亚音速边界求解了典型的二维矩形管道的内部流场,并将计算结果与DS2V软件的结果和解析解进行了比较,对比结果验证了压力入口、恒流量入口、压力出口和泵出口等亚音速边界条件的准确性。另外还计算了具有放气壁面的二维矩形管道的流场,并将其与DS2V的结果进行对比,对比结果验证了放气壁面条件的准确性。然后通过研究形状和边界条件更复杂的T形管内部流场,进一步验证了集成的亚音速边界条件的准确性。最后计算了相同入口流量下在滑移流态和自由分子流态下不同长度管道内部的流场。仿真结果在验证了边界条件准确性的同时,还显示改进的恒流量入口边界可以更精确的实现流量的设定值,并可起到减小流场计算区域的作用,而改进的泵出口边界则可在只设定抽速的条件下达到预期的出口压力。研究了稀薄气体流态下毛细管出口分子速度的角度分布。首先使用试验粒子Monte Carlo方法计算了自由分子流态下圆柱形毛细管出口分子速度的角度分布,在此基础上提出了增加圆锥形末端以减少出口分子束发散的方法,并得到了管道长径比、锥形部分的几何尺寸和CL反射模型的调节系数对出口分布的影响规律。然后使用DSMC方法和压力入口边界条件研究了更宽流域范围内圆柱形毛细管和带锥形末端毛细管出口分子速度的角度分布,并分析了气体的稀薄程度以及CL反射模型和VSS碰撞模型系数对出口速度角度分布的影响。仿真结果显示在宽的流域范围内圆锥形末端在不同的反射模型系数和碰撞模型系数下都具有减小出口分子速度角度分布发散的作用。使用DSMC方法和恒流量入口边界研究了不同入口流量下圆柱形毛细管出口外分子密度的角度分布及分布的半峰全宽(FWHM),并将计算结果与已有实验结果和理论值进行了比较。最后还研究了圆锥形末端和测量位置对出口外分子密度角度分布的影响。仿真结果与实验结果的对比表明使用DSMC方法和恒流量入口边界可以准确计算圆柱形毛细管出口分子密度的角度分布;圆锥形末端虽然会增加出口分子密度角度分布的FWHM,但是却可以在矩离出口较远的位置增大小极角处的分子密度,并减小大极角处的分子密度,达到减小出口分子束发散的作用。