Skip to main content
QUICK REVIEW

[论文解读] Second order parameter-uniform convergence for a finite difference method for a singularly perturbed linear parabolic system

Franklin Victor, Joseph Paramasivam Mathiyazhagan|arXiv (Cornell University)|Aug 14, 2010
Differential Equations and Numerical Methods参考文献 5被引用 26
一句话总结

本文提出了一种在Shishkin分段均匀网格上求解具有不同小参数的奇异摄动线性抛物系统有限差分法。该方法在时间方向证明了一阶收敛性,在空间方向证明了二阶收敛性,且对所有摄动参数均保持一致收敛,通过解的精细分解和离散最大值原理分析,实现了空间方向二阶参数一致收敛。

ABSTRACT

A singularly perturbed linear system of second order partial differential equations of parabolic reaction-diffusion type with given initial and boundary conditions is considered. The leading term of each equation is multiplied by a small positive parameter. These singular perturbation parameters are assumed to be distinct. The components of the solution exhibit overlapping layers. Shishkin piecewise-uniform meshes are introduced, which are used in conjunction with a classical finite difference discretisation, to construct a numerical method for solving this problem. It is proved that the numerical approximations obtained with this method are first order convergent in time and essentially second order convergent in the space variable uniformly with respect to all of the parameters.

研究动机与目标

  • 为具有多个不同小参数的奇异摄动抛物系统开发一种鲁棒的数值方法。
  • 解决具有不同摄动参数的线性抛物反应-扩散系统中重叠边界层的问题。
  • 在时间和空间方向均实现参数一致收敛,尤其在空间方向实现二阶收敛。
  • 将已知的标量和两组分情况的结果推广至一般的n组分系统。
  • 证明该方法在参数大小差异显著时仍能保持一致收敛。

提出的方法

  • 在Shishkin分段均匀网格上应用经典的有限差分格式,该网格在边界层附近自适应加密。
  • 将解分解为光滑部分和奇异部分,以分别分析误差。
  • 为离散算子建立离散最大值原理,以确保稳定性和收敛性。
  • 利用对解的光滑部分和奇异部分的精确边界,分析局部截断误差。
  • 使用数学归纳法推导奇异部分的精确估计,尤其在引理7中。
  • 采用比较原理和与网格相关的边界,控制不同区域中离散解的误差。

实验结果

研究问题

  • RQ1在Shishkin网格上,有限差分法能否对具有不同参数的一般n组分奇异摄动抛物系统实现空间方向二阶收敛?
  • RQ2收敛速率是否对所有小摄动参数εi均保持一致有界?
  • RQ3重叠边界层如何影响收敛行为,能否通过分段均匀网格有效处理?
  • RQ4在给定假设下,离散最大值原理能否推广至具有非对角扩散矩阵的系统?
  • RQ5时间和空间方向的确切收敛阶是多少,能否在所有参数值下统一证明?

主要发现

  • 该方法在时间方向实现了对所有摄动参数一致的一阶收敛。
  • 该方法在空间方向实现了对所有参数一致的二阶收敛,包含一个对数因子。
  • 误差界为 ||U - u|| ≤ C N⁻² (ln N)³,证明了参数一致收敛。
  • 分析依赖于通过数学归纳法导出的奇异部分的精确边界。
  • 离散最大值原理和比较原理成功推广至具有非对称但类似M-矩阵结构的系统情形。
  • 由于采用适当的网格分级和参数依赖的估计,即使最小参数εn远小于其他参数,该方法仍能保持一致收敛。

更好的研究,从现在开始

从论文设计到论文写作,大幅缩短您的研究时间。

无需绑定信用卡

本解读由 AI 生成,并经人工编辑审核。