Danke dass du mich daran erinnerst dass dein Zwilling schon wieder eine ganze Seite
hat. Damit sich niemand durch euren Müll blättern muss um etwas zum Thema zu lesen zurück zum Anfang:
Nachdem
die kritische Elite schon seit Monaten keine
Wurfparabel liefern kann bleibt es mal wieder an den Relativisten hängen.
Wir nehmen also eine Erde und komprimieren sie zum schwarzen Loch. Im Abstand von 10 Schwarzschildradien zum Schwerpunkt platzieren wir einen Werfer, der einen Ball mit der neutonischen Orbitalgeschwindigkeit wirft; im 1. Beispiel in einem Winkel von 0°, und im 2. Beispiel mit 45°. Links wird nach Newton geplottet, und rechts mit den selben Startbedingungen die relativistische Wurfparabel nach Einstein. Der innere Kreis repräsentiert den Ereignishorizont der Masse, und der äußere die Schale auf der der Werfer steht.
Beispiel 1 (0°)
Beispiel 2 (45°)
RechnungCode:
(* relativistische Wurfparabel | yukterez.net 2016 | Syntax: Mathematica *)
G = 678*^-13; M = 6*^24; c = 3*^8; rs = 2 G M/c^2;
v0 = Sqrt[G M/r0]; vr0 = v0/Sqrt[2]; r0 = 10 rs; vθ0 = v0/r0/Sqrt[2]; θ0 = 0;
T = 7*^-9; step = T/500;
sol =
NDSolve[{
r''[t] == -((G M)/r[t]^2) + r[t] θ'[t]^2 - (3 G M)/c^2 θ'[t]^2,
r'[0] == vr0,
r[0] == r0,
θ''[t] == -((2 r'[t] θ'[t])/r[t]),
θ'[0] == vθ0,
θ[0] == θ0,
τ'[t] == Sqrt[c^2 r[t] + r[t] r'[t]^2 - c^2 rs + r[t]^3 θ'[t]^2 - r[t]^2 rs θ'[t]^2]/(c Sqrt[r[t] - rs] Sqrt[1 - rs/r[t]]),
τ[0] == 0
}, {r, θ, τ}, {t, 0, T + step},
WorkingPrecision -> 40, MaxSteps -> Infinity];
x[t_] := (Sin[Evaluate[θ[t] /. sol]] Evaluate[r[t] /. sol])[[1]]
y[t_] := (Cos[Evaluate[θ[t] /. sol]] Evaluate[r[t] /. sol])[[1]]
Do[Print[
Rasterize[Show[
Graphics[{Circle[{0, 0}, r0], Circle[{0, 0}, rs]}, Frame -> True, ImageSize -> 400, PlotRange -> 14 rs],
Graphics[{Point[{x[т], y[т]}]}],
ParametricPlot[{x[t], y[t]}, {t, 0, т}, PlotStyle -> {Opacity[0.3], Gray}],
ParametricPlot[{x[t], y[t]}, {t, Max[Re[т - 1*^-10], 0], т}, PlotStyle -> {Red}]]]],
{т, step, T, step}]
Vergleichend,
Anmelden
Registrieren
FAQ
Suche
