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

jh140042-NA23

次世代降着円盤シミュレータの開発
松元亮治（千葉大学）
概要 重力を及ぼす天体のまわりに形成される降着円盤の時間発展を計算する３次元磁気流体シミュレータを改訂して、その適用範囲を拡大することを目的として、共同研究を実施した。まず、HLLD 法に基づく高次精度の磁気流体コードを実装したシミュレー
ションソフトウェア CANS+を整備・公開し、シミュレーションスクール等で使用した。
また、一般相対論的な３次元輻射磁気流体コードを作成してブラックホール近傍のシミュレーションを可能にし、高い降着率を持つブラックホール降着流におけるダイナモ計算やジェット形成に適用した。輻射輸送方程式を各光線方向について解く直接解法の実装、シミュレーション結果に基づく輻射スペクトル計算モジュールの作成も進めている。
自己重力計算モジュールを実装した解適合格子法に基づく流体シミュレーションコードSFUMATO の公開を進め、利用者からのフィードバックを加えた改訂を行っている。

1.

共同研究に関する情報

町田真美（九州大学・大学院理学研究院）(1) 共同研究を実施した拠点名

２温度降着円盤の大局的 3 次元磁気流体数

東京大学

値実験の実施、スペクトル計算(2) 共同研究分野

松本洋介（千葉大学・大学院理学研究科）■ 超大規模数値計算系応用分野

シミュレーションコードのチューニング、

超大規模データ処理系応用分野

コード整備・公開

超大容量ネットワーク技術分野



超大規模情報システム関連研究分野

花輪知幸（千葉大学・先進科学センター）
自己重力磁気流体コードの並列性能向上(3) 参加研究者の役割分担

中村賢仁（九州産業大学・工学部）

松元亮治（千葉大学、代表者）、全体統括、陰的差分を用いた熱伝導モジュールの並列

降着円盤シミュレーション実施

性能向上、降着円盤への適用

片桐孝洋（東京大学・情報基盤センター、川島朋尚（国立天文台）副代表）、陰解法に適用できる並列性能の高超臨界降着流の輻射スペクトル計算

い連立１次方程式解法の実装

朝比奈雄太（千葉大学・大学院理学研究科）

廣瀬重信（海洋研究開発機構、副代表）

一般相対論的輻射磁気流体コードの実装、陰的差分を用いた輻射磁気流体コードの並降着円盤、ジェット形成への適用

列性能向上、降着円盤への適用

工藤祐己（千葉大学・大学院理学研究科）

大須賀 健（国立天文台・理論研究部）

陰的差分を用いた熱伝導モジュールの並列輻射輸送方程式解法の実装、降着円盤シミ性能向上、降着円盤への適用ュレーションへの適用

谷田部紘希（千葉大学・大学院理学研究科）

高橋博之（国立天文台・天文シミュレーショ

降着円盤状態遷移シミュレーションの実施

ンプロジェクト）中島研吾（東京大学・情報基盤センター）

一般相対論的輻射磁気流体コードの実装、陰解法に適用できる並列性能の高い連立１降着円盤への適用次方程式解法の実装、並列性能向上

松本倫明（法政大学・人間環境学部）

松本正晴（東京大学・情報基盤センター）

AMR 法に基づく自己重力磁気流体コードの並列性能向上、星形成領域への適用

陰解法に適用できる並列性能の高い連立１1次方程式解法の実装、並列性能向上

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

2.

研究の目的と意義に対応する降着率（臨界降着率）よりも高い降着

本共同研究の目的は、重力を及ぼす天体の周り率でブラックホールにガスが落下する超臨界降着に形成される降着円盤の時間発展を計算する３次流のシミュレーションに適用し、輻射圧で加速さ元磁気流体シミュレータを改訂して、その適用可れ、磁場によって絞られたジェットが噴出するこ能範囲を拡大することにある。降着円盤では磁気とを示した。また、一般相対論的な輻射磁気流体ストレスに起因する角運動量輸送や輻射と物質のコードのテスト計算を実施した。

相互作用などを考慮する必要があり、磁場と輻射を考慮した輻射磁気流体シミュレーションが必要

5. 今年度の研究成果の詳細になる。我々は、これまでの共同研究によって輻

（１）CANS+コードの整備

射輸送方程式を角度方向に積分したモーメント式

次世代降着円盤シミュレータの基盤となる磁気を差分法によって解く解法を実装してきた。

流体コードの実装、並列計算におけるコードの最本研究では一般相対論化した輻射磁気流体コー適化を進めた。それらを整備し、典型的な天体磁ドをシミュレータに実装し、ブラックホール近傍気流体現象の初期条件、境界条件等をあらかじめのシミュレーションに適用可能にする。次世代の準備してシミュレーションを簡便に実行できるよエクサフロップス計算機を視野に入れて吸収・散うにした基本課題モジュール群と磁気流体方程式

乱を考慮した輻射輸送方程式を角度積分せずに各の時間発展を差分法によって数値的に解くシミュ

光線方向について解く解法を実装すること、シミレーションエンジンをパッケージ化した CANS＋ュレーション結果をもとに輻射スペクトルを計算

http://www.astro.phys.s.chiba-u.ac.jp/cans/

するモジュールを実装することも目的とする。

を公開した。CANS+の特徴として、•近似リーマン解法 HLLD をベースとした有限体積法（Miyoshi & Kusano, 2005）

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

輻射の時間スケールが力学的時間スケールに比MP5 法を用いた高次精度（空間 5 次精度）の数値補間法（Suresh & Huynh, 1998）

べて短い系では輻射場の計算に陰解法を用いるた•め、大規模な連立方程式を解く必要がある。また、数値的単極磁場を取り除くための hyperbolic星形成や活動銀河中心核円盤等、自己重力を考慮

divergence cleaning 法（Dedner et al., 2002）する必要がある場合はポアソン方程式の解法の効が挙げられる。これらを実装することにより降着率を高めることも必要である。本研究では、並列円盤で重要な物理過程である圧縮性磁気流体乱流効率の高い陽的差分法に基づくコード開発に加え、差分法向けに最適化された行列解法ライブラリを開発している東京大学情報基盤センターの中島研吾氏らとの共同研究によって、輻射場、重力場の陰的解法の並列効率を高めることを目指す。

4. 前年度までに得られた研究成果の概要
昨年度まで実施していた共同研究「天体活動現象の輻射磁気流体シミュレーション」では、輻射
図１：CANS+コードによる降着円盤の磁気流体シミュレーション結果。左から右に向かって時間発展している。色はトロイダル磁場成分を表す。上の段は空間 5 次精度の MP5 法、下の段は2 次精度の MUSCL 法による計算結果を示す。

輸送方程式を角度方向に積分した 0 次、1 次のモーメント式を解く 1 次モーメント法（M1 法）に基づく特殊相対論的輻射磁気流体コードを実装して球対称降着流の上限光度であるエディントン光度2学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月を減衰させることなく降着円盤の長時間発展を追輻射４元速度、メトリックである。輻射とガスはうことが可能になった。

吸収・散乱過程により輻射４元力を通じてエネ

図１は降着円盤の時間発展を MP5 法と２次精ルギー・運動量をやり取りする。従ってこれらの度の MUSCL 法で比較したものである。ともに、式を解くことでガスと輻射を無矛盾に決定する事ダイナモ効果による大局的磁場成分の反転（赤→ができる。重力を表すメトリックは時間に依存せ青）が得られているが MP5 法では長時間経過後ず、シュバルツシルト半径で発散しないも円盤構造と磁気乱流が維持されている。

Kerr-Schild メトリックを用いた。
このコードを降着円盤の大局的シミュレーションに適用した結果を図２に示す。ボリュームでン（２）一般相対論的３次元輻射磁気流体コード
我々はこれまで特殊相対論的輻射磁気流体コー

ダリングした色(青-赤)は密度を表し、線は磁力線ドを用いて球対称降着流の上限光度（エディントを表す。初期条件として太陽の 10 倍の質量の無回ン光度）に対応する臨界降着率より高い降着率で転ブラックホール近傍にトーラスを置き、円盤内ブラックホールに質量が供給される超臨界円盤かにはポロイダル磁場が配位されている。円盤内赤らのアウトフローの生成と加速について研究を行

道面の密度は 0.1g/cm3 程度である。計算開始後、ってきた。しかしこの計算ではブラックホールに軸対称モードの磁気回転不安定性が成長し、角運よる重力はニュートン近似を用いて記述していた。 動量を円盤外縁へと輸送する事でガス降着が始まそのためブラックホール近傍の物理、特に円盤のる。計算開始後 t=0.02s 程度で降着円盤はブラッ内縁半径がどこにあるのか、そしてブラックホークホールへ達し、その後準定常的にガスが降着す

ルの回転が円盤構造に影響を与えるのか、といっる。ガス降着率は臨界降着率の 1000 倍程度である。

た疑問には答えられていなかった。また軸対称を円盤内は輻射エネルギーが卓越し、輻射圧、磁気仮定して計算していたため、ガス降着を考える上圧、ガス圧の比は 100:10:1 である。また円盤からで必須である円盤内乱流を正しく扱う事ができなはアウトフローが形成され、その速度は 0.3c 程度いという問題があった。となる（図２, 上下にのびる青い領域）。これらのこれらの問題を解決すべく我々は３次元一般相

結果はこれまでの特殊相対論的輻射磁気流体シミ対論的輻射磁気流体コードを作成した。基礎方程ュレーション結果と矛盾は無い事が確かめられた。

式は以下の通りである。

上から順に質量保存、磁気流体・輻射のエネルギー運動量保存式、誘導方程式である。

,はメトリックの行列式とクリストッフェル記号で、重力やブラックホールの回転の効果を表す。
輻射場は輻射強度を運動量空間で積分したモーメント量によって定式化され、以下の M1 クロージャーを用いて方程式は閉じる

ここで

図 2：一般相対論的輻射磁気流体計算結果。
青—赤のレンダリングデータは密度、線は磁力線、上下にのびる青いボリュームデータはアウトフロー速度を表す。

はそれぞれ輻射エネルギー密度、3

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

また、これまで取り扱う事ができなかった効果された軟 X 線が高温ガスによって逆コンプトン散として３次元効果がある。図３は r=10rg での（rg乱されることにより、観測される硬 X 線を作り出はシュバルツシルト半径）トロイダル磁場の時間すと考えられる。

変化を示す。トロイダル磁場の極性が時間と共にこのように、磁気的な角運動量輸送率増加に伴う変化していることがわかる。これは円盤内でダイガス温度と輻射温度のずれに起因する高温プラズナモ機構が働き、定常的な磁場増幅が起きているマによって硬 X 線起源を説明するモデルはこれまことを示している。このようなダイナモ機構が一でに無く、ガスだけでなく輻射、磁場、重力と行般相対論効果を含めた輻射磁気流体計算でわかっった全ての物理を無矛盾に扱うことで得られた新たのは初めてである。

しい知見である。

図 3: r=10rg でのトロイダル磁場の時間進化。オ
レンジと青はそれぞれケプラー回転方向に対して順方向と逆方向の磁場を表す。

円盤内はガス密度が高く、吸収・散乱に対して光学的に厚いためにガスと輻射が熱平衡状態にあり、その温度は１億度程度となる。ブラックホール近傍ではガス温度と輻射温度の間にズレが生じ、ガス温度は 100 億程度まで上がる(図４左上)。こ
図 4: 赤道面のガス温度分布。左は無回転、右は高の輻射温度とガス温度がずれる半径が円盤の内縁速回転するブラックホールの結果。上下は初期磁半径と見なせる。従来の理論的考察ではこの半径場配位の違い(上はポロイダル磁場、下はトロイダはガスがブラックホールの周りをケプラー回転でル磁場)を示す。

き る 限 界 の 半 径 (Innermost Stable CircularOrbit, ISCO)で決まっていると考えられていた。

（３）自己重力 AMR コード SFUMATO の試験公開

今回の計算の結果、円盤内縁半径は ISCO よりも解適合格子法（AMR）に基づく自己重力流体コー外側にあり、磁場形状に依存する事がわかった。

ド SFUMATO をコミュニティの一部に試験的に公開

子午面内の磁場が強い場合、磁気回転不安定性にした。このコードにはポアソン方程式を、マルチよって方位角磁場が作られる。この方位角磁場とグリッド法を用いて解くモジュールが実装されて子午面磁場によって強い角運動量輸送が働き、ガおり、ガス系の自己重力が重要になる星形成や活スはブラックホールへと素早く落ちることができ動銀河中心核、ブラックホールによる恒星の潮汐る。その結果、ガスと輻射は熱平衡にならず、ブ破壊事象のシミュレーション等に適用することがラックホール近傍に高温ガス雲が形成される。一できる。ユーザの一部では既に科学的な成果が出方トロイダル磁場が優勢な場合、磁気トルクは比始めている。また、開発者が想定していない計算

較的弱く、ガスと輻射は熱平衡になりやすいため、モデルに利用されることにより、SFUMATO コード

円盤内縁半径は内側へと移動する事が解った。の仕様に不具合があることが報告された。これら

高効率の角運動量輸送によって高温ガス雲がブのユーザからのフィードバックにもとづいて、ラックホール近傍に形成されると、円盤から放射SFUMATO コードの仕様の変更を検討した。
4学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月

SFUMATO コードを用いて実施した集団的星形成が高く、降着時間スケールが短くなる場合、ブラのシミュレーションでは、ベースグリッドが 512^3ックホール近傍においてガス温度が輻射温度よりメッシュという比較的大規模な計算においても、もかなり高くなる領域が出現することが明らかに現実的な計算が可能であることを示した。また、なった。ブラックホールの回転の効果を取り入れシンク粒子法を実装し、シンク粒子数が 1000 個程たコードも実装され、ブラックホールの回転によ度の場合でも現実的な計算時間で計算を実施可能

る有効ポテンシャルの変化がジェットやアウトフにした。分子雲コアの分裂のシミュレーションでローのダイナミクスに及ぼす影響を調べることがは、乱流をともなった分子雲コアが重力収縮中に可能になった。

分裂して多重星系を形成し、多重星の周囲に擾乱シミュレーション結果に基づく輻射スペクトルが伝播してアーク状の構造が形成されることを確計算について、現在、効率の良い疑似乱数の生成認した。
このようなアーク状の構造は、最近の ALMAアルゴリズムである Mersenne twister 法をこれ望遠鏡による観測で発見された。このシミュレーまでわれわれが用いてきたモンテカルロ法に基づションでは自己重力 AMR の特性を活かして高解像

く特殊相対論的輻射スペクトル計算コードに実装を可能にしただけではなく、前述のシンク粒子法中である。また、輻射磁気流体シミュレーションによって長時間の進化を追跡可能にした。

結果を用いて輻射スペクトルを計算する際に、シ連星のガス降着のシミュレーションでは、連星ンクロトロン放射を考慮する計算モジュールを実を取り囲む円盤である周連星円盤を調べた。連星装してテスト中である。今後は、このコードをに向かって落下するガスは角運動量を持っているCANS+の解析用モジュールとして整備し、公開してため、直接連星に落下せずに、一旦周連星円盤にいく予定である。

落下する。連星からの重力トルクによって、周連ブラックホールの近傍では一般相対論的効果を星円盤には渦状腕が形成される。このような渦状考慮した輻射スペクトル計算をすることが必要不腕は、最近の ALMA 望遠鏡による観測で発見された。 可欠である。今年度は一般相対論的輻射スペクト
ル計算コードの開発に着手し、測地線方程式を 4

このシミュレーションでは AMR の特性を活かし、従来よりも高解像なシミュレーションを実現した。 次精度 Runge-Kutta 法で解く計算モジュールを作成した（図５参照）。今後、より効率的に正確な

6.

解 を 求 め る た め に 、 adaptive stepsize control

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

開発した高次精度の磁気流体コードは CANS+

algorithm により積分の stepsize を変化させながとして整備し、公開されている。また、本年度千ら計算する embedded Runge-Kutta method に葉大で開催した「宇宙磁気流体・プラズマシミュ変更する。この embedded Runge-Kutta method

レーションサマースクール」で CANS+を使ったに は い く つ か 種 類 が あ る の で (Fehlberg 1970,

演習を行い、シミュレーション初心者にも使いや

Dormand & Prince 1980, Cash & Karp 1990)、

すいコードとして利用されている。今後は、本研実際に測地線方程式を解いて性能評価を行い、最究グループで開発した輻射磁気流体方程式のソル適な方法を採用する。今後はこの計算モジュール

バー等を CANS+に実装して公開することにより、をこれまでに用いてきた輻射スペクトル計算コーコミュニティの幅広いニーズに応えていく予定でドに実装し、一般相対論的輻射スペクトル計算コある。

ードを完成させブラックホール降着流の輻射スペ昨年度から開発してきた一般相対論的輻射磁磁クトルを計算する予定である。

気流体コードを用いて３次元計算を行うことで、輻射輸送方程式を各光線方向について解く直接ブラックホール超臨界降着円盤の３次元構造がわ

解法コードの実装とテスト計算も進展している。

かってきた。特に、磁場による角運動量輸送効率5学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月

7.

研究成果リスト(1) 学術論文
・H.R. Takahashi, K. Ohsuga, ‘Radiation Drag
Effects in Black Hole Outflows fromSuper-critical Accretion Disks via Special
Relativistic Radiation
Magnetohydrodynamics Simulations’, Publ.
Astron. Soc. Japan 印刷中・K. Hashizume, K. Ohsuga, T. Kawashima, and M.
Tanaka, ‘Radiation Hydrodynamics
Simulations of Wide-angle Outflows fromSuper-critical Accretion Disks around Black
Holes’, Publ. Astron. Soc. Japan 印刷中・S. Hirose, ‘Magnetic turbulence and

図 5：反時計回りに回転するブラックホール近傍の光子軌道のテスト計算結果。長さの単位はブラックホール質量。回転ブラックホールによる時空の引きずりにより光子は回転しながらブラックホールに吸い込まれていく。

thermodynamics in the inner region of
protoplanetary discs’, Mon. Not. Roy.
Astron. Soc. 448, 3105-3120, 2015
・T. Matsumoto, T. Onishi, K. Tokuda, S.

自己重力 AMR コード SFUMATO について、今後はInutsuka,‘An origin of arc structures

磁場を考慮したシミュレーションを遂行し、星形

deeply embedded in dense molecular cloud

成における磁場の役割について定量的に調べる。

cores’, Mon. Not. Roy. Astron. Soc. Letters,

磁場はすでに SFUMATO コードに実装済である。分

449, L123-L127, 2015

子雲スケールのシミュレーションである集団的星・T. Matsumoto, K. Dobashi, T. Shimoikura,

形成のシミュレーションでは、とくに両極性拡散‘Star Formation in Turbulent Molecularによる磁場の散逸が重要になる。そこで SFUMATO

Clouds with Colliding Flow’, Astrophys. J.

への両極性拡散の実装についても今後検討する。

801, article id. 77 (21pages), 2015

SFUMATO コードの公開にともなうユーザからのフ・Y. Asahina, T. Ogawa, T. Kawashima, N.

ィードバックを検討し、SFUMATO コードをさらにFurukawa, R. Enokiya, H. Yamamoto, Y. Fukui,発展させるべく、開発を継続する。

and R. Matsumoto, ‘Magnetohydrodynamic

以上のように陽的差分に基づく磁気流体コード、Simulations of a Jet Drilling an HI Cloud:

ポアソンソルバーを実装した自己重力流体・磁気

Shock Induced Formation of Molecular Clouds

流体 AMR コード、回転ブラックホールに適用可能

and Jet Breakup’, Astrophys. J. 789, article

な一般相対論的輻射磁気流体コード、シミュレー

id. 79 (11pages), 2014

ション結果に基づく輻射スペクトル計算コードの・S. Hirose, O. Blaes, J.H. Krolik, M.S.B.

開発と降着円盤シミュレータへの実装は順調に進Coleman, T. Sano,‘Convection Causes

んでいる。今後は、super time stepping 法を用

Enhanced Magnetic Turbulence in Accretion

いた陽的な熱伝導モジュールの実装、輻射輸送方

Disks in Outburst’, Astrophys. J. 787,

程式の直接解法の実装と、陰的解法の高速化につarticle id 1(14pages), 2014

いての研究を加速する。

・K. Dobashi, T. Matsumoto, T. Shimoikura, H.

Saito, K. Akisato, K. Ohashi, K. Nakagomi,
6

学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月‘Colliding Filaments and a Massive Dense(4) 国内会議発表

Core in the Cygnus OB 7 Molecular Cloud’,・川島朋尚、大須賀 健、松元亮治「 巨大ブラッAstrophys. J. 797, article id. 58 (22pages),

クホールによる潮汐破壊事象の磁気流体・輻射2014

流体数値実験 」地球惑星科学連合 2014 年大会、・S. Takakuwa, M. Saito, K. Saigo, T. Matsumoto,パシフィコ横浜、2014 年 5 月

J. Lim, T. Hanawa , P.T.P. Ho,, ‘Angular

・廣瀬重信、「熱対流によって強められる磁気乱流」、

Momentum Exchange by Gravitational Torques

乱流研究のフロンティア(京大 数理解析研究所

and Infall in the Circumbinary Disk of theRIMS 研究集会)、京都大学数理解析研究所, 2014

Protostellar System L1551 NE’, Astrophys. J.

年 7 月 23-25 日
・高橋博之、'Relativistic Magnetic

796, article id. 1 (17pages), 2014

Reconnection'、 超新星/ガンマ線バースト研究(2) 国際会議プロシーディングス

会 2014 、理化学研究所、2014 年 8 月 25 日-27
日 (invited)(3) 国際会議発表

・高橋博之、大須賀健、「超臨界降着円盤の一般相・T. Kawashima, T. Ogawa, T. Ono, H. Oda , R.

対論的輻射磁気流体シミュレーション」、日本天Matsumoto. 'X-ray Flares Induced by the

文学会秋季年会、山形大学、2014 年 9 月 11 日

Magnetic

Reconnection

inBlackHole

・松本倫明 , 土橋 一仁, 下井倉ともみ、「乱流とAccretion Flows', MR2014, Tokyo, Japan, May,衝突をともなう高密度分子雲における集団的星

20-22, 2014

形成」
、日本天文学会 2014 年秋季年会、山形大･ R. Matsumoto, Y. Kudoh, Y. Asahina, M.

学，2014 年 9 月

Machida, T. Kawashima, ‘Global・川島朋尚、松元亮治、「 ブラックホール降着流

Magnetohydrodynamic Simulations of Diskを通過するガス雲による降着率変化」日本流体Dynamos’, MR2014, Tokyo, Japan, May 20-22,

力学会年会 2014、東北大学、2014 年 9 月
・高橋博之、大須賀健、‘Numerical Study of2014
・T. Kawashima, T. Ogawa, T. Ono, H. Oda,

Supercritical Accretion Disks’, Plasma

R. Matsumoto, 'Three-Dimensional MHD

Conference 2014, 朱鷺メッセ、新潟、2014 年

simulations of Disk Dynamos', 4th East Asian

11 月 20 日

School and Workshop on Laboratory, Space,・川島朋尚、松元亮治、「ブラックホール降着流とAstrophysical Plasmas, Harbin, China, July

ガス雲の相互作用による質量降着率と磁場の時

28-August 1, 2014間変化」, Plasma Conference 2014, 朱鷺メッ・ S. Hirose, ‘Radiation MHD Simulations ofthAccretion Disks’, in 4

セ、新潟、2014 年 11 月 21 日

East-Asian School・高橋博之、大須賀健、関口雄一郎、川島朋尚、Space,「ブラックホール降着円盤の３次元一般相対論

Astrophysical Plasmas, Harbin Institute of的輻射磁気流体シミュレーション」、第 27 回 理

Technology, July 28-August 1, 2014

論懇シンポジウム、2014 年 12 月 24 日-26 日andWorkshop

onLaboratory,・ S. Hirose, ‘Magnetic turbulence in innerradiiofprotoplanetarydisks’,・松本倫明, 土橋 一仁, 下井倉ともみ、「乱流とin

衝突をともなう高密度分子雲における集団的星

Non-ideal MHD, Stability, and Dissipation in形成 II」、日本天文学会 2015 年春季年会、大阪

Protoplanetary Disks, Niels Bohr Institute,大学、2015 年 3 月
・ 松本倫明, 徳田一起, 大西利和, 犬塚修一郎、

August 4-8, 20147

学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月「分子雲コアの中心に存在するアーク構造の起源」日本天文学会 2015 年春季年会、大阪大学、2015 年 3 月
・廣瀬重信, FU Ori アウトバーストと原始惑星系
円盤の熱不安定について, 日本天文学会春季年会, 大阪大学, 2015 年 3 月 19 日
・川島朋尚、松元亮治、「ガス雲の通過に伴う銀河中心ブラックホール高温降着流の光度変化」日
本天文学会 2015 年春季年会、大阪大学、2015
年 3 月 20 日
・高橋博之、大須賀健、関口雄一郎、川島朋尚、「一般相対論的輻射磁気流体計算で探る高降着率円盤構造」、日本天文学会春季年会、大阪大学、2015 年 3 月 21 日
(5) その他（特許，プレス発表，著書等）（プレス発表）双子の赤ちゃん星を育むガスの渦巻き、2014 年
12 月 4 日、高桑繁久、斎藤正雄、西合一矢、松本倫明、Jeremy Lim, 花輪知幸、Paul T.P. Ho、国立天文台

8