Register to reply

Orbitals of celestial bodies

by eNathan
Tags: bodies, celestial, orbitals
Share this thread:
eNathan
#1
Mar11-13, 11:19 PM
P: 352
Hi all,

I have a question about calculating the orbitals of bodies of mass in space (newton's basic laws). I am writing a program to simulate the orbitals of bodies in space -- basically, you define the object's mass, size, location, and initial velocities and watch how they interaction via gravitational attraction.

So we know the following equations that govern the motion of these bodies:

Magnitude of force due to gravitational attraction:
[itex]F=G\frac{m1*m2}{r^2}[/itex]
The direction is, at any instant in time, points in the same direction as displacement between the centers of their mass (points in direction of the other body of mass).

We know the magnitude and direction of the force, and from this, we would say that the acceleration due to gravity that object 1 undergoes is given by:
[itex]a=\frac{F}{m1}[/itex]
The direction of which is the same direction as the force vector.

Now suppose I asked the following question: what is the total displacement that object 1 undergoes given any arbitrary [itex]Δt[/itex] (and initial velocity is known)? Currently, my program just changes the velocity vector based on the direction of [itex]a[/itex] at whatever moment the refresh was called, and from that changes the object's position, but I know this is not entirely accurate (e.g., this is comparable to finding the area under a curve by diving it into tiny rectangles), because the direction of [itex]a[/itex] is constantly changing.

What is the calculus method of doing this?

Much appreciated!
Phys.Org News Partner Physics news on Phys.org
Physicists discuss quantum pigeonhole principle
First in-situ images of void collapse in explosives
The first supercomputer simulations of 'spin?orbit' forces between neutrons and protons in an atomic nucleus
A.T.
#2
Mar12-13, 05:51 AM
P: 3,914
Quote Quote by eNathan View Post
What is the calculus method of doing this?
You use this:
http://en.wikipedia.org/wiki/Euler_method
Better are:
http://en.wikipedia.org/wiki/Midpoint_method
http://en.wikipedia.org/wiki/Runge%E...3Kutta_methods

There's lots of code examples online for the specific task you are doing.


Register to reply

Related Discussions
Motion of electrons in orbitals and shape of orbitals Chemistry 5
Looking for specific data about thrust, hydrogen, and celestial bodies Classical Physics 6
Finding the equilibrium bwetween two celestial bodies Astronomy & Astrophysics 2
Distribution of stars and other celestial bodies Astronomy & Astrophysics 0
Constructing orthogonal orbitals from atomic orbitals Atomic, Solid State, Comp. Physics 8