论文部分内容阅读
在奇异摄动边值问题的高阶导数项中,扰动参数ε的存在导致解在边界处产生宽度依赖于ε的边界层。用标准有限差分法或有限元法求得的数值解在边界层内一般会产生剧烈的震荡,造成逼近结果在极值点附近失真,尤其是对薄边界层问题。随着奇异摄动问题出现在时下的很多物理与工程建模中,寻求关于薄边界层问题的数值解法也变得尤为重要。 在有限差分法的基础上,Beckett等人提出了基于等距分布网格的有限差分格式,这种等距分布网格具有一阶一致收敛性,而且可以等量地分布边界层和外层内的剖分网点。借助于Chebyshev配置法,Tang Tao等人引入了m-SINE变换,提出了专门求解奇异摄动问题的Chebyshev拟谱法。这一坐标拉伸技术达到了使边界层内有足够多配置点的目的。鉴于奇异摄动边值问题的多尺度特征,本文在等距分布网格和m-SINE函数的基础上,利用区域分割技术提出了用Chebyshev拟谱法分别求解边界层子问题和外层子问题、并在分断点处满足一阶连续的分片C1格式。所做的主要工作是: 本文在第二章构建了适用于含对流项的反应扩散方程的等距分布网格;第三章给出了分片C1格式的具体实现以及分断点、迭代参数m的选取策略;在第四章,通过数值实验,我们用分片C1格式求解不同类型的奇异摄动边值问题,验证该方法的逼近效果,主要优点有:(1)分片C1格式适用于求解同时含对流项的反应扩散方程、具有不同宽度的边界层问题,尤其是薄边界层问题;(2)我们给出了仅通过奇异摄动问题的形式就能先验选取迭代参数m的值的结论;(3)利用等距分布网格剖分点求解分断点,能够更加精确地逼近问题的极值点,从而减小了极值点附近的逼近误差,提高了整体的逼近效果。