对流换热边界条件下球体内固液相变过程移动热源法求解.pdf

第20 卷 第 1期 中 国科 学 院研 究 生 院学 报 Vol.20　No.1 2003 年3 月 Jour nal of the Graduate School of the Chine e A cademy of Science M ar.2003 对流换热边界条件下球 内固液相变过程 移动热源法求解＊ 白凤武　卢文强 (中国科学院研究生院物理系, 北京 100039) (2002 年5 月9 日收稿;2002 年6 月 12 日收修改稿) 摘　要　应用移动热源法求解了对流换热条件下球 内具有明显相变界面的固液相变问题, 分 析了球 半径、相变材料潜热、外界对流换热条件对相变过程的影响, 并进一步考虑实际相变过 程存在的变对流边界条件, 构造了一实用解法.为求解相变问题的数值方法提供了一校核其可 靠性的新途径, 还为实验结果提供了相互比较的数据.该方法及其结果对相变材料的分析与设 计具有重要的实用价值. 关键词　移动热源法, 固液相变, 对流换热 中图分类号　O551.3 　　相变材料(PCM )是工业上广泛应用的一种热控材料, 由于具有在相变过程中吸收或者放出相变潜 热, 而温度却变化很小的特点, 使其在航天器的热控制、电子设备冷却等领域发挥着重要的作用. 　　大多数纯质有机相变材料的相变传热过程具有确定移动界面, 按照移动边界的形成原因, 有固 的 熔化、液 的凝固及一些固 的热结和烧蚀等.由于相变传热过程存在一个随时间变化的相界面, 使得传 热区域被界面分成两个具有不同相态的部分.而且在相界面上具有潜热的吸收或释放, 造成界面两侧热 流不连续, 使得此类问题在数学上呈现强非线性, 解的叠加原理不能使用.目前仅对少数一维半无限大、 [1] 无限大区域且具有简单边界条件的理想化情形能够精确求解 .文献[2]应用移动热源法, 对球形相变材 料定壁温边界条件下的相变过程进行了求解, 并使其作为与实验结果相互比较的依据.而对于相变材料 在对流换热情况下的相变过程, 还没有见到理论解析结果和数值计算结果的报道, 使得该类问题的实验 和数值计算结果的分析缺少理论结果的比较. 　　本文就对流换热下的球 相变材料的换热过程进行了分析, 根据换热过程的具 特点, 结合格林函 数法和移动热源法, 对相变过程进行了解析求解, 很好地计算了相变界面随时间的推进情况, 分析了诸因 素对相变过程的影响, 为实际相变过程提供了可靠的理论结果, 并对变对流换热边界条件进行了分析. 1　问题分析及假设 　　在相变材料的诸多应用中, 往往需要考虑对流换热情况下液 的凝固或者固 的熔化过程.如 :潜热 型功能流 的强化换热过程中, 相变材料通过与周围流 间的对流换热, 进行凝固或者熔化过程, 进而放 [3, 4] 出或者吸收潜热, 达到强化换热的目的 .由于在功能流 换热过程中, 相变材料是封装在一些直径只 有几十到几百微米的微胶囊之中, 可以认为微胶囊在相变过程为球形.微胶囊处在外围流 的包围之中, 其外界对流换热系数以及流 温度都是随着换热过程发生变化的, 形成了一个典型的变对流换热边界条 ＊国家自然科学基金资助项目(No;国家重点基础研究发展规划项目(No.G200026300) 　1　　　4 中国科学院研究生院学报 第20 卷 件的相变问题.为了说明问题方便, 本文仅对凝固过程进行分析, 换热过程可以做如下描述. 　　初始温度高于凝固点温度的液 球, 置于低于凝固点温度的外界流 环境之中, 由于对流换热作用, 液 球被逐渐冷却, 当球 外表面降到液 的凝固点温度后, 发生凝固, 在进一步冷却下, 达到完全凝固, 直至与外界流 达到换热平衡, 即冷却到外部流 温度. [5] 　　解析求解相变问题的方法有