2023 Fiscal Year Research-status Report

数値シミュレーションを用いた土星の環のカッシーニ間隙形成過程の研究

Research Project

Project/Area Number	20K04047
Research Institution	Kyoto Women's University
Principal Investigator	道越秀吾京都女子大学, データサイエンス学部, 准教授 (60572229)
Project Period (FY)	2020-04-01 – 2025-03-31
Keywords	太陽系 / 小惑星 / リング / 数値シミュレーション
Outline of Annual Research Achievements	土星探査機カッシーニの観測により、土星の環の構造が驚くほど多様であることが明らかになった。中でも、カッシーニ間隙と呼ばれる領域には、未だ解明されていない数多くの構造が存在している。本研究の目的は、数値シミュレーションを用いてカッシーニ間隙の形成過程を検証し、その中に存在する構造の起源を解き明かすことである。土星の環を完全にシミュレーションするためには、膨大な数の粒子と非常に長い時間スケールが必要となる。現在の計算機の性能では、不可能であるため、シミュレーションコードを十分に最適化し、現実の物理的性質を損なわないように適切なスケール変換を行って考察する必要がある。 2020年度には、シミュレーションコードの改良と最適化に取り組み、2021年度には、土星の環と衛星の相互作用に関する数値実験を行った。これまでにも土星の環と衛星の相互作用に関するシミュレーションは行われてきたが、自己重力を考慮した3次元N体シミュレーションは本研究が初めてであった。本研究では、3次元N体シミュレーションを実施し、カッシーニ間隙の形成に関係すると考えられるリンドブラッド共鳴による密度波の再現に成功した。また、シミュレーション結果を線形摂動解析の理論研究と比較したところ、非定常構造やエヴァネッセント領域（波動が伝播しない領域）の発生などの定性的な理論予測と一致することが確認された。しかし、十分な空間精度で長時間のシミュレーションを行うことは現状では困難であり、さらなるコードの高速化が必要であることが明らかになった。 2022年度には、粒子衝突処理と時間積分のアルゴリズムを工夫することで、シミュレーションの高速化に取り組んだ。2023年度は、アルゴリズムの理論的考察をさらに進め、定式化の正当性を確認した。現在は、さらなるテスト計算による検証を進め、論文の執筆を行っている段階である。
Current Status of Research Progress	Current Status of Research Progress 3: Progress in research has been slightly delayed. Reason 本研究の進捗状況は、当初の計画と比較するとやや遅れが生じている。その主な理由は以下の通りである。初年度においては、計算コードの開発とリファクタリングを行い、拡張性と高速化を実現することができた。また、カッシーニ間隙とリンドブラッド共鳴に関する先行研究を精査し、現実的なタイムスケールでの間隙形成シミュレーションを行うためのパラメータを検討した。選定したパラメータを用いてシミュレーションを実施した結果、リンドブラッド共鳴による密度波形成と粘性過安定の形成を再現することができた。しかしながら、一部の結果は理論と整合するものの、理論的には予想されていない自己重力による微細構造が観察され、詳細な検証が必要であることが明らかになった。この詳細な検証を行うためには、さらなる大規模かつ長時間のシミュレーションが不可欠であることが判明した。そのため、研究計画を一度変更し、シミュレーションのさらなる高速化を実現するためのアルゴリズムとコード開発に注力することとなった。この計画変更により、研究全体の進捗がやや遅れることとなった。昨年度は、開発したアルゴリズムの正当性を理論的に検証し、概ねその作業を完了することができた。現在は、その成果をまとめた論文を執筆中である。それと並行して、開発したコードを用いて当初予定していたシミュレーションを行うために、最適なパラメータの探索を進めている。しかしながら、コードの高速化・最適化については、今後の課題として残されている。以上のように、本研究では、当初の計画通りに進捗させることができなかった部分があり、研究全体としてはやや遅れが生じていると言える。しかしながら、この過程で得られた知見や開発されたコードは、今後の研究を前進させるものと期待される。残された課題に着実に取り組み、研究の進捗を加速させていきたい。
Strategy for Future Research Activity	シミュレーションコードの開発においては、一昨年度より衝突アルゴリズムの改良に取り組んでおり、高速化を目指している。昨年度は、新しいアルゴリズムの理論的な検証をほぼ終え、現在論文を執筆中である。今後は、様々なテスト問題でコードの検証を継続的に行い、実際のシミュレーション研究に適用できるよう、検証とパラメータ探索を進めていく予定である。さらなる長時間計算を可能にするために、重力計算や衝突計算の高速化が課題となっている。このため、今年度は高性能な計算サーバを導入し、コード開発を加速させる計画である。 2021年度では、密度波の基本的なシミュレーションを実施し、線形摂動解析などの理論研究との比較を行った。これはテスト計算であったため、粒子数は実際よりもかなり少なく設定していた。今年度は、改良されたシミュレーションコードを用いて、より現実に近いパラメータ設定での大規模なシミュレーションを実施することを目標としている。これにより、理論や観測との整合性を確認し、研究成果としてまとめていく予定である。今後の研究推進においては、以下の点に重点を置く。 (1) 新しいアルゴリズムの理論的検証を完了し、論文化を進める。また、様々なテスト問題でコードの検証を継続的に行い、実際のシミュレーション研究への適用を目指す。(2) 重力計算や衝突計算の高速化により、より長時間のシミュレーションを可能にする。このため、高性能な計算サーバを導入し、コード開発を加速させる。(3) 改良されたシミュレーションコードを用いて、より現実的なパラメータ設定での大規模なシミュレーションを実施する。これにより、理論や観測との整合性を確認し、研究成果としてまとめる。これらの取り組みを通じて、シミュレーションコードの性能向上と、より現実に即した大規模なシミュレーションの実現を目指す。
Causes of Carryover	当初の研究計画では、計算サーバおよび解析サーバを年度内に導入する予定であったが、昨年度は研究代表者が学内業務が多忙となる諸事情があり、研究に必要不可欠なコード開発が遅延したため、予定していた計算サーバの導入には至らなかった。その結果、研究費の一部が未使用となり、次年度使用額が生じることとなった。昨年度の研究では、当初の計画通りにはいかなかったものの、アルゴリズム開発において一定の進捗があった。このアルゴリズムを用いて大規模なシミュレーションを行うためには、大容量のメモリや高速なCPUが利用可能な高性能な計算サーバが必要不可欠である。そこで、次年度の研究では、前年度の未使用分の研究費を活用し、当初計画していた計算サーバの導入を実施する予定である。この計算サーバを用いることで、開発したアルゴリズムの性能評価や、大規模データを用いた解析を効率的に行うことが可能となる。これにより、研究のさらなる進展と、当初の研究目的の達成が見込まれる。以上のように、次年度使用額が生じた理由は、昨年度の研究の遅延によるものであるが、その分の研究費を有効活用し、研究の進展を図るための使用計画を立てている。