Forum der Vereinigung der Sternfreunde

Forum of the German Amateur Astronomy Association
Aktuelle Zeit: 14. Mai 2021, 08:47:58 AM

Alle Zeiten sind UTC+02:00




Ein neues Thema erstellen  Auf das Thema antworten  [ 1 Beitrag ] 
Autor Nachricht
BeitragVerfasst: 20. April 2021, 15:04:44 PM 
Offline
Meister
Benutzeravatar

Registriert: 20. Januar 2013, 20:03:54 PM
Beiträge: 1110
Wohnort: Leipzig
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:
Code:
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")


Dateianhänge:
1505.jpg
1505.jpg [ 15.05 KiB | 564 mal betrachtet ]

_________________
Uwe Pilz, Fachgruppen Kometen und Astrophysik/Algorithmen.
Oft benutzte Instrumente: Swarovski SLC 7x50B, Fujinon 16x70 FMT-SX-2, TMB Apo 105/650, Skywatcher Evostar 120/900 ED, Ninja Dobson 320/1440
Nach oben
   
Beiträge der letzten Zeit anzeigen:  Sortiere nach  
Ein neues Thema erstellen  Auf das Thema antworten  [ 1 Beitrag ] 

Alle Zeiten sind UTC+02:00


Wer ist online?

Mitglieder in diesem Forum: 0 Mitglieder und 2 Gäste


Du darfst keine neuen Themen in diesem Forum erstellen.
Du darfst keine Antworten zu Themen in diesem Forum erstellen.
Du darfst deine Beiträge in diesem Forum nicht ändern.
Du darfst deine Beiträge in diesem Forum nicht löschen.
Du darfst keine Dateianhänge in diesem Forum erstellen.

Suche nach:
Gehe zu:  
Powered by phpBB® Forum Software © phpBB Limited
Deutsche Übersetzung durch phpBB.de