Das Geopotenzialmodell von Wikipedia gibt die Form der verknüpften Potenzialkomponente an als
und die Beschleunigungskomponenten (Gradient dieses Potentials) sind:
Aus dem Geopotentialmodell und dieser Antwort .
Ich habe geplant Und unten für die Flugzeug, wo ist die Rotationsachse der Erde. Die Verteilung in drei Dimensionen wird einfach um diese gedreht Achse.
Frage: ist der dominierende Teil des Quadrupolmoments der Erde, aber es ist nur eine Komponente.
import numpy as np
import matplotlib.pyplot as plt
Re = 6378137. # meters
J2 = 1.7555E+25 # m^5/s^2
hw = 400
x = np.linspace(-2.5*Re, 2.5*Re, 2*hw+1)
x, z = np.meshgrid(x, x)
y = np.zeros_like(x)
x2, y2, z2 = [w**2 for w in (x, y, z)]
r2 = x2 + y2 + z2
r = np.sqrt(r2)
r7 = r**7
u = J2 * r**-5 * 0.5 * (3*z2 - r**2)
ax = J2 * (x/r7) * (6*z2 - 1.5*(x2+y2))
ay = J2 * (y/r7) * (6*z2 - 1.5*(x2+y2))
az = J2 * (z/r7) * (3*z2 - 4.5*(x2+y2))
for thing in (u, ax, ay, az):
thing[r<Re] = np.nan
axz = np.sqrt(ax**2 + az**2)
if True:
plt.figure()
plt.subplot(1, 2, 1)
umax = np.nanmax(np.abs(u))
plt.imshow(u, cmap='seismic', vmin=-umax, vmax=umax)
plt.colorbar()
plt.title('specific potential energy (J/kg)', fontsize=16)
plt.plot([hw, hw], [0.3*hw, 1.7*hw], '-k')
plt.xlim(0, 2*hw+1)
plt.ylim(0, 2*hw+1)
plt.subplot(1, 2, 2)
plt.imshow(axz, cmap='inferno')
plt.colorbar()
plt.title('acceleration magnitude (m/s^2)', fontsize=16)
plt.plot([hw, hw], [0.3*hw, 1.7*hw], '-k')
plt.xlim(0, 2*hw+1)
plt.ylim(0, 2*hw+1)
plt.suptitle("Earth's geopotential's J2 component", fontsize=16)
plt.show()
Repräsentiert J 22 den Rest des Quadrupolmoments der Erde?
Mit einer gebräuchlicheren Schreibweise l Grad , m Ordnung gibt es die zonalen l 2 m 0 , tesseralen l 2 m 1 und sektoralen l 2 m 2 sphärischen Harmonischen, aber Sie kennen sie wahrscheinlich bereits.
Wie sieht J22 aus ? Wie ist seine Form und Symmetrie?
Hier ist, was ich für l 2 m 2 bekomme :
Zum Vergleich, hier ist, was ich für Ihr Diagramm "Beschleunigungsgröße" (l 2 m 0 ) bekomme (es ist nicht klar, was 800 darstellt):
Wikipedia geht tatsächlich detailliert auf die Komponenten des Gravitationskraftfeldes der Erde ein.
In der Wikipedia-Referenz haben wir nicht für . Vielmehr wird dieser Term in eine Cosinus-Komponente zerlegt und eine Sinuskomponente . Diese Tabelle enthält Ergebnisse von einem aus einer Reihe von Modellen, die als JGM-Modelle bezeichnet werden und vom Goddard Space Center entwickelt wurden:
Zonale Koeffizienten
Tesserale Koeffizienten [sic; Fälle mit sind sektoral]
Hier bedeutet Multiplizieren mit einer Zehnerpotenz, deren Exponent also folgt bedeutet multiplizieren mit . Die absoluten Werte der obigen Zahlen stellen die relative Stärke des Gravitationspotentials dar, das sich aus der sphärischen harmonischen Komponente ergibt, in einem Abstand, der gleich dem durchschnittlichen Erdradius ist, und mit der isotropen Komponente "sphärische Masse" (was wäre in dieser Nomenklatur) skaliert auf .
Nehmen wir nun an, Sie bringen einen Satelliten in eine niedrige Erdumlaufbahn, in einem radialen Abstand, der im Wesentlichen gleich der oben verwendeten Referenz ist. Unter den Quadrupolkomponenten eindeutig ist am größten, mit Und viele Größenordnungen kleiner, und Und dazwischen fallen. Sie können den Quadrupol daher vertretbar nur mit dem beschreiben ( ) Und ( , ) Komponenten, wobei die sehr kleinen ignoriert werden Stücke.
Aber das ist nicht die ganze Geschichte. Während die verbleibende Quadrupolkomponente in der Nähe der Erdoberfläche winzig ist, einige polare Komponenten höherer Ordnung ( ) sind nicht. Wenn Sie die einschließen Komponenten des Quadrupols zur Modellierung Ihrer erdnahen Bahnbewegung sollten Sie auch mindestens die ähnlich großen einbeziehen , Und Begriffe, um Ihre Genauigkeit gegenüber dem rotationssymmetrischen Quadrupol wirklich zu verbessern. Die Abweichungen von einem rituell symmetrischen Quadrupol sind komplexer als nur eine zusätzliche Quadrupolkomponente.
Der Grund dafür ist, dass die Abweichungen von einem rotationssymmetrischen Quadrupol durch die Verteilung von großräumigen geografischen Merkmalen wie Gebirgszügen, Becken und tektonischen Subduktionsgrenzen bestimmt werden. Solche Merkmale auf der Erde sind zu komplex, um von einer oder zwei sphärischen harmonischen Komponenten erfasst zu werden. Wenn Sie die geologische Komplexität der Erde erfassen müssen, um Ihren Satelliten zu modellieren, benötigen Sie dazu viele sphärische harmonische Terme.
Sie können die Dinge vereinfachen, wenn Sie stattdessen Ihren Satelliten in eine große Höhe bringen, sagen wir eine geostationäre Umlaufbahn. Die Komponenten mit größer Werte verblassen schneller mit der Entfernung als solche mit kleineren Werten, sodass Sie diese je nach der von Ihnen benötigten Genauigkeit möglicherweise auslassen Und Komponenten. Die geografischen Warzen, die wir aus der Nähe sehen, werden durch eine weiter entfernte Betrachtung geglättet. Das ist, wo Sie am ehesten Erfolg haben werden , , Und .
Die Erde ist nicht der einzige Körper mit solchen Gravitationskomplexitäten. Jeder scheinbar hydrostatisch ausgeglichene, aber solide Körper wird wahrscheinlich seine eigenen geografischen Unebenheiten haben. Insbesondere Mars hat den großen Höhenunterschied zwischen Nord und Süd, das Hellas-Becken und Olympus Mons. Wenn wir einen Mars-Orbiter mit hoher Genauigkeit verfolgen wollen, sollten wir besser unsere sphärischen Harmonischen auffrischen.
äh
jet
. Ich erinnere mich, dass ich einen Tweet mit einem alten Haiku über Jet gelesen habe, aber ich kann es jetzt nicht finden. Ich habe "Don't use Jet!" in diesem und diesem Plot und entschuldigte sich dafür auch in einem von ihnen. Ich werde nach einer besseren cmap suchen und updaten ;-)äh
Gerrit
äh