So berechnen Sie die Temperatur des planetaren Gleichgewichts in binären Systemen

Ich arbeite an einem Programm, das die grundlegenden Statistiken eines terrestrischen Planeten in Doppelsternsystemen generiert. Ich bin nicht der Beste in dieser Art von Mathematik, also habe ich Probleme, die planetare Gleichgewichtstemperatur zu berechnen.

Die Formel, auf die ich mich gerade beziehe, wurde vor fast 4 Jahren in einem Thread gepostet: Wie man die erwartete Oberflächentemperatur eines Planeten berechnet .

Formel

Es scheint ziemlich genau zu sein, da (wenn ich mich recht erinnere) die Verwendung eine gute Antwort für Erde und Mars liefert. Dies gilt jedoch für ein System mit einem strahlenden Körper und nicht für ein System mit zwei strahlenden Körpern (Doppelsterne).

Mein Instinkt ist, die beiden Sterne als einen großen Stern zu behandeln und die Wattzahl zu addieren ( w / M 2 ) der Planet empfängt von den verschiedenen Sternen. Also Gesamtwatt = Watt Stern 1 + Watt Stern 2. Das erscheint mir vernünftig, zumindest für die Zwecke der Berechnung der Planetentemperaturen.

Ich habe bereits daran gearbeitet, die von den Planeten empfangenen Watt aus den verschiedenen Quellen zu berechnen. Ich habe diese Gleichung verwendet und Werte für beide Sterne eingefügt.

F ( L / S 2 )

  • L = Leuchtkraft (in Bruchteil der Sonne), Sterne 1 & 2 (jeweils): 0.0159&0.7758
  • S = Große Halbachse (in AU), Sterne 1 & 2 (jeweils): 0.1228&3.6000
  • F = Mittlerer Sonnenfluss der Erde ( w / M 2 ) 1362

Die obige Berechnung ergibt also einen mittleren Sonnenfluss für den Beispielplaneten von 1344 w / M 2 von Stern 1 und 92 w / M 2 von Star 2. Wenn ich diese zusammenzähle, erhalte ich eine mittlere globale Wattzahl von 1436 w / M 2 .

Ich bin mir jedoch nicht sicher, wie ich die ursprüngliche Gleichung ändern soll, um den oben berechneten Sonnenfluss zu verwenden (idealerweise ohne auf die ursprünglichen Variablen der Leuchtkraft und der großen Halbachse verweisen zu müssen). Wie ändert man die Formel richtig?

Jede Hilfe wäre sehr willkommen!

Antworten (1)

Die vom Planeten abgestrahlte Leistung unter der Annahme, dass er auf seiner gesamten Oberfläche die gleiche Temperatur hat:

P R A D = ϵ σ T e Q 4 4 π R P l 2

Wo ϵ ist der Emissionsgrad (um der Formel in der Frage zu entsprechen, setzen Sie diese auf 1), σ ist die Stefan-Boltzmann-Konstante , T e Q ist die Gleichgewichtstemperatur, und R P l ist der Radius des Planeten.

Die vom Planeten von den Sternen absorbierte Leistung, vorausgesetzt, er ist ausreichend weit von jedem Stern entfernt, so dass es eine vernachlässigbare Beleuchtung von diesem Stern jenseits von 90 ° vom substellaren Punkt gibt:

P A B S = ich F ich ( 1 A ich ) π R P l 2

Wo F ich ist der Fluss aus der ich Stern, und A ich ist die planetare Albedo in Bezug auf die ich Stern (dies kann unterschiedlich sein, da das Reflexionsvermögen wellenlängenabhängig ist und die Sterntemperatur beeinflusst, welche Wellenlängen emittiert werden).

Setzen Sie also abgestrahlte und absorbierte Leistung gleich und ordnen Sie neu an:

T e Q = [ ich F ich ( 1 A ich ) 4 ϵ σ ] 1 / 4

Die Einschränkung hier ist, dass die Flüsse aufgrund der Entfernungsänderungen des Planeten von den Sternen zeitvariabel sein werden (dasselbe gilt für einen Planeten in einer exzentrischen Umlaufbahn um einen einzelnen Stern), und es ist nicht trivial, herauszufinden, wie das geht mittel das aus. Der reale Planet wird Zeit brauchen, um sich aufzuheizen und abzukühlen, und in dieser Zeit werden sich die Entfernungen und Flüsse ändern.

Nehmen Sie also die Zahlen, die Sie aus dieser Berechnung der "Gleichgewichtstemperatur" erhalten, mit einer Prise Salz.


Wie gewünscht, hier einige Beispielrechnungen. Zunächst einmal der Plausibilitätscheck: Erde. Nimmt man 1361 W/m 2 als Wert der Solarkonstante und verwendet man eine Albedo von 0,3 und einen Emissionsgrad von 1, ergibt sich die Temperatur wie folgt

T e Q = [ 1361   W M 2 × ( 1 0,3 ) 4 × 1 × ( 5.6704 × 10 8   W M 2 K 4 ) ] 1 / 4 255   K

Ein bisschen kühl, aber nicht so schlimm: Der Treibhauseffekt hält die tatsächliche Erde wärmer, was durch die Verwendung eines Emissionsgrads berücksichtigt werden könnte ϵ < 1 .

Nun zu dem Beispiel in Ihrer Frage mit Flüssen von 1344 W/m 2 und 92 W/m 2 . Nehmen wir außerdem an, dass, während der erste Stern sonnenähnlich ist (also werde ich den Wert von 0,3 für die Albedo beibehalten), der zweite Stern kühler als die Sonne ist. Dieser Stern emittiert mehr von seinem Licht im Rot, wo der Planet weniger reflektiert. Ich werde für diesen Stern eine Albedo von 0,25 verwenden.

Die Berechnung funktioniert wie folgt (Weglassen der Einheiten im Berechnungsschritt aus Platzgründen):

T e Q = [ 1344 × ( 1 0,3 ) + 92 × ( 1 0,25 ) 4 × 1 × ( 5.6704 × 10 8 ) ] 1 / 4 258   K

Hoffentlich hilft das!

Danke für diese Antwort. Ich habe ein bisschen Probleme, einen Sinn zu finden, wie man die Gleichung anwendet. Könnten Sie vielleicht ein Beispiel geben, indem Sie Zahlen für zwei Planeten eintragen? Wenn das zu viel verlangt ist, keine Sorge!
@n_bandit - sicher, keine Sorge.
Danke! Dies ist eine großartige Antwort. :) Die einzige verbleibende Frage, die ich habe, betrifft den 0,25-Exponenten. Würde sich dieser Exponent unter allen denkbaren Bedingungen ändern? Hat es zum Beispiel eine Beziehung zur Rotationsgeschwindigkeit des Planeten?
@n_bandit - die Potenz von 1/4 ist auf das Stefan-Boltzmann-Gesetz für Schwarzkörperstrahlung zurückzuführen. Diese geht in die erste Gleichung für die Strahlungsleistung des Planeten ein, die den Term σT⁴ enthält. Eine Änderung dieses Exponenten würde bedeuten, dass ein ganz anderer physikalischer Prozess stattfand (und würde auch einige Änderungen an den beteiligten physikalischen Konstanten erfordern, damit die Einheiten richtig funktionieren).