Bu çalışma, yakın uzay koşullarında azalan hava yoğunluğu ve Reynolds sayısının yüksek basınç oranlı transonik bir santrifüj kompresörde performansı nasıl düşürdüğünü ve kayıpların çark ile kanatlı difüzör arasında nasıl yeniden dağıldığını incelemiştir. Araştırmacılar, NASA CC3 santrifüj kompresörünün çark ve kanatlı difüzörünü birlikte içeren üç boyutlu, sürekli rejimli bir sayısal model kurmuş; 0, 10, 20 ve 25 km irtifaya karşılık gelen giriş koşullarını aynı düzeltilmiş dönme hızında karşılaştırmıştır.

Referans Reynolds sayısı deniz seviyesindeki 5,70 × 10⁵ değerinden 25 km’de 2,28 × 10⁴ değerine düşmüştür. Aynı aralıkta hesaplanabilen en yüksek toplam basınç oranı yaklaşık 4,50’den 3,92’ye, en yüksek izentropik verim 0,87’den 0,76’ya ve en yüksek düzeltilmiş kütlesel debi yaklaşık 4,90 kg·s⁻¹ değerinden 4,45 kg·s⁻¹ değerine gerilemiştir. En yüksek izentropik verimdeki göreli düşüş %12,9’dur.

Kayıp dağılımı irtifayla birlikte değişmiştir. Deniz seviyesinde toplam boyutsuz entropi artışının %58’i difüzörden, %42’si çarktan kaynaklanmıştır. Reynolds sayısı azaldıkça çark kaybı daha hızlı artmış; 20 km’de iki bileşen toplam kayba eşit katkı verir hâle gelmiştir. Bununla birlikte 10–25 km aralığında difüzör verimi %5,72, çark verimi ise %4,05 azalmıştır. Bu sonuç, difüzörün çark çıkışındaki bozulmuş akış koşullarına daha duyarlı olduğunu göstermektedir.

Türkiye açısından değerlendirme: Bulgular, Türkiye’de yüksek irtifa uzun havada kalışlı hava araçları, yakın uzay platformları, yüksek irtifa pistonlu motor destek sistemleri, turboşarjlar ve kompakt hava sıkıştırma makineleri üzerinde çalışan araştırma ekipleri için yararlı bir mekanizma çerçevesi sunmaktadır. Türkiye’ye özgü uygulama için modelin yerli kompresör geometrileriyle yeniden kurulması, gerçek giriş bozulmaları ve atmosfer koşullarının eklenmesi, zamana bağlı tam halka hesaplamaları yapılması ve yüksek irtifa test düzeneği veya düşük basınç odasında deneysel doğrulama gerçekleştirilmesi gerekir. Araştırmadan Türkiye’de kullanılan herhangi bir motorun basınç oranı, irtifa sınırı, kararlılık payı, yakıt tüketimi veya görev performansı doğrudan çıkarılamaz.

Çalışmanın ortaya koyduğu temel kayıp zinciri şöyledir: Reynolds sayısındaki azalma sınır tabakasının türbülansa geçişini geciktirmekte, özellikle muhafazaya yakın emme yüzeyinde düşük momentumlu akışkanın birikmesini kolaylaştırmakta, yerel ayrılma ve girdap çekirdeği yayılımını artırmakta, çark çıkışındaki akış açısını bozmakta ve bu bozulmuş akış difüzörün giriş kayıplarını büyütmektedir.

Araştırmanın temel sorusu nedir?

Yakın uzay ve yüksek irtifa koşullarında atmosfer basıncı ile yoğunluk büyük ölçüde azalır. Kompresörün geometrisi ve karakteristik hızı değişmese bile yoğunluktaki düşüş, akış içindeki atalet kuvvetlerinin viskoz kuvvetlere oranını ifade eden Reynolds sayısını azaltır. Bu değişim, kanat yüzeyindeki sınır tabakasını, uç açıklığı akışını, çark çıkışındaki hız dağılımını ve difüzör basınç geri kazanımını etkileyebilir.

Araştırma yalnızca “kompresör verimi ne kadar azalır?” sorusunu sormamaktadır. Asıl amaç, performans kaybının çark ile difüzör arasında nasıl dağıldığını ve çark içinde başlayan akış bozulmasının difüzör girişine hangi mekanizmayla taşındığını belirlemektir.

Çalışmanın ana mekanizma soruları şunlardır:

• Düşük Reynolds sayısı, kanat yüzeyindeki laminer-türbülanslı geçişi nasıl değiştirmektedir?
• Düşük momentumlu bölgeler hangi açıklık konumlarında yoğunlaşmaktadır?
• Muhafazaya yakın emme yüzeyinde ayrılma ve yerel geri dolaşım gelişmekte midir?
• Uç kaçağı ve girdap çekirdeklerinin uzamsal kapsamı nasıl değişmektedir?
• Çark çıkışındaki bozulma difüzör giriş açısını ve difüzör kaybını nasıl etkilemektedir?

NASA CC3 kompresörü neden seçilmiştir?

NASA CC3, yüksek basınç oranlı transonik santrifüj kompresörlerin sayısal incelenmesinde kullanılan, geometrisi ve deneysel performans verileri kamuya açık bir referans sistemidir. Çalışmada kanat profilleri, meridyen kanal koordinatları ve temel tasarım özellikleri NASA tasarım raporlarından yeniden oluşturulmuştur.

Şekil 1, kompresörün üstten, önden ve eğik görünümlerini sunmaktadır. Geometri; ana kanatları, ara kanatları, çark çıkışını ve çevresel olarak yerleştirilen kanatlı difüzörü içermektedir. Şekil 2’deki hesaplama alanı ise tek bir periyodik çark geçidi ile buna karşılık gelen kanatlı difüzör geçidini göstermektedir.

Reynolds sayısı nasıl tanımlanmıştır?

İrtifalar arasında viskoz etkileri karşılaştırmak için referans Reynolds sayısı şu bağıntıyla hesaplanmıştır:

\[ Re=\frac{\rho u_2 b_2}{\mu} \]

Burada \(\rho\) giriş havasının yoğunluğunu, \(u_2\) çark çıkışındaki çevresel hızı, \(b_2\) çark çıkışı kanat yüksekliğini ve \(\mu\) havanın dinamik viskozitesini göstermektedir. Çalışmada \(b_2=17\) mm alınmıştır.

Yoğunluk ideal gaz bağıntısıyla hesaplanmıştır:

\[ \rho=\frac{p}{RT} \]

Basınç ve sıcaklık irtifayla birlikte değiştiği için yoğunluk hızla azalmıştır. 25 km’deki yoğunluk deniz seviyesinin yaklaşık %3,2’sine, Reynolds sayısı ise yaklaşık %4’üne düşmüştür.

Neden aynı gerçek dönme hızı kullanılmamıştır?

Farklı irtifalarda giriş sıcaklığı değiştiği için kompresörler doğrudan aynı mekanik devirde karşılaştırılsaydı çevresel Mach sayıları aynı olmayacaktı. Bu nedenle araştırmacılar aynı düzeltilmiş dönme hızını korumuş ve gerçek mekanik hızı her irtifanın giriş sıcaklığına göre değiştirmiştir.

Düzeltilmiş dönme hızı:

\[ N_{\mathrm{cor}} = \frac{N}{\sqrt{T^*/T_0}} \]

Düzeltilmiş kütlesel debi ise:

\[ \dot{m}_{\mathrm{cor}} = \dot{m} \frac{\sqrt{T^*/T_0}}{P^*/P_0} \]

biçiminde ifade edilmiştir. Burada \(T_0=288{,}15\) K ve \(P_0=101.325\) Pa referans koşullarıdır.

Bu düzenleme, çalışmanın “aynı gerçek motor devri” değil, aynı aerodinamik düzeltilmiş hız koşullarını karşılaştırdığı anlamına gelmektedir.

Sayısal model nasıl kurulmuştur?

Üç boyutlu, sürekli, sıkıştırılabilir Reynolds ortalamalı Navier–Stokes denklemleri ANSYS CFX ile çözülmüştür. Türbülans için SST \(k-\omega\) modeli kullanılmıştır. Bu model, duvara yakın akışları ve ters basınç gradyanı altındaki ayrılmaları çözmeye yönelik bir yapıya sahiptir.

Düşük Reynolds sayısında laminer-türbülanslı geçişin değişimini hesaba katmak için SST modeli, \(\gamma-Re_\theta\) geçiş modeliyle birleştirilmiştir. Model iki ek büyüklüğü izlemektedir:

• Aralıklılık katsayısı \(\gamma\): Sıfıra yaklaştığında laminer veya geçiş hâlindeki, bire yaklaştığında büyük ölçüde türbülanslı akışı ifade etmektedir.
• Momentum kalınlığı Reynolds sayısı \(Re_\theta\): Sınır tabakası geçişinin başlaması ve gelişmesiyle ilişkilidir.

Girişte toplam basınç, toplam sıcaklık ve eksenel akış yönü; difüzör çıkışında statik basınç tanımlanmıştır. Karakteristik eğrileri oluşturmak için çıkış geri basıncı kademeli olarak yükseltilmiştir. Bütün katı yüzeyler kaymasız ve adyabatik kabul edilmiştir.

Çark ve difüzör arasında veri nasıl aktarılmıştır?

Dönen çark alanı ile sabit difüzör alanı arasında Frozen Rotor arayüzü kullanılmıştır. Bu yaklaşım çark ve difüzörün göreli çevresel konumunu sabit tutarak hız, basınç ve sıcaklık düzensizliklerini arayüzden geçirir. Böylece çark çıkışındaki çevresel düzensizliğin bir kısmı difüzör girişine aktarılabilmektedir.

Bununla birlikte Frozen Rotor ve sürekli RANS yaklaşımı aşağıdaki süreçleri fiziksel zaman içinde çözmemektedir:

• Kanat geçiş frekansına bağlı rotor–stator etkileşimi,
• uç kaçağı girdabının zamana bağlı hareketi ve parçalanması,
• dönen stall hücreleri,
• surge salınımları,
• kararlılık sınırı yakınındaki geçici akış yapıları.

Bu nedenle çalışmadaki düşük debili son yakınsamış noktalar gerçek surge sınırını temsil etmemektedir. Araştırmacıların kullandığı “hesaplanabilen en yüksek basınç oranı”, yalnızca sürekli çözümlerin yakınsadığı aralıktaki en yüksek değerdir.

Ağ bağımsızlığı nasıl değerlendirilmiştir?

Çark ve difüzör için yapılandırılmış ağlar ANSYS TurboGrid ile üretilmiştir. Ön kenar, arka kenar, uç açıklığı, emme yüzeyi, difüzör girişi ve difüzör kanadı ön kenarı gibi yüksek gradyan beklenen bölgeler yerel olarak inceltilmiştir.

Şekil 3’te çark yüzeylerindeki \(y^+\) değerlerinin büyük ölçüde birin altında kaldığı görülmektedir. Bu, geçiş ve SST modellerinin duvara yakın çözüm gereksinimini desteklemektedir.

Şekil 4 ve Tablo 2’ye göre toplam hücre sayısı 494.000’e ulaştıktan sonra daha sık ağların basınç oranına etkisi sınırlı kalmıştır. Bu nedenle diğer hesaplamalarda 233.000 çark ve 261.000 difüzör hücresinden oluşan ağ kullanılmıştır.

Model NASA deneyleriyle ne kadar uyumludur?

Model, deniz seviyesindeki basınç oranı ve izentropik verim karakteristikleri bakımından NASA deney verileriyle karşılaştırılmıştır. Şekil 5 ve Şekil 6’da sayısal eğriler deneysel eğilimleri genel olarak izlemektedir.

• Toplam basınç oranındaki ortalama bağıl hata: %8,00
• İzentropik verimdeki ortalama bağıl hata: %5,56
• Basınç oranı eğrisindeki ortalama kütlesel debi sapması: %2,58
• Basınç oranı eğrisindeki en yüksek kütlesel debi sapması: %4,87
• Verim eğrisindeki ortalama kütlesel debi sapması: %2,06
• Verim eğrisindeki en yüksek kütlesel debi sapması: %5,17

En büyük noktasal farklar boğulmaya yakın dik eğri kolunda görülmüştür. Küçük bir boğulma sınırı kayması, aynı debide büyük basınç oranı ve verim farkı oluşturabilmektedir. Bu nedenle model ana çalışma bölgesindeki göreli değişimleri incelemek için kullanılmış; gerçek boğulma veya kararlılık sınırının yüksek doğrulukla tahmini için yeterli kabul edilmemiştir.

İrtifa arttığında kompresör performansı nasıl değişmiştir?

Şekil 7 ve Şekil 8, irtifa arttıkça basınç oranı ve verim karakteristiklerinin daha düşük düzeltilmiş debilere ve daha düşük performans düzeylerine kaydığını göstermektedir.

Deniz seviyesinden 25 km’ye çıkıldığında en yüksek izentropik verim 0,112 azalmış, yani yaklaşık %12,9 düşmüştür. Yüksek verimli çalışma aralığı da daralmıştır. Bu durum, düşük Reynolds sayısının yalnızca basınç üretme kapasitesini değil, kayıpların düşük tutulabildiği çalışma aralığını da sınırladığını göstermektedir.

Entropi üretimi nasıl hesaplanmıştır?

İdeal gaz için özgül entropi değişimi:

\[ \Delta s = c_p\ln\left(\frac{T_{\mathrm{out}}}{T_{\mathrm{in}}}\right) - R\ln\left(\frac{p_{\mathrm{out}}}{p_{\mathrm{in}}}\right) \]

bağıntısıyla ifade edilmiştir. Araştırmacılar kaybı karşılaştırmak için üstel boyutsuz entropi parametresini tanımlamıştır:

\[ s^* = \exp\left(\frac{\Delta s}{R}\right) = \frac{ \left(T_{\mathrm{out}}/T_{\mathrm{in}}\right)^{\kappa/(\kappa-1)} }{ p_{\mathrm{out}}/p_{\mathrm{in}} } \]

İzentropik ideal süreçte \(s^*=1\)’dir. Viskoz yayınım, şoklar, ayrılma ve karışma arttıkça değer birin üzerine çıkar.

Üstel parametre doğrusal olarak toplanamadığı için difüzör kaybı logaritmik entropi artışı üzerinden ayrıştırılmıştır:

\[ \ln s^*_{\mathrm{diff}} = \ln s^*_{\mathrm{stage}} - \ln s^*_{\mathrm{imp}} \]

Kademe boyutsuz entropi parametresi 1,246’dan 1,412’ye çıkarak %13,3 artmıştır. Çark parametresi %7,9 yükselmiştir. Difüzör parametresi önce 10 km’de hafifçe azalmış, ardından 20 ve 25 km’de belirgin biçimde artmıştır.

Kayıp çark ve difüzör arasında nasıl dağılmıştır?

Bileşenlerin toplam entropi artışındaki payları:

\[ C_{\mathrm{imp}} = \frac{\ln s^*_{\mathrm{imp}}} {\ln s^*_{\mathrm{stage}}} \times 100 \]

\[ C_{\mathrm{diff}} = \frac{\ln s^*_{\mathrm{diff}}} {\ln s^*_{\mathrm{stage}}} \times 100 \]

bağıntılarıyla hesaplanmıştır.

Şekil 9’daki yığılmış sütunlar, deniz seviyesinde difüzör kaybının baskın olduğunu; Reynolds sayısı azaldıkça çark kaybının hızla büyüdüğünü göstermektedir. 20–25 km’de iki bileşen birbirine yakın paylara ulaşmıştır.

Entropi payı ile bileşen veriminin düşme oranı aynı kavram değildir. 10–25 km arasında:

• Kademe verimi %9,55,
• Çark verimi %4,05,
• Difüzör verimi %5,72

azalmıştır. Şekil 10’daki sonuç, difüzörün performansının bozulan giriş akışına çarktan daha duyarlı olduğunu göstermektedir.

Sınır tabakası geçişi nasıl değişmiştir?

Şekil 11’de çark geçidinin açıklığın %10, %50 ve %90 konumlarındaki türbülans aralıklılığı gösterilmektedir. Deniz seviyesinde geçidin büyük kısmı yüksek \(\gamma\) değerlerine sahiptir ve sınır tabakası erken biçimde türbülanslı duruma geçmektedir.

İrtifa yükseldikçe düşük aralıklılık bölgeleri kanat yüzeyleri boyunca ve akış yönünde genişlemektedir. En belirgin genişleme muhafazaya yakın, açıklığın %90 konumunda görülmüştür. Bu durum düşük Reynolds sayısının geçişi geciktirdiğini ve laminer veya geçiş hâlindeki bölgelerin alanını artırdığını göstermektedir.

Düşük aralıklılık alan oranı şu şekilde tanımlanmıştır:

\[ R_{\gamma<0.5} = \frac{A(\gamma<0.5)} {A_{\mathrm{span}}} \times100 \]

Şekil 12’ye göre geçiş baskılanması yalnızca kanat ucunda oluşmamış, 25 km’de açıklığın tamamına yayılmıştır. Bununla birlikte 20–25 km arasındaki en hızlı artış muhafazaya yakın kesitte gerçekleşmiştir.

Düşük hızlı akış bölgeleri nasıl genişlemiştir?

Şekil 13, muhafazaya yakın span = 0,9 kesitindeki Mach sayısı dağılımlarını göstermektedir. Deniz seviyesinde düşük hızlı akış büyük ölçüde duvar çevresinde ve geçit çıkışında sınırlı kalmıştır. 20 km’de emme yüzeyinin orta ve arka bölümlerinde daha sürekli bir düşük hızlı şerit gelişmiş; 25 km’de bu şerit arka kenardaki düşük hızlı alanla birleşmiştir.

Düşük Mach sayılı alan oranı:

\[ R_{Ma<0.4} = \frac{A(Ma<0.4)} {A_{\mathrm{span}}} \times100 \]

olarak tanımlanmıştır.

Şekil 14’te muhafazaya yakın kesitin bütün irtifalarda en yüksek düşük hızlı alan oranına sahip olduğu görülmektedir. 20 km sonrasında bu kesitte alan oranı yaklaşık sabit kalmış, ancak akışın biçimi ve ayrılma yapısı bozulmaya devam etmiştir. Bu nedenle yalnız alan yüzdesi akışın şiddetini tam olarak açıklamamaktadır.

Yerel ayrılma nerede gelişmiştir?

Şekil 15’te span = 0,9 kesitinin akım çizgileri karşılaştırılmıştır. Deniz seviyesinde çizgiler büyük ölçüde düzenli ve süreklidir. İrtifa yükseldikçe geçidin orta ve arka bölümündeki akım çizgileri emme yüzeyine doğru toplanmış ve yanal sapma artmıştır.

Şekil 16’daki büyütülmüş karşılaştırmada 0 km’de çizgiler kavisli olmakla birlikte kapalı bir geri dolaşım yapısı oluşturmamaktadır. 25 km’de ise muhafaza tarafındaki emme yüzeyi yakınında akım çizgileri kıvrılarak yerel girdap benzeri bir ayrılma ve geri dolaşım bölgesi meydana getirmiştir.

Bu bulgu, mekanizma zincirinin ilk üç aşamasını desteklemektedir:

1. Türbülanslı geçişin gecikmesi,
2. Duvara yakın momentum taşınımının zayıflaması,
3. Ters basınç gradyanı altında düşük momentumlu akışkanın birikerek ayrılması.

Üç boyutlu ana akış ve girdap sistemi nasıl değişmiştir?

Şekil 17’deki üç boyutlu akım çizgileri, deniz seviyesinde ana akışın geçit boyunca daha düzenli olduğunu göstermektedir. 20 ve 25 km’de düşük hızlı bölge genişlemiş, akım çizgileri emme yüzeyi ve muhafazaya doğru daha fazla sapmış ve yerel dolaşıklık artmıştır.

Şekil 18’de Mach sayısı alanı ile Q-kriteri girdap çekirdekleri birlikte gösterilmiştir. Deniz seviyesinde girdap yapıları daha çok kanat ucu ve muhafaza yakınında, nispeten sürekli şeritler hâlindedir. Yüksek irtifada girdap çekirdekleri geçidin iç kısmına doğru yayılmış ve daha dağınık bir geometri oluşturmuştur.

Normalize edilmiş girdap çekirdeği alanı:

\[ A_Q^* = \frac{A_Q} {A_{Q,0\mathrm{km}}} \]

bağıntısıyla hesaplanmıştır.

Şekil 19, girdap çekirdeği izoyüzeylerinin uzamsal kapsamının Reynolds sayısı azaldıkça büyüdüğünü göstermektedir. Bu ölçüt girdabın dönme şiddetini veya vortisite büyüklüğünü değil, aynı Q eşiğinde tanımlanan izoyüzey alanını karşılaştırmaktadır. Ayrıca sonuçlar sürekli RANS çözümündeki zaman ortalamalı girdap topolojisini temsil etmektedir.

Çark içindeki bozulma difüzöre nasıl aktarılmıştır?

Çark çıkışı ile difüzör girişi arasındaki akış açısı, radyal ve çevresel hız bileşenlerinden hesaplanmıştır. Pozitif açı normal dışa doğru radyal taşınımı, negatif açı ise radyal hız bileşeninin tersine döndüğü anormal sapma veya yerel geri akışı göstermektedir.

Şekil 20’de 25 km koşulunda negatif açılı bölgenin muhafaza boyunca belirgin biçimde genişlediği görülmektedir. Alan oranı 6,60 yüzde puan artarken en düşük açı çok daha kuvvetli bir ters sapmaya ulaşmıştır.

Kanatlı difüzörde giriş akış açısı, gelen akış ile difüzör kanadı ön kenarının uyumunu belirler. Çevresel ve açıklık yönünde farklılaşan giriş açıları yerel hücum açısını, tıkanmayı, ayrılmayı ve karışma kaybını artırabilir. Bu nedenle çark içindeki düşük hızlı ayrılmanın etkisi çark çıkışında sona ermemekte, difüzör kaybına aktarılmaktadır.

Çalışmanın ortaya koyduğu birleşik kayıp mekanizması

1. İrtifa yükseldikçe yoğunluk ve Reynolds sayısı azalır.
2. Kanat yüzeyindeki laminer-türbülanslı geçiş gecikir.
3. Duvara yakın momentum alışverişi zayıflar.
4. Özellikle muhafaza tarafındaki emme yüzeyinde düşük hızlı akışkan birikir.
5. Ters basınç gradyanı altında yerel ayrılma ve geri dolaşım gelişir.
6. Uç kaçağı, ikincil akış ve ayrılma bölgeleri arasındaki etkileşim güçlenir.
7. Girdap çekirdekleri geçit içine yayılır ve ana akış düzeni bozulur.
8. Çark çıkışında negatif radyal hız ve akış açısı düzensizliği artar.
9. Difüzör giriş uyumu kötüleşir; ayrılma, karışma ve basınç geri kazanım kayıpları büyür.
10. Çark ile difüzör kayıpları birlikte artarak kademe verimini ve çalışma aralığını düşürür.

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

• Reynolds sayısının 5,70 × 10⁵ değerinden 2,28 × 10⁴ değerine düşmesi, sayısal modelde basınç oranını, verimi ve debi kapasitesini azaltmıştır.
• En yüksek izentropik verim 0,87’den 0,76’ya düşmüştür.
• Kademe boyutsuz entropi parametresi %13,3 artmıştır.
• Deniz seviyesinde difüzör kaybı baskınken yüksek irtifada çark kaybının payı büyümüştür.
• 25 km’de geçiş baskılanması çark açıklığının tamamında belirgin hâle gelmiştir.
• Muhafazaya yakın düşük Mach sayılı alan oranı yaklaşık %30’dan %42’ye çıkmıştır.
• Muhafaza tarafındaki emme yüzeyinde yerel ayrılma ve geri dolaşım gelişmiştir.
• Normalize edilmiş girdap çekirdeği alanı 1,000’dan 1,206’ya yükselmiştir.
• Difüzör girişindeki negatif akış açılı alan %9,32’den %15,92’ye çıkmıştır.
• En düşük giriş açısı −20,3°’den −64,5°’ye düşmüştür.
• Çark çıkışı akış bozulması ile artan difüzör kaybı arasında tutarlı bir mekanizma ilişkisi kurulmuştur.

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

• Modelin 20 veya 25 km’de gerçek bir kompresörün mutlak performansını aynı doğrulukla tahmin ettiği deneysel olarak gösterilmemiştir.
• Karakteristik eğrilerin düşük debili son noktaları gerçek surge sınırı değildir.
• Düşük Reynolds sayısının zamana bağlı stall hücrelerini ve rotor–stator etkileşimini nasıl değiştirdiği çözülmemiştir.
• Girdap çekirdeği alanındaki artış, girdap şiddetinin aynı oranda arttığı anlamına gelmemektedir.
• Difüzör kaybındaki bütün artışın yalnızca negatif giriş açılarından kaynaklandığı doğrudan kanıtlanmamıştır.
• Farklı kompresör geometrilerinin aynı kayıp paylarını veya eşik irtifalarını göstereceği belirlenmemiştir.
• Gerçek uçuş giriş bozulmaları, nem, parçacıklar, buzlanma ve sıcaklık düzensizlikleri incelenmemiştir.
• Kanat veya muhafaza ısı transferi modellenmemiştir.
• Yapısal gerilme, yorulma, rotor dinamiği ve mekanik güvenlik değerlendirilmemiştir.
• Önerilen hassas bölgeler için yeni bir optimize edilmiş geometri veya deneysel akış kontrolü geliştirilmemiştir.

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

Araştırmanın güçlü yönü, genel performans eğrilerini tek başına yorumlamak yerine kademe, çark ve difüzör kayıplarını aynı modelde birlikte incelemesidir. Entropi artışının bileşenlere ayrılması, düşük Reynolds sayısının yalnızca difüzörü veya yalnızca çarkı etkilediği biçimindeki tek bileşenli açıklamaları sınırlandırmaktadır.

İkinci güçlü yön, kayıp artışının birbirini tamamlayan çok sayıda akış göstergesiyle izlenmesidir. Geçiş aralıklılığı, düşük Mach sayılı alan, akım çizgileri, Q-kriteri girdap çekirdekleri ve difüzör giriş açısı aynı mekanizma zincirinde değerlendirilmiştir.

Üçüncü güçlü yön, NASA deneyleriyle yapılan temel doğrulamanın yalnızca görsel eğri benzerliğine bırakılmaması; basınç oranı, verim ve kütlesel debi sapmalarının sayısal olarak raporlanmasıdır.

Başlıca sınırlılıklar

Model yüksek irtifalarda deneysel olarak doğrulanmamıştır. 0 km doğrulamasındaki ortalama basınç oranı hatası %8,00 ve verim hatası %5,56’dır. Bu hata düzeyi, yüksek irtifadaki küçük performans farklarının mutlak hassasiyetle yorumlanmasını sınırlandırmaktadır.

Çalışma sürekli RANS ve tek geçitli Frozen Rotor modeline dayanmaktadır. Zamana bağlı kanat geçiş etkileri, tam halka düzensizlikleri, stall başlangıcı ve surge dinamiği doğrudan çözülememektedir.

NASA CC3 geometrisi yayımlanmış koordinatlardan yeniden kurulmuştur. Gerçek deney donanımıyla küçük geometrik farklılıklar, uç açıklığı belirsizlikleri ve sınır koşulları doğrulama farklarına katkı sağlayabilir.

Bütün duvarlar adyabatik kabul edilmiştir. Yüksek irtifa uçuşunda duvar sıcaklığı ve ısı transferinin geçiş davranışına etkisi bu çalışma kapsamında değerlendirilmemiştir.

Giriş koşulları standart atmosfer değerlerine dayanmaktadır. Gerçek hava aracının hava alığı, manevrası, sınır tabakası yutması veya çevresel giriş bozulması modele eklenmemiştir.

Metnin giriş bölümündeki bazı literatür atıfları “Error! Reference source not found.” biçiminde bozulmuştur. Bu durum çalışma bulgularını doğrudan değiştirmemekle birlikte kaynak izlenebilirliği ve yayın hazırlığı açısından düzeltilmesi gereken bir biçimlendirme sorunudur.

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

Sayısal yöntem özeti

Ana performans değişimleri

Şekillerin bilimsel mesajları

• Şekil 1: NASA CC3 kompresörünün üstten, önden ve eğik geometrik görünümlerini göstermektedir.
• Şekil 2: Ana kanat, ara kanat, uç açıklığı, periyodik yüzeyler ve kanatlı difüzörü içeren hesaplama alanını sunmaktadır.
• Şekil 3: Çark duvarlarındaki \(y^+\) değerlerinin yakın duvar çözümü için düşük tutulduğunu göstermektedir.
• Şekil 4: Toplam hücre sayısı yaklaşık 494.000’e ulaştıktan sonra basınç oranının ağ inceltmesine sınırlı duyarlılık gösterdiğini ortaya koymaktadır.
• Şekil 5: Sayısal ve deneysel basınç oranı karakteristiklerini, ana çalışma ve boğulmaya yakın bölgeleri karşılaştırmaktadır.
• Şekil 6: Sayısal ve deneysel izentropik verim eğrilerinin ana çalışma bölgesindeki uyumunu göstermektedir.
• Şekil 7: İrtifayla basınç oranı eğrilerinin aşağı ve daha düşük debilere kaydığını göstermektedir.
• Şekil 8: En yüksek verimin ve yüksek verimli çalışma aralığının irtifayla azaldığını göstermektedir.
• Şekil 9: Entropi artışının çark ile difüzör arasındaki dağılımını ve irtifayla değişen katkı oranlarını göstermektedir.
• Şekil 10: Kademe, çark ve difüzör verimlerinin 10–25 km arasında düştüğünü; en büyük bileşen verim kaybının difüzörde oluştuğunu göstermektedir.
• Şekil 11: Açıklığın üç konumunda düşük aralıklılık bölgelerinin irtifayla kanat yüzeyi boyunca genişlediğini göstermektedir.
• Şekil 12: Düşük aralıklılık alan oranının özellikle 20 km sonrasında hızla arttığını göstermektedir.
• Şekil 13: Muhafazaya yakın kesitte düşük Mach sayılı alanın emme yüzeyi ve geçit çıkışı boyunca büyüdüğünü göstermektedir.
• Şekil 14: Düşük hızlı alan oranının bütün açıklık konumlarında arttığını ve en yüksek değerin muhafaza yakınında kaldığını göstermektedir.
• Şekil 15: İrtifayla akım çizgilerinin düzeninin bozulduğunu ve emme yüzeyi yönünde yanal sapmanın arttığını göstermektedir.
• Şekil 16: Deniz seviyesinde bulunmayan belirgin yerel geri dolaşımın 25 km’de geliştiğini göstermektedir.
• Şekil 17: Çark içindeki üç boyutlu ana akışın yüksek irtifada daha dağınık ve düzensiz hâle geldiğini göstermektedir.
• Şekil 18: Düşük hızlı alanlarla Q-kriteri girdap çekirdeklerinin daha fazla örtüştüğünü ve girdapların geçit içine yayıldığını göstermektedir.
• Şekil 19: Normalize edilmiş girdap çekirdeği alanının irtifayla 1,000’dan 1,206’ya yükseldiğini göstermektedir.
• Şekil 20: 25 km’de difüzör girişindeki negatif akış açılı bölgelerin daha geniş ve daha şiddetli olduğunu göstermektedir.

Yöntemsel yorum

Çalışmanın bulguları, düşük Reynolds sayısı etkisinin basit bir verim düzeltme katsayısı olarak değerlendirilmemesi gerektiğini göstermektedir. Viskozite etkisi önce sınır tabakasını, ardından geçit içindeki üç boyutlu akış yapısını ve son olarak bileşen eşleşmesini değiştirmektedir.

Özellikle muhafaza tarafındaki emme yüzeyi, çark çıkışı ve difüzör girişi kritik bölgeler olarak belirlenmiştir. Gelecekteki tasarım çalışmaları yalnızca çarkın tek başına verimini artırmak yerine şu amaçları birlikte değerlendirmelidir:

• Muhafazaya yakın düşük hızlı alanı sınırlandırmak,
• Emme yüzeyi arka bölgesindeki ayrılmayı geciktirmek,
• uç kaçağı girdaplarının geçit içine yayılımını azaltmak,
• çark çıkış akış açısını difüzör kanatlarıyla daha iyi eşleştirmek,
• çark ve difüzörü bütünleşik biçimde optimize etmek.

Kaynak ve Yöntem Notu

Çalışmanın tam özgün adı:Low-Reynolds-number effects on a transonic centrifugal compressor under near-space conditions: Performance degradation and coupled impeller–diffuser loss evolution

Yazarlar ve sıraları: Shaohua Lv; Yunxiao Wang; Shen Gao; Siwei Li; Yanxing Zhao; Maoqiong Gong.

Eş birinci yazar veya eş katkı bilgisi: Çalışmada yer almamaktadır.

Sorumlu yazar: Yanxing Zhao.

Kurumsal bağlantılar:

1. State Key Laboratory of Cryogenic Science and Technology, Technical Institute of Physics and Chemistry, Chinese Academy of Sciences, Pekin, Çin.
2. School of Advanced Interdisciplinary Sciences, University of Chinese Academy of Sciences, Pekin, Çin.

DOI: 10.2139/ssrn.7201739.

Yayın platformu: SSRN.

Dergi: Çalışmada kabul edilmiş veya yayımlanmış bir hakemli dergi adı yer almamaktadır.

Özgün yayınevi: Hakemli nihai yayın bulunmadığından doğrulanmış bir makale yayınevi yoktur. SSRN, Elsevier bünyesinde işletilen bir preprint ve erken araştırma paylaşım platformudur.

Yayın yılı: 2026.

Kaynak türü: NASA CC3 deney verileriyle deniz seviyesinde doğrulanan üç boyutlu, sürekli RANS ve geçiş modeli tabanlı sayısal araştırma preprinti.

Hakemlik durumu: Çalışma hakem değerlendirmesinden geçmemiştir. Bulgular erken araştırma sonuçları olarak değerlendirilmelidir.

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

Yazar katkıları: Shaohua Lv ilk taslak, doğrulama, yazılım ve veri düzenleme; Yunxiao Wang inceleme-düzenleme, biçimsel analiz ve kavramsallaştırma; Shen Gao inceleme-düzenleme ve yöntem; Siwei Li inceleme-düzenleme ve yöntem; Yanxing Zhao inceleme-düzenleme, kaynaklar, biçimsel analiz ve kavramsallaştırma; Maoqiong Gong danışmanlık ve kaynaklar görevlerini üstlenmiştir.

Fon desteği: Çalışma Çin Ulusal Doğa Bilimleri Vakfının 52322602 ve 52227811 numaralı destekleriyle finanse edilmiştir.

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

Bu Türkçe açıklama yüklenen çalışmanın 65 sayfası; başlık ve yazar bilgileri, denklemler, altı tablo, Şekil 1–20, eklenen yüksek çözünürlüklü akış alanları, performans grafikleri, ağ yapıları, aralıklılık dağılımları, Mach sayısı alanları, akım çizgileri, girdap çekirdekleri ve difüzör giriş açısı haritaları incelenerek hazırlanmıştır. Bilimsel içerik yalnızca çalışmanın yöntemleri, sayısal sonuçları ve kendi değerlendirmelerine dayandırılmıştır.

Dış kaynak kontrolü yalnızca özgün başlık, yazarlar, SSRN kaydı, DOI ve preprint durumunun bibliyografik doğrulanması amacıyla kullanılmıştır. Çalışmada bulunmayan dış aerodinamik sonuçlar açıklamaya eklenmemiştir.

Başlıca yöntemsel sınırlılıklar; yüksek irtifa koşullarında bağımsız deneysel doğrulama bulunmaması, sürekli RANS ve Frozen Rotor yaklaşımının zamana bağlı rotor–stator etkileşimini çözememesi, gerçek stall ve surge sınırlarının belirlenememesi, tek geçitli periyodik model kullanılması, duvarların adyabatik kabul edilmesi ve gerçek uçuş giriş bozulmalarının hesaba katılmamasıdır.

Çalışmadaki “hesaplanabilen en yüksek basınç oranı” gerçek surge sınırındaki basınç oranı olarak yorumlanmamalıdır. Benzer şekilde normalize edilmiş Q-kriteri alanı girdap çekirdeklerinin uzamsal kapsamını ifade etmekte, girdap şiddetini doğrudan ölçmemektedir.

Giriş bölümündeki bazı literatür bağlantıları “Error! Reference source not found.” biçiminde bozulmuştur. Bu kaynak biçimlendirme sorunu, nihai hakemli sürüm hazırlanırken düzeltilmelidir.