1. 引言
实际问题中,一个微分方程为反映一定的物理意义,往往会含有很小的正参数
,且它们的解和解的导数当
趋近零时出现不连续极限,因此这些方程被称为奇异摄动方程,
称为摄动参数。奇异摄动问题出现于科技和工程的各种分支领域之中,例如河网水质模拟问题、大Péclet数对流热传输问题、运动介质中的电磁场问题、半导体器件模型的漂移扩散方程以及最著名的流体流动问题中的高雷诺数Navier-Stokes方程。因此,奇异摄动问题的数值解法成为计算数学和科学计算学科的一个重要的研究课题。我们考虑如下形式的非守恒型线性奇异摄动两点边值问题:
(1)
其中
是一个充分小的正摄动系数,p和q是充分光滑的函数,且对任意
,
分别满足
,
,且
,
,其中,
都是固定常数,
,
为一类正常数,且在不同处代表不同的值。满足上述条件时,问题(1)只含有唯一解u,且u在
处存在一个指数边界层(见 [1] [2] [3])。
Linß在 [4] 中研究了一类奇异摄动拟线性两点边值问题,提出了一种边界值问题数值处理的自适应方法,并进行了数值试验。Kumar和Srinivasan在 [5] 中利用二阶中心差分格式,提出了一种新的求解对流占优对流扩散奇异摄动问题的自适应网格方法。Kumar和Srinivasan在 [6] 中研究了一种新的自适应网格策略,用于求解具有规则边界层的抛物面对流扩散型奇异摄动问题。该方法利用引入一个熵样变量的辅助不等式也可以作为一个非常有效的适应指标的观点,是对 [5] 中方法的一种拓展。对于拟线性守恒型方程
,
,
,
。郑权,叶福林在 [7] 中采用以弧长函数作为控制函数的基于等分布原理的自适应方法,得到了一阶收敛的误差估计。对于线性非守恒型方程
,
,
,
。杨继明,陈艳萍在 [8] 中采用移动网格方法,得到了一阶收敛的误差估计。对于上述问题陈艳萍在 [9] 中根据近似解的弧长等分布控制函数达到了一阶收敛速度。与以前对完全离散化方法的一些分析不同,该问题不需要此边值问题的保守形式。
本文研究了自适应网格上的一种迎风格式,利用等分布数值解来迭代网格,最终构建成等分布数值解弧长的网格来求解非守恒型线性问题(1),证明了该格式的一致一阶收敛性,并通过数值实验验证了其有效性和误差估计。
2. 稳定性估计
根据
,我们得到
的伴随算子
。伴随算子
的格林函数
满足
(2)
其中
是狄拉克
函数。然后我们可以得到以下形式的解:
(3)
其中
。
通过比较原则我们可以得到格林函数的界:
(4)
定义范数:
我们需要更进一步。对每个
,设辅助函数
满足:
。通过分部积分我们可以得到:
(5)
因为
并且G满足(2)及其边界条件,我们可以得到:
(6)
从(2)中,我们可以得到:
根据 [2] 中§1.1.2,我们设:
(7)
再通过(7)式我们可以得到:
(8)
所以:
。显然,
,因此我们可得:
当
时,上式在
上的积分产生:
所以
,即:
在
上单调减少。
对于
,在
上对(2)式上积分,同理我们得到:
在
上单调增加。
下一步是利用上述单调性和(4)的
界来约束
的
范数。
首先,我们可得:
(9)
由(5)和(9),就有:
因此我们可得:
,即:
(10)
定理1:设
且满足
,
,那么,
证明:由(10)可得:
又由T是线性算子可得:
证毕。
3. 后验误差估计
构造问题(1)的差分格式:
(11)
引入二次插值,对于每个
。
(12)
我们使用以下恒等式:
(13)
由(11)知对于
,
对于任意
,由(1),(12)可知
(14)
积分上面等式
(15)
由(13),(15)知,对于
,
(16)
为
的离散化
的分段函数:
,
。
由定理1,(1),(13)及(16),及模的定义,可知
(17)
由引言知
(18)
由(18),对
进行估计,进而对
进行估计:
(19)
同理由(18),(19)对
进行估计:
(20)
同理由(18),对
进行估计:
(21)
由(18),对
进行估计:
(22)
由(18),对
进行估计:
(23)
由(17),(19),(20),(21),(22),(23)知:
(24)
由
,得:
(25)
由(24),(25)得:
(26)
定理2:设u是问题(1)的准确解,
是问题(1)在差分格式(11)下任意网格上计算所得的数值解,那么,存在常数C,使得:
证明:由无穷模范数定义,及(26)得:
证毕。
4. 自适应网格上的一阶一致收敛性
设
表示
通过外积
联系的伴随算子。根据
因此,
(27)
定义离散格林函数
如下:
且
,
,其中
通过离散比较原则,我们可以得到
,
;
。
对于任何网格函数,我们都可得到以下关系:
(28)
易证,在
时,
是i的单调递增函数,在
时,
是i的单调递减函数。
因此对每个
,我们有:
所以我们得到:
(29)
其中
是离散解
的弧长。
对
有:
(30)
所以,由(29),(30)我们有:
(31)
定理3:设
是(1)在迎风格式(11)下满足弧长等分布的自适应网格
上的解,则:
证明:根据Brouwer不动点定理,存在均匀分布弧长的自适应网格,且对
,从点
到点
的弧长:
由定理2和(31),就得到:
证毕。
5. 数值实验
在BakhvalovShishkin网格和Shishkin网格和自适应网格上比较了迎风格式求解线性非守恒问题(1)的数值结果. 数值收敛速度由下式计算:
另外,
是在迎风格式下BS网格上计算得到的数值解的最大节点误差。
是BS网格上解的数值收敛速度。
是在迎风格式下S网格上计算得到的数值解的最大节点误差。
是S网格上解的数值收敛速度。
是在迎风格式下自适应网格上计算得到的数值解的最大节点误差。
是自适应网格上解的数值收敛速度。
算例1:我们研究下列线性奇异摄动两点边值问题:
该问题的准确解为
其中
,
。
BS网格和S网格和自适应网格上迎风格式的数值结果如表1、表2和表3所示。
Table 1. Maximum errors and the corresponding rate of convergence in example for ε = 10−2
表1. 对于ε = 10−2,算例1中数值解的最大节点误差
Table 2. Maximum errors and the corresponding rate of convergence in example for ε = 10−5
表2. 对于ε = 10−5,算例1中数值解的最大节点误差
Table 3. Maximum errors and the corresponding rate of convergence in example1 for ε = 10−8
表3. 对于ε = 10−8,算例1中数值解的最大节点误差
6. 结论与展望
本文研究了自适应网格上的迎风格式,利用等分布数值解来迭代网格,最终构建成等分布数值解的网格来求解非守恒型线性问题(1),证明了该格式的一致一阶收敛性。是对 [9] 的更一般化的推广,加入了一般项。并且数值实验进行了分别与层适应网格BS网格,S网格的对比,实验结果也说明了本文的自适应网格误差更小,效果更好。但仍有改进的空间,今后我将主要研究通过中点迎风格式,进一步地缩小误差。
致谢
感谢郑权老师的悉心教导。
NOTES
*通讯作者。