微通道内给定壁面热流的气体换热研讨.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编廿．：063095 微通道内给定壁面热流的气体换热研究半 索晓娜闰孝红曾敏王秋旺徐进良 西安交通大学动力工程多相流国家重点实验室，西安，710049 摘要本文利用实施给定热流边界条件的DSMC方法，对短通道内给定壁面热流边界条件下的气体换 热情况逍行了模拟。结果表明，壁面热流密度增大导致通道内压力分布非线性明显。随着热流密度的 增大，截面速度分布趋于甲缓，壁面处速度梯度增大，滑移速度增大。给定热流密度的通道内壁面温 度和通道截面气流平均温度的差值沿程增大，温度梯度沿程下降，并且稀薄性增大换热降低。 关键词DSMC，给定热流边界，换热特性，微通道 l引言 微尺度下的对流传热问题广泛地存在于仪器科学、微电子学、生物工程及一些先进 微能量系统中IlJ。直接模拟蒙特卡罗(DSMC)方法对于研究稀薄气体的流动传热问题 是一种有效的工具12J。由于MEMS运行时的自身发热，相应壁面处给定热流的边界条件 在微通道问题中是普遍存在的，因此研究微通道内稀薄气体在给定壁面热流边界条件下 的流动与换热特性有着十分重要的意义。贺群武等在文ffa天131中采用逆温度抽样算法发展 了DSMC方法下给定热流边界条件的实施方法，计算结果表明，该方法能够准确抽样 得到壁面沿程温度，同时也能准确计算得到流场内压力、速度分布。本文将采用这一方 法分析研究壁面处给定热流密度的微通道内气体换热情况。 ． 2三种壁面边界条件的统一处理 为了研究不同壁面边界条件下的换热特性，将三种壁面边界条件进行统一。初始化 时假设壁面温度等于外部对流流体温度％运行一段时间流动达到稳态，然后根据壁面 处热鼍守恒式(1)更新壁面温度得到式(2)，假定忽略轴向导热。 q。=办，(弓一瓦)=忍。(乙一乙) (1) 瓦=弓一吼／乃 (2) 上式中外部对流流体温度乃，外部对流换热系数幻均为已知常数。当办，较大时，用 式(2)更新壁面温度，其中壁面热流‰由分子运动参数统计，并随着壁面温度不断更新。 当胁较小时，需要先假定壁面温度，然后由(1)式计算出壁面热流量，之后采用文献【3】 提到的ITS方法抽样出更新的壁面温度，直到达到稳定值。计算表明当以小于105时， 可采用(2)式进行计算。给定壁面温度的边界条件与绝热的边界条件分别是外部对流边界 条件的两个特例。当厅，趋于0时，外部对流边界变成绝热边界；当所趋于无穷大时，对 流边界变成给定壁温的边界。 为了研究方便，取算例通道为短通道，尺寸为5×1I．tm2，出El压力为肛larm，进出 口乐比为P，n／Po=|．2，通道壁面温度与进口流体温度均为300K。计算：【质为氮气，分子 之间作用势采用VHS模型，分子与壁面的碰撞采用完全热适应漫反射模型。 对丁微通道内的气体换热，改变外部对流换热系数时计算所得壁面温度分布及通道 基金项目：困家自然科学媾金资助项目． 599 截面平均气体温度如图l所示。由图可观察到以下现象：1)当外部对流换热系数办≯106 W／(m2-K)时，壁面温度与外部环境温度乃非常接近，通道内气体温度沿程分布与给定壁 温边界条件下的气体沿程温度分布几乎一致。这与前面所述的当向，趋于无穷大时，对流 边界转变成给定壁温的边界的理论相符。2)给定壁温边界条件下，通道内的气体截面 平均温度上升最快，但整体换热效果最差。且在通道后半段，气流温度与壁面温度相当， 传热量极小，近似绝热流动【4．5】。3)当缸106W／(m2．K)时，气流温度和壁面温度同时沿 着流动方向缓慢升高，在本文计算工况下，通道长度不足以使得壁面温度升高至外部环 境温度。当以非常小，如矗尸100W／(m2K)，流动近似为绝热流动，通道内的气体在压差 驱动下流动。壁面温度和气流温度沿着流动方向温度上升非常缓慢。4)各种计算工况 下都出现了在通道