これらのモデルに現れる典型的パターンとして、1次元では進行パルス、2次元ではリングパターンが挙げられる.これらのパターンの特徴は衝突消滅性である.この現象はBZ反応の実験にも観察されていたことから、2変数単安定興奮性反応拡散方程式系の進行パルスは衝突消滅することがこれまで信じられてきた.
しかしながら、近年、方程式系(1),(2)とは異なった形の興奮性反応拡散方程式系において、進行パルスが必ずしも衝突消滅しないことが報告された.Petrov(1994)らは、3点平衡点を持つ自己触媒モデルに対して、また、Krisher,Mikhailov(1994)、Kawaguchi,Mimura(1997)はglobal coupling項を持つ実質3変数の方程式系(1),(2)に対して、1次元進行パルスの反射、2次元スポットパターンの反射が起こることを報告している.
以上のことから、必ずしも衝突消滅は興奮性反応拡散方程式系の進行パルスの特徴的な現象でないことがわかってきた.しかしながら、FHNモデルは衝突消滅性を持つことから、いかなるタイプの興奮性反応拡散方程式が衝突消滅するのか、あるいはしないのかは依然として問題となっていた.
本申請論文は以上の問題を考察することである.すなわち系(1)に対して
(A-1)f(u,)=0は3次関数的非線形項とし,g(u,)=0は単調関数とする.そして、それらは交点を1点だけ持つ、(1)の拡散を無視した系では、小さいに対して、興奮性を持つ.
(A-2)(1)は安定進行パルス解を持つ.
を満たす2成分反応拡散方程式系に対して次の問題を提起する.
仮定(A-1),(A-2)のもとで、果たして方程式系(1)は衝突非消滅現象を起こすか?
この問題に対して、次の非線形項
を持つ1次発熱反応拡散方程式系(1)を考える.ここで、k(u)=exp(u/(1+u/c))である.ただし、(1)の変数u(t,x),v(t,x)はそれぞれ、時刻t場所xにおける温度と化学反応物質濃度に対応する.a,h,c,c1は正定数とする.
パラメータa,h,vc,c,を適当に固定すると、(3)は(A-1)を満たす.更に、dを小さくとると、安定進行パルス解が存在することがわかり、(A-2)を満たす.このとき、2つの1次元進行パルス及び2次元リングパターンは衝突消滅する.一方、安定進行パルス解が存在する((A-2)が成立する)範囲でdを大きくとると、進行パルスの速度は遅くなり、2つの進行パルスが反発反射する(図2).また、2次元ではリングパターンが崩壊し複雑なパターンが出現する.このことは、たとえ(A-1),(A-2)を満たしたとしても、衝突非消滅現象を起こし、その結果複雑なパターンが出現する興奮性反応拡散方程式が存在することを示唆している.ここではこのようなパターンダイナミクスの理解の第一歩として次の問題を考える.
(P-1)衝突非消滅現象の本質は何か?
本論文の目的は問題(P-1)に対して数値計算そして有限次元に縮約された系の解析から接近することである.論文の内容は(I)計算機支援による衝突非消滅現象の考察、(II)衝突非消滅現象の解析的考察、(III)数値実験による衝突非消滅性に起因する2次元パターンの考察、から成っている.
(I)計算機支援による衝突非消滅現象の考察 ここでは数値実験から問題(P-1)を明らかにすることが目的である.最初に、(1),(3)に対して、衝突非消滅性を持つ進行パルス解の特徴を明らかにするために、安定進行パルス解や安定振動進行パルス解の存在範囲を調べた.その結果、安定な速度の遅い進行パルスの存在が数値計算より示唆された.何故速度の遅い進行パルス解が存在するのかを明らかにするために、定常パルス解の解構造と定常パルス解からの分岐を数値計算によって調べた.数値計算によって、定常パルス解からのsuper-critical分岐による速度の遅い進行パルス解の出現し、衝突非消滅現象が出現することが数値計算によって示唆された.この分岐によって現れた進行パルス解は、dが大きいことから、の側方抑制性を持っていると考えられる.
側方抑制性とは次のことである:(1),(3)において、の拡散がuのそれより大きいとき、uが拡散と反応によって拡がるより、uの高い場所でが消費され、その周囲のが消費された部分に速く拡散してしまい、周囲のは少なくなってしまう.その結果uの拡がりが抑制され、進行パルスは速度が遅くなる.このような現象を引き起こす性質を側方抑制性という.
そこで、この側方抑制性を持つ進行パルスが存在すれば反射するのかを考えた.(3)とは異なるタイプの非線形項を持つGray-Scottモデルにも、衝突非消滅現象が出現することが知られているので、速度の遅い進行パルスが存在するか否かを調べることにした.その結果、定常パルス解からのsuper-critical分岐よって、速度の遅い進行パルス解は出現すること、その進行パルスは反射することが数値計算によって示唆された.(図3)この進行パルスもまた、側方抑制性をもっていることがわかる.
これらのことから、次を予想する.
(C-1)反応拡散系のタイプによらず、もしも側方抑制性を持つ速度の遅い進行パルス解が存在するならば、衝突非消滅現象の性質を持つのではないだろうか?
以下では、側方抑制性を持つ速度の遅いパルスを単に速度の遅い進行パルスと呼ぶ.
(II)衝突非消滅現象の解析的考察 予想(C-1)を解析的に論じるために速度の遅い進行パルス解が存在する系として、次の双安定FHN方程式を考える.
ここでf(u)=u(1-u)(u-a).ただし、,,d,a,は正定数,十分小さい時定数、0=0とu-=0は3点平衡点P,Q,Rを持ち、PとRは安定平衡点、Qは不安定平衡点となっている.
系(4)において、非常に速度の遅い進行パルス解が存在し、その解は反射することが数値計算によって示された(図4).そこで、進行パルス解の反射を議論するために、パルスの位置を↓0の界面Front、Back(1(t)>2(t))として定義し、(4)の極限方程式を導出する.この極限方程式に対する界面の相互作用問題において、数値計算によって界面の反射現象が確認できた(図5).
次に、界面だけで記述できる方程式を導出するために、f(u)に対して区分線形f(u)=-u+H(u-a)を仮定する.ただし、H(x)=1,x>0, H(x)=0,x<0である.ここで、界面の速度が非常に遅いと仮定して、界面の時間変化を記述する4次元常微分方程式系を導出することができた.この系は、数値計算によって、図5と定性的に類似した界面の反射現象を示すことが確認された(図6).また、速度の遅い界面は幅が非常に広いことがわかる.そこで幅が非常に広い場合、すなわち1-2≫0の場合には、界面の反射は本質的に、Front同士の反射、Back同士の反射、Front同士の反射の3つの過程から成ることが明らかになった.
こうして、この常微分方程式は速度の遅い進行パルス解の反射機構を理解するためのよい近似になっていることが示唆された.しかし、この正当化と界面反射の厳密な証明は今後に残された課題である.
(III)衝突非消滅性に起因する2次元パターンの考察 最後に、衝突非消滅現象に起因して出現する2次元時間空間パターンについて考察する.ここでは、(1),(3)に現れる複雑な2次元パターンの生成機構の本質は何であるかを数値計算によって考察した.最初に、プラナー進行パルスの安定性を調べた.その結果、速度の遅い進行パルスが存在するとき、プラナー進行パルスは不安定であることが示唆された.この結果が一般に成り立つかどうかを検証するために、Gray-Scottモデルに対して、プラナー進行パルスの性質を調べた結果、予想したように速度の遅い進行パルスはプラナー不安定であることが示唆された.この事実から、一般的に速度の遅い進行パルスは、衝突非消滅性だけでなくプラナー不安定性も持っていることが示唆された.そして、(1),(3)とGray-Scottモデルでは、このような状況下でのリングパターンは衝突非消滅性とプラナー不安定性によって崩壊し、複雑なパターンを生成することが示唆された.以上の結果は、反応拡散系のタイプによらず、もし速度の遅い進行パルスが存在するならば、2次元では複雑なパターンを生成するのではないか?という期待を抱かせる.
図1 (2)のf,gのヌルクライン.図2 (1),(3)の進行パルスの反射現象.図3 Gray-Scottモデルの進行パルスの反射現象.図4 (4)の進行パルスの反射現象.図5 極限方程式の界面の反射現象.図6 4次元常微分方程式での界面の反射現象.