[論文レビュー] Can a traveling wave connect two unstable states? The case of the nonlocal Fisher equation

Grégoire Nadin, Benoı̂t Perthame+1 more|arXiv (Cornell University)|Nov 20, 2010

Mathematical and Theoretical Epidemiology and Ecology Models参考文献 6被引用数 6

ひとこと要約

本稿は、非局所的フィsher方程式における進行波解の性質を調査している。ここでは、定常状態 $ u=0 $ は単安定性のため、$ u=1 $ はチューリング不安定性のため、両者とも動的に不安定である。著者らは、有限領域における人工的境界条件を用いた新しい数値アルゴリズムを開発し、$ u=1 $ がチューリング不安定であっても、単調な進行波が $ 0 $ から $ 1 $ へと接続可能であることを示した。これは、非局所性に起因する不安定性の下で進行波が存在しないはずという古典的予想に反する。

ABSTRACT

This note investigates the properties of the traveling waves solutions of the nonlocal Fisher equation. The existence of such solutions has been proved recently in \cite{BNPR} but their asymptotic behavior was still unclear. We use here a new numerical approximation of these traveling waves which shows that some traveling waves connect the two homogeneous steady states $0$ and $1$, which is a striking fact since $0$ is dynamically unstable and $1$ is unstable in the sense of Turing.

研究の動機と目的

• 非局所反応拡散系において、2つの不安定な定常状態を接続する進行波が存在するかという未解決問題を解明すること。
• 定常状態 $ u=1 $ がチューリング不安定である場合の非局所的フィッシャー方程式における進行波解の漸近的挙動を明確にすること。
• 非局所相互作用および不安定な境界条件の下でも、高い精度で進行波解を計算可能な数値手法の開発と検証を行うこと。
• $ u=1 $ がチューリング不安定である場合、波が $ u=1 $ に接続するか、安定な周期的状態に接続するかを、特に $ \mu $ が大きい場合に特定すること。
• 非局所効果によって進行波の単調性が失われる臨界閾値 $ \mu_c \approx 8.9 $ を特定すること。

提案手法

• 有限領域 $[x_l, x_r]$ における進行波方程式を解くための数値アルゴリズムを考案。人工的境界条件として $ v(x_l) = 1 $, $ v(x_r) = 0 $, および $ v(0) = \epsilon $ を設定。
• 計算領域を $ I_1 = [x_l, 0] $ と $ I_2 = [0, x_r] $ に分割し、それぞれにディリクレ境界条件を課した2つの楕円型問題を別々に解く。
• 畳み込み項 $ \phi \star v $ は、$ v $ を $ (-\infty, x_l) $ で 1 に、$ (x_r, \infty) $ で 0 に拡張することで、非局所的評価を可能にする。
• 解の連続性を確保するため、$ x=0 $ における微分のジャンプを最小化するように、波の位相速度 $ \sigma^L $ をシューティング法で調整する。
• 収束性は、空間グリッド幅 $ \Delta x $ の細分化と領域長 $ L $ の増加によって評価され、誤差は $ x=0 $ における $ |E(0)| $ で測定される。
• 本手法は、一様核 (5)、非負で compact な台を持つ核 (6a)、(6b)、ガウス核 (7) の複数のカーネルに対してテストされ、手法の頑健性とカーネル依存性が評価された。

実験結果

リサーチクエスチョン

• RQ1 非局所的フィッシャー方程式において、動的に不安定な状態 $ u=0 $ とチューリング不安定な状態 $ u=1 $ を接続する進行波解が存在するか？
• RQ2 $ u=1 $ にチューリング不安定性が存在する場合、単調な進行波の存在が妨げられるのか、それとも依然として存在可能なのか？
• RQ3 進行波の単調性が失われる臨界的な $ \mu $ の値は何か？
• RQ4 特にフーリエ変換が負となる核（例：$ \sin \xi / \xi $）を用いた場合、波の形状や安定性にどのような影響を与えるか？
• RQ5 初期データが compact な台を持つ初期値問題の解は、進行波に吸引されるのか、それとも脈動的前線のようなより複雑な構造と共存するのか？

主な発見

• $ \mu=1000 $ および $ \mu=2500 $ の場合、核 (5) を用いた数値的検証により、$ u=1 $ がチューリング不安定であっても、$ u=0 $ から $ u=1 $ へと接続する進行波解が存在することが示された。
• $ \mu=10 $ の場合、波はほぼ単調であり、最大値が $ 1.0028 $ に達するなど、弱い非単調性を示す。$ \mu=10 $ 時の最大値はわずかに 1 を超える。
• 単調性の喪失の臨界閾値は、$ \mu_c \approx 8.9 $ として数値的に確認され、$ x \to -\infty $ 近傍での線形化から得られる解析的推定値と整合的である。
• $ \mu $ が増加するにつれ、波の尾部に減衰する振動が現れ、その振幅は領域長 $ L $ の増加に伴いゆっくり減少する。これは、$ 1 $ から $ 0 $ へと接続する波に収束している可能性を示唆する。
• 一様核 (5)、負のフーリエ変換を持つ核（例：$ \sin \xi / \xi $）およびガウス核 (7) を含む、すべてのテストしたカーネルに対して、進行波解は $ u=1 $ に接続されており、周期的状態には接続されていない。
• 数値的手法により、$ \Delta x $ を小さくし、$ L $ を大きくするにつれ、誤差 $ |E(0)| $ が減少し、$ \mu=64 $ の場合に $ \sigma^L \to 15.367 $ に収束する。これは理論的波速 $ 2\sqrt{64} = 16 $ に近い。