学際大規模情報基盤共同利用・共同研究拠点 平成 27 年度共同研究 最終報告書 2016 年 5 月

jh150005-NA04

核融合プラズマ研究のための超並列粒子シミュレーションコード開発とその可視化
大谷寛明（核融合科学研究所）
概要 核融合プラズマシミュレーション研究において、粒子モデルは第一原理的な計算が可能なため粒子の運動論的効果を解析するために活用されている。しかし、装置全体を扱う系への適用は、膨大な数の粒子運動の計算が必要なため困難である。本研究課題
は超並列に対応したプラズマ粒子コードを開発して、微視的物理が装置全体に影響を与えるような現象に対して、第一原理からその発生機構の解明を目指す。この微視的なシ
ミュレーションで得られた知見から微視的現象を表すモデルを開発して磁気流体力学モデルのような巨視的モデルに組み込み、経験的モデルより予測性の高い巨視的シミュレ
ーションの実現を目指す。このような大規模系のポストプロセスとしての可視化も困難なため、シミュレーションと同時に可視化も行う In-situ 可視化の研究開発も行う。

1. 共同研究に関する情報ンとの比較，微視的物理のモデル化(1) 共同研究を実施した拠点名

小川智也（北里大）：流体シミュレーションと

本課題は、核融合科学研究所（核融合研）、神の比較，微視的物理のモデル化

戸大学、情報通信通信研究機構（情報通信研）、

中島 浩（京大）：動的負荷分散ライブラリの北里大学、京都大学、名古屋大学の複数の研開発・実装支援究者が参加する共同研究であり、名古屋大学

深沢圭一郎（京大）：流体シミュレーションとの大規模計算資源を利用した。

の比較，微視的物理のモデル化(2) 共同研究分野


荻野正雄（名大）： FX100 向け最適化・実行

超大規模数値計算系応用分野

支援

2. 研究の目的と意義(3) 参加研究者の役割分担

大谷寛明（核融合研）：総括・粒子コード開発・

核融合プラズマのシミュレーション研究

可視化研究

では、磁気流体力学(MHD)モデルやジャイロ運大 野 暢 亮 （ 兵 庫 県 立 大 ）： そ の 場 可 視 化

動論モデル、粒子モデルなどが使われている。

(In-situ 可視化)の研究開発

粒子モデルは荷電粒子の運動方程式と宇佐見俊介（核融合研）：粒子コードの移植・

Maxwell 方程式を直接解くため、MHD モデルや高速化

ジャイロ運動論モデルでは取り扱えない微視沼波政倫（核融合研）：粒子コードの移植・高的物理を解くことができ、第一原理的なシミ

速化ュレーションが可能である。しかし、巨視的長谷川裕記（核融合研）：粒子コードの高速化

な時空間スケールを扱うには、膨大な数の粒三浦英昭（核融合研）：流体シミュレーション

子運動を計算する必要がある (具体的には数との比較，微視的物理のモデル化

十億～数百億)。

堀内利得（核融合研）：問題設定，結果の評価

トカマク型装置で発生しているディスラプ

臼井英之（神戸大）：動的負荷分散ライブラリション現象や鋸歯状振動などはプラズマ閉じの活用支援

込めの崩壊を招いたり閉じ込め性能を悪化さ三宅洋平（神戸大）：動的負荷分散ライブラリせたりするため、今後の核融合プラズマ研究の実装支援の課題の一つと考えられている。これらの現田 光江（情報通信研）：流体シミュレーショ

象は、粒子の運動論的効果が重要な役割を果1

学際大規模情報基盤共同利用・共同研究拠点 平成 27 年度共同研究 最終報告書 2016 年 5 月

たす磁気リコネクションが引き金となって発

この課題の研究対象である磁気リコネクシ

生していると考えられ、微視的物理を正確にョン現象は太陽フレアや地球磁気圏サブスト

取り扱いつつ装置全体にわたる巨視的な現象

ームのような爆発現象、降着円盤やブラック

をシミュレーションする必要がある。

ホール形成に伴うジェットなどでも重要な役

本研究課題では次のような課題を挙げて割を果たすと考えられている。本研究課題に

粒子モデルを基礎としたシミュレーションコ

より、ディスラプションや鋸歯状振動の発生

ード群の整備を行い、トーラス系プラズマに機構の解明が進めばプラズマ閉じ込め性能のおける微視的物理が装置全体に影響を与える

改善につながり核融合プラズマ研究の飛躍的ような現象に対して、第一原理からその発生な発展が期待でき、さらには天体現象といっ機構の解明を目指し、核融合プラズマの閉じた多岐多様にまたがる多くの分野に貢献する

込め性能改善への寄与を図る。効率的なシミことができる。

ュレーションコードの開発からモデリング、

本課題による粒子シミュレーションの計算可視化までをひとつのパッケージとした、総結果を使えば、第一原理に基づく微視的物理合的な研究基盤の構築を目指す。

現象のモデルを提供することができ、それを①超並列に対応した粒子コードの開発：核融組み込むことでより精密な流体シミュレーシ

合プラズマでは粒子分布が非一様なので、ションが実行できると期待される。このように

ミュレーション領域を単純に領域分割して並モデリングという面からも、宇宙天気予報、

列化しても計算効率は上がらない。超並列に銀河形成、核融合プラズマなどで行われてい対応したアルゴリズム開発・研究を行う。

る流体シミュレーションに対する非常に大き②大規模シミュレーションデータ解析のためな波及効果が期待できる。

の「In-situ 可視化(その場可視化)」の研究

3. 当拠点公募型共同研究として実施した意義

開発：大規模シミュレーションではストレー

本研究課題の推進には、計算科学分野・計ジの制限などから計算したデータを全て保存

算機科学分野・可視化情報学分野の協調的相することは不可能であり、また、解析用マシ

補的な研究体制が必要である。具体的には、

ンのメモリ不足から全データを可視化解析す①では OhHelp の開発活用に向けた計算科学

るなどのポストプロセスも困難である。その分野と計算機科学分野の連携、②では何をど

解決方法の一つとして計算実行とともに可視のように可視化するかを議論するための計算化作業を行う In-situ 可視化の研究・開発を科学分野と可視化情報学分野の連携、また、行う。

In-situ 可視化ライブラリの最適化のための③巨視的シミュレーションに活用するための、

計算機科学分野との連携、③では粒子コード

大規模粒子シミュレーションで得られた知見による研究者と流体コードによる研究者の連に基づく微視的現象のモデル化の研究：たと携である。

え超並列化した粒子コードによっても装置全

このように幅広い研究分野の研究者を含む体などのグローバルな系を扱うことは現在・研究体制を構築するため当該公募型共同研究

近未来のスーパーコンピュータでも不可能では不可欠である。この課題は超大規模計算機ある。粒子コードで得られた知見から微視的と超大容量ストレージを有する情報基盤を必現象を表すようにモデル化(パラメータ化)し要としたグランドチャレンジ的な問題でもあて MHD モデルのような巨視的モデルに組み込る。以上の理由から、多数の先導的研究者がむことで、経験的モデルより予測性の高い巨在籍する JHPCN 公募型共同研究で実施する必視的シミュレーションの実現を目指す。

要がある。
2

学際大規模情報基盤共同利用・共同研究拠点 平成 27 年度共同研究 最終報告書 2016 年 5 月

4. 前年度までに得られた研究成果の概要

コネクション発生領域で発生して排出される

新規課題なので、該当しない。

という現象が地球磁気圏など様々な系で観測

5. 今年度の研究成果の詳細されている。２次元 PASMO コードを使って、本研究課題で超並列化のための開発を行って駆動型磁気リコネクションにおけるプラズモ

いる粒子コードは、電磁粒子コードの PASMO

イドの役割について調べた。

と静電粒子コードの p3bd である。この超並列化のための開発(３次元コードの PASMO と
p3bd、MHD コード)は次節で報告する。ここでは PASMO コード(２次元コード)と p3bd コード
によるプロダクトランの成果についてまとめる。続いて、波動―粒子相互作用に基づく異常抵抗モデル（粒子シミュレーションの結果に基づく微視的モデル）をグローバル MHD シ
ミュレーションに組み込んだ場合の異常抵抗モデルの妥当性を調べた結果について報告する。また、粒子コードの計算結果をモデル化するだけでなく、粒子コードと巨視的コード

図 1：電磁場エネルギーからイオンと電子へを直接接続する多階層モデルの開発も進めてのエネルギー遷移の時間変化。multi-islands

いる。その結果についても報告する。最後にと no-island はそれぞれ流入幅が広がる速さ

In-Situ 可視化ライブラリ VISMO の研究・開が速い場合と遅い場合に対応する。

発について報告する。
初期条件は一様な縦磁場を加えた温度一定○PASMO コードによる成果
PASMO コード(PArticle Simulation codeの 1 次元平衡解で与え、上流境界から粒子をfor Magnetic reconnection in an Open

シミュレーション領域内へと流入させるため

system)は、磁気リコネクション研究のためのに上流境界に駆動電場を加えた。プラズマは、開発された、陽解法による 3 次元電磁粒子シ

最初に上流境界の中心から流入し始めて、そミ ュ レ ー シ ョ ン コ ー ド (H.Ohtani &の後、その流入する幅を広げ最終的には上流

R.Horiuchi : PFR2009)であり、上流境界から境界全体から一様に流入する。この流入するプラズマと磁場フラックスを流入させること幅が広がる速度を変えることでリコネクショ

で、駆動型磁気リコネクションを研究するこンのダイナミクスが大きく変化することが分とができる。下流領域が開放境界になっていかった。幅が広がる速度が遅い場合、プラズ

るので粒子が自由に出入りができ、開放系にモイドが発生しないが、速い場合、プラズモ

おける長時間のダイナミクスを解析すること

イドがリコネクション発生領域で繰り返し発が可能である。

生して下流方向へ排出されることが分かった。

PASMO コードには３次元と２次元のコード

このときの電磁場からイオンや電子へのエネがある。３次元コードは６節で述べるように

ルギー遷移を調べてみると(図 1)、流入する

FX100 への移植等を進めている。２次元コー幅が広がる速さが速い場合(multi-islands)

ドについてはいくつかのプロダクトランが行

は、遅い場合(no-island)に比べ、エネルギー

われ、成果を挙げた。

遷移がより強く行われていることが分かった。
また、別の詳細な解析で、multi-island の場

プラズモイドと呼ばれるプラズマの塊がリ
3

学際大規模情報基盤共同利用・共同研究拠点 平成 27 年度共同研究 最終報告書 2016 年 5 月

合と no-island の場合のリコネクション率を比較すると、multi-islands の場合の方が、より大きなリコネクション率となることが分かった。以上のように、プラズモイドの発生がエネルギー遷移やリコネクション率と大きな相関があることが分かった。
○静電コード p3bd コードによる成果
p3bd コードは、プラズマブロブ(磁力線に沿ったフィラメント状のプラズマコヒーレン
ト構造)による非拡散的動径方向輸送を解明するために開発された 3 次元静電粒子シミュ

図３：ブロブ内部における電子の磁力線方向

レ ー シ ョ ン コ ー ド で あ り (S.Ishiguro &速度分布。電位が高い領域では電子温度が高H.Hasegawa:JPP2006)、このコードを用いて周くなっている。

辺領域のプラズマ輸送現象における運動論的○微視的現象のモデル化と巨視的モデルへのダイナミクスを解析している。磁力線方向の組み込み

両端は固体壁（磁場閉じ込め核融合装置のダ

磁気圏サブストームにある多階層性をグロ

イバータ板に相当）を再現して、荷電粒子がーバル MHD の中で無矛盾に扱う研究を行って吸収される境界条件を設定している。磁場閉いるが、今年度は波動―粒子相互作用に基づじ込め核融合装置周辺領域のスクレイプ・オ

く異常抵抗モデルの評価を行うため、サブス

フ・レイヤープラズマを模擬することができ

トームのグローバル MHD シミュレーションをる。このコードを使って、ブロブ内部に自己行った。MHD に適用する際に、磁気リコネク

無撞着に形成される電流系（図２）、電位構造

ションの周辺磁場がダイナミカルに変化するが確認され、さらに、電位構造に起因する温ことを考慮して、異常電気抵抗モデルに含ま度構造などが明らかになった（図３）

れプラズマ密度と磁場は MHD に従って変化す(Hasegawa and Ishiguro:PoP2015)。

るとした。

図４：サブストームオンセット後の磁場の双極子化の様子。横軸は時間、縦軸は磁気圏磁場の南北成分を表す。左右の図では異常抵抗モデルの大きさが異なるが、グローバル磁場
図２：ブロブ内部に形成された電流系（左手の変化はおおよそ違いがない。

前側の流線）。右奥の等値面は電子密度分布を示している。

グローバルな磁気圏のふるまいとして、サブストームオンセット後の磁場の双極子化の4

学際大規模情報基盤共同利用・共同研究拠点 平成 27 年度共同研究 最終報告書 2016 年 5 月

シミュレーション結果を示したのが図４であＸ点を持ち間欠的なリコネクションが継続する。図からもわかるとおり、異常抵抗モデル

る場合である。

の大きさにほぼ依存しないことが分かった。
これはグローバル MHD シミュレーションでは○In-situ 可視化ライブラリ VISMO の開発

解像度が不十分であり、波動―粒子相互作用

In-situ 可視化ライブラリ VISMO の可視化に関しては、サブストームにおける異常抵抗

手法はソフトウェアレンダリングであり、開モデルの妥当性の検討が困難であることを意

発言語はシミュレーション研究者になじみの味している。このことから、グローバルでは深い Fortran90 を採用している。可視化手法

なく局所モデルを用いた粒子モデルと MHD のとして、粒子・等値面・スライス・流線・矢直接の比較が必須であることが示された。

印を実装している。In-situ 可視化の難点は、○PASMO コードの応用（PASMO と MHD コードの視点を変えるなど異なる可視化を行おうとす連結）

ると、シミュレーションを再度実行する必要

PASMO コードは MHD シミュレーションコーがあることである。その難点を克服する一つ

ドと結合させて、多階層シミュレーションモの手段として、点群データを活用して対話性

デルの開発にも活用されている。これは、磁を持たせるデータ可視化手法を開発して、

気リコネクションの階層横断現象解明のため

VISMO に組み込んだ。図５はテストデータにのモデルである。シミュレーションの実空間

ついてその場可視化を行い、専用ビューワでを分割し、マクロ領域は MHD、ミクロ領域は表示した例である。

粒子手法(PASMO)により計算するモデルである。このモデルを使って、マクロ階層のダイ
ナミクスが磁気リコネクションのミクロな物
理に及ぼす影響について、世界で初めての解析を行った。その結果、マクロ階層のプラズ
マフローのパターンによって、図５(a)のよう
図５：VISMO によるテストデータの可視化。

に単一のＸ点を持ち定常なリコネクションに緩和する場合、(b)のように複数のＸ点を持ち

6. 今年度の進捗状況と今後の展望

間欠的なリコネクションが継続する場合に分○PASMO コードについて類できることが分かった。さらに詳しく調べたところ、このリコネクション物理の違いに

３次元 PASMO コードはこれまで日立製作所

関して定量的に体系づけることに成功した。

製 SR16000 で稼働していたが、今年度は富士
通製 FX100 に移植するための最適化を進めた。
コンパイラによる自動並列ができなかった箇所について OpenMP でスレッド並列を行うことで、最適化前と比べ７倍の性能向上を実現した。しかし、最もコストの高い箇所(粒子運
動方程式を解き、電流密度を計算するルーチン)の構造が巨大すぎるため、SIMD 化やソフ

図５：多階層モデルによる磁気リコネクショ

トウェアパイプラインによる最適化がほとん

ンの解析結果。(a)は単一のＸ点を持ち定常などできておらず、今後の課題である。また、

リコネクションに緩和する場合、(b)は複数のOhHelp を 導 入 する た めの 検 討 を 進 め た 。
5

学際大規模情報基盤共同利用・共同研究拠点 平成 27 年度共同研究 最終報告書 2016 年 5 月

OhHelp ライブラリ(H.Nakashima,et al: Proc.

適化を進める。

Intl. Conf. Supercomputing, pp.90-99, June

○In-situ 可視化ライブラリ VISMO について

2009)は、粒子プラズマシミュレーション向け

In-situ 可視化ライブラリ VISMO では点群の動的負荷分散アルゴリズムを適用したライ

データによる疑似的な対話性を持たせた可視

ブラリで、非一様な粒子密度分布の核融合プ化の実装ができた。今後は、PASMO コードにラズマ粒子シミュレーションにおいて、ロー実装して FX100 での実行を計画している。

ドバランスが悪くなる単純一様領域分割コー○MHD と PIC の比較について

ドに対して、良好な負荷分散が期待される。

磁気リコネクションモデル構築のために、○p3bd コードについて

AMR-MHD コードグローバル MHD コードに電気

FX100 への移行作業などを中心に進めた。

抵抗モデルを導入する作業を進めた。また、

これまで、日立製作所製 SR16000 での計算実

AMR の粗い方のグリッドとグローバル MHD と績があったが、同規模の計算で比較した場合、の接続に着手した。

SR16000 及び富士通製 FX10 より 1.5 倍の速度

○自己評価向上が見られた。SIMD 化なども試みているが、FX100 への最適化では、粒子コード・流体

それ以上の速度向上が見られないため、アルコードともに高速化に成功したことから良好

ゴリズムの再検討を含めて、更なる高速化のな成果を挙げることができた。更なる大規模作業を進めた。一方、本コードでは、これま計算機に向けた最適化には、OhHelp の導入もで、初期の粒子分布として、ブロブ構造しか含めアルゴリズム自体の再検討も必要である

与えることができなかったが、コードの改良

ことがわかり、一定の成果を挙げることがでをおこない、ホール構造（背景プラズマよりきたといえる。また、In-Situ 可視化の研究

密度の低い、磁力線に沿ったフィラメント状では、In-Situ 可視化の最大の弱点である対構造）も扱えるようになった。これにより、

話性の欠如に対して疑似的な対話性を導入す

非拡散的動径方向輸送に関する研究の幅が、ることができ、非常に大きなインパクトを与非常に大きく拡がった。

える成果を得ることができた。巨視的コード○MHD コードについてのための微視的モデルの研究では、波動―粒

我々がこれまで使用してきた、核融合用非子相互作用に基づく異常抵抗モデルの評価は直交曲線座標系表示に基づく電磁流体力学

グローバル MHD シミュレーションでは困難で(MHD) シ ミ ュ レ ー シ ョ ン コ ー ド 、 MHD Inあることがわかり、今後の研究の方向性を得

Non-Orthogonal System (MINOS) コ ー ド のることができて一定の成果を挙げることがで

FX100 への最適化および試験計算を行った。

きたといえる。

このシミュレーションコードは、スイッチによる切り替えにより、物理モデルを圧縮性

7.

MHD 方 程 式 、 拡 張 MHD 方 程 式 、 圧 縮 性(1) 学術論文

Navier-Stokes 方程式など複数の物理モデル

・ H.Ohtani,の中から選択可能などの特徴をもつ。現段階

S.Ishiguro,

で は 、 圧 縮 性 MHD 方 程 式 、 圧 縮 性“Visualization of dust particle data with

Navier-Stokes 方程式を直角格子で解く場合

plasmaについて最適化および試験計算が進行し、最virtual-reality

適化前に比べて 40%程度の高速化に成功した。

Contributions to Plasma Physics (2015).

今後は、非直交格子で解く場合について、最・H.Miura, K.Araki, F.Hamba: “Hall effects
6

研究成果リスト
M.Shoji,N.Ohno,A.Kageyama,
simulation
system”,

Y.Suzuki,
Y.Tamura:

resultsbysubmittedto学際大規模情報基盤共同利用・共同研究拠点 平成 27 年度共同研究 最終報告書 2016 年 5 月

andsub-grid-scalemodeling

inNumerical Simulation of Plasmas (ICNSP2015),

turbulence

August 12-14, 2015, Golden, Colorado, USA.

simulations”, to appear in Journal of・S. Usami, R. Horiuchi, H. Ohtani, M. Den:

Computational Physics (2016).

“Development

・ H.Hasegawa and S.Ishiguro, "Study ofSimulation

Model

Based

on

MHD-PIC

Self-Consistent Particle Flows in a Plasmainterlocking

andStudieson

Magnetic

Blob with Particle-in-Cell Simulations",

Reconnection”, JIFT Workshop on "Extended

Phys. Plasmas, 22, pp.102113 (1-8) (2015).

MHD and MHD simulations for magnetizedmagnetohydrodynamicofaMulti-Hierarchy(2) 国際会議プロシーディングス

plasmas”, August 10-11, 2015, Golden,(3) 国際会議発表

Colorado, USA.

・ H.Ohtani,

M.Shoji,Y.Suzuki,・H.Miura, K. Araki, A.Ito: “Short-wave Hall

Y.Tamura:and gyro-viscous effects and their modelling“Visualization of dust particle data within extended magnetohydrodynamic simulations

plasmaof

S.Ishiguro,

N.Ohno,A.Kageyama,
simulation

virtual-reality

results

system”,

byInternational

instability

International

andturbulence”,

Conference

on

24th

Numerical

Workshop on Plasma Edge Theory in Fusion

Simulation of Plasmas (ICNSP2015), August

Devices, September 9-11, 2015, Nara, Japan.

12-14, 2015, Golden, Colorado, USA.

・R.Horiuchi,H.Ohtani, S.Usami:“Macro-Micro
Interlocked
PlasmaSimulation

Phenomena

Reconnection”,

for

M.Nunami,“Hierarchical

H.Usui,

Multiscale

like

Magnetic

Rayleigh-Taylor instability on extended MHD

Internationalwith AMR framework and visualization”, 24thSpace

H.Miura:A.M.Ito,

for

12th

School/Symposium

・ T.Hatori,International

Simulations

simulation

Conferenceonof

Numerical(ISSS-12), July 3-10, 2015, Prague, Czech.

Simulation of Plasmas (ICNSP2015), August・ R.Horiuchi: “Basis of Particle-In-Cell

12-14, 2015, Golden, Colorado, USA.

(PIC) simulation and its applications”, 5th・H. Miura, T. Hatori, and K. Araki: “Scale

East-Asia School and Workshop on Laboratory,

Hierarchy

Space, and Astrophysical plasmas, August

Turbulence”, The 15th European Turbulence

17-22, 2015, POSTECH, Korea.

Conference, August 25-28, 2015 , Delft,

・R. Horiuchi, S. Usami, H. Ohtani: “Role of

Netherlands.

magneticislandsin・ H. Hasegawa and S. Ishiguro: “Study ofprocessofcollisionlessreconnection”,
Conferenceonenergy24th
Numerical

conversion

Plasmadrivenin

Homogeneous

Coherent

Structure

Hall

MHD

Transport

International

Dynamics with Particle-in-Cell Simulation”,

Simulation

25th

of

International

Toki

Conference,

Plasmas (ICNSP2015), August 12-14, 2015,

November 3-6, 2015, Toki, Gifu, Japan.

Golden, Colorado, USA.

・ H. Hasegawa and S. Ishiguro: “ Kinetic

・ S.Usami, H.Miura, R.Horiuchi, H.Ohtani,

Effects on Plasma Blob Dynamics with Plasma

M.Den: “Improvement of the Multi-Hierarchy

Sheath” , 57th Annual Meeting of the APS

Model for Magnetic Reconnection Studies –Division of Plasma Physics, November 16-20,

Interlocking between PIC and Extended MHD

2015, Savannah, Georgia, USA.

-”,・ M.Den, R.Horiuchi, S.Usami, Y.Kubota,

24thInternational

Conferenceon7学際大規模情報基盤共同利用・共同研究拠点 平成 27 年度共同研究 最終報告書 2016 年 5 月

T.Tanaka, and H.Ohtani: "Comparative Study

ミュレーションモデルのさらなる発展」プラズ

of Effective Resistivity Models Evaluated in

マシミュレータシンポジウム 2015、土岐。

Kinetic Studies Using MHD Simulation",・宇佐見俊介、大谷寛明、堀内利得、田光江：
「実

Fall

空間分割型の多階層シミュレーションによる

Meeting, December 14-18, 2015, Sanfransisco,、日本物理
駆動型磁気リコネクション研究 II」

USA.

学会 2015 年秋季大会、吹田。

・ T.Ogawa, S.Usami, R.Horiuchi, M.Den,and「プラズマコヒーレン・長谷川裕記, 石黒静児：

K.Yamashita:"Multi-scaleof

ト構造の輸送ダイナミクスの 粒子シミュレー

using

ション」プラズマシミュレータシンポジウム

AmericanGeophysicalmagnetic

Union

(AGU)

simulation

reconnection

particle-in-cell with magnetohydrodynamics

2015、土岐。

on adaptively refined mesh" P2-59, The 25th・長谷川裕記、石黒静児：「プラズマコヒーレン

International Toki Conference, Toki, Gifu.

ト構造の輸送ダイナミクスの粒子シミュレー

・ T.Ogawa,ション」、プラズマ・核融合学会第 32 回年会、

S.Usami,R.Horiuchi,

M.Den,

K.Yamashita, and T.Tanaka: "Simulation of

名古屋市千種区。

Magnetic Reconnection in Magnetotail by・長谷川裕記：「連結階層モデルを用いた大規模Interlockingandシミ ュレーシ ョン」、分子科 学研究所 Open

Magnetohydrodynamics on Hierarchical Mesh",seminar「局所と大域：データサイエンスと階SM51B-2565, AGU 2015 Fall Meeting, December

層性の創発」、愛知県岡崎市。

Particle-in-Cell「第・長谷川裕記、KINDSP プロジェクトグループ：

14-18, Sanfransisco, USA, 2015.
・S.Usami, H.Miura, R.Horiuchi, H.Ohtani, and一原理的手法による周辺プラズマダイナミク

M.Den: "Improvement of a Multi-Hierarchy

スの研究」、 第 5 回 IFERC-CSC 研究会、東京都

Simulation Model for Magnetic Reconnection

港区。

Studies",Toki・長谷川裕記、石黒静児：「ホール構造伝播の粒

Conference, Ceratopia Toki, Toki, Japan,

子シミュレーション」、日本物理学会第 71 回年

November 3-6, 2015.

次大会、仙台市泉区。

・ S.Usami,

25thInternational

R.Horiuchi,

andH.Ohtani:

・S.Usami, R.Horiuchi, and H.Ohtani："Studies"Simulation studies on effective ion heatingonin the reconnection downstream", US-Japan

Downstream of Magnetic Reconnection", 「自Workshop on Magnetic Reconnection, Andaz

然科学における階層と全体」シンポジウム, 国Napa Hyatt, Napa, CA, USA, March 7-11, 2016.

立天文台（東京）。

・ N.Ohno and H.Ohtani, “Development of・宇佐見俊介、堀内利得、大谷寛明：リコネク

In-Situ Visualization tool VISMO”, The 6thション下流における実効的イオン加熱のシミ

AICSュレーション研究、磁気リコネクション研究のInternational

Symposium,

February

22-23, 2016, (Kobe, Japan)

the

Ion

Heating

Mechanisminthe

最前線と今後の展望、国立天文台（東京）。

(4) 国内会議発表(5) その他（特許，プレス発表，著書等）

・大谷寛明、堀内利得、宇佐見俊介：「開放系電・核融合科学研究所プレスリリース「核融合研

磁粒子シミュレーションコード PASMO の開発」、究が更に進展プラズマシミュレータシンポジウム 2015、土岐。

融合し、バーチャルリアリティで表示－」(平・宇佐見俊介、三浦英昭、堀内利得、大谷寛明、―実験データと計算データを

成 28 年 4 月 5 日)、

田光江：「拡張 MHD と PIC 連結による多階層シ

http://www.nifs.ac.jp/press/160405.html
8