Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
河村 拓馬; 井戸村 泰宏
Journal of Visualization, 23(4), p.695 - 706, 2020/08
被引用回数：0 パーセンタイル：100(Computer Science, Interdisciplinary Applications)粒子ベースボリュームレンダリングに基づくIn-situ可視化システムは、多変量ボリュームレンダリングに基づく高度にスケーラブルで柔軟な視覚解析環境を提供する。従来のCPUプラットフォームでは優れた計算パフォーマンスを示したが、最新のメニーコアプラットフォームでの高速計算により、関数パーサーと粒子データI/Oに関連する性能のボトルネックが明らかになった。本論文では、新しいSIMD対応の関数パーサーと、タスクベースのスレッド並列化に基づく非同期データI/Oを開発する。8208基のIntel Xeon Phi7250 (Knights Landing)プロセッサーで構成されるOakforest-PACS上の数値実験では、強スケーリングを約100kコアに維持しながら、一桁高い処理速度を実証した。
河合 智賀*; 井戸村 泰宏; 小川 雄一*; 山田 弘司*
Physics of Plasmas, 27(8), p.082302_1 - 082302_11, 2020/08
被引用回数：0弱磁気シアにおける大域的ジャイロ運動論モデルに基づいてトロイダル電子温度勾配駆動(ETG)乱流を調べた。大域的分布効果のために高トロイダルモード数nのトロイダルETGモードは外側の磁気面で励起され、強い線形分散をもたらす。この結果得られる非等方な波-乱流境界とエネルギー逆カスケードが帯状流の自己組織化を生成する。これは大域的ジャイロ運動論モデル特有の機構である。この自己組織化はランダムノイズによって初期化した減衰乱流とトロイダルETG乱流の両方で確認された。また、イオン電子温度比と乱流強度が決める臨界パラメータによってこの自己組織化過程が帯状流と等方的渦を生成することも示した。
長谷川 雄太; 小野寺 直幸; 井戸村 泰宏
第25回計算工学講演会論文集(CD-ROM), 4 Pages, 2020/06
都市部における風況と汚染物質拡散のアンサンブル計算を行うため、局所細分化格子ボルツマン法を用いたGPUペースのCFDコードを開発した。本コードはPascalまたはVolta世代のGPUアーキテクチャ向けに最適化されており、数km四方の計算領域・数m解像度の格子において実時間で風況解析を行うことができる。開発したコードを用いて、オクラホマシティで行われたフィールド実験JU(Joint Urban)2003の再現計算を行った。計算では、風況は実験値とよく一致し、また、トレーサ濃度のアンサンブルの平均値および最大値がFactor2(計算値と実験値の比が1/2倍
2倍)の条件を満たすことを確認した。
小野寺 直幸; 井戸村 泰宏; Ali, Y.*; 下川辺 隆史*; 青木 尊之*
第25回計算工学講演会論文集(CD-ROM), 4 Pages, 2020/06
原子力機構では3次元多相流体解析手法としてJUPITERを開発している。本研究では、JUPITERの圧力Poisson方程式解法として、適合細分化格子(AMR)を用いたマルチグリッド前提条件付き共役勾配法(P-CG)を開発した。計算の高速化として、全ての計算カーネルはCUDAを用いて実装すると共に、GPUスーパーコンピュータ上にて高い性能を発揮する様に最適化した。開発したマルチグリッド圧力Poisson解法は、オリジナルのP-CG法と比較して約1/7の反復回数で収束することが確認された。また、TSUBAME3.0上で8から216GPUまでの強スケーリング性能測定により、更なる3倍の高速化が達成された。
小野寺 直幸; 井戸村 泰宏; 上澤 伸一郎; 山下 晋; 吉田 啓之
Mechanical Engineering Journal (Internet), 7(3), p.19-00531_1 - 19-00531_10, 2020/06
福島第一原子力発電所の廃炉の方法の一つとして、空気のみで冷却する乾式法が挙げられる。日本原子力研究開発機構(JAEA)では、空冷性能のCFD評価手法としてJUPITERコードを開発している。しかしながら、JUPITERコードにおいて、複雑な原子炉内構造物を捉えた解析を実施するためには、非常の多くの計算資源が必要となる。このような問題に対して、本研究ではGPUスーパーコンピュータおよび適合格子細分化(AMR)法を適用した格子ボルツマン法に基づくCityLBMコードを開発している。CityLBMにて乾式法を模擬したJAEAの実験に対して検証計算を行なった結果、AMR格子の一様格子への収束性、および、実験値の再現を確認した。同じ解像度および同数の並列数にて計算速度の比較を行った結果、4台のGPU(NVIDIA Tesla V100)を用いたCityLBM法は、36台のCPU(Intel Xeon E5-2680v3)を用いたJUPITERの6.7倍の速度にて解析が可能であることが示された。
朝比 祐一*; Latu, G.*; Bigot, J.*; 前山 伸也*; Grandgirard, V.*; 井戸村 泰宏
Concurrency and Computation; Practice and Experience, 32(5), p.e5551_1 - e5551_21, 2020/03
被引用回数：0 パーセンタイル：100(Computer Science, Software Engineering)2つのジャイロ運動論コード、GYSELA, GKVを最新のアクセラレータ環境、Xeon Phi KNL, Tesla P100 GPUに移植した。一台のSkylakeプロセッサーに比べ、KNLにおけるGYSELAの逐次計算カーネルは1.3x、P100 GPUにおけるGKVの逐次計算カーネルは7.4x高速化された。GYSELAとGKVのスケーリングテストをそれぞれ16-512 KNLおよび32-256 P100 GPUで実施し、GYSELAのセミラグランジアンカーネルおよびGKVの畳み込みカーネルにおけるデータ転置通信が主要なボトルネックとなることがわかった。この通信コストを削減するために、これらのコードにパイプライン法およびタスク並列法に基づく通信オーバーラップを実装した。
井戸村 泰宏
Physics of Plasmas, 26(12), p.120703_1 - 120703_5, 2019/12
被引用回数：0 パーセンタイル：100(Physics, Fluids & Plasmas)本論文では大域的full-fジャイロ運動論シミュレーションを用いてイオン温度勾配駆動(ITG)乱流が支配的な水素(H)および重水素(D)のLモードプラズマにおける水素同位体と規格化ジャイロ半径
の影響を示す。断熱応答電子を用いたイオン加熱数値実験では、エネルギー閉じ込め時間はBohm的なスケーリングで決まり、水素同位体質量依存性は見られない。運動論的電子を用いた電子加熱数値実験は電子からイオンへの衝突性エネルギー移行に対する同位体効果に起因する明確な同位体質量依存性を示し、約1.4倍の加熱パワー比でH, Dプラズマは同様のイオン,電子温度分布を示す。H, Dプラズマの規格化無衝突イオンジャイロ運動論方程式はが同じ場合に等価となり、ITG乱流に対する衝突効果も弱い。このため、水素同位体質量依存性は主にスケーリングと加熱源によって決まる。
Ali, Y.*; 小野寺 直幸; 井戸村 泰宏; 伊奈 拓也*; 今村 俊幸*
Proceedings of 10th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA 2019), p.1 - 8, 2019/11
被引用回数：2 パーセンタイル：3.15大規模線形問題の反復法ソルバはCFDコードで共通に用いられる。前処理付共役勾配(P-CG)法は最も広く用いられている反復法の一つである。しかしながら、P-CG法では、特に演算加速環境において、大域的集団通信が重要なボトルネックとなる。この問題を解決するために、省通信版のP-CG法がますます重要になっている。本論文では多相CFDコードJUPITERにおけるP-CG法と前処理付チェビシェフ基底省通信CG(P-CBCG)法を最新のV100GPUに移植する。全てのGPUカーネルは高度に最適化され約90%のルーフライン性能を達成し、ブロックヤコビ前処理はGPUの高い演算性能を引き出すように再設計し、さらに残された袖通信のボトルネックは通信と計算のオーバーラップによって回避した。P-CG法とP-CBCG法の全体性能は大域的集団通信と袖通信の省通信特性によって左右され、GPUあたりのノード間通信帯域が重要となることが示された。開発したGPUソルバはKNLにおける以前のCPUソルバの2倍に加速され、Summitにおいて7,680GPUまで良好な強スケーリングを達成した。
井戸村 泰宏; 渡邉 智彦*; 藤堂 泰*
シミュレーション, 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
被引用回数：0 パーセンタイル：100(Physics, Fluids & Plasmas)Full-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
被引用回数：2 パーセンタイル：34.07(Computer Science, Interdisciplinary Applications)陰的時間積分スキームを用いてSugama衝突演算子をジャイロ運動論的ブラゾフコードGKVに実装した。新手法はオペレータスプリッティング、陰的時間積分、クリロフ部分空間反復法ソルバを用いており、線形化衝突演算子の詳細に依存しない汎用的な手法となっている。数値計算テストでは衝突項が制限する時間ステップ幅を超えて安定な計算が実証できた。データ転置を用いることで、反復計算中に通信が発生しない分散メモリシステムにおける効率的実装を実現した。これにより、本手法は計算効率向上と計算コスト削減を同時に達成し、アプリケーションの全体性能を大きく加速する。
小野寺 直幸; 井戸村 泰宏; Ali, Y.*; 下川辺 隆史*
Proceedings of 9th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA 2018) (Internet), p.9 - 16, 2018/11
被引用回数：2 パーセンタイル：12.32計算の高速化に向けて適合細分化格子(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
被引用回数：0 パーセンタイル：100多相流体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
被引用回数：6 パーセンタイル：41.93(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
被引用回数：6 パーセンタイル：4.68本研究では、局所細分化格子を適用した格子ボルツマン法を開発した。計算コードは、東京工業大学の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.472.39倍加速され、1,280ノード処理におけるデータ縮約通信は全体コストの513%から約1%に削減された。
