Skip to main content
QUICK REVIEW

[论文解读] A parameter uniform fitted mesh method for a weakly coupled system of two singularly perturbed convection-diffusion equations

Saravana Sankar Kalaiselvan, John J. H. Miller|arXiv (Cornell University)|Sep 24, 2018
Differential Equations and Numerical Methods参考文献 7被引用 5
一句话总结

本文提出了一种参数一致的拟合网格方法,用于求解具有不同扰动参数的弱耦合二维奇异摄动对流-扩散方程组。在分段均匀的Shishkin网格上采用迎风有限差分格式,该方法在最大范数下对两个扰动参数均实现了第一阶收敛,且收敛性与参数无关。该方法得到了严格的误差估计分析和数值实验的支持,展示了在不同参数取值下均具有鲁棒的收敛性。

ABSTRACT

In this paper, a boundary value problem for a singularly perturbed linear system of two second order ordinary differential equations of convection- diffusion type is considered on the interval [0, 1]. The components of the solution of this system exhibit boundary layers at 0. A numerical method composed of an upwind finite difference scheme applied on a piecewise uniform Shishkin mesh is suggested to solve the problem. The method is proved to be first order convergent in the maximum norm uniformly in the perturbation parameters. Numerical examples are provided in support of the theory.

研究动机与目标

  • 为求解具有不同小扰动参数的两个奇异摄动对流-扩散方程组,开发一种鲁棒的数值方法。
  • 解决由于弱耦合和不同扰动参数导致解分量中边界层宽度不等的挑战。
  • 确保收敛性对两个扰动参数均一致,尤其在层丰富的解存在时。
  • 为解的分解和导数估计提供严格的理论框架,相较于先前工作采用更精确的定义,实现改进。

提出的方法

  • 在分段均匀的Shishkin网格上应用迎风有限差分格式,对两个二阶奇异摄动对流-扩散方程组进行离散化。
  • 构建Shishkin网格以在x=0处解析边界层,网格的密铺程度根据较小的扰动参数ε₁和较大的ε₂进行调整。
  • 将解分解为正则部分和层部分,利用最大值原理和稳定性估计分析层函数。
  • 为离散问题建立了最大值原理,确保数值格式的正性和稳定性。
  • 证明了该方法在最大范数下对ε₁和ε₂均以第一阶速率一致收敛。
  • 通过在基准问题上进行数值实验,验证了理论收敛性,实验中变化了ε₁和ε₂的取值。

实验结果

研究问题

  • RQ1在Shishkin网格上采用拟合有限差分格式,能否实现对具有不同扰动参数的弱耦合二维奇异摄动对流-扩散方程组的参数一致收敛?
  • RQ2在存在耦合的情况下,解分量中的边界层结构如何依赖于各自的扰动参数ε₁和ε₂?
  • RQ3从先前工作中移除限制性假设ε₁ ≤ C N⁻¹ 和 ε₂ ≤ C N⁻¹ 对收敛性分析和方法设计有何影响?
  • RQ4如何改进解的分解方法,以更准确地捕捉每个扰动参数对层行为的独立贡献?
  • RQ5实际观测到的收敛速率和误差常数是多少?当ε₁和ε₂趋于零时,这些量是否保持稳定?

主要发现

  • 所提出的在Shishkin网格上的迎风有限差分格式在最大范数下对两个扰动参数ε₁和ε₂均实现了第一阶收敛,且收敛性与参数无关。
  • 计算得到的收敛阶稳定在约0.483,表明在所有测试参数取值下均具有鲁棒的收敛性。
  • 计算得到的(ε₁,ε₂)-一致误差常数在N ≥ 128时稳定在约2.755,证实了参数一致性。
  • 数值结果表明,随着N增大,误差减小,并在ε₁和ε₂减小时趋于稳定,验证了理论收敛行为。
  • 该方法成功捕捉了预期的层结构:在u₁中表现为ε₁层,在两个分量中均表现为ε₂层,具体取决于边界层条件。
  • 对解的正则部分与层部分的分解进行了改进并正确定义,为理解解的结构提供了更清晰的洞察。

更好的研究,从现在开始

从阅读论文到最终审阅,大幅缩短您的研究时间。

无需绑定信用卡

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