用泡函数法求解抛物型强对流-扩散方程.pdfVIP

维普资讯 第7卷 第6期 2007年3月 科 学 技 术 与 工 程 Vo1．7 No．6 Mar．2007 1671-1819(2007)6-0959—05 ScienceTechnologya·ndEngineering @ 2007 Sci．Tech．Engng． 用泡函数法求解抛物型强对流-扩散方程 李文涛 陈焕贞 (山东师范大学教学科学学院，济南250014) 摘 要 利用标准的Galerkin方法解决强对流扩散问题时，格式稳定性差并伴有强烈的数值震荡现象。为克服上述缺陷，人 们提出了借助于泡函数构造数值方法以得到稳定化的数值解 ，并成功地应用于求解椭圆型强对流扩散问题。文 中利用泡函 数构造数值 ，求解抛物型强对流扩散问题的稳定化方法，并给 出收敛性分析。 关键词 泡函数法 抛物型强对流扩散问题 收敛性 中图法分类号 0242．2 文献标识码 A ． 我们知道，强对流扩散方程在工程上被广泛地 角剖分或矩形剖分，为剖分单元 c (力)为分 应用，有着重要的现实意义；然而，用标准的Galerkin 片多项式空间，当 为三角剖分时，不妨记 = 方法在解决该问题时，求解过程中出现运算量过大， ：={∈HO(力)：l，∈P()}， 为矩形剖分时， 计算误差累积过快，最终结果不稳定等问题。20世 记 = ：={∈no(力)：l，∈Q ()}，其中，P 纪90年代出现一种名为泡函数法的新方法，它不但 ()和Q (T)分别表示三角单元和矩形单元 上的 具有 良好的逼近性，而且具有很好的稳定性。 k次多项式空间。不同的是在每个单元 上构造 考虑如下椭圆型强对流扩散问题 ： ()的一个有限维子空间 ，，并令 =Q ，为泡 ㈩ 函数空间，记 = + ，则问题 (1)泡函数方法可 定义为：求 ^= ￡+B∈ ，满足 其中，Lu=一eAu+口· ，力∈R 为有界多角形区 0(^，^)=(f，^)，V ^∈ (2) 域，tl,=(tl,】，tl,2，…，tl,)满足div0=0，E为常数，并且 有 ∈《l0l·diam21厂∈L(∈．)。 (2)式中，0(，)=∈IIn Vvdx+IIn(0· )vdx。 由于 ∈《l0l·diam21，问题 (1)表现为强对流 占优，利用标准的Galerkin方法求解该问题时，边界 由文献 [2]，若令N，=dim{ ∈ }，{i，i 层会导致近似解出现震荡，数值格式不稳定。为克 = 1，2，…，Ⅳ，}表示 在单元 上的一组基，则 服标准Galerkin方法的数值不稳定性，文献 [1]提出 Br=span{br．1，br．2，…，b ， }，其中b ∈Ho 了利用泡函数构造稳定数值方法的思路，即在标准 ()， ∈Ho()分别满足： Galerkin方法的基础上，通过在每个小单元 上引 入 内非0，OT和 外恒为0的泡函数，来增加离散 f =一 ，in， (3) b【 =0，onOT 空间 基的个数 ，从而解决近似解的震荡问题。 同标准Galerkin方法相同，令 表示力上的三 f珥=，，in (4) ‘ 【=0，onOT 2006