dxStefan数相变传热研究.pdfVIP

1992年第 5期 低 温 工 程 NO．5 1992 总第69期 CRYOGENICS Sum No．69 C d ^ q ， dxStefan数相变传热研究 ， 曲 鹏 ● _ — — — 一 T 毛Li (北京航空航无大学) 摘要 讨论 了受对流冷却 的相变传热 问题 应用近似分析方法得 出了小 Stefan 数下凝 目时间的近似解，并与精确解和摄功解作 了}E较。方法简捷、直观，所得结 果具有一定精度 进一步计论表 明，该分析方法叶任意边界条件下的相变 问题 有普 遍适用性，可解决更 奇复杂 的相变 问题 主题词 相变，移动边界 ，冷冻 。传热 Tt 相变温度[℃] 符 襄 号 T ， 流体温度[℃] 温度波动幅度[℃] Ta Tn的常数项[℃] 比热 [kTkg ℃ ] T =0表面温度[℃] 板厚度Im3 对流换热系数EkJm～s ℃ ] 希 腊 符 号 导热系数[kJr玎 S ℃ ] 口 导温系数Im S ] 相变潜热EkJkg ] 6 凝固层厚度Era3 边界热流 [kJtTt S’] p 密度Ikgm ] 时间[s] 目 凝 固时间Is3 温度 [℃] 角速度 [s ] 0 引言 固液问物态变化 (相变)的传热问题广泛存在于工程实践中，这一领域包括航天热控制， 飞机防冰、气动烧蚀，热能贮存、铸造冶金、化工、食品及生物组织冻结。在相变过程中， 存在分隔两相的固一液界面，其位置是时问的函数 ，且伴随着热量的吸收和 释 放，因而这一 问题是强非线性的，精确解只在极少数情况下才能得到 “。对复杂的几何形状、边界条件和 初始条件的相变 问题，只能通过近似分析或数值计算求解I”“。 对于工程问题，相变传热问题的简单分析解变得十分实际，吸引了大量研究者的注意力。 文献[4]回顾了这些工作。 在大围范内，相变问题的传热特性 由Stefan数确 定 (Ste=c(T，一r)／L)，dxStefan — l8 一 数代表了一大类工程实际问题 。在这些问题中，相变潜热在相变过程 中起主要作用．许多研 究者便基于这一特点构造相变传热的近似分析解 为此，一些研究者基于数值计算去推 测 “I 或通过将控制方程转化为一新独立变量方程求解，而使控制方程 更 为复 杂 “J或应用摄动 法，直接将解展开成 SLelan数或其它小参数的级数来求解”。小Stefan数的研究可参见文 献 [7，8]。但由于上述研究中分析求解过程复杂，在应用上受到了较大限制。 本文利用最基本的热传导近似解理论来分析相变传热问题，并将着重解决 小 Stefan效 的构造问题 。这一分析法是基于假设相变过程 中温度分布是稳态的，因而较其它方法简单。 实际上，对于小 Scefan数相变传热过程，这一假设是合理的。此分析方法的求解对任何外 部边界条件均有效，且可解决其它外形及初始温度不为相变温度的更为复杂的相变问题 。 本文通过一相当普遍、并不复杂的相变问