改进的奇异摄动自适应移动网格方法.docVIP

改进的奇异摄动自适应移动网格方法 　　摘要：用等分弧长函数来控制网格剖分，用迎风有限差分格式来求解一类奇异摄动两点边值问题的自适应算法。本文用了的数值试验证明了算法的可行性和高效性。 　　关键词:奇异摄动；自适应网格；迎风有限差分格式；等分原则 　　中图分类号：O241文献标识码：A文章编号：1009-3044(2008)06-10ppp-0c 　　 　　A Improved Adaptive Grid Method for a Singularly Perturbed Problem 　　LI Li 　　(Information department of Hunan business college, Changsha 410205,China) 　　Abstract: In this paper, the discrete solution are generated by an upwind finite difference scheme and the grid is formed by equidistributing a monitor function based on arc-length. A improved numerical experiment proved that the algorithm is feasible and efficient. 　　Key words: singular perturbation problems; adaptive mesh; upwind finite difference; equidistribution principle. 　　 　　1 引言 　　 　　近年来，研究带边界层对流占优的对流扩散问题移动网格算法收敛性分析问题的网格构造主要有两种：特殊网格（如B-type meshes, S-type meshes）、自适应网格。对比之下，在自适应网格方面的工作甚少。 　　我们考虑如下奇异摄动的两点边值问题的数值近似解法： 　　 　　这里0ε1,ε是一个很小的正的参数。本文考虑的是更特殊的情形，即p(x)=1。 　　对于(1)形式的对流扩散方程问题，我们已有一些自适应网格方面的结论。这些都是从控制函数的等分上所得到的结果[1,4,5,6]。Qiu, Sloan 和Tang[6]研究基于半离散方法的收敛速度问题，这种方法指出精确解被应用于控制函数中。这就使得分析变得简化，并且能够在解区间上给出一个清晰的网格剖分结构。他们证明了对于任意给定的γ∈[0,1]，总存在一个不依赖于ε和N的正常数C(γ)，使得 　　 　　其中N是剖分节点的总数，并且足够大，ll03.tif 是数值近似值。本文的主要目的是用一种改进的数值试验方法证明利用o(N-1lnN)替代（2）式的右端从而提高结论（2）的一致收敛性的可行性。其改进之处在于剖分节点构造上的工作有了大大的简化。 　　 　　2 问题描述 　　 　　我们考虑的是一种特殊的奇异摄动两点边值问题的数值近似解法。方程如下所示： 　　 　　其中ε是一个很小的正数，且0ε1。令[0,1]上一个任意的剖分：ΩN={xj|x0x1…xN=1}，在ΩN上我们有如下定义： 　　 　　其中所用到的算子如下面定义： 　　 　　在 [0,1]上，数值网格由等分弧长函数M(x)=??1+(u(x))2给出，并且产生一个映射x=x(ξ)： 　　 　　其中L是u在[0,1]上的弧长。这种方法叫做半离散方法。在这篇文章中我们运用数学试验讨论这种数值方法的误差分析等。 　　在论文中，C代表一个普通的独立于ε以及网格剖分的常数，并且在不同的地方可以取不同的值。这些工作的主要结果由以下定理给出[7]。 　　定理1令u(x)为（3）的一个精确解，且令ll08.tif为由（6）所定义的网格上的有限差分格式（4）和（5）所得到的，那么存在一个不依赖于ε和N的正常数C，使得：|u(xi)-uNi|≤CN-1lnN，0iN 　　 　　3 数值试验 　　 　　3.1 求真解表达式? 　　微分方程 　　 　　的真解表达式为：ll10.tif 　　 　　3.2 用迎风格式进行离散? 　　我们运应迎风格式得到： 　　 　　此时系数矩阵A以及矩阵B和U分别为： 　　 　　那么方程（9）化为： 　　AU=B， 　　3.3 构造剖分节点 　　要求出方程的解，关键在于构造出剖分节点，从而得到hi。它不同于我们平时一般所采取的方法：平均剖分。它是由一个控制函数（6）来给出。其中ξ是计算坐标，x是物理坐标。这个式子所代表的意义就是：把计算坐标上等分弧长的点映射到物理坐标上，而得到的这些物理