Next: 4 -周期解
Up: 単独保存則方程式の周期解の数値解析
Previous: 2 差分近似解
(PDF ��������: kiyou1997.pdf)
3 周期解への収束
周期解の存在証明[3] では, 外力の
-方向の全変動が
影響を与えていたので, 実験に用いた時間周期的外力は, 最大値を保
ったままで全変動を変えられるように次のような関数とした.
の全変動は
,
の最大値は
と
なる.
ただし, 以下の実験では, この外力の変動は変化させず,
すべて
,
と固定し,外力については
,
のみを変化
させて数値計算を行った.
また,
の初期値
は定数
で与えた.
-方向の
分割幅
は分割数
に対して
となるように定め,
-方向の分割幅
は, CFL 条件と,
がちょうど
の
整数倍になるように決定したが, 粘性効果が大きくなりすぎないように
なるべく大きな値に取った.
この外力の周期
と分割数
を変えて数値実験してみた結果を以
下に示す.
図 1:
:
,
,
=100,101,102,...,110
(
,
)
|
図 2:
:
,
,
=100,101,102,...,110
(
,
)
|
Fig. 1,2 は
のときの
の
グラフを
から
まで, すなわち外力の
周期目から
周期目までのグラフを重ねて書いたもので,
Fig. 1 は
, Fig. 2 は
のときの様子
である.
Fig. 1 のグラフはほとんど一本のグラフに見え, 関数
がほぼ同じものであることがわかる.
これはこの解がこの時点でほぼ
-周期解になっていることを意味する.
一方 Fig. 2 のグラフは一本のグラフにまとまっているとは
いいがたい.
の経過にともなうグラフのまとまり具合を, 外力の周期
を変化させて
ながめてみると, この二つの例のように非常に速くまとまる
とまとまりの
悪い
があり, しかもこれらが
の変化に対して入れ替わりにあらわれる.
このまとまりの悪い
では何らかの不安定さが働いているのではないかと
予想される.
Next: 4 -周期解
Up: 単独保存則方程式の周期解の数値解析
Previous: 2 差分近似解
Shigeharu TAKENO
2001年 7月 20日