数学论文.doc

《数学论文.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、基于高斯函数形式的热平衡积分法AGaussianFunctionBasedHeatBalanceIntegralMethod袁国静，令锋1.内蒙古工业大学理学院，内蒙古呼和浩特010051；2.肇庆学院计算机学院，广东肇庆526061投寄：潍坊学院学报收稿日期：2010-06-22；修改日期：2010-基金项目：广东省自然科学基金资助项目（04011600）作者简介：袁国静（1983-），女，山东潍坊人，内蒙古工业大学与肇庆学院联合培养硕士研究生。通讯作者：令锋（1963-），男，陕西岐山人，教授，博士，l

2、ingf@zqu.edu.cn基于高斯函数形式的热平衡积分法袁国静，令锋（1.内蒙古工业大学，内蒙古呼和浩特010051；2.肇庆学院，广东肇庆526061）摘要：采用高斯函数形式的温度分布，应用热平衡积分法及热平衡积分法的细化方法求得单相融化问题的近似解。通过讨论得知，当数的倒数小于1时，热平衡积分法的细化可有效提高精度。关键词：热平衡积分法；细化；高斯函数中图分类号：O241.82文献标志码：A0引言热平衡积分法（HBIM）是1958年由Goodman最早提出的一类可以用来生成由微分方程控制的热传导方程

3、的近似解的半解析的有效方法。但是，HBIM的计算结果对选取的温度分布函数具有敏感性，Goodman最先提出选用二次多项式作为温度分布函数，同时也简单提及三次函数的情况，有算例表明，选择三次函数作为温度分布函数的精度却比二次多项式的结果更差，为了减少计算结果对温度分布函数的依赖性，提高计算结果精度，Noble提出了HBIM的细化（Refinement）方法。Bell据此法求解了伴有相变的一维问题，获得较好的结果。Mosally等人以单相融化问题为模型，采用与数值实验结果相对比的方法，给出了当数等于1时单相融化

4、问题HBIM细化解的收敛速度。徐湘田和令锋通过理论分析而不依赖数值试验，分别给出了单相融化问题和Neumann问题的热平衡积分细化解收敛性的证明。Gauss函数在Hermite多项式的定义中起着重要作用。热传导问题导热过程的特征促使人们设想指数函数可能是适当的温度分布函数的选取形式，但计算表明，选择指数函数作为温度分布函数的精度反而比选取二次多项式的结果要差，所以，选取高斯函数作为温度分布函数可能是一种有益的尝试。本文以单相融化问题为模型，以高斯函数为温度分布形式，应用热平衡积分法及热平衡积分法的细化方法求

5、解单相融化问题，从而得到细化可提高精度的条件。1模型方程研究的定解问题采用半无限大介质中固体融化问题的无量纲形式，初始温度为它的融化温度，数学描述如下：，，（1），（2）（3），，（4），，（5）其中，其中称为数。方程式（1）—（5）的解析解为，，（6），（7）其中是超越方程（8）的根，称为融化参数。2高斯函数形式的热平衡积分解方程式（1）两边关于在内应用热平衡积分法，得到热平衡积分方程（9）把式（5）代入上式得（10）研究高斯函数形式热平衡积分方程近似解，令（11）参量为常数，由式（11）得，（12）把式

6、（12）计算结果代入式（10）得（13）由式（11）得（14）由式（14）和式（5）得（15）由式（3）（4）得，（16）联立式（13）（15）（16）得（17）由式（15）（16）得（18）（19）当，取不同值时，融化界面的精确解、高斯函数形式热平衡积分解、二次多项式形式热平衡积分解作比较，如表1所示表1精确解，高斯函数形式热平衡积分解和二次多项式形式热平衡积分解比较精确解高斯函数形式二次多项式形式数值解相对误差数值解相对误差0.22.11942.0698-0.02342.21520.04520.41.7

7、2481.7079-0.00981.79540.04090.61.50241.4942-0.00551.55500.03500.81.35171.3470-0.00341.39260.030311.24021.2372-0.00241.27300.02641.21.15191.1508-0.00181.18000.0235由表1可见，高斯函数形式热平衡积分解比二次多项式形式解有更高的精度。3高斯函数形式的热平衡积分细化解将区间等分，得到个长度为的子区间，，，令：，，(20)由上式可知：未知参数为个。（1）据

8、式（3）（4）得，（21）（2）据在节点上的连续性得当时，节点上温度近似为当时，节点上温度近似为所以，（22）（3）同理，据在节点上导数的连续性得，（23）（4）当时（24）由式（24）和式（5）得（25）（5）在内对积分得个热平衡积分方程，（26）由式（20）得，，（27）把式（27）结果代入式（26）得（28）令式（28）等于式（25）得代数方程（29）（6）在区间上（30）由式（20）得，（31）把式（31