论文部分内容阅读
稀薄气体流动广泛存在于航空航天飞行和微/纳机电系统中,高空稀薄效应对飞行器气动力和气动热的计算影响严重,对飞行器的控制、防热有着重要意义;而微/纳尺度下的大表面积比和表面力主导的效应,与微/纳电子元器件的散热和阻力计算都有着密切关系,需要给予极大的关注和研究力度。以Kn数为标准,可以将稀薄气体流动分为滑移流区、过渡流区和自由分子流区,由于每个流区中气体分子间相互作用和气体分子与壁面分子间作用的所占比重不同,所以需要用不同方法进行建模处理和求解,由此导致了目前稀薄气体流动研究中的问题和难点所在——滑移边界条件的形式不统一、Kn层内流动规律尚不明晰、分子与表面作用模型不准确,研究清楚这些问题对于稀薄气体流动的精确建模有着很大的帮助。分子动力学方法是一种确定性的粒子模拟方法,不同于DSMC,它从最基本的物理规律出发,应用牛顿运动定理来求解每个粒子的位置、速度和轨迹,并用统计的方法得出系统的宏观物理量,因此更适宜对稀薄气体的流动机理进行模拟研究和分析。基于上述背景,在本文中,我们采用分子动力学方法进行建模,用开源的LAMMPS程序搭建了计算平台,对微/纳通道中的Couette流动进行了模拟,研究了Kn数、剪切速度、温度和壁面参数变化对稀薄气体流动的影响,计算了壁面滑移速度随这些参数变化的规律,得出了滑移速度与Kn数、温度成反比,与剪切速度成正比,而壁面参数的影响相对较弱,且比较复杂,但对微/纳流动的影响也是不可忽略的结论。