Interesting animations by Yukterez, unfortunately it's a bit beyond me to calculate relativistic orbits at the moment
Below is my 2nd attempt at a simulation using more stars
I probably should at least put a few more equations into make things a bit easier. Plummer's density is given by:
\rho(r) = \frac{3M}{4\pi R^3}\left( 1+ \left(\frac{r}{R} \right)^2 \right)^{-\frac{5}{2}}
where
M = total mass of star cluster
R = magic scaling parameter
Now let X_{n} be a...
Hi all,
I refer to the following pdf document, in particular the appendix:
http://adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=1974A%26A...37..183A&link_type=ARTICLE&db_key=AST&high=
End goal is to distribute N stars each having mass m (looks like equal mass stars is the easiest scenario)...
Thanks for the link Loren I will check it out and yes would be to compare. I am going to try and create a more 3 dimensional perspective like Janus's and get the camera zoom around with some smoke/fog effects if I can ever figure out how to get it working. The smoke/fog is turning out to be a...
After great difficulty I managed to get something uploaded to youtube.
Below is my cubes/stars orbiting the black hole
Not quite as good as Janus's but I plan to improve in subsequent simulations
Ok I think I found the problem. Looks like there were 2 issues, the first one was user error on my behalf. That x-value discrepancy we were having was because I accidentally gave you the value for S14 rather than S13. The other issue was in the jungle of + and - signs I got one of them around...
Nice camera work there Janus - you did well to extract data for 10 stars.
tony873004 I take it the calculation for E is still cactus in your code01.html calculator? I thought my E=2.11 was pretty close
I got those instructions from:
https://downloads.rene-schwarz.com/download/M001-Keplerian_Orbit_Elements_to_Cartesian_State_Vectors.pdf
Thanks for the new calculator I will play with the numbers and try to see what's going on.
For the new data can I just set M=0 for all stars?
I just dug up an old dusty book from the garage and found it gave the following for S1
\mathit{a} = 20.5 \text{ mpc} \\
\mathit{e} = 0.496 \\
\mathit{i} = 120.82 \text{ deg} \\
\Omega = 341.61 \text{ deg} \\
\omega = 115.3 \text{ deg} \\
P = 132 \text{ years}
The trouble with this data...
I was hoping it would be a simple fix of adding .0 to the integer values given in the table for S13 as it is the only star with integer values but doesn't look to be the case unfortunately, especially since all other stars are producing orbits.
I am using 3.61 \times 10^6 M_{\odot} \approx...
Unless I have miscalculated somewhere, the S13 star doesn't look to be in a relatively stable orbit. Based on that data it will fall into the black hole (in my simulation it gets flung into infinity after a few seconds). Noted that the uncertainty in the measurements for that particular star...
Thanks for the detailed information, you probably saved me a few weeks by mentioning that semi-major axis trap.
What I did was use the above info to calculate the initial position and velocity of the star S1 and then used Euler's method to calculate subsequent positions using 3.6 million solar...
That's beautiful thanks tony873004. If I understand correctly, using the information below for the S1 star
we have that
M_{0} = 360 - \frac{360}{T} (t_{0}-P)
using D=P
which means
M_{0} \approx 360 - \frac{360}{94.1}(2002.6 - 2016.25)
Does that sound reasonable?
Sounds great if you could find more complete and current data.
Yes it does look like I need help converting to Cartesian coordinates. My main trouble is with the mean anomaly which isn't provided in the table.
We have a\text{, } e\text{, } \omega\text{, } \Omega\text{, } i but no M_{0} =...
Excellent thanks tony873004 that animation should it extremely easy for me to do the animation since it has the nice Cartesian coordinates. Were you able to do that animation using the orbital parameters in this thread or was it already done?
Thanks Janus that table is looking better. Looks...