数学论文.doc.doc数学论文.doc.doc

收藏 分享

资源预览需要最新版本的Flash Player支持。
您尚未安装或版本过低,建议您

基于高斯函数形式的热平衡积分法AGAUSSIANFUNCTIONBASEDHEATBALANCEINTEGRALMETHOD袁国静2,1,令锋21内蒙古工业大学理学院,内蒙古呼和浩特010051;2肇庆学院计算机学院,广东肇庆526061投寄潍坊学院学报收稿日期20100622;修改日期2010基金项目广东省自然科学基金资助项目(04011600)作者简介袁国静(1983),女,山东潍坊人,内蒙古工业大学与肇庆学院联合培养硕士研究生。通讯作者令锋(1963),男,陕西岐山人,教授,博士,LINGFZQUEDUCN基于高斯函数形式的热平衡积分法袁国静2,1,令锋2(1内蒙古工业大学,内蒙古呼和浩特010051;2肇庆学院,广东肇庆526061)摘要采用高斯函数形式的温度分布,应用热平衡积分法及热平衡积分法的细化方法求得单相融化问题的近似解。通过讨论得知,当STEFAN数的倒数小于1时,热平衡积分法的细化可有效提高精度。关键词热平衡积分法;细化;高斯函数中图分类号O24182文献标志码A0引言热平衡积分法(HBIM)是1958年由GOODMAN最早提出的一类可以用来生成由微分方程控制的热传导方程的近似解的半解析的有效方法1。但是,HBIM的计算结果对选取的温度分布函数具有敏感性,GOODMAN最先提出选用二次多项式作为温度分布函数,同时也简单提及三次函数的情况,有算例表明,选择三次函数作为温度分布函数的精度却比二次多项式的结果更差32,为了减少计算结果对温度分布函数的依赖性,提高计算结果精度,NOBLE提出了HBIM的细化(REFINEMENT)方法。BELL据此法求解了伴有相变的一维STEFAN问题54,获得较好的结果。MOSALLY等人以单相融化问题为模型,采用与数值实验结果相对比的方法,给出了当STEFAN数等于1时单相融化问题HBIM细化解的收敛速度6。徐湘田和令锋通过理论分析而不依赖数值试验,分别给出了单相融化问题7和NEUMANN问题8的热平衡积分细化解收敛性的证明。GAUSS函数在HERMITE多项式的定义中起着重要作用。热传导问题导热过程的特征促使人们设想指数函数可能是适当的温度分布函数的选取形式,但计算表明,选择指数函数作为温度分布函数的精度反而比选取二次多项式的结果要差9,所以,选取高斯函数作为温度分布函数可能是一种有益的尝试。本文以单相融化问题为模型,以高斯函数为温度分布形式,应用热平衡积分法及热平衡积分法的细化方法求解单相融化问题,从而得到细化可提高精度的条件。1模型方程研究的定解问题采用半无限大介质中固体融化问题的无量纲形式,初始温度为它的融化温度,数学描述如下622XUTU,0TSX,0T(1)00,XU,0X(2),1,0TU0T(3)0,TXU,TSX,0T(4)DTDSXU,TSX,0T(5)其中10MUUCLSTE,其中STE称为STEFAN数。方程式(1)(5)的解析解为6ERFTXERFTXU21,,0TSX,0T(6)TTS2,0T(7)其中是超越方程12EERF(8)的根,称为融化参数。2高斯函数形式的热平衡积分解方程式(1)两边关于X在,0S内应用热平衡积分法,得到热平衡积分方程00XTSXSXUXUDXTU(9)把式(5)代入上式得00XSXUDTDSDXTU(10)研究高斯函数形式热平衡积分方程近似解,令2,SXCESXBATXV(11)参量CBA,,为常数,由式(11)得TV,0XXV(12)把式(12)计算结果代入式(10)得CCEEBBCDTDSSCC2122(13)由式(11)得2CCTSXCEESBXV(14)由式(14)和式(5)得21CBEDTDSSC(15)由式(3)(4)得1A,0CBEA(16)联立式(13)(15)(16)得CECECCC211221(17)由式(15)(16)得TTCS22212(18)221C(19)当1T,取不同值时,融化界面S的精确解、高斯函数形式热平衡积分解、二次多项式形式热平衡积分解作比较,如表1所示表1精确解,高斯函数形式热平衡积分解和二次多项式形式热平衡积分解比较精确解高斯函数形式二次多项式形式数值解相对误差数值解相对误差02211942069800234221520045204172481707900098179540040906150241494200055155500035008135171347000034139260030311240212372000241273000264121151911508000181180000235由表1可见,高斯函数形式热平衡积分解比二次多项式形式解有更高的精度。3高斯函数形式的热平衡积分细化解将区间,0SN等分,得到N个长度为NS的子区间,10V,0NV,NSIXI令2SXCIIIIESXBAV,IIXXX1,NI,2,120由上式可知未知参数为13N个。(1)据式(3)(4)得11A,NCNNEBA(21)(2)据在节点IX上的连续性得当,1IIXXX时,节点IX上温度近似为2NICIIIIENIBAV当,1IIXXX时,节点IX上温度近似为2111NICIIIIENIBAV所以21211NICINICIIIIIEBEBNIAA,1,2,1NI,(22)(3)同理,据在节点IX上导数的连续性得21212121212NICENICEBBINICINICIIII,1,2,1NI,(23)(4)当SXXN,时21NCNTSXCESBXVN(24)由式(24)和式(5)得21NCNCEBDTDSSN(25)(5)在IIXX,1内对X积分得N个热平衡积分方程11IIIIXXXXXXXVXVDXTV,NI,2,1(26)由式(20)得IXXXV,1IXXXV,TV(27)把式(27)结果代入式(26)得2222221121221211122121212NICIINICINICIINICIINICIINICIIIIIIIIEBCNIEBEBCNIEBNICEBNICEBCDTDSS(28)令式(28)等于式(25)得代数方程21211121212221NICEBNICEBCINICIINICIIII2222121212221NICIINICINICIINICINCNIIIINEBCNIEBEBCNIEBCEB1,2,1NI,(29)(6)在区间SXN,上11NNXXTSXSXXVXVDXTV(30)由式(20)得DXTVSXN1,1NXXXV(31)把式(31)结果和式(5)代入式(30)得第N个热平衡积分方程同(5)做法令第N个热平衡积分方程等于式(25)得代数方程2112122NNCEBCNNNCNNNNNNNCNNCNNNCNCNNCNCNNCEBCEBEBEBCEBNNNNN2112122121122NI(32)由式(21)(22)(23)(25)(29)(32)可求得13N未知参数。4结果与讨论为了比较讨论,取2N为例,图1给出了当1T代入不同的值后高斯函数形式的热平衡积分细化解和热平衡积分解与精确解的误差比较。01020304050607080911112000200400600801012013变量误差热平衡积分细化解热平衡积分解图1近似解与准确解的误差随值变化情况比较由图1可看出,当STEFAN数的倒数小于1时,热平衡积分细化解与精确解的误差明显小于热平衡积分解与精确解的误差即当STEFAN数的倒数小于1时,细化可有效提高精度。参考文献1GOODMANTRTHEHEATBALANCEINTEGRALANDITSAPPLICATIONTOPROBLEMSINVOLVINGACHANGEOFPHASEJTRANSASMEJOURNALOFHEATTRANSFER,1958,803353422LANGFORDDTHEHEATBALANCEINTEGRALMETHODJINTJHEATMASSTRANSFER,1973,16242424283CRANKJTHEMATHEMATICSOFDIFFUSIONMOXFORDCLARENDONPRESS,1964,3103254BELLGEAREFINEMENTOFTHEHEATBALANCEINTEGRALMETHODAPPLIEDTOAMELTINGPROBLEMJINTJHEATMASSTRANSFER,1978,21135713625BELLGESOLIDIFICATIONOFALIQUIDABOUTACYLINDRICALPIPEJINTJHEATMASSTRANSFER,1979,22168116866MOSALLYF,WOODAS,ALFHAIDAONTHECONVERGENCEOFTHEHEATBALANCEINTEGRALMETHODJAPPLIEDMATHEMATICALMODELLING,2005,299039127徐湘田,令锋单相融化问题热平衡积分细化解的收敛性J内蒙古大学学报自然学版,2009,40(2)1281318徐湘田,令锋NEUMANN问题热平衡积分细化解的收敛性J工程数学学报,2010,27(2)1151249MOSALLYF,WOODAS,ALFHAIDAANEXPONENTIALHEATBALANCEINTEGRALMETHODJAPPLIEDMATHEMATICSANDCOMPUTATION,2002,13087100AGAUSSIANFUNCTIONBASEDHEATBALANCEINTEGRALMETHODYUANGUOJING2,1,LINGFENG2(1INNERMONGOLIAUNIVERSITYOFTECHNOLOGY,HOTHOT010051,CHINA2ZHAOQINGUNIVERSITY,ZHAOQING526061,CHINA)ABSTRACTGAUSSIANFUNCTIONISSELECTEDASTHEAPPROXIMATETEMPERATUREPROFILE,ANAPPROXIMATESOLUTIONFORONEPHASEMELTINGPROBLEMISPRESENTEDBYUSINGHEATBALANCEINTEGRALMETHODANDITSREFINEDMETHODITSHOWSTHAT,WHENTHERECIPROCALOFTHESTEFANNUMBER1,THEREFINEMENTOFHEATBALANCEINTEGRALMETHODCANEFFECTIVELYIMPROVEPRECISIONKEYWORDSHEATBALANCEINTEGRALMETHOD,REFINEMENT,GAUSSIANFUNCTION收稿日期20100622基金项目广东省自然科学基金资助项目(04011600)作者简介袁国静(1983),女,山东潍坊人,内蒙古工业大学与肇庆学院联合培养硕士研究生。研究方向相变热传导问题近似解通讯作者令锋(1963),男,陕西岐山人,教授,博士,LINGFZQUEDUCN
编号:201402160857000906    类型:共享资源    大小:264.50KB    格式:DOC    上传时间:2014-02-16
  
10
关 键 词:
数学论文
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
关于本文
本文标题:数学论文.doc.doc
链接地址:http://www.renrendoc.com/p-260906.html

当前资源信息

4.0
 
(2人评价)
浏览:38次
BCEAAEE88296385EDA2815A44814E61A上传于2014-02-16

官方联系方式

客服手机:17625900360   
2:不支持迅雷下载,请使用浏览器下载   
3:不支持QQ浏览器下载,请用其他浏览器   
4:下载后的文档和图纸-无水印   
5:文档经过压缩,下载后原文更清晰   

相关搜索

精品推荐

相关阅读

人人文库
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

网站客服QQ:2846424093    人人文库上传用户QQ群:460291265   

[email protected] 2016-2018  renrendoc.com 网站版权所有   南天在线技术支持

经营许可证编号:苏ICP备12009002号-5