| 研究概要 |
本年度は前年度までに求まった近似方程式の高精度検証を行った.まず蔵本・シバシンスキー方程式(K-S方程式)の数値計算法について計算スキームの高精度化をおこない,乱流状態のスペクトルを計算した結果,定波数から高波数までスペクトル法と完全に一致する高精度解が得られた.この計算コードを用いて前年度に行った有効粘性係数をカルマン・ブーシーフィルタを用いて計算した結果,その値はほぼ一致し,ノイズレベルの検討を新たに行った結果,この値がZ可測推定であると同時に,この有効粘性においてほぼノイズレベルが最小になっていることが明らかとなった.これはこの係数が物理的にも有効であることが確認されたことを意味する.
さらに大偏差統計を用いて揺動力項のクラメル関数を求めた.前年度までに得られた分布は粗視化パラメータTの値の妥当性の検証がなされていなかったが,本年度はこれをもパラメータにして適切なTを検討した結果,このTが前述の有効粘性係数を求めるのに必要なTのほぼ最小値となることが明らかとなった.以上のように本年度は得られた近似方程式が伊藤過程(ノイズを付加したバーガース方程式)で近似でき,、本研究で得られた係数が最適であることが明らかとなった.
一方,本年度の研究の過程でK-S方程式のエネルギー散逸構造を維持するためにはこのようにして得られた近似方程式を単純に数値積分(数値伊藤積分)するだけでは不十分であること,さらにそもそもこのような数値積分法は発展途上であり,和書は存在しないことが明らかとなった.これは新たな(数値)確率積分が必要であることを意味しており,文献サーベイも含め来年度の課題である.
|
