【2017年整理】数学论文.doc

基于高斯函数形式的热平衡积分法 A Gaussian Function Based Heat Balance Integral Method 袁国静，令锋 1.内蒙古工业大学 理学院，内蒙古 呼和浩特 010051； 2.肇庆学院 计算机学院，广东 肇庆 526061 投寄：潍坊学院学报 收稿日期：2010-06-22； 基金项目：广东省自然科学基金资助项目 作者简介：袁国静（1983-），女，山东潍坊人，内蒙古工业大学与肇庆学院联合培养硕士研究生。 通讯作者：令锋（1963-），男，陕西岐山人，教授，博士，lingf@ 基于高斯函数形式的热平衡积分法 袁国静，令锋 （1.内蒙古工业大学，内蒙古 呼和浩特 010051；2.肇庆学院，广东 肇庆 526061） 摘 要：采用高斯函数形式的温度分布，应用热平衡积分法及热平衡积分法的细化方法求得单相融化问题的近似解。 通过讨论得知， 当数的倒数小于1时，热平衡积分法的细化可有效提高精度。 关 键 词：热平衡积分法；细化；高斯函数 中图分类号 ：O241.82 文献标志码：A 0 引言 热平衡积分法（HBIM）是1958年由Goodman最早提出的一类可以用来生成由微分 方程控制的热传导方程的近似解的半解析的有效方法。 但是，HBIM的计算结果对选取的温度分布函数具有敏感性，Goodman最先提出选用二次多项式作为温度分布函数，同时也简单提及三次函数的情况，有算例表明，选择三次函数作为温度分布函数的精度却比二次多项式的结果更差，为了减少计算结果对温度分布函数的依赖性，提高计算结果精度，Noble提出了HBIM的细化（Refinement）方法。 Bell据此法求解了伴有相变的一维问题，获得较好的结果。 Mosally等人以单相融化问题为模型，采用与数值实验结果相对比的方法，给出了当数等于1时单相融化问题HBIM细化解的收敛速度。徐湘田和令锋通过理论分析而不依赖数值试验，分别给出了单相融化问题和Neumann问题的热平衡积分细化解收敛性的证明。 Gauss函数在 HYPERLINK /zh-hans/%E5%9F%83%E5%B0%94%E7%B1%B3%E7%89%B9%E5%A4%9A%E9%A1%B9%E5%BC%8F \o 埃尔米特多项式 Hermite多项式的定义中起着重要作用。热传导问题导热过程的特征促使人们设想指数函数可能是适当的温度分布函数的选取形式，但计算表明，选择指数函数作为温度分布函数的精度反而比选取二次多项式的结果要差，所以，选取高斯函数作为温度分布函数可能是一种有益的尝试。 本文以单相融化问题为模型，以高斯函数为温度分布形式，应用热平衡积分法及热平衡积分法的细化方法求解单相融化问题，从而得到细化可提高精度的条件。 1 模型方程 研究的定解问题采用半无限大介质中固体融化问题的无量纲形式，初始温度为它的融化温度，数学描述如下： ，， （1） ， （2） （3） ，， （4） ，， （5） 其中，其中称为数。 方程式（1）—（5）的解析解为 ，， （6） ， （7） 其中是超越方程 （8） 的根，称为融化参数。 2 高斯函数形式的热平衡积分解 方程式（1）两边关于在内应用热平衡积分法，得到热平衡积分方程 （9） 把式（5）代入上式得 （10）研究高斯函数形式热平衡积分方程近似解，令