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

jh150004-NA03

超並列宇宙プラズマシミュレーションの研究
三宅洋平（神戸大学・課題代表者）、中島浩（京都大学）、加藤雄人（東北大学）、
小路真史（名古屋大学）、臼井英之（神戸大学）、大村善治（京都大学）
岩下武史（北海道大学）、深沢圭一郎（九州大学）、木倉佳祐（神戸大学）
概要 本研究の目的は、Particle-in-Cell 方式に基づく宇宙プラズマ粒子シミュレーションを大規模並列計算機システムへ効率的に実装し、宇宙プラズマ理工学の重要問題に適用することである。今年度は、今後さらなる普及が予想されるメニーコアプロセッサの効率的利用のため、プラズマ粒子カーネル計算部分の SIMD 最適化手法の検討を進める。
これに加え、OhHelp 実装済の全粒子プラズマシミュレーションにより、「月」に見られる日照－日陰境界付近のプラズマダイナミクスと表面帯電特性を解明する。また地球放射線帯の相対論的高エネルギー電子の加速機構、そして非線形イオンプラズマ波動励起過程に関する粒子シミュレーションを実施し、宇宙プラズマ理工学分野における先駆的な研究成果を創出する。

1.

共同研究に関する情報

大な数のプラズマ代表粒子と空間格子点上の電磁(1) 共同研究を実施した拠点名

場の相互作用計算がシミュレーションの大部分を本課題は、神戸大学・京都大学・東北大学・名占める。個々のプラズマ粒子はシミュレーション

古屋大学・北海道大学・九州大学の複数の研究者空間の任意の位置に存在しかつ時間的に移動する。

が参加する共同研究であり、京都大学の大規模計

このため、旧来の均等型領域分割方式による MPI算資源を利用する。

並列では、プロセッサ間の負荷バランスが崩れる。

(2) 共同研究分野：超大規模数値計算系応用分野

そこで計算機科学分野と応用分野双方の研究者に(3) 参加研究者の役割分担

より共同開発した、動的負荷分散技法 OhHelp の参加研究者の役割分担は前年度の体制に、大学適用により、この問題を解決することに成功した。

院生の木倉氏を加えた構成である。詳細は H26 年②は、PIC カーネル計算の共有メモリ並列化時度最終報告書、表 1 を参照願いたい。宇宙プラズにメモリアクセス競合を引き起こしうる。したがマ理工学分野の中でも粒子モデル (PIC) シミュってこれを回避するために適切なスケジューリン

レーションを専門に扱う研究者と流体計算の専門グが必要となる。また、アクセスパターンが粒子家、そして計算機科学分野研究者・技術者が有機位置に依存して決まる、いわゆる間接アクセスを的に協力して行う共同研究体制となっている。

有したままだと、コンパイル時に SIMD 演算命令が発行されない。このため、本報告書で検討する

2.

研究の目的と意義

ように、コード最内ループのプログラム構造やデ本研究の目的は、Particle-in-Cell (PIC) 方式にータ構造を大幅に改変することが必要となる。

基づく宇宙プラズマシミュレーションを対象に、メニーコアプロセッサを基盤とする将来の超並次世代の大規模並列計算システムへの高効率実装列システムにプラズマ粒子計算を高効率に実装すを実現することにある。この目標に向けて解決するため、上記 2 点を同時に解決するコードを開発べき点は、①不均一な粒子分布に起因する負荷不することが、当課題の大きな目標である。

均衡、および②粒子位置座標に依存して決定され実問題への応用に関しては、OhHelp 適用済のる格子点配列へのランダムアクセス、である。

全粒子シミュレーションコード EMSES により、①について、PIC 粒子シミュレーションは、膨月などの天体の固体表面と宇宙プラズマの相互作
1

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

用現象を再現し、将来の着陸探査計画に貢献可能

流体ハイブリッドコードへの OhHelp 実装が完了な定量データを取得する。またプラズマ波動の 1し、その有効性が確認された。

種であるホイッスラーモード・コーラス放射と光応用分野における成果としては、衛星プラズマ速の 99 %を超える相対論的高エネルギー電子と環境シミュレータ EMSES を用いた大規模 3 次元の相互作用に関する大規模計算機実験を実施し、シミュレーション解析により、地球磁気圏の希薄高エネルギー粒子生成過程を明らかにする。

なプラズマ流中におかれた衛星周辺に形成される上述の課題を遂行することは OhHelp 手法の有大規模なウェイク構造を再現した。また当該コー効性の実証という数値アルゴリズム的観点、およドを用いた大規模シミュレーション解析は、太陽び超大規模モデルによる宇宙プラズマ粒子シミュ

近傍環境探査衛星と太陽コロナプラズマ間の電磁レーションの実現と最先端研究成果の創出という的相互作用の研究にも適用された。

宇宙科学的観点の双方において大きな意義がある。

理学的応用としては、電磁イオンサイクロトロ
ントリガード放射のハイブリッドシミュレーショ

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

ンにより、宇宙プラズマ中のプロトン分布関数が本研究は、アルゴリズム開発、コード開発、大大きく歪められることが発見された。また、放射規模シミュレーション実行の全ての段階で、複数線帯での相対論的電子生成機構において重要な役の大学の計算科学者と同センターの計算機科学者割を果たすと考えられているプラズマ波動・コーとの緊密な連携によって実施されている。このよラス放射の発生過程を再現する計算機実験を実施

うな研究体制の重要性を端的に表す象徴とも言うし、その波動特性が背景磁場強度の空間構造の変べき中核技術が動的負荷分散技法 OhHelp である。 化に応じて大きく変化することを明らかにした。
OhHelp は、粒子と格子という 2 つの異なる並列計算の効率性を同時に達成する一般的手法の一つ

5. 今年度の研究成果の詳細として捉えることができ、またそのアルゴリズム(1) 月縦孔周辺プラズマ環境の大規模全粒子シミは、木構造を用いた負荷均衡状況の再帰的なチェュレーション

ックなど、計算機科学的かつ非数値的な発想に基「かぐや」衛星による地形カメラ観測により、づいている。したがって、非数値的な並列計算ア月面に特徴的な縦孔構造が発見されている。このルゴリズムによる並列数値シミュレーションへの縦孔や地下空洞を直接探査する UZUME 計画が貢献を表す好例となっている。

立案されている中、将来の着陸探査への影響を知る上でも縦孔周辺プラズマおよび帯電ダスト環境

4. 前年度までに得られた研究成果の概要を事前評価することは重要である。そこで

領域分割方式の超並列計算が困難とされてきたOhHelp 適用済の全粒子プラズマシミュレーショ

プラズマ粒子 PIC シミュレーションに関して、プンコード EMSES を、月面昼側の月縦孔周辺プラロセス間の負荷不均衡を回避できる動的負荷分散ズマ環境の定量評価に適用した。

技法 OhHelp を宇宙プラズマシミュレーションに

本解析では、Marius の丘に発見されている直径適用した結果、優れた並列計算効率が得られるこ50 m、深さ 45 m の月縦孔に着目し、その上空もとが実証された。また PIC のカーネル計算部分に含めて、おおよそ 200 m×200 m×1 km の空間をついて、マルチカラー法を用いたマルチスレッド計算領域とし、京都大学 CRAY XC30 の最大 128計算技法を提案し、Xeon Phi 上で 102 並列度までノードを利用した計算により、①縦孔近傍のプラ良好なスケーラビリティが維持されることが確認

ズマ分布および電気環境、②その中での帯電ダスされた。これに加えて高エネルギー電子を粒子、トの挙動、を明らかにした。

背景電子を流体としてモデル化する独自の粒子・

一般的にプラズマに固体表面が接している状態
2

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

図 1:月縦孔近傍の電位分布および帯電ダスト挙動の大規模シミュレーション解析．
では、プラズマ荷電粒子が表面に捕捉されること積しているという仮定の下での計算であるが、静により、帯電が生じる。図 1a に太陽風（太陽光）電気力の作用により、ダストが浮遊していること入射角度=30のケースにおける月縦孔を含めたが図から確認できる。ダストの持つ電荷は日向で表面電位分布を示す。図から明らかなように、縦は正、日陰では負であるが、一方で日向および日孔内部では、場所によって帯電値が大きく異なり、陰の月面からはそれぞれ上向きおよび下向きの電大きな電位差（40 V 程度）が生じている。これは、場が発生しているため、ダストが浮かび上がる向

日照面では先述した光電子放出により、正に帯電きに静電気力が作用する。この静電気力と下向きするのに対し、日陰面からは光電子は放出されず、の重力との釣り合いにより、ダストが上下に振動周辺を浮遊する電子の一部が表面に付着することする運動が確認されている。現在はダストの直径で負に帯電するためである。この電位差により、と質量をある一定値に固定した計算であるが、今日照－日陰境界付近には大きな電場が生じている後はダストの直径に分布があることを想定した計ことが示唆され、着陸探査においては留意すべき算を実施し、平均的なダスト密度分布を算出する点であるといえる。この結果は学術論文として、ことを目指す。

Icarus にて発表済である。
帯電は月面のみならず、月表層に堆積している(2) 月磁気異常‐太陽風プラズマ相互作用の大規微小な塵であるダストに対しても発生する。帯電模全粒子シミュレーション

したダストは、月縦孔周辺の特異な電気環境のも上述の解析は、月面上空およそ 1 km までの描とで浮遊し、特徴的なダイナミクスを示すことが像であるが、より上空では太陽風プラズマと月に予想されるが、こうした帯電ダストが月面探査用局所的に存在する磁気異常の間の電磁的相互作用のローバーに付着することにより、機器の故障なが電気環境の決定に大きな役割を果たす。そこで、どの障害を引き起こす恐れがある。そこでダスト

Reiner Gamma と呼ばれる月磁気異常をモデルケの挙動を明らかにすることを目的に、テスト粒子ースとして採用し、プラズマ粒子シミュレーショシミュレーション解析を実施した。

ンによる解析を実施した。ダイポール中心から磁図 1b に縦孔周辺でのダスト分布のスナップシ

気圧と太陽風動圧が釣り合う点までの距離を磁気ョットを示す。初期状態ではダストは月面上に堆異常の代表長 L とすると、磁気異常ではその代表
3

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

図 2:月磁気異常と太陽風プラズマの相互作用によって形成される大規模電流構造．
長 L が太陽風電子のジャイロ半径よりも十分大き(3) 惑星磁気圏における波動粒子相互作用に関す

く、イオンのジャイロ半径より小さい。このようる大規模ハイブリッドシミュレーションな条件下では、太陽風プラズマと磁気異常磁場間放射線帯は光速の 99 %を超える相対論的高エの相互作用により小型磁気圏が形成される。またネルギー粒子が地球を取り巻くように分布してい磁気圏境界層における太陽風電子とイオンの磁場る領域である。放射線帯の相対論的粒子のフラッに対する応答差に起因する静電界が磁気圏形成にクス量が、いつ・どこで・どのように変動するか

重要な役割を果たす。

を明らかにすることは、宇宙天気分野での最重要

本シミュレーションでは、特に磁気圏境界層に課題の一つである。本研究の目的は OhHelp 適用

おける電子ダイナミクスによって決定づけられる済の電子ハイブリッドコードを用いた大規模数値電流構造に着目して、詳細解析を実施した。シミ解析により、放射線帯電子を作り出す過程で重要ュレーション結果から、境界層低緯度領域においな役割を果たすとされるプラズマ波動、ホイッスては朝側から夕方側に向かう電子ドリフト運動がラーモード・コーラス放射の励起過程を再現し、顕著であるが、中高緯度領域においては逆に夕方相対論的高エネルギー電子が生み出される物理過側から朝側に向かう電子の流れが主に見えた。こ程を定量的に理解することである。のような二層の電子ドリフト運動の構造により磁今年度は、発生するコーラス放射の波動特性が、気圏層間側、すなわち磁気異常上空の南北両半球背景磁場の空間スケールならびに高エネルギー電において電子電流の渦構造が見られた（図 2）。特子の初期分布に対してどのように変化するかを究に赤道面における電子電流には、電子－イオン間明している。放射線帯は固有磁場を持つ惑星の磁の荷電分離によって発生する電場と磁気異常磁場気圏に普遍的に存在するが、コーラス放射も同様による E×B ドリフトと、磁気圏境界層付近でのに惑星磁気圏に共通したプラズマ波動である。本

電子圧力の急速な減少に伴う反磁性ドリフト電流研究の成果から、地球・木星・水星磁気圏でのコが大きな役割を果たしていることが明らかとなっーラス放射発生条件の考察が可能となるデータセた。これらの結果は小型磁気圏の構造、ひいてはットが得られており、学術論文への投稿を準備中その周辺でのイオンのダイナミクスにも大きく影である。

響を及ぼすという点で、月プラズマ科学上、きわ6.

めて重要な成果ということができる。

今年度の進捗状況と今後の展望(1) 当初の目標に対する達成度の自己評価：90%4

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

スする必要がある。ここで vp とt はそれぞれ粒子(2) 進捗の詳細 1：PIC カーネル計算の SIMD 最適

速度と時間ステップ幅を表し、(xpvpt)により粒実装の検討
Xeon Phi における SIMD 演算機構は、NEC SX子移動前の位置を示している。ここでコンパイラシリーズに代表されるようなベクトル計算機におは関数によって計算されるインデックスに規則けるパイプライン処理に概念上は類似しているが、 性を見つけることができないため、8 つの粒子そ
Xeon Phi における SIMD 演算器はリストベクト

れぞれに対する E((xp))、B((xp))、J((xpvpt))、ルのような強力な機構を持たないことや、ベクトそして J((xp))から構成される SIMD ベクトルルレジスタの数が少ないなどの理由により、複雑48+24=72 本をレジスタに登録することで SIMDなデータアクセスパターンを有するアプリケーシ命令発行を試みる。しかし、現実には Xeon Phi

ョンへの適用性は必ずしも高くない。また SIMDで登録可能な最大ベクトル数は 32 であるため最命令発行の可否は最終的にはコンパイラの判断に終的には SIMD 命令発行が不可であると判断され委ねるしかないため、カーネル計算部のプログラてしまうのである。
この問題を解決するためには、次の 2 つの対策

ム構造やデータ構造を変更しながら、SIMD 命令が発行可能な実装を見つけ出す、といった試行錯をとればよいと考えられる。

誤を繰り返すこととなる。

A) 参照書き込みを行う格子点を 8 つの粒子の間

このレベルの最適化においては、PIC 計算のカで共通化する．

ーネル部分である粒子速度位置更新、電流計算をB) A の事実が「コンパイラから見て」明らかに含む粒子ループの SIMD-aware 実装が焦点とななるよう、電磁場や電流の参照書き込みにスる。特に電流計算では、各粒子がつくる電流要素カラー変数を用いる．

を、それらの位置座標に隣接する格子点上の電流

A の対策は、粒子の所属セルをキーとしたソーテ格子点配列に配分する。この粒子位置は実行時まィングで実現される。ここで、粒子ループ中ではで未知数であるため、格子点配列へのアクセスは粒子位置更新が行われるため、その際に所属セルしばしばランダムになる。

が変化する可能性があり、その場合は A が満たさ

Xeon Phi は 512 bit 幅の SIMD ベクトルを持つれなくなる点に注意する必要がある。そこで従来ため、原理上は 8 つの粒子を同時に処理することの長大な粒子ループを、速度更新、位置更新前のが当面の目標であるが、従来のプログラムにおい電流計算、位置更新後の電流計算という 3 つのルては、ループ内のプログラム構造が複雑すぎるとープに分割することで、各ループの中で常に A をして、Intel コンパイルによる SIMD 命令発行が満たすことが保証される。

なされていなかった。この主な要因は、電磁場やB については、これまで 1 重の粒子ループであ電流配列に対する間接的なアクセスにある。粒子ったものを、セルに関する外側ループとそのセルループの各繰り返しにおいて、速度更新の際にはに所属する粒子の内側ループから成る、2 重ルー粒子に隣接する 8 格子点分の電磁場ベクトル（6プ構造に変更する。これによりセルループの内側成分）、すなわち 48 成分の E((xp))や B((xp))を参では、(xp)は共通であるから、E((xp))、B((xp))

照する。ここで xp は浮動小数点であらわされた粒を e や b、J((xp))を j といったスカラー変数にコ子位置、関数はそれを周辺の 8 格子点に整数化すピーすることが可能である。内側での粒子ループる処理を表す。同様に電流計算においては、粒子

では、元の配列ではなく、e、b、そして j を対象が存在するセルの辺上に定義される電流ベクトルに計算を実施する。これらを踏まえ、変更前後の12 成分に参照書き込みが行われるが、それを粒子ループ構造を表したのが、図 3 および 4 である。

移動前後それぞれについて行う必要があるため、このような変更を施すことにより、PIC 計算のカ

計 24 成分の J((xpvpt))および J((xp))にアクセーネル部分である粒子速度位置更新、電流計算に5

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

ついて、コンパイルによる SIMD 命令発行が可能

E[][][]/B[][][]となることを確認した。

for(p in a block)

次に粒子のデータ構造に関して、①各粒子の持48

48

Update vp
Update xp
Scatter J((xp-vp))
Scatter J((xp))

つ座標成分のメモリ上での配置、および②多数の粒子の格納方法、の変更を実施する。まず①につUpdate vp
Update xp
Scatter J((xp-vp))
Scatter J((xp))

48

いて、粒子配列のレイアウトを構造体配列48（ Structure of Array: SoA ） か ら 配 列 構 造 体

J[][][]

図 3:SIMD 最適化前の PIC カーネルループ。

（Array of Structure: AoS）へ変更する。この変換の目的は、複数粒子の同一座標成分がメモリ上for(each cell)で連続して並ぶようにすることにある。これによE[][][]/B[][][]り、図 4 に示した同一セル内粒子のループにおけ48

Scalarized e/b

るベクトルレジスタへの粒子データの load/store

for(part.in cell)

Update vp

に関して、SIMD 命令が発行されることになり、

Update vp

for(part.in cell)
Scatter J((xp))
Update xp

演算効率の向上が見込める。
②に関しては、従来のように複数粒子のデータ

Scatter J((xp))
Update xp
24

Scalarized jを粒子配列上に前から順に詰めて配置するのではなく、図 5 のように同一所属セルの粒子毎に固めJ[][][]

Sort particles
for(each cell)for(part.in cell)て配置することとし、所属セルが異なる粒子群の間には適当な数のマージン（空き領域）を設ける。

Scatter J((xp))

これは次に述べる「On-the-fly 方式の粒子ソーテ

Scatter J((xp))
24Scalarized jィング」の実装のために必要な処置である。ここで所属セルが異なる粒子群の間のマージンの幅は、J[][][]

図 4:SIMD 最適化後の PIC カーネルループ。

各セルに存在する粒子数の増減に従って変化し、より、可能となっている。また一般的に各時間スあるセルに粒子が集中する状況においては、あらテップでセル境界を超える粒子の割合は比較的小

かじめ定められたセル毎の格納可能個数に対してさいと想定されるため、本方針はソーティングコオーバーフローを引き起こすことも想定される。

ストの削減において非常に有効である。

このような場合には、例外処理としてオーバーフ本方式に基づく粒子ソーティングの動作例を図

ローが発生したセル用のバッファ長を増加させる

5 に示す。図中セル c1 に属する 3 番地の粒子が隣処理を行う。この例外処理はメモリ上での大規模のセル c2 に移動したとする。まず移動する粒子のな粒子移動を必要とするため、慎重な実装が必要情報をセル c2 のマージンの先頭 78 番地に移動すとなるが、本稿で詳述することは避ける。

る。これによりセル c2 に所属する粒子について、(3) 進捗の詳細 2：On-the-fly 方式にもとづく低49 から 78 番地まで粒子データが隙間なく格納されることが確認できる。このとき、セル c1 についコスト粒子ソート手法の検討
本方式による粒子ソーティングの要点は、メモては 3 番地に空きが生じた状態となっている。こリ空間上で移動対象となる粒子の個数を、実際にの空き番地を埋めるにあたっては、セル c1 に所属セル間を移動した粒子個数に抑えることで、ソーする粒子群の末尾である 32 番地のデータを当該番地に移動させる。これにより移動元のセルにおティングのコストを削減することにある。これはいても、粒子データが先頭から隙間なく格納され同一セル内の粒子に関して、格納順序を考慮するることを保証する。このような手順の繰り返しに必要が無い、という PIC アルゴリズム上の特徴により粒子ソーティングにかかるオーバーヘッドコ
6

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

図 5:粒子データの格納方法と On-the-fly 方式に基づくソーティング手法．
ストを最小化することが可能である。
前節までに述べた最適化を実装した結果、プラズマ粒子計算のカーネル計算部である粒子速度更新および電流計算の最内ループがベクトル化されること、および計算結果が正しいことを確認した。
以下では本コードを用いた性能評価結果を紹介する。15  120  120 グリッドの計算空間内に、電子とイオンをそれぞれ 2,073,600 個ずつ一様に分布させ、100 時間ステップ分計算を実施した。粒子の初速度は x、y、z のいずれかランダムな成分
図 6:SIMD 最適化後の PIC 計算性能評価．に 0.01 grid/step を与えた。また、粒子ソーティ
ング用に導入した同一セル所属セルの間のマージ

現在取り組んでいる SIMD 最適化の先には、こンの大きさは cell あたりの粒子数（本評価では 96）れまで取り組んできたノード間並列やノード内並の 75%と設定した。

列技法との組み合わせにより、メニーコアアーキ

図 6 に使用コア数を変化させて測定した性能評

テクチャを基盤とした超並列システムに対応した価の結果を示す。60 コアかつハイパースレッド数

プラズマ粒子シミュレーション技法の確立を目指

HT=3 の 180 スレッドで実行した際に 414  106す。これは簡単な作業ではなく、各アーキテクチparticles/sec の最高性能を得た。また台数効果にャ階層の最適化実装の間の両立性を丁寧に検討しついては 1 コア実行時の 53.1 倍の速度向上が得らていかなくてはならない。H28 年度以降、こうしれた。図には示されていないが、SIMD 最適化をた課題にも取り組んでいく。

施す前のプログラムでは、180 スレッド実行での性能は 62.7  106 particles/sec にとどまっており、今回の SIMD 最適化により 6.6 倍の性能向上を達

また高度化された計算コードによる大規模シミ
ュ―レーションを、引き続き宇宙プラズマ理工学の重要問題に適用していく。特に今後は、日本で成できたことになる。本稿で述べたような包括的計画される宇宙探査ミッション（例えば ERG 等）な最適化の必要性を示す結果と言える。と密接に関連した形で、宇宙プラズマ物理過程に関する大規模数値解析を実施する。

(4) 今後の展望
7

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

of Whistler Mode Emissions, 12th ISSS,

7.

Prague, 2015.

研究成果リスト5.

(1) 学術論文
1.

2.

Miyake, Y. and M.N. Nishino, Electrostaticand H. Nakashima, Simulation study ofEnvironment near Lunar Vertical Hole: 3D

dependencies of the generation process of

Plasma Particle Simulations, Icarus, 260,

whistler-mode

301-307, 2015.

emissions on temperature anisotropy of

Miyake, Y., Y. Nishimura, and Y. Kasaba,energetic electrons, 12th AOGS, Singapore,

Asymmetric

2015.

Electrostatic

Environment

around Spacecraft in Weakly Streaming

3.

6.

Miyake, Y., and M.N. Nishino, Full-Particle
Environment near Lunar Vertical Holes,

Nakashima, H., Manycore challenge in

AGU Fall Meeting, San Francisco, 2015.
7.

on

Electrostatic

PlasmaMiyake, Y., and M.N. Nishino, Full-Particle

TFlops peak performance for simulation

Simulations

codes

Environment near Lunar Vertical Holes,

withirregularcomputation,on

12th81-94, 2015.

Conference, Woodlands, 2016.

1.

onIn-SpaceElectric

Nakashima,Dependencies

Processof

2.

of

3.

the加藤雄人, 深沢圭一郎, 木星磁気圏におけるッド・MHD 連成計算, 第 138 回 SGEPSS
総会および講演会, 東京, 2015.

4.

Nakashima, H., K. Kikura, and Y. Miyake,PrototypeMeeting, Tokyo, 2015.

Implementation

andItsFundamental Performance Evaluation of aM.N. Nishino, Fullon

三宅洋平, 西野真木, 月縦孔プラズマ静電環

コーラス放射発生過程に関する電子ハイブリ

Whistler-Mode

Magnetosphere, URSI-Japan Radio Science

Simulations

and

59 回宇宙科学技術連合講演会, 鹿児島, 2015.

Energetic Electrons in the Earth's Inner

Particle

System境に関する大規模粒子シミュレーション, 第Emissions on Temperature Anisotropy ofandComputing

Infrastructure, つくば, 2015.（招待講演）

FieldKato, Y., Y. Omura, Y. Miyake, H. Usui, andMiyake, Y.,

Science三宅洋平, 超並列宇宙プラズマ粒子シミュレ
AdvancedMeeting, Tokyo, 2015.（招待講演）

Generation

Planetaryー シ ョ ン の 研 究 , Annual Meeting on

Measurements, URSI-Japan Radio Science

H.

andPlasma(4) 国内会議発表

Miyake, Y., Spacecraft-Plasma Interaction
Effects

Lunar

Electrostatic

Computers & Electrical Engineering, 46,(3) 国際会議発表

Manycore-Aware

PlasmaOhHelp’ed

PIC

Simulation Code, 第 153 回ハイパフォーマ

Environment around Lunar Vertical Hole,

ンスコンピューティング研究発表会, 松山,

12th ISSS, Prague, 2015.
4.

hiss-like2015.

該当なし

3.

andSimulations(2) 国際会議プロシーディングス2.

chorusPlasmas, J. Geophys. Res., 120, 6357–6370,

particle-in-cell simulation: How to exploit 1

1.

Katoh, Y., Y. Omura, Y. Miyake, H. Usui,2016.

Katoh, Y., Y. Omura, Y. Miyake, H. Usui, H.
Nakashima, Electron-hybrid Simulations of(5) その他（特許，プレス発表，著書等）

the Dependencies of the Generation Process

該当なし8