×

用迎风差分格式分析奇异摄动微分方程的稳定性。 (英语) Zbl 1312.65174号

作者将逆风差分格式应用于以下具有边界条件的奇摄动微分方程(SPDE):\[-\varepsilon\Delta u+\alpha u_x+\beta u_y+cu=f\,\,\text{in}\,S,\quad u=g\,\,\text{on}\,\Gamma,\]其中,\(Delta=\frac{\partial^2}{\parial x^2}+\frac}\partial y^2}),\(S\)是矩形:\(S=\{(x,y)|0<x<\pi\},\,\Gamma=\partial-S\)为其边界,参数\(\varepsilon\)非常小,即,}\(\valepsilon\ll 1\)。
本文利用线性算子的条件数分析了迎风差分格式的数值稳定性。最初,它们从使用微分算子及其共轭算子的特征值计算的传统条件数开始。但是,如果参数\(\varepsilon\leq 10^{-8}\),那么传统的条件数将是巨大的,这将给SPDE问题方案的稳定性分析带来问题。然后,作者引入了两个新的条件数{即}有效条件数和实际条件数,这两个条件数可以用极大值原理进行计算。上述两个条件数都提供了由舍入误差、截断误差和离散化误差引起的解误差的较小界。
很明显,数值格式的稳定性分析与线性算子的条件数有关。在这里,作者首先使用微分算子及其共轭算子的特征值来计算传统的条件数,然后发现对于(varepsilon \leq 10^{-8}),得到的条件数很大。然后,他们利用最大值原理计算有效条件数,该原理取决于非齐次项和边界条件。实际条件数由离散最大值原理计算。
在比较了所有三个条件数之后,他们发现有效条件数和实际条件数的最坏情况也小于线性算子的传统条件数。作者还证明了当离散极大值原理存在时,有效条件数和实际条件数都能保证格式的良好调节性。根据有效条件编号和实际条件编号,即使对于\(\varepsilon\rightarrow 0\),也不满足方案的ill-conditioning。这澄清了巨大的传统条件数所造成的困境。

MSC公司:

65N12号 含偏微分方程边值问题数值方法的稳定性和收敛性
65号06 含偏微分方程边值问题的有限差分方法
35B25型 偏微分方程背景下的奇异摄动
35J25型 二阶椭圆方程的边值问题
65层35 矩阵范数、条件、缩放的数值计算
35磅50英寸 PDE背景下的最大原则
65奈拉 涉及偏微分方程的边值问题的误差界
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] J.T.H.Miller、E.O’Riordian和G.I.Shishkin,奇异摄动问题的拟合数值方法,世界科学,新加坡,1996年·Zbl 0915.65097号
[2] H.G.Roos、M.Stynes和T.Tobiska,奇异摄动微分方程的数值方法,对流扩散和流动问题,Springer,柏林,1991年。
[3] S.Wang和Z.C.Li,奇异摄动对流扩散方程的有限元和体积方法与各向异性网格加密的非协调组合,数学计算72(2003),1689-1709·Zbl 1024.65109号
[4] F.Mazzia和D.Trigante,二阶奇异摄动问题的数值方法,计算数学应用23(1992),81-89·Zbl 0757.65094号
[5] H.O.Kreiss、N.K.Nichols和D.L.Brown,两点边值问题的数值方法,SIAM J Numer Anal23(1986),325-368·Zbl 0608.65049号
[6] 李振中,黄洪涛,数值偏微分方程的有效条件数,数值线性代数应用15(2008),575-594·Zbl 1212.65410号
[7] 李振中,黄海涛,陈建涛,魏耀伟,有效条件数及其应用,计算89(2010),87-112·Zbl 1198.65085号
[8] Z.C.Li、H.T.Huang、Y.Wei和A.H.‐D。程,数值偏微分方程的有效条件数,科学出版社,北京,2013。
[9] J.R.Rice,《矩阵计算和数学软件》,麦格劳-希尔图书公司,纽约,1981年·Zbl 0459.65008号
[10] T.F.Chan和D.E.Foulser,有效良好条件线性系统,SIAM J Stat Compute9(1988)963-969·兹伯利0664.65041
[11] S.Christiansen和P.C.Hansen,有效条件数应用于某些边界方法的误差分析,J Comput Appl Math54(1994),15-36·Zbl 0834.65033号
[12] J.M.Banoczi、M.Chiu、G.E.Cho和C.F.Ipsen,《右手边对线性系统解精度影响的缺乏》,SIAM科学计算杂志20(1998),203-227·Zbl 0914.65047号
[13] J.T.Chen、S.R.Lin和K.H.Chen,用边界元法求解拉普拉斯方程时的退化尺度问题及其处理,国际数值方法工程62(2005),233-261·Zbl 1179.74161号
[14] 李振聪,黄洪涛,使用局部网格细化的有限元方法的有效条件数,应用数值数学59(2009),1779-1795·Zbl 1180.65150号
[15] M.Krizek和P.Neitaanmaki,《关于线性三角形元素梯度的全局超收敛性》,《计算应用数学杂志》18(1987),221-233·Zbl 0602.65084号
[16] 李振聪,黄洪涛,黄春生,卢涛,方奇,涉及边界奇异性的泊松方程有限差分方法的有效条件数,数值函数分析优化21(2011),659-681·Zbl 1416.65404号
[17] T.Yamamoto,Q.Fang和X.Chen,Dirichlet问题Shortley‐Weller近似的超收敛和非超收敛,数值函数分析优化22(2001),455-470·Zbl 0991.65108号
[18] J.M.Ortega,《数值分析,第二课程》,SIAM出版社,费城,1970年。
[19] P.G.Ciarlet,《椭圆问题的基本误差估计》,有限元方法(第一部分),P.G..Ciarlet和(编辑)J.L.Lions(编辑),编辑,荷兰北部,阿姆斯特丹,1991年,第17-352页·Zbl 0712.65091号
[20] G.H.Golub和C.F.van Loan,《矩阵计算》,第二版,约翰霍普金斯大学出版社,巴尔的摩和伦敦,1989年·Zbl 0733.65016号
[21] J.H.Wilkinson,《代数特征值问题》,克拉伦登出版社,牛津,1965191·Zbl 0258.65037号
[22] F.Cucker、H.Diao和Y.Wei,《关于Moore-Penrose逆和线性最小二乘问题的混合和分量条件数》,《数学计算》76(2007),947-963·Zbl 1115.15004号
[23] Z.C.Li、H.S.Tsai、S.Wang和J.J.H.Miller,二维边界层奇摄动微分方程的精确和近似解析解,计算数学应用55(2008),2602-2622·Zbl 1142.65385号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。