数学论文[精编]

《数学论文[精编]》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

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

2、静1'2,令锋2（1.内蒙古工业大学,内蒙古呼和浩特010051；2.肇庆学院,广东肇庆526061）摘要：采用高斯函数形式的温度分布，应用热平衡积分法及热平衡积分法的细化方法求得单相融化问题的近似解。通过讨论得知，当Stefan数的倒数0小于1时，热平衡积分法的细化可有效提高精度。关键词：热平衡积分法；细化；高斯函数中图分类号:0241.82文献标志码：A0引言热平衡积分法（HBIM）是1958年由GoodmanM早捉出的一类可以川來生成由微分方程控制的热传导方程的近似解的半解析的有效方法山。但是，HB1M的计算结果对选取的温度分布函数具有敏感性，Goodman最先提出选用

3、二次多项式作为温度分布函数，同时也简单捉及三次函数的情况，有算例表明，选择三次函数作为温度分布函数的精度却比二次多项式的结果更差〔27,为了减少计算结果对温度分布函数的依赖性，提高计算结果棉度，Noble提出了HBIM的细化（Refinement）方法。Bell据此法求解了伴有相变的Stefan问题片讥获得较好的结果。Mosally等人以单相融化问题为模型，采用与数值实验结果相对比的方法，给出了当Stefan数等于1时单相融化问题HBIM细化解的收敛速度⑹。徐湘III和令锋通过理论分析而不依赖数值试验，分别给出了单相融化问题⑺和Neumann问题⑻的热平衡积分细化解收敛性的证

4、明。Gauss函数在Hennite多项式的定义屮起着重要作用。热传导问题导热过程的特征促使人们设想指数函数可能是适当的温度分布函数的选取形式，但计算表明，选择指数函数作为温度分布函数的精度反而比选取二次多项式的结果要差⑼，所以，选取高斯函数作为温度分布函数可能是一种有益的尝试。本文以单相融化问题为模型，以高斯函数为温度分布形式，应用热平衡积分法及热平衡积分法的细化方法求解单相融化问题，从而得到细化町提高精度的条件。1模型方程研究的定解问题采用半无限人介质中固体融化问题的无最纲形式，初始温度为它的融化温度，数学描述如下⑹：dud2udtdx200(1)u(x

5、,0)=0,x>Q(2)u(0j)=1,r>0(3)u(x,t)=0,x=5(r),f〉0(4)ducds-=-/>—,x=5(r)，/>0(5)oxdt其中"士二不土T其中氐称为皿数。方程式（1）一（5）的解析解为⑹00(6)(7)英中2是超越方程(8)的根，称为融化参数。2高斯函数形式的热平衡积分解方程式（1）两边关于兀在［0,$］内应用热平衡积分法，得到热平衡积分方程du~dxdug（f）dx.t=0(9)把式（5）代入上式得qdsdu—0矿丟(10)x=0研究高斯函数形式热平衡积分方程近似解，令参量abc为常数，由式（11）得dv把式(12)计算结

6、果代入式（10）得-2bcclth(ec-2cec-1)+2cf3由式(11)dv=-(ec+2cec)x=s(t)S由式(14)和式（5）得ds_阳(l+2c)~dt~山式(3)(4)得d=1,a--bec=0联立式(13)(15)(16)得(1+2c)[2cec(l+0)-d+1]=2c0由式(15)(16)得(13)(14)(15)(16)(17)(18)(19)当t=lf0取不同值时，融化界而$的和确解、高斯函数形式热平衡积分解、二次多项式形式热平衡积分解作比较，如表1所示表1精确解，高斯函数形式热平衡积分解和二次多项式形式热平衡积分解比较0精确解高斯函数形式二次多

7、项式形式数值解相对误差数值解相对误差0.22」1942.0698-0.02342.21520.04520.41.72481.7079-0.00981.79540.04090.61.50241.4942-0.00551.555()0.035()0.81.35171.3470-0.00341.39260.030311.24021.2372-0.00241.27300.02641.21.15191.1508-0.00181.18000.0235由表］可见，高斯函数形式热平衡积分解比二次多项式形式解冇更高的精度