Araştırma, 300 km başlangıç yüksekliğinden yeniden girişe başlayan bir uydu için atmosferik sürükleme, yerçekimi ve itki etkilerinin zamanla yüksekliği, hızı ve uçuş yolu açısını nasıl değiştirdiğini incelemektedir. Modelde 500 kg uydu kütlesi, 1 m² kesit alanı, 2,2 sürükleme katsayısı, 1000 N sabit itki ve 8500 m atmosferik ölçek yüksekliği kullanılmış; diferansiyel denklemler 0–2000 s aralığında MATLAB ode45 ile çözülmüştür. Simülasyon, irtifa azaldıkça üstel olarak artan atmosfer yoğunluğunun sürükleme kuvvetini keskin biçimde yükselttiğini ve alt atmosferde güçlü yavaşlamanın baskın mekanizmasına dönüştüğünü göstermektedir. Bununla birlikte model iki boyutlu, sabit sürükleme katsayılı ve termal etkileri içermeyen sadeleştirilmiş bir simülasyondur; ayrıca yüzeyde sonlandırma koşulu bulunmadığı için sayısal çözüm Dünya yüzeyinin yaklaşık 33,6 km altına kadar devam etmektedir.

Kaynakta son hesaplanan durum yaklaşık −33.639,09 m irtifa, 9,20 m/s hız ve −90,00° uçuş yolu açısı olarak verilmektedir. Ancak negatif irtifa fiziksel bir yeniden giriş yüksekliği değildir; makalenin kendi açıklamasına göre ode45 entegrasyonunun yüzey sınırından sonra devam etmesinden kaynaklanan sayısal aşımı temsil eder.

Araştırmanın temel problemi nedir?

Çalışmanın çıkış noktası, uydu yeniden girişinde farklı fiziksel kuvvetlerin aynı anda etkili olmasıdır. Uydu irtifa kaybettikçe yerçekimi, atmosferik yoğunluk ve aerodinamik sürükleme değişmekte; varsa itki kuvveti de hareket üzerinde ek kontrol olanağı sağlamaktadır. Bu etkileşimler doğrusal olmayan diferansiyel denklemlerle ifade edildiğinden, yeniden giriş hareketinin zaman içinde sayısal olarak çözülmesi gerekir.

Yazarlar, gelişmiş yeniden giriş ve enkaz değerlendirme yazılımlarına kıyasla daha erişilebilir bir MATLAB çerçevesi oluşturmaya çalışmaktadır. Literatür bölümünde yörünge mekaniği, aerodinamik sürükleme, yelken ve tether tabanlı de-orbit teknikleri, güneş etkinliğine bağlı atmosfer değişimi, yeniden giriş yönlendirmesi ve mevcut enkaz değerlendirme yazılımları tartışılmakta; önerilen MATLAB yaklaşımı özellikle eğitim ve erken aşama görev incelemeleri için daha sade bir araç olarak konumlandırılmaktadır.

Modelde kullanılan fiziksel parametreler

Hareket denklemleri nasıl kurulmuştur?

Kaynak, yeniden giriş hareketini öncelikle yükseklik h, hız v ve uçuş yolu açısı γ üzerinden tanımlar. Birinci denklem yükseklik değişim hızını verir:

\[ \frac{dh}{dt}=v\sin\gamma \]

Burada h uydu yüksekliğini, v hız büyüklüğünü ve γ uçuş yolu açısını temsil eder. Negatif bir uçuş yolu açısı, modelin işaret konvansiyonunda alçalma hareketine karşılık gelmektedir.

Kaynakta verilen hız denklemi şöyledir:

\[ \frac{dv}{dt}=\frac{F_{\mathrm{thrust}}-D}{m} \]

F_thrust uygulanan itki kuvvetini, D aerodinamik sürükleme kuvvetini ve m uydu kütlesini gösterir. Önemli bir kaynak içi ayrıntı, bu gösterilen denklemde yerçekiminin hız doğrultusundaki bir bileşeninin açıkça bulunmamasıdır. Buna rağmen sonuç tartışmasında net ivmenin itki, sürükleme ve yerçekiminin birleşik etkisini içerdiği belirtilmektedir. Makale bu iki gösterimi tam olarak uzlaştıran ek bir denklem vermemektedir.

Uçuş yolu açısının evrimi için verilen denklem:

\[ \frac{d\gamma}{dt}=-\frac{g}{v}\cos\gamma \]

Kontrol girdisi eklenmiş biçim ise:

\[ \frac{d\gamma}{dt}=-\frac{g}{v}\cos\gamma+u(t) \]

Buradaki u(t), yazarların manuel veya otomatik olabileceğini belirttiği kontrol fonksiyonudur. Ancak çalışma belirli bir u(t) kontrol yasası, optimizasyon amaç fonksiyonu veya optimizasyon algoritması tanımlayıp sonuçlarını raporlamamaktadır. Bu nedenle başlıktaki “optimizasyon” ifadesine karşın sunulan hesaplamanın ana gövdesi bir yeniden giriş simülasyonudur.

Aerodinamik sürükleme nasıl hesaplanmaktadır?

Kaynakta aerodinamik sürükleme kuvveti klasik karesel hız bağımlılığıyla ifade edilmektedir:

\[ D=\frac{1}{2}\rho(h)v^2 C_d A \]

Burada ρ(h) yüksekliğe bağlı atmosfer yoğunluğunu, v uydu hızını, C_d sabit alınan sürükleme katsayısını ve A referans kesit alanını ifade eder. Bu yapı nedeniyle yoğunluk ve hız, sürükleme kuvvetinin büyüklüğünü birlikte belirlemektedir.

Yerçekimi yüksekliğe göre nasıl değişmektedir?

Yerçekimi ivmesi Dünya merkezinden uzaklığa bağlı ters-kare ilişkisiyle hesaplanmaktadır:

\[ g=\frac{GM}{(R_{\mathrm{earth}}+h)^2} \]

G evrensel kütleçekim sabiti, M Dünya kütlesi, R_earth Dünya yarıçapı ve h yüzey üzerindeki yüksekliktir. Şekil 6'da uydu Dünya'ya yaklaştıkça yerçekimi kuvvetinin nispeten küçük fakat düzenli bir artış gösterdiği görülmektedir.

Atmosfer yoğunluğu modeli

Çalışma atmosfer yoğunluğunu tek ölçek yüksekliğine sahip üstel bir fonksiyonla temsil etmektedir:

\[ \rho(h)=\rho_0\exp\left(-\frac{h}{H}\right) \]

Burada \(\rho_0=1.225\,\mathrm{kg/m^3}\) deniz seviyesi yoğunluğu, \(H=8500\,\mathrm{m}\) ise ölçek yüksekliğidir. Bu denklem modelin aerodinamik sürüklemeyi irtifaya bağlayan temel mekanizmasıdır. İrtifa azaldıkça \(\rho(h)\) hızla yükselmekte ve sürükleme kuvveti büyümektedir.

Ancak çalışma bu tek üstel atmosfer modelini 300 km başlangıç yüksekliğine kadar uygulamaktadır. Kaynakta zamana bağlı güneş etkinliği, jeomanyetik indisler, atmosferik bileşim değişimi veya çok katmanlı atmosfer modeli kullanılmamaktadır. Dolayısıyla “yüksekliğe bağlı atmosfer yoğunluğu” içerilmekle birlikte atmosferin gerçek değişkenliğinin tamamı modellenmemektedir.

Başlangıç ve çözüm koşulları

• Başlangıç yüksekliği: 300.000 m.
• Başlangıç hızı: Kaynak metinde “7.500 m” biçiminde yazılmıştır; Şekil 3 hız ekseni ise m/s kullanmaktadır.
• Başlangıç uçuş yolu açısı: −5°, radyana dönüştürülmüştür.
• Çözüm aralığı: 0–2000 s.
• Sayısal çözücü: MATLAB ode45.

Yazarlar −5° değerini, kaynakta belirtilen yaklaşık −1° ile −10° aralığı içerisinde “orta derecede sığ” bir yeniden giriş açısı olarak kullanmaktadır. Ancak bu seçimin ayrı bir optimizasyon analiziyle elde edildiği gösterilmemiştir.

Şekil 1: Simülasyon iş akışı

Şekil 1, yöntemi beş temel adım halinde göstermektedir: başlangıç, iki boyutlu yörünge mekaniği modelinin kurulması, uydunun başlangıç dairesel yörüngesine yerleştirilmesi, bozucu kuvvetlerin eklenmesi ve ode45 ile sayısal çözüm. Bu akış, çalışmanın esasının diferansiyel denklem tabanlı zaman entegrasyonu olduğunu açık biçimde göstermektedir.

Verianla Live: Görselleştirme, bu makaledeki görünür bilimsel veri tablosundan tarayıcıda oluşturulur. Tablo bilimsel kaynak-of-truth olarak korunur.

Şekil 2: İrtifa zamanla nasıl değişmektedir?

Şekil 2'de uydu yaklaşık 300 km başlangıç irtifasından aşağı doğru ilerlemektedir. Eğrinin ilk bölümündeki hızlı irtifa kaybını daha sonra yüzey referansının altına uzanan sayısal çözüm izlemektedir. Makalenin daha sonraki metni, son irtifanın −33.639,09 m olmasını fiziksel bir yeraltı hareketi olarak değil, entegratörün \(h=0\) sınırında durdurulmaması nedeniyle oluşan sayısal aşım olarak açıklamaktadır.

Şekil 3: Hız profilinin ana mesajı

Şekil 3'te hız başlangıçta yüksek yörüngesel seviyedeyken, yeniden girişin erken bölümünde çok keskin biçimde düşmektedir. Yazarlar bu değişimi atmosfer yoğunluğunun artmasıyla güçlenen aerodinamik sürüklemeye bağlamaktadır. Grafiğin dikey ekseni m/s'dir. Bu nokta önemlidir; çünkü yöntem metnindeki başlangıç hızı satırında birim “m” olarak yazılmıştır.

Şekil 4: Uçuş yolu açısı

Başlangıçtaki −5° uçuş yolu açısı zamanla daha negatif değerlere ilerlemekte ve simülasyon sonunda −90°'ye ulaşmaktadır. Kaynak bunu hareketin giderek daha dik ve sonunda neredeyse düşey bir alçalma durumuna geçmesi olarak yorumlamaktadır.

Şekil 5: Sürükleme kuvveti neden kısa süreli büyük bir tepe oluşturuyor?

Şekil 5, sürükleme kuvvetinin yeniden girişin belirli bir bölümünde keskin bir tepe oluşturduğunu ve ardından hızın güçlü biçimde azalmasıyla yeniden düştüğünü göstermektedir. Kaynağın temel yorumu, uydu daha yoğun atmosfere girdiğinde aerodinamik sürüklemenin ana yavaşlatıcı etkiye dönüşmesidir. Makale grafiğin ham veri noktalarını tablo halinde vermediğinden, yalnız görselden okunabilecek yaklaşık değerleri kesin veri gibi yeniden üretmek uygun değildir.

Şekil 6: Yerçekimi kuvveti

Yerçekimi kuvveti sürükleme grafiğine kıyasla çok daha yumuşak değişmektedir. Şekilde yaklaşık 4,4–5,0 kN bandında yukarı yönlü bir eğilim görülür; ancak kaynak metninin bilimsel olarak güvenli ana sonucu yalnız uydu irtifa kaybettikçe ters-kare bağıntısı nedeniyle yerçekiminin bir miktar artmasıdır.

Şekil 7: Net ivme

Net ivme grafiğinde yeniden girişin yoğun aerodinamik bölümüne karşılık gelen güçlü bir negatif tepe vardır. Yazarlar bu durumu sürükleme etkisinin baskın hale gelmesiyle açıklamaktadır. Bununla birlikte makalenin gösterdiği \(dv/dt\) denkleminde yerçekimi açık bir terim olarak bulunmadığı için, Şekil 7'nin “itki + sürükleme + yerçekimi” birleşimini hangi tam matematiksel formülle kullandığı kaynakta yeterince açık değildir.

Şekil 8: Son yeniden giriş koşulları

Başlıktaki “optimizasyon” ne ölçüde gerçekleştirilmiştir?

Makale, uçuş yolu açısına \(u(t)\) kontrol girdisi eklenebileceğini göstermekte ve itkinin yeniden giriş karakteristiklerinin kontrolünde kullanılabileceğini tartışmaktadır. Bununla birlikte kaynakta optimize edilen açık bir amaç fonksiyonu, karar değişkeni, kısıt, optimizasyon algoritması veya farklı çözümler arasından seçilmiş optimum sonuç bulunmamaktadır. Dolayısıyla çalışmadan doğrudan desteklenen sonuç, bir MATLAB yeniden giriş simülasyon çerçevesi oluşturulduğudur; formal yörünge optimizasyonunun gerçekleştirildiği sonucu daha güçlü bir iddia olur.

Türkiye açısından ne ifade edebilir?

Çalışmanın verileri Türkiye'ye ait değildir ve model Türkiye'ye özgü bir uzay görevi üzerinde doğrulanmamıştır. Bununla birlikte sade MATLAB yapısı; yörünge mekaniği, yeniden giriş dinamiği ve sayısal diferansiyel denklem çözümünün eğitim veya erken aşama mühendislik çalışmaları açısından anlaşılmasına yardımcı olabilecek bir örnek sunmaktadır. Türkiye'de gerçek bir uydu yeniden giriş planlamasına uygulanması istenirse atmosfer modeli, araç geometrisi, aerodinamik katsayılar, termal yükler, üç boyutlu yörünge, kontrol sistemi, gerçek görev verileri ve güvenlik sınırları gibi unsurların ayrıca doğrulanması gerekir.

Çalışma ne söylüyor?

• Verilen parametre setiyle MATLAB ode45 kullanılarak sadeleştirilmiş bir yeniden giriş simülasyonu kurulabileceğini gösteriyor.
• Model içinde irtifa azaldıkça atmosfer yoğunluğu ve aerodinamik sürükleme etkisinin belirginleştiğini gösteriyor.
• Simülasyonda hızın güçlü biçimde azaldığını ve uçuş yolu açısının −90°'ye doğru dikleştiğini gösteriyor.
• Kaynak, alt atmosferde sürüklemeyi temel yavaşlatıcı mekanizma olarak yorumluyor.
• Yüzey sonlandırma koşulu bulunmadığında sayısal çözümün fiziksel sınırı aşabileceğini kendi sonuçlarıyla gösteriyor.

Çalışma neyi göstermiyor?

• Gerçek bir uydunun yeniden giriş uçuş testini veya saha doğrulamasını göstermiyor.
• Üç boyutlu yeniden giriş yörüngesi sunmuyor.
• Atmosferik kaldırmayı modellemiyor.
• Aerodinamik katsayıların Mach sayısı, sıcaklık veya yönelime göre değişimini modellemiyor.
• Aerodinamik ısınmayı veya termal koruma sistemi yüklerini hesaplamıyor.
• Parçalanma, ablation veya enkazın hayatta kalma analizini içermiyor.
• Belirsizlik veya hassasiyet analizi vermiyor.
• Başlığa rağmen formal bir optimizasyon problemi çözüp optimum yörünge raporlamıyor.
• −33,6 km son irtifayı fiziksel yeniden giriş konumu olarak doğrulamıyor; bu değer sayısal aşım olarak açıklanıyor.

Çalışmanın Yöntemi ve Bulguları

Sayısal yöntem

Çalışma, yeniden giriş hareketini adi diferansiyel denklemler sistemi biçiminde kurmuş ve MATLAB'ın değişken adımlı sayısal ODE çözücüsü ode45 ile 0–2000 saniye boyunca çözmüştür. Kaynakta tolerans ayarları, maksimum adım boyutu, olay fonksiyonu veya hata kontrol parametreleri ayrıca bildirilmemiştir. Özellikle \(h=0\) için bir sonlandırma olayı tanımlanmaması, çözücünün Dünya yüzeyini geçerek negatif irtifaya devam etmesine yol açmıştır.

Model zinciri

1. 300 km başlangıç irtifası ve başlangıç durumu tanımlanmıştır.
2. Yerçekimi irtifaya bağlı ters-kare bağıntısıyla hesaplanmıştır.
3. Atmosfer yoğunluğu üstel \(\rho(h)\) modeliyle hesaplanmıştır.
4. Sürükleme kuvveti \(D=\frac12\rho v^2 C_dA\) ile hesaplanmıştır.
5. Hareket denklemleri ode45 ile zamana göre bütünleştirilmiştir.
6. İrtifa, hız, uçuş yolu açısı, sürükleme kuvveti, yerçekimi kuvveti ve net ivme profilleri grafiklenmiştir.
7. Simülasyon sonunda son irtifa, hız ve uçuş yolu açısı raporlanmıştır.

Ana nicel koşullar ve bulgular

Modelin güçlü yönleri

Çalışmanın başlıca güçlü yanı, yeniden giriş fiziğinin temel bileşenlerini az sayıda açık denklem ve MATLAB'ın yaygın bir sayısal çözücüsüyle bir araya getirmesidir. İrtifa, hız, uçuş yolu açısı ve kuvvet profillerinin ayrı grafiklerle gösterilmesi, model içindeki neden-sonuç zincirinin eğitim açısından izlenmesini kolaylaştırmaktadır. Ayrıca yazarların negatif son irtifayı gizlemek yerine sayısal aşım olarak açıkça tartışması, çözücü sınır koşullarının önemini görünür hale getirmektedir.

Temel sınırlılıklar

• Yörünge iki boyutlu düzlemsel modelle sınırlandırılmıştır.
• Sürükleme katsayısı sabittir.
• Kaldırma kuvveti ihmal edilmiştir.
• Termal ısınma modellenmemiştir.
• Üç boyutlu yönelim ve atmosferik uçuş dinamiği yoktur.
• Aktif kontrol stratejisi sonuçlarla doğrulanmamıştır.
• Atmosfer tek üstel yoğunluk modeliyle temsil edilmiştir.
• Ham simülasyon verisi veya kod listesi makalede verilmemiştir.
• Yüzey temasını durduran olay koşulu kullanılmamıştır.
• Formal optimizasyon prosedürü tanımlanmamıştır.

Bu nedenle çalışma, ayrıntılı görev sertifikasyonu veya gerçek yeniden giriş güvenliği tahmini yerine, sadeleştirilmiş sayısal modelleme ve eğitim/erken tasarım bağlamında değerlendirilmelidir.

Kaynak ve Yöntem Notu

Özgün başlık: Simulation and Optimization of Satellite Re-Entry Trajectories Using MATLAB

Yazarlar: Mohamed Shuaib.A; Dr. S. Charulatha

Sorumlu yazar: Mohamed Shuaib.A

Kurum: Department of Aerospace, Hindustan Institute of Technology and Science, Chennai, Tamilnadu, India

Dergi: International Journal of Advanced Engineering Research and Science (IJAERS)

Cilt / Sayı: 13 / 3

Sayfalar: 15–21

Yayın tarihi: Mart 2026; çevrimiçi erişim tarihi kaynakta 17 Mart 2026 olarak verilmiştir.

DOI: 10.22161/ijaers.133.2

Yayınevi: AI Publication

Kaynak türü: Hakemli dergide yayımlanan MATLAB tabanlı sayısal simülasyon/modelleme çalışması.

Lisans: Creative Commons Attribution 4.0 International (CC BY 4.0).

Finansman: İncelenen sürümde ayrı finansman beyanı bulunmamaktadır.

Veri erişilebilirliği: İncelenen sürümde ayrı veri erişilebilirliği beyanı bulunmamaktadır.

Çıkar çatışması: İncelenen sürümde ayrı çıkar çatışması beyanı bulunmamaktadır.

Yazar katkıları: İncelenen sürümde ayrı CRediT/yazar katkı tablosu bulunmamaktadır.

Resmî DOI bağlantısı: https://doi.org/10.22161/ijaers.133.2

Bilimsel anlatım yüklenen çalışmanın metnine, denklemlerine, Tablo 1'e ve Şekil 1–8'e dayanmaktadır. Dış kaynaklardan yeni deneysel sonuç veya yeni bilimsel bulgu eklenmemiştir. Çalışma gerçek bir yeniden giriş uçuş deneyi değil, sadeleştirilmiş iki boyutlu bir sayısal simülasyondur. Sabit sürükleme katsayısı, basitleştirilmiş atmosfer modeli, kaldırma ve termal etkilerin ihmal edilmesi, üç boyutlu dinamiklerin bulunmaması ve yüzey sonlandırma koşulunun kullanılmaması sonuçların yorum sınırını belirlemektedir.