Détermination de l'effet d'une petite force variable sur la précession planétaire du périhélie

14

Existe-t-il une technique analytique pour déterminer l'effet d'une petite accélération transversale variable sur le taux de précession des aspides (strictement pas une précession mais la rotation de la ligne d'aspides) d'une planète en orbite autour du Soleil dans un plan 2D selon la loi de gravité newtonienne ?

J'ai modélisé ces effets dans un modèle informatique réitératif et je voudrais vérifier ces mesures.

La formule d'accélération transversale est

UNEt=(K/c2)VrVtUNEr.

Où:-

c est la vitesse de la lumière,

K est une constante de grandeur entre 0 et +/- 3, de telle sorte que K/(c2)<<1 .

Ar est l'accélération de la planète vers le Soleil due à l'influence gravitationnelle newtonienne du Soleil, ( UNEr=gM/r2 ).

Vr est la composante radiale de la vitesse de la planète par rapport au Soleil (+ = mouvement loin du Soleil)

Vt est la composante transversale de la vitesse de la planète par rapport au Soleil (+ = direction du mouvement vers l'avant de la planète le long de sa trajectoire orbitale). Vectorialement Vt = V - Vr où V est le vecteur vitesse instantanée totale de la planète par rapport au Soleil.

Supposons que la masse de la planète est petite par rapport au Soleil

Aucun autre organisme n'est dans le système

Tous les mouvements et accélérations sont confinés au plan bidimensionnel de l'orbite.

MISE À JOUR

La raison pour laquelle cela m'intéresse est qu'une valeur de K = +3 dans mon modèle informatique produit des valeurs anormales (non newtoniennes) de taux de rotation périaphés très proches d'environ 1% de celles prédites par la relativité générale et de quelques pour cent de ceux observés par les astronomes (Le Verrier, mis à jour par Newcomb).

Formule (Einstein, 1915) pour la rotation péri-dérivée dérivée de GR (radians par orbite) de http://en.wikipedia.org/wiki/Apsidal_precession

ω=24.π3.une2.T-2.c-2.(1-e2)-1

MISE À JOUR 4

J'ai accepté la réponse de Walter. Non seulement a-t-il répondu à la question initiale (Existe-t-il une technique ...?) Mais aussi son analyse produit une formule qui non seulement confirme les effets simulés par ordinateur de la formule d'accélération transversale (pour K = 3) mais qui aussi (de manière inattendue pour moi) est essentiellement équivalent à la formule d'Einstein 1915.

du résumé de Walter (dans la réponse de Walter ci-dessous): -

: (à partir d'une analyse de péturbation du premier ordre) le demi-grand axe et l'excentricité sont inchangés, mais la direction du périapse tourne dans le plan de l'orbite au rythme où est la fréquence orbitale et avec axe semi-majeur. Notez que (pour ) cela correspond au taux de précession de la relativité générale (GR) à l'ordre (donné par Einstein 1915).Ωvc=ΩaaK=3v 2 c /c2

ω=Ωvc2c2K1-e2,
Ωvc=ΩuneuneK=3vc2/c2
steveOw
la source
Cherchez-vous toujours une réponse?
Walter
@Walter. Oui. J'ai posé une question similaire sur physics.stackexchange.com/questions/123685/… mais aucune réponse solide n'a encore été reçue.
steveOw
@Walter. J'ai également demandé à math.stackexchange.com/questions/866836/… .
steveOw
Oui, il existe des méthodes analytiques approximatives (théorie des perturbations), valables dans la limite de . Vous pouvez peut-être clarifier un peu votre question. Quelle est la direction de l'accélération transversale (je comprends que «transverse» signifie perpendiculaire à la vitesse instantanée, mais il n'est pas clair si l'accélération est dans le plan de l'orbite ou perpendiculaire ou un mélange). K1
Walter
Il y a une différence entre votre question ici et celle sur les mathématiques (et la physique): ici l'accélération transversale est proportionnelle à l'accélération radiale et est un nombre sans dimension, là l'accélération radiale n'a aucun effet sur l'accélération transversale et K doit être un accélération (bien que vous parliez d'un «nombre»). KK
Walter

Réponses:

5

Vous voudrez peut-être utiliser la théorie des perturbations . Cela vous donne seulement une réponse approximative , mais permet un traitement analytique. Votre force est considérée comme une petite perturbation de l'orbite elliptique Keplerian et les équations résultantes du mouvement sont étendues en puissances de . Pour la théorie de la perturbation linéaire, seuls les termes linéaires dans K sont conservés. Cela conduit simplement à intégrer la perturbation le long de l'orbite d'origine non perturbée. En écrivant votre force comme vecteur, l'accélération perturbatrice est a = K G MKK avecvr=v r la vitesse radiale (v ˙ r ) et vt=(v - r (v r ))la composanterotation devitesse (la vitesse complète moins la vitesse radiale). Ici, le point ci-dessus indique une dérivée temporelle et un chapeau le vecteur unitaire.

a=KGMr2c2vrvt
vr=vr^vr˙vt=(v-r^(vr^))

Maintenant, cela dépend de ce que vous entendez par « effet ». Étudions les changements de l'axe semi-principal orbital , de l'excentricité e et de la direction du périapse.unee


Pour résumer les résultats ci - dessous : le demi-grand axe et l'excentricité sont inchangés, mais la direction du périapse tourne dans le plan de l'orbite à la vitesse Ωest la fréquence orbitale etvc=Ωaavecunaxe semi-majeur. Notez que (pourK=3) celacorrespond au taux de précession de la relativité générale (GR)à l'ordrev 2 c /c2(donné par Einstein 1915 mais non mentionné dans la question d'origine).

ω=Ωvc2c2K1e2,
Ωvc=ΩaaK=3vc2/c2

changement d'axe semi-majeur

De la relation (avec E = 1a=GM/2El'énergie orbitale) que nous avons pour le changement deadû à une accélération externe (non Keplerian) ˙ a =2a2E=12v2GMr1a En inséranta(notez quevvt=h2/r2avec le vecteur de moment angulairehrv), on obtient ˙ a =2a2Kh2

a˙=2a2GMva.
avvt=h2/r2hrv Etant donnéla moyenneorbitevrf(r)=0pour toute fonctionf(voir cidessous), ˙ un=0.
a˙=2a2Kh2c2vrr4.
vrf(r)=0fa˙=0

changement d'excentricité

De , on trouve e ˙ e = - h ˙ hh2=(1e2)GMune Nous savons déjà que ˙ un=0, donc seulement besoin de considérer le premier terme. Ainsi, e ˙ e =-(rv)(ra)

ee˙=-hh˙gMune+h2une˙2gMune2.
une˙=0 où j'ai utilisé l'identité (ab)(cd)=ac
ee˙=-(rv)(rune)gMune=-r2vunegMune=-Kh2unec2vrr2,
et le fait r a p = 0 . Encorefoisv r / r 2= 0 et donc˙ e= 0 .(uneb)(c)=unecb-unebcrunep=0vr/r2=0e˙=0

changement de direction du périapse

evh/GMr^e

e˙=a(rv)+v(ra)GM=2(va)r(rv)aGM=2Kc2h2vrrr4Kc2vr2vtr
a(bc)=(ac)b(ab)cra=0e˙=ωe
ω=ΩKvc2c2(1e2)1h^.
ω=|ω|ee˙=ee˙=0

K(vc/c)20aeaeK

mentions légales Aucune garantie que l'algèbre est correcte. Vérifie ça!


Annexe: moyennes orbitales

vrF(r)F(r)F(r)F(r)F=F

vrF(r)=1T0Tvr(t)F(r(t))t=1T[F(r(t))]0T=0
T

e˙

r=une((cosη-e)e^+1-e2péchéηk^)etr=une(1-ecosη)
ek^h^e^ehη=ηesinη,d=(1ecosη)dη
=(2π)102πd=(2π)102π(1ecosη)dη.
˙=Ω=gM/une3r
v=vc1-e2cosηk^-péchéηe^1-ecosη
vcΩune=gM/uneunevr=r^v=vcepéchéη(1-ecosη)-1
vt=vc1-e2(cosη-e)k^-(1-e2)péchéηe^(1-ecosη)2.

h2vrrr4=Ωvc2k^e(1e2)3/22π02πsin2η(1ecosη)4dη=Ωvc2e2(1e2)k^vr2vtr=Ωvc2k^e2(1e2)1/22π02πsin2η(cosηe)(1ecosη)4dη=0,
e^
2h2vrrr4vr2vtr=Ωvc2ek^(1e2)
Walter
la source
Les commentaires ne sont pas pour une discussion approfondie; cette conversation a été déplacée vers le chat .
called2voyage