[論文レビュー] Mathematical modeling of pattern formation in sub- and supperdiffusive reaction-diffusion systems

V. V. Gafiychuk, Bohdan Datsko+1 more|ArXiv.org|Nov 2, 2006

Nonlinear Dynamics and Pattern Formation参考文献 25被引用数 7

ひとこと要約

本稿は、分数階微分（Caputo型、0 < α < 2）を用いた反応拡散系におけるパターン形成を調査する。分数階指数αは分岐パラメータとして機能する：α < α₀のとき、系は古典的RD系と類似した安定で定常状態の散逸的構造を示す。一方、α > α₀では、周期的またはカオス的空間時間パターンが出現し、均一な振動が高次のαでパターンを破壊する。臨界定数α₀は、振動不安定性の発生を示し、解析的に特定され、数値シミュレーションでも確認されている。

ABSTRACT

This paper is concerned with analysis of coupled fractional reaction-diffusion equations. It provides analytical comparison for the fractional and regular reaction-diffusion systems. As an example, the reaction-diffusion model with cubic nonlinearity and Brusselator model are considered. The detailed linear stability analysis of the system with Cubic nonlinearity is provided. It is shown that by combining the fractional derivatives index with the ratio of characteristic times, it is possible to find the marginal value of the index where the oscillatory instability arises. Computer simulation and analytical methods are used to analyze possible solutions for a linearized system. A computer simulation of the corresponding nonlinear fractional ordinary differential equations is presented. It is shown that the increase of the fractional derivative index leads to the periodic solutions which become stochastic at the index approaching the value of 2. It is established by computer simulation that there exists a set of stable spatio-temporal structures of the one-dimensional system under the Neumann and periodic boundary condition. The characteristic features of these solutions consist in the transformation of the steady state dissipative structures to homogeneous oscillations or space temporary structures at a certain value of fractional index.

研究の動機と目的

• 分数階時間微分が反応拡散系におけるパターン形成に与える影響を分析すること。
• 分数階指数α₀を特定すること。これは、サブ拡散およびスーパー拡散系における振動不安定性が発生する臨界値である。
• 古典的（整数階）系と比較して、分数階反応拡散系のダイナミクスを検討すること。
• αの増加に伴い、定常状態の散逸的構造から空間時間的振動への遷移を調査すること。
• ノイマン境界条件および周期的境界条件下での空間時間的構造の頑健性と多様性を調査すること。

提案手法

• Caputo時間微分（α ∈ (0,2)）を用いて、連立分数階反応拡散方程式を定式化する。
• 均一平衡状態のまわりで線形化し、得られた固有値問題を解析することで、分数階系の線形安定性分析を実施する。
• 特徴的時間の比と分数階指数に依存する臨界定数α₀を解析的に導出する。これは、振動不安定性が発生する閾値としての役割を果たす。
• 線形化系の数値シミュレーションを用いて、振動的およびカオス的解の出現を分析する。
• ノイマンおよび周期的境界条件下で、三次非線形性およびブルッセルレーター・モデルを含む非線形分数階ODEおよびPDEをシミュレートする。
• スペクトル法を用いて空間モードを分析し、異なるα値における摂動の時間的変化を追跡する。

実験結果

リサーチクエスチョン

• RQ1 分数階微分指数αは、反応拡散系の安定性およびパターン形成にどのように影響を与えるか？
• RQ2 振動不安定性が発生する臨界値α₀は何か？また、これは系のパラメータによってどのように決定されるか？
• RQ3 分数階系における空間時間的パターンは、古典的（α=1）反応拡散系のそれとどのように比較できるか？
• RQ4 αがα₀を超えて増加する際、どのようなダイナミクス的挙動—周期的、カオス的、または均一な振動—が出現するか？
• RQ5 なぜ高次のαで散逸的構造が均一な振動に崩壊するのか？これは系の非線形性にどのように依存するか？

主な発見

• 臨界定数α₀は分岐パラメータとして機能し、α < α₀では安定で減衰する振動（damped oscillations）、α > α₀では不安定で振動的またはカオス的ダイナミクスに分岐する。
• α < α₀のとき、系は古典的反応拡散系と同一の形の定常状態の散逸的構造を示し、パターンの形状に顕著なずれは認められない。
• α = α₀で、系は持続的振動に遷移し、三次非線形性モデルでは均一な振動が空間パターンを破壊する。
• αがα₀を超えて増加するにつれ、系は周期的振動から準カオス的空間時間的構造に進化する。これはブルッセルレーター・モデルで特に顕著に観察される。
• ブルッセルレーター・モデルでは、散逸的構造の振幅が大きく、滑らかでないプロファイルであるため、高次のαでも非均一な振動が持続し、均一な振動に収束しない。
• 選択的分数階シミュレーションにより、第一変数の分数階微分が振動の出現に主として寄与することが示された。