太陽将来計画ワークショップ
2006年6月29日 三鷹

太陽数値シミュレーションの将来
横山 央明 ・ 磯部 洋明（東京大学地球惑星）

太陽シミュレーションの歴史 浮上磁場を例に 1/2
前スペース観測時代
Shibata (1980)

数十点

数＋Mflops

1次元流体

1万点
「ひのとり」時代
Shibata et al. (1989) 2次元MHD(Cf. P4 4G 数Gflops)「ようこう」時代「Solar-B」さきがけ
Matsumoto et al. (1998) 3D

Shibata et al. (1992) TY& Shibata (1995) リコネクション

600Mflops 10万点

30Gflops 5十万点

太陽シミュレーションの歴史 浮上磁場を例に 2/2「SOLAR-B」直前

Isobe et al. (2005) 3Dリコネクション 地球シミュレータ
2億点
数＋Tflops
Miyagoshi et al. (2006)

3千万点

Tflops
Notoya et al. (2006) 3D プラズマ放出

Tflops3千万点

スーパーコンピュータ環境について• これまで– 数十Mflopsから数十Tflopsへ（6桁）– 数十点（1980）から数億点（現在）へ7桁••••

国立天文台 VPP 数百GFlops →2008年更新予定
宇宙研SX6 TFlops
地球シミュレータ（ES）数十TFlops
理研 特定高速計算機（「京速計算機」） ペタFlops

議論されたおもなテーマ
CME

コロナ加熱・太陽風

松本

櫻井・鈴木
フレア
リコネクション
草野・新田・横山
粒子加速今田

エネルギー蓄積・トリガ
草野・宮腰・井上・能登谷
塩田・西田
真柄・野澤・宮腰・清水

浮上磁場

放射流体
廣瀬・川道
工藤

彩層・スピキュール
黒点
高橋・磯部

光球・対流 磯部

議論されたおもなテーマ
CME

コロナ加熱・太陽風

松本

櫻井・鈴木
フレア
リコネクション
草野・新田・横山
粒子加速今田

エネルギー蓄積・トリガ
草野・宮腰・真柄・井上・能登谷
塩田・西田・清水
宮腰・野澤

浮上磁場

放射流体
廣瀬・川道
工藤

彩層・スピキュール
黒点
高橋・磯部

光球・対流 磯部

浮上磁場と対流との相互作用（磯部・宮腰）• 今の浮上磁場の３次元計算はどれも小さい、時間も短い。
• 対流が入るとあっという間に壊れる。(Isobe& Shibata 2004、Cheung et al. 2005)
深く大きな計算領域で計算

黒点構造のモデリング、黒点形成（磯部・高橋）

光球下の音速分布、速度場(Kossovichev)

対流コラプス

光球下の磁場構造は？

radiative cooling

Photosphere
cold
downflow
flux tube
Inhibition of convection

intense flux tube
Thomas et al. 2002 半暗部の構造の起源:pumping?

議論されたおもなテーマ
CME

コロナ加熱・太陽風松本

櫻井・鈴木
フレア
リコネクション
草野・新田・横山
粒子加速今田

エネルギー蓄積・トリガ
草野・宮腰・井上・能登谷
塩田・西田真柄・野澤・宮腰・清水

浮上磁場

放射流体
廣瀬・川道
工藤

彩層・スピキュール
黒点
高橋・磯部

光球・対流 磯部

磁気対流とコロナ加熱（磯部）•Katsukawa & Tsuneta 2005高温ループの方がフィリングファクターが小さい。磁束管が自由に動けるから？

Isobe (2006) 2次元シミュレーション

Steiner et al. 1998

3次元シミュレーションに拡張して、ねじりAlfven波発生の効果を取り入れるコロナ加熱・太陽風（鈴木）議論されたおもなテーマ
CME

コロナ加熱・太陽風

松本

櫻井・鈴木
フレア
リコネクション
草野・新田・横山
粒子加速今田

エネルギー蓄積・トリガ
草野・宮腰・井上・能登谷
塩田・西田
真柄・野澤・宮腰・清水

浮上磁場

放射流体
廣瀬・川道
工藤

彩層・スピキュール
黒点
高橋・磯部

光球・対流 磯部• 京大グループでの役割分担– 浮上磁場、エネルギー蓄積 (清水)– フレア・CMEのトリガー、エネルギー解放(西田)（塩田）– 太陽風、惑星間空間の伝播 (松本)• 最終的にはこれらを統合して計算する• 下層での計算結果を、上層の初期条件・境界条件として計算する• 同時にすべてをとく必要はない

フィラメント放出過程（能登谷）
リコネクションした磁力線は周囲の磁力線を引き伸ばして上昇していく。

リコネクションに伴って電流シートが散逸される。

データ駆動型モデリング• データ同化が基づく数値シミュレーション
数値シミュレーション（MHD)
観測データ
ベクトル
マグネトグラム

初期境界条件
初期境界条件
初期境界条件

計算結果
計算結果
計算結果

観測データ
X、紫外、Hα• 同化手法の開発– 3次元コロナ磁場の同化（初期条件）– 光球面電場の同化（境界条件）• シミュレーションモデルの開発– 活動領域から全球コロナへのmulti-resolutionモデル

議論されたおもなテーマ
CME

コロナ加熱・太陽風

松本

櫻井・鈴木
フレア
リコネクション
草野・新田・横山
粒子加速今田

エネルギー蓄積・トリガ
草野・宮腰・井上・能登谷
塩田・西田
真柄・野澤・宮腰・清水

浮上磁場

放射流体
廣瀬・川道
工藤彩層・スピキュール
黒点
高橋・磯部

光球・対流 磯部

Tanuma et al. (2001)「乱流」リコネクション
MHDマクロスケールと運動論的ミクロスケール「大域的電流シート」の中間スケール
リコネクションの間欠性・非定常性~104 km
Tajima &Shibata (1997)~10 mzy
横山(2001)xBx

連結階層リコネクションモデル（杉山、草野）E  EMHD  (1   )EPIC

電磁流体（MHD）モデル
一般化されたオームの法則
EMHD  U  B  0

B   E
t
1 EPIC   B  0 J PIC
2c t

マックスウェル方程式
粒子（PIC）モデル

議論されたおもなテーマ
CME

コロナ加熱・太陽風

松本

櫻井・鈴木
フレア
リコネクション
草野・新田・横山
粒子加速今田

エネルギー蓄積・トリガ
草野・宮腰・井上・能登谷
塩田・西田
真柄・野澤・宮腰・清水浮上磁場放射流体
廣瀬・川道
工藤

彩層・スピキュール
黒点
高橋・磯部

光球・対流 磯部輻射流体と光球・彩層のダイナミクス（川道ほか京大グループ）• フレアループの1次元輻射流体シミュレーション– 分光・偏光観測との比較により、非熱的粒子、彩層ダイナミクスの診断– White light flare– 恒星フレアへの応用• 彩層ダイナミクス・加熱– 浮上磁場、Ellerman bomb
– MHD波動、リコネクション• MHD＋部分電離– 光球・彩層リコネクション– 浮上磁場 (Leake & Arber 2006)– 原始星降着円盤への応用• (M)HD＋電離非平衡イオン– リコネクションによる加熱(加速)時間＜電離平衡に達する時間– 電離状態が異なる＝＞電離平衡を過程した分光診断に影響あり（Solar-B/EISの観測の解釈）– 1Dフレア／コロナループ(Bradshaw & Mason 2003)

「太陽シミュレーションの今後」キーワード• 浮上磁場・黒点形成・コロナ波動加熱• フレア・CMEエネルギー蓄積 データ連携シミュレーション• 階層間結合 リコネクションなど• マルチ物理輻射、電離非平衡など

東京大学横山グループの今後 1/2• Solar-Bを見据えた研究– 波動による加熱• 磁気対流による波動発生• 波動によるコロナループ加熱の多次元シミュレーション– 磁気対流・浮上磁場• 黒点・ポアの形成・崩壊• デルタ型黒点（活動領域）の形成– フレア・プラズマ放出• 浮上磁場によるコロナ磁場の非平衡・不安定化– 磁気リコネクションの基礎研究• 「乱流」リコネクション

「太陽シミュレーションの今後」キーワード• 浮上磁場・黒点形成・コロナ波動加熱• フレア・CMEエネルギー蓄積 データ連携シミュレーション• 階層間結合 リコネクションなど• マルチ物理輻射、電離非平衡など

ここまでは、現在のサイエンス（とSOLAR-B）の延長
何か欠けている、、、

太陽物理学の重要課題• フレア• コロナ加熱• ダイナモ！！！！！

太陽物理学の重要課題• フレア済ようこう– エネルギー解放ーリコネクション– エネルギー蓄積・トリガ– 粒子加速• コロナ加熱済？ Solar-B• ダイナモ

東京大学横山グループの今後 2/2• Solar-Bの先– 「より深く」「より遠くへ」をめざす• 対流からダイナモへ– 対流層中での浮上磁場– 回転勾配層– 天体ダイナモ 銀河、降着円盤• 他天体へ「太陽物理学は天文学」– 特殊相対論・一般相対論– 放射流体力学 はじめはFLD近似• 粒子加速– Fokker-Planckシミュレーション

Miesch et al. 2006

Brummel et al. 2002

MHD＋輻射以外のマルチフィジクスを含むシミュレーション• MHD＋部分電離– 光球・彩層リコネクション– 浮上磁場 (Leake & Arber 2006)– 原始星降着円盤への応用• (M)HD＋電離非平衡イオン– リコネクションによる加熱(加速)時間＜電離平衡に達する時間– 電離状態が異なる＝＞電離平衡を過程した分光診断に影響あり（Solar-B/EISの観測の解釈）– 1Dフレア／コロナループ(Bradshaw & Mason 2003)

3D シミュレーション• Gudiksen and Nordlund (2005)– model granulation as a driverFuture Work
今までの解析結果より、厳密なForce Free解であるならば、RelaxationMethodを用いてほぼ収束する事はだいたいわかった。

実際の太陽コロナ（１） 本当にForce-Freeなの？？ 速度場の影響は？？
（２） そもそも有限β効果を無視してもいいのか？？ 得られる観測磁場データは光球面上のもので、ここはβ=0の近似は使えない。
厳密なForce-Free解（Low & Lou）での検証
この場合は基本的に収束する。我々も検証済み。

数値シミュレーション
結果を用いての検証（非平衡度合い、有限βの効果など）
このSTEPは難しいと思う。

光球面磁場データを用いての検証
この場合は収束しない。
上の理由もあるし、そもそもデータが粗いので、Force-Free解かどうかが怪しい。

フラックスロープ構築
Amari et al. 2003
bipolar potential field + twist motion (bottom boundary)Realistic simulationの最近の成果• 深い計算領域でマルチスケール対流(Benson, Stein,Nordlund)• 対流＋磁束管浮上(Cheung, Schuessler, MorenoInsertis)

q=0.005 (捻れ弱)Miyagoshi et al. (2006)
コロナ上空およそ14000kmへ浮上後の磁力線構造（俯瞰図）

・計算の空間時間スケール(～2万km、１時間)を現実の活動領域形成（5万km～、数日）に近づける---太陽表面磁場の時間発展を観測と詳細比較

・熱的プロセス（エネルギー式：断熱）をリアリスティックに---浮上磁束によるコロナループ形成---放射冷却による黒点形成

高磁気レイノルズ数領域のリコネクション

MHD性• 多重リコネクション• 乱流リコネクション• Hall MHDリコネクション– Whistler wave
• Collision-less リコネクション– Meandering motion• 異常抵抗– LHDI, Ion Acoustic Instability, etc. 非MHD性マルチスケール問題としてのリコネクション

Instability
非MHD
Slow shock
Fast shock

乱流
meandering

Hall MHD
Whistler wave

LHDI, Ion Acoustic Inst., etc.

MHD