流线扩散有限元方法在分层网格上收敛性分析.pdfVIP

2015年 2月 计 算 数 学 第 37卷第 1期 Feb．，2015 MATHEMATICANUMERICA SINICA Vol，37．No．1 流线扩散有限元方法在分层网格上的 收敛性分 1) 尹云辉 祝 鹏 (嘉兴学院数理与信息工程学院浙江嘉兴314001) 杨宇博 (嘉兴学院南湖学院，浙江嘉兴314001) 摘 要 本文在分层网格上分析了采用线性元的流线扩散有限元方法求解一维对流扩散型奇异摄动问 题的一致收敛性．在E N 的前提下，可以证明在 印 范数下的一致误差估计为O(N (1og1)) 在数值算例部分对理论结果进行了验证． 关键词：奇异摄动问题；流线扩散有限元；分层网格 MR (2000)主题分类：65M12，65N12 1．引 言 奇异摄动问题广泛存在于物理工程，流体力学，化学力学等应用领域，这类问题的解常常 出现局部剧烈变化，用通常的有限差分法和有限元方法等数值方法求解时往往在急剧变化地 带产生非物理振荡，得不到令人满意的结果．即通常的数值格式对小参数 问题是不稳定的． 1979年 Hughes和 Brooks将标准 Galerkin有限元法和最小二乘方法相结合，提出流线 扩散有限元 (SDFEM)|1]j此方法是具有稳定性强、精度高的特点。流线扩散有限元方法对所 有的扰动参数都是稳定的，对抑制非物理振荡有很好的效果．为了获得精确的计算结果，通常 采用强稳定性的数值方法 (诸如 SDFEM)在层自适应网格上求解奇异摄动问题．文献 [2】在 Shishkin网格上给出了无穷范数下的误差估计；文献 [3】给出了在流线扩散范数 (SD范数)下 的误差分析，但是由于 ∑ lIb·V(u—UI)lI。r1无法估计，只分析了ll 一钆lIsD的误差；为 j Cn0 了解决这个 问题，文献 [4,5】在离散的SD范数下分析了 一Uh的误差． 2006年DurAn首次提出了分层网格_6J．该网格的优点是相邻单元的尺寸变化平缓，不用 考虑粗细网格的转变点，计算效果也不再依赖于网格构造中出现的参数．本文试图在分层网 格上讨论 ∑ lIb·(一UI) 。f1，并给出流线扩散有限元方法在 sD范数下的一致收敛性 K C 【22 分析． 本文研究的是一维对流扩散型的奇异摄动问题，其方程为 l{uLu，(n0)、：z一+6，+c( ， in(0)1) (1．1) = (1)=0 、 2014年 3月 19日收到 1)基金项 目：浙江省 自然科学基金 (LQ12A01014)和浙江省教育厅科研项 目(Y201431793)资助 计 算 数 学 其中￡为一任意小的正常数 (设E N-1)，6()，c()，f(x)均为光滑函数，且b(x)一 ， 0， c(x)0．在这些假设条件下，(1．1)存在唯一解，并且解通常在流入边界处X=0有一个宽度 为0(slog )的边界层． 2．分层网格 将 2【=[0，1】区域分为细网格部分 2【1=[0，C1Clog]和粗网格部分Q2=[CIClog1，1] 其中Cl 吉，如图l所示： 0q 1。g l 图 1区域划分 给定网格参数h0，区域 [0，1]上网格剖分 {)0定义如下：