Bu çalışma, gaz türbini motorlarının mimarisini, termodinamik çalışma koşullarını ve toplam ağırlığını aynı matematiksel problem içinde birlikte optimize etmektedir. Araştırmacılar; ara soğutucu, ısı geri kazanım sistemi, döner detonasyonlu yeniden ısıtma bölümü ve üçüncü akış kanalı gibi isteğe bağlı bileşenleri ikili karar değişkenleriyle temsil eden bir karışık tamsayılı doğrusal olmayan programlama modeli geliştirmiştir. Model, 10, 30 ve 50 kN seyir itkisi sınıflarında küresel verimlilik ile güç sistemi ağırlığı arasındaki dengeyi Pareto cepheleri üzerinden araştırmıştır. Önerilen yeni motor mimarileri 30 ve 50 kN sınıflarında optimize edilmiş geleneksel motorlara göre yaklaşık %15 daha yüksek küresel verimlilik ve %18-22 daha düşük ağırlık üretmiştir. Bununla birlikte bulgular, hakem değerlendirmesinden geçmemiş, seyir tasarım noktasına odaklanan ve gerçek motor ya da uçuş deneyi içermeyen bir ön tasarım modeline dayanmaktadır.

Araştırmanın ikinci katkısı, fiziksel olarak mümkün olmayan veya verimsiz mimari dallarını erken aşamada eleyen Çok Başlangıçlı Fizik Bilgili Dal ve Sınır algoritmasıdır (Multi-Start Physics-Acquainted Branch and Bound, MS-PABnB). Tam faktöriyel örnekleme, AdaBoost fizibilite sınıflandırıcısı, LSBoost sıcaklık tahmini ve yakıt-hava oranı belirleme aşamalarının birlikte kullanılması, uygulanabilir başlangıç noktalarının oranını %5,83'ten %88,6'ya yükseltmiştir. MS-PABnB, aynı çözüm ortamındaki geleneksel Dal ve Sınır yaklaşımından genel olarak daha az düğüm araştırmış, daha az başlangıç noktasına ihtiyaç duymuş ve daha iyi ölçeklendirilmiş amaç fonksiyonu değerleri üretmiştir.

Seçilen seyir görevlerinde model, 30 ve 50 kN sınıfları için referans motorlara kıyasla yakıt tüketimi ve karbondioksit salımında yaklaşık %28'e varan azalmalar hesaplamıştır. Ancak bu sonuçlar doğrudan ölçülmüş uçuş sonuçları değildir. Görev hesabı, referans motorların itki profilleri, hesaplanan özgül yakıt tüketimleri ve ampirik emisyon ilişkileri kullanılarak yapılmıştır. Özellikle 10 kN tasarımında yakıt tüketimi azalmasına rağmen azot oksit salımı referans motora göre %86,86 daha yüksek hesaplanmıştır.

Araştırmanın temel sorusu nedir?

Araştırmanın temel sorusu, bir uçak gaz türbini motorunda yalnızca basınç oranı, hava debisi veya yanma sıcaklığı gibi sürekli çalışma değişkenlerinin değil, motorun hangi bileşenlerden oluşacağının da aynı anda optimize edilip edilemeyeceğidir. Geleneksel tasarım yaklaşımında motor mimarisi çoğunlukla önceden belirlenir; ardından bu sabit mimari içindeki termodinamik değişkenler ayarlanır. Böyle bir sıra, henüz başlangıçta seçilmeyen bileşen kombinasyonlarının değerlendirme dışında kalmasına neden olur.

Çalışma bu sınırlamayı, bileşen kurulumlarını sıfır veya bir değerini alan ikili değişkenlerle tanımlayarak aşmayı amaçlamaktadır. Optimizasyon böylece yalnızca “Bu motor hangi basınç oranında çalışmalı?” sorusunu değil, “Ara soğutucu kurulmalı mı?”, “Isı geri kazanımı kullanılmalı mı?”, “Çekirdek akışın bir bölümü döner detonasyon odasında yeniden ısıtılmalı mı?” ve “Üçüncü bir hava akışı bulunmalı mı?” sorularını da çözmektedir.

Çalışma neden önemlidir?

Motor verimliliğini artıran her bileşen uçuş sırasında net kazanç sağlamaz. Bir ısı eşanjörü egzoz enerjisinin bir bölümünü geri kazanabilir; ancak eşanjörün ağırlığı ve oluşturduğu basınç kaybı bu kazancı azaltabilir. Daha büyük bir fan itki verimliliğini yükseltebilir; fakat fan çapı, motor kaportası, aktarma sistemi, ağırlık ve hava aracı entegrasyonu üzerinde yeni kısıtlar oluşturur. Bu nedenle yalnızca termodinamik verimin en yüksek olduğu çevrim, uçak açısından en uygun motor olmak zorunda değildir.

Araştırmacılar bu çatışmayı iki amacı aynı anda değerlendirerek ele almaktadır:

• Küresel motor verimliliğinin mümkün olduğunca yükseltilmesi,
• Fan, kompresörler, yanma odası, türbinler, eşanjörler, karıştırma odası, yapısal elemanlar ve yardımcı sistemler dâhil toplam güç sistemi ağırlığının mümkün olduğunca düşürülmesi.

Bu yaklaşım tek bir “en iyi motor” yerine, verimlilik artışı karşılığında ne kadar ağırlık kabul edildiğini gösteren bir Pareto cephesi üretmektedir. Pareto cephesindeki her nokta, bir hedefin daha fazla iyileştirilmesinin diğer hedefte kayıp oluşturduğu bir tasarımı temsil eder.

Önerilen motor mimarisi nasıl çalışmaktadır?

Çalışmadaki motor akış şeması, klasik iki şaftlı turbofan yapısını dört isteğe bağlı unsurla genişletmektedir:

• Ara soğutucu: Düşük basınç ve yüksek basınç kompresörleri arasındaki havayı soğutarak sonraki sıkıştırma işi üzerinde etkili olur.
• Isı geri kazanım sistemi: Çekirdek egzozundaki ısının daha soğuk bir akışa aktarılmasını sağlar.
• Döner detonasyonlu yanma bölümü: Yüksek basınç türbininden sonraki çekirdek akışın belirli bir bölümünü yeniden ısıtarak basınç kazançlı yanmayı yaklaşık olarak temsil eder.
• Üçüncü akış kanalı: Fan sonrasında ayrılan havanın bir bölümünü ısı eşanjörlerinden geçirip ikincil akışla yeniden birleştirir.

Motor şemasında ana çekirdek akış, fan akışı ve üçüncü akış ayrı yollarla gösterilmektedir. Düşük basınç kompresörü ve yüksek basınç kompresörü sayısal olarak dört sanal bölüme ayrılmıştır. Bu sanal bölümler gerçek kompresör kademeleri değildir; gaz özelliklerini kademeli biçimde güncellemek ve türbin soğutma havasının hangi kompresör noktasından alınacağını seçmek için kullanılan hesaplama düğümleridir.

Döner detonasyon bölümü etkin olduğunda çekirdek akışın tamamı değil, optimizasyon tarafından belirlenen bir bölümü yeniden ısıtılır. Yeniden ısıtılmış akış, basınç salınımlarının düzenlenmesi amacıyla bir karıştırma odasından geçirilir ve düşük basınç türbini bölümlerine yönlendirilir. Model aynı zamanda türbin soğutmasının gerekli olup olmadığını ve soğutma havasının yüksek basınç kompresörünün hangi sanal kademesinden çekileceğini belirlemektedir.

Matematiksel modelin ana yapısı nedir?

Model iki sıkı bağlı alt sistemden oluşmaktadır. Termodinamik modül her istasyondaki kütle, enerji ve momentum dengelerini çözerek itkiyi, sıcaklıkları, basınçları, akış hızlarını ve verimleri hesaplar. Ağırlık modülü ise ortaya çıkan debi, çap, kademe sayısı, yüzey alanı ve malzeme parametrelerinden bileşen ağırlıklarını tahmin eder. İki modül ardışık biçimde değil, aynı optimizasyon problemi içinde eş zamanlı çözülür.

Bir motor istasyonundaki kütlesel debinin genel biçimi şu ilişkiyle verilmektedir:

\[ \dot{m}_s = \dot{m}_{\pi}\left(FM_s-\sum_{i\in\mathcal{I}}\sum_{s'\in\mathcal{X}_s}SB_{s',i}+\sum_{i\in\mathcal{I}_s}GB_i+\sum_{c\in\mathcal{Z}_s}FAR_c\right)RF_s \]

Burada \(\dot{m}_s\), ilgili motor istasyonundaki kütlesel debiyi; \(\dot{m}_{\pi}\), çekirdek hava debisini; \(FM_s\), çekirdek, bypass veya üçüncü akışa bağlı debi çarpanını; \(SB_{s',i}\), kompresörden çekilen soğutma havasını; \(GB_i\), türbine yeniden verilen soğutma havasını; \(FAR_c\), yakıt-hava oranını ve \(RF_s\), yeniden ısıtılan akış oranını temsil eder. Denklem, motor boyunca hava çekilmesi, yakıt eklenmesi ve soğutma akışının yeniden katılması nedeniyle debinin sabit kalmadığını ifade etmektedir.

Ardışık motor istasyonları arasındaki toplam sıcaklık ve toplam basınç değişimleri iki genel oranla tanımlanmaktadır:

\[ T_{0s}=T_{0s'}TR_{c,s'} \]

\[ P_{0s}=P_{0s'}PR_{c,s'} \]

\(TR\), bileşenin sıcaklık oranını; \(PR\), basınç oranını gösterir. Kompresör, türbin, yanma odası, eşanjör, kanal veya karışım işlemi için uygun fiziksel denklem bu iki orana bağlanmaktadır. Böylece motorun tamamı istasyondan istasyona izlenebilen ortak bir yapıda modellenmektedir.

Döner detonasyonlu yeniden ısıtma nasıl modellenmiştir?

Döner detonasyon odasının basınç kazancı, idealize edilmiş sabit hacimli ısı ekleme sürecine benzetilmektedir. Kullanılan temel basınç oranı şu biçimdedir:

\[ PR_{RDC,\pi42}=\frac{T^{CV}_{0\pi43}}{T_{0\pi42}}\left(\frac{T_{0\pi43}}{T^{CV}_{0\pi43}}\right)^{\frac{\gamma_{\pi43}}{\gamma_{\pi43}-1}} \]

\(PR_{RDC,\pi42}\), döner detonasyon odasının tahmini basınç oranını; \(T_{0\pi42}\), odaya giriş toplam sıcaklığını; \(T_{0\pi43}\), gerçek çıkış toplam sıcaklığını; \(T^{CV}_{0\pi43}\), eşdeğer sabit hacimli süreç sıcaklığını ve \(\gamma\), özgül ısı oranını ifade eder. Model gerçek detonasyon dalgalarını zamana bağlı üç boyutlu akış çözümüyle hesaplamamaktadır. Bunun yerine detonasyonlu yanmanın çevrim seviyesindeki olası basınç ve sıcaklık etkisini yaklaşık bir termodinamik bağıntıyla temsil etmektedir.

Net itki ve motor verimi nasıl hesaplanmıştır?

Net itki, çekirdek ve ikincil nozullardaki momentum değişimleri ile nozul çıkış basınçlarının ortam basıncından farklı olmasının oluşturduğu basınç itkisinin toplamıdır:

\[ F_{Net}=\dot{m}_{\pi}\left[V_{\pi9}\left(1+\sum_{c\in\mathcal{Z}}FAR_c\right)-V_{\infty}\right]+A_{\pi9}(P_{\pi9}-P_{\infty})+\dot{m}_{\pi}(\beta+\sigma)(V_{\beta9}-V_{\infty})+A_{\beta9}(P_{\beta9}-P_{\infty}) \]

\(V_{\pi9}\) ve \(V_{\beta9}\), çekirdek ve ikincil nozul çıkış hızlarını; \(V_{\infty}\), uçuş hızını; \(\beta\), bypass oranını; \(\sigma\), üçüncü akış oranını; \(A\), nozul çıkış alanını ve \(P\), basınçları göstermektedir. Optimizasyonda net itki, hedeflenen 10, 30 veya 50 kN değerine yaklaşık \(10^{-12}\) mertebesindeki toleransla eşitlenmektedir.

Küresel verimlilik, uçağa aktarılan yararlı itki gücünün yakıtla sağlanan enerjiye oranıdır:

\[ \eta_{glob}=\frac{F_{Net}V_{\infty}}{HV\dot{m}_{\pi}\sum_{c\in\mathcal{Z}}FAR_c}=\eta_{th}\eta_{prop} \]

\(\eta_{glob}\), küresel verimliliği; \(HV\), yakıtın alt ısıl değerini; \(\eta_{th}\), ısıl verimliliği ve \(\eta_{prop}\), itki verimliliğini temsil eder. Çalışmada Jet-A1 benzeri yakıt için ısıl değer \(44\times10^6\) J/kg olarak alınmıştır. Eşitlik, yakıt enerjisinin önce termodinamik çevrimde mekanik ve jet enerjisine, ardından bu jet enerjisinin yararlı itkiye dönüştürüldüğünü ifade eder.

Ağırlık modeli neleri kapsamaktadır?

Toplam güç sistemi ağırlığı şu şekilde birleştirilmiştir:

\[ W_{PwP}=W_{add}+W_{struct}+\sum_{c\in\mathcal{C}}W_c \]

\(W_c\), ayrı motor bileşenlerinin ağırlığını; \(W_{struct}\), yapısal ağırlığı; \(W_{add}\), kablolar, bağlantı elemanları ve doğrudan modellenmeyen yardımcı sistemleri ifade eder. Fan, kompresör, yanma odası ve türbin ağırlıkları çap, debi, eksenel hız, kademe sayısı ve motor gücü gibi değişkenlerden türetilmiştir.

Isı eşanjörü ağırlığı yüzey alanına bağlı olarak hesaplanmaktadır:

\[ W_c=2K_{hex}S_c\lambda_c \]

\(K_{hex}=3,6\) kg/m², eşanjör çekirdeğinin alan yoğunluğunu; \(S_c\), eşanjör alanını; \(\lambda_c\), eşanjörün kurulup kurulmadığını gösteren ikili değişkeni temsil eder. İki katsayısı; kapak, kanal ve bağlantı borularının eşanjör çekirdeğine yaklaşık eşit ek kütle oluşturduğu varsayımından gelmektedir.

Döner detonasyon odasının ağırlığı ise yeniden ısıtılan akış oranıyla ana yanma odası ağırlığının doğrusal ölçeklenmesi üzerinden tahmin edilmiştir:

\[ W_{RDC}=\chi W_{CC} \]

\(\chi\), yeniden ısıtılan çekirdek akış oranını; \(W_{CC}\), ana yanma odasının tahmini ağırlığını göstermektedir. Bu basit ilişki, olgunlaşmış bir döner detonasyon motoru ağırlık veri tabanı bulunmadığı için kullanılan ön tasarım varsayımıdır ve sonuçların önemli belirsizliklerinden biridir.

Verimlilik ile ağırlık arasındaki denge nasıl kurulmuştur?

Verimlilik ve ağırlık farklı birimlere ve sayısal ölçeklere sahip olduğu için önce tek amaçlı optimizasyonlarla iki referans değer belirlenmiştir: ulaşılabilen en yüksek küresel verimlilik ve ulaşılabilen en düşük güç sistemi ağırlığı. Daha sonra iki amaç normalize edilmiştir:

\[ \Theta_{Eff}=\frac{\eta^{Anch}_{glob}-\eta_{glob}}{\eta^{Anch}_{glob}} \]

\[ \Theta_{We}=\frac{W_{PwP}-W^{Anch}_{PwP}}{W^{Anch}_{PwP}} \]

Birleştirilmiş amaç fonksiyonu şu biçimdedir:

\[ \Theta_{Opt}=\alpha_{We}\Theta_{We}+\alpha_{Eff}\Theta_{Eff} \]

\[ \alpha_{We}+\alpha_{Eff}=1 \]

\(\alpha_{Eff}\) büyüdükçe optimizasyon verimliliğe, \(\alpha_{We}\) büyüdükçe düşük ağırlığa daha fazla öncelik verir. Araştırmacılar katsayıları 0 ile 1 arasında 0,025 aralıklarla değiştirerek Pareto cephesini yaklaşık olarak oluşturmuştur. Ağırlıklı toplam yöntemi, doğrusal olmayan ve dışbükey olmayan bir problemde yalnızca desteklenen Pareto çözümlerini yakalayabildiği için ortaya çıkan cephe bütün olası baskın olmayan çözümlerin eksiksiz kümesi değildir.

MS-PABnB algoritması neden fizik bilgili olarak tanımlanmaktadır?

Standart Dal ve Sınır algoritması bütün ikili değişkenleri genel matematiksel kararlar gibi ele alır. MS-PABnB ise gaz türbininin fiziksel bağımlılıklarını kullanarak bazı dalların denenmesini gereksiz kılar. Çalışmada dört temel eleme kuralı uygulanmıştır:

1. Bir eşanjör kurulmamışsa o eşanjörün ısı kapasitesi oranını belirleyen alt kararlar araştırılmaz.
2. İki eşanjörün ısı kapasitesi oranı seçicilerinin farklı değerleri fiziksel olarak uygulanamaz çözüm oluşturuyorsa ikinci değişken ilkine bağlanır.
3. Bir türbinin soğutması etkin değilse bütün soğutma havası çekiş noktaları sıfırlanır. Bir çekiş noktası seçildiğinde aynı türbin için diğer noktalar araştırılmaz.
4. Bir dalda uygulanabilir ve daha iyi bir mevcut çözüm erken bulunursa gereksiz alt dalların araştırılması durdurulur.

Modelde 16 ikili değişken bulunmaktadır. Bunların hiçbir ilişki kullanılmadan bütün kombinasyonları araştırılsaydı yalnızca ikili yapı için \(2^{16}=65.536\) olasılık ortaya çıkacaktı. Fiziksel kurallar, gerçek dışı veya birbirini tekrar eden kombinasyonları çözüm ağacının erken aşamalarında kaldırmaktadır.

Başlangıç noktaları nasıl oluşturulmuştur?

Doğrusal olmayan optimizasyonlar başlangıç değerlerine duyarlı olduğu için araştırmacılar üç aşamalı bir başlangıç noktası üretim sistemi geliştirmiştir:

1. Tam faktöriyel deney tasarımı: Toplam hava debisi, üçüncü akış oranı, yeniden ısıtılan akış oranı, toplam basınç oranı ve eşanjör alanları belirlenen sınırlar içinde örneklenmiştir.
2. Vekil modelle eleme: AdaBoost sınıflandırıcısı bir aday çevrimin hedef itkiyi sağlayıp sağlayamayacağını tahmin etmiştir. Uygulanabilir görülen adayların türbin giriş sıcaklığı LSBoost regresyon modeliyle kestirilmiştir.
3. Yakıt-hava oranı belirleme: Elemeden geçen adaylarda hedef itkiyi karşılayacak yakıt-hava oranı hesaplanmış; termodinamik kısıtları sağlayamayan noktalar çıkarılmıştır.

Yalnızca deney tasarımı kullanıldığında 360 adayın 21'i uygulanabilir olmuş ve fizibilite oranı %5,83'te kalmıştır. Vekil model elemesiyle 160 aday uygulanabilir bulunmuş ve oran %44,4'e çıkmıştır. Yakıt-hava oranı belirleme aşamasından sonra 319 adayın uygulanabilir hâle gelmesiyle oran %88,6'ya ulaşmıştır.

Model nasıl doğrulanmıştır?

Termodinamik model, üç referans motor sınıfında GASTURB yazılımıyla karşılaştırılmıştır. Kullanılan referanslar yaklaşık 10 kN sınıfında ALF502L, 30 kN sınıfında PW2040 ve 50 kN sınıfında CFM56-2C1 motorlarıdır. Net itki ve özgül yakıt tüketimindeki farklar bütün sınıflarda %0,1'in altında kalmıştır.

Bununla birlikte iç motor değişkenlerinde farklar daha büyüktür. Örneğin 50 kN sınıfında düşük basınç türbini çıkış basıncı GASTURB modelinde 70,69 kPa, geliştirilen modelde 65,75 kPa; çekirdek nozul çıkış hızı sırasıyla 478,04 ve 469,48 m/s hesaplanmıştır. Bu nedenle doğrulama, toplam itki ve yakıt tüketimi bakımından güçlü uyum gösterse de bütün iç istasyon değerlerinin aynı kesinlikte eşleştiği şeklinde yorumlanmamalıdır.

Ağırlık modeli, üç motor büyüklüğünde ticari motor veri tabanlarıyla karşılaştırılmıştır. Ortalama mutlak yüzde hata yaklaşık 10 kN sınıfında %10,8, 30 kN sınıfında %5,6 ve 50 kN sınıfında %6,3 olarak bildirilmiştir. Aynı sınıflardaki hata standart sapmaları sırasıyla %6,0, %4,9 ve %4,6'dır. Bazı tekil motorlarda hata %15-16 düzeyine çıkmaktadır.

Algoritma karşılaştırması ne göstermiştir?

Model, GAMS 50.4.1 ve LINDO API 14.0 ortamında Genelleştirilmiş İndirgenmiş Gradyan yöntemi kullanılarak çözülmüştür. Hesaplamalarda yaklaşık 50 çoklu başlangıç noktası, \(10^{-12}\) fizibilite toleransı ve her çözüm için 3.600 saniyelik üst zaman sınırı kullanılmıştır.

Model istatistikleri tablosu 483 kısıt bildirmektedir: 451 eşitlik, 17 küçük-eşit ve 15 büyük-eşit kısıtı. Bunların 351'i doğrusal olmayan, 132'si doğrusaldır. Toplam 471 değişken içinde 16 ikili değişken bulunmaktadır. Ancak sonuç bölümünde “471 sürekli değişken ve 16 ayrık karar” ifadesi kullanılmıştır. Bu ifade, 471 toplam değişken gösteren tabloyla tutarlı değildir ve yüklenen sürümde açıklanmayan bir yazım veya sınıflandırma farklılığıdır.

Başlangıç noktası duyarlılığı grafikleri, üç itki sınıfında da MS-PABnB'nin referans en iyi çözüme geleneksel çok başlangıçlı Dal ve Sınır yönteminden daha az başlangıç noktasıyla yaklaştığını göstermektedir. 50 kN sınıfında düğüm başına hesaplama maliyeti bazı bölgelerde yükselse de, gerekli başlangıç noktası sayısının azalması toplam süre avantajını korumuştur.

Düğüm sayısı grafikleri, MS-PABnB'nin Pareto cephesi boyunca daha az düğüm araştırdığını ve çeyrekler arası dağılımının daha dar olduğunu göstermektedir. Özellikle 30 kN sınıfında verimlilik ağırlığının yaklaşık 0,5-0,6 olduğu bölgede geleneksel yöntemin düğüm sayısı belirgin biçimde yükselirken fizik bilgili yaklaşımda aynı sıçrama görülmemiştir. Normalize edilmiş amaç fonksiyonu grafiklerinde de MS-PABnB üç sınıfta genel olarak daha düşük ortalama değerler üretmiştir.

Ana motor tasarımı bulguları nelerdir?

Seçilen temsilî tasarımlarda ara soğutucu ve ısı geri kazanım eşanjörü hiçbir itki sınıfında kullanılmamıştır. Eşanjörlerin oluşturduğu ağırlık ve basınç kaybı, seçilen verimlilik-ağırlık dengesinde sağlayabilecekleri termodinamik kazancı aşmıştır. Araştırmacılar yalnızca verimliliğe çok yüksek ağırlık verilen başka Pareto noktalarında ısı geri kazanımının avantajlı olabildiğini belirtmektedir.

30 ve 50 kN sınıflarındaki temel kazanç; çekirdek motorun küçültülmesi, akışın bir bölümünün döner detonasyon bölümünde yeniden ısıtılması ve bypass ile üçüncü akış oranlarının yükseltilmesinden doğmuştur. Yeniden ısıtma çekirdek akışın güç yoğunluğunu artırmış; daha küçük çekirdekle hedef itkinin korunmasına ve fan çevresinden daha fazla hava geçirilmesine imkân vermiştir.

30 kN tasarımında ONAED, OCA'ya göre küresel verimliliği %15,40 artırmış ve ağırlığı %22,38 azaltmıştır. 50 kN tasarımında verimlilik artışı %15,73, ağırlık azalması %18,23'tür. 10 kN sınıfında iki mimari neredeyse aynı sonucu vermiştir: verimlilik farkı yalnızca %0,38, ağırlık farkı %0,19'dur.

10 kN sınıfındaki sınırlı kazancın nedeni, hem yeni hem geleneksel mimarinin çekirdek debisi için belirlenen 17,57 kg/s alt sınıra ulaşmasıdır. Çekirdeğin daha fazla küçültülememesi, yeni bileşenlerin bypass oranını artırmak için sağlayabileceği ek tasarım serbestliğini sınırlandırmıştır.

Bütün seçilen tasarımlarda toplam basınç oranı 33,1 olan üst sınıra ulaşmıştır. Yanma odası çıkış sıcaklığı 10 kN tasarımında yaklaşık 1.659-1.668 K, 30 ve 50 kN tasarımlarında 1.860 K üst sınırındadır. Bu sonuçlar, optimizasyonun verilen kısıtlar içinde yüksek basınç oranı ve yüksek tepe sıcaklığını tercih ettiğini göstermektedir; ancak malzeme ömrü, ayrıntılı soğutma tasarımı ve azot oksit oluşumu bütün fiziksel ayrıntılarıyla modellenmemiştir.

Pareto grafikleri nasıl yorumlanmalıdır?

Pareto grafiklerinde ONAED çözüm uzayı, OCA çözüm uzayını tamamen içermektedir. OCA, yeni mimarideki isteğe bağlı bileşenlerin tamamının kapatıldığı özel bir alt problemdir. Bu nedenle ONAED'in OCA'dan daha kötü olmaması, yalnızca yeni teknolojilerin fiziksel başarısından değil, optimizasyon probleminin matematiksel olarak daha geniş bir tasarım uzayına sahip olmasından da kaynaklanmaktadır.

30 ve 50 kN sınıflarında iki Pareto cephesi arasındaki ayrım belirgindir. 10 kN sınıfında ise cepheler büyük ölçüde çakışmaktadır. Hiperhacim göstergesi yeni mimarinin referans motora göre baskın olduğu amaç uzayını ölçmektedir. Mutlak hiperhacim kazancı 30 kN sınıfında, ONAED'in OCA'ya göre oransal üstünlüğü ise 50 kN sınıfında daha belirgindir.

Seyir görevi ve emisyon sonuçları nelerdir?

Seçilen tasarımlar üç temsilî seyir görevinde değerlendirilmiştir. Yaklaşık 10 kN sınıfı için 2,5 saatlik Canouan-Fort Lauderdale, 30 kN sınıfı için 4 saatlik Lihue-Los Angeles ve 50 kN sınıfı için 8 saatlik São Paulo-Washington seyir bölümleri kullanılmıştır. Görevler gerçek uçuş rotalarından türetilmiş olsa da yeni motorlar gerçek uçaklara yerleştirilerek uçurulmamıştır.

30 kN sınıfında ONAED'in yakıt tüketimi OCA'ya göre %19,62, referans PW2040 modeline göre %27,64 daha düşüktür. 50 kN sınıfında azalma OCA'ya göre %21,23, referans motora göre %27,58'dir. Karbondioksit sonuçları kullanılan sabit yakıt-CO₂ dönüşüm katsayısı nedeniyle yakıt tüketimindeki değişimi aynı yüzdeyle izlemektedir.

Azot oksit sonucu her sınıfta aynı yönde değildir. 30 kN ONAED tasarımı referans motordan %8,90, 50 kN tasarımı %22,68 daha düşük NOx üretmiştir. Buna karşılık 10 kN ONAED tasarımında NOx, referans ALF502L modelinden %86,86 daha yüksek hesaplanmıştır. Araştırmacılar bu artışı yüksek yanma sıcaklığı ve bu sınıfta yakıt tasarrufunun sınırlı kalmasıyla ilişkilendirmektedir.

CO₂, yakılan her kilogram Jet-A1 için 3,159 kg CO₂ katsayısıyla hesaplanmıştır. NOx değerleri Boeing Fuel Flow Method 2 temelli basınç ve sıcaklık ölçeklemesine dayanmaktadır. CO ve yanmamış hidrokarbonlar yüksek güçlü seyirde sabit artık emisyon değerleriyle modellenmiştir. Dolayısıyla emisyon tablosu doğrudan egzoz gazı ölçümü değil, görev boyunca hesaplanan yakıt akışı ve ampirik emisyon ilişkilerinin sonucudur.

Çalışmanın güçlü yönleri nelerdir?

• Motor mimarisi, termodinamik çalışma noktası ve ağırlık aynı optimizasyon probleminde birleştirilmiştir.
• Yeni mimari ile geleneksel mimari aynı kısıtlar ve aynı model altında karşılaştırılmıştır.
• Termodinamik model üç motor sınıfında GASTURB sonuçlarıyla karşılaştırılmıştır.
• Ağırlık modeli çok sayıda gerçek motor verisiyle ayrı olarak değerlendirilmiştir.
• Başlangıç noktası üretimi, vekil makine öğrenmesi modelleri ve fizik tabanlı dal eleme kuralları nicel olarak karşılaştırılmıştır.
• Verimlilik-ağırlık dengesi tek bir tasarım yerine bütün Pareto cephesi boyunca incelenmiştir.
• Yakıt ve emisyon etkileri üç farklı itki ve seyir süresi sınıfında değerlendirilmiştir.

Çalışmanın sınırlılıkları nelerdir?

• Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.
• Gerçek motor prototipi, test hücresi deneyi veya uçuş doğrulaması bulunmamaktadır.
• Optimizasyon yalnızca seyir tasarım noktasına odaklanmaktadır; kalkış, tırmanış, alçalma ve değişken çalışma koşulları birlikte optimize edilmemiştir.
• Kompresör ve türbin performans haritaları, geçici rejimler ve motor kontrol dinamikleri açık biçimde modellenmemiştir.
• Motor kaportası sürüklemesi, hava aracı entegrasyonu, montaj etkileri ve fan çapının bütün aerodinamik sonuçları modele dâhil değildir.
• Akustik kısıtlar ve döner detonasyon odasının ayrıntılı basınç salınımları değerlendirilmemiştir.
• Döner detonasyon odası, karıştırma bölümü ve eşanjör ağırlıkları yaklaşık korelasyonlara dayanmaktadır.
• Isı eşanjörlerinin ayrıntılı kanatçık, kanal ve akış geometrileri çözülmemiştir.
• Ağırlıklı toplam yöntemi dışbükey olmayan Pareto cephesinin bütün baskın olmayan noktalarını bulmayı garanti etmemektedir.
• Algoritmanın karşılaştırıldığı referans optimum, 150 başlangıç noktasıyla bulunan en iyi bilinen çözümdür; yüklenen metinde bağımsız bir küresel optimalite sertifikası sunulmamaktadır.
• ONAED çözüm uzayı OCA'yı içerdiğinden ONAED'in Pareto üstünlüğünün bir bölümü problem yapısından kaynaklanmaktadır.
• Yüklenen dosyada atıf yapılan 55-116 numaralı denklemler ile 12-17 numaralı tabloları içeren ekler bulunmamaktadır.

Çalışma neyi desteklemektedir?

Çalışma, termodinamik çevrim ve ağırlık tahmini birlikte optimize edildiğinde, özellikle 30 ve 50 kN seyir itkisi sınıflarında yeniden ısıtmalı ve üçüncü akışlı motor mimarilerinin model içinde geleneksel iki akışlı tasarımlardan daha iyi verimlilik-ağırlık dengeleri üretebildiğini desteklemektedir. Ayrıca fiziksel bağımlılıkları kullanan dal eleme kurallarının aynı çözüm ortamında geleneksel Dal ve Sınır yaklaşımından daha az düğümle ve daha kararlı sonuçlarla çalışabildiğini göstermektedir.

Çalışma neyi kanıtlamamaktadır?

Çalışma, önerilen motorların üretilebilir, sertifikalandırılabilir veya mevcut hava araçlarına doğrudan takılabilir olduğunu kanıtlamamaktadır. Hesaplanan ağırlık azalmaları gerçek üretim motoru tartımları değildir. Yakıt ve emisyon değerleri uçuşta ölçülmemiştir. Döner detonasyon sisteminin mekanik ömrü, gürültüsü, titreşimi, soğutma ihtiyacı ve güvenilirliği doğrulanmamıştır. Seyir noktasındaki yüksek bypass ve küçük çekirdek tercihinin kalkış ve tırmanış performansını karşılayacağı da gösterilmemiştir.

Türkiye açısından yöntemsel değeri nedir?

Çalışmanın Türkiye açısından esas değeri, hazır bir motor tasarımı sunmasından çok, farklı motor bileşenlerinin aynı sayısal ortamda seçilmesini ve boyutlandırılmasını sağlayan yöntemsel çerçevesidir. Gaz türbini, turbofan, yardımcı güç ünitesi veya başka enerji dönüşüm sistemlerinde çalışan araştırma ekipleri; bileşen kurulumu gibi ayrık kararlarla basınç oranı, debi, sıcaklık ve geometri gibi sürekli değişkenleri ortak bir MINLP modeli içinde değerlendirebilir.

Bununla birlikte bu yaklaşımın gerçek bir mühendislik programına aktarılabilmesi için yüksek doğruluklu hesaplamalı akışkanlar dinamiği, malzeme ve ömür analizleri, bileşen haritaları, test hücresi verileri, üretim kısıtları ve bütün görev profiliyle doğrulanması gerekir. Çalışma bu aşamaların yerine geçmemekte; ön tasarım uzayını daraltabilecek bir hesaplama yöntemi önermektedir.

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

Araştırma tasarımının teknik özeti

Model sabitleri ve çalışma koşulları

• İdeal gaz sabiti: 287,05 J/(kg K).
• Yakıt ısıl değeri: \(44\times10^6\) J/kg.
• Karıştırma odası için IN-738 yoğunluğu: 8.110 kg/m³.
• Isı eşanjörü alan yoğunluğu: 3,6 kg/m².
• Toplam basınç oranı aralığı: 15-33,1.
• Yeniden ısıtılan akış oranı aralığı: 0-0,5.
• Ara soğutucu ve geri kazanım eşanjörü yüzey alanı: 0-750 m².
• Deniz seviyesine düzeltilmiş toplam debi örnekleme aralığı: 17,57-802,4 kg/s.

Termodinamik doğrulama sonuçları

En önemli nicel bulgular

• Başlangıç noktası fizibilitesi %5,83'ten %88,6'ya yükselmiştir.
• MS-PABnB üç itki sınıfında geleneksel yönteme göre daha az arama düğümü kullanmıştır.
• 30 kN ONAED tasarımında küresel verimlilik 0,442, OCA tasarımında 0,383'tür.
• 30 kN ONAED ağırlığı 1.634 kg, OCA ağırlığı 2.104 kg'dır.
• 50 kN ONAED verimliliği 0,412, OCA verimliliği 0,356'dır.
• 50 kN ONAED ağırlığı 2.425 kg, OCA ağırlığı 2.966 kg'dır.
• 30 kN seyir görevinde ONAED yakıt tüketimi referans motordan %27,64 daha düşüktür.
• 50 kN seyir görevinde ONAED yakıt tüketimi referans motordan %27,58 daha düşüktür.
• 10 kN sınıfında ONAED'in NOx sonucu referans motordan %86,86 daha yüksektir.

Yöntemsel yorum sınırı

Bu sonuçlar, tanımlanan değişken sınırları, sabit bileşen verimleri, yaklaşık ağırlık bağıntıları ve seyir tasarım noktası için geçerlidir. Modelin daha düşük ağırlık üretmesi, fiziksel motorun aynı ağırlıkla üretilebileceğini tek başına göstermez. Benzer şekilde seyir yakıtındaki azalma, kalkıştan inişe kadar bütün görev yakıtının aynı oranda azalacağını kanıtlamaz.

Kaynak ve Yöntem Notu

• Çalışmanın tam özgün adı: A Physics-Acquainted Branch-and-Bound Algorithm for Multi-Objective Mixed-Integer Nonlinear Optimization of Gas Turbine Engines
• Yazarlar ve sıraları: Carlos Ávila Catalán; Jorge Saavedra; Luis Cadarso
• Eş birinci yazarlık: Eş katkı veya eş birinci yazarlık bilgisi belirtilmemiştir.
• Sorumlu yazar: Carlos Ávila Catalán
• Sorumlu yazar e-postası: carlos.avila@urjc.es
• Kurum: Universidad Rey Juan Carlos, Camino del Molino 5, Fuenlabrada, 28943 Madrid, İspanya
• DOI: 10.2139/ssrn.6945142
• Hakemli dergi: Hakemli bir dergi adı doğrulanmamıştır. Dosyada çalışma “Preprint submitted to Elsevier” olarak tanımlanmaktadır.
• Yayın platformu: SSRN / SSRN Electronic Journal
• Özgün yayınevi veya platform işletmecisi: Elsevier bünyesindeki SSRN
• Yayın yılı: 2026
• Kaynak türü: Sayısal modelleme ve çok amaçlı optimizasyon araştırması niteliğinde preprint
• Hakemlik durumu: Bu çalışma bir preprinttir ve hakem değerlendirmesinden geçmemiştir.
• Resmî preprint bağlantısı:SSRN çalışma sayfası
• DOI bağlantısı:10.2139/ssrn.6945142

Bu Türkçe açıklama, yüklenen çalışmanın metni, formülleri, tabloları, motor mimarisi şeması, algoritma akış şeması, Pareto grafikleri ve emisyon tabloları incelenerek hazırlanmıştır. Çalışmada bulunmayan yeni bir bilimsel sonuç, motor performansı veya uygulama başarısı eklenmemiştir. Dış kaynak kontrolü yalnızca yazar, kurum, DOI, platform ve yayın durumu gibi bibliyografik kimlik bilgilerinin doğrulanması amacıyla kullanılmıştır.

Yüklenen dosya 36 sayfada sona ermekte, ancak iç sayfa numaralandırması toplam 55 sayfaya işaret etmektedir. Ana metinde atıf yapılan ek denklemler, model kümeleri, kalibrasyon parametreleri ve 12-17 numaralı tablolar dosyada bulunmamaktadır. Bu eksiklik, modelin bağımsız olarak bütünüyle yeniden kurulmasını ve bütün katsayıların denetlenmesini engellemektedir.

Çalışmanın temel sınırlılıkları; seyir tasarım noktasına odaklanması, deneysel motor doğrulaması içermemesi, geçici ve tasarım dışı çalışma koşullarını modellememesi, yeni bileşenlerin ağırlığını yaklaşık bağıntılarla hesaplaması, ayrıntılı motor-hava aracı entegrasyonunu dışarıda bırakması ve emisyonları doğrudan ölçmek yerine ampirik yöntemlerle tahmin etmesidir. Sonuçlar bu sınırlar içinde bir ön tasarım ve hesaplamalı yöntem çalışması olarak değerlendirilmelidir.