[论文解读] Asymptotic Behaviours of Solutions for Finite Difference Analogue of the Chipot-Weissler Equation

Houda Hani, Moez Khenissi|arXiv (Cornell University)|Feb 10, 2015

Differential Equations and Numerical Methods参考文献 4被引用 3

一句话总结

本文分析了Chipot-Weissler方程的有限差分数值逼近，该方程是一个具有反应项与梯度相关扩散项组合的非线性抛物型PDE。在参数和初始数据受限的条件下，建立了数值解向连续解的收敛性，但揭示了数值解与解析解在爆破集上的差异，突显了在数值上捕捉渐近爆破行为的局限性。

ABSTRACT

This paper deals with nonlinear parabolic equation for which a local solution in time exists and then blows up in a finite time. We consider the Chipot-Weissler equation. We study the numerical approximation, we show that the numerical solution converges to the continuous one under some restriction on the initial data and the parameters of the non linearity. Moreover, we study the numerical blow up sets and we show that although the convergence of the numerical solution is guaranteed, the numerical blow up sets are sometimes different from that of the PDE

研究动机与目标

• 使用有限差分格式研究Chipot-Weissler方程解的数值行为。
• 在初始数据和参数p与q的特定约束下，研究数值解向真实解的收敛性。
• 将数值爆破集与连续PDE的实际爆破集进行比较。
• 分析连续与离散设定下解的渐近行为。

提出的方法

• 对Chipot-Weissler方程应用半隐式有限差分格式：$ u_t = u_{xx} + u^p - |u_x|^q $，定义在 $ x \in (-1,1), t > 0 $ 上。
• 采用均匀网格对空间和时间进行离散化，空间导数通过中心差分逼近，时间导数通过后向欧拉法逼近。
• 在初始数据及参数 $ p > 1 $，$ 1 \leq q \leq \frac{2p}{p+1} $ 的限制下，证明了数值解向连续解的收敛性。
• 通过随时间与空间追踪解的最大值，对爆破集进行数值分析。
• 理论分析结合能量估计与比较原理，以控制解的有界性并调控爆破时间。
• 通过数值实验比较离散解与连续解的爆破集。

实验结果

研究问题

• RQ1 在初始数据和参数p与q的适当条件下，有限差分格式是否收敛于Chipot-Weissler方程的真实解？
• RQ2 数值爆破集与连续PDE的实际爆破集相比如何？
• RQ3 数值格式能否准确捕捉解的渐近行为，特别是有限时间爆破？
• RQ4 参数p与q在数值解的收敛性与爆破动力学中起什么作用？

主要发现

• 在初始数据及参数 $ p > 1 $，$ 1 \leq q \leq \frac{2p}{p+1} $ 的适当限制下，数值解收敛于连续解。
• 尽管存在收敛性，数值爆破集可能与PDE的实际爆破集不同，表明离散逼近未能保持爆破结构。
• 当梯度项 $ |u_x|^q $ 较强时，数值解的爆破时间可能不收敛于真实爆破时间。
• 理论分析证实局部解存在且在有限时间内爆破，与PDE行为一致。
• 数值模拟表明，即使解在范数意义下收敛，爆破位置在离散化下仍可能发生偏移。
• 爆破集的差异源于反应项与梯度项之间的非线性耦合，该耦合在有限差分格式中未能被充分保留。