Pversion间断有限元方法求解奇异摄动问题.pdfVIP

第33卷第1期 云南师范大学学报 V01．33No．1 ofYunnanNormal 2013年1月 Journal University Jan．2013 P—version间断有限元方法求解奇异摄动问题+ 陶利娜1， 谢资清1’2 (1．贵州师范大学数学与计算机科学学院，贵州贵阳550001； 2．湖南师范大学数学与计算机科学学院，湖南长沙410081) 摘要： 主要讨论用tversion间断有限元方法(IX3)求解一维奇异摄动反应扩散问题。证明了数值 解的存在唯一性，并给出了相应的数值算例．数值结果表明，与h-versionIX}方法相比较，为了达到同 样数量级的误差，p-version IX}方法的自由度大大下降，并且具有指数收敛性。 关键词：p-version间断有限元；奇异摄动反应扩散问题；指数收敛 中图分类号：0242．21文献标志码：A 文章编号： 1007--9793(2013)01—0039一07 1 引 言 在科学和工程计算中，人们经常遇到带小参数的奇异摄动问题，由于受奇异摄动参数的影响，真解 往往会产生所谓的边界层或内部层现象，即解的梯度在某些狭窄的区域内产生剧烈变化，而边界层或内 部层的存在给数值计算带来了极大的影响。事实上这类问题广泛出现在量子力学、流体力学、弹性力 学、最优控制和化学反应等领域，因此研究其高效的数值方法具有重要的理论意义和应用价值。 本文考虑带Dirichlet边界条件的奇异摄动反应扩散问题： z∈(0，1) re2∥+饿_厂， (1) 【“(0)一Uo，z￡(1)一“1． 其中厂和C(fo)是给定的函数，￡为很小的参数。 由于边界层的存在，如果用传统的有限元方法或差分方法进行计算，往往会产生非物理振荡．当参 数非常小时，为了达到一定的精度，需要强制性地使网格的步长非常小，然而过度的加密网格又会使计 算量剧增，这样的网格令人难以接受．因此，数值求解这类奇异摄动问题极富挑战性。目前对于奇异摄 动问题的边界层计算，存在两种基本的处理方法：第一种是采用局部加密的方法，如等级网格或Shishkin 网格[1]，将这些网格与某些数值方法结合往往能得到拟最优一致误差估计。特别是近年来，人们开始关 注在局部加密网格下用间断有限元方法求解奇异摄动问题及其相应的一致收敛性分析[2_6]，另一种方 DG或者hp-versionDG有限元方法，所谓p-version有限元方法是一种纯粹提高多项式 法是P—version p-version有限元方法的收敛速度是h-version有限元方法的两倍以上。 * 收稿日期：2012—1l一28 基金项目：国家自然科学基金资助项目；贵州省科学技术基金资助项目(LKSF2010]05)． 作者简介：陶利娜(1987一)，女，贵州遵义人，硕士研究生，主要从事有限元方面研究． 通信作者：谢资清(1966一)，女，教授，博士生导师，主要从事间断有限元方面研究． 万方数据 ·40· 云南师范大学学报(自然科学版) 第33卷 受文[4—8]的启发，根据间断有限元方法具有高效、稳定、高精度以及p-version有限元方法具有指 DG 方法的数值格式；第三节证明p