Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
井戸村 泰宏; 渡邉 智彦*; 藤堂 泰*
シミュレーション, 38(2), p.79 - 86, 2019/06
次世代核融合実験炉ITERにおける炉心プラズマ性能の評価・予測、および、より高性能な炉心運転シナリオの探索に向けて、ポスト「京」をターゲットとするエクサスケール核融合プラズマシミュレーションの研究開発を推進している。本稿では、これまでに開発してきたエクサスケール計算技術や既存実験装置を対象とする実証研究の成果を紹介し、ポスト「京」における核融合プラズマシミュレーションの展望を議論する。
朝比 祐一*; Grandgirard, V.*; Sarazin, Y.*; Donnel, P.*; Garbet, X.*; 井戸村 泰宏; Dif-Pradalier, G.*; Latu, G.*
Plasma Physics and Controlled Fusion, 61(6), p.065015_1 - 065015_15, 2019/05
パーセンタイル：100Full-fジャイロ運動論コードGYSELAを用いて輸送過程へのポロイダル対流セルの影響を調べた。この目的のために、対流セルのフィルタを適用し、フィルタ有無のシミュレーション結果を比較した。フィルタを適用することで磁気ドリフトに駆動されるエネルギー束が半減することがわかった。対流セルの周波数スペクトは乱流レイノルズ応力テンソルの周波数と対応し、対流セルが乱流によって駆動されることを示した。この対流セルの効果は乱流と新古典のダイナミクスの相互作用と考えられる。
小野寺 直幸; 井戸村 泰宏; 河村 拓馬; 上澤 伸一郎; 山下 晋; 吉田 啓之
Proceedings of 27th International Conference on Nuclear Engineering (ICONE-27) (Internet), 6 Pages, 2019/05
福島第一原子力発電所の廃炉の方法の一つとして、乾式法が挙げられる。日本原子力研究開発機構(JAEA)では、空冷性能のCFD評価手法としてJUPITERコードを開発している。しかしながら、JUPITERコードにおいて、複雑な原子炉内構造物を捉えた解析を実施するためには、非常の多くの計算資源と計算時間が必要となる。このような問題に対して、本研究ではGPUスーパーコンピュータに適した格子ボルツマン法に基づくCityLBMコードを開発している。CityLBMにてDry methodを模擬したJAEAの実験に対して検証計算を行なった結果、JUPITERコードと同様の結果が得られることが示された。また、同じ解像度および同数の並列数にて計算速度の比較を行った結果、GPUを用いたCityLBM法は、CPUを用いたJUPITERの1/6の計算時間にて解析が行えることが示された。以上の結果より、CityLBMは熱流動解析コードの有効な手法の一つであることが示された。
前山 伸也*; 渡邉 智彦*; 井戸村 泰宏; 仲田 資季*; 沼波 政倫*
Computer Physics Communications, 235, p.9 - 15, 2019/02
パーセンタイル：100(Computer Science, Interdisciplinary Applications)陰的時間積分スキームを用いてSugama衝突演算子をジャイロ運動論的ブラゾフコードGKVに実装した。新手法はオペレータスプリッティング、陰的時間積分、クリロフ部分空間反復法ソルバを用いており、線形化衝突演算子の詳細に依存しない汎用的な手法となっている。数値計算テストでは衝突項が制限する時間ステップ幅を超えて安定な計算が実証できた。データ転置を用いることで、反復計算中に通信が発生しない分散メモリシステムにおける効率的実装を実現した。これにより、本手法は計算効率向上と計算コスト削減を同時に達成し、アプリケーションの全体性能を大きく加速する。
小野寺 直幸; 井戸村 泰宏; Yussuf, A.*; 下川辺 隆史*
Proceedings of 9th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA 2018) (Internet), p.9 - 16, 2018/11
計算の高速化に向けて適合細分化格子(AMR)法を適用した格子ボルツマン法(LBM)に対して、通信削減マルチタイムステップ法(CRMT)を提案した。本手法はテンポラルブロッキング法に基づく定式化を行うことで、GPU計算で大きなボトルネックとなる通信回数の削減が可能となる。東京工業大学のTSUBAMEおよび東京大学のReedbushスーパーコンピュータにて性能測定を実施した結果、通信コストが64%に削減され、200GPUまでの弱および強スケーリング結果が改善された。以上の高速化により、2km四方の計算領域に対して1m解像度の風速5msの実時間解析が可能であることが示された。
井戸村 泰宏; 伊奈 拓也*; 山下 晋; 小野寺 直幸; 山田 進; 今村 俊幸*
Proceedings of 9th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA 2018) (Internet), p.17 - 24, 2018/11
多相流体CFDコードJUPITERの圧力ポアソン方程式に省通信マルチグリッド前処理付共役勾配(CAMGCG)法を適用し、省通信クリロフ部分空間法と計算性能と収束特性を比較した。JUPITERコードにおいてCAMGCGソルバ問題サイズによらずロバーストな収束特性を有し、通信削減と収束特性向上を両立することから、通信削減のみを実現する省通信クリロフ部分空間法に対する優位性が高い。CAMGCGソルバを
億自由度の大規模多相流体CFDシミュレーションに適用して反復回数を前処理付CG法の
に削減し、Oakforest-PACSにおける8,000ノードまでの良好な強スケーリングとCG法の
倍の性能向上を達成した。
小野寺 直幸; 井戸村 泰宏
Proceedings of 26th International Conference on Nuclear Engineering (ICONE-26) (Internet), 7 Pages, 2018/07
放射性物質の拡散計算は、核セキュリティの観点から非常に重要である。近年、低消費電力および演算性能の高さよりGPUによる科学技術計算が注目されている。本研究では、GPUを用いた格子ボルツマン法による計算コード手法を構築することで、高精度かつ高速な風況解析手法を構築した。東京工業大学のスーパーコンピュータTSUBAMEによる弱スケーリングの性能測定を行った結果、4から144台のGPUに対して良いスケーリングが得られるとともに、CPUとの比較ではノードあたり30倍程度の性能を達成した。産業技術総合研究所の風洞実験に対する検証計算では、AMR法を適用することで全てを細かい格子で計算するのと比較して10%程度の格子点数で解析が可能となるとともに、実験による観測と良く一致している結果が得られた。
松岡 清吉; 井戸村 泰宏; 佐竹 真介*
Physics of Plasmas, 25(2), p.022510_1 - 022510_10, 2018/02
被引用回数：1 パーセンタイル：53.13(Physics, Fluids & Plasmas)プラズマ輸送シミュレーションにおいて、ジャイロ運動論方程式をスケール分離を仮定せずプラズマ全体について第一原理に基づいて解く大域的full-fシミュレーション手法が注目されている。本研究では、複雑な3次元磁場形状をもつステラレータ型プラズマを大域的full-fジャイロ運動論シミュレーションコードであるGT5Dに適用するため、GT5Dの差分計算手法の拡張、及び3次元磁場平衡構築コードVMECとのインターフェースの開発を行った。開発したコード群の妥当性検証を目的として輸送に関する標準的なベンチマーク計算を行った。その結果、GT5Dが局所理論解析モデルや他の大域的新古典輸送コードの結果とよく一致することを示した。
井戸村 泰宏; 伊奈 拓也*; 真弓 明恵; 山田 進; 今村 俊幸*
Lecture Notes in Computer Science 10776, p.257 - 273, 2018/00
前処理付チェビシェフ基底省通信共役勾配(P-CBCG)法を多相熱流体CFDコードJUPITERにおける圧力ポアソン方程式に適用し、8,208台のKNLプロセッサを搭載したOakforest-PACS上で計算性能と収束特性を前処理付共役勾配(P-CG)法や前処理付省通信共役勾配(P-CACG)法と比較した。P-CBCG法は収束特性のロバースト性を維持しつつ集団通信回数を削減する。このロバースト性向上により、P-CACG法と比べて一桁以上大きい省通信ステップ数を実現する。2,000プロセッサを用いた場合、P-CBCG法はP-CG法, P-CACG法と比べてそれぞれ1.38倍, 1.17倍高速であることを示した。
小野寺 直幸; 井戸村 泰宏
Lecture Notes in Computer Science 10776, p.128 - 145, 2018/00
本研究では、局所細分化格子を適用した格子ボルツマン法を開発した。計算コードは、東京工業大学のGPUベースのスーパーコンピュータTSUBAME3.0を用いて開発を行い、最新のPascalアーキテクチャに対して最適化を行なった。1から36ノードを用いた弱スケーリングの性能測定では、NVIDIA TESLA P100を用いたGPU計算がBroadwellによるCPU計算の10倍以上の高速化が達成された。
山下 晋; 伊奈 拓也*; 井戸村 泰宏; 吉田 啓之
第31回数値流体力学シンポジウム講演論文集(DVD-ROM), 7 Pages, 2017/12
原子力機構では、過酷事故時の炉内溶融物移行挙動進展を現象論的に評価するために、数値流体力学的手法に基づく3次元多相多成分熱流動解析手法(JUPITER)の開発を行っている。原子炉過酷事故では、燃料棒1本スケール(数cm)から圧力容器スケール(十数メートル)まで非常に広範囲に渡るマルチスケール・マルチフィジクス現象を含むため、溶融物の移行挙動の詳細な予測のためには、高性能な大規模計算が必要不可欠となる。このような問題を解決するため、大規模並列計算に適した計算スキームの適用だけでなく、計算負荷の多くを占める圧力Poissonソルバーに対する反復解法を超並列計算に適したソルバーを適用した。これら手法を用いて、原子炉過酷事故時の炉内状況予測に向けた溶融物移行挙動解析手法の有効性評価と大規模計算への対応のため、実験解析及び予備解析を実施し、妥当性並びに大規模計算における良好な並列性能を有することを確認した。
井戸村 泰宏; 伊奈 拓也*; 真弓 明恵; 山田 進; 松本 和也*; 朝比 祐一*; 今村 俊幸*
Proceedings of 8th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA 2017), p.7_1 - 7_8, 2017/11
ジャイロ運動論的5次元オイラーコードGT5Dに省通信一般化最小残差(CA-GMRES)法を適用し、一般化共役残差(GCR)法を用いたオリジナルコードとの性能比較をJAEA ICEX(Haswell)、Plasma Simulator(FX100)、Oakforest-PACS(KNL)において実施した。CA-GMRES法はGCR法に比べて約3.8倍の演算密度となることから、メモリとネットワークの帯域が制限された将来のエクサスケールアーキテクチャに適合する。性能評価の結果、GCR法に比べて計算カーネルは1.47
2.39倍加速され、1,280ノード処理におけるデータ縮約通信は全体コストの513%から約1%に削減された。
山下 晋; 伊奈 拓也; 井戸村 泰宏; 吉田 啓之
Nuclear Engineering and Design, 322, p.301 - 312, 2017/10
被引用回数：4 パーセンタイル：18.09(Nuclear Science & Technology)過酷事故時の沸騰水型原子炉における溶融物の詳細な挙動について大きな注目を集めている。この挙動を明らかにするために、原子力機構では3次元多相多成分熱流動解析コードJUPITERを開発している。本論文では、JUPITERの妥当性を確認するために計算手法の基礎的妥当性検証及び実験結果との比較を実施し、良好な一致を得ることができた。加えて新たに開発したハイブリッド並列Poissonソルバーを導入することによって劇的に性能が向上した。そして、スーパーコンピュータ「京」において20万コアまでのストロングスケーリングを達成した。これらJUPITERの物理的、計算機的能力は、過酷事故時の各種溶融現象の評価を可能にするものと言える。
朝比 祐一*; Grandgirard, V.*; 井戸村 泰宏; Garbet, X.*; Latu, G.*; Sarazin, Y.*; Dif-Pradalier, G.*; Donnel, P.*; Ehrlacher, C.*
Physics of Plasmas, 24(10), p.102515_1 - 102515_17, 2017/10
パーセンタイル：100(Physics, Fluids & Plasmas)トカマクプラズマにおける熱流駆動型のイオン温度勾配乱流を計算するために2つの大域的full-Fジャイロ運動論コードのベンチマークを行う。この目的のために、full-Fジャイロ運動論方程式を現実的な熱流束固定条件で計算するセミ・ラグランジアンコードGYSELA、および、オイラーコードGT5Dを採用する。時空間特性に注目して雪崩的な輸送現象を評価した。自己組織化臨界現象(SOC)的な振舞いを議論するために統計解析を実施し、両方のコードで高周波側で
スペクトルから
スペクトルへの遷移を確認した。このベンチマークに基づき、SOC的な振舞いは数値計算法に依存しないロバーストな特徴であることを検証した。
松岡 清吉; 井戸村 泰宏; 佐竹 真介*
Physics of Plasmas, 24(10), p.102522_1 - 102522_9, 2017/10
被引用回数：1 パーセンタイル：79.99(Physics, Fluids & Plasmas)軸対称磁場閉じ込め装置であるトカマク型プラズマにおいて、誤差磁場等に起因する3次元的な非軸対称成分を持つ摂動磁場の効果が閉じ込め性能向上や不安定性制御の観点から注目されている。近年、3次元摂動磁場が駆動する衝突性粘性について、従来より広く用いられてきたバウンス平均理論モデルと大域的運動論シミュレーションの両者で得られた結果が一致しないことが指摘され、その物理機構の解明が課題となっていた。本研究では、二つの異なるモデルに基づいた大域的運動論シミュレーションを用いて詳細な解析を行い、上記の差異が生じる物理機構について検討した。その結果、理論モデルと大域的シミュレーションの不一致が、(1)粒子軌道の効果によって粒子軌道の共鳴構造が失われ、粒子が実効的に感じる摂動磁場強度が弱まること、(2)粒子軌道に沿った位相混合により速度空間構造内に微細構造が生成され輸送量の減衰が起きること、の二つが原因であることを明らかにした。
Mantica, P.*; Bourdelle, C.*; Camenen, Y.*; Dejarnac, R.*; Evans, T. E.*; G
rler, T.*; Hillecheim, J.*; 井戸村 泰宏; Jakubowski, M.*; Ricci, P.*; et al.
Nuclear Fusion, 57(8), p.087001_1 - 087001_19, 2017/08
被引用回数：1 パーセンタイル：79.99(Physics, Fluids & Plasmas)本会議報告は2016年9月5日8日にスイス、レザンにて開催されたthe 21st Joint EU-US Transport Task Force Workshopにおける報告と議論をまとめたものである。ワークショップは以下の8つのトピックから構成された: full-F運動論的乱流シミュレーションの進展; 高Z、低Z不純物の輸送、制御、および、プラズマ閉じ込めへの影響; 炉心、周辺輸送への3次元効果(MHD、外部磁場、ステラレータを含む); 予測性の高いプラズマ設計; 電子熱輸送とマルチスケール統合; Scrape-Off Layer(SOL)におけるパワー減衰長の理解; L-H遷移におけるSOLの役割; 乱流計測に対する乱流基礎特性の実証研究。
井戸村 泰宏
Physics of Plasmas, 24(8), p.080701_1 - 080701_5, 2017/08
被引用回数：3 パーセンタイル：47.07(Physics, Fluids & Plasmas)大域的full-fジャイロ運動論モデルに基づく電子加熱変調数値実験から、電子加熱によるイオン温度勾配駆動(ITG)乱流から捕捉電子モード(TEM)乱流への遷移が密度勾配の急峻化とプラズマ回転の変化をもたらすことを示した。回転変化時のトロイダル角運動量バランスをトロイダル角運動量保存則の直接観測によって明らかにし、イオン系の乱流応力に加え、イオン系の新古典応力、径方向電流、イオン系と電子系のトロイダル電場応力が重要となることを示した。ITGフェーズとTEMフェーズにおけるトロイダルトルクの反転はイオン系方向電流の反転によるものであり、これは粒子輸送と運動量輸送の相互作用であることを明らかにした。イオン系と電子系の径方向電流は両極性条件を満たすようにバランスし、電子系の径方向電流は電子系のトロイダル電場応力によって打ち消される。これは電子系のトロイダル電場応力が間接的にトロイダルトルクに影響していることを示す。
山田 進; 伊奈 拓也*; 佐々 成正; 井戸村 泰宏; 町田 昌彦; 今村 俊幸*
Proceedings of 2017 IEEE International Parallel & Distributed Processing Symposium Workshops (IPDPSW) (Internet), p.1418 - 1425, 2017/08
本発表では、BaileyのDouble-Doubleアルゴリズムを利用した4倍精度基本線形代数演算ライブラリ「BLAS」の高速化と、それを4倍精度固有値計算ルーチンに適用した際の計算性能の向上について発表する。特筆すべき成果は、積和演算の中間結果を高精度で保持して計算できるFMA命令を利用することで、我々がこれまでに開発した4倍精度BLASと比較し、代表的なルーチンで2045%の高速化を実現したことである。さらに、4倍精度固有値ソルバ「QPEigenK」が利用している4倍精度BLASをFMA命令を利用して開発したものに置き換えるだけで、原子力機構のICE Xおよび東京大学のFX10の2つのスパコンの大規模並列計算において1020%の高速化を実現した。この成果は、計算機の大規模化に伴って必要になりつつある高精度計算の高速な実行に資する成果である。
朝比 祐一*; Latu, G.*; 伊奈 拓也; 井戸村 泰宏; Grandgirard, V.*; Garbet, X.*
IEEE Transactions on Parallel and Distributed Systems, 28(7), p.1974 - 1988, 2017/07
被引用回数：1 パーセンタイル：60.51(Computer Science, Theory & Methods)セミ・ラグランジュ法における間接メモリアクセス、有限差分法におけるストライドメモリアクセスといった複雑なメモリアクセスパターンを有する核融合プラズマ乱流コードの高次元ステンシル計算をGPGPUやXeon Phiプロセッサ等の演算加速器上で最適化した。どちらのデバイスでも、Array of Structure of Array (AOSOA)データレイアウトが連続的なメモリアクセスに有効である。Xeon Phiでは時空間データ局所性の向上によるローカルキャッシュの効率的利用が必要不可欠である。GPGPUではテクスチャメモリの利用がセミ・ラグランジュ法の間接メモリアクセス性能を向上する。これらの最適化により、アクセラレータ用核融合カーネルはCPU用カーネルに比べてSandy Bridge (CPU)用最適化コードに比べて1.4x - 8.1x高速化した。
河村 拓馬; 野田 智之; 井戸村 泰宏
Supercomputing Frontiers and Innovations, 4(3), p.43 - 54, 2017/07
最新のメニーコア・プラットフォーム上で、In-Situ粒子ベースボリュームレンダリング(In-Situ PBVR)に基づくIn-Situデータ探索フレームワークのパフォーマンスを検証した。In-Situ PBVRはコストの高い順序計算を回避することができ、並列モンテカルロサンプリングによって、大規模なボリュームデータを小さな可視化用粒子データに変換する。この特徴により、大規模シミュレーションのボトルネックである、1ノードあたりのメモリ制限や、計算速度とノード間通信速度のパフォーマンスギャップを回避することができる。さらに、サイズの小さな粒子データをクライアントPCに転送することによる自由な視点位置のボリュームレンダリング、そして非同期のファイルベースの制御シーケンスによるバッチ処理実行時の可視化パラメータを変更という、In-Situデータ探索が可能になった。In-Situ PBVRは、8208個のIntel Xeon Phi7250(Knights Landing)プロセッサで構成されたOakforest-PACS上で、約10万コアまでの強スケーリング(問題サイズを固定し、CPU数を増加させて処理を加速させるスケーリング)と省メモリ特性を示した。
