[论文解读] Finite difference approximations for the first-order hyperbolic partial differential equation with point-wise delay

Paramjeet Singh, Kapil K. Sharma|arXiv (Cornell University)|Dec 5, 2010

Stability and Controllability of Differential Equations参考文献 23被引用 3

一句话总结

本文提出了一种二阶精度的有限差分格式——Lax-Friedrichs 和 Leap-Frog——用于求解具有点延迟的一阶双曲PDE，确保了相容性、在CFL条件下的稳定性以及收敛性。该方法通过构建一种特殊网格，使延迟自变量与网格点对齐，从而能够处理小延迟和大延迟，经由一维和二维数值算例验证，空间和时间方向均表现出二阶收敛性。

ABSTRACT

Explicit numerical methods based on Lax-Friedrichs and Leap-Frog finite difference approximations are constructed to find the numerical solution of the first-order hyperbolic partial differential equation with point-wise delay or advance, i.e., shift in space. The differential equation involving point-wise delay and advance models the distribution of the time intervals between successive neuronal firings. We construct higher order numerical approximations and discuss their consistency, stability and convergence. The numerical approximations constructed in this paper are consistent, stable under CFL condition, and convergent. We also extend our methods to the higher space dimensions. Some test examples are included to illustrate our approach. These examples verify the theoretical estimates and shows the effect of point-wise delay on the solution.

研究动机与目标

• 开发适用于具有点延迟或提前的一阶双曲PDE的更高阶数值方法，扩展先前的一阶格式。
• 确保在非恒定、光滑系数条件下，数值格式在CFL条件下的稳定性和收敛性。
• 构建一种专门的有限差分网格，使延迟的空间自变量与网格点对齐，从而实现对延迟值的精确逼近。
• 将该方法推广至二维及更高维空间，同时保持二阶精度。
• 通过数值实验案例展示点延迟对解行为的影响。

提出的方法

• 将Lax-Friedrichs和Leap-Frog有限差分格式应用于具有点延迟的双曲PDE，采用包含延迟空间域 $[-\alpha, 0]$ 的交错网格。
• 构建一种特殊网格，使得延迟自变量 $x - \alpha$ 恰好映射到一个网格点，从而确保离散解可直接访问延迟值。
• 通过使用二阶精度的有限差分近似替换时间与空间导数，推导出格式，保持相容性。
• 在CFL条件下分析稳定性，该条件限制时间步长相对于空间步长，以维持数值稳定性。
• 通过在 $x$ 和 $y$ 方向应用相同的有限差分模板，将方法扩展至二维空间，边界条件定义在 $[-\alpha, 0] \times [-\beta, 0]$ 上。
• 通过初始数据和边界数据迭代计算数值解，延迟值从存储的网格点中检索。

实验结果

研究问题

• RQ1 能否为具有点延迟的一阶双曲PDE构造出二阶精度的有限差分格式，以确保相容性和收敛性？
• RQ2 CFL条件如何影响Lax-Friedrichs和Leap-Frog等显式格式在存在延迟时的稳定性？
• RQ3 点延迟在多大程度上改变了解的形态？该方法能否准确捕捉这种影响？
• RQ4 该方法能否在保持二阶精度和稳定性的同时，推广至二维及更高维空间？
• RQ5 数值误差如何随网格细化而变化？该方法是否在空间和时间方向均实现二阶收敛？

主要发现

• 所提出的Lax-Friedrichs和Leap-Frog格式对于具有点延迟的双曲PDE具有相容性，且在CFL条件下无条件稳定，并具有收敛性。
• 数值实验结果证实了在空间和时间方向上的二阶收敛性，误差表显示当网格间距减半时，误差减少约四倍。
• 对于示例1（$\alpha = 0.05$），Lax-Friedrichs格式的最大绝对误差从 $\Delta x = 1/100$ 时的 0.053623 降低至 $\Delta x = 1/800$ 时的 0.004520，表明实现了二阶收敛。
• 对于示例2（$\alpha = 0.1$），Leap-Frog格式的均方根误差也表现出类似降低，从 0.041308 降至 0.000612，对应相同的网格细化序列。
• 在 $\alpha = \beta = 0.5$ 的二维数值实验中，验证了该方法在高维空间中的适用性，$t = 0.5$ 时刻的解稳定且精确，可视化结果清晰。
• 在图中可明显观察到延迟对解行为的影响，表现为相位偏移和振幅调制，与理论预期一致。