Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
前山 伸也*; 渡邉 智彦*; 仲田 資季*; 沼波 政倫*; 朝比 祐一; 石澤 明宏*
Nature Communications (Internet), 13, p.3166_1 - 3166_8, 2022/06
被引用回数：19 パーセンタイル：95.69(Multidisciplinary Sciences)乱流輸送は、磁場閉じ込め核融合プラズマを閉じ込めるための重要な物理過程である。最近の理論的,実験的研究により 小さい(電子)スケールと大きい(イオン)スケールの乱流の間にクロススケールの相互作用が存在することが明らかにされている。従来の乱流輸送モデルではクロススケール相互作用が考慮されていないため、将来の核燃焼プラズマ実験においてクロススケール相互作用を考慮する必要があるかどうかを明らかにする必要がある。核燃焼プラズマ実験では、核融合で生まれたアルファ粒子によって高い電子温度が維持されるため、プラズマの性質が今まで実験されてきたものと大きく異なると予測される。本論文では、スーパーコンピュータを用いたシミュレーションにより、高電子温度プラズマにおける電子スケールの 電子温度プラズマの乱流は、電子だけでなく燃料や灰の乱流輸送にも影響を与えることを明らかにした。電子スケールの乱流は、イオンスケールの微小的不安定性の原因である共鳴電子の軌道を乱し、大きなスケールの乱流揺らぎを抑制する。同時に、イオンスケールの乱流渦も電子スケールの乱流を抑制する。これらの結果は 異なるスケールの乱流が互いに排他的であることを示す。また、クロススケール相互作用により、熱流束が減少する可能性を示す。
前山 伸也*; 渡邉 智彦*; 井戸村 泰宏; 仲田 資季*; 沼波 政倫*
Computer Physics Communications, 235, p.9 - 15, 2019/02
被引用回数：6 パーセンタイル：47.95(Computer Science, Interdisciplinary Applications)陰的時間積分スキームを用いてSugama衝突演算子をジャイロ運動論的ブラゾフコードGKVに実装した。新手法はオペレータスプリッティング、陰的時間積分、クリロフ部分空間反復法ソルバを用いており、線形化衝突演算子の詳細に依存しない汎用的な手法となっている。数値計算テストでは衝突項が制限する時間ステップ幅を超えて安定な計算が実証できた。データ転置を用いることで、反復計算中に通信が発生しない分散メモリシステムにおける効率的実装を実現した。これにより、本手法は計算効率向上と計算コスト削減を同時に達成し、アプリケーションの全体性能を大きく加速する。
菅野 龍太郎*; 沼波 政倫*; 佐竹 真介*; 松岡 清吉; 高丸 尚教*
Nuclear Fusion, 58(1), p.016033_1 - 016033_7, 2018/01
被引用回数：0 パーセンタイル：0.00(Physics, Fluids & Plasmas)摂動磁場による磁気面の破壊を伴う乱れた磁場構造をもつトーラスプラズマ中の電子熱輸送について研究を行った。本論文では、摂動磁場に沿った粒子の平行方向軌道に起因する熱輸送が、捕捉粒子軌道の効果によって低減されることを明らかにした。
前山 伸也*; 渡邉 智彦*; 井戸村 泰宏; 仲田 資季*; 石澤 明宏*; 沼波 政倫*
Nuclear Fusion, 57(6), p.066036_1 - 066036_10, 2017/05
被引用回数：17 パーセンタイル：66.93(Physics, Fluids & Plasmas)電子およびイオン温度勾配(ETG/ITG)モードを含むマルチスケールプラズマ乱流を電磁的ジャイロ運動論シミュレーションによって調べた。非線形モード結合の3波移行解析により、電子スケールとイオンスケールのスケール間相互作用を明らかにした。この相互作用の1つは、イオンスケール乱流による電子スケール乱流の抑制である。ここで、ITGが駆動する短波長の渦がせん断流れのように作用し、ETG乱流を抑制する。もう一つのスケール間相互作用は、電子スケール乱流の存在下でのイオンスケール乱流の増大である。これは、短波長の帯状流によって引き起こされる。この帯状流はITG乱流における運動論的通過電子の応答によって生成され、剪断によってITG乱流を抑制し、ETG乱流によって減衰される。どちらの場合も、電子スケールとイオンスケールの間のサブイオンスケールの構造が、スケール間相互作用において重要な役割を果たす。
菅野 龍太郎*; 沼波 政倫*; 佐竹 真介*; 松岡 清吉; 高丸 尚教*
Contributions to Plasma Physics, 56(6-8), p.592 - 597, 2016/08
被引用回数：2 パーセンタイル：10.88(Physics, Fluids & Plasmas)共鳴摂動磁場と径電場の効果を含む準定常状態の衝突性輸送を計算するための
法に基づくドリフト運動論コードを開発した。本論文では、コードの妥当性を複数のテスト問題によって確認した。その結果、径電場がない場合に磁力線のカオス的振る舞いによって強く影響を受ける電子輸送が、正の径電場によって低減されることを明らかにした。
石澤 明宏*; 井戸村 泰宏; 今寺 賢志*; 糟谷 直宏*; 菅野 龍太郎*; 佐竹 真介*; 龍野 智哉*; 仲田 資季*; 沼波 政倫*; 前山 伸也*; et al.
プラズマ・核融合学会誌, 92(3), p.157 - 210, 2016/03
幅広いアプローチ協定に基づいて国際核融合エネルギー研究センター(IFERC)の計算機シミュレーションセンター(CSC)に設置された高性能計算機システムHeliosは、2012年1月に運用を開始し、日欧の磁気核融合シミュレーション研究に供用され、高い利用率の実績を示すとともに、炉心プラズマ物理から炉材料・炉工学にわたる広い分野で多くの研究成果に貢献している。本プロジェクトレビューの目的は、国内の大学や研究機関においてHeliosを利用して進められているシミュレーション研究プロジェクトとその成果を一望するとともに、今後予想される研究の進展を紹介することである。はじめにIFERC-CSCの概要を示した後、各研究プロジェクト毎にその目的、用いられる計算手法、これまでの研究成果、そして今後必要とされる計算を紹介する。
前山 伸也; 渡邉 智彦*; 井戸村 泰宏; 仲田 資季; 沼波 政倫*; 石澤 明宏*
Parallel Computing, 49, p.1 - 12, 2015/11
被引用回数：7 パーセンタイル：51.36(Computer Science, Theory & Methods)Optimization techniques of a plasma turbulence simulation code GKV for improved strong scaling are presented. This work is motivated by multi-scale plasma turbulence extending over multiple spatio-temporal scales of electrons and ions, whose simulations based on the gyrokinetic theory require huge calculations of five-dimensional (5D) computational fluid dynamics by means of spectral and finite difference methods. First, we present the multi-layer domain decomposition of the multi-dimensional and multi-species problem, and segmented MPI-process mapping on 3D torus interconnects, which fully utilizes the bi-section bandwidth for data transpose and reduces the conflicts of simultaneous point-to-point communications. These techniques reduce the inter-node communication cost drastically. Second, pipelined computation-communication overlaps are implemented by using the OpenMP/MPI hybrid parallelization, which effectively mask the communication cost. Thanks to the above optimizations, GKV achieves excellent strong scaling up to 
600 k cores with high effective parallelization rate 99.99994% on K, which demonstrates its applicability and efficiency toward a million of cores. The optimized code realizes multi-scale plasma turbulence simulations covering electron and ion scales, and reveals cross-scale interactions of electron- and ion-scale turbulence.
前山 伸也*; 井戸村 泰宏; 渡邉 智彦*; 仲田 資季*; 矢木 雅敏; 宮戸 直亮; 石澤 明宏*; 沼波 政倫*
Physical Review Letters, 114(25), p.255002_1 - 255002_5, 2015/06
被引用回数：98 パーセンタイル：95.25(Physics, Multidisciplinary)実イオン・電子質量比かつ実
値のマルチスケールジャイロ運動論的乱流シミュレーションをはじめて実現した。ここで、値はプラズマ圧力と磁気圧の比によって与えられ、微視的不安定性に対する電磁効果を特徴付ける。電子・イオンスケールのスケール間相互作用を明らかにするために両スケールにおける数値解析を行った。実質量比のスケール分離があっても、イオンスケール乱流が電子スケールストリーマを消去し、イオンだけでなく電子の熱輸送も支配する。また、イオンスケールモードが有限効果によって安定化される場合には乱流輸送への電子スケールダイナミクスの寄与は無視できなくなり、イオンスケール乱流輸送を増大させることがわかった。
渡邉 智彦*; 井戸村 泰宏; 前山 伸也; 仲田 資季; 洲鎌 英雄*; 沼波 政倫*; 石澤 明宏*
Journal of Physics; Conference Series, 510, p.012045_1 - 012045_11, 2014/05
被引用回数：2 パーセンタイル：48.33(Computer Science, Interdisciplinary Applications)配位空間と速度空間を含む位相空間において分布関数と電磁場の揺らぎを伴うプラズマ乱流が発展する。磁場核融合プラズマにおける分布関数の詳細構造を超並列計算機上で実行されるジャイロ運動論シミュレーションによって明らかにした。ドリフト波乱流のジャイロ運動論シミュレーションによって位相空間におけるエントロピー移送、ヘリカル磁場による帯状流の増大と輸送低減を示した。最先端の高性能数値計算をイオンスケールと電子スケールを含むマルチスケール乱流シミュレーション、および、ITER規模の大域的乱流輸送シミュレーションに適用した。
仲田 資季; 松山 顕之; 相羽 信行; 前山 伸也; 沼波 政倫*; 渡邉 智彦*
Plasma and Fusion Research (Internet), 9, p.1403029_1 - 1403029_12, 2014/04
局所ジャイロ運動論乱流シミュレーションコードGKVへ、MHD平衡コードMEUDASから計算される実トカマクMHD平衡データが新たに組み込まれた。拡張されたコードは他の計算コードとのベンチマークテストにより実装の健全性検証が行われ、ITG不安定性の分散関係の良い一致が確認された。さらに、本コードを用いて、JT-60SAトカマク装置で期待されている2種類のプラズマ平衡[ITER-like平衡およびHighly-shaped平衡]に対するITG不安定性解析および残留ゾーナルフロー解析を実施した。計算結果の詳細な比較により、ITER-like平衡に比べHighly-shaped平衡の方が、相対的に低いITG不安定性成長率と高い残留ゾーナルフローレベルを示すことが明らかとなり、より高い安定性特性を示すことが確認された。
井戸村 泰宏; 仲田 資季; 山田 進; 町田 昌彦; 今村 俊幸*; 渡邉 智彦*; 沼波 政倫*; 井上 晃*; 堤 重信*; 三吉 郁夫*; et al.
International Journal of High Performance Computing Applications, 28(1), p.73 - 86, 2014/02
被引用回数：17 パーセンタイル：75.01(Computer Science, Hardware & Architecture)A plasma turbulence research based on five dimensional (5D) gyrokinetic simulations is one of the most critical and demanding issues in fusion science. To pioneer new physics regimes both in problem sizes and in time scales, an improvement of strong scaling is essential. In this work, the strong scaling and the parallel efficiency of a gyrokinetic toroidal 5D Eularian code GT5D is dramatically improved by novel developing communication overlap techniques. On the K-computer, excellent strong scaling is achieved beyond 100 k cores with keeping the sustained performance of 10% ( 307 TFlops using 196,608 cores), and simulations for next generation large-scale fusion experiments are significantly accelerated.
前山 伸也; 渡邉 智彦*; 井戸村 泰宏; 仲田 資季; 沼波 政倫*; 石澤 明宏*
Plasma and Fusion Research (Internet), 8, p.1403150_1 - 1403150_8, 2013/11
We extend our gyrokinetic Vlasov simulation code GKV to run efficiently on peta-scale supercomputers. A key numerical technique is the parallel Fast Fourier Transform (FFT) required for parallel spectral calculations, where masking of the cost of inter-node transpose communications is essential to improve strong scaling. To mask communication costs, computation-communication overlap techniques are applied for FFTs and transpose with the help of the hybrid parallelization of message passing interface and open multi-processing. Integrated overlaps including whole spectral calculation procedures show better scaling than simple overlaps of FFTs and transpose. The masking of communication costs significantly improves strong scaling of the GKV code, and makes substantial speed-up toward multi-scale turbulence simulations.
flows in helical plasmas渡邉 智彦*; 沼波 政倫*; 洲鎌 英雄*; 佐竹 真介*; 松岡 清吉*; 石澤 明宏*; 前山 伸也; 田中 謙治*
Proceedings of 24th IAEA Fusion Energy Conference (FEC 2012) (CD-ROM), 8 Pages, 2012/10
Gyrokinetic simulation of ion temperature gradient turbulence and zonal flows for helical plasmas has been validated against the Large Helical Device experiments with high ion temperature, where a reduced modeling of ion heat transport is also considered. It is confirmed by the entropy transfer analysis that the turbulence spectrum elongated in the radial wavenumber space is associated with successive interactions with zonal flows. A novel multi-scale simulation for turbulence and zonal flows in poloidally-rotating helical plasmas has demonstrated strong zonal flow generation byturbulence, which implies that turbulent transport processes in non-axisymmetric systems are coupled to neoclassical transport through the macroscopic flows determined by the ambipolarty condition for neoclassical particle fluxes.
井戸村 泰宏; 仲田 資季; 山田 進; 町田 昌彦; 今村 俊幸*; 渡邉 智彦*; 沼波 政倫*; 井上 晃*; 堤 重信*; 三吉 郁夫*; et al.
Proceedings of 31st JSST Annual Conference; International Conference on Simulation Technology (JSST 2012) (USB Flash Drive), p.234 - 242, 2012/09
A plasma turbulence research based on 5D gyrokinetic simulations is one of the most critical and demanding issues in fusion science. To pioneer new physics regimes both in problem sizes and in time scales, an improvement of strong scaling is essential. Overlap of computations and communications is a promising approach to improve strong scaling, but it often fails on practical applications with conventional MPI libraries. In this work, this classical issue is clarified, and resolved by developing communication overlap techniques with mpi_test and communication threads, which work even on conventional MPI libraries and network hardwares. These techniques dramatically improve the parallel efficiency of a gyrokinetic Eularian code GT5D on K and Helios, which adopt dedicated and commodity networks.
西原 功修*; 砂原 淳*; 佐々木 明; 沼波 政倫*; 田沼 肇*; 藤岡 慎介*; 島田 義則*; 藤間 一美*; 古河 裕之*; 加藤 隆子*; et al.
Physics of Plasmas, 15(5), p.056708_1 - 056708_11, 2008/00
被引用回数：132 パーセンタイル：97.65(Physics, Fluids & Plasmas)Extreme ultraviolet (EUV) radiation from laser-produced plasma (LPP) has been thoroughly studied for application in mass-production of the next generation semiconductor devices. One critical issue for realization of a LPP-EUV light source for lithography is the conversion efficiency (CE) from incident laser power to EUV radiation of 13.5 nm wavelength (within 2% bandwidth). Another is solving a problem of damage caused by debris reaching a EUV collecting mirror. We here present an improved power balance model, which can be used for the optimization of laser and target conditions to obtain high CE. An integrated numerical simulation code has been developed for target design. The code is shown to agree well with experimental results not only for the CE but also for detailed EUV spectral structure. We propose a two pulse irradiation scheme for high CE and reduced ion debris using a carbon dioxides laser and a droplet or punch-out target.
前山 伸也; 渡邉 智彦*; 井戸村 泰宏; 仲田 資季; 沼波 政倫*; 石澤 明宏*
no journal, ,
To analyze multi-scale plasma turbulence in magnetically confined plasma, we have extended the gyrokinetic Vlasov simulation code GKV to run efficiently on the K computer. A key numerical technique is the speed up of the parallel two-dimensional (2D) Fast Fourier Transform (FFT) algorithms. MPI/OpenMP hybrid parallelization with 4D domain decomposition is newly implemented to use a large number of cores. Perpendicular dynamics is solved by using the spectral method with 2D FFTs, where the data transpose often degrades scalability. The MPI process mapping is optimized for the 3D torus inter-node network on the K computer, so that the transpose communications are performed among the neighboring nodes located in a 3D box shape. Additionally, overlaps of FFTs and data transpose are implemented by using OpenMP directives. Thanks to the developed methods, GKV achieves almost linear speed-up beyond 100k cores with high parallel efficiency about 99.9998 per cent.
前山 伸也; 渡邉 智彦*; 井戸村 泰宏; 仲田 資季; 矢木 雅敏; 宮戸 直亮; 沼波 政倫*; 石澤 明宏*
no journal, ,
国際熱核融合実験炉ITERでは核燃焼プラズマ実験が計画されており、核融合生成高エネルギーアルファ粒子により主に加熱される電子の熱輸送機構の早急な解明が求められている。本研究では、電子熱輸送解析においてこれまで十分に考慮されてこなかった微視的乱流のマルチスケール性について、ジャイロ運動論的シミュレーションに基づく解析を行った。これにより、従来スケール分離が仮定されてきたイオン温度勾配駆動(ITG)乱流と電子温度勾配駆動(ETG)乱流の間で、乱流の非線形結合を介したマルチスケール相互作用が存在すること、および、それによる電子熱輸送スペクトルが大きく変化することを明らかにした。
仲田 資季; 本多 充; 吉田 麻衣子; 浦野 創; 前山 伸也; 沼波 政倫*; 渡邉 智彦*
no journal, ,
ITERやDEMOに代表される大規模な燃焼プラズマの輸送特性の解析・予測においては、核融合プラズマの第一原理モデルであるジャイロ運動論に基づいた乱流輸送シミュレーションが有効であり、その定量予測性能の検証が近年の課題となっている。本研究では、乱流シミュレーションコードGKV-Jへ実トカマク形状効果を新たに導入し、JT-60Uトカマクに対する第一原理的乱流輸送シミュレーションを実現した。イオン種に加え、通過/捕捉電子も運動論的に取り扱うことにより、半径位置 に応じた微視的不安定性の同定が可能となり、それらが駆動する乱流中のゾーナルフロー形成過程や熱・粒子に関する乱流輸送特性が明らかとなった。拡張されたGKV-Jによる乱流シミュレーション結果は、従来の流体輸送モデルでは定量性を著しく欠いていた領域に対しても、イオン/電子熱輸送レベルに関する実験結果をよく再現することに成功している。加えて、JT-60SA配位での乱流輸送における形状度の影響を同定することにも成功した。
仲田 資季; 本多 充; 吉田 麻衣子; 浦野 創; 前山 伸也; 沼波 政倫*; 渡邉 智彦*
no journal, ,
ITERやDEMOに代表される次世代の燃焼プラズマの乱流輸送特性の解明と予測には、核融合プラズマの第一原理モデルであるジャイロ運動論に基づいた乱流輸送シミュレーションが不可欠である。本研究では、著者らが開発している第一原理乱流シミュレーションコードGKV-Jを駆使して、実際のJT-60Uトカマク平衡を用いたITG-TEM駆動乱流輸送シミュレーションを実現した。まず、微視的線形安定性解析により、プラズマ内側領域ではITGモード、中間領域ではITG-TEMモード、外側領域ではTEMモードが支配的となることを同定した。次に非線形乱流輸送シミュレーションを実施し、GKV-Jのシミュレーション結果が従来の流体輸送モデルでは定量性を著しく欠いていた領域に対しても、イオン/電子熱輸送レベルに関する実験結果をよく再現することを確認した。また、詳細な乱流データの解析により、イオン/電子熱輸送、粒子輸送といったそれぞれの輸送チャンネルに対して異なる乱流/ゾーナルフロー強度の依存性が存在することが明らかになった。
井戸村 泰宏; 仲田 資季; 山田 進; 町田 昌彦; 今村 俊幸*; 渡邉 智彦*; 沼波 政倫*; 井上 晃*; 堤 重信*; 三吉 郁夫*; et al.
no journal, ,
A plasma turbulence research based on 5D gyrokinetic simulations is one of the most critical and demanding issues in fusion science. To pioneer new physics regimes both in problem sizes and in time scales, an improvement of strong scaling is essential. Overlap of computations and communications is a promising approach to improve strong scaling, but it often fails on practical applications with conventional MPI libraries. In this work, this classical issue is clarified, and resolved by developing communication overlap techniques with mpi_test and communication threads. These techniques dramatically improve the parallel efficiency of a gyrokinetic Eularian code GT5D on the K-computer, and ITER-size simulations are significantly accelerated.
