蒙特卡罗方法中分子作用模型的设计及应用精选.pdf

第 26 卷 第 1 期 计 　　算 　　物 　　理 Vol . 26 ,No . 1 2009 年 1 月 　CHINESE JOURNAL OF COMPUTATIONAL PHYSICS 　 Jan . , 2009 　 ( ) [文章编号] 　1001246X 2009 蒙特卡罗方法中分子作用模型的设计及应用 1 2 2 1 王学德 , 　伍贻兆 , 　夏 　健 , 　林晓宏 ( 1 南京理工大学动力工程学院 , 江苏 南京 　210094 ; 2 南京航空航天大学航空宇航学院 , 江苏 南京 　210016) [摘 　要] 　利用LennardJones (6 - 12) 作用势 , 设计 GSS3 分子碰撞模型 , 通过与已有的分子模型及实验数值 对比 , 该模型能够给出正确的气体输运性质. 将 GSS3 模型应用到非结构网格 DSMC 中 , 模拟二维超声速平板绕 流和圆柱绕流. 数值计算结果初步验证该分子模型的合理性与有效性. [ 关键词] 　蒙特卡罗方法 ; 分子作用模型 ; DSMC ; 稀薄气体 [ 中图分类号] 　V21125 [文献标识码] 　A 0 　引言 ( ) [ 1 ,2 ] 在求解稀薄气体动力学和微尺度流动问题中,直接物理模拟蒙特卡罗方法 DSMC 发挥了极为重要 的作用 ,也是迄今为止唯一获得巨大成功的方法. 其基本思想是用有限个模拟分子代替真实气体分子 ,在计 算机中存储模拟分子的位置坐标 、速度分量以及内能 ,其值随模拟分子的运动与边界的作用以及模拟分子之 间的碰撞而不断变化 ,通过对计算网格内模拟分子的运动状态取样统计达到求解真实气体流动的目的. 由于 DSMC 方法是直接以分子为研究对象的统计方法 ,因此为了获得气体真实的物理性质 ,研究者必须明确分子 之间的相互作用势 ,即分子作用模型. 本文首先以LennardJones 势为作用模型 ,在前人工作的基础上 ,设计了 GSS3 分子碰撞模型 ,该模型同时具有吸引和排斥的分子作用力 ,能够给出正确的气体输运性质. 然后给出 了该分子作用模型在非结构网格 DSMC 中的具体应用 ,验证该分子模型的合理性与有效性. 1 　GSS3 分子作用模型的设计 11 　DSMC 常用分子模型 ( ) [ 1] 在 DSMC 中 ,最早采用的是硬球模型 HS ,它不变的分子直径与碰撞截面简化了理论分析与数值计 ( ) [2 ] 算 ,但与物理事实不符 ,为此 ,Bird 提出了变径硬球模型 VHS ,VHS 模型假设总碰撞截面 、偏转角与粘性 ξ ξ ξ - 2 - + 12 σ σ ε χ ( ) μ ε 系数满足 T