Kepler yörüngesi

küçükresim|300px|Kepler Orbitand eksentriklerini çeşitli biçimlerde gösteren diyagram. Mavi bir hiperbolik yörünge (e> 1) 'dir. Yeşil bir parabolik yörünge (e = 1) 'dir. Kırmızı eliptik yörünge (0 r(\theta) = \frac{a(1-e^2)}{1+e\cos(\theta)} *r mesafedir *a yörüngenin boyutunu tanımlayan yarı büyük eksendir *e yörüngenin şeklini tanımlayan dış merkezliktir *\theta yörüngedeki nesnenin mevcut konumu ile yörüngedeki merkez cisme en yakın olduğu (enberi noktası olarak adlandırılan) konum arasındaki açı olan gerçek ayrıklıkdır. Alternatif olarak, denklem şu şekilde ifade edilebilir: : r(\theta) = \frac{p}{1+e\cos(\theta)} Burada p, eğrinin yarı özkirişi olarak adlandırılır. Denklemin bu formu, özellikle yarı büyük ekseni sonsuz olan parabolik yörüngelerle uğraşırken kullanışlıdır. Bu yasaları gözlemleriyle geliştiren Kepler, hiçbir zaman bu hareketleri açıklayacak bir teori geliştirememiştir.Bate, Mueller, White. s. 177-181 === Isaac Newton === Isaac Newton, 1665 ve 1666 yılları arasında hareket, kütleçekim ve diferansiyel hesap ile ilgili çeşitli kavramlar geliştirdi. Bununla birlikte bu kavramlar, 1687 yılında hareket yasalarını ve evrensel kütleçekim yasasını ana hatlarıyla belirttiği Principia'ya kadar yayınlanmadı. Üç hareket yasasından ikincisi şöyle der: > Bir cismin ivmesi, cisme etki eden net kuvvetle paralel ve doğru > orantılıdır. Net kuvvet yönündedir ve cismin kütlesi ile ters orantılıdır: : > \mathbf{F} = m\mathbf{a} = m\frac{d^2\mathbf{r}}{dt^2} *\mathbf{F} kuvvet > vektörüdür *m kuvvetin etki ettiği cismin kütlesidir *\mathbf{a} ivme > vektörüdür, \mathbf{r} konum vektörünün ikinci zaman türevidir. Kurallara bakılırsa, denklemin bu şekli yalnızca aşağıda yapılan basitleştirici varsayımlara dayalı olarak doğru olan sabit kütleli bir nesne için geçerlidir. Newton'un kütleçekim yasası şöyle der: > Her nokta kütle, her iki noktayı kesen doğru boyunca işaret edilen bir > kuvvetle diğer nokta kütleleri çeker. Kuvvet, iki kütlenin çarpımı ile > orantılı ve nokta kütleler arasındaki uzaklığın karesi ile ters orantılıdır: > :F = G \frac{m_1 m_2}{r^2} * F iki nokta kütlesi arasındaki kütleçekim > kuvvetinin büyüklüğüdür. * G kütle çekimi sabitidir * m_1 birinci nokta > kütlesinin kütlesidir * m_2 ikinci nokta kütlesinin kütlesidir * r iki nokta > kütlesi arasındaki mesafedir Newton, hareket yasaları ve evrensel kütleçekim yasasından astronomide yörünge hareketine özgü Kepler yasalarını türetebildi. Kepler yasaları gözlem verileriyle iyi bir şekilde desteklendiğinden, bu tutarlılık Newton'un genelleştirilmiş teorisinin ve birleşik göksel ve sıradan mekaniğin geçerliliğine güçlü bir destek sağladı. Bu hareket yasaları, 20. yüzyılın başlarında Albert Einstein özel ve genel görelilik kavramlarını ortaya koyana kadar modern gök mekaniğinin temelini oluşturdu. Çoğu kullanımda Kepler hareketi, gezegenlerin ve uyduların hareketlerini nispeten yüksek doğruluk derecelerine yaklaştırır ve astronomi ve astrodinamikte yaygın olarak kullanılır. === Sadeleştirilmiş iki cisim problemi === İki cisimli bir sistemde bir cismin hareketini çözmek için iki sadeleştirici varsayım yapılabilir: # Cisimler küresel olarak simetriktir ve nokta kütleler olarak ele alınabilir. # Cisimlere karşılıklı kütleçekimi dışında etki eden iç ve dış kuvvetler yoktur. Büyük gök cisimlerinin şekilleri kürelere yakındır. Simetri ile, bir kütle noktasını homojen bir küreye doğru çeken net kütleçekim kuvveti, merkeze doğru yönlenmelidir. Kabuk teoremi (Isaac Newton tarafından da kanıtlanmış olan), kürenin yoğunluğu derinliğe göre değişse bile (çoğu gök cisimleri için olduğu gibi) bu kuvvetin büyüklüğünün, kütlenin tamamı kürenin ortasında toplanmış gibi aynı olduğunu belirtir. Bundan hemen sonra iki homojen küre arasındaki çekim, her ikisinin de kütlesi merkezinde yoğunlaşmış gibidir. Asteroitler veya uzay araçları gibi daha küçük nesneler genellikle bir küreden güçlü bir şekilde ayrılan bir şekle sahiptir. Ancak bu düzensizlikler tarafından üretilen kütleçekim kuvvetleri, merkezi cismin kütleçekimi ile karşılaştırıldığında genellikle önemsizdir. Düzensiz bir şekil ile mükemmel bir küre arasındaki fark da mesafelerle azalır ve yörünge mesafelerinin çoğu, yörüngedeki küçük bir cismin çapıyla karşılaştırıldığında çok büyüktür. Bu nedenle bazı durumlar için şekil düzensizliği, kesinlik üzerinde önemli bir etki olmaksızın ihmal edilebilir. Bu etki, özellikle alçak yörüngelerde bulunan yapay Dünya uyduları için oldukça belirgindir. Gezegenler değişen hızlarda dönerler ve bu merkezkaç kuvveti nedeniyle hafif basık bir şekil alabilirler. Böyle yassı bir şekille kütle çekimi, homojen bir küreninkinden bir miktar sapacaktır. Daha büyük mesafelerde bu yassılığın etkisi ihmal edilebilir hale gelir. Güneş Sistemindeki gezegensel hareketler, nokta kütleler olarak ele alınırsa yeterli hassasiyetle hesaplanabilirler. Kütleleri m_1 ve m_2 olan iki nokta kütle nesnesi ve bazı eylemsiz referans çerçevesine göreceli konum vektörleri \mathbf{r}_1 ve \mathbf{r}_2, kütleçekim kuvvetlerine maruz kalır: : m_1 \ddot{\mathbf{r}}_1 = \frac{-G m_1 m_2}{r^2} \mathbf{\hat{r}} : m_2 \ddot{\mathbf{r}}_2 = \frac{ G m_1 m_2}{r^2} \mathbf{\hat{r}} burada \mathbf{r}, kütle 1'in kütle 2'ye göre göreceli konum vektörüdür ve şu şekilde ifade edilir: : \mathbf{r} = \mathbf{r}_1 - \mathbf{r}_2 \mathbf{\hat{r}}, o yöndeki birim vektördür ve r, bu vektörün uzunluğudur. Kendi kütlelerine bölmek ve ikinci denklemi birinciden çıkarmak, birinci cismin ikinciye göre ivmesi için hareket denklemini verir: = - \frac{\alpha}{r^2} \mathbf{\hat{r}}|}} \alpha, yerçekimi parametresidir ve eşittir : \alpha = G(m_1 + m_2) Birçok uygulamada, üçüncü bir sadeleştirici varsayım yapılabilir: # * Merkezi cisim ile karşılaştırıldığında yörüngedeki cismin kütlesi önemsizdir. Matematiksel olarak, m1 >> m2, yani 'dir. Genellikle \mu = G\,M olarak gösterilen bu tür standart kütleçekim parametreleri, yörüngedeki uydularından çok daha büyük kütlelere M sahip olan Güneş, büyük gezegenler ve Ay için yaygın olarak kullanılır. Bu varsayım, sadeleştirilmiş iki cisim problemini çözmek için gerekli değildir, ancak özellikle Dünya yörüngesindeki uydular ve Güneş'in yörüngesindeki gezegenlerle ilgili hesaplamaları basitleştirir. Jüpiter'in kütlesi bile Güneş'inkinden 1047 kat daha azdır, bu da α değerinde %0,096'lık bir hata oluşturur. Dikkate değer istisnalar arasında Dünya-Ay sistemi (kütle oranı 81,3), Plüton-Charon sistemi (kütle oranı 8,9) ve ikili yıldız sistemleri sayılabilir. Bu varsayımlar altında iki cisim problemi için diferansiyel denklem matematiksel olarak tamamen çözülebilir ve Kepler'in gezegensel hareket yasalarını takip ederek sonuçta oluşan yörüngeye "Kepler yörüngesi" denir. Tüm gezegenlerin yörüngeleri yüksek doğrulukta Güneş etrafındaki Kepler yörüngeleridir. Küçük sapmalar, gezegenler arasındaki çok daha zayıf kütleçekim etkilerinden ve Merkür durumunda genel görelilikten kaynaklanmaktadır. Dünya etrafındaki yapay uyduların yörüngeleri, makul bir tahminle Güneş, Ay ve Dünya'nın yassılığından dolayı küçük düzensizliklere sahip Kepler yörüngeleridir. Hareket denkleminin bütün kütleçekim ve kütleçekim olmayan kuvvetlerin (güneş radyasyon basıncı ve atmosfer direnci gibi) hesaba katılmasıyla sayısal olarak entegre edilmesi gereken yüksek doğruluktaki kullanımlarda, Kepler yörünge kavramları çok önemlidir ve yoğun olarak kullanılır. === Kepler yörünge öğeleri === küçükresim|300px|sağ|Kepler yörünge öğeleri Kepler yörüngesinin altı parametre ile tanımlanabilirliğinden söz etmek gerekir. Üç boyutlu uzayda hareket eden bir nesnenin hareket bir konum vektörü ve bir hız vektörü ile karakterize edilir. Her vektörün üç bileşeni vardır, bu yüzden uzayda bir yörünge tanımlamak için gerekli değerlerin toplam sayısı altıdır. Bir yörünge (Kepler Elementleri olarak bilinen) ve pozisyon ve hızına bağlı olarak ölçülebilen altı elementten oluşur ki bunların üçü daha önce gösterilmiştir. Bu altı uygun elementten beş tanesi yörünge için değişmezdir. (iki sürekli değişen vektöre bir tezat niteliğinde olmak üzere). Kendi yörüngesinde içindeki bir nesnenin gelecekteki konumu tahmin edilebilir ve yeni konum ve hız kolayca yörünge elemanları ile elde edilebilir. Bunlardan ikisi eksenin boyutunu ve şeklin belirler: * Yarımajör eksen (a\,\\!) * eksantriklik (e\,\\!) Üç tanesi yörünge alanının oryantasyonunu belirler: * eğim (i\,\\!)yörünge düzlemi ile referans düzlemi arasındaki açıyı tanımlar. * Yükselen nodun uzunluğu (\Omega\,\\!) Referans yönü ve referans düzlemi üzerinde yörüngeye yukarı kesişme (yükselen nod) arasındaki açıyı tanımlar. * Periapsis argüman(\omega\,\\!) artan nod ve periapsis arasındaki açıyı belirler. Ve son olarak: * Gerçek anomali ( u)periapsis ölçülen yörünge boyunca yörüngedeki cismin konumunu tanımlar. Çeşitli başka değerlerin yerine gerçek bir anomali kullanılabilir, bunlardan en sık kullanılanları mean anomoli olan M ve periapsisden yana geçen süre olan T’dir. i, \Omega ve \omega sadece yörünge düzlemi içinde nesnenin hareketini tartışırken referans çerçevesinde yörünge yönünü tanımlayan açısal ölçümlerdir ve kesin olarak gerekli değildir. Burada tamamlanmışlık adına belitilmişlerdir ancak aşağıdaki ispatlar için gerekli değillerdir. === Yukarıdaki diferansiyel denklemin () matematiksel çözümü === Herhangi bir merkez kuvveti altında yani r bir kuvvet paraleli ve hareket için, belirli bağıl açısal momentumu \mathbf{H} = \mathbf{r} \times {\dot{\mathbf{r}}} sabit kalır: \dot {\mathbf{H}} = \frac{d}{dt}\left(\mathbf{r} \times {\dot{\mathbf{r}}}\right) = \dot{\mathbf{r}} \times {\dot{\mathbf{r}}} + \mathbf{r} \times {\ddot{\mathbf{r}}} =\mathbf{0} + \mathbf{0} = \mathbf{0} Pozisyonun çapraz ürünü vektör ve onun hızı olduğu için sabit kaldığından, aynı düzelmde bulunmak zorundadırlar. ( \mathbf{H} ’ye ortogonal). Bu vektör, fonksiyonun, bir düzlem eğrisi olduğu anlamına gelir. Denklemin kökeni etrafında simetri vardır, çünkü kutupsal koordinatlarda çözmek kolaydır. Bununla birlikte, bu denklem dikkat etmek önemlidir doğrusal ivmeye refere eder \left (\ddot{\mathbf{r}} \right ), açısalın karşıtı olarak \left (\ddot{\theta} \right ) yahut radyal \left (\ddot{r} \right ) ivmelenme. Bu nedenle, denklemi değiştiren kişinin dikkatli olması gerekir.Şimdi kartezyen koordinat sistemine (\hat{\mathbf{x}} \ , \ \hat{\mathbf{y}}) ve kutupsal birim vektörlere (\hat{\mathbf{r}} \ , \ \hat{\boldsymbol\theta}) düzlemi octoganelinde bakalım \mathbf{H} : \hat{\mathbf{r}}=\cos(\theta)\hat{\mathbf{x}} + \sin(\theta)\hat{\mathbf{y}} \hat{\boldsymbol\theta}=-\sin(\theta)\hat{\mathbf{x}} + \cos(\theta)\hat{\mathbf{y}} Şimdi vektör \mathbf{r} fonksiyonunu yeniden yazabiliriz ve diferansiyali şöyledir as: \mathbf{r} =r ( \cos\theta \hat{\mathbf{x}} + \sin \theta \hat{\mathbf{y}}) = r\hat{\mathbf{r}} \dot{\mathbf{r}} = \dot r \hat {\mathbf r} + r \dot \theta \hat {\boldsymbol{\theta}} \ddot{\mathbf{r}} = (\ddot r - r\dot\theta^2)\hat{\mathbf{r}} + (r \ddot\theta + 2 \dot r \dot\theta) \hat{\boldsymbol\theta} (Vektör calculusunu inceleyiz). Bunları yerine yazarsak şunu buluruz: () (\ddot r - r\dot\theta^2)\hat{\mathbf{r}} + (r \ddot\theta + 2 \dot r \dot\theta) \hat{\boldsymbol\theta} = \left (-\frac{\mu}{r^2}\right )\hat{\mathbf{r}} + (0)\hat{\boldsymbol\theta} Bu sıradan olmayan polar diferansiyel denklem verir: ^2 = - \frac {\mu} {r^2}|}} Bu denklemi çözmek için, öncelikle her zaman diferansiyel lineerini ortadan kaldırmak gerekir. Şuna ulaşırız: H = |\mathbf{r} \times {\dot{\mathbf{r}}}| = |(r\cos(\theta), r\sin(\theta), 0) \times (\dot{r}\cos(\theta)-r\sin(\theta)\dot{\theta}, \dot{r}\sin(\theta)+r\cos(\theta)\dot{\theta}, 0)| = |(0,0,r^2\dot\theta)| = r^2\dot\theta ()’ün zaman diferansiyalinii aldığımız zaman, şuna ulaşırız: {r^3}|}} () ve (}) denklemleri ’ın diferansiyellerini elememize izin verir.. In order to eliminate the time derivatives of \theta’ın zaman diferansiyellerini ellimine etmek için uygun kısaltmaları bulmak bunun için de zincir kuralı kullanmamız gerekir. ^2 + \frac {dr} {d\theta} \cdot \ddot {\theta}|}} \ddot{r} - r {\dot{\theta}}^2 = - \frac {\mu} {r^2} \frac {d^2r} {d\theta^2} \cdot {\dot {\theta}}^2 + \frac {dr} {d\theta} \cdot \ddot {\theta} - r {\dot{\theta}}^2 = - \frac {\mu} {r^2} \frac {d^2r} {d\theta^2} \cdot \left (\frac{H}{r^2} \right )^2 + \frac {dr} {d\theta} \cdot \left (- \frac {2 \cdot H \cdot \dot{r}} {r^3} \right ) - r\left (\frac{H}{r^2} \right )^2 = - \frac {\mu} {r^2} Diferansiyel denklem () değişken değişimi ile analitik çözülebilir Bu verilenleri kullanarak, in ()’deki tüm zaman diferansiyelleri yoksayılabilir, ve bizi r \theta\,. \ddot{r} - r {\dot{\theta}}^2 = - \frac {\mu} {r^2} () ve ()’da belirtilenleri \frac {d^2r} {d\theta^2} ve \frac {dr} {d\theta} için kullanarak genel bir sonuca ulaşılır. e ve \theta_0\,’un başlangıçtaki s ve \frac {ds} {d\theta}’n başlangıç değerlerinin integral sabitleri olursa, integralinin sabitini kullanmak yerine birim vektörleri kullanıma alınır ve seçilen yörünge alanındaki koordinat sistemini belirler, böylece sıfır değerini alır ve e is pozitiftir. Bu şu anlama gelir: ’in maksimum olduğu yerde sıfırdır ve böylece minimumdur. p parametresini \frac {H^2}{\mu} olarak tanımlarsak şu ortaya çıkar. ==== Alternatif türevi ==== Bu denklemi polar diferansiyel denklemini kullanmadan çözmenin bir diğer yolu şöyledir: \mathbf{u}‘yu bir birim vektör olarak tanımlayın, örneğin,\mathbf{r} = r\mathbf{u} and \ddot{\mathbf{r}} = -\frac{\mu}{r^2}\mathbf{u} gibi. Burdan yola çıkarak \mathbf{H} = \mathbf{r} \times \dot{\mathbf{r}} = r\mathbf{u} \times \frac{d}{dt}(r\mathbf{u}) = r\mathbf{u} \times (r\dot{\mathbf{u}}+\dot{r}\mathbf{u}) = r^2(\mathbf{u} \times \dot{\mathbf{u}}) + r\dot{r}(\mathbf{u} \times \mathbf{u}) = r^2\mathbf{u} \times \dot{\mathbf{u}} Şimdi şunu değerlendiriniz \ddot{\mathbf{r}} \times \mathbf{H} = -\frac{\mu}{r^2}\mathbf{u} \times (r^2\mathbf{u} \times \dot{\mathbf{u}}) = -\mu\mathbf{u} \times (\mathbf{u} \times \dot{\mathbf{u}}) = -\mu[(\mathbf{u}\cdot\dot{\mathbf{u}})\mathbf{u}-(\mathbf{u}\cdot\mathbf{u})\dot{\mathbf{u}}] (Üçlü vektör ürünü). Şunu dikkate alınız \mathbf{u}\cdot\mathbf{u} = |\mathbf{u}|^2 = 1 \mathbf{u}\cdot\dot{\mathbf{u}} = \frac{1}{2}(\mathbf{u}\cdot\dot{\mathbf{u}} + \dot{\mathbf{u}}\cdot\mathbf{u}) = \frac{1}{2}\frac{d}{dt}(\mathbf{u}\cdot\mathbf{u}) = 0 Bu verileri bir önceki denkleme yerleştirdiğimizde \ddot{\mathbf{r}}\times\mathbf{H}=\mu\dot{\mathbf{u}} İki tarafın da integralini alırsak \dot{\mathbf{r}}\times\mathbf{H}=\mu\mathbf{u} + \mathbf{c} Burada c sabit vektördür. Bunu r' ile birleştirmek ortaya ilginç bir sonuç çıkarır \mathbf{r}\cdot(\dot{\mathbf{r}}\times\mathbf{H})=\mathbf{r}\cdot(\mu\mathbf{u} + \mathbf{c}) = \mu\mathbf{r}\cdot\mathbf{u} + \mathbf{r}\cdot\mathbf{c} = \mu r(\mathbf{u}\cdot\mathbf{u})+rc\cos(\theta)=r(\mu + c\cos(\theta)) Burada \theta \bar{r} ve \bar{c} arasındaki açıdır. r’ye göre çözersek r = \frac{\mathbf{r}\cdot(\dot{\mathbf{r}}\times\mathbf{H})}{\mu + c\cos(\theta)} = \frac{(\mathbf{r}\times\dot{\mathbf{r}})\cdot\mathbf{H}}{\mu + c\cos(\theta)} = \frac{|\mathbf{H}|^2}{\mu + c\cos(\theta)} Dikkat ediniz ki (r,\theta) vektör fonksiyonunun polar koordinatlarıdır. Verilenleri yerine koyarsak ve, şu denkleme ulaşırız Bu fokal noktada merkezi olan bir konik bölmenin polar koordinatlarının denklemidir. \theta\, Argümanına gerçek anomali denir. === Yörünge denkleminin özellikleri === e\ =\ 0\, ise bu çapı p olan bir dairedir. 0\ < e\ <\ 1\, ise bu aşağıdakine sahip bir elipstir. = a \cdot \sqrt{1-e^2}|}} e\ =\ 1\, ise fokal uzunluğu \frac {p}{2} olan bir dairedir. e\ >\ 1\, ise bir aşağıdakine sahip bir hiperboldür. = a \cdot \sqrt{e^2-1}|}} Sonraki görsel, bir elipsi (kırmızı), bir paraboolü (yeşil) ve bir hiperbolü (mavi) gösterir. küçükresim|Dış merkezi 0.7 olan bir elips Kepler yörüngesi, Parabol Kepler yörüngesi ve dış merkezi 1.3 olan bir Kepler yörüngesi. Fokal noktaya olan uzaklık, denklemle de gösterildiği gibi polar açının yatay çizgiye yakın açısının bir fonksiyonudur() Yatay çizgi üzerinde nokta odak noktasından sağa doğru çıkıyorsa nokta \theta = 0\, odak noktasına olan minimum uzaklık \frac {p}{1 + e} değerini alır (pericentre) .Elips için de odak uzaklığı maksimum değerini aldığı bir apocentre vardır\frac {p}{1 - e}. Hiperbol için \theta\, ’nın aralığı: :\left [ -\cos^{-1}\left(-\frac{1}{e}\right) < \theta < \cos^{-1}\left(-\frac{1}{e}\right)\right ] Bir parabol için aralık; :\left [ -\pi < \theta < \pi \right ] Diferansiyelin zinci kuralını kullanarak (), () denklemi ve p’nin tanımı \frac {H^2}{\mu} olmak üzere radyal hız elemanı şudur. \cdot e \cdot \sin \theta|}} ve teğetsel bileşen olan (hız bileşeni V_r)’ye dik ) \cdot (1 + e \cdot \cos \theta)|}} Polar argümanı \theta\, ve zaman t arasındaki bağlantı elips ve hiperbolik yörüngeler için küçük değişiklikler gösterebilir. Elips şeklindeki bir yörünge için, dış merkezli anomali olan E’ye geçeriz ve burdan ve sonuç olarak. açısal momentum H olmak üzere Zamana göre entegrasyon olan t ile şu sonuca ulaşılır. zaman varsayımını yaparak integral sabiti 0 olur. p nin tanımı gereği şu şekilde yazılabiir ( E - e \cdot \sin E)|}} Hiperbolik bir yörünge için parametre yaratmak için hiperbolik fonksiyon kullanılır. bu durumda şuna ulaşılır açısal momentum H olmak üzere Zamana göre entegrasyon olan t ile şu sonuca ulaşılır. örneğin (e \cdot \sinh E-E)|}} t’nin hangi zamanda belirli bir gerçek anomaliye (\theta\,) ne ulaştığını bulmak için, buna uyan parametre olan E’nin () ile zamana olan bağını hesaplaması, elips için () ile bağını ve hiperbolik yörünge için. Dikkat ediniz ki () ve () bu aralıkta bir yer belirler. \left [ -\infin < t < \infin\right ] \longleftrightarrow \left [-\infin < E < \infin \right ] === Bazı ilave formüller === Ayrıca bakınız: Merkez denklemi - analitik açılımlar ()ve ()’den edilinen elips yörünge için ve buna göre ()’dan devam ettiğimizde şu sonuç çıkar : \tan^2 \frac{\theta}{2} = \frac{1-\cos \theta}{1+\cos \theta}= \frac{1-\frac{\cos E-e}{1-e \cdot \cos E}}{1+\frac{\cos E-e}{1-e \cdot \cos E}}= \frac{1-e \cdot \cos E - \cos E+e}{1-e \cdot \cos E + \cos E-e}= \frac{1+e}{1-e} \ \cdot\ \frac{1-\cos E}{1+\cos E}= \frac{1+e}{1-e} \ \cdot\ \tan^2 \frac{E}{2} Dış merkezli anomalinin geometrik yapısını belirlemek için, açıkça vektörler(\ \cos E\ ,\ \sin E\ ) ve (\ \cos \theta\ ,\ \sin \theta\ )‘nin x aksisiyle aynı yanda olmaları gerekir. Buradan yola çıkarak \left( \cos\frac{E}{2}\ ,\ \sin\frac{E}{2} \right) ve \left( \cos\frac{\theta}{2}\ ,\ \sin\frac{\theta}{2} \right) aynı çeyrek daireye sahiptir. Burdan yola çıkarak \cdot \tan \frac{E}{2}|}} ve (\ x\ ,\ y\ ) vektörünün polar argümanlarınınve n "\operatorname{arg}(x\ ,\ y)" örneğinde olduğu gibi seçilmiş olması gerekir. \operatorname{arg}(x\ ,\ y)nın sayısal hesaplaması için, standart ATAN2(y,x) fonksiyonunun, ( ya ikili DATAN2(y,x)), FORTRAN gibi programlama dillerinde kullanılması mümkündür. Dikkate alınız ki bu aralıklar arasında alan vardır. :\left [ -\infin < \theta < \infin\right ] \longleftrightarrow \left [-\infin < E < \infin \right ] Hiperbolik bir yörünge için () ve ()dan yola çıkarak and therefore that ve sonuç olarak : \tan^2 \frac{\theta}{2} = \frac{1-\cos\theta}{1+\cos \theta}= \frac{1-\frac{e-\cosh E}{e \cdot \cosh E-1}}{1+\frac{e-\cosh E}{e \cdot \cosh E-1}}= \frac{e \cdot \cosh E - e +\cosh E}{e \cdot \cosh E + e -\cosh E}= \frac{e+1}{e-1}\ \cdot\ \frac{\cosh E-1}{\cosh E+1}= \frac{e+1}{e-1}\ \cdot\ \tanh^2 \frac{E}{2} \tan \frac{\theta}{2} ve \tanh \frac{E}{2} aşağıdaki ifadeyi sağlayan aynı işarete sahiptir. \cdot \tanh \frac{E}{2}|}} Bu ilişki gerçek anomali ve E parametresi arasında geçiş yapmak için uygundur. (burada E ()’e zaman yönünden bağlıdır.) Aralık arasındaki alana dikkat ediniz. :\left [ -\cos^{-1}\left(-\frac{1}{e}\right) < \theta < \cos^{-1}\left(-\frac{1}{e}\right)\right ] \longleftrightarrow \left [-\infin < E < \infin \right ] ve aşağıdaki bağla açıklanabilir. :\tanh ^{-1}x=\frac{1}{2}\ln \left( \frac{1+x}{1-x} \right) ()’ye bakarak, P’nin elips bir yörünge için yörünge periyodu |}} Kuvvet alanına potansiyel enerjinin ilgisi () :-\frac {\mu} {r} (), (), () ve ()’ kadar kinetik ve potansiyel enerjinin özeti :\frac{{V_r}^2+{V_t}^2}{2}-\frac {\mu} {r} Elips yörünge için (), (), () ve ()’den sonra kinetik ve potansiyel enerjinin hiperbol yörünge için :\frac{{V_r}^2+{V_t}^2}{2}-\frac {\mu} {r} Eylemsizlik koordinat sistemiyle ilişkisi and from (), (), olan bir yörünge düzleminde pericentre’ye doğru () ve ()’dan yapılan hız elementleri çıkarımı === Belirli bir başlangıç durumuna karşılık Kepler yörüngesinin belirlenmesi === : \hat{x} \ , \ \hat{y} Bu diferansiyel denklemi () için "başlangıç değeri problemi" dir, ki 6 boyutlu bir durum vektründe ilk denklemdir ve şöyle yazılır \cdot \sin \theta|}} \cdot (e +\cos \theta)|}} Başlangıç değeri vektörünün herhangi bir değeri için, bu başlangıç değeri probleminin çözümüne uygun Kepler Yörüngesi aşağıdaki algoritmayla bulunur Ortogonal birim vektörleri (\hat{r}\ ,\ \hat{t}) tanımalamak gerekirse ve r > 0 and V_t > 0 (), (), () ve ()’dan sonra sırasını takip eder. # {{NumBlk|:|p = \frac{{(r \cdot V_t)}^2}{\mu }|}} e \ge 0 and \theta‘yı tanımlarken aşağıdakinin olduğu yerde # |}} aynı r, V_r and V_t değerlerine sahip olan gerçek anomaliye sahip bir Kepler yörüngesi bulunur ve bu () ve ()’de gösterilmiştir. Eğer bu Kepler yörüngesi aynı vektörlerine sahipse ve gerçek anomalisi () ve ()de tanımladığı gibi ise, Kepler yörüngesinin (\hat{r}\ ,\ \hat{t}) vektörleri istenen değer olan (\ \bar{r_0}\ ,\bar{v_0}\ )’yi gerçek anomalisi \theta olmak üzere alır. Standart atalet sabit koordinat sistemi (\hat{x}\ ,\ \hat{y}) yörüngesel düzlemde (\hat{x} pericentre'nin homojen küresindin merkezinden gelecek şekilde) konik alanın oryantasyonunu belirliyorsa (elips, parabol, hiperbol), bu ilişkiyle hesaplanabilir. # # () ve ()’ün V_r=0 olduğu durumda bağlantılı olduğuna dikkat ediniz ve :V_t=V_0=\sqrt{\frac{\mu}{p}}=\sqrt{\frac{\mu}{\frac{{(r \cdot V_t)}^2}{\mu }}} örneğin # |}} Diaresel bir yörüngeye uyan başlangıç durumunu (\ \bar{r_0}\ ,\bar{v_0}\ ) gösterir. === Oskülatör kepler yörüngesi === (\bar{r} , \bar{v} ) gibi herhangi bir durum vektörü için, bu durumdaki Kepler yörüngesi aşağıda tanılmlanan algoritma ile hesaplanabilir. İlk parametreler p , e , \theta ile belirlenir ve \hat{x} , \hat{y} yörüngesel düzelemde ortogonal birim vektörleri de () ve () arasındaki ilişkiyi kullanarak bulunur. Şimdi, eğer hareket denklemi aşağıdaki gibi ise # = \operatorname{\bar{F}}(\bar{r},\dot {\bar{r}},t)|}} aşağıdakinin :\operatorname{\bar{F}}(\bar{r},\dot {\bar{r}},t) bir fonksyion olup :-\mu \cdot \frac {\hat{r}} {r^2} aşağıdaki parametlere sahip olmaması gerekir. p ,\, e ,\, \theta,\, \hat{x} ,\, \hat{y} \bar{r},\dot {\bar{r}} tarafından tanımlanmış, zaman içinde değişecek olan durumlar, Keppler Yörüngesinin tersine yalnızca \theta’da değişecektir. T zamanında bu şekilde "hareket denkleminin" çözümü olarak aynı "durum vektör" olan () 'da hesaplanan Kepler yörüngesi bu durumda "oskülatör" olduğu söylenebilir. Bu konsept, şu durumda yararlıdır. :\operatorname{\bar{F}}(\bar{r},\dot {\bar{r}},t)=-\mu \cdot \frac {\hat{r}} {r^2}+\operatorname{\bar{f}}(\bar{r},\dot {\bar{r}},t) olmak üzere :\operatorname{\bar{f}}(\bar{r},\dot {\bar{r}},t) küçük bir “bozucu güç”tür çünkü, bir örnekle açıklamak gerekirse diğer gök cisimlerinden sönük bir yer çekimi bu durumda örnek gösterilebilir. Oskülasyon halindeki Kepler yörüngesinin parametreleri, ancak o zaman yavaşça değişecek ve oskülasyon halindeki Kepler yörüngesi iyi bir tahminle oskülasyonun öncesinde ve sonrasında gerçek yörüngeye hatrı sayılır bir zaman diliminde girecektir. Bu konsept aynı zamanda bir roket uçuşu için de faydalı olabilir çünkü itme kapatıldığı zaman hangi Keppler yörüngesinde devam edebileceği belirlenmiş olur. “Daire olmaya yakın” bir yörünge konsepti için, dışmerkezli vektörü şeklinde tanımlamak yararlıdır. (), () ve ()’dan yola çıkarak: {V_0}|}} Örneğin ( \bar{r} ,\bar{v} ) durum vektörlerinin diferansiyel fonskyonlarıdır \bar{e}, ayrıca bu durum dairesel bir yörüngeye karşılık geliyor ise de geçerlidir. == Ayrıca bakınız == * İki cisim problemi * Kütleçekimsel iki cisim problemi * Kepler problemi * Kepler'in gezegensel hareket yasaları * Eliptik yörünge * Hiperbolik yörünge * Kaçış yörüngesi ==Kaynakça== * El'Yasberg "Theory of flight of artificial earth satellites", Israel program for Scientific Translations (1967) * * Kategori:Yörüngeler Kategori:Johannes Kepler