DSMC方法下恒热流边界条件研究.pdfVIP

中国工程热物理学会 传热传质学 第十一届学术会议 编号：053097 DSMC方法下恒热流边界条件8 贺群武王秋旺索晓娜 (西安交通大学动力工程多相流国家重点实验室，陕西西安，10049) 摘要本文详细讨论了DsMc方法中流场温度、能量及边界热流的统计方法，发展了～种从边界热流 求得与壁面碰撞分子的反射速度的方法。该方法被称为ITs(InverseTcmpera【uresamplin只)方法。在此 基础上，本文发展，DsMc方法中壁面处恒热流边界条件的实施方法。计算结果表明：lTs方法能准 确抽样反射分子的特征温度，进而求得分子反射速度。基于该方法的给定热流边界条件可以准确求得 壁商处温度分布。以及流场内压力、速度。 关键词DsMc．ITs方法，给定热流边界，微通道 l引言 simu】ationofMonte 直接模拟蒙特卡罗DsMc(Direct cano)方法【1瞧研究稀薄气体 流动传热问题的一种有效。】：具。对于DsMc方法在微通道流动传热的模拟，人们开展 了J“泛的研究，也取得了明显的进展，但其研究大多没有考虑通道内的换热特性。在实 际应用中，由于MEMs运行时的自身发热，相对应壁面处给定热流的边界条件在微通 道问题中是普遍存在的【“，但目前为止尚未见到有关DsMc方法中恒热流边界条件处理 方法的报道。 2InVerse 7Ibmperature Sampling(ITS) 对单一气体，一定体积内气体宏观速度％为该体积内所有分子速度的平均值： co=√“02+v02+w02=√玎2+矿2+面2(1) 式中，“o、K、wo为流体宏观运动速度％分量。设c7为分子热运动速率，“’、v’、Ⅵ／ 为其三个分量，“、v、w为分子绝对运动速度c分量，则它们之间的关系为： c。=c—c0 co=c (2) 通常分子振动的特征温度都很高(如氨气分子的振动特征温度为3337K)，大多数工况 的温度不足以激发工质分子的振动能量，而且单原子分子气体没有转动自由度，因此对 于单原子分子气体，其能量仅还有平动能，即： 孑=_， (3) 亏，=去矾c“=去m(∥2+v。2+w。) (4) 当气体分子与固体壁面碰撞时，气体分子与固体壁面碰撞的充分度和碰撞能量均与 壁面性质有关，反射分子的特征温度可按下式进行调整11】： 吒=器 ㈣ 口。称为物面热协调系数，哎∈【o，1】，蓐、巨分别表示入射分子与从物面反射回气体 基金项目：国家自然科学基金资助项目 1553 流场的分子所具有的平均能量，E表示当反射回气体流场中的分子的温度与实际物面 温度r相等时气体分子所具有的平均能量。两种极限情况即：a，=l，称为完全热协调； 口，=o，称为完全不协调，此时E，=E，分子在反射前后温度不变。与此对应，根据 分子与壁面碰撞后分子反射速度不同，反射模型可以分为镜面反射(specujarRenection) 模型和漫反射(DimlseRenection)模型。对于镜面反射，分子与固壁碰撞后其速度平行于 同壁的切向分量保持不变，而法向分量的大小不变，方向与碰撞前相反：对于漫反射模 型，气体分子与同壁碰撞后的速度与碰撞前的速度完全没有关系，反射分子作为一个整 体在速度上满足MaxweIijan速度分布…： ／ 、； 腓n