双安定反応拡散方程式の 有限領域内の界面の持続時間

Presentation on theme: "双安定反応拡散方程式の 有限領域内の界面の持続時間"— Presentation transcript:

1 双安定反応拡散方程式の 有限領域内の界面の持続時間
堀川　洋 香川大学工学部情報工学科

2 1. 問題 1次元双安定反応拡散方程式の1次元有限領域における不安定界面の持続時間 (安定定常状態までの過渡時間)
1次元双安定反応拡散方程式の1次元有限領域における不安定界面の持続時間　(安定定常状態までの過渡時間) 界面の移動速度は領域の長さに対して指数関数的に小さくなる。 → 界面の持続時間は領域の長さに対して指数関数的に増加する。 本研究では以下の2点について調べた。 ･ランダムな初期値から生じる界面の持続時間の分布 → 界面の持続時間は逆べき型の分布に従う。 ･時間空間的雑音の界面の持続時間への影響 → 一定の強度の雑音によって平均持続時間が長くなる。 2 2

3 2. モデル 1次元有限領域内での双安定反応拡散方程式の初期境界値問題 安定一様解：u = -1, 1 (不安定一様解：u = 0)
ランダムな初期条件：E{φ(x)φ(x’)} = δ(x – x’) → 過渡的な界面パターン　・・・ 極めて長い持続時間 → 安定一様解(u(x) = -1 or 1) Demonstration 3 3

4 界面の位置(u(x0, t) = 0)：x0 = l/2 - l+ l+：界面の左側の領域の長さ
界面の方程式[1, 2] 界面の位置(u(x0, t) = 0)：x0 = l/2 - l+ l+：界面の左側の領域の長さ 領域長：l → 大 界面の速度 → 極めて小 過渡時間(空間一様解に収束するまでの時間)：T → 領域長：lに対して指数関数的に増加[2] T ∝ exp(l+), exp(l - l+) l l+ -l/ l/2 x0 4 4

5 3. 界面の持続時間 3.1. 矩形パターンからの持続時間 初期値：l+(0) = l0 (< l/2)
例えば、y = exp(αl+)とおく → 式(2)の解： l+(T) = 0 → 界面の持続時間：T(l0) 領域の長さ：無限大(l → ∞) → 簡単な式： 持続時間T･･･初期長さl0に対して指数関数的に増加 5 5

6 3.2. ランダム初期値からの持続時間の分布 初期の界面の長さ ～ 領域内(0, l)に一様分布 → 界面の持続時間：Tの分布：h(T)
カットオフ時間：Tc = exp(αl/2)/(αk) ･T < Tc　逆べき型の分布 ･T > Tc　指数分布 カットオフよりも長い持続時間の割合：Prob{T > Tc} ≈ 21/2exp(-2)/l ≈ 0.19/l　 → lの逆べきで減少 6 6

7 図2 ランダムな初期値から発生する界面の持続時間の分布：h(T) (領域長：l = 15， Tc ≈ 3.4×107)
持続時間の平均、分散 → 指数関数的に増加 h(T) ∝ 1/T 図2　ランダムな初期値から発生する界面の持続時間の分布：h(T) 　 (領域長：l = 15， Tc ≈ 3.4×107) 図3　領域の長さ：lに対するランダムな初期値からの界面の平均持続時間：m(T) 7 7

8 3.3. 雑音の影響 式(1)の右辺に空間･時間的な強度σxの白色雑音が加わった場合 Demonstration
界面の持続時間・・・式(2)の右辺に白色雑音を加えた次式の初通過時間問題 持続時間：T l+(0) = l0 → l+(T) = 0 or l T0 T1 T3 T4 8 8

9 界面の片側の領域の長さl+：l0 → 0 への初通過時間(FPT)の平均と分散[3]：
l+ = 0：吸収境界、l+ = l/2：反射境界 界面の片側の領域の長さl+：l0 → 0 への初通過時間(FPT)の平均と分散[3]： 平均持続時間は非単調に変化し、ある強度の雑音によって最大となる。 ただし、初期長さl0につれてピークの位置も移動するため、ランダム初期値から生じる界面の分布は単調減少となる。 図4　界面の平均持続時間の雑音強度(標準偏差)に対する変化　 (l = 15, l0 = 4) 9 9

10 3.4. 空間2次元の場合 Demonstration
ランダムな初期値からの過渡時間は、1次元での表式(式(6) – (9))で良く近似される。 図A2　2次元領域(-l/2 < x, y < l/2)におけるランダムな初期値からの平均過渡時間：m(T) 図A1　2次元領域(-5 < x, y < 5 (l = 10))におけるランダムな初期値からの過渡時間の分布：h(T) 10 10

11 4. 反応項の非対称性の影響 l(0) lc 1 - δ 分岐点：lc l+ → 0 (l(0) < lc)
l+ → l (l(0) > lc) l(0) lc -l/ l/2 11 11

12 矩形パターンからの持続時間 初期長さ：l0に対する持続時間：Tの増加は指数関数的なものから線形なものになる。
図A4　平衡点のずれ：δに対する界面の持続時間：T (l = 15, l0 = 4) 図A3　初期長さ：l0に対する界面の持続時間：T 12 12

13 5. まとめ 双安定反応拡散方程式の有限領域における界面の持続時間について調べた。
　双安定反応拡散方程式の有限領域における界面の持続時間について調べた。 界面の移動速度は領域の長さに対して指数関数的に小さくなる。 → 界面の持続時間は領域の長さに対して指数関数的に増加する[1, 2]。 　本研究では、 ･ランダムな初期値から生じる界面の持続時間の分布 → 持続時間の分布にはカットオフが存在 　　カットオフ以下では持続時間は逆べき型の分布に従う 　　カットオフ以上では指数分布に従うが、その割合は界面の長さの逆数に比例して 　　減少する ･時間空間的雑音の界面の持続時間への影響 → 一定の強度の雑音によって平均持続時間が長くなる。 　　Noise-sustained patterns 13 13

14 関連事項 ･領域長が大きい場合 → 複数(多数)の界面からなるパルス状パターン 隣接する界面は衝突と合体を繰り返して空間一様解に収束
　　　　　　　　 　　　　　　　　隣接する界面は衝突と合体を繰り返して空間一様解に収束 Demonstration 界面の速度 ～ 界面間の距離の指数関数で近似 各衝突までの時間 ～ 距離に対して指数関数的に増加[1, 2] その裏返しとして、界面間の平均距離は時間の対数に比例してゆっくりと増加[4] ･環状発振器型の神経回路網モデルの過渡振動の持続時間 → 回路内の素子数に対して指数関数的に増加[5]。 いずれの場合も、不安定解(界面と周期解) → 双安定な空間一様解のセパラトリクスを成す 　　ただ1次元の不安定固有空間を持つ 不安定空間の固有値は系の大きさ(領域の長さと素子数)に対して指数関数的に減少[2, 6] 14 14