Notebook

これは日々の作業を通して学んだことや毎日の生活で気づいたことをを記録しておく備忘録である。

HTML ファイル生成日時: 2025/11/03 14:57:45.082 (台灣標準時)

Rebound を使い軌道積分する際の注意点

Rebound を使って軌道積分をする際に注意すべき点があると気づいたでござる。

太陽系の三次元的な構造を視覚的に理解するために、太陽、八個の惑星、冥王星、それから数万個の小惑星の軌道積分を Rebound を使って行い、その軌道積分の結果を Matplotlib の 3D plot を使って可視化し、アニメーションを作成してみたでござる。その際、質量ゼロとして軌道積分する小惑星の数が数千くらいのときと、数万にしたときで、軌道積分にかかる時間が大きく異なることに気がついたでござる。小惑星の数が 5000 のときと、小惑星の数が 50000 のときで、計算にかかる時間が 10 倍ではなく、それよりももっと長いようでござる。 sim.add (m=0.0, a=...) としているので、 massless な particle として扱われていると思っていたのでござるが、他の particle に与える重力の影響が計算されているようでござる。それで、計算量が膨大になっていたようでござる。

Rebound の API Documentation を読むと、 simulation object に対して、 .N_active という量を設定できて、他の天体に重力を与える天体の数を決めることができるようでござる、ということがわかったでござる。つまり、太陽、八個の惑星、冥王星についてだけ、他の天体に重力を及ぼす、ということにして、小惑星はこの十個の天体から重力を受けるけれども、小惑星は太陽、惑星、冥王星に重力を及ぼさないし、小惑星同士も重力を及ぼし合わないということにしたい場合には、 sim.N_active = 10 とすればよいようでござる。つまり、こうすると、太陽、惑星、冥王星の十個は、他の全ての天体に重力を及ぼし、それ以外の小惑星は十個の天体から重力を受けるけれども、自分以外の天体に重力を及ぼさない、ということになるようでござる。

太陽、惑星、冥王星は他の全ての天体に重力を及ぼし、小惑星は太陽、惑星、冥王星に重力を及ぼすことにするが、小惑星同士に働く重力は考慮しない、というように設定することもできるようでござる。そのためには、 sim.N_active = 10 に加えて、 sim.testparticle_type = 1 とすればよいようでござる。

デフォルトの設定だと、以下のようになっているようでござる。

parameter	value
sim.N_active	0
sim.testparticle_type	0

Variables の Particles についての部分に関連する記述があったでござる。

やはり documentation をきちんと読み込まないといけぬ、ということを改めて実感したでござる。

About this article:

author: daisuke
file: 20231221_00.html
category: Science___Astronomy
title: Rebound を使い軌道積分する際の注意点
mode: public
last modified: 2023/12/21 09:31:48 (Taiwan Standard Time)
html generated: 2025/11/03 14:57:45.082 (Taiwan Standard Time)

Frequently accessed files

Misc___Taiwan/20240207_00.html

1229 page views
title: 台灣から台灣の外に EMS で荷物を発送する方法

Misc___Taiwan/20240819_00.html

965 page views
title: 住所から台灣の郵便番号を調べる方法

Computer___TeX/20231107_00.html

793 page views
title: LaTeX での counter の利用方法について

Computer___NetBSD/20230119_00.html

596 page views
title: NetBSD でバイナリーパッケージを利用する方法

Computer___TeX/20230726_01.html

567 page views
title: Beamer の Metropolis テーマで block の色を変える方法

Misc___Japan/20240610_00.html

555 page views
title: NHK ラジオの「聞き逃し」の変更点 (2024 年 06 月)

Computer___NetBSD/20250301_01.html

533 page views
title: yt-dlp で YouTube の動画をダウンロードするときのこと

Computer___Network/20240130_00.html

501 page views
title: Google Colaboratory で Python 3.12 を使う方法

Computer___Network/20230516_00.html

481 page views
title: OpenVPN 2.6 を使い VPN Gate に接続するときの注意点

Misc___Taiwan/20240903_01.html

476 page views
title: 台湾の郵便局で EMS を利用して荷物を海外に送る方法

Computer___TeX/20240411_00.html

476 page views
title: LuaTeX を使って PDF ファイルを作成する方法

Computer___NetBSD/20240805_03.html

419 page views
title: NetBSD 10 のインストール

Computer___FreeBSD/20220621_0.html

400 page views
title: FreeBSD での X.org の設定の仕方

Computer___Hardware/20240820_00.html

371 page views
title: Raspberry Pi 5 の演算性能

Computer___NetBSD/20240810_00.html

361 page views
title: Raspberry Pi 5 に NetBSD-10 をインストールしてみた

Misc___Taiwan/20250728_01.html

359 page views
title: 「我那小小多山的國家」という言葉が流行っている

Computer___Debian/20230102_00.html

350 page views
title: Debian GNU/Linux で KVM を使い仮想機械を動かす手順の記録

Computer___NetBSD/20241102_00.html

344 page views
title: ImageMagick の convert コマンドについて

Computer___TeX/20240414_00.html

342 page views
title: LuaTeX での fontspec を使ったフォントの設定の方法

Computer___NetBSD/20250409_00.html

341 page views
title: Raspberry Pi 5 に NetBSD-10.1 をインストールする手順

Computer___Python/20250330_00.html

332 page views
title: Jupyter Notebook で Matplotlib を使うときのこと

Computer___TeX/20240410_00.html

331 page views
title: pdfTeX を使って PDF ファイルを作る方法

Food___Taiwan/20230615_01.html

331 page views
title: 「郭記牛肉食堂」の豬肉炒飯と紅燒牛肉湯 (2023 年 01 月下旬)

Book___Chinese/20240424_00.html

328 page views
title: 曹操の詩

Computer___NetBSD/20250307_00.html

326 page views
title: uim と mozc による日本語入力の設定について

Computer___TeX/20240414_01.html

324 page views
title: LuaTeX を使って日本語と中国語の漢字が混在した PDF ファイルを作成する方法

Computer___NetBSD/20220818_1.html

324 page views
title: Emacs の markdown-mode について

Computer___Network/20230726_00.html

322 page views
title: git の SSL certificate problem の解決方法

Computer___TeX/20230503_00.html

322 page views
title: LaTeX CJK で日本語や中国語を取り扱うための準備について

Computer___NetBSD/20240804_00.html

317 page views
title: NetBSD で NVMM を使って仮想機械を動かす方法

HTML file generated by Kinoshita Daisuke.