1 / 13

12/11/16 数値演習 Gravitational calculation in RK 細野七月

12/11/16 数値演習 Gravitational calculation in RK 細野七月. Introduction. ◆地惑で頻出の数値計算応用例. 重力. Gravity. particle. N. two body problem. interaction. Physics. in Physical Unit M=1kg, L=1m, T=1sec => G=6.67e-11. in Planetary Unit M=1Solar Mass, L=1AU, G=1 => T=1.0/2π year. L5.

coen
Download Presentation

12/11/16 数値演習 Gravitational calculation in RK 細野七月

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 12/11/16 数値演習Gravitational calculation in RK細野七月

  2. Introduction ◆地惑で頻出の数値計算応用例 重力

  3. Gravity particle N two body problem interaction

  4. Physics in Physical Unit M=1kg, L=1m, T=1sec => G=6.67e-11 in Planetary Unit M=1Solar Mass, L=1AU, G=1 => T=1.0/2π year

  5. L5 Physical Application L3 L1 L2 ラグランジュポイント 検索 L4

  6. Trojan ISAS/JAXA

  7. Simulation 5.0-5.4 AU Gaussian distribution massless particles Sun 5.2 AU eccentricity=0 9.6 AU eccentricity=0

  8. Simulation z = 0

  9. Memory Structure

  10. Result

  11. Result t = 300 yr

  12. Summary ●RK法で重力計算をして、Trojanの再現を試みた。 ⇒それっぽいものはできた！ 　⇒それっぽくないのもできた… ●どうすればそれっぽくないのを取り除けるか？ ├1)そもそも300年という計算時間は短すぎ。 │　もっと長く計算するべき。 ├2)他の惑星からの摂動も考えるべき。 ├3)離心率や軌道傾斜角を無しにしてしまっている。 │　実際には円軌道では無いのだから考慮するべき。 └4)そもそも彼らはいつできたのか？ 　　現在の太陽系と同じセッティングでいいのか？

  13. % HOW 2 COMPILE %>gfortran runge.f90 xrunge f90 % % HOW 2 PLOT %>gnuplot -e “i = 0” plot.plt % % HOW 2 MAKE %>Makefile %>Makefile is not configured... % % quit % %see you again! :) % Usage

More Related