[論文レビュー] 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

ひとこと要約

本稿は、反応項と勾配に依存する拡散項を併せ持つ非線形放物型PDEであるChipot-Weissler方程式の有限差分数値近似を分析している。パラメータおよび初期データの制約の下で、数値解が連続解に収束することを確立しているが、数値的・解析的 Blow-up 集合の間には乖離が生じており、漸近的 Blow-up 挙動を数値的に捉えることの限界を示している。

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 に特定の制約を課した下で、数値解が真の解に収束するかを調査すること。
• 離散的解のBlow-up集合と連続PDEの実際のBlow-up集合を比較すること。
• 連続的および離散的設定における解の漸近的挙動を分析すること。

提案手法

• Chipot-Weissler方程式 $ u_t = u_{xx} + u^p - |u_x|^q $ に半implicitな有限差分スキームを適用する。$ x \in (-1,1), t > 0 $ で定義される。
• 空間と時間に一様グリッドを用い、空間微分は中央差分で、時間微分は後退Euler法で近似する。
• 初期データおよびパrameter $ p > 1 $, $ 1 \leq q \leq \frac{2p}{p+1} $ に制限を課した下で、数値解が連続解に収束することを証明する。
• Blow-up集合は、時間と空間全域における解の最大値を追跡することで数値的に分析する。
• 理論的分析ではエネルギー推定と比較原理を組み合わせ、解の境界を定め、Blow-up時間を制御する。
• 離散的および連続的解のBlow-up集合を比較するための数値実験を実施する。

実験結果

リサーチクエスチョン

• RQ1 初期データおよびパラメータ p と q に適切な条件を課した場合、有限差分スキームはChipot-Weissler方程式の真の解に収束するか？
• RQ2 数値的Blow-up集合は、連続PDEの実際のBlow-up集合とどのように異なるか？
• RQ3 特に有限時間 Blow-up を含む解の漸近的挙動は、数値スキームで正確に捉えることができるか？
• RQ4 パラメータ p と q は、数値解の収束性およびBlow-upダイナミクスにどのような役割を果たすか？

主な発見

• 初期データおよびパラメータ $ p > 1 $, $ 1 \leq q \leq \frac{2p}{p+1} $ に適切な制限を課した下で、数値解は連続解に収束する。
• 収束はしているが、数値的Blow-up集合とPDEの実際のBlow-up集合には乖離が生じており、離散近似においてBlow-up構造が保存されていないことを示唆している。
• 特に勾配項 $ |u_x|^q $ が強い場合、数値解のBlow-up時間は真のBlow-up時間に収束しないことがある。
• 理論的分析により、局所解が存在し、有限時間でBlow-upすることが確認され、PDEの挙動と整合的である。
• 数値シミュレーションでは、解のノルム収束が見られる中で、Blow-up位置が離散化によってずれることが示された。
• Blow-up集合の乖離は、反応項と勾配項の非線形的結合が有限差分スキームで完全に保存されていないことに起因している。