Diferansiyel Denklemlerle Yörünge Çözümü
Kütleçekim altında hareket eden bir cismin yörüngesi, Newton denklemlerinin diferansiyel formda çözülmesiyle elde edilir. Bu yaklaşım, iki cisim probleminden başlayarak daha karmaşık sistemlere kadar genişletilebilir.
Temel fikir şudur: kuvvet biliniyorsa hareket denklemi yazılır, diferansiyel denklem çözülerek yörünge bulunur.
Newton Hareket Denklemi
İki cisim problemi, indirgenmiş kütle kullanılarak tek cisim problemine indirgenebilir.
Newton yasası:
$$ m \frac{d^2 \vec{r}}{dt^2} = – \frac{G M m}{r^3} \vec{r} $$
Sadeleştirilirse:
$$ \frac{d^2 \vec{r}}{dt^2} = – \frac{G M}{r^3} \vec{r} $$
Bu ikinci mertebeden, doğrusal olmayan bir diferansiyel denklemdir.
Açısal Momentumun Korunumu
Merkezcil kuvvet altında açısal momentum korunur:
$$ \vec{L} = \vec{r} \times m \vec{v} $$
Büyüklük olarak:
$$ L = m r^2 \dot{\theta} $$
Buradan:
$$ r^2 \dot{\theta} = \frac{L}{m} $$
Bu sonuç problemi iki boyuta indirger.
Polar Koordinatlara Geçiş
Konum vektörü polar koordinatlarda yazılır:
$$ \vec{r} = r \hat{e}_r $$
Radyal denklem:
$$ \ddot{r} – r \dot{\theta}^2 = – \frac{G M}{r^2} $$
Açısal momentum bağıntısı yerine konursa:
$$ \ddot{r} – \frac{L^2}{m^2 r^3} = – \frac{G M}{r^2} $$
Bu denklem yörüngenin temel diferansiyel formudur.
Değişken Dönüşümü
Çözümü kolaylaştırmak için:
$$ u = \frac{1}{r} $$
tanımı yapılır.
Türev dönüşümü sonucunda şu denklem elde edilir:
$$ \frac{d^2 u}{d\theta^2} + u = \frac{G M m^2}{L^2} $$
Bu lineer ikinci mertebe diferansiyel denklemdir.
Genel Çözüm
Denklemin genel çözümü:
$$ u(\theta) = \frac{G M m^2}{L^2} \left( 1 + e \cos\theta \right) $$
Dolayısıyla:
$$ r(\theta) = \frac{L^2}{G M m^2} \frac{1}{1 + e \cos\theta} $$
Bu ifade konik kesit denklemdir.
Yörünge Türleri
Eksantriklik ( e ) yörünge tipini belirler:
- ( 0 \le e < 1 ) → Elips
- ( e = 1 ) → Parabol
- ( e > 1 ) → Hiperbol
Yarı latus rektum:
$$ p = \frac{L^2}{G M m^2} $$
Yörünge denklemi:
$$ r = \frac{p}{1 + e \cos\theta} $$
Enerji ile Bağlantı
Toplam enerji:
$$ E = \frac{1}{2} m v^2 – \frac{G M m}{r} $$
Eksantriklik enerji ve açısal momentuma bağlıdır:
$$ e = \sqrt{1 + \frac{2 E L^2}{m (G M)^2}} $$
Buradan enerji işareti ile yörünge tipi arasındaki ilişki açıkça görülür.
Sayısal Çözüm
Analitik çözüm yalnızca ters kare kuvvet için mümkündür. Genel durumda diferansiyel denklem sayısal olarak çözülür:
$$ \frac{d^2 \vec{r}}{dt^2} = \vec{a}(\vec{r}) $$
Yaygın yöntemler:
- Euler yöntemi
- Runge-Kutta yöntemleri
- Leapfrog integratörü
Enerji korunumu kontrolü için:
$$ \frac{dE}{dt} = 0 $$
yaklaşık olarak sağlanmalıdır.
Özet
Yörünge problemi şu diferansiyel denklemle başlar:
$$ \frac{d^2 \vec{r}}{dt^2} = – \frac{G M}{r^3} \vec{r} $$
Değişken dönüşümü ile:
$$ \frac{d^2 u}{d\theta^2} + u = sabit $$
Çözüm:
$$ r = \frac{p}{1 + e \cos\theta} $$
Bu sonuç, Newton mekaniğinin diferansiyel denklemlerle çözüldüğünde konik kesit yörüngeler verdiğini gösterir. Diferansiyel denklemler yaklaşımı, hem analitik hem sayısal gök mekaniğinin temelini oluşturur.