Forum der Vereinigung der Sternfreunde

Liebe Sternfreunde, im aktuellen Heft habe ich vorgerechnet, welchen Weg ein Photon nimmt, das tangential an einem Schwarzen Loch vorübergeht. Der minimale Abstand ist 1,5 Schwarzschildradien (r0=1.5). Weiter innen liegende Bahnen führen in das Schwarze Loch hinein. Ihr könnt also gern mit verschiedenen r0 experimentieren. Ich zumindest fand das sehr anschaulich. Die Abbildung zeigt r0=1,505 – das Photon kreuzt die eigene Bahn zwei Mal. Ganz nahe an 1,5 heran funktioniter das Priogramm übrigens nicht mehr, für diesen Grenzfall ist es nicht ausgelegt.

Ich werde übrigens zu diesem Thema auf der virtuell stattfindenden Würzburger Frühjahrstagung sprechen, am kommenden Samstag 13 Uhr bis 13:30 Uhr. Der Einwahllink liegt auf der Webseite unseres Vereins.

Hier der Programmcode aus dem Heft, auf der Webseite und auch komplett:

from turtle import *
from random import *
from math import *

def plot(x,y): # einen Punkt setzen
    penup()
    goto(50*x,50*y)
    pendown()
    dot(2)
    hideturtle()

def mycircle():
    penup()
    goto (0,-50)
    begin_fill()
    circle(50)
    end_fill()

        
# Hauptprogramm
tracer(0,0) # Turtle aus
mycircle()
r0=float(input("r0: "))
r=1.00000001*r0
phi=0
Delta=1e-3
while r<10:
    for i in range (10):
        drdphi=r*r*r*r/(r0*r0*r0)
        drdphi=drdphi*(r0-1)-r*(r-1)
        if (drdphi<0):
            print("Fehler: r0 ist ungueltig")
        drdphi=sqrt(drdphi)
        r=r+Delta*drdphi
        phi=phi+Delta
        x=r*cos(phi)
        y=r*sin(phi)            
        plot(x,y)
        # auch nach unten zeichnen
        x=r*cos(-phi)
        y=r*sin(-phi)            
        plot(x,y)
    update()
z=input("Enter")

_________________
Uwe Pilz, Fachgruppen Kometen und Astrophysik/Algorithmen.
Oft benutzte Instrumente: Fujinon 16x70 FMT-SX-2, TMB Apo 105/650, Skywatcher Evostar 120/900 ED, Ninja Dobson 320/1440, Smartscope Dwarf-3