水中シルトフェンスの変形・移動を可能とする3次元モデリング; 流体・シルトフェンス連成計算の実現

山田 進; 町田 昌彦; 田中 みのり*; 関 克己*; 有川 太郎*

日本応用数理学会論文誌, 31(1), p.20 - 43, 2021/03

https://doi.org/10.11540/jsiamt.31.1_20

原子力規制庁の受託研究に基づき開発した水路等に設置したたわみのあるシルトフェンスを流れ場に合わせて変形させるアルゴリズムを報告する。さらに、このアルゴリズムを用いて、水路内の流れ場とシルトフェンスの連成シミュレーションを実施しその変形を評価するとともに、実際の水路実験の結果との比較を行い、提案したアルゴリズムの妥当性を報告する。

Status of study of long-term assessment of transport of radioactive contaminants in the environment of Fukushima (FY2018) (Translated document)

長尾 郁弥; 新里 忠史; 佐々木 祥人; 伊藤 聡美; 渡辺 貴善; 土肥 輝美; 中西 貴宏; 佐久間 一幸; 萩原 大樹; 舟木 泰智; et al.

JAEA-Research 2020-007, 249 Pages, 2020/10

JAEA-Research-2020-007.pdf:15.83MB

2011年3月11日に発生した太平洋三陸沖を震源とするマグニチュード9.0の東北地方太平洋沖地震とそれに伴って発生した津波により、東京電力(現東京電力ホールディングス)福島第一原子力発電所の事故が発生し、その結果、環境中へ大量の放射性物質が放出された。この事故により放出された放射性核種は、その大部分が森林に沈着している。これに対し、面積が広大であり大量の除去土壌などが生じる、多面的な森林の機能が損なわれる可能性があるなどの問題があり、生活圏近傍を除き、汚染された森林の具体的な除染計画はない。そのため、未除染の森林から放射性セシウムが流出し、既に除染された生活圏に流入することで空間線量率が上がってしまうのではないか(外部被ばくに関する懸念)、森林から河川に流出した放射性セシウムが農林水産物に取り込まれることで被ばくするのではないか、規制基準値を超えて出荷できないのではないか(内部被ばくに関する懸念)などの懸念があり、避難住民の帰還や産業再開の妨げとなる可能性があった。日本原子力研究開発機構では、環境中に放出された放射性物質、特に放射性セシウムの移動挙動に関する「長期環境動態研究」を2012年11月より実施している。この目的は、自治体の施策立案を科学的側面から補助する、住民の環境安全に関する不安を低減し、帰還や産業再開を促進するといった点にある。本報告書は、原子力機構が福島県で実施した環境動態研究におけるこれまでの研究成果について取りまとめたものである。

大規模多相流体解析向け省通信型マルチグリッド前処理付き共役勾配法

井戸村 泰宏; 小野寺 直幸; 山田 進; 山下 晋; 伊奈 拓也*; 今村 俊幸*

スーパーコンピューティングニュース, 22(5), p.18 - 29, 2020/09

多相多成分熱流動解析コードJUPITERの圧力ポアソン方程式に省通信型マルチグリッド前処理付き共役勾配(CAMGCG)法を適用し、従来のクリロフ部分空間法と計算性能と収束特性を比較した。CAMGCGソルバは問題サイズによらずロバーストな収束特性を示し、通信削減と収束特性向上を両立することから、通信削減のみを実現する省通信クリロフ部分空間法に対する優位性が高い。CAMGCGソルバを900億自由度の大規模多相流体解析に適用し、前処理付共役勾配法ソルバと処理性能を比較した。このベンチマークにおいて、反復回数は約1/800に削減され、Oakforest-PACS上で8,000ノードに至る良好な強スケーリングを維持しつつ約11.6倍の性能向上を達成した。

High performance eigenvalue solver for Hubbard model; Tuning strategies for LOBPCG method on CUDA GPU

山田 進; 町田 昌彦; 今村 俊幸*

Parallel Computing; Technology Trends, p.105 - 113, 2020/00

https://doi.org/10.3233/APC200030

本発表は科学研究費補助金(科研費)に従い実施した強相関ハバードモデル計算に現れる固有値問題に対する高性能計算に関するものである。具体的には、ハバードモデルの計算に現れる固有値計算に固有値計算ソルバの1つであるLOBPCG法を適用した際の高速化についての発表である。特筆すべき成果は、現在主流のプロセッサの1つであるGPUのアーキテクチャに合わせたデータの格納方法を提案し、実際に行列計算を高速化したことである。さらに、複数の線形計算をまとめて実行することで、データへのアクセス回数を減らすことができ、さらなる高速化も実現した。これらの高速化により、これまでの方法と比較し全体で約1.4倍の高速化を実現した。なお、この成果は科研費研究「エクサスケール計算機を想定した量子モデルシミュレーションに対する並列化・高速化」の研究成果である一方、GPUを利用した高性能計算にも資する成果である。

疎行列の固有値計算

山田 進

固有値計算と特異値計算; 計算力学レクチャーコース, p.77 - 104, 2019/12

本書「固有値計算と特異値計算」は一般の研究者向けに固有値計算及び特異値計算の基礎理論から応用例までを解説する書籍である。発表者は本書の第3章「疎行列の固有値計算」を執筆し、疎行列向きの固有値計算方法として頻繁に用いられているLanczos法やLOBPCG法などのアルゴリズムを解説している。また、実際にこれらの解法を用いて固有値問題を計算した結果を示すとともに、それぞれの解法の有効性や問題点などを紹介している。

福島における放射性セシウムの環境動態研究の現状(平成30年度版)

長尾 郁弥; 新里 忠史; 佐々木 祥人; 伊藤 聡美; 渡辺 貴善; 土肥 輝美; 中西 貴宏; 佐久間 一幸; 萩原 大樹; 舟木 泰智; et al.

JAEA-Research 2019-002, 235 Pages, 2019/08

JAEA-Research-2019-002.pdf:21.04MB

2011年3月11日に発生した太平洋三陸沖を震源とするマグニチュード9.0の東北地方太平洋沖地震とそれに伴って発生した津波により、東京電力福島第一原子力発電所の事故が発生し、その結果、環境中へ大量の放射性物質が放出され、その大部分が森林に沈着している。これに対し、面積が広大であり大量の除去土壌等が生じる、多面的な森林の機能が損なわれる可能性があるなどの問題があり、生活圏近傍を除き、汚染された森林の具体的な除染計画はない。そのため、未除染の森林から放射性セシウムが流出し、既に除染された生活圏に流入することに対する懸念があり、避難住民の帰還や産業再開の妨げとなる可能性があった。原子力機構では、環境中に放出された放射性物質、特に放射性セシウムの移動挙動に関する「長期環境動態研究」を2012年11月より実施している。この目的は、自治体の施策立案を科学的側面から補助する、住民の環境安全に関する不安を低減し、帰還や産業再開を促進するといった点にある。本報告書は、原子力機構が福島県で実施した環境動態研究におけるこれまでの研究成果について取りまとめたものである。

Communication avoiding multigrid preconditioned conjugate gradient method for extreme scale multiphase CFD simulations

井戸村 泰宏; 伊奈 拓也*; 山下 晋; 小野寺 直幸; 山田 進; 今村 俊幸*

Proceedings of 9th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA 2018) (Internet), p.17 - 24, 2018/11

https://doi.org/10.1109/ScalA.2018.00006

多相流体CFDコードJUPITERの圧力ポアソン方程式に省通信マルチグリッド前処理付共役勾配(CAMGCG)法を適用し、省通信クリロフ部分空間法と計算性能と収束特性を比較した。JUPITERコードにおいてCAMGCGソルバ問題サイズによらずロバーストな収束特性を有し、通信削減と収束特性向上を両立することから、通信削減のみを実現する省通信クリロフ部分空間法に対する優位性が高い。CAMGCGソルバを億自由度の大規模多相流体CFDシミュレーションに適用して反復回数を前処理付CG法のに削減し、Oakforest-PACSにおける8,000ノードまでの良好な強スケーリングとCG法の倍の性能向上を達成した。

High performance LOBPCG method for solving multiple eigenvalues of Hubbard model; Efficiency of communication avoiding Neumann expansion preconditioner

山田 進; 今村 俊幸*; 町田 昌彦

Lecture Notes in Computer Science 10776, p.243 - 256, 2018/00

https://doi.org/10.1007/978-3-319-69953-0_14

本発表は科学研究費補助金(科研費)研究の一環として実施した量子問題計算に現れる大規模な対称疎行列であるハミルトニアン行列の複数の固有値および固有ベクトルを計算する際に用いる省通信ノイマン展開前処理の有効性についての報告である。具体的な成果は、ハミルトニアン行列の複数の固有値とそれに対応する固有ベクトルをLOBPCG法を用いて反復計算する際の収束性を向上させるための新規の前処理方法を提案し、実際の計算から既存の方法よりも短時間で計算できることを示したことである。さらに、問題の物理的性質を利用して前処理計算の通信回数を減少させるアルゴリズムを提案し、原子力機構のスパコンSGI ICEXおよび理化学研究所の京スパコンの数千コアを用いた並列計算において、高速化が実現することを確認した。この成果は、量子計算の高速計算だけでなく、今後大規模化していく並列計算機の有効利用に資する成果でもある。

Communication avoiding Neumann expansion preconditioner for LOBPCG method; Convergence property of exact diagonalization method for Hubbard model

山田 進; 今村 俊幸*; 町田 昌彦

Parallel Computing is Everywhere, p.27 - 36, 2018/00

https://doi.org/10.3233/978-1-61499-843-3-27

本発表は科研費研究の一環として実施した固有値計算ソルバLOBPCG法の前処理に関する研究についての発表である。LOBPCG法は反復解法であり適切な前処理を用いることで収束性が向上することが知られているが、量子問題に表れるハバードモデルの固有値を計算する際に物理パラメータによっては既存の前処理では収束性が向上しないことがある。そこで、線形方程式の反復解法で前処理として使われているノイマン展開を利用した前処理を適用し、そのような問題に対しても収束性が向上することを見出した。さらに、問題の物理的性質を考慮し、演算回数は若干増加するが、通信回数を減少させる方法を提案した。この方法で開発したコードを、原子力機構の並列計算機SGI ICEXを用いて並列シミュレーションしたところ、これまでの方法よりも約20%高速に計算できることを確認した。この成果は、固有値問題に有効な前処理方法を見出しただけではなく、ネットワーク構造が複雑化してきている並列計算機での性能向上にも資する成果である。

Quadruple-precision BLAS using Bailey's arithmetic with FMA instruction; Its performance and applications

山田 進; 伊奈 拓也*; 佐々 成正; 井戸村 泰宏; 町田 昌彦; 今村 俊幸*

Proceedings of 2017 IEEE International Parallel & Distributed Processing Symposium Workshops (IPDPSW) (Internet), p.1418 - 1425, 2017/08

https://doi.org/10.1109/IPDPSW.2017.42

本発表では、BaileyのDouble-Doubleアルゴリズムを利用した4倍精度基本線形代数演算ライブラリ「BLAS」の高速化と、それを4倍精度固有値計算ルーチンに適用した際の計算性能の向上について発表する。特筆すべき成果は、積和演算の中間結果を高精度で保持して計算できるFMA命令を利用することで、我々がこれまでに開発した4倍精度BLASと比較し、代表的なルーチンで2045%の高速化を実現したことである。さらに、4倍精度固有値ソルバ「QPEigenK」が利用している4倍精度BLASをFMA命令を利用して開発したものに置き換えるだけで、原子力機構のICE Xおよび東京大学のFX10の2つのスパコンの大規模並列計算において1020%の高速化を実現した。この成果は、計算機の大規模化に伴って必要になりつつある高精度計算の高速な実行に資する成果である。