Matrix-free defect correctionを用いた陰的FR法の安定性解析

摘要

近年,複雑形状周りの実用LESを実現する次世代の流体ソルバ一の研究が盛hに行われており，非構造格子を用いながら高い空間精度を達きることが出来る，高次精度非構造格子法が注目を集めている.代表的な手法は不連続ガレルキン(Discontinuous Galerkn, DG)法であるが，最近では構築がより容易な流束再構築(Flux Reconstruction, FR)法が提案され，関連する研究が多数報告されている.これらの手法はセル内に高次精度化のための自由度を持ち，セル境界で近似リーマン解法を用い数値流束を評価することが特徴である.高次精度化のために必要な情報がセル内に集まっているため，领域分割計算のコア間の通信に必要なデータがセル界面のみのため，高い並列性を期待することが出来る.FR法をベースとしたオープンソースコードであるPyFRは18,000個のCPUを使用して，13.7PFLOP/s(Peak比58.0%)という非常に高い演算性能を自している.