Bu çalışma, kompresör içindeki akışkan ile gövde, kanat ve disk gibi katı parçalar arasındaki ısı transferinin ölçülen ve hesaplanan kompresör performansını nasıl etkilediğini NASA Rotor 37 transonik eksenel kompresörü üzerinde incelemiştir. Araştırmada adyabatik duvar, 288,15 K sıcaklıktaki izotermal gövde, yalnızca gövdeyi içeren eşlenik ısı transferi ve gövde-kanat-disk katı bölgelerini içeren eşlenik ısı transferi olmak üzere dört sayısal yapılandırma karşılaştırılmıştır. Isıl sınır koşulları toplam basınç oranını yalnızca sınırlı ölçüde değiştirirken, özellikle göbek ve dış gövde yakınındaki toplam sıcaklık dağılımını ve geleneksel yöntemle hesaplanan verim eğrilerini belirgin biçimde değiştirmiştir. Bununla birlikte çalışma, sıcaklık farkından hesaplanan verim değişiminin tamamının aerodinamik kayıp değişimi olarak yorumlanamayacağını göstermektedir.

En önemli termodinamik bulgu, akıştan katı bölgelere geçen ısı eşdeğer adyabatik verim hesabına dâhil edildiğinde dört yapılandırmanın verim eğrilerinin büyük ölçüde birbirine yaklaşmasıdır. Bu sonuç, geleneksel verim eğrileri arasındaki farkların önemli bir bölümünün akışın tersinmezliğindeki gerçek bir değişimden değil, enerji bilançosunda ısı akısının hesaba katılmamasından kaynaklandığını düşündürmektedir. Fiziksel duruma en yakın seçenek olarak değerlendirilen gövde-kanat-disk eşlenik ısı transferi yapılandırması, geleneksel verim hesabında yalnızca yaklaşık %0,3 düzeyinde bir iyileşme üretmiştir; çalışmada bu değerin mutlak yüzde puanı mı yoksa göreli değişim mi olduğu ayrıca açıklanmamıştır.

Akış alanı sonuçları, ısı transferinin tamamen etkisiz olmadığını da göstermektedir. Gövde-kanat-disk eşlenik çözümünde göbek yakınındaki kanat işi ve toplam basınç oranı küçük fakat tutarlı biçimde yükselmiş, kanadın emme yüzeyindeki arka kenar köşe ayrılmasının başlangıcı gecikmiştir. Buna karşılık izotermal gövde koşulu deneysel verim eğrisine daha yakın görünmesine rağmen gerçek eşlenik ısı transferi çözümünden çok daha fazla soğutma oluşturmuş ve uç duvar sıcaklık profilini fiziksel süreçle uyuşmayacak ölçüde değiştirmiştir. Dolayısıyla yalnızca deneysel toplam verim eğrisine yakınlık, seçilen duvar sınır koşulunun fiziksel doğruluğunu kanıtlamamaktadır.

Sonuçlar yalnızca tek kademeden yalıtılmış NASA Rotor 37 geometrisinin, ağırlıklı olarak tasarım hızının %80’indeki sayısal analizine dayanmaktadır. Çalışma bir preprinttir; doğrudan katı sıcaklığı veya ısı akısı ölçümleriyle doğrulanmış bir eşlenik ısı transferi deneyi, çok kademeli kompresör analizi ya da gerçek motor koşullarında saha doğrulaması sunmamaktadır. Bu nedenle bulgular, belirli model varsayımları ve sınır koşulları içinde değerlendirilmelidir.

Türkiye açısından olası önemi

Çalışmanın Türkiye’ye ilişkin doğrudan bir veri seti veya uygulama sonucu bulunmamaktadır. Bununla birlikte yöntemsel sonuç, Türkiye’de yürütülebilecek gaz türbini, uçak motoru ve yüksek hızlı kompresör araştırmaları açısından önemlidir: Bir hesaplamalı akışkanlar dinamiği modelinin deneysel verim eğrisine yaklaşması, kullanılan ısıl sınır koşulunun mutlaka fiziksel olarak doğru olduğu anlamına gelmemektedir. Bulguların Türkiye’deki tasarım çalışmalarına aktarılabilmesi için gövde, kanat ve disk sıcaklıklarının ölçüldüğü deneylerle eşlenik modellerin doğrulanması; farklı hızlar, kademeler ve çalışma noktalarının incelenmesi; ayrıca ilk kademeden çıkan radyal sıcaklık profilinin sonraki kademelerin eşleşmesine etkisinin değerlendirilmesi gerekir. Çalışma, herhangi bir yerli motorun performansını, güvenilirliğini veya endüstriyel başarısını kanıtlamamaktadır.

Araştırmanın temel sorusu nedir?

Araştırmanın temel sorusu, kompresörde akışkan ile katı yüzeyler arasında gerçekleşen ısı transferinin iki farklı sonucu birbirinden nasıl ayırdığıdır: Birincisi kompresörün gerçek aerodinamik kayıpları ve kanat işi, ikincisi ise giriş-çıkış sıcaklıkları üzerinden hesaplanan görünür termodinamik verimdir. Kompresör hesaplamalarında kanat, göbek ve dış gövde yüzeylerinin çoğu zaman adyabatik kabul edilmesine karşın gerçek deney düzenekleri ve motorlar çevreleriyle ısı alışverişinde bulunur. Bu nedenle çıkış toplam sıcaklığındaki bir değişiklik yalnızca akışa aktarılan mekanik işi veya aerodinamik kaybı değil, katı yüzeylerden geçen ısıyı da içerebilir.

Araştırmacılar özellikle şu ayrımı incelemiştir: İzotermal bir dış gövde koşulunun deneysel verim eğrisine daha iyi uyması, bu koşulun gerçek fiziksel ısı transferini doğru temsil ettiği anlamına mı gelmektedir; yoksa uyum, yalnızca çıkış sıcaklığının yapay biçimde düşürülmesinden mi kaynaklanmaktadır?

Çalışma neden önemlidir?

Bir kompresörün toplam basınç oranı, çıkış sıcaklığı, verimi ve kararlılık sınırı motor veya gaz türbini tasarımının temel performans göstergeleridir. Duvarlardan dışarı aktarılan ısı, çıkış sıcaklığını düşürdüğünde geleneksel adyabatik verim denklemi daha yüksek bir verim gösterebilir. Ancak bu yükseliş, kanat geçidindeki şok dalgalarının, sınır tabakaların, uç kaçak akışının veya ayrılma bölgelerinin gerçekten iyileştiği anlamına gelmeyebilir.

Bu ayrım çok kademeli kompresörlerde daha da önemlidir. Bir kademenin çıkışındaki radyal sıcaklık profili, sonraki kademenin giriş koşulunu oluşturur. Toplam basınç oranı çok az değişse bile göbek veya dış gövde yakınındaki sıcaklık dağılımının değişmesi, kademeler arası eşleşmeyi etkileyebilir. Çalışma bu çok kademeli etkiyi doğrudan hesaplamamış, gelecekte araştırılması gereken bir konu olarak tanımlamıştır.

İncelenen kompresör ve çalışma koşulları

Araştırma nesnesi, NASA Glenn Research Center tarafından gelişmiş ve sekiz kademeli bir çekirdek kompresörün giriş kademesi olarak tasarlanmış NASA Rotor 37’dir. Rotor 37, transonik rejimde çalışan, düşük en-boy oranlı 36 adet çoklu dairesel yay biçimli kanada sahip yalıtılmış bir eksenel kompresör rotorudur.

Şekil 1, rotor kanadı, göbek, dış gövde, giriş ve çıkış bölgelerini hem üç boyutlu hem meridyen düzleminde göstermektedir. Meridyen görünümünde hesaplama bölgesinin eksenel uzunluğu yaklaşık 14,86 cm olarak verilmiştir.

Karşılaştırılan dört ısıl yapılandırma

Yöntem ve Şekil 2’deki tanımlara göre doğru yapılandırma eşleşmesi bu şekildedir. Buna karşın sonuç bölümündeki bazı cümlelerde yapılandırma 2 ile 3’ün açıklamaları karışmıştır. Ayrıca sonuç kısmında yapılandırma 4 zaman zaman “gövde-disk CHT” olarak kısaltılmış olsa da yöntem bölümünde kanat katı bölgesinin de çözüme dâhil edildiği açıkça belirtilmiştir.

Sayısal çözüm nasıl gerçekleştirilmiştir?

Akışkan bölgesindeki üç boyutlu Reynolds ortalamalı Navier-Stokes denklemleri ve katı bölgelerdeki ısı transferi denklemleri CFX 2025 yazılımıyla çözülmüştür. Karmaşık sıkıştırılabilir akış için kayma gerilmesi taşınımı türbülans modeli (shear stress transport, SST) ve otomatik duvar işlemi kullanılmıştır. Taşınım ve türbülans terimleri yüksek çözünürlüklü şemayla ayrıklaştırılmıştır.

Rotorun çevresel olarak tekrarlanan geometrisi nedeniyle tek kanat geçidi ve dönel periyodiklik koşulu kullanılmıştır. Girişte toplam basınç 101,325 kPa, toplam sıcaklık 288,15 K ve eksenel akış tanımlanmış; çıkışta statik basınç uygulanmıştır. Dış gövde duvarı için karşıt dönme duvar koşulu, akışkan-katı arayüzlerinde ise sıcaklık ile normal ısı akısı sürekliliğini sağlayan frozen-rotor yaklaşımı kullanılmıştır. Bütün katı bölgelerin malzemesi çelik kabul edilmiştir.

Yakın tepe verimi ve kararsızlık sınırına yakın seçilmiş çalışma noktalarında zamana bağlı hesaplamalar yapılmıştır. Araştırmacılar bu hesaplamaların kararlı çözümlere göre biraz daha düşük toplam basınç oranı ve daha yüksek verim verdiğini, böylece deneysel sonuçlara yaklaştığını bildirmiştir. Ancak zaman adımı, toplam fiziksel çözüm süresi ve zamansal ağ bağımsızlığı gibi ayrıntılar metinde verilmemiştir.

Katı bölge enerji denklemi

Katı bölgedeki enerji taşınımı aşağıdaki denklemle ifade edilmiştir:

\[ \frac{\partial(\rho h)}{\partial t}+\nabla\cdot(\rho U_s h)=\nabla\cdot(\lambda\nabla T)+S_E \]

Burada ρ yoğunluğu, h özgül entalpiyi, t zamanı, U_s katı bölgenin hareket hızını, λ ısıl iletkenliği, T sıcaklığı ve S_E hacimsel enerji kaynağını göstermektedir. Sol taraftaki terimler katıda enerji birikimini ve katı hareketiyle enerji taşınımını; sağ taraftaki iletim terimi sıcaklık gradyanı nedeniyle katı içinde yayılan ısıyı temsil eder. Hacimsel bir ısı kaynağı bulunması hâlinde bunun etkisi S_E ile hesaba katılır.

Dış yüzey taşınım denklemi

Dış gövde veya disk yüzeyinden çevreye geçen ısı akısı şu ilişkiyle hesaplanmıştır:

\[ q=h_c(T_o-T_w) \]

q, birim yüzey alanından geçen ısı akısıdır ve W/m² birimiyle ifade edilir. h_c taşınımla ısı transfer katsayısıdır ve W/(m²·K) birimindedir. T_o dış ortam sıcaklığını, T_w ise duvar sıcaklığını gösterir. İşaret, dış ortam ile duvar arasındaki sıcaklık farkına bağlı olarak ısının hangi yönde aktığını belirtir.

Ana metinde doğal taşınım için dış gövde katsayısı 25 W/(m²·K) olarak belirtilirken Şekil 2(b)’de 15 W/(m²·K) yazılmıştır. Bu iki değer arasındaki uyuşmazlık metinde açıklanmamıştır. Dış ortam sıcaklığı her iki durumda da 288,15 K’dir. Gövde-kanat-disk çözümünü gösteren Şekil 2(c), dış disk yüzeyi için yaklaşık 1100 W/(m²·K) değerini göstermektedir.

Dönen disk için Nusselt bağıntısı

Dış disk yüzeyindeki taşınım katsayısının belirlenmesinde sabit ve dönen diskler için verilen ampirik bağıntı kullanılmıştır:

\[ Nu_r=0{,}0253(Re_r)^{0{,}8} \]

Nu_r radyal konuma bağlı Nusselt sayısını, Re_r ise yerel Reynolds sayısını ifade eden boyutsuz büyüklüklerdir. Nusselt sayısı, taşınımla ısı geçişinin yalnızca iletimle gerçekleşecek ısı aktarımına oranını temsil eder. Reynolds sayısının 0,8 kuvveti, dönme ve akış hızının artmasıyla taşınımın doğrusal olmayan biçimde güçlendiğini ifade eder.

Geleneksel adyabatik verim

Deneysel ve sayısal performansın karşılaştırılmasında kullanılan geleneksel adyabatik verim denklemi şöyledir:

\[ \eta_{ad}=\frac{\pi^{(\gamma-1)/\gamma}-1}{T_{t,2}/T_{t,1}-1} \]

η_ad adyabatik verimi, π toplam basınç oranını, γ özgül ısı oranını, T_t,1 giriş toplam sıcaklığını ve T_t,2 çıkış toplam sıcaklığını gösterir. Denklem, belirli bir basınç artışı için ideal izentropik sıcaklık artışını gerçek toplam sıcaklık artışıyla karşılaştırır.

Bu tanım duvarlardan ısı geçmediği durumda termodinamik olarak uygundur. Akışkan katı yüzeylere ısı verdiğinde çıkış sıcaklığı yalnızca kompresör işi nedeniyle belirlenmez. Bu nedenle ısı kaybı, paydadaki sıcaklık artışını küçülterek gerçek aerodinamik kayıplar değişmese bile görünürde daha yüksek verim üretebilir.

Eşdeğer adyabatik verim

Isı akısını enerji bilançosuna geri eklemek için eşdeğer çıkış toplam sıcaklığı ve eşdeğer adyabatik verim tanımlanmıştır:

\[ \eta'_{ad}=\frac{\pi^{(\gamma-1)/\gamma}-1}{T'_{t,2}/T_{t,1}-1} \]

\[ T'_{t,2}=T_{t,2}+\frac{Q}{\dot{m}C_p} \]

T′_t,2, ısı transferi etkisi düzeltilmiş eşdeğer çıkış toplam sıcaklığıdır. Q akış alanından katı bölgelere geçen toplam ısı transferini, ṁ kütlesel debiyi ve C_p sabit basınçtaki özgül ısı kapasitesini gösterir. Q/(ṁC_p) terimi, birim akışkan kütlesinden katıya geçen enerjinin sıcaklık karşılığıdır.

Gündelik bir benzetmeyle, geleneksel hesap kompresörden çıkan havanın sıcaklığını yalnızca makinenin yaptığı işin sonucu kabul eder. Eşdeğer hesap ise havanın gövdeye bıraktığı ısıyı yeniden enerji bilançosuna ekleyerek, kompresör işini duvar soğutmasının etkisinden ayırmaya çalışır.

Toplam entalpi artışı ve Euler kanat işi

Akışkanın toplam entalpi artışı şu denklemle hesaplanmıştır:

\[ \Delta h_t=C_p(T_{t,2}-T_{t,1}) \]

Δh_t toplam entalpi artışıdır ve J/kg veya kJ/kg birimindedir. Bu değer, giriş ile çıkış toplam sıcaklık farkının özgül ısı kapasitesiyle çarpılmasıyla elde edilir. Isı transferinin bulunduğu durumda bu artış hem rotorun yaptığı işi hem de akışkanın katı yüzeylerle enerji alışverişini içerir.

Rotor kanatlarının akışkana aktardığı özgül iş Euler türbomakine denklemiyle hesaplanmıştır:

\[ L_u=U_2V_{t,2}-U_1V_{t,1} \]

L_u birim kütle başına kanat işini, U yerel kanat çevresel hızını ve V_t akışın çevresel hız bileşenini ifade eder. 1 ve 2 indisleri giriş ve çıkış konumlarını belirtir. Adyabatik durumda Δh_t ile L_u arasındaki fark duvar sürtünmesiyle ilişkilendirilirken, ısı transferli yapılandırmalar arasındaki ek fark katı yüzeylerden geçen ısıyı yansıtır.

Zorlanmış taşınım bağıntısı

Turbofan motorunda dış gövdenin fan çıkışındaki hızlı akışa maruz kalabileceği senaryo için düz levha korelasyonu kullanılmıştır:

\[ Nu=(0{,}03Re^{0{,}8}-871)Pr^{1/3} \]

Nu Nusselt sayısını, Re Reynolds sayısını ve Pr Prandtl sayısını ifade eder. Prandtl sayısı momentum yayılımı ile ısıl yayılım arasındaki oranı gösterir. Araştırmada dış akış hızı 100 m/s kabul edilmiş ve bu bağıntıyla dış gövde için 250 W/(m²·K) taşınım katsayısı elde edilmiştir.

Ağ bağımsızlığı ve deneysel karşılaştırma

Akışkan ve katı bölgeler için ayrı ağ bağımsızlığı incelemeleri yapılmıştır. Akışkan bölgesinde beş ağ seviyesi değerlendirilmiş; sabit 105.000 Pa çıkış statik basıncında kütlesel debi, toplam basınç oranı ve izentropik verim izlenmiştir. Duvara bitişik ilk hücre yüksekliği bütün yapılandırmalarda 0,003 mm seçilmiş ve y⁺ değerinin 3’ün altında tutulduğu bildirilmiştir.

Akışkan ağı 1,05 milyondan 1,22 milyon elemana çıkarıldığında kütlesel debi değişimi 2×10⁻³ kg/s, izentropik verim değişimi ise 4×10⁻⁴ içinde kalmıştır. Gövde katı ağı 65 binden 79 bin elemana çıkarıldığında verim değişimi 1×10⁻⁵, bütünleşik ısı akısı değişimi 7×10⁻⁴ W içinde kalmıştır.

Hesaplama maliyeti ile doğruluk arasındaki denge için akışkan bölgesinde yaklaşık 1,05 milyon, gövde katı bölgesinde 65 bin ve disk katı bölgesinde 79 bin eleman seçilmiştir. Nihai ağda kanat yüksekliği boyunca 87 hücre ve uç açıklığı içinde 17 radyal hücre bulunmaktadır. Şekil 5, kanat ucu ve akışkan-katı arayüzlerinde ağın yoğunlaştırıldığını göstermektedir.

Şekil 4’te adyabatik çözüm deneysel NASA Rotor 37 verileriyle karşılaştırılmıştır. Basınç oranı eğrisinin genel biçimi yakalanmış, ancak sayısal değerler çalışma aralığının büyük bölümünde deneysel basınç oranından biraz yüksek kalmıştır. Verim eğrisinin eğilimi yakalanmasına rağmen sayısal verim deneysel değerlerin altında kalmıştır. Metin, tepe izentropik verimindeki sapmanın %3 ve diğer bir sapmanın %1,5 içinde olduğunu bildirmektedir; ancak ikinci oranın toplam basınç oranına mı yoksa başka bir ölçüte mi ait olduğu cümlede açıkça belirtilmemiştir.

Genel performans eğrileri ne göstermektedir?

Şekil 6’nın temel mesajı, ısı transferi koşullarının toplam basınç oranını çok az, geleneksel denklemle hesaplanan verimi ise görünür biçimde değiştirdiğidir. Dört basınç oranı eğrisi çalışma aralığı boyunca neredeyse üst üste gelmiştir. Son kararlı hesaplama noktaları da birbirine yakın kalmış; yalnızca gövde eşlenik ısı transferi yapılandırmasında küçük bir kararlılık payı avantajı bildirilmiştir.

İzotermal gövde yapılandırması, adyabatik çözüme göre daha yüksek geleneksel verim üretmiş ve deneysel verim eğrisine daha fazla yaklaşmıştır. Buna karşılık yalnız gövdeyi veya gövde-kanat-diski gerçek katı iletimiyle çözen eşlenik modeller çok daha sınırlı bir termal etki göstermiştir. Fiziksel açıdan en gerçekçi kabul edilen yapılandırma 4’te verim artışı yaklaşık %0,3 olarak bildirilmiştir.

Bu karşılaştırma, izotermal koşulun deneyle daha iyi eşleşmesinin tek başına fiziksel doğrulama sayılamayacağını göstermektedir. İzotermal duvar, gövde sıcaklığını bütün yüzey boyunca 288,15 K’ye zorlayarak gerçek katı iletimi ve dış taşınımla oluşandan çok daha yüksek bir soğutma sağlamıştır.

Eşdeğer verim neden farklı bir sonuç vermektedir?

Şekil 7’de ısı akısı eşdeğer çıkış sıcaklığına eklendiğinde dört eşdeğer adyabatik verim eğrisi büyük ölçüde üst üste gelmektedir. Bu, geleneksel verim eğrileri arasındaki belirgin farkın benzer büyüklükte bir aerodinamik tersinmezlik farkını temsil etmediğini göstermektedir. Başka bir ifadeyle izotermal gövdenin yüksek görünen verimi, büyük ölçüde akıştan dışarı çekilen ısının kompresör işi gibi değerlendirilmemesinden kaynaklanmaktadır.

Çalışma bütün aerodinamik kayıpların değişmediğini kesin olarak kanıtlamamaktadır. Göbek yakınında küçük fakat ölçülebilir akış alanı değişimleri bulunmuştur. Bununla birlikte toplam verim eğrisindeki büyük görünür fark ile aerodinamik değişimin büyüklüğü aynı değildir.

Radyal basınç ve sıcaklık profilleri

Şekil 8’de tepe verimi koşulundaki toplam basınç oranı profilleri, kanat yüksekliğinin büyük bölümünde deneysel sonuçlarla benzer eğilim göstermiş ve dört ısıl yapılandırma arasında çok az fark sergilemiştir. Göbek yakınındaki büyütülmüş görünümde yapılandırma 4’ün küçük bir toplam basınç oranı avantajı görülmektedir.

Toplam sıcaklık oranı profilleri ise ısıl sınır koşuluna daha duyarlıdır. İzotermal gövde yapılandırması dış gövdeye yakın üst açıklık bölgesinde en büyük sıcaklık düşüşünü oluşturmuştur. Gövde CHT ve gövde-kanat-disk CHT sonuçları deneysel sıcaklık eğilimine izotermal sonuçtan daha yakın kalmıştır. Göbeğe en yakın yaklaşık %10’luk bölgede yapılandırma 4, disk ve kanat üzerinden gerçekleşen ısı transferi nedeniyle adyabatik göbekli yapılandırmalardan farklı bir sıcaklık profili üretmiştir.

Şekil 9’daki kararsızlık sınırına yakın durumda benzer davranış sürmüş, ancak dış gövde yakınındaki ısı transferinden etkilenen radyal bölge genişlemiştir. İzotermal yapılandırma yine dış gövde yakınında en büyük toplam sıcaklık azalmasını oluşturmuştur. Göbek bölgesinde yapılandırma 4’ün toplam basınç oranı ve entalpi dağılımındaki küçük farkları korunmuştur.

Entalpi artışı, kanat işi ve ısı akısı

Şekil 10 ve 11, toplam entalpi artışı ile Euler kanat işini ayrı ayrı göstermektedir. Toplam entalpi artışı, özellikle göbek ve dış gövde yakınında ısıl sınır koşuluna göre değişmiştir. Buna karşılık kanat işi profilleri açıklığın büyük bölümünde birbirine çok yakın kalmıştır. Bu sonuç, sıcaklık ve entalpi profilindeki değişimin tamamının rotorun akışkana aktardığı mekanik işteki değişimden kaynaklanmadığını desteklemektedir.

Yapılandırma 4’te göbek yakınında kanat işinin küçük ölçüde yükseldiği görülmüştür. Bu artışa aynı bölgede hafif bir toplam basınç oranı yükselişi eşlik etmiştir. Yapılandırma 2 çok daha büyük ısı akısı oluşturmasına rağmen benzer bir kanat işi veya basınç oranı kazancı üretmemiştir. Dolayısıyla yalnızca daha fazla soğutma, daha iyi aerodinamik performans anlamına gelmemektedir.

Şok dalgaları ve akış ayrılması

Şekil 12’de tepe verimi koşulunda emme yüzeyindeki ayrılmanın arka kenar yakınında başladığı ve göbeğe doğru güçlenerek köşe ayrılması niteliği kazandığı görülmektedir. Kanat yüksekliğinin yaklaşık %20’sinin üzerinde, eksenel kordun yaklaşık 0,4’ü civarında düşük duvar kayma gerilmesi bölgesi bulunmasına rağmen akış tepe verimi koşulunda yeniden yüzeye bağlı kalmaktadır.

Şekil 13’te kararsızlık sınırına yaklaşıldığında ayrılma daha belirgin hâle gelmekte, yerel bir ayrılma kabarcığı oluşmakta ve akış daha ileride yeniden yüzeye bağlanmaktadır. Şekil 14’te %90 açıklık kesitindeki ön kenar şokunun emme yüzeyine çarparak sınır tabakayla etkileştiği ve ayrılmayı başlattığı görülmektedir. Şekil 15’te %10 açıklık kesitindeki şok daha zayıftır ve aynı tür bir şok kaynaklı ayrılma oluşturmamaktadır.

Genel renk haritalarında dört yapılandırma arasındaki fark küçük olduğu için Şekil 16’da belirli radyal kesitlerde duvar kayma gerilmesi eğrileri karşılaştırılmıştır. %50 açıklıkta ayrılma ve yeniden bağlanma konumları dört çözümde yaklaşık aynıdır. %10 açıklıkta yapılandırma 4, arka kenar köşe ayrılmasının başlangıcını diğer yapılandırmalara göre geciktirmiştir. %90 açıklıkta izotermal yapılandırma şok kaynaklı ayrılmayı en fazla geciktiren çözüm olarak görünmektedir; ancak bu sonuç, izotermal koşulun gerçek dışı derecede güçlü soğutma oluşturduğu sınırlılığıyla birlikte değerlendirilmelidir.

Şekil 17’de kararsızlık sınırına yakın %10 açıklık kesitindeki statik basınç dalgalanmalarının kök ortalama kare değerleri karşılaştırılmıştır. Yapılandırma 4 en düşük dalgalanma alanını üretmiş ve araştırmacılar bunu göbek yakınında daha kararlı bir akış alanına işaret eden bir bulgu olarak yorumlamıştır. Bu sonuç bütün kompresör için deneysel bir kararlılık garantisi değildir; yalnızca incelenen sayısal kesit ve çalışma noktası için geçerlidir.

Dış gövdedeki zorlanmış taşınımın etkisi

Doğal taşınım yerine 100 m/s dış akışa karşılık gelen 250 W/(m²·K) zorlanmış taşınım katsayısı uygulandığında dış gövde üzerinden geçen ısı akısı belirgin biçimde artmıştır. Şekil 18, zorlanmış taşınım yapılandırmasının bazı gövde bölgelerinde doğal taşınıma göre daha yüksek mutlak ısı akısı oluşturduğunu göstermektedir.

Buna karşın Şekil 19’daki dış gövde duvar kayma gerilmesi dağılımı görünür ölçüde değişmemiştir. Araştırmacılar genel kompresör performansında da belirgin bir değişiklik tespit etmemiştir. Bu sonuç, dış yüzey soğutmasının ısı akısını artırmasının akış geçidindeki aerodinamik yapıyı aynı ölçüde değiştirmediğini göstermektedir.

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

• Adyabatik, izotermal ve iki farklı eşlenik ısı transferi yaklaşımı aynı kompresör ve çalışma koşulları altında karşılaştırılmıştır.
• Geleneksel verim ile ısı akısı düzeltilmiş eşdeğer verim ayrılarak termodinamik değerlendirme hatası görünür hâle getirilmiştir.
• Genel performans eğrilerinin yanında radyal profiller, kanat işi, entalpi artışı, duvar kayma gerilmesi, Mach sayısı ve basınç dalgalanmaları incelenmiştir.
• Akışkan ve katı bölgeler için ağ bağımsızlığı çalışmaları yapılmıştır.
• Adyabatik temel çözüm NASA Rotor 37 deneysel performans verileriyle karşılaştırılmıştır.
• Kararsızlık sınırına yakın seçilmiş noktalar zamana bağlı çözümlerle ayrıca incelenmiştir.
• Dış gövde doğal taşınımına ek olarak zorlanmış taşınım senaryosu değerlendirilmiştir.

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

• Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.
• Analiz yalıtılmış tek bir rotorla sınırlıdır; statorlar ve çok kademeli kompresör eşleşmesi hesaplanmamıştır.
• Ayrıntılı inceleme ağırlıklı olarak tasarım hızının %80’ine odaklanmaktadır.
• Katı bölgelerin tamamı çelik kabul edilmiştir; gerçek parça malzemeleri, temas dirençleri ve kaplamalar modellenmemiştir.
• Dış gövde taşınım katsayısı ampirik kabullerle belirlenmiştir ve doğal taşınım değeri metin ile Şekil 2 arasında tutarsızdır.
• Gövde, kanat ve disk sıcaklıklarının veya ısı akılarının doğrudan deneysel ölçümüyle CHT doğrulaması sunulmamıştır.
• Tek geçitli dönel periyodiklik yaklaşımı çevresel olarak asimetrik akış yapılarını temsil edemez.
• SST-RANS yaklaşımı bütün geçici türbülans ve ayrılma ölçeklerini doğrudan çözmez.
• Akışkan-katı arayüzünde frozen-rotor yaklaşımı kullanılmıştır.
• Zamana bağlı hesapların zaman adımı, örnekleme süresi ve zamansal yakınsama ayrıntıları verilmemiştir.
• Sayısal kararsızlık sınırı, çıkış basıncının 100 Pa artırılmasının çözümü ıraksatmasıyla belirlenmiştir; bu tanım deneysel döner stall veya surge ölçümüyle aynı değildir.
• Yalnızca ağ bağımsızlığı ayrıntılı olarak sunulmuş; malzeme özellikleri, dış taşınım katsayısı ve türbülans modeli için kapsamlı belirsizlik analizi verilmemiştir.

Çalışmanın desteklediği sonuçlar

• Isıl sınır koşulları çıkış toplam sıcaklığını ve özellikle uç duvarlara yakın radyal sıcaklık profilini belirgin biçimde etkileyebilir.
• Isı akısı hesaba katılmadan hesaplanan geleneksel kompresör verimi, aerodinamik kayıp değişimini olduğundan büyük gösterebilir.
• Isı akısı eşdeğer adyabatik verime dâhil edildiğinde dört yapılandırma arasındaki verim farkı büyük ölçüde azalır.
• Toplam basınç oranı dört ısıl yapılandırma arasında yalnızca küçük değişiklikler göstermiştir.
• Gövde-kanat-disk CHT çözümünde göbek yakınındaki ayrılma başlangıcı gecikmiş, kanat işi ve basınç oranı küçük ölçüde artmıştır.
• İzotermal gövdenin deneysel toplam verime daha yakın sonuç vermesi, sınır koşulunun fiziksel olarak doğru olduğunu tek başına kanıtlamamaktadır.
• Dış gövdedeki zorlanmış taşınım ısı akısını artırmış, ancak incelenen koşulda genel performansı görünür biçimde değiştirmemiştir.

Çalışmanın kanıtlamadığı sonuçlar

• Çalışma, izotermal gövde koşulunun gerçek bir kompresörde doğru sınır koşulu olduğunu kanıtlamamaktadır.
• Yaklaşık %0,3 verim değişimi gerçek motor ölçeğinde yakıt tüketimi veya görev performansı iyileşmesi anlamına gelmemektedir.
• Sayısal ayrılma gecikmesi, deneysel stall marjı veya surge sınırı artışını doğrudan kanıtlamamaktadır.
• Tek rotor sonuçları çok kademeli yüksek basınç kompresörlerine doğrudan genellenemez.
• Çalışma belirli bir motor, uçak veya enerji tesisinde saha başarısı göstermemektedir.
• Zorlanmış dış soğutmanın her hızda, her kompresörde veya daha yüksek sıcaklık farklarında etkisiz olduğu sonucuna varılamaz.
• Sayısal model uzun dönem dayanım, termal gerilme, yorulma veya parça ömrü değerlendirmesi sunmamaktadır.
• Çalışma ticari bir kompresör tasarımının güvenlik veya performans garantisini sağlamamaktadır.

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

Teknik yöntem özeti

Ana bulguların karşılaştırmalı özeti

Genel bilimsel değerlendirme

Araştırmanın en güçlü sonucu, kompresör veriminin yorumlanmasında termal enerji alışverişi ile aerodinamik kaybın birbirinden ayrılması gerektiğini göstermesidir. İzotermal gövde koşulu deneysel verim değerine daha yakın sonuç verse de ayrıntılı sıcaklık profilleri ve eşlenik çözümler, bu uyumun fiziksel ısı transferini doğru temsil etmediğini göstermektedir.

Yerel akış alanı açısından ısı transferi tamamen pasif değildir. Kanat ve disk katı bölgelerinin çözüme dâhil edilmesi göbek yakınındaki sıcaklık, ayrılma başlangıcı, basınç dalgalanması ve kanat işi üzerinde küçük etkiler oluşturmuştur. Ancak bu yerel iyileşmeler, geleneksel verim eğrilerinde görülen bütün farkı açıklayacak büyüklükte değildir.

Çalışmanın çok kademeli kompresörlere ilişkin sonucu bir öngörüdür: Uç duvar sıcaklık profili belirgin biçimde değiştiği için sonraki kademelerin giriş sıcaklığı ve kademe eşleşmesi etkilenebilir. Bu mekanizma mevcut araştırmada doğrudan test edilmemiştir. Gelecekte çok kademeli geometri, daha büyük gaz-katı sıcaklık farkları ve ölçülmüş katı sıcaklıklarıyla doğrulanmış CHT modelleri gereklidir.

Kaynak ve Yöntem Notu

Çalışmanın özgün adı: The Effect of Heat Transfer on the Compressor Performance

Yazarlar ve sıraları: Xuegao Wang; Hongwei Yan; Yuze Liu; You Zuo; Xinqian Zhen; Baotong Wang.

Eş birinci yazar: Çalışmada eş birinci yazarlık veya eşit katkı bilgisi belirtilmemiştir.

Sorumlu yazar: Hongwei Yan. Metinde sorumlu yazar e-posta adresi lepidum@outlook.com olarak verilmiştir.

Yazar-kurum eşleşmeleri:

• Xuegao Wang: National Key Laboratory of Transient Physics, Nanjing University of Science and Technology, Nanjing, Çin.
• Hongwei Yan: School of Vehicle and Mobility, Tsinghua University, Pekin, Çin; SinoHawk Power Technology Co., Ltd, Neijiang, Çin.
• Yuze Liu: State Key Laboratory of Multimodal Artificial Intelligence Systems, Institute of Automation, Chinese Academy of Sciences, Pekin, Çin; School of Artificial Intelligence, University of Chinese Academy of Science, Pekin, Çin.
• You Zuo: Digital Twin Research Center, Institute of Engineering Thermophysics, Chinese Academic of Science, Pekin, Çin; National Key Laboratory of Science and Technology on Advanced Light-Duty Gas-Turbine, Chinese Academic of Science, Pekin, Çin.
• Xinqian Zhen: School of Vehicle and Mobility, Tsinghua University, Pekin, Çin.
• Baotong Wang: Institute for Aero Engine, Tsinghua University, Pekin, Çin.

Kurum adları, çalışmada verildiği biçime mümkün olduğunca sadık tutulmuştur. Metinde bazı kurumlar için “Chinese Academy of Sciences”, bazıları için ise “Chinese Academic of Science” yazımı kullanılmıştır; bu yazım farkı sessizce tek biçime dönüştürülmemiştir.

DOI: 10.2139/ssrn.7201740

Kaynak türü: Sayısal modelleme ve simülasyon çalışması niteliğinde preprint araştırma makalesi.

Yayın platformu: SSRN.

Dergi: Hakemli bir dergi adı bu sürümde yer almamaktadır ve yayımlanmış dergi sürümü doğrulanamamıştır.

Yayınevi: Ayrı bir özgün dergi yayınevi bulunmamaktadır. SSRN, Elsevier’e ait bir preprint ve erken araştırma platformudur; bu kayıt hakemli nihai dergi yayını değildir.

Yayın yılı: 2026.

Hakemlik durumu: Bu çalışma bir preprinttir ve hakem değerlendirmesinden geçmemiştir. Metnin her sayfasında “Preprint not peer reviewed” uyarısı bulunmaktadır.

Resmî kaynak bağlantısı: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=7201740

Veri kullanılabilirliği: Yazarlar, verilerin talep üzerine sağlanacağını bildirmiştir.

Çıkar çatışması: Yazarlar, çalışmayı etkileyebilecek bilinen bir mali çıkar veya kişisel ilişki bulunmadığını beyan etmiştir.

Yazar katkısı notu: Makalenin yazar satırında altı araştırmacı bulunmasına rağmen CRediT katkı beyanında yalnızca Xuegao Wang ve Hongwei Yan’ın görevleri açıklanmıştır. Xuegao Wang için veri düzenleme, araştırma ve ilk taslak yazımı; Hongwei Yan için biçimsel analiz ile gözden geçirme ve düzenleme görevleri belirtilmiştir. Yuze Liu, You Zuo, Xinqian Zhen ve Baotong Wang’ın katkı rolleri bu sürümde ayrıca açıklanmamıştır.

Bu Türkçe açıklama, yüklenen 15 sayfalık çalışma baştan sona incelenerek hazırlanmıştır. Bilimsel yöntem, sayısal değerler, denklemler, şekiller ve sonuçlar yalnızca çalışmada sunulan içerikten aktarılmıştır. Dış kaynaklar yalnızca DOI, SSRN kaydı, platform niteliği ve hakemlik durumunun bibliyografik doğrulanması amacıyla kullanılmış; makalenin bilimsel bulgularına dışarıdan yeni sonuç eklenmemiştir.

Çalışmanın başlıca sınırlılıkları; hakem değerlendirmesinden geçmemesi, tek ve yalıtılmış rotor geometrisine dayanması, ayrıntılı analizlerin tasarım hızının %80’inde yoğunlaşması, gerçek katı sıcaklıkları ve ısı akılarıyla doğrudan deneysel CHT doğrulaması sunmaması, çok kademeli eşleşmeyi hesaplamaması ve doğal taşınım katsayısı ile bazı yapılandırma adlandırmalarında metin içi tutarsızlıklar bulunmasıdır. Sonuçlar bu sınırlılıklar dikkate alınarak okunmalıdır.