热辐射计内稀薄气体流动的蒙特卡罗模拟.pdfVIP

　第 22 卷 　第 4 期 空 气 动 力 学 学 报 Vol . 22 , No . 4 　 　2004 年 12 月 ACTA AERODY NAMICA SINICA Dec . ,2004 　 ( ) 文章编号 : 2004 热辐射计内稀薄气体流动的蒙特卡罗模拟 吕章德1 ,2 ,王绍民1 ( 1. 浙江大学物理系 ,浙江 杭州 　310027 ;2 . 绍兴文理学院物理系 ,浙江 绍兴 　312000) ( ) 摘 　要 :本文应用直接仿真蒙特卡罗 DSMC 方法 ,研究了 Crookes 热辐射计内由于叶片两表面的温度不同而引起的 ( ) 稀薄气体的流动 ,给出了热辐射计叶片高温面温度 900K ,低温表面温度为 300 K , 四种不同克努森数 Kn 下流速矢 量场 、温度场和叶片两面压力差的仿真结果 。并与热辐射计在不同 Kn 数下的转速实验测量结果进行了比较 ,两者 结果基本一致 。 关键词 :稀薄气体动力学 ; DSMC 方法 ; 热辐射计 ; Boltzmann 方程 中图分类号 :V211. 3 　　　文献标识码 : A 0 　引　言 性决定的 ,例如系统的结构 ,渐近理论不能描述在较 大 Knudsen 数下的气流行为 。本文应用 DSMC[ 15 ] 方 　　Crookes 热辐射计是由英国化学物理学家 Willian 法 ,讨论了热辐射计中铝箔叶片两表面具有不同的均 Crookes 于 1872 年发明的[ 1] ,在物理上主要用于演示 匀温度分布下、在较大 Knudsen 数范围内的气体分子 和定性测量光辐射的能量 。它是一个密闭的具有一 速度流场 、温度场和辐射计叶片表面的压力差 。并对 定真空度的玻璃球形容器 ,球形腔内有一个固定四片 不同 Kn 下的稀薄气体流动的物理机理作了解释 。 方形叶片的转子 ,转子被顶在一根被固定在球形腔内 ( ( ) ) 1 　基本方程和边界值条件 支架上的针尖上 参见图 5 ,叶片由铝箔制成 ,铝箔 的一面涂黑 , 另一面则涂银 。当辐射计受到光照射 ( ) 　　本文用限制在二维容器 如无限长的矩形管 的 时 ,转子会朝黑面后退的方向转动 。 静止稀薄气体来仿真热辐射计 ,其横截面如图 1 所 　　对热辐射计的工作原理的定性分析已由Maxwell ( 示 ,0 X 1 L , 0 X 2 L , 其中 X i 是直角坐标 X 3 于 1