How to scale down solar system data to simulatable values

Question

Okay, I am seriously ashamed for asking this especialy when I geniualy study on physics but there is something that bugs me with the simulation I am working on.

I am re-creating the solar system in an n-body simulation that i programmed before. And I've a problem with scaling down the solar data. So even if I use kg and km for the metric units, the values are far bigger than the variables can hold in programming. Also as some of you know, bigger the value is, bigger the floating point error it makes. (error noise in data) Also it makes it slower to process.

I decided to scale down the data with a reference point, and for that, I took the earth's radius as 1 unit. And scaled down every other distance and radius according to it. (So a unit is 6371 km just to be clear)

But I am not sure whether if I should scale down the mass or not. My common sense says that I should scale down the mass so density of each body should remain same. So I took the density, and calculated a new mass value for each body, with the new scaled down radius. But I am somehow couldn't conviced myself about If it's true or not. So here I am, asking to you :) Should I also scale down the mass?

PS.1: I used using $F = GMm/r^2$ equation for the calculactions as usual. (Iterating it through each body pairs)

If there are other programmers like me interested in making a simulation like this, how did you accomplish this data size problem? Are there any better solutions than scaling down the values?

PS. I have created an excel file that does the scale conversion. So I am sharing the sheet in OneDrive. (http://1drv.ms/1NIekGo) If you can check my calculations and values, that I'd also be really helpful to me. Thanks for any help.

Answer 1

As barrycarter mentioned in his comment, you should be more concerned with units and less with scale.

Generally, it's best to stick to conventional units that people recognize. (It will keep your head on straight and make it easier for others to review your work.) In astronomy, these are a little different than the standard SI units, since things are—shall we say—big, and numbers get out of hand quickly (as you have noted). Here are some units suggested in the Astronomical system of units Wikipedia page:

• Time: Day. This is probably not too important, so feel free to use seconds here. (After all, it's only a 5 order of magnitude difference.)
• Mass: There are several conventions, but the one I've personally seen most is solar masses. Since you're modeling our solar system, it makes a convenient reference.
• Length: I would definitely use astronomical units. (What an appropriate name!)

If you stick to these units, you should be able to limit error due to rounding and reduce the computational load of large numbers.

Answer 2

Consider a nearly circular orbit. On average, $ v = 2 \pi r / T $, and the gravitational force balances the centrifugal force: $$ {G M m \over r^2} = {m v^2 \over r} = {4 \pi^2 m r \over T^2} $$ Solve for $G$ and substitute values for Earth's orbit around the Sun: $$ G = {4 \pi^2 r^3 \over M T^2 } = 4 \pi^2 {\mathrm{AU}^3 \over M_{\odot} \ \mathrm{y}^2} $$ If you want day units instead, use a conversion factor to cancel the years: $$ G = 39.5 {\mathrm{AU}^3 \over M_{\odot} \ \mathrm{y}^2} \left( {1 \ \mathrm{y} \over 365.25 \ \mathrm{d}} \right)^2 = 2.96 \times 10^{-4} {\mathrm{AU}^3 \over M_{\odot} \ \mathrm{d}^2} $$ The factor-label method also works if you start with the MKS value: $$ G = 6.67 \times 10^{-11} \mathrm{m^3 \over kg \ s^2} \left( {1 \ \mathrm{AU} \over 1.496 \times 10^{11} \ \mathrm{m}} \right)^3 \left( {1.99 \times 10^{30} \ \mathrm{kg} \over 1 \ M_{\odot}} \right) \left( {86400 \ \mathrm{s} \over 1 \ \mathrm{d}} \right)^2 = {2.96 \times 10^{-4} {\mathrm{AU}^3 \over M_{\odot} \ \mathrm{d}^2}} $$ Validation tests will tell whether your simulation works correctly. Set up a variety of simple systems where you have a clear idea of what should happen, and then check what actually does happen.

Answer 3

I programmed one by myself some time ago. I used full SI units in combination with doubles (64-bit floating numbers). They work great for the scale of our solar system and are still extremly precise.

        var sun = new Star();
        sun.Position = new Vector3D(0,0,0);
        sun.Velocity = new Vector3D(0,0,0);
        sun.Mass = 1.998855e30;

        var earth = new Planet();
        earth.Mass = 5.9722e24;
        earth.Position = new Vector3D(0, 149.6e9, 0);
        earth.Velocity = new Vector3D(29780,0,0);

        system.World.Objects.Add(sun);
        system.World.Objects.Add(earth);

If you are intersted, my code for this simulation is on github: https://github.com/RononDex/Simulation

The simulation itself is in the subfolder Simulation.Testing