奇异摄动问题在Bakhvalov-Shishkin网格上流线扩散有限元逼近.pdfVIP

2013年 11月 计 算 数 学 第 35卷第 4期 NOV．，2013 MATHEMATICA NUMERICA SINICA Vo1．35，No．4 奇异摄动问题在 Bakhvalov-Shishkin网格 上的流线扩散有限元逼近木1) 尹云辉 祝 鹏 (嘉兴学院数理与信息工程学院，浙江嘉兴314001) 杨宇博 (嘉兴学院南湖学院，浙江嘉兴314001) 摘 要 本文采用线性插值的流线扩散有限元在 Bakhvalov—Shishkin网格上求解一维对流扩散型的 奇异摄动问题．在 E N 的前提下，可以得到，关于扰动参数 E是一致收敛的．在离散的 SD 范数下，其 11,一Ul的误差阶提高到N_。，u—Uh的误差阶达到 Ⅳ (1nⅣ) ．最后，通过数值 算例，验证了理论分析． 关键词：奇异摄动问题；流线扩散有限元；Bakhvalov—Shishkin网格 MR (2000)主题分类：65M12，65N12 1．引 言 奇异摄动问题是一类广泛存在于物理工程，流体力学，化学力学等研究领域的重要问题， 这类 问题的解常常出现局部剧烈变化，如含有边界层或内层等．这类 问题用通常的数值方法， 如用有限差分法和有限元方法求解时往往在急剧变化地带产生非物理振荡，得不到令人满意 的结果．即通常的数值格式对小参数问题是不稳定的．人们感兴趣 的是对所有的扰动参数都 有效的方法，也就通常所说的一致收敛的格式． 1979年Hughes和Brooks将标准Galerkin有限元法和最小二乘方法相结合，提出流线扩 散有限元 (SDFEM)LlJ，此方法是通过对检验函数的适当选取，施加了一个沿流场方向的人工 粘性项，具有稳定性强、精度高的特点．流线扩散有限元方法对所有的扰动参数都是稳定的， 对抑制非物理振荡有很好的效果．目前，在文献 [2,3，4】给出了流线扩散有限元在Shishkin网格 上的一致性分析，其误差阶为 Ⅳ (InⅣ)0．本文将给出流线扩散有限元在 Bakhvalov—Shishkin 网格上的误差分析． 本文研究的是一维对流扩散型的奇异摄动问题，其方程为 ! -- EU I{~,(：，+b()扎+c()=， in0(．1) (1．1) nh=-- u (1)= 0， 其中 ￡为一任意小的正常数 (设E N-1)，6(z)，c()，f(x)均为光滑函数，且 b(x) 0， c(x)0．在这些假设条件下，(1．1)存在唯一解，并且解通常在流出边界处 =1有一个宽度 为 O(~InI)的边界层． 2012年 12月24日收到． )基金项 目：浙江省 自然科学基金 (LQ12A01014)以及嘉兴学院科研启动基金 资助 366 计 算 数 学 2013住 2．Bakhvalov—Shishkin(B—S)网格 网格参数=min(1，1nJ7v)，将Q=0【，1】区间分为粗网格部分2【=0[，1一和细 网格部分 2【2=[1一 ，1]，如图1所示： 0 l一 l 图 1Ⅸ问划分 粗网格区间 2【划分为 N等分，粗网格步长 H= ，则网格点为