天体MHD数値シミュレーションの魅力と魔力.ppt

天体ＭＨＤ数値シミュレーションの魅力と魔力 柴田一成 京大理花山天文台 本講演の内容 １．はじめに ２．天体ＭＨＤ現象のおもしろさ ３．天体ＭＨＤシミュレーションは超困難 ４．天体ＭＨＤシミュレーションの魅力 ５．天体ＭＨＤシミュレーションの魔力。はまると危険（落とし穴の数々） ６．むすび：ノーベル賞課題（超難問）に挑戦せよ １．はじめに 　天体ＭＨＤシミュレーションとは？ 　　 　　　　天体現象への応用を目的として、 　　　　時間に依存する電磁流体（ＭＨＤ） 　　　　方程式を数値的に解くこと。 　 　　シミュレーション＝数値実験 　　　 　　　　　　　　 天体ＭＨＤシミュレーションを取り巻く迷信?偏見の数々 　１）シミュレーションは何でも解決してくれる 　２）シミュレーションは、境界条件をうまく選　べば必ず望む結果を出してくれる 　３）シミュレーションは、簡単である。頭が悪い 　４）シミュレーションは、きれいな絵やムービー　を見せて人をだまそうとする手段である 迷信にまどわされるな １）ＭＨＤシミュレーションは、万能ではない。 ２）ＭＨＤシミ ュレーションでは、境界条件をコントロールするのは至難のワザ。 ３）ＭＨＤシミュレーションは、流体シミュレーションの１０倍くらい困難。精神力、体力が必要。 ４）シミュレーション?ムービーは、実におもしろくて役に立つ。人（または自分）を「だます」のではなく、「教育する」のである。 私の場合 大学院に入りたての頃、コンピュータにも、シミュレーションにも全く興味はなかった。 　ＦＯＲＴＲＡＮも知らなかった。 そんな私が、いかにして、天体ＭＨＤシミュレーションに、はまっていったか、本日はその魅力と魔力について語りたい ２．天体ＭＨＤ現象のおもしろさ 天体ＭＨＤ現象の代表例 　　　宇宙ジェット（ＡＧＮ、連星系、原始星） 　　　太陽フレア?ジェット　　 まずは、魅惑あふれる画像とムービーを見よう 宇宙ジェット（活動銀河核ジェット：電波銀河 Cyg A　with VLA） 宇宙ジェット（原始星ジェット：　ＨＨ１－２with HST） 宇宙ジェット（連星系ジェット：ＳＳ４３３ with ASCA） 宇宙ジェットと降着円盤（活動銀河核ジェット） 宇宙ジェットと降着円盤（原始星ジェット） 太陽ジェット（Ｈα）（京大飛騨天文台： Kurokawa et al. 1988） 太陽コロナ?フレア?ジェット（ようこう軟Ｘ線） 太陽フレアからの衝撃波（京大飛騨天文台：Ｈα） 太陽コロナ質量放出（ＳＯＨＯ／ＬＡＳＣＯ　白色光） かに星雲中の波（ＨＳＴ） ３．天体ＭＨＤシミュレーションは超困難　（でも、だから、やりがいがある） 電磁流体力学の方程式がそもそも複雑　８変数の非線形偏微分方程式　　ｃｆ）流体力学の方程式は５変数 電気抵抗が入るとさらに複雑　　磁気リコネクションの物理は未だ解明　　されていない 電磁流体（ＭＨＤ）方程式 　　未知数８個：　密度（ρ）、速度ベクトル（ｖ）、 　　　　　　　　　　　　磁束密度ベクトル（Ｂ）、圧力（ｐ） 　　方程式８個：　非線形連立偏微分方程式 　　 　 数値計算の方法 差分法：　微分を差分で近似。 有限の格子点上　　　　　　で計算 粒子法：　超粒子の運動方程式を解く 　　　　　　　　ラグランジュ法。ＭＨＤには不向き。 　　　　　　　　（磁場を解く必要があるので、　　　　　　　　　　　　　格子点はどうしても必要） 　　 電磁流体波（スローモード）も奇妙 天体ＭＨＤシミュレーションはどこが難しいか？ 境界条件（ＭＨＤでは自由境界は、きわめ　　　　　　　て困難） 密度成層（アルフベン速度が大きく変化） 　　　例）太陽光球密度／コロナ密度＝ 初期条件（初期平衡解を求めることさえ、　　　　　　　　　大仕事） 　　＝＞数値不安定の嵐！ 太陽ジェットの磁気リコネクション?モデル（Yokoyama and Shibata 1995)密度分布 ４．天体ＭＨＤシミュレーションの魅力　　 とにかく、ムービーはおもしろい。 「人をだます」などと言わずに、素直に感動しよう。 １）　横山：　太陽フレア?ジェットの磁気リコネ　　　　　　　クション?モデル（天文学会奨励賞） ２）　　林：　　原始星フレアの磁気リコネクション?モデル（SGI コンピュータグラフィックス賞） ３）　工藤：　宇宙ジェットのＭＨＤモデル ４）　田沼：　銀河における磁気リコネクション　 太陽ジェットの磁気リコネクション?モデル（Yokoyama and Shibata 1995)温度分布 太陽コロナＸ線ジェット 足元でマイクロ?フレア 長