Wie findet man die Dichte eines Planeten und seines Kerns unter Berücksichtigung der Gravitationskompression in ihnen?

Stellen Sie sich einen Planeten vor P des Radius R größer als die Erde, deren Masse ist R 3 . Nach dem Quadratwürfelgesetz hat dieser Planet die gleiche Dichte wie die Erde ... vorausgesetzt, dass die Gravitationskompression auf ihm die gleiche ist wie auf unserem Planeten, aber das ist falsch, weil der Planet daher eine größere Schwerkraft haben wird , wird die Gravitationskompression darin intensiver sein, wodurch es kleiner und dichter als erwartet wird (ist diese Argumentation richtig?). Aber was mir wichtig ist, ist nicht nur die Dichte des Planeten, sondern auch die seines Kerns, der ebenfalls zunehmen würde.

Der Erdkern hat eine Dichte von 11.000 kg/m³, mehr als zu erwarten wäre, wenn man ihn nur als Gemisch aus zwei Stoffen (vereinfacht Fe und Ni) ohne Berücksichtigung der Gravitationskompression berechnen würde. Wie kann ich die Dichte des Kerns berechnen? P (und von sich aus) unter Berücksichtigung der Gravitationskompression? Ist das ohne großen Aufwand möglich?

Es gibt keine praktische Gleichung für die Beziehung zwischen diesen Dingen, einfach weil reale Materialien keine konstante Dichte oder Temperatur haben, wenn sich der Radius (auf dem Planeten) ändert. Hier ist ein Link zu einer Seite mit Links zu verschiedenen Planetenmodellen, die Ihnen eine Vorstellung von der Komplexität geben. Es gibt eine (nicht leicht verständliche) Anwendung, die über Wolfram Player verwendet werden kann, wenn Sie es ausprobieren möchten. Ich vermute, dass dieses Material zu komplex für Sie sein wird, aber schauen Sie es sich an.
Eine heikle Frage wegen der Phasenumwandlung: Zumindest Eisen kann offenbar unter hohem Druck von flüssig wieder in kristallinen Feststoff übergehen. phys.org/news/2017-02-theory-earth-core-solid-extreme.html . Die Phase kann definitiv die Dichte beeinflussen. Ich schlage das vor? wäre gut für den planetaren Wissenschaftsstapel, da es definitiv harte Wissenschaft ist, aber ich weiß nicht, wie ich es migrieren soll.
Vielleicht die Art und Weise, wie die Dichte der Erde kartiert wurde? Dazu gehörten viele Dinge wie das Verfolgen von Erdbebenwellen, das Messen exakter Schwerkraftwerte usw. Den Wert im Kern zu kennen, ist ohne Messungen eine chaotische Sache.

Antworten (1)

Die Technik, um dies zu tun, ähnelt derjenigen, die beim Konstruieren von Sternenmodellen verwendet wird. Sie kennen einige der Eigenschaften Ihres Objekts - in diesem Fall sehen sie aus wie Masse und Radius. Sie möchten die innere Struktur des Planeten herausfinden, einschließlich der zentralen Dichte und des Drucks sowie der Dichte- und Druckprofile als Funktionen des Radius. Der beste Weg, dies zu tun, ist wie folgt:

  1. Bestimmen Sie eine Beziehung zwischen Dichte und Druck - was wir eine Zustandsgleichung nennen. Wir schreiben dies als ρ = F ( P ) .
  2. Schätzen Sie den zentralen Druck - und damit die zentrale Dichte - anhand der obigen Beziehung.
  3. Integrieren Sie numerisch die Differentialgleichungen, die die Struktur des Planeten bestimmen. Beginnen Sie im Kern und gehen Sie nach außen.
  4. Sobald Ihre Gleichungen anzeigen, dass Sie die Oberfläche des Planeten erreicht haben, sehen Sie sich die vorhergesagte Masse und den vorhergesagten Radius an.
  5. Passen Sie angesichts des Ergebnisses von Schritt 4 Ihre Schätzungen für den zentralen Druck und die Dichte an und wiederholen Sie die Schritte 3 und 4.

Die beiden Gleichungen, die Sie neben der Zustandsgleichung benötigen, sind die Gleichung des hydrostatischen Gleichgewichts ( 1 ) und die Massenkontinuitätsgleichung ( 2 ) :

(1) D P D R = G M R ρ R 2
(2) D M D R = 4 π R 2 ρ
Hier, P ist Druck, ρ ist Dichte, R ist Radius, G ist die Gravitationskonstante, und M R ist die in einem Radius enthaltene Masse R :
M R = 0 R 4 π R ' 2 ρ ( R ' ) D R '
Sagen wir das bei einem bestimmten Radius R ich , Sie kennen den Druck, die Dichte und die Masse, die in diesem Radius eingekapselt sind - P ich , ρ ich , Und M R , ich . Angenommen, Ihre Integration verwendet eine gewisse Schrittgröße Δ R . Dann die Werte am Radius R ich + 1 Sind
R ich + 1 = R ich + Δ R
P ich + 1 = P ich + Δ P = P ich + D P D R Δ R = P ich G M R , ich ρ ich R ich 2 Δ R
ρ ich + 1 = F ( P ich + 1 )
M R , ich + 1 = M ich + Δ M = M ich + D M D R Δ R = M ich + 4 π R ich 2 ρ ich Δ R
Fahren Sie damit fort, bis Sie die Oberfläche des Planeten erreichen, dh bei P = 0 , und wiederholen Sie dann die Schritte 3 bis 4 mit einer neuen Schätzung, bis Sie mit der vorhergesagten Masse und dem Radius des Modells zufrieden sind.

Ich habe dieses Verfahren in meiner Antwort auf eine andere Frage implementiert , obwohl ich die Schritte 3 bis 4 nur einmal ausführen musste - ich hatte einen vernünftigen zentralen Druck, der alles herausgefunden hatte. Ich habe eine Zustandsgleichung von Seager et al. 2008 , entwickelt für felsige Planeten.

Sie brauchen Informationen, die Sie nicht haben. Du schreibst ρ ( R ) aber es ist auch eine Funktion der Zusammensetzung. Das weißt du nicht.
@puppetsock Das OP hat vermutlich eine Zusammensetzung gewählt (z. B. Eisen, Perowskit usw.). Das verlinkte Papier gibt Beispielparameter für verschiedene Zusammensetzungen (siehe Tabelle 3), und es ist ganz einfach, je nach Zusammensetzung des Planeten einfach die richtigen auszuwählen (oder vielleicht Zwischenwerte zu wählen, wenn der Planet eine heterogenere Zusammensetzung hat). Es ist nicht schwer, eine richtige zu finden ρ ( P ) zusammensetzungsabhängige Zustandsgleichung.
Es ist nicht schwer, eine richtige zu finden ρ ( P ) Zusammensetzungsabhängige Zustandsgleichung Es mag für Sie nicht schwer sein, aber ich habe den Eindruck, dass dies für das OP sehr schwierig ist.
@StephenG Ich habe einen Artikel mit mehreren verschiedenen EOSs verlinkt, die die meisten felsigen Planeten abdecken. . . Ich bin mir nicht sicher, was die verbleibenden Einwände sind. Wenn das OP eine bestimmte Zusammensetzung hat, kann ich natürlich nach einer geeigneten suchen, wenn sie nicht in der Arbeit behandelt wird, aber im Moment bin ich mir nicht sicher, ob die aktuellen Optionen unzureichend sind. :-/
Ich werde die Antwort innerhalb von ein paar Tagen mit einem detaillierteren Abschnitt zu EOSs aktualisieren. Ich denke, zukünftige Leser könnten das über das OP hinaus nützlich finden, daher werde ich nach Möglichkeit näher darauf eingehen.