会员注册 | 登录 | 微信快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

   首页 人人文库网 > 资源分类 > DOC文档下载

数学论文.doc.doc

  • 资源星级:
  • 资源大小:264.50KB   全文页数:8页
  • 资源格式: DOC        下载权限:注册会员/VIP会员
您还没有登陆,请先登录。登陆后即可下载此文档。
  合作网站登录: 微信快捷登录 支付宝快捷登录   QQ登录   微博登录
友情提示
2:本站资源不支持迅雷下载,请使用浏览器直接下载(不支持QQ浏览器)
3:本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

数学论文.doc.doc

基于高斯函数形式的热平衡积分法AGaussianFunctionBasedHeatBalanceIntegralMethod袁国静2,1,令锋21.内蒙古工业大学理学院,内蒙古呼和浩特0100512.肇庆学院计算机学院,广东肇庆526061投寄潍坊学院学报收稿日期20100622修改日期2010基金项目广东省自然科学基金资助项目(04011600)作者简介袁国静(1983),女,山东潍坊人,内蒙古工业大学与肇庆学院联合培养硕士研究生。通讯作者令锋(1963),男,陕西岐山人,教授,博士,lingfzqu.edu.cn基于高斯函数形式的热平衡积分法袁国静2,1,令锋2(1.内蒙古工业大学,内蒙古呼和浩特0100512.肇庆学院,广东肇庆526061)摘要采用高斯函数形式的温度分布,应用热平衡积分法及热平衡积分法的细化方法求得单相融化问题的近似解。通过讨论得知,当Stefan数的倒数小于1时,热平衡积分法的细化可有效提高精度。关键词热平衡积分法细化高斯函数中图分类号O241.82文献标志码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精确解,高斯函数形式热平衡积分解和二次多项式形式热平衡积分解比较精确解高斯函数形式二次多项式形式数值解相对误差数值解相对误差0.22.11942.06980.02342.21520.04520.41.72481.70790.00981.79540.04090.61.50241.49420.00551.55500.03500.81.35171.34700.00341.39260.030311.24021.23720.00241.27300.02641.21.15191.15080.00181.18000.0235由表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代入不同的值后高斯函数形式的热平衡积分细化解和热平衡积分解与精确解的误差比较。0.10.20.30.40.50.60.70.80.911.11.200.020.040.060.080.10.120.13变量误差热平衡积分细化解热平衡积分解图1近似解与准确解的误差随值变化情况比较由图1可看出,当Stefan数的倒数小于1时,热平衡积分细化解与精确解的误差明显小于热平衡积分解与精确解的误差.即当Stefan数的倒数小于1时,细化可有效提高精度。参考文献1GoodmanTR.TheheatbalanceintegralanditsapplicationtoproblemsinvolvingachangeofphaseJ.TransASMEJournalofHeatTransfer,1958,80335342.2LangfordD.TheheatbalanceintegralmethodJ.IntJHeatMassTransfer,1973,1624242428.3CrankJ.ThemathematicsofDiffusionM.OxfordClarendonPress,1964,310325.4BellGE.ArefinementoftheheatbalanceintegralmethodappliedtoameltingproblemJ.IntJHeatMassTransfer,1978,2113571362.5BellGE.SolidificationofaliquidaboutacylindricalpipeJ.IntJHeatMassTransfer,1979,2216811686.6MosallyF,WoodAS,AlFhaidA.OntheconvergenceoftheheatbalanceintegralmethodJ.AppliedMathematicalModelling,2005,29903912.7徐湘田,令锋.单相融化问题热平衡积分细化解的收敛性J.内蒙古大学学报自然学版,2009,40(2)128131.8徐湘田,令锋.Neumann问题热平衡积分细化解的收敛性J.工程数学学报,2010,27(2)115124.9MosallyF,WoodAS,AlFhaidA.AnexponentialheatbalanceintegralmethodJ.AppliedMathematicsandComputation,2002,13087100.AGaussianFunctionBasedHeatBalanceIntegralMethodYUANGuojing2,1,LINGFeng2(1.InnerMongoliaUniversityofTechnology,Hothot010051,china2.ZhaoqingUniversity,Zhaoqing526061,china)AbstractGaussianfunctionisselectedastheapproximatetemperatureprofile,anapproximatesolutionforonephasemeltingproblemispresentedbyusingheatbalanceintegralmethodanditsrefinedmethod.Itshowsthat,whenthereciprocalofthestefannumber1,therefinementofheatbalanceintegralmethodcaneffectivelyimproveprecision.Keywordsheatbalanceintegralmethod,refinement,gaussianfunction收稿日期20100622基金项目广东省自然科学基金资助项目(04011600)作者简介袁国静(1983),女,山东潍坊人,内蒙古工业大学与肇庆学院联合培养硕士研究生。研究方向相变热传导问题近似解通讯作者令锋(1963),男,陕西岐山人,教授,博士,lingfzqu.edu.cn

注意事项

本文(数学论文.doc.doc)为本站会员()主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网(发送邮件至[email protected]或直接QQ联系客服),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。

[email protected] 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5