2次元乱流は3次元乱流と異なる性状の代表例であるが、本研究では特に特異性の発生に焦点を絞って数値的に研究した。ここでは非粘性準地衡近似方程式に従う速度場を選び、大規模数値シミュレーションによって場の特異性の発生を調べた。この方程式は速度場と温度場の関係がLaplacianの平方根で結ばれている点で、2次元乱流と3次元乱流の中間に位置しており、94年にConstantinらによる数値実験(1024^2モード)から特異性が発生することが示唆されているものである。また彼らは解析的に、発散が生じるとすればそれは温度場の微分によるものであることを示している。 そこで本研究では、2048^2モードを用いて数値計算を実行した。その結果数値精度の範囲内において、非粘性問題における温度場の勾配について、その最大値は時間的に二重指数関数的に増大すること、またそのL^2ノルムは時間的に指数的に増大することを見い出した。また粘性問題においては、温度場のL^2ノルムの減衰率εについてε〜R^<-0.45>が成り立つこと(R:Reynolds数)、および温度勾配のL^2ノルムが最大になる時刻t_cについてt_c〜loglogRが成り立つことを見い出した。これらの結果はいずれも、温度勾配について、有限時間の発散を支持せず、従って、Contantinらによる有限時間の発散を支持する数値結果は精度的問題を含むことが結論される。
|