| 研究概要 |
第一段階として、減衰および外力の付加により定常状態にある非圧縮性一様等方性乱流のDNSデータの生成を行った.格子点数は512×512×512とし、計算には大阪大学サイバー・メディア・センター設置のSX-5を利用し、研究経費その他中の計算機使用料は、この使用料に充てられた.第二段階として,この両DNSデータを用いて、渦層から渦管への変換過程を解析し、渦管が2葉の渦層の相互作用により形成され、2葉の渦層上の渦度ベクトルの配置によって3種類の過程が存在する事を示した.特に、両ベクトルが直交する場合、渦層上の伸張方向のひずみ速度に圧縮が発生し、伴って形成される渦管は渦層上の渦度ベクトルにたいし横断方向に軸を持つ事、又、この配置による変換過程の発生頻度の増加につれて顕著なエネルギーのカスケードが起きる事を示した.第三段階として、Oldroyd構成方程式を用いて変換過程の生成にたいする粘弾性効果を解析した.その結果、変換過程の発生と乱流の生成の抑止効果は、Oldroyd-A構成方程式がB構成方程式に比べて大きい事を示した.これは、B方程式では粘弾性エネルギーが主に渦管に沿って分布し、粘弾性応力による圧力が中心部で極大値をとる事により渦管の生成を抑止するのにたいし、A方程式では主に渦層に沿って分布し、渦管への変換を抑止するためである事を明らかにした.次に、両Oldroyd構成方程式の定常解を求め、LESの非線形モデルとの類似点と相違点を示し、粘弾性効果を取り入れたSGSモデルをA方程式に基づいて提案した.海外共同研究者のHolm研究員は、Navier-Stokes-αmodelと上記定常解の関係を明らかにし、N.S.-αmodelが定常解とは異なる非対称なテンソル・モデルを与える事を示した.
|
