资源描述：
《边界条件随时间变化的热传导问题的一种近似解法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、边界条件随时间变化的热传导问题的一种近似解法第33卷第2期Vo1.33No.2菏泽学院JournalofHezeUniversity2011年3月Mar.20ll文章编号:1673—2103(2011)02—0048—05边界条件随时间变化的热传导问题的一种近似解法袁国静一,令锋(1.内蒙古212业大学理学院,内蒙古呼和浩特010051;2.肇庆学院计算机学院,广东肇庆526061)摘要:在热平衡积分法和改进的热平衡积分法的基础上,提出了一种求解热传导问题的新方法.应用此方法求解边界条件随时间变化的热传导问题,结果表明:此方法
2、计算过程简单且计算精度高.关键词:热传导;固定温度边界条件;线性温度边界条件中图分类号:0241.82文献标志码:A引言热平衡积分法(HeatBalanceIntegralMethod;HBIM)是Goodman…提出的解决热传导问题的一种解析方法.即通过对导热微分方程在f时刻的摄动深度(￡)关于空问变量积分得到一个积分方程,选取满足边界条件的函数代人积分方程,求得整个区域的温度分布.用简单函数来近似描述温度场,计算难度低于分析解法,计算速度快于数值解法,且计算精度较高,这使其在导热问题的求解中被广泛应用.改进的热平衡积分法(
3、RefinedIntegralMethod;RIM)是Sadoun等人提出的一种半解析方法.即通过对导热微分方程关于空间变量积分2次并结合传统的热平衡积分方程得到一个新的积分方程,选取满足边界条件的函数代入这个积1叮'分方程,求得整个区域的温度分布.该方法建立在前人提到的改进积分基础上J,避免了在=0处产生ox的误差,在一定Stefan数范围内该方法比其他解析方法显示出更高的准确性2J.近似函数的选择是应用热平衡积分法及改进的热平衡积分法需要考虑的关键问题,大多数学者采用的形式为多项式函数和指数函数一J,Goodman最先提出
4、选用二次函数,同时也提及三次函数的情况J,Antic和Hill在2003年研究谷物储藏的扩散问题中采用的三次函数.Lang~rd介绍了测量误差函数E(￡)=—,,T'2,2f,I警一ldx,使E(f)最小化得到n的值.在相同的条件下由测量误差函数确定,t值得到的解较先前n…ua值得到的解显示出更高的精确性J,但是求解n值的代数计算过程较繁琐,首先须对假设的近似温度分布函数求得对t和的偏导数,然后代人E(t)取平方后在区间[0,]上积分求得关于6,6和凡的复杂的表达式,把通过热平衡积分法和改进的热平衡积分法求得的摄动深度代人求得
5、的表达式,最后对(f)求最小化得到n的值.本文在热平衡积分法和改进的热平衡积分法的基础上,提出了一种求解热传导问题的新方法,应用此方法求解了边界条件随时间变化的热传导问题,结果表明:此方法计算精度高且计算过程简单.48收稿日期:2011—1—19基金项目:广东省自然科学基金资助项目(04011600)作者简介:袁国静(1983一),女,山东潍坊人,在渎硕士研究生,研究方向:相变热传导问题.令锋(1963一),男,陕西岐山人,教授,博士,硕士研究生导师,研究方向:冰川冻土2011年袁国静,等:边界条件随时间变化的热传导问题的一种
6、近似解法第2期1模型方程边界条件随时间变化热传导问题的半无限大区域模型方程如下:～OT:磐0<<∞afaT=0T:(t)塑一0a对任意西(t)解析解可表示为:t=0=0—}∞(1)(2)(3)(4)():f业孚dA(5)21T如(t—A)丁由于对随时间变化的边界条件咖(t),式(5)的求解过程繁琐甚至很多情况下的结果无法求得,故考虑求解过程相对简单的近似解.2近似解法介绍为保证近似解与温度场光滑融合,近似温度分布函数需满足T(6,t)=0,(6,t)=0,假设温度近似函数为:=(￡)(1一(6)2.1热平衡积分法将
7、式(1)关于空间变量积分,有:~O—T,dx=Ox(7)—,,整理得:一乳.㈩将式(6)代人式(8)得:dt【1】=6(9)【n+J2.2改进的热平衡积分法将式(1)关于空间变量积分两次,有:(=(一:.)整理得:6r一广:一f一6一OTI(11)JoOt』.OtJ:oc3xJ:0将式(1)代人式(11)左边第一项,整理得:UoTd:r(o,)(12)将式(6)代人式(12)得:201I丘菏泽学院第2期f】=(13i2)d【(n十)(n+)J丫将不同的温度边界条件代人式(9),式(13)求得摄动深度,令两式求得的结果相等得/7
8、,值.3边界条件随时间变化条件下的近似解Sadoun在文献[2]中研究了固定温度边界条件和线性温度边界条件两种情形,本文以这两种情形为例研究热传导问题.3.1固定温度边界条件(t)=I将条件(t)=l代人式(13)并积分得:=(14)式(9)中假定n是常数,将条件()=l代人
