2023 Fiscal Year Research-status Report
Revealing the Solar Wind and Corona Heating Problem through Extreme Ultraviolet Emission from plasma in non-equilibrium ionization
| Project/Area Number |
23K03456
|
|---|
| Research Institution | Nagoya University |
|---|
Principal Investigator |
松本 琢磨 名古屋大学, 宇宙地球環境研究所, 特任助教 (50728326)
|
|---|
| Project Period (FY) |
2023-04-01 – 2026-03-31
|
| Keywords | 太陽遷移層 / 太陽コロナ / 紫外線放射 / 非平衡電離 / 解適合格子 |
|---|
| Outline of Annual Research Achievements |
太陽の彩層とコロナの間には、厚さが約0.1Mm以下で、遷移層と呼ばれる極めて薄い層が存在する（コロナの圧力スケール長は約100Mm）。この層では、温度帯によっては、C IVやN V、O VIなど、LiやNaに似たイオンからの紫外線輝線が放射される。太陽コロナからの紫外線放射については、基本的には電離平衡を仮定した理論を使用して観測をよく再現できるが、遷移層からの輝線に関しては、数倍の過大または過小評価がなされることが知られている。
この知見を踏まえ、2023度は、彩層とコロナの間で起こる蒸発や凝縮過程に伴う電離平衡からの乖離に注目し、太陽大気のMHDシミュレーションと同時に、電離の時間発展（移流、電離、再結合）を解くシミュレーションコードを開発することを目標とした。初めに、非平衡電離の効果を導入した1次元のMHDシミュレーションを行った結果、蒸発に伴いLi-やNa-likeイオンの異常な増光が観察された（Matsumoto 2024, ApJ）。さらに、この計算から、非平衡電離の効果を正確に取り扱うには、大局的なコロナの動的現象（約100 kmスケール）と、遷移層での移流（約1km以下）を同時に解くことが必要であることが明らかになった。したがって、最初に想定していた単純な格子を用いた2D MHDシミュレーションでは、遷移層での非平衡電離を解くのが困難であると判断した。
そのため、解適合格子を使用して遷移層構造を空間的に分解する方針に転換し、研究を継続することにした。現在は、解適合格子上で1Dから3DまでのMHDシミュレーション、一様重力、grey近似輻射輸送、潜熱を考慮した状態方程式、Super Time Stepを用いたSpitzer熱伝導、電離の時間発展式を実装し、そのデバッグとテスト計算を進めている。
|
| Current Status of Research Progress |
Current Status of Research Progress
2: Research has progressed on the whole more than it was originally planned.
Reason
2023度は、太陽遷移層からの異常な紫外線放射を解明する糸口として、非平衡電離を伴う1次元MHDシミュレーションコードの開発と数値実験を行い論文として発表した（Matsumoto 2024 ApJ)。この研究では、太陽光球で対流運動によって励起されるアルフベン波が太陽大気を100万度まで加熱しコロナを形成するという既存のシミュレーションモデル(Moriyasu & Shibata 2004) を採用した。このモデルに新たに非平衡電離を計算するモジュールを追加し、コロナの形成・崩壊に伴い、電離平衡からずれたイオンが放射する紫外線がどのように観測できるかを議論した。モデル中には長期的な蒸発・凝縮過程や、衝撃波と伴う短期的な温度・密度変動が含まれる。その結果、蒸発時には観測と整合的なLi-, Na-likeイオンからの電離平衡で予測されるよりも異常な増光がみられた。また、逆に凝縮時にはLi-, Na-likeイオンからの異常な減光がみられた。しかしながら、短期的な変動である衝撃波に伴う密度・温度変動では、短期的には電離平衡からのずれがみられるものの、時間平均すると平衡からのずれはほぼ消えてしまうことも判明した。これらのことから、Li-, Na-likeイオンの異常な増光・減光は彩層とコロナとの間の質量輸送に関連していることが示唆される。したがって本研究を進めることで、Li-, Na-likeイオンの紫外線観測から、質量輸送の情報が間接的に得られることが期待できる。
|
| Strategy for Future Research Activity |
2023年度に行った1次元計算から、非平衡電離計算には非常に薄い遷移層を十分に空間分解して移流計算を解くことが必要であることが明らかになった。しかしながら、分解すべき最小空間スケールが0.1km程度以下であるのに対して、太陽コロナの構造は数万km程度となっている。したがって、当初想定していた単純な格子を用いた多次元 MHDシミュレーションでは、遷移層での非平衡電離を解くのは計算コストの観点から困難であることが判明した。
そのため本研究では、必要な個所に必要な空間分解能を動的に割り当てることが可能な解適合格子を使用して遷移層構造を空間的に分解する方針に転換した。過去の研究では1次元計算に対して解適合格子を適用し計算の実行速度を向上させる試みがなされてきたが、これを多次元に拡張した例は少ない。現在は、解適合格子上で1Dから3DまでのMHDシミュレーション、一様重力、grey近似輻射輸送、潜熱を考慮した状態方程式、Super Time Stepを用いたSpitzer熱伝導、電離の時間発展式を実装し、そのデバッグとテスト計算を進めている。
2024年度ではこのコードを用いて表面対流層からコロナ上空までを含めた2次元MHDシミュレーションを行い、遷移層を空間分解しつつ非平衡電離を解くことに挑戦したい。多次元化することで遷移層のうねりに伴う視線方向の重なりなどの効果によって、どれくらい紫外線強度が変化するのかを定量的に評価したい。また、2028年打ち上げ予定のSolar-Cの観測に向けて、紫外線放射synthesisの研究も進めていく予定である。
|
| Causes of Carryover |
2023年度に行った1次元シミュレーション計算の結果から、多次元計算の実行には、当初予定していたシミュレーションに新たに解適合格子の実装することが必要であることが判明した。したがって2023年度では国際会議での発表や国外研究者との打ち合わせよりも、シミュレーションコードの開発に注力した結果、渡航費分の当該助成金が生じた。現在はコードのデバッグおよびテスト計算を行っており、2024年度ではこのコードを用いた初期結果をもとにして、国際会議での発表や国外での打ち合わせを実施したい。
|
