I've managed to find a way to correct for the obleateness, and the results are pretty remarkable with regards to to Gallian moons, The absolute error on Europa's position has been reduced to about 32000 km, a vast improvement.
However I am at a loss what is causing the huge errors on Saturn...
Yes I've noticed Saturns error is very high which I can't explain.
The numbers I posted are the largest deviations that occur in a 3 year period. Europa deviates by about 100 degrees while with Io, it's deviation is >180 degrees, it it remains it the same orbit as the reference.
I get your...
I don't experience any significant changes in accuracy when changing using 15 steps/hour compared to 10 steps/hour. And as most bodies are simulated quite accurately, I don't expect it to be a issue with the integrator. I am using Velocity Verlet,
Below is a exhaustive list of the bodies that...
Thanks, can you point me to some resources on how to do that?
I am making a solar system simulator, and the errors I get when simulating the larger moons of jupiter are about the diameter of their orbits when simulating three years. I have made some animations and can see that they fall behind...
I have been working on implementing a solar system system simulator in C++ - but am getting incorrect results on the order of 10^10 km, and it seems that the planets are moving directly away from the sun. I suspect that there is a mistake in the integrator (Velocity Verlet) that I have posted...
I am working on a simulation of the solar system, but I am having problems with the orbits with Europa and Io shown in the image below, the orbits of Jupiter, Ganymede and Callisto are consistent with NASA's simulations. For all bodies I have used identical parameters in generating the initial...
Thank you all for your suggestions. @Mark44 Is correct that I am used to Python with only a passing familiarity with C++.
The project I am working on is N-Body simulator for predicting Solar System orbits in Python, and I would like to speed
up the computations without loosing the graphical...
Is there a computer friendly list with the known bodies in the Solar System astroid belt, that contains their physical properties, position and velocities?
Errh.. I had accidentally chosen a another coordinate center than the one I intended. It simulates beautifully now!
So thank you very much, all of you, for your help.
Hi,
Below are the starting values, each line is one body and the order is
Sun
Mercury
Venus
Earth
Mars
Jupiter
Saturn
Uranus
Neptune
Pluto
The format is x, y, z, vx, vy, vz and the units are meters and meters/second.
X, Y, Z...
Hi all.
Before I tried Verlet, I had written a RK4 integrator where I got the same results.
I have rewritten the program and (still) get the same results. Would it be helpful for me to post the code for the new integrator (still Verlet) here, or should I post it on the programming forum?
The initial values for the simulations, x, y, z position and velocities in all three dimensions, are the first entries in the two diagnostic data sets. The simulator does not estimate the velocities, even though the diagnostic data set does include velocities. The two diagnostic data sets where...
The initial values are from the same dataset as the diagnostics. The main loop is given below:
def gen_coords():
if not os.path.exists('trajectories/'):
os.makedirs('trajectories/')
# Generate coordinates and save to file
# NumPy arrays for coordinates
trajectory_earth =...
Hi, as a side project I am making a solar system simulator in Python, but I am getting extreme inaccuracies.
As the image shows - I am comparing the simulation to data obtained from NASA's Horizon.
I have written a a basic Verlet integrator in the Python class posten below:
import numpy as...
I am currently working on problem 32-15 in Calculus by Spivak, and in question (b)
in the bottom line there is a relation
[\phi_1'(b)\phi_2(b) -\phi_1'(a)\phi_2(a)]+[\phi_1(b)\phi_2'(b)-\phi_1(a)\phi_2'(a)]>0
But I can only get it to work out if
[\phi_1'(b)\phi_2(b)...
I am reading §22 of Topology by Munkres, in Theorem 22.2 the function g is said to be constant on each set p^(-1)({y}). However the only explicit property given in Corollary 22.3 to the function g is that it is continuous and surjective, but Theorem 22.2 to g in the proof. Is it implied that g...
The meaning of "different" in Munkres' Topology
Hi, I'm working on problem 20.8(b) (page 127f) in Munkres' "Topology", the problem is to show that four topologies are "different". Does different in this context mean that they are unequal - in which case one can contain the other, or...
Hello, I am studying mathematics and physics on my own. However, it is the mathematical side of physics that intrest me.
I would like some recommendations on which books would be beneficial for me to read, to gain an understanding of Special and General Relativity, and ultimatly string...