球内凝固过程的近似解分析.pdfVIP

维普资讯 第 19卷第 4期 青 岛大 学 学报 (工 程 技 术 版 ) JOURNALOFQINGDAOUNIVERSITY (ET) Vo1．19NO．4 2004年 12月 Dec．200 4 文章缩号 ：1006—9798(2004)04—0080—06 球内凝固过程的近似解分析 杨 娟 ，杨启容 ，吕剑锋 (1．青岛大学机 电工程学院，山东青岛266071；2．青岛朗讯通讯设备有限公司，山东青岛266071) 摘要：对球 内凝固相变过程作了较为系统和全面的理论分析 ，建立了能够揭示在第三类 边界条件下球体内相变传热过程物理机制的数学模型，用 Lighthill奇异摄动法给出了球 体内对称凝固问题的近似分析解 ，并以Ba(OH) ·8HO为相变材料，得到了相变界面、 相变介质温度随时间和热流密度随相变位置的变化规律 ，以及球囊半径和对流换热系数 对凝固过程的影响。得到结论：对一般工程背景，球囊半径为 1．5～2．0cm，h1000 w／(m ·K)较为合适，从而为以球为相变贮能单元的贮能装置的优化设计提供了一定的 理论参考。 关键词 i球 内凝固；奇异摄动法；相变贮能 中图分类号iTK123 文献标识码：A 相变贮能换热器有多种形式，如管簇式、球形和平板式、堆积床式。然而 ，对球形贮能换热器而言，它有 结构简单 、安装维修方便、可靠性强 、阻力小、换热性能好 、价格便宜及适于系统扩充等特点。所以，对球形相 变贮能换热器的研究是较有意义的。 对于球形相变贮能换热器，相变材料密封于球体内，其固液界面的移动无法直接观测到。因此 ，对于球 形相变贮能器的测控要比盘管式的系统困难，所 以，很有必要从传热学的角度对球体内部的相变过程进行理 论分析。可采取数值法和分析解法去求解相变传热问题 ，以获得相界面的移动规律 、球内温度分布以及相关 参数对其换热的影响，为系统的优化设计提供理论依据 。但采用数值分析只能获得离散点处的数值解，并且 数值方法对于工程实际来说过于繁琐 ，同时存在着调试复杂、不易确定各参数间的关系等特点。近似分析法 虽有许多种 ，但无论哪种方法其共同的思路就是抓住问题的主要方面而忽略其次要方面，具有简捷 、可靠的 优点，可清晰地描绘相界面与时间的函数关系，得到球内凝 固时的温度分布。近几十年来，摄动法在解决量 子力学、流体力学 、天体物理以及工程技术中出现的大量非线性问题时取得 了很大的成功，并显示出优越性。 与数值法相比，摄动法可以提供各种参数间近似分析关系式 ，与其它近似方法相 比，它又有较大的灵活性 。 因此 ，近年来有不少人试用这一方法解 Stefan问题 ，通常是选择 Ste数为摄动量 ，将温度及界面移动位置按 Ste数展开 ，利用 Ste：0，即：忽略显热时的分析解——准稳态解作为零阶近似，求得它的逐项近似解 。Ste 数表示材料在相变过程中显热与潜热之比，亦即表面上所得到的加热量用于加热液体部分相对大小的一种 量度 。对于一般工程 问题，Ste1。 对于圆柱或圆球，直接用正则摄动法得到的解在界面移动到中心时不是一致有效的，本文采用 Lighthill 奇异摄动法[1]，引入 2个新的变量 ，以消除中心点的奇性 ，得到适用于整个凝固过程的一致有效解。 1 模型的建立 1．1 物理模型的建立 在实际工程中，许多装满相变材料的小球堆积在蓄热罐中，放热时，蓄热系统 中的低温液体流经球体与 收稿 日期 ：2004一O6一O1 作者简介 ：杨娟 (1980一)，女 ，在读硕士研究生，研究方向为车辆 电子技术 。 维普资讯 第 4期 杨 娟 ，等 ：球内凝固过程的近似解分析