Kleines Rechenprogramm zum Lichtsegler
02.01.2013 um 01:16Anzeige
Z. schrieb:Weil wir nicht gerade.. beste Freunde sind...etc.. und ich keine Lust habe mich wieder zu ärgern!Das Video ist doch nicht von dir, wo also ist das Problem?
Z. schrieb:Das Problem ist nicht das Video. Videos bereiten keine Probleme.Du postest den Link nicht, weil ich dir deswegen Probleme bereiten könnte? Was ist so besonders an diesem Link? Versteh die Geheimnistuerei nicht.
"Menschen" bereiten Probleme.
Z. schrieb:anstatt mein Frage zu respektieren und mir eine ehrliche Begründung oder die Motivation zu nennen.Meine Güte Z. Du hast gesagt die Simulation ist nicht von dir. Ist in Foren doch üblich nach einer Quellenangabe zu fragen. Würde die anderen hier vielleicht auch interessieren.
// Syntax für Zeit t und Geschwindigkeit v
// als Funktion von Strecke R beim Sonnensegel.
// Code & Plot by Yukterez
// Variablen
r1 := 15*AU:
A := (1e10*m)^2:
pA := 0.05*kg/m^2:
MX := 2000*Msol:
Msail := 100*kg+A*pA:
// Parameter
E0 := 1368*kg/sek^3:
EX := E0*(MX/Msol)^3.5*(AU/r1)^2:
Msol := 2e30*kg:
G := 6.67384e-11*m^3/kg/sek^3:
c := 299792458*m/sek:
AU := 1.5e11*m:
z := sqrt((1+v/c)/(1-v/c))-1:
F := EX*A/c/(1+z):
g := MX*G/r1^2:
a := F/Msail-g:
r2 := r1+R:
// Evaluierung
ai := int(a/(r/r1)^2, r=r1..r2)/R:
t := sqrt(2)*sqrt(R)/sqrt(ai):
// Plot
kg:=1: m:=1: sek:=1:
V := plot::Implicit2d(-v+ai*t,
R=0..5*AU, v=0..c): plot(V);
// Startfähigkeit
aS := EX*A/c/Msail-g:
OK = sign(int(aS/(r/r1)^2, r=(r1+1e-2)..(r1+1e+2))/1e2);
reset()
// ____________________________________________________________
// Syntax für Zeit t und Geschwindigkeit v
// als Funktion von Parameter x und Strecke R.
// Code & Rechnung by Yukterez
// ANIMATIONSPARAMETER: (x)
// Variablen
r1 := 15*AU:
A := (1e10*m)^2:
pA := 0.05*kg/m^2/(x):
MX := 400*Msol*(x):
Msail := 100*kg+A*pA:
// Parameter
E0 := 1368*kg/sek^3:
EX := E0*(MX/Msol)^3.5*(AU/r1)^2:
Msol := 2e30*kg:
G := 6.67384e-11*m^3/kg/sek^3:
c := 299792458*m/sek:
AU := 1.5e11*m:
z := sqrt((1+v/c)/(1-v/c))-1:
F := EX*A/c/(1+z):
g := MX*G/r1^2:
a := F/Msail-g:
r2 := r1+R:
// Evaluierung
ai := int(a/(r/r1)^2, r=r1..r2)/R:
t := sqrt(2)*sqrt(R)/sqrt(ai):
// Plot
kg:=1: m:=1: sek:=1:
V := plot::Implicit2d(-v+ai*t,
R=0..5*AU, v=0..1.1*c, x=1..5, Frames=200): plot(V);
reset()
Ein Quadratkilometer wiegt somit 757 g.Da ich, was natürlich Fehlerhaft sein kann?, bei einer Entfernung zur Sonne von 3 Sr, und entsprechenden, der Berechnung zu entnehmenden Vorgaben, Gewicht, Fläche, Effektivität etc, auf ca. 4943 ms anfängliche Beschleunigung komme, habe ich mir erlaubt dies zunächst als "stetig" anzunehmen, um Zeit und Strecke zu umreissen, bis 99% c anliegen. Wie gesagt Borduhr.
Besonders aber macht´s mir einen Riesenspass, mal wieder die theoretischen Antriebe zu Sichten und zu berechnen. Im Moment bin ich gerade nochmal am prüfen wie schnell ein Lichtsegel im Orbit um die Sonne theoretisch werden kann. ......Und es geht auch nicht um v > c.
Das ist a priori mein Postulat zum Thema Sonnensegeln!
Bisher werden von der ESA nur 25% von c angenommen und selbst wenn wir hier nicht auf nahe c kommen können, wird dieser Versuch hier zeigen können, das es schneller geht. ;) Derer festen Überzeugung ist er, seine/r eine/r. Thats it.
Z. schrieb:Von daher denke ich, das es sich bei Resultat um eine einigermassen korrekte Näherung hält.Darf ich noch mit einer Antwort rechnen?
Wie ist Ihr Ergebnis, mit Ihrem gezeigtem Rechenweg?
yukterez schrieb:Selbst wenn wir das gelten liessen wäre die Formel dafür eine andereAuf welches Ergebnis Sie (ms Segel) mit Ihrer dieser Formel bekommen wenn Sie die gleichen Daten zugrunde legen? Die Daten haben Sie ja.