边界节点法计算二维瞬态热传导问题.PDF

　 ， ３５ ２ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１４年２月１５ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３５，Ｎｏ．２，Ｆｅｂ．１５，２０１４ 文章编号：１０００⁃０８８７（２０１４）０２⁃０１１１⁃１０ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 边界节点法计算二维瞬态热传导问题∗ 师晋红，　 傅卓佳，　 陈　 文 （１． 河海大学 力学与材料学院，南京 ２１００９８； ２． 河海大学 水文水资源与水利工程科学国家重点实验室，南京 ２１００９８） （我刊编委陈文来稿） 摘要：　 采用边界节点法（ＢＫＭ）结合双重互易法（ＤＲＭ）求解二维瞬态热传导问题．采用差分格式 处理时间变量，可将原瞬态热传导方程转化为一系列非齐次修正的Ｈｅｌｍｈｏｌｔｚ方程．随后，方程的解 可分为特解和齐次解两部分计算，引入双重互易法在区域内部配点求解方程的特解，采用边界节 点法仅需边界配点求解方程的齐次解．给出的数值算例显示该方法计算精度高，适用性好，具有很 好的稳定性和收敛性，适合求解瞬态热传导问题． 关　 键　 词：　 瞬态热传导；　 边界节点法；　 双重互易法；　 差分格式；　 修正的Ｈｅｌｍｈｏｌｔｚ方程 中图分类号：　 Ｏ２４１．８２；Ｏ３４３．６　 　 　 文献标志码：　 Ａ ｄｏｉ：１０．３８７９／ ｊ．ｉｓｓｎ．１０００⁃０８８７．２０１４．０２．００１ 引　 　 言 ［１］ 近年来，功能梯度材料、形状记忆合金等许多新兴材料被广泛用于高温涡轮发动机 及 ［２］ 伸缩翼系统 的热障涂层等航天工程中．这些材料经常需要在高温环境中长期连续工作，所以 很有必要了解这些材料及相关结构随时间演变的热传导性能，因此亟需发展计算瞬态热传导 问题的数值算法，以便为航天工程材料的研发设计提供参考与指导． ［３⁃９］ 目前，很多数值算法 已用于求解瞬态热传导问题，例如，有限差分法、有限元法、边界元 方法等．有限差分法和有限元法在求解问题时要求在整个求解域上进行离散，在处理复杂几何 区域等问题时存在计算量过大、网格生成困难等难题．边界元方法是一种边界型数值算法，将 计算域降低了一维，然而边界元法选用含奇异性的基本解作为插值基函数，不可避免地需要处 理费时费力的奇异积分计算问题． ［１０］ 最近几十年发展起来的无网格方法 ，基于点的近似，克服网格依赖的缺陷，减少因网格 ［５⁃７，１１⁃１２］ ［１３］ 畸变引起的计算困难，已广泛用于求解瞬态热传导问题 ．其中，基本解方法 是与边界 ∗ 收稿日期：　 ２０１３⁃１０⁃０７；修订日期：　 ２０１３⁃１１⁃１９ 基金项目：　 国家重点基础研究发展计划 （９７３ 计划）（２０１０ＣＢ８３２７０２）； 国家杰出青年基金 ；国家自然科学基金１１３０２０６９）；高等学校学科创新引智计划 （“１１１”计划）（Ｂ１２０３２） 作者简介：　 师晋红（１９８９—），女，山西人，硕士生（Ｅ⁃ｍａｉｌ：ｓｈｉｊｈｈｈｕ＠１６３．ｃｏｍ）； 陈文（１９６７—），男，江苏人，教授，博士，博士生导师（通讯作者． Ｅ⁃ｍａｉｌ：ｃｈｅｎｗｅｎ＠ｈｈｕ． ｅｄｕ．ｃｎ）． １１１ １１２