Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
山田 進
固有値計算と特異値計算; 計算力学レクチャーコース, 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
被引用回数：3 パーセンタイル：17.89(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
本発表は科学研究費補助金(科研費)研究の一環として実施した量子問題計算に現れる大規模な対称疎行列であるハミルトニアン行列の複数の固有値および固有ベクトルを計算する際に用いる省通信ノイマン展開前処理の有効性についての報告である。具体的な成果は、ハミルトニアン行列の複数の固有値とそれに対応する固有ベクトルを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
本発表では、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
被引用回数：21 パーセンタイル：16.01(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法ソルバに比べて高い性能を達成した。
佐々 成正; 山田 進; 町田 昌彦; 今村 俊幸*
Nonlinear Theory and Its Applications, IEICE (Internet), 7(3), p.354 - 361, 2016/07
フールエ変換を連続使用した場合の丸め誤差の累積誤差について議論を行った。偏微分方程式の数値計算において、フールエ変換を繰り返し連続使用した場合に丸め誤差が累積することを確認した。計算精度の毀損を補うため4倍精度高速フーリエ変換を用いた計算を行って有効性を確認した。
山田 進; 今村 俊幸*; 町田 昌彦
Parallel Computing; On the Road to Exascale, p.361 - 369, 2016/00
本発表では簡単な演算を高速に実行できるアクセラレータ(補助演算装置)の1つであるGPUを用いてハバードモデルの計算に現れる固有値問題を高速に計算する手法について発表する。特筆すべき成果は、ハバードモデルの物理的性質、およびGPUの構造を考慮してデータの格納形式や計算方法を考案したところ、GPUで一般的に利用されている計算方法よりも約2倍高速に計算できることを実際のGPUを利用した計算から示したことである。さらに、6コアのCPUを用いた並列計算と比較しても、23倍高速に計算できることを確認した。この成果は、今後高性能計算において主流になると考えられているアクセラレータを利用した計算機の有効利用に資する成果である。
湯山 貴裕; 石堀 郁夫; 倉島 俊; 吉田 健一; 石坂 知久; 千葉 敦也; 山田 圭介; 横山 彰人; 薄井 絢; 宮脇 信正; et al.
Proceedings of 12th Annual Meeting of Particle Accelerator Society of Japan (インターネット), p.302 - 304, 2015/09
日本原子力研究開発機構のイオン照射施設TIARAでは4台の加速器により、材料・バイオ技術の研究開発への利用を主として、広範囲のエネルギー及び多様なイオン種のビームを提供している。本発表では2014年度のTIARAの稼働状況、保守・整備及び技術開発を報告する。保守・整備及び技術開発の主要な内容を以下に示す。サイクロトロンの高周波系において、ショート板用接触子に焼損が発生したため、接触子の交換及び焼損箇所の研磨を行うことで復旧させた。原因調査の結果、経年劣化によりフィードバックケーブルが断線しかかっていたため、不必要な高電圧が印加されたことが原因と判明した。サイクロトロン制御システムに関して、サポートが停止されたWindows XPをWindows 7に変更し、これに伴い制御システムを更新するとともに、トレンドグラフのログデータ保存機能、操作画面上の制御対象一括選択機能の付加など、各種機能を向上させた。C
イオンビームの計測に関して、複雑な二次荷電粒子を生成するCイオンビームの正確な電流測定のために、サプレッサー電極の構造を改良することで二次荷電粒子を十分捕集するファラデーカップを開発した。
操上 広志; 北村 哲浩; 山田 進; 町田 昌彦
Proceedings of 23rd International Conference on Nuclear Engineering (ICONE-23) (DVD-ROM), 6 Pages, 2015/05
福島環境中の放射性セシウムの移動に係る様々な時間・空間スケールの課題に対応するため、複数の数値モデルが準備された。本論文はいくつかのケーススタディとともに、環境中での放射性セシウムの移動に係る諸課題に対する原子力機構のモデル化アプローチの一部について記述する。
山田 進; 北村 哲浩; 操上 広志; 町田 昌彦
Proceedings of Joint International Conference on Mathematics and Computation, Supercomputing in Nuclear Applications and the Monte Carlo Method (M&C + SNA + MC 2015) (CD-ROM), 9 Pages, 2015/04
本発表は福島長期環境動態研究(F-TRACE)の一環として実施した2次元河川シミュレーションコードNays2Dの並列化に関する報告である。放射性セシウムは粒径の小さい土砂に吸着し、河川など流水の影響を受け土砂と一緒に移動する。そのため、河川による土砂の移動をシミュレーションすることでセシウムの振る舞いを予測することができる。河川の流れは地形の影響を大きく受けるため、詳細な土砂の堆積位置を知るためには2次元モデルのシミュレーションが必要になるが、計算量が多く広い範囲のシミュレーションを実施することが困難である。そこで、2次元河川シミュレーションコードNays2Dを並列化した。特筆すべき成果は、この並列化により、原子力機構のBX900を用いた並列計算で10倍以上の高速化を実現し、実際に、大柿ダム内の粘土の振る舞いのシミュレーションを実施し、現実的な時間でシミュレーションすることに成功したことである。この結果は、並列化により、現実的な時間で詳細な土砂の振る舞いを予測することが可能になったことを示しており、計算科学だけでなく、環境動態研究にも資する成果である。
Cs transport and accumulation in the Ogaki Dam of eastern Fukushima山田 進; 北村 哲浩; 操上 広志; 山口 正秋; Malins, A.; 町田 昌彦
Environmental Research Letters, 10(1), p.014013_1 - 014013_9, 2015/01
被引用回数：15 パーセンタイル：34.93(Environmental Sciences)本論文は福島長期環境動態研究(F-TRACE)の一環として実施した2次元河川シミュレーションコードNays2Dによる福島県の大柿ダムでの放射性セシウムの流動及び堆積分布予測に関する研究についての報告である。現状のシミュレーションコードでは放射性物質の移動は計算できないが、放射性セシウムは粘土のような粒径の小さい土砂に吸着しやすいことに着目し、土砂の移動を計算しそれを放射性セシウムの移動とみなすことにした。実際に、2013年9月の洪水時や平均的な洪水時の流量等を境界条件として用いて、大柿ダム内の土砂の流れのシミュレーションを原子力機構のBX900で実施し、土砂の移動を計算した。特筆すべき成果は、粒径の小さいシルト及び粘土の振る舞いはダムの排水口の高さに大きく依存しており、排水口の位置を高くすることでダム内にシルト及び粘土が多くとどまるようになることを示したことである。この結果は、ダムの水位を適切に調整することで下流への土砂に付着したセシウムの移動を防げる可能性があることを示しており、住民の被ばくの低減に資する成果である。
北村 哲浩; 操上 広志; 山口 正秋; 小田 好博; 齋藤 龍郎; 加藤 智子; 新里 忠史; 飯島 和毅; 佐藤 治夫; 油井 三和; et al.
Nuclear Science and Engineering, 179(1), p.104 - 118, 2015/01
被引用回数：6 パーセンタイル：35.56(Nuclear Science & Technology)福島第一原子力発電所事故に伴い環境に放出されその後地表に降下した放射性物質の分布を予測することは重要で、速やかに進めて行く必要がある。このような予測を行うために、放射性物質として特に放射性セシウムに着目し、現在複数の数理モデルを開発している。具体的には、土壌の表層流出に伴う放射性セシウムの移行については土壌流亡予測式を用いた流出解析、河川における核種移行については河川解析コードTODAM・iRICを用いた移行解析、河口域における土砂堆積については3次元解析コードROMS等を応用した堆積解析を行っている。また、セシウムと土壌の吸着メカニズムについては分子原子レベルの分子挙動計算法を用いた解析を開始しており、最終目標として吸着係数等の把握を目指している。
