研究課題
気液二相流計算の気液界面を捕獲スキームとしてVOF(Volume of Fluid)値を保存させるスキームとLevel Set 法を組み合わせた手法を試行した。界面の曲率が格子間隔に近づくと非保存形は体積を失い（小さな気泡などが消える）、保存形では位相誤差が大きくなるため存在しない気液界面（分裂）が発生することを確認した。この問題を解決するため、気液界面捕獲のためだけに流体を解く格子より細かいSub-gridを準備し、VOF値をそのSub-grid上で保存形を用いて計算する方法を検討した。Navior-Stokes方程式の非線形項や気液界面の曲率を精度よく求めるために、VOF値からLevel Set値を再構築により求める。そこでの反復計算には、高次精度のWENOスキームを導入した。これらは計算量が多く、GPU計算で効率良く計算するためには、GT200コアではGPUのShared Memoryの使い方やレジスタを有効に使う必要があったが、FermiコアのGPUではL1キャッシュに頼ることで十分高速計算できることが明らかになった。気液二相流の数値計算では、気体から液体に連続的に遷移する2～3メッシュの層を設け、ここでの1000倍の密度変化がPoisson方程式の微分演算子の内側に入り込み、疎行列計算を困難にする。BiCGStab法だけでなく前処理法が必要であり、不完全LU分解をスムーザとするMG(Multi-Grid)法を導入した。さらにGPU用のMG前処理には、隣接Thread間の依存性を排除するためにRed & Black法を導入した。
2: おおむね順調に進展している
年度当初に予定していたGPUのリリースが遅れるなどの理由で予算を繰り越したが、平成24年度でほぼ予定した研究計画を達成した。
平成25年度は当初の研究計画の通り進めることができる。
すべて 2013 2012 2011 その他
すべて 雑誌論文 (10件) (うち査読あり 7件) 学会発表 (32件) (うち招待講演 13件) 図書 (3件)
油空圧技術
巻: Vol. 51, No.3 ページ: 7-13
Journal of Computational Physics
巻: Volume 230, Issue 14 ページ: 5787-5805
日本機械学会論文集B編
巻: Vol.77, No.781 ページ: 1699-1714
日本流体力学会雑誌「ながれ」
巻: Vol.30 ページ: 75-78
Journal of Crystal Growth
巻: Vol.318 ページ: 40-45
Computer Methods in Biomechanics and Biomedical Engineering
巻: Volume 15, Issue 7, ページ: 771-778
Parallel Computing
巻: Volume 37, issue 9 ページ: 521-535
日本気象学会誌「天気」
巻: 新用語解説, Vol.58, No.7 ページ: 657-660
日本シミュレーション学会誌「シミュレーション」
巻: Vol.30, No.3, ISSN 0285-9947 ページ: 163-172
Progress in Nuclear Science and Technology
巻: Vol.2 ページ: 639-642
