基于Sigmoid函数光滑化的等效热容和有限元法求解相变传热问题.pdfVIP

第l9卷 5期 应用基础与工程科学学报 V01．19．No．5 2011年10月 JOURNALOFBASICSCIENCEAND ENGINEERING October2011 文章编号：1005-0930(2011)05-0817-013 中图分类号 ：0241 文献标识码：A doi：10．3969／j．issn．1005-0930．2011．05．015 基于 Sigmoid函数光滑化的等效热容和有 限元法求解相变传热 问题 何宜谦 ， 杨海天2 (1．大连理工大学建设工程学部土木工程学院，辽宁 大连 116024；2．大连理工大学工业装备结构分析国家重点 实验室，辽宁 大连 116024) 摘要：为避免等效热容法中因等效热容函数阶跃变化对计算精度的影响，提出了 一 个基于Sigmoid函数的光滑焓函数，通过焓对温度的导数，可得到光滑和连续 的等效热容函数，并结合有限元法与有限差分法建立了求解相变传热问题的数 值模型．通过数值算例，对所提方法进行了验证，探讨了有限元网格密度、热容阵 的数值积分格式以及Sigmoid函数相关参数对计算结果的影响，与解析解及 以 往的平均化方法相比，结果令人满意． 关键词：传热；相变；等效热容法；Sigmoid函数；有限元方法 在 自然界和各种工业生产过程中，很多问题与相变传热密切相关．例如，冰的融化和 凝固…，铸造过程中金属构件的凝固 J，晶体的生长_3J，相变材料的应用_4剖等． 相变传热问题中随时间移动的两相界面是造成相变传热问题求解困难的主要原因． 由于这类问题一般难于解析求解，大量工作集中于数值求解方法的研究，这些数值方法主 要可分为前沿跟踪法 (Front—trackingmethods)和固定区域法 (Fixed-domainmethods)L61． 前沿跟踪法连续不断地追踪相变界面位置，具有较高的精度，但为了追踪移动界面， 需要连续地进行网格重构、坐标变换或采用时空有限元技术 J，因此前沿追踪法的计算 量一般相对较大而且操作相对复杂．另外，前沿追踪法不适于求解非等温相变(具有相变 温度区间)或相变界面形状较为复杂的问题 J． 固定区域法将两相区域看作一个连续的区域，将相变界面的边界条件作为一个另外 的物理量(如焓，等效热容等)隐含在整个区域的控制方程中 J．和前沿追踪法相比，固定 区域法不需要连续地追踪移动界面和进行网格重构，因此计算量较小，操作较简单，并能 够处理具有相变温度区间的非等温相变问题 1-2,6]．根据潜热处理方法的不同，固定区域 法可进一步细分为焓法、等效热容法、热流法等 J． 在固定区域法中，等效热容法是被研究得最早和应用最为广泛的一种数值方法 J． 收稿 日期：2010-03-24；修订 日期：2011~04-20 基金项 目：国家自然科学基金项 目10721062；国家重点基础研究规划项 目[2005CB321704， 2010CB832703]；辽宁省中青年学术带头人基金 ；中央高校基本科研业务费专项资金资助 作者简介：何宜谦(1983一)，博士后．E-mail：heyiqian@hotmail．Corn 通讯作者：杨海天(1958一)，教授 ，博士生导师．Tel：0411E．mail：haitian@dlut．edu．ca 8l8 应用基础与工程科学学报 Vo1．19 其早期研究工作可见Vollerj『，Lewis等 J，Dalhuijsen和SegalL8等学者的综述性文章．近 年来也有不少的相关文献报道：高志华等_9将等效热容法应用于青藏铁路地基相变过程 的数值模拟中，Chang等 叫将等效热容法应用于硅薄膜镭射结晶过程的数值模拟中， 谗壤较 ＼ Li 将等效热容法应用于冷冻外科手术的数值分析中，Purlis等 将等效热容法应用于 鲫 ∞ ∞ ∞ ∞ 加 m O