Timeline for Mass distribution in NDSolve
Current License: CC BY-SA 3.0
16 events
| when toggle format | what | by | license | comment | |
|---|---|---|---|---|---|
| Jul 24, 2016 at 16:16 | vote | accept | Feyre | ||
| Jul 19, 2016 at 16:35 | comment | added | Feyre |
@march Oh, I see, well that shouldn't be too hard then. I was actually already thinking about converting to AU, but if this is a standard measure I'll use that.
|
|
| Jul 19, 2016 at 16:34 | answer | added | Siav Josep | timeline score: 6 | |
| Jul 19, 2016 at 16:30 | comment | added | march | Well, scaling quantities in a reasonable set of units is a good skill to have in any situation, really. You want as many of the parameters to be on the order of 1 as possible. | |
| Jul 19, 2016 at 16:29 | comment | added | Feyre | @march I'll look into it, though it may be above my pay-grade, I'm not an astronomer, doing this for fun. | |
| Jul 19, 2016 at 15:48 | comment | added | Feyre | @m_goldberg I meant that the result is accurate to that many digits for Earth and Mars. I'm using higher digit precision in the actual calculation. I'm still refining the model, the asteroid belt is one avenue of refinement. | |
| Jul 19, 2016 at 15:48 | comment | added | Feyre |
@KraZug Well I just input the result of a generalised calculation. In the actual notebook I have a formula which calculates the acceleration vectors for n bodies.
|
|
| Jul 19, 2016 at 15:22 | comment | added | march | I would recommend working in something like N-body units. Keeping around all those factors of ten can be a problem, precision-wise, I think. | |
| Jul 19, 2016 at 15:10 | history | tweeted | twitter.com/StackMma/status/755419506172915712 | ||
| Jul 19, 2016 at 14:58 | comment | added | Kuba | related by the topic: Only final result from NDSolve | |
| Jul 19, 2016 at 14:52 | comment | added | m_goldberg | In a computation giving only 3-4 significant digits over a year, I strongly suspect the asteroid belt objects have no detectable effect on the major bodies of the solar system. | |
| Jul 19, 2016 at 14:49 | comment | added | SPPearce | You could also label the vectors with an index, e.g. X[1][t_] or X["sun"][t_] etc, and that lets you iterate over those. | |
| Jul 19, 2016 at 14:48 | history | edited | m_goldberg | CC BY-SA 3.0 |
Minor clean-up
|
| Jul 19, 2016 at 14:47 | comment | added | SPPearce | While it doesn't answer your mass question, you really should be in vector form here: X1[t_] = {x1[t], y1[t], z1[t]};X2[t_] = {x2[t], y2[t], z2[t]}; eqns = {X1''[t] == (g m2 (X2[t] - X1[t]))/Norm[(X2[t] - X1[t])]^3, X2''[t] == (g m1 (X1[t] - X2[t]))/Norm[(X2[t] - X1[t])]^3}; s = NDSolve[{eqns, vel, pos}, {X1[t], X2[t]}, {t, 0, tmax}]; | |
| Jul 19, 2016 at 14:28 | history | edited | J. M.'s missing motivation♦ |
edited tags
|
|
| Jul 19, 2016 at 14:11 | history | asked | Feyre | CC BY-SA 3.0 |