Help in runge kutta problem to determine period of a comet

Click For Summary
SUMMARY

The discussion focuses on using the Runge-Kutta methods, specifically order 2 and 4, to determine the orbital period of Halley's comet, with perihelion at 0.586 AU and aphelion at 35.1 AU. The period is calculated using Kepler's third law, where the semimajor axis is derived from the average of perihelion and aphelion distances. Participants suggest setting up a system of linear differential equations based on Newton's laws and the inverse square law of gravitation, followed by programming these equations to simulate the comet's trajectory and calculate the orbital period.

PREREQUISITES
  • Understanding of Runge-Kutta methods, specifically order 2 and 4
  • Familiarity with Kepler's laws of planetary motion
  • Knowledge of Newton's laws of motion and gravitational theory
  • Basic programming skills for numerical simulations
NEXT STEPS
  • Implement Runge-Kutta methods in a programming language of choice to solve differential equations
  • Study Kepler's laws in detail, focusing on their application in celestial mechanics
  • Explore numerical methods for solving ordinary differential equations (ODEs)
  • Research the principles of conservation of energy and angular momentum in orbital mechanics
USEFUL FOR

Astronomy students, physicists, and engineers interested in celestial mechanics, numerical methods, and the dynamics of cometary orbits will benefit from this discussion.

Gvilla1000
Messages
1
Reaction score
0
Hi everyone:

I have trouble in an assignment I was given. The instructions are to utilize Runge-Kutta order 2 and 4 to determine the period of the comet Halley (perihelion = 0.586 AU and aphelion = 35.1 AU).

The period of the orbit is given by Kepler's third law:
Period = Sm3/2

where Sm is the Semimajor axis acquired from:
Sm= (P+A)/2

where P is the perihelion and A is the aphelion

My problem is...if this is such a straightforward answer...how can I use Runge-Kutta in a problem like this? what would be my function that I would need to use Runke Kutta on?
 
Technology news on Phys.org
You could use Newton's laws of motion and the inverse square law of gravitiation. Or you could set up some equations using conservation of energy and angular momentum. You will also need the initial conditions, e.g. the speed of the planet at periohelion that gives you the correct aphelion distance.

It's hard to say what is the "best" way to set this up without knowing the bigger picture of how this course fits into your studies. A computer scientist who knows about numerical methods but not much physics, a mech engineer who knows about dynamics but not much astronomy, and an astronomer, would probably do it three different ways.
 
Last edited:
If I had given this problem to my students, I would have meant it as follows:

1. Set up a system of linear differential equations to describe the motion of the comet in the Sun's gravitational field, starting with F = ma and Newton's law of gravitation.

2. Write a program to solve those equations using the specified methods, to generate the path of the comet given its initial position and velocity.

3. For the initial position, use the known perihelion for Halley's comet. Adjust the initial velocity so that the aphelion is the known value for Halley's comet.

4. Trace the comet through one complete orbit and find how much time it takes.

But your professor might have meant something else, so you had best ask him if this is the sort of thing he is looking for.
 

Similar threads

  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 65 ·
3
Replies
65
Views
9K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 4 ·
Replies
4
Views
8K
  • · Replies 1 ·
Replies
1
Views
3K
Replies
8
Views
15K
  • · Replies 19 ·
Replies
19
Views
19K
Replies
11
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 13 ·
Replies
13
Views
4K