• 検索結果がありません。

核融合プラズマ乱流シミュレーションの階層的並列化

N/A
N/A
Protected

Academic year: 2021

シェア "核融合プラズマ乱流シミュレーションの階層的並列化"

Copied!
1
0
0

読み込み中.... (全文を見る)

全文

(1)2011年ハイパフォーマンスコンピューティングと計算科学シンポジウム High Performance Computing Symposium 2011. HPCS2011 2011/1/18. 核融合プラズマ乱流シミュレーションの階層的並列化 井戸村 泰宏、 ジョリエ セバスチャン 日本原子力研究開発機構、東京都台東区東上野 6-9-3 核融合プラズマにおける乱流現象の第一原理モデルであるジャイロ運動論モデルに基づく 5 次元格子コード GT5D にお ける階層的並列化手法について議論する。GT5D では方程式系の物理的特性や対称性に基づいて 5 次元位相空間の階層 的な多次元領域分割を行う。これを複数の MPI コミュニケータと SMP を組み合わせた階層的通信構造によって実装す ることにより高いスケーラビリティを実現し、BX900 システム 16384 コアを用いて実効性能 19.4TFlops を達成した。. 散化し、①に関してはハミルトン方程式で与えられる非圧 縮乱流場による移流の保存則を満たす無散逸保存型差分 によって精度と数値安定性を両立し、②については中心差 分を適用する。時間積分には半陰的ルンゲ・クッタ法を使 用し、1010 次元の連立一次方程式として与えられる陰的ス テップを一般化共役残差法に基づく反復法によって解く。 ③のポアソン方程式はトーラス方向にフーリエモード展 開を適用し、トーラス断面を 2 次元有限要素近似で取り扱 うが、乱流場の非等方構造を反映したフィルターを適用し て有限要素演算子サイズを圧縮することにより、直接法で. 乱流場を解く。これらの差分・有限要素演算子に演算子の 物理的特性や対称性に応じて設計した階層的領域分割を 適用し、3 つの MPI コミュニケータと SMP によって構成 される階層的ハイブリッドネットワーク構造上で並列化 する。これによって 1 つの集団通信を行う MPI プロセス 数を百程度以下に抑制し、1 万コアを超える領域まで高い スケーラビリティを実現した。ここで、SMP によるハイ ブリッド並列化は MPI プロセス数削減による通信コスト 削減のみならず、MPI バッファ等によるメモリ使用量の削 減にも大きく寄与する。図 1 では複数の超並列プラットフ ォーム上で GT5D のストロングスケーリングを計測し、 BX900 シ ス テ ム 16384 コ ア の 実 効 性 能 と し て ~19.4TFlops(対ピーク性能比~10.1%)を達成した。講演 では演算カーネル性能やスケーラビリティの各プラット フォーム間の比較等についても議論する。 Sustained Performance (GFlops/total). 乱流輸送現象は核融合炉心プラズマの閉じ込め性能を 支配するため核融合エネルギー開発における重要な課題 となっている。中性流体乱流と比較して、核融合プラズマ 乱流は、①電子流体とイオン流体という移動度や特性長が 異なる流体成分が電磁場を介して相互作用する、②閉じ込 め磁場構造によって非等方な乱流構造が形成される、そし て、③約 1 億度という高温の核融合プラズマでは衝突周波 数が非常に小さく、波と粒子の共鳴現象や粒子軌道効果と いった運動論的効果が顕著となる、といった特徴がある。 3 次元空間の中性流体乱流シミュレーションと異なり、核 融合プラズマ乱流シミュレーションは、特に、③の運動論 的効果のために 5 次元位相空間(3 次元空間+2 次元速度) における粒子分布の取り扱いを必要とし、膨大な計算量に なる。このため、現実的なパラメータでシミュレーション を行うことが長年困難であったが、近年急速に発展してい る超並列計算機の利用によって実装置規模の 5 次元ジャイ ロ運動論シミュレーションが盛んに行われるようになっ てきた[1,2,3]。しかしながら、現在の計算資源では ITER (国際熱核融合実験炉)のような大型装置の計算は難しく、 将来のペタ超級の計算資源の利用、および、それに向けた シミュレーションの超並列最適化が必要とされている。 ジャイロ運動論方程式系は①乱流場による 5 次元位相空 間の移流、②粒子衝突効果による速度空間の移流・拡散、 ③乱流場を決定するポアソン方程式から構成される。数値 計算法としては粒子モデル、および、格子モデルという 2 つの異なるアプローチが存在するが、本研究で開発してい る GT5D[4,5]は格子モデルを採用している。トーラス状の 炉心プラズマの 5 次元位相空間を 1010 程度の構造格子で離. BX900 SR16k. 10000 . FX1 Altix T2K 1000 . BG/P. 100  100. 1000. 10000 Cores. 図 1:GT5D コードのストロングスケーリング 参考文献 [1] 渡邉, 洲鎌, 井戸村, プラズマ・核融合学会誌 81 (2005), pp.534-546, pp. 581-592, pp. 686-702. [2] Y. Idomura, T.-H. Watanabe, and H. Sugama, C. R. Physique 7, 650 (2006) [3] X. Garbet, Y. Idomura, L. Villard, T.-H. Watanabe, Nucl. Fusion 50, 043002 (2010). [4] Y. Idomura et al., Comput. Phys. Commun. 179, 391 (2008).; Nucl. Fusion 49, 065029 (2009). [5] S. Jolliet et al., submitted to J. Comput. Phys.. 58. ⓒ 2011 Information Processing Society of Japan.

(2)

参照

関連したドキュメント

自ら将来の課題を探究し,その課題に対して 幅広い視野から柔軟かつ総合的に判断を下す 能力 (課題探究能力)

心臓核医学に心機能に関する標準はすべての機能検査の基礎となる重要な観

しかしながら生細胞内ではDNAがたえず慢然と合成

ここで融合とは,バンカーが伝統的なエリートである土地貴族のライフスタ

線遷移をおこすだけでなく、中性子を一つ放出する場合がある。この中性子が遅発中性子で ある。励起状態の Kr-87

議論を深めるための参 考値を踏まえて、参考 値を実現するための各 電源の課題が克服さ れた場合のシナリオ

 福島第一廃炉推進カンパニーのもと,汚 染水対策における最重要課題である高濃度

入学願書✔票に記載のある金融機関の本・支店から振り込む場合は手数料は不要です。その他の金融機