Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
田中 みのり*; 渡部 真史*; 町田 昌彦; 山田 進; 榎本 容太*; 郡司 滉太*; 有川 太郎*
土木学会論文集,B2(海岸工学)(インターネット), 76(2), p.I_103 - I_108, 2020/11
原子力規制庁の受託研究に基づき実施した水路に2枚のシルトフェンスを設置した際の汚濁物質の拡散防止効果の評価を目標とした水理模型実験の結果報告である。シルトフェンスを適切な間隔で2枚設置すると1枚目のシルトフェンスの影響で2枚目のふかれ量が少なくなり、1枚だけの場合よりもより多くの汚濁物質が捕捉できるようになることを報告する。また、平面水槽実験においても同様の結果が得られたことも報告する。
長尾 郁弥; 新里 忠史; 佐々木 祥人; 伊藤 聡美; 渡辺 貴善; 土肥 輝美; 中西 貴宏; 佐久間 一幸; 萩原 大樹; 舟木 泰智; et al.
JAEA-Research 2020-007, 249 Pages, 2020/10
JAEA-Research-2020-007.pdf:15.83MB
2011年3月11日に発生した太平洋三陸沖を震源とするマグニチュード9.0の東北地方太平洋沖地震とそれに伴って発生した津波により、東京電力(現東京電力ホールディングス)福島第一原子力発電所の事故が発生し、その結果、環境中へ大量の放射性物質が放出された。この事故により放出された放射性核種は、その大部分が森林に沈着している。これに対し、面積が広大であり大量の除去土壌などが生じる、多面的な森林の機能が損なわれる可能性があるなどの問題があり、生活圏近傍を除き、汚染された森林の具体的な除染計画はない。そのため、未除染の森林から放射性セシウムが流出し、既に除染された生活圏に流入することで空間線量率が上がってしまうのではないか(外部被ばくに関する懸念)、森林から河川に流出した放射性セシウムが農林水産物に取り込まれることで被ばくするのではないか、規制基準値を超えて出荷できないのではないか(内部被ばくに関する懸念)などの懸念があり、避難住民の帰還や産業再開の妨げとなる可能性があった。日本原子力研究開発機構では、環境中に放出された放射性物質、特に放射性セシウムの移動挙動に関する「長期環境動態研究」を2012年11月より実施している。この目的は、自治体の施策立案を科学的側面から補助する、住民の環境安全に関する不安を低減し、帰還や産業再開を促進するといった点にある。本報告書は、原子力機構が福島県で実施した環境動態研究におけるこれまでの研究成果について取りまとめたものである。
山田 進; 町田 昌彦; 操上 広志
Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.140 - 146, 2020/10
本発表は、2011年の福島第一原子力発電所の事故により環境中に放出され、地上に降下した放射性セシウムのダム内での振る舞いの研究についての発表である。大きなダムでは、夏期に上層の水温が高く下層が低い状態が保たれる水温の成層化と呼ばれる状態が観測されており、成層化しているダム内での詳細な振る舞いを計算機シミュレーションにより評価するためには、その状態を考慮した機能を追加することが必須である。そこで、これまでに著者らが開発してきた3次元流体シミュレーションコード「3D-Sea-SPEC」に水平方向の移動に関する補正機能を追加することで、温度の成層化を再現することを可能にした。本発表では、「3D-Sea-SPEC」の紹介を行うとともに、成層化を再現するための補正方法について報告する。また、水温を変化させたシミュレーションを実施し、ダム内が成層化する条件を調査するとともに、水温の成層化の有無による放射性セシウムの振る舞いの違いを評価し、その結果を報告する。
井戸村 泰宏; 小野寺 直幸; 山田 進; 山下 晋; 伊奈 拓也*; 今村 俊幸*
スーパーコンピューティングニュース, 22(5), p.18 - 29, 2020/09
多相多成分熱流動解析コードJUPITERの圧力ポアソン方程式に省通信型マルチグリッド前処理付き共役勾配(CAMGCG)法を適用し、従来のクリロフ部分空間法と計算性能と収束特性を比較した。CAMGCGソルバは問題サイズによらずロバーストな収束特性を示し、通信削減と収束特性向上を両立することから、通信削減のみを実現する省通信クリロフ部分空間法に対する優位性が高い。CAMGCGソルバを900億自由度の大規模多相流体解析に適用し、前処理付共役勾配法ソルバと処理性能を比較した。このベンチマークにおいて、反復回数は約1/800に削減され、Oakforest-PACS上で8,000ノードに至る良好な強スケーリングを維持しつつ約11.6倍の性能向上を達成した。
町田 昌彦; 山田 進; 岩田 亜矢子; 乙坂 重嘉; 小林 卓也; 渡辺 将久; 船坂 英之; 森田 貴己*
Journal of Nuclear Science and Technology, 57(8), p.939 - 950, 2020/08
被引用回数：0 パーセンタイル：100(Nuclear Science & Technology)2011年4
5月にかけて発生した東京電力ホールディングス・福島第一原子力発電所2号機及び3号機からの汚染水の海洋への直接流出以後、神田は相対的に小さいが連続的な放射性物質の流出が引き続き起こっていることを指摘している。しかし、その期間は2012年9月までであり、その後の流出量の推定についての報告はない。そこで、本論文では、その後を含めて2011年4月から2018年6月までの7年間に渡り
Csの流出量を推定した結果を報告する。報告のない時期、国・東京電力ホールディングスは、流出を抑制するための努力を続け、港湾内海水の放射性核種濃度は徐々に減少している。われわれは、一月単位でCsの流出量を二つの手法、一つは神田の提案した手法だがわれわれの改良を加えた手法とボロノイ分割によるインベントリー評価法を使い評価した。それらの結果から、前者の手法は常に後者の手法と比べて保守的だが、前者の後者に対する比は1桁の範囲内であることが分かった。また、それらの推定量から簡単に沿岸域に対するインパクトを評価し、特に魚食による内部被ばく量を推定したところ、福島第一原子力発電所(1F)の海洋流出量に基づく内部被ばく分は極めて小さいことが分かった。
山田 進; 町田 昌彦; 今村 俊幸*
Parallel Computing; Technology Trends, p.105 - 113, 2020/00
本発表は科学研究費補助金(科研費)に従い実施した強相関ハバードモデル計算に現れる固有値問題に対する高性能計算に関するものである。具体的には、ハバードモデルの計算に現れる固有値計算に固有値計算ソルバの1つであるLOBPCG法を適用した際の高速化についての発表である。特筆すべき成果は、現在主流のプロセッサの1つであるGPUのアーキテクチャに合わせたデータの格納方法を提案し、実際に行列計算を高速化したことである。さらに、複数の線形計算をまとめて実行することで、データへのアクセス回数を減らすことができ、さらなる高速化も実現した。これらの高速化により、これまでの方法と比較し全体で約1.4倍の高速化を実現した。なお、この成果は科研費研究「エクサスケール計算機を想定した量子モデルシミュレーションに対する並列化・高速化」の研究成果である一方、GPUを利用した高性能計算にも資する成果である。
松本 和也*; 井戸村 泰宏; 伊奈 拓也*; 真弓 明恵; 山田 進
Journal of Supercomputing, 75(12), p.8115 - 8146, 2019/12
被引用回数：0 パーセンタイル：100(Computer Science, Hardware & Architecture)ジャイロ運動論的トロイダル5次元オイラーコードGT5Dにおける反復法線形ソルバの性能向上に向けて省通信一般化最小残差法(CA-GMRES)をCPU-GPUハイブリッドクラスタで実装した。CA-GMRESに加え、計算量を削減するために我々が提案した修正版CA-GMRES(M-CA-GMRES)の実装と評価も行った。本研究から、集団通信回数の最小化と密行列積演算による高効率演算というCA-GMRESの利点が実証された。性能評価は1ノードあたりNVIDIA Tesla P100 GPU4台を搭載したReedbush-L GPUクラスタで実施した。この結果、M-CA-GMRESによりCA-GMRES, 一般化共役残差法(GCR), GMRESに比べてそれぞれ1.09x, 1.22x, 1.50xの高速化が示された。
山田 進
固有値計算と特異値計算; 計算力学レクチャーコース, p.77 - 104, 2019/12
本書「固有値計算と特異値計算」は一般の研究者向けに固有値計算及び特異値計算の基礎理論から応用例までを解説する書籍である。発表者は本書の第3章「疎行列の固有値計算」を執筆し、疎行列向きの固有値計算方法として頻繁に用いられているLanczos法やLOBPCG法などのアルゴリズムを解説している。また、実際にこれらの解法を用いて固有値問題を計算した結果を示すとともに、それぞれの解法の有効性や問題点などを紹介している。
2018年6月までの7年間に渡る月間流出量の推定町田 昌彦; 山田 進; 岩田 亜矢子; 乙坂 重嘉; 小林 卓也; 渡辺 将久; 船坂 英之; 森田 貴己*
日本原子力学会和文論文誌, 18(4), p.226 - 236, 2019/12
2011年45月にかけて発生した東京電力HD・福島第一原子力発電所2号機及び3号機からの汚染水の海洋への直接流出以後、神田は相対的に小さいが連続的な放射性物質の流出が引き続き起こっていることを指摘している。しかし、その期間は2012年9月までであり、その後の流出量の推定についての報告はない。そこで、本論文では、その後を含めて2011年4月から2018年6月までの7年間に渡りCs-137の流出量を推定した結果を報告する。報告のない時期、国・東京電力HDは、流出を抑制するための努力を続け、港湾内海水の放射性核種濃度は徐々に減少している。われわれは、一月単位でCs-137の流出量を二つの手法、一つは神田の提案した手法だがわれわれの改良を加えた手法とボロノイ分割によるインベントリー評価法を使い評価した。それらの結果から、前者の手法は常に後者の手法と比べて保守的だが、前者の後者に対する比は1桁の範囲内であることが分かった。また、それらの推定量から簡単に沿岸域に対するインパクトを評価し、特に魚食による内部被ばく量を推定したところ、1Fの海洋流出量に基づく内部被ばく分は極めて小さいことが分かった。
長尾 郁弥; 新里 忠史; 佐々木 祥人; 伊藤 聡美; 渡辺 貴善; 土肥 輝美; 中西 貴宏; 佐久間 一幸; 萩原 大樹; 舟木 泰智; et al.
JAEA-Research 2019-002, 235 Pages, 2019/08
JAEA-Research-2019-002.pdf:21.04MB
2011年3月11日に発生した太平洋三陸沖を震源とするマグニチュード9.0の東北地方太平洋沖地震とそれに伴って発生した津波により、東京電力福島第一原子力発電所の事故が発生し、その結果、環境中へ大量の放射性物質が放出され、その大部分が森林に沈着している。これに対し、面積が広大であり大量の除去土壌等が生じる、多面的な森林の機能が損なわれる可能性があるなどの問題があり、生活圏近傍を除き、汚染された森林の具体的な除染計画はない。そのため、未除染の森林から放射性セシウムが流出し、既に除染された生活圏に流入することに対する懸念があり、避難住民の帰還や産業再開の妨げとなる可能性があった。原子力機構では、環境中に放出された放射性物質、特に放射性セシウムの移動挙動に関する「長期環境動態研究」を2012年11月より実施している。この目的は、自治体の施策立案を科学的側面から補助する、住民の環境安全に関する不安を低減し、帰還や産業再開を促進するといった点にある。本報告書は、原子力機構が福島県で実施した環境動態研究におけるこれまでの研究成果について取りまとめたものである。
Wang, H.*; 大津 秀暁*; 千賀 信幸*; 川瀬 頌一郎*; 武内 聡*; 炭竃 聡之*; 小山 俊平*; 櫻井 博儀*; 渡辺 幸信*; 中山 梓介; et al.
Communications Physics (Internet), 2(1), p.78_1 - 78_6, 2019/07
被引用回数：4 パーセンタイル：25.43(Physics, Multidisciplinary)陽子(あるいは中性子)過剰核の効率的な生成経路を探索することは、原子核反応研究の主な動機のひとつである。本研究では、
Pdに対する核子当たり50MeVの陽子および重陽子入射による残留核生成断面積を逆運動学法によって測定した。その結果、重陽子入射ではAgやPd同位体の生成断面積が大きくなることを実験的に示した。また、理論計算による解析から、この生成断面積の増大は重陽子の不完全融合反応に起因することを示した。これらの結果は、陽子過剰核の生成において重陽子のような弱束縛核の利用が有効であることを示すものである。
井戸村 泰宏; 伊奈 拓也*; 山下 晋; 小野寺 直幸; 山田 進; 今村 俊幸*
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 16th International Conference on Accelerator and Large Experimental Physics Control Systems (ICALEPCS 2017) (Internet), p.1470 - 1473, 2018/01
J-PARC加速器制御用のネットワークシステムは、2005年以来10年以上運用されてきた。加速器制御ネットワークは、J-PARCの3加速器(リニアック, 3Gevシンクロトロン, MR)と3実験施設(物質・生命科学実験施設,ニュートリノ実験施設,ハドロン実験施設)をカバーし、エッジスイッチの総数は約250である。その全体構成(光ネットワーク網)、冗長機能、VLAN構成、について現状を説明する。次に、オフィスネットワーク(JLAN)と加速器制御ネットワークとの関係について説明する。2つのネットワークの間にFirewallを導入して直接通信を禁止する一方で、双方から制限付きで通信できる別ネットワークを設定した。この手法で、アクセス制限(セキュリティ)と接続性(利便性)の両立を図っている。また、J-PARC加速器制御システムのセキュリティ運用について説明する。オフィスネットワークからのウィルス感染を防ぐための複数のアクセス制限や、アンチウィルスソフトの端末への適用状況について解説する。
山田 進; 今村 俊幸*; 町田 昌彦
Lecture Notes in Computer Science 10776, p.243 - 256, 2018/00
被引用回数：0 パーセンタイル：100本発表は科学研究費補助金(科研費)研究の一環として実施した量子問題計算に現れる大規模な対称疎行列であるハミルトニアン行列の複数の固有値および固有ベクトルを計算する際に用いる省通信ノイマン展開前処理の有効性についての報告である。具体的な成果は、ハミルトニアン行列の複数の固有値とそれに対応する固有ベクトルをLOBPCG法を用いて反復計算する際の収束性を向上させるための新規の前処理方法を提案し、実際の計算から既存の方法よりも短時間で計算できることを示したことである。さらに、問題の物理的性質を利用して前処理計算の通信回数を減少させるアルゴリズムを提案し、原子力機構のスパコンSGI ICEXおよび理化学研究所の京スパコンの数千コアを用いた並列計算において、高速化が実現することを確認した。この成果は、量子計算の高速計算だけでなく、今後大規模化していく並列計算機の有効利用に資する成果でもある。
井戸村 泰宏; 伊奈 拓也*; 真弓 明恵; 山田 進; 今村 俊幸*
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倍高速であることを示した。
山田 進; 今村 俊幸*; 町田 昌彦
Parallel Computing is Everywhere, p.27 - 36, 2018/00
本発表は科研費研究の一環として実施した固有値計算ソルバLOBPCG法の前処理に関する研究についての発表である。LOBPCG法は反復解法であり適切な前処理を用いることで収束性が向上することが知られているが、量子問題に表れるハバードモデルの固有値を計算する際に物理パラメータによっては既存の前処理では収束性が向上しないことがある。そこで、線形方程式の反復解法で前処理として使われているノイマン展開を利用した前処理を適用し、そのような問題に対しても収束性が向上することを見出した。さらに、問題の物理的性質を考慮し、演算回数は若干増加するが、通信回数を減少させる方法を提案した。この方法で開発したコードを、原子力機構の並列計算機SGI ICEXを用いて並列シミュレーションしたところ、これまでの方法よりも約20%高速に計算できることを確認した。この成果は、固有値問題に有効な前処理方法を見出しただけではなく、ネットワーク構造が複雑化してきている並列計算機での性能向上にも資する成果である。
井戸村 泰宏; 伊奈 拓也*; 真弓 明恵; 山田 進; 松本 和也*; 朝比 祐一*; 今村 俊幸*
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%に削減された。
山田 進; 伊奈 拓也*; 佐々 成正; 井戸村 泰宏; 町田 昌彦; 今村 俊幸*
Proceedings of 2017 IEEE International Parallel & Distributed Processing Symposium Workshops (IPDPSW) (Internet), p.1418 - 1425, 2017/08
被引用回数：3 パーセンタイル：21.83本発表では、BaileyのDouble-Doubleアルゴリズムを利用した4倍精度基本線形代数演算ライブラリ「BLAS」の高速化と、それを4倍精度固有値計算ルーチンに適用した際の計算性能の向上について発表する。特筆すべき成果は、積和演算の中間結果を高精度で保持して計算できるFMA命令を利用することで、我々がこれまでに開発した4倍精度BLASと比較し、代表的なルーチンで2045%の高速化を実現したことである。さらに、4倍精度固有値ソルバ「QPEigenK」が利用している4倍精度BLASをFMA命令を利用して開発したものに置き換えるだけで、原子力機構のICE Xおよび東京大学のFX10の2つのスパコンの大規模並列計算において1020%の高速化を実現した。この成果は、計算機の大規模化に伴って必要になりつつある高精度計算の高速な実行に資する成果である。
小林 恵太*; 奥村 雅彦; 山田 進; 町田 昌彦; 青木 秀夫*
Physical Review B, 94(21), p.214501_1 - 214501_7, 2016/12
被引用回数：24 パーセンタイル：16.9(Materials Science, Multidisciplinary)平坦バンドが超伝導性を示す可能性を探るために、平坦バンドを形成する最も単純な準一次元系の一つであるダイアモンド鎖上の斥力相互作用するフェルミオン系について調べた。厳密対角化法と密度行列繰り込み群法を用いて調べた結果、フェルミエネルギーに近い空の平坦バンドと相互作用する分散バンドが満たされる1/3フィリングよりも少しだけ小さなフィリングで、長い相関距離を持つクーパー対が有意な束縛エネルギーを持つことがわかった。さらに、この対相関関数は、ダイアモンド鎖の外側のサイトに存在するフェルミオン対によるものであることを明らかにした。また、厳密に1/3フィリングの時、系は絶縁体になり、ダイアモンド鎖の外側のサイトに存在するフェルミオンがトポロジカルに区別可能なエンタングル状態を形成していることがわかった。
真弓 明恵; 井戸村 泰宏; 伊奈 拓也; 山田 進; 今村 俊幸*
Proceedings of 7th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA 2016) (Internet), p.17 - 24, 2016/11
左前処理省通信共役勾配(LP-CA-CG)法を多相数値流体力学コードJUPITERの圧力Poisson方程式に適用した。LP-CA-CG法の演算密度を分析し、内積処理と三項間漸化式処理のループ分割を行うことで演算密度を大きく向上した。ブロックヤコビ前処理及びアンダーラップ前処理を適用した2つのLP-CA-CGソルバを開発した。京コンピュータ上では局所的な1対1通信のスケールが良好であることと、アンダーラップ前処理を適用すると収束性が悪くなることにより、ブロックヤコビ前処理ソルバにより良好な性能が得られた。このソルバは3万ノードまで良好な強スケーリングを示し、大域的集団通信のコストを69%削減することにより従来のCG法ソルバに比べて高い性能を達成した。
