一类热传导方程移动边界问题的数值求解.pdf

第28卷第12期 商丘师范学院学报 V01．28No．12 OF NORMALUNIVERSITY Dec． 2012 2012年12月 JOURNAL SHANGQIU 一类热传导方程移动边界问题的数值求解 李 玉 山 (甘肃政法学院计算机科学学院，甘肃兰州730070) 摘要：给出了数值求解具有多层区域热方程移动边界问题的算法，算法是基于有限差分格式的，给出处理界 面条件和曲边的方法，分析并且给出了最佳的处理方案；最后用数值例子说明我们的方法具有较高的精度． 关键词：热传导方程；多层区域；有限差分 中图分类号：0241．82 文献标识码：A Numericalmethodof foraclassofheatconduction movingboundaryproblem equation LIYushan ofPoliticalScience of Institute and (School Science，Gansu Law，Lanzhou730070，China) Computer thefinitedifferencemethodtosolvetheheatconductionwith use Abstract：We equationmuhilayerdomain，give methodstodealwiththeinterfaceconditionand andobtainitsbest movingboundary works numericalresultsshowthatour method effectively． proposed conduction difference words：heat Key equation；muhilayerdomain；finite 0引言 热传导方程在实际应用中，得到了很多的关注，这方面的技术广泛地应用在钢铁、电子和材料等方面．在 其中有一类重要且困难的问题就是复合材料中经常出现的多层问题．这类问题通常会和高温环境下的非线 性边界条件有关¨。2j，这给问题带来了一定的难度． 关于热传导方程的研究，已经有了非常丰富的经典理论结果∞“J．但是针对不同的模型，还是要不同的 解决方法．特别是当出现不同的非线性边界条件时候，就要有特殊的解决方法．在多层热传导问题中，经常出 现一类重要的边界条件就是线性条件，针对这方面的研究，大多是关于一层模型，本文考虑具有线性连接条 件的一维三层热传导模型．数值求解热传导方程有很多的数值方法，主要的方法有有限差分法，有限元和边 界元等等．针对本文的问题，选取了其中相对简单的有限差分法，并且对特殊的连接条件和移动边界，分别采 用中心差商近似和插值近似来处理．数值结果表明，本文采用的数值方法是有效的． 1问题的提出 具有三层区域热方程移动边界的问题可表示如下： 在每一层区域内的温度分布满足方程 警㈠小n；≥㈠叭舢一_1，2，3， (1) 边界条件和连接条件是 收稿日期：2012—07—19；修回日期：2012-08—04 作者简介：李玉山(198