非定常扩散方程的虚边界元解法.pdfVIP

第 26卷 第 6期 北京工商大学学报 (自然科学版) Vo1．26No．6 2008年 11月 JournalofBeijingTechnologyandBusinessUniversity(NaturalScienceEdition) Nov．2008 73 文章编号：1671—1513(2008)06—0073．04 非定常扩散方程的虚边界元解法 汪学海 ， 吴胤林2 (1．平顶山工学院 数理系，河南 平顶山 467000；2．重庆大学 数理学院，重庆 400030) 摘 要：对于二维齐次和非齐次非定常扩散方程问题，利用与时间有关的基本解，基于单层位势的 延拓，建立虚边界积分方程，然后用虚边界元法求解．最后，给出了数值算例验证 了虚边界元法求 解非定常扩散方程问题的可行性和有效性 ． 关键词 ：非定常扩散方程 ；虚边界元；单层位势 中图分类号 ：O241．82 文献标识码 ：A 非定常扩散方程是一类典型、简单但是常用的 其中 是R 中的单连通域，I1：an是光滑 的，边 抛物型偏微分方程，应用广泛，存在有很多区域型和 界 r=r】+-／2．时间域为 工：￡0≤t≤T+。。，系 边界型的数值求解方法_1。]．用边界元求解抛物型 数 k对不同的物理方程有不同的解释．若 f(，t)= 方程存在多种方法，传统的方法涉及到对时间和空 0，(1)式表示的就是齐次扩散方程 的混合初边值 问 间变量的积分，其中包含有指数积分函数，有奇异积 题 ． 分． 在边界r外的一定距离处选取边界r ，称其为 虚边界元算法在求解区域之外假想的虚边界上 虚边界．并且假定在虚边界上分布有与时间有关的 存在虚拟密度函数，通过 已知边界条件推出这个虚 虚拟密度函数 q(，r)，而且 在穿越边界时是连 拟源后再计算原问题的解．由于场 点和源点分别位 续的．利用扩散方程 问题的与时间有关的基本解， 于虚、实边界上，避免了传统边界元方法在边界以及 通过格林公式_6j，即得到虚边界积分方程为： 靠近边界 的地方存在 的积分奇异性 问题 ．同时，该 算法还保持了传统边界元方法所具有的将问题降低 (，t)=I 0()“(，t；，t0) + Jn 一 维、计算精度高等特点．通常的虚边界元方法只 见于应用到椭圆型 问题_5J，本文将其推广到抛物型 愚』q，()“(；)r+ 问题的求解，并以非定常扩散方程为例，讨论了具体 jlIjInf(，r) (，；，r)d~dr 的数值求解过程． 1 齐次和非齐次非定常扩散方程的虚 = ．n『“。()旦 + 边界元积分方程 志j’，口(，r)旦—=。 sr+ 考虑如下二维非齐次非定常扩散方程的混合初 边值问题 ： jIjn 0 H ． z f =Au(，)+(，)(，)E21×J