Akademik araştırmalar, anlaşılır dil

Verianla | Akademik Araştırmalardan Türkçe Ekonomi ve Bilim İçerikleri

27 Eylül 2026, Pazar
VERİANLABağımsız bilim yayıncılığı
Menüyü aç veya kapat
...
Home / Uygulamalı Bilimler / Mühendislik / Hipersonik Yüzey Sürtünmesi Ölçümlerinde Çift Hata Kontrolü: Sensör Aralığı ve Montaj Sapmaları Nasıl Düzeltiliyor?
Mühendislik

Hipersonik Yüzey Sürtünmesi Ölçümlerinde Çift Hata Kontrolü: Sensör Aralığı ve Montaj Sapmaları Nasıl Düzeltiliyor?

Bu çalışma, hipersonik akışlarda duvar kayma gerilmesini doğrudan ölçmek için kullanılan yüzer elemanlı yüzey sürtünme terazilerindeki iki temel hata kaynağını aynı yöntem içinde ele almaktadır.

30/07/2026  Veri Anla 30 görüntüleme
Hipersonik Yüzey Sürtünmesi Ölçümlerinde Çift Hata Kontrolü: Sensör Aralığı ve Montaj Sapmaları Nasıl Düzeltiliyor?

Bu çalışma, hipersonik akışlarda duvar kayma gerilmesini doğrudan ölçmek için kullanılan yüzer elemanlı yüzey sürtünme terazilerindeki iki temel hata kaynağını aynı yöntem içinde ele almaktadır. Araştırmacılar, yüzer eleman ile çevresindeki duvar arasındaki 0,3 mm yapısal aralığın oluşturduğu sistematik hatayı 20–40 km irtifa ve Mach 5–7 aralığındaki 25 sayısal akış koşulunda modellemiş; montaj sırasında oluşan X ve Y yönlü eş eksenlilik sapmaları ile H yönündeki yüzey hizası sapmasını Tip B belirsizlik yöntemiyle değerlendirmiştir. En uygun sistematik hata modeli ikili dördüncü dereceden polinom olarak belirlenmiş, önerilen ±50 µm montaj toleransı altında tüm uçuş zarfındaki en yüksek göreli genişletilmiş belirsizlik ±%10,46 olarak hesaplanmıştır. Bununla birlikte çalışma hakem değerlendirmesinden geçmemiş bir preprinttir; sonuçlar gerçek bir hipersonik rüzgâr tünelinde aynı sensör üzerinde deneysel olarak doğrulanmamış ve düzeltme polinomunun katsayıları yayımlanan metinde verilmemiştir.

Model sonuçları, ölçüm elemanının çevresindeki aralığın akışı boşluğa çekerek çevrintili bir kavite akışı oluşturduğunu ve yüzer elemanın kenarında halka biçiminde yüksek duvar kayma gerilmesi meydana getirdiğini göstermektedir. 25 km irtifa ve Mach 6,5 koşulunda 0,3 mm aralık, ölçüm noktasında duvar kayma gerilmesini 40,38 Pa artırarak belirgin bir fazla tahmine yol açmıştır. Montaj sapmaları arasında en kritik bileşen ise yüzer elemanın duvar düzleminin üstünde veya altında kalmasını temsil eden ΔH sapmasıdır; bu bileşenin toplam rastgele belirsizliğe katkısı incelenen koşullarda %94,44 ile %98,46 arasında değişmiştir.

Araştırmanın temel sorusu nedir?

Araştırmanın temel sorusu, hipersonik yüzey sürtünme terazilerindeki yapısal aralık hatası ile montaj sapmalarından doğan ölçüm belirsizliğinin tek bir hassasiyet yönetimi çerçevesinde kontrol edilip edilemeyeceğidir. Çalışma iki farklı hata türünü birbirinden ayırmaktadır:

  • Sistematik hata: Yüzer ölçüm elemanı ile sabit duvar arasındaki yapısal aralığın akış alanını sürekli ve belirli bir yönde değiştirmesi.
  • Rastgele hata veya belirsizlik: Sensör montajında ortaya çıkan X, Y ve H yönlü konum sapmalarının ölçüm sonucunu değiştirmesi.

Sistematik hata, uygun bir düzeltme değeri çıkarılarak telafi edilebilir. Montaj kaynaklı sapmalar ise her kurulumda aynı değerde olmayabileceğinden, çalışmada doğrudan “düzeltilen” bir hata olarak değil, toleransla sınırlandırılan ve belirsizlik bütçesine aktarılan bir bileşen olarak ele alınmıştır.

Duvar kayma gerilmesi neden önemlidir?

Duvar kayma gerilmesi (wall shear stress, WSS), akışkanın araç yüzeyine teğetsel yönde uyguladığı sürtünme kuvvetinin birim alana düşen değeridir ve pascal cinsinden ifade edilir. Hipersonik araçlarda bu büyüklüğün doğru belirlenmesi; sürtünme sürüklemesi, yüzey ısınması, sınır tabakası davranışı, menzil tahmini ve sürükleme azaltma tasarımları açısından önem taşır.

Dolaylı ölçüm yöntemleri hız profili veya ampirik bağıntılardan WSS çıkarmaya çalışırken, yüzer elemanlı sürtünme terazisi yüzey üzerindeki teğetsel kuvveti doğrudan mekanik yapı üzerinden algılar. Bununla birlikte ölçüm elemanının hareket edebilmesi için çevresinde bırakılan küçük boşluk, ölçülmek istenen akışın kendisini değiştirir. Bu nedenle sensör, yalnızca akışı ölçen pasif bir eleman değil, yerel akış alanına müdahale eden bir geometridir.

Literatürde hangi boşluk hedeflenmiştir?

Önceki çalışmalar aralık büyüklüğünün, yüzey hizası sapmasının veya eş eksenlilik hatasının etkilerini çoğunlukla ayrı ayrı incelemiştir. Bazı araştırmalar yüzer elemanın duvarın üzerinde kalmasının ölçüm hatasını %100’ün üzerine çıkarabildiğini, bazı çalışmalar ise aralık genişledikçe kavite emişinin ve ölçüm sapmasının arttığını bildirmiştir.

Bu çalışma, belirli bir hipersonik sürtünme terazisi geometrisi için iki eksikliği hedeflemektedir. Birincisi, yapısal aralığın oluşturduğu sistematik hatayı yalnızca tek bir akış koşulunda değil, tam bir irtifa–Mach zarfında düzeltmektir. İkincisi, X ve Y yönlü eş eksenlilik sapmaları ile H yönündeki yüzey hizası sapmasını birlikte değerlendirerek montaj toleransı ve genişletilmiş belirsizlik belirlemektir.

Sayısal rüzgâr tüneli nasıl kurulmuştur?

Üç boyutlu hesaplama alanının akış doğrultusundaki uzunluğu 5000 mm, genişliği ve yüksekliği ise 1000 mm olarak belirlenmiştir. Sürtünme terazisi giriş sınırından 2500 mm ileride ve montaj plakasının ön kenarından 500 mm geride konumlandırılmıştır. Bu yerleşim, sınır tabakasının sensöre ulaşmadan önce gelişmesi amacıyla seçilmiştir.

Yüzer elemanın çapı 20 mm, terazi yapısının toplam yüksekliği 70 mm ve eleman ile duvar arasındaki aralık 0,3 mm’dir. Sistematik hata hesabında aynı geometri iki biçimde modellenmiştir: fiziksel aralığın bulunduğu 0,3 mm’lik model ve ideal aralıksız model. Bu iki çözüm arasındaki WSS farkı, aralık kaynaklı düzeltme miktarı olarak kullanılmıştır.

Montaj sapmaları üç doğrultuda tanımlanmıştır:

  • ΔX: Yüzer elemanın akış doğrultusundaki yatay konum sapması.
  • ΔY: Akışa dik yatay doğrultudaki eş eksenlilik sapması.
  • ΔH: Elemanın duvar yüzeyinin üstünde veya altında kalması; başka bir ifadeyle yüzey hizası sapması.

Çalışmanın bazı bölümlerinde ΔH için “flatness”, bazı bölümlerinde “flushness” terimi kullanılmaktadır. Şemalar ve fiziksel açıklamalar, bu değişkenin yüzer elemanın duvar düzlemine göre düşey konumunu temsil ettiğini göstermektedir.

Hangi uçuş koşulları incelenmiştir?

Temel uçuş zarfı 20, 25, 30, 35 ve 40 km irtifalar ile Mach 5, 5,5, 6, 6,5 ve 7 hızlarının birleşiminden oluşan 25 koşulu kapsamaktadır. Atmosferik basınç ve sıcaklık değerleri NRLMSIS 2.0 atmosfer modelinden alınmıştır.

İrtifaStatik basınçStatik sıcaklık
20 km5357,25 Pa211,6 K
25 km2433,95 Pa219,1 K
30 km1138,74 Pa231,1 K
35 km544,61 Pa235,1 K
40 km269,47 Pa243,2 K

Duvar kaymama koşuluyla ve 300 K sabit sıcaklıkta tanımlanmış, çıkış sınırı basınç çıkışı olarak belirlenmiştir. Çözümler kararlı durumlu örtük bir çözücüyle yürütülmüş, sıcaklığa bağlı viskozite Sutherland yasasıyla hesaplanmıştır.

Türbülans ve ağ yapısı nasıl ele alınmıştır?

Hipersonik sınır tabakası ve WSS dağılımı, k–ω SST türbülans modeliyle hesaplanmıştır. Duvara en yakın ilk ağ tabakasının yüksekliği 2 × 10−5 m olarak seçilmiş ve boyutsuz duvar uzaklığı yaklaşık \(y^+ \approx 1\) düzeyinde tutulmuştur. Bu değer, duvar yakınındaki hız gradyanının duvar fonksiyonlarına aşırı bağımlı olmadan çözülmesi amacıyla kullanılmıştır.

Aralık bölgesindeki keskin hız ve basınç değişimlerini yakalamak için yapılandırılmamış ağ ve yerel ağ inceltmesi uygulanmıştır. Ek Şekil S1’de ağ hücresi sayısı artırıldığında aralıklı ve aralıksız modellerdeki WSS değerlerinin belirli bir yoğunluktan sonra küçük değişimler gösterdiği görülmektedir. Araştırmacılar, bu plato bölgesindeki ağı hesaplama maliyeti ile sonuç kararlılığı arasında uygun çözüm olarak seçmiştir.

Referans sıcaklık yöntemiyle doğrulama

Aralıksız modelin WSS tahmini, yüksek hızlı sınır tabakaları için kullanılan referans sıcaklık yöntemiyle karşılaştırılmıştır. Çalışmada duvar veya geri kazanım sıcaklığı şu şekilde hesaplanmaktadır:

\[ T_{\omega}=T_e\left(1+r\frac{\gamma-1}{2}Ma^2\right) \]

Burada \(T_{\omega}\) duvar veya geri kazanım sıcaklığını, \(T_e\) serbest akış statik sıcaklığını, \(r\) geri kazanım katsayısını, \(\gamma\) özgül ısı oranını ve \(Ma\) Mach sayısını temsil eder. Çalışmada \(r\) için yaklaşık 0,85, hava için \(\gamma\) değerinde yaklaşık 1,4 kullanılmıştır.

Referans sıcaklık bağıntısı şöyledir:

\[ \frac{T^*}{T_e}=1+0.039Ma^2+0.5\left(\frac{T_{\omega}}{T_e}-1\right) \]

\(T^*\), yoğunluk ve viskozite gibi özelliklerin değerlendirildiği referans sıcaklıktır. Referans Reynolds sayısı şu şekilde tanımlanmıştır:

\[ Re_L^*=\frac{\rho^*U_eL}{\mu^*} \]

\(\rho^*\) referans sıcaklıktaki yoğunluk, \(U_e\) serbest akış hızı, \(L\) karakteristik uzunluk ve \(\mu^*\) referans sıcaklıktaki dinamik viskozitedir.

Referans sürtünme katsayısı:

\[ C_f^*=\frac{0.664}{\sqrt{Re_L^*}} \]

Çalışmada WSS için kullanılan son ifade ise şöyledir:

\[ \tau_{\omega}=\frac{1}{2}\frac{\rho^*}{\rho}\frac{\mu^*}{\mu}C_f^*\rho U_e^2 \]

Burada \(\tau_{\omega}\) pascal cinsinden duvar kayma gerilmesini, \(\rho\) ve \(\mu\) gerçek akış yoğunluğu ile dinamik viskoziteyi temsil eder.

25 km irtifa ve Mach 6 koşulunda sayısal model 62,94 Pa, referans sıcaklık yöntemi ise 64,93 Pa WSS üretmiştir. Aradaki yaklaşık %3,05 fark, aralıksız model için temel doğrulama göstergesi olarak sunulmuştur.

Sürtünme katsayısı karşılaştırması

Farklı serbest akış dinamik basınçlarının etkisini azaltmak için boyutsuz yüzey sürtünme katsayısı kullanılmıştır:

\[ C_f=\frac{\tau_{\omega}}{\frac{1}{2}\rho U_e^2} \]

25 km ve Mach 6 koşulunda hesaplanan değer \(C_f=1.56\times10^{-3}\) olmuştur. Çalışmanın ek tablosunda benzer türbülanslı hipersonik koşullar için bildirilen sonuçlar \(1.16\times10^{-3}\), \(1.69\times10^{-3}\) ve yaklaşık \(1.7\times10^{-3}\) düzeyindedir. Mevcut hesap bu aralığın içinde kalmaktadır.

Makalede Denklem 7 altında \(q_e\) değişkeni “dinamik basınç” olarak tanımlanmasına rağmen yazılan matematiksel ifade, yukarıdaki \(C_f\) formülünü aynen tekrar etmektedir:

\[ q_e=\frac{\tau_{\omega}}{\frac{1}{2}\rho U_e^2} \]

Bu ifade boyutsuz bir oran üretirken dinamik basıncın biriminin pascal olması gerekir. Bu nedenle Denklem 7 ile onu izleyen metin arasında açık bir sembol veya denklem tutarsızlığı bulunmaktadır. Çalışmayı sessizce düzeltmek yerine, bu teknik sorun sonuçların yorumunda dikkate alınmalıdır.

Aralık sistematik hatayı nasıl oluşturmuştur?

25 km ve Mach 6,5 koşulundaki akım çizgileri, ana akışın bir bölümünün basınç farkı nedeniyle 0,3 mm aralığa çekildiğini göstermektedir. Akış, yüzer elemanın altındaki kavitede çevrintili bir yapı oluşturmakta ve karşı taraftan tekrar ana akışa katılmaktadır.

Bu kavite akışı, elemanın çevresinde halka biçiminde yüksek WSS bölgesi üretmiştir. Şekil 2’de halkanın en yüksek kesimleri kırmızı ve sarı renklerle, elemanın merkezindeki daha düşük gerilme alanı ise mavi tonlarla gösterilmektedir. Aralıklı ve aralıksız çözüm karşılaştırıldığında ölçüm konumundaki fark 40,38 Pa’dır. Bu pozitif fark, sensörün gerçek duvar sürtünmesini olduğundan yüksek ölçmesine neden olmaktadır.

Sistematik hata nasıl tanımlanmıştır?

Aralık kaynaklı mutlak hata şu şekilde tanımlanmıştır:

\[ \Delta WSS=WSS_{\mathrm{meas}}-WSS_{\mathrm{true}} \]

\(WSS_{\mathrm{meas}}\), 0,3 mm aralıklı sensör modelinden; \(WSS_{\mathrm{true}}\) ise ideal aralıksız modelden elde edilen değerdir. Araştırmacılar hata oranı yerine pascal cinsinden mutlak farkı modellemiştir. Bunun nedeni, gerçek WSS’nin irtifa ve Mach sayısıyla büyük ölçüde değişmesi ve küçük gerçek değerlerde göreli hatanın yapay biçimde büyüyebilmesidir.

Neden dördüncü dereceden polinom seçilmiştir?

İrtifa ve Mach sayısının aralık hatası üzerindeki doğrusal olmayan ortak etkisini modellemek için ikili üçüncü, dördüncü ve beşinci dereceden polinomlar karşılaştırılmıştır.

ModelTerim sayısıR²Düzeltilmiş R²RMSEAICBIC
İkili üçüncü derece100,99140,98621,5987 Pa43,4655,65
İkili dördüncü derece150,99710,99310,9212 Pa25,9044,18
İkili beşinci derece210,93220,59324,4804 Pa116,99142,58

Dördüncü dereceden model en yüksek R² ve düzeltilmiş R² değerlerini, en düşük RMSE’yi ve en düşük bilgi kriterlerini üretmiştir. Beşinci dereceden modelde terim sayısı artmasına rağmen düzeltilmiş R²’nin 0,5932’ye düşmesi ve RMSE’nin 4,4804 Pa’ya çıkması, aşırı uyum işareti olarak değerlendirilmiştir.

Model, ana 25 koşula ek olarak 22 km–Mach 5,2; 26 km–Mach 5,8; 33 km–Mach 6,2 ve 37 km–Mach 6,7 koşullarındaki dört doğrulama noktasıyla sınanmıştır. Şekil 3’te doğrulama noktaları yüzeye yakın, artıklar ise sıfır çevresinde gösterilmektedir.

Bununla birlikte yüklenen sürüm, ikili dördüncü dereceden polinomun 15 katsayısını vermemektedir. Dolayısıyla okuyucu yalnızca irtifa ve Mach sayısını kullanarak bildirilen düzeltme yüzeyini bağımsız biçimde yeniden hesaplayamaz. Sunulan başarı ölçütleri model seçimini desteklese de katsayı eksikliği yeniden üretilebilirliği sınırlar.

Montaj sapmaları için ölçüm modeli

Montaj kaynaklı değişimler şu genel modelle ifade edilmiştir:

\[ \tau=f(X,Y,H)=\tau_0+C_XX+C_YY+C_HH+C_{XY}XY+C_{XYH}XYH \]

\(\tau_0\), montaj sapması olmayan durumdaki WSS’dir. \(C_X\), \(C_Y\) ve \(C_H\), ilgili doğrultulardaki birim konum değişimine karşı WSS duyarlılıklarını temsil eder. \(C_{XY}\) ile \(C_{XYH}\) ise sapmaların birlikte bulunması durumundaki etkileşim terimleridir.

Metin, bütün etkileşim katsayılarının birimini “Pa/mm” olarak belirtmektedir. Ancak \(XY\) terimi mm², \(XYH\) terimi mm³ boyutundadır. Bu nedenle \(C_{XY}\) ve \(C_{XYH}\) katsayılarının boyutlarının ayrı ayrı açıklanması gerekir; mevcut sürümde bu boyutsal ayrım verilmemiştir.

Tip B standart belirsizlik nasıl hesaplanmıştır?

Her montaj sapmasının \(-a\) ile \(+a\) arasında eş olasılıklı, düzgün dağılım gösterdiği varsayılmıştır. Standart belirsizlik:

\[ u(X)=u(Y)=u(H)=\frac{a}{\sqrt{3}} \]

Burada \(a\), mikrometre cinsinden montaj toleransının yarı genişliğidir. Örneğin ±50 µm toleransta \(a=50\) µm alınır.

Duyarlılık katsayıları kısmi türevlerle tanımlanmıştır:

\[ C_X=\frac{\partial \tau}{\partial X},\qquad C_Y=\frac{\partial \tau}{\partial Y},\qquad C_H=\frac{\partial \tau}{\partial H} \]

\[ C_{XY}=\frac{\partial^2\tau}{\partial X\partial Y} \]

\[ C_{XYH}=\frac{\partial^3\tau}{\partial X\partial Y\partial H} \]

Tek yönlü standart belirsizlik bileşenleri:

\[ u_X(\tau)=|C_X|u(X),\qquad u_Y(\tau)=|C_Y|u(Y),\qquad u_H(\tau)=|C_H|u(H) \]

Etkileşim bileşenleri ise:

\[ u_{XY}(\tau)=|C_{XY}|u(X)u(Y) \]

\[ u_{XYH}(\tau)=|C_{XYH}|u(X)u(Y)u(H) \]

Birleşik standart belirsizlik kareler toplamının kareköküyle hesaplanmıştır:

\[ u_c(\tau)=\sqrt{u_X^2+u_Y^2+u_H^2+u_{XY}^2+u_{XYH}^2} \]

Yaklaşık %95 kapsama düzeyi için \(k=2\) alınmış ve genişletilmiş belirsizlik:

\[ U=k\,u_c(\tau) \]

Sonuç şu biçimde raporlanmıştır:

\[ \tau=\tau_0\pm U\qquad(k=2) \]

30 km ve Mach 6 koşulunda hangi sapma en etkilidir?

ΔH sapması, X ve Y sapmalarından çok daha güçlü etki göstermiştir. Pozitif ΔH, yüzer elemanı duvarın üzerine çıkararak sınır tabakasının ön kenarında bir çıkıntı oluşturur. Akış bu çıkıntı çevresinde hızlanır ve ölçülen WSS belirgin biçimde artar. Negatif ΔH ise elemanı duvar düzleminin altına indirerek küçük bir çukur oluşturur; yerel yavaşlama ve çevrinti WSS’yi azaltabilir.

Negatif ΔX sapması, elemanın ön kenarını akışa göre ileri taşıdığı için sınır tabakasını daha erken bozar ve pozitif ΔX’e göre daha büyük hata oluşturur. Y yönündeki sapmanın etkisi, geometrik simetri nedeniyle pozitif ve negatif yönlerde aynı kabul edilmiştir.

Tek yönlü hatayı yaklaşık ±%10 içinde tutmak amacıyla bütün doğrultular için ±50 µm montaj toleransı önerilmiştir. Bu koşulda 30 km ve Mach 6 için genişletilmiş belirsizlik ±4,94 Pa, göreli genişletilmiş belirsizlik ise ±%9,63 olarak hesaplanmıştır.

Sapmalar gerçekten bağımsız kabul edilebilir mi?

30 km ve Mach 6 koşulunda X–Y eş eksenlilik etkileşiminin en yüksek bağlaşım gücü %8,75; eş eksenlilik ile H yönlü yüzey hizası arasındaki en yüksek bağlaşım ise %6,61 olarak bildirilmiştir. Araştırmacılar bu seviyeleri zayıf bağlaşım olarak yorumlamış ve belirsizlik bileşenlerini doğrusal duyarlılıklara dayalı kareler toplamı yöntemiyle birleştirmiştir.

Bu sonuç, etkileşimin sıfır olduğu anlamına gelmez. Çalışmanın iddiası, incelenen tolerans ve uçuş koşullarında etkileşim terimlerinin ana tek yönlü etkilerden küçük olduğudur. Özellikle daha büyük montaj hatalarında veya farklı sensör geometrilerinde aynı varsayımın geçerli olacağı gösterilmemiştir.

Tam uçuş zarfındaki belirsizlik

Önerilen ±50 µm tolerans; düşük irtifa–yüksek hız, orta koşul ve yüksek irtifa–düşük hız olmak üzere üç temsilî durumda ayrıntılı olarak değerlendirilmiştir.

Koşulİrtifa ve hızHesaplanan WSSGenişletilmiş belirsizlikGöreli belirsizlik
120 km, Mach 7204,67 Pa±21,41 Pa±%10,46
230 km, Mach 651,32 Pa±4,94 Pa±%9,63
340 km, Mach 518,06 Pa±1,52 PaYaklaşık ±%8,42

Bütün koşullarda tek yönlü montaj hatasının en yüksek değeri %8,83 olmuştur. En yüksek toplam göreli belirsizlik, dinamik basıncın daha yüksek ve sınır tabakasının daha ince olduğu 20 km–Mach 7 koşulunda ortaya çıkmıştır.

Sınır koşullarındaki bağlaşım gücü düşük irtifa–yüksek hız durumunda en fazla %15,41, yüksek irtifa–düşük hız durumunda ise en fazla %9,52 olarak bildirilmiştir. Düşük irtifa ve yüksek hızda küçük geometrik sapmaların yerel akış alanını daha güçlü değiştirmesi nedeniyle bağlaşım artmıştır.

Belirsizlik bütçesine hangi yön hâkimdir?

Uçuş koşuluΔX katkısıΔY katkısıΔH katkısı
20 km, Mach 7%2,40%3,16%94,44
30 km, Mach 6%1,73%0,53%97,74
40 km, Mach 5%1,49%0,05%98,46

Bu dağılım, montaj sırasında en kritik kontrolün yüzer elemanın duvarla aynı düzlemde tutulması olduğunu göstermektedir. X ve Y eş eksenlilik ayarları ihmal edilmemelidir; ancak incelenen geometri ve tolerans aralığında toplam belirsizliğin büyük bölümü ΔH bileşeninden kaynaklanmıştır.

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

  • Sistematik aralık hatası ile montaj belirsizliğini tek bir ölçüm doğruluğu çerçevesinde birleştirmiştir.
  • 20–40 km ve Mach 5–7 aralığında 25 koşulu kapsayan parametrik bir uçuş zarfı kullanılmıştır.
  • Aralıklı ve aralıksız geometriler doğrudan karşılaştırılmıştır.
  • Duvar yakınında \(y^+\approx1\) sağlayan ağ, yerel ağ inceltmesi ve ağ bağımsızlığı incelemesi uygulanmıştır.
  • Sayısal sonuçlar hem referans sıcaklık yöntemi hem de yayımlanmış türbülanslı rüzgâr tüneli ve CFD sonuçlarıyla karşılaştırılmıştır.
  • Polinom derecesi yalnızca R² ile değil; düzeltilmiş R², RMSE, AIC ve BIC ölçütleriyle seçilmiştir.
  • Belirsizlik bileşenlerinin yönlere göre katkısı ayrı ayrı hesaplanarak montaj önceliği belirlenmiştir.

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

  • Bu çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.
  • Önerilen sistem aynı sensör üzerinde fiziksel hipersonik rüzgâr tüneli deneyiyle doğrulanmamıştır.
  • Düzeltme modeli yalnızca 0,3 mm aralığa, 20 mm yüzer elemana ve çalışmada kullanılan özel geometriye aittir.
  • İkili dördüncü dereceden polinomun katsayıları verilmediğinden düzeltme fonksiyonu bağımsız olarak yeniden kullanılamamaktadır.
  • Denklem 7’de dinamik basınç tanımı ile yazılan matematiksel ifade birbiriyle uyuşmamaktadır.
  • Çarpım terimlerinin duyarlılık katsayıları için verilen birimler boyutsal olarak ayrıntılandırılmamıştır.
  • Duvar 300 K’de izotermal kabul edilmiştir; aerodinamik ısınma ile değişen duvar sıcaklığının ölçüm hatasına etkisi değerlendirilmemiştir.
  • Kararlı RANS ve k–ω SST modeli kullanılmıştır; geçiş, şok–sınır tabakası etkileşimi, zamana bağlı titreşimler ve model-form belirsizliği ayrı bir belirsizlik bileşeni olarak hesaba katılmamıştır.
  • Tip B hesap, montaj sapmalarının düzgün dağılım gösterdiği varsayımına dayanır; gerçek imalat ve kurulum dağılımları ölçülmemiştir.
  • Tekrarlı fiziksel ölçümlerden elde edilen Tip A istatistiksel belirsizlik değerlendirilmemiştir.
  • Literatür karşılaştırmaları benzer ancak tamamen aynı sensör, geometri ve serbest akış koşullarından oluşmamaktadır.

Çalışma neyi desteklemektedir?

Çalışma, belirli bir yüzer elemanlı hipersonik sürtünme terazisinde 0,3 mm çevresel aralığın WSS’yi sistematik biçimde fazla tahmin ettirebildiğini desteklemektedir. Simülasyonlar ayrıca ±50 µm montaj toleransının incelenen uçuş zarfında tek yönlü montaj hatalarını yaklaşık ±%10 içinde tutabildiğini ve belirsizlik bütçesinin ana bileşeninin ΔH yüzey hizası sapması olduğunu göstermektedir.

Çalışma aynı zamanda, sistematik hata düzeltmesinin irtifa ve Mach sayısına bağlı doğrusal olmayan bir yüzey olarak modellenebileceğini ve rastgele montaj belirsizliğinin Tip B yaklaşımıyla ölçüm sonucuna eklenebileceğini göstermektedir.

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

Araştırma, önerilen ±50 µm toleransın bütün yüzey sürtünme sensörleri için evrensel olduğunu kanıtlamamaktadır. Sonuçlar farklı eleman çaplarına, aralık büyüklüklerine, kavite geometrilerine, laminer akışlara veya Mach 7 dışındaki koşullara doğrudan aktarılamaz.

Çalışma ayrıca gerçek bir rüzgâr tüneli montajında toplam ölçüm hatasının kesin olarak ±%10,46’dan küçük kalacağını göstermemektedir. Bildirilen değer; sayısal model, düzgün dağılım varsayımı, belirlenen duyarlılık katsayıları ve yalnızca incelenen montaj bileşenleri altında hesaplanan genişletilmiş belirsizliktir. Sensör elektroniği, sıcaklık sürüklenmesi, titreşim, kalibrasyon, malzeme davranışı ve tekrarlanabilirlik gibi başka hata kaynakları mevcut hesapta yer almamaktadır.

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

Çalışmanın yöntemi, Türkiye’de hipersonik rüzgâr tünelleri, yüksek hızlı aerodinamik test sistemleri, mikro kuvvet sensörleri, MEMS ölçüm elemanları ve füze veya uzay aracı yüzey testleri geliştiren araştırma grupları açısından aktarılabilir bir yaklaşım sunmaktadır. En önemli aktarılabilir sonuç belirli bir ±50 µm değeri değil, sensör geometrisinden doğan sistematik hatanın sayısal olarak haritalanması ve montaj toleransının belirsizlik bütçesiyle birlikte belirlenmesidir.

Yerel bir sensör tasarımında aynı yöntemin uygulanabilmesi için Türkiye’de geliştirilen gerçek sensör geometrisinin, gerçek aralık değerinin ve rüzgâr tüneli çalışma zarfının yeniden modellenmesi gerekir. Sayısal düzeltme yüzeyi, fiziksel kalibrasyon deneyleriyle doğrulandıktan sonra ölçüm yazılımına veya test sonrası veri işleme zincirine aktarılabilir.

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

Yöntem bileşeniUygulanan koşulTemel sonuçYorum sınırı
Hesaplama alanı5000 × 1000 × 1000 mm üç boyutlu alanSensör girişten 2500 mm, plaka ön kenarından 500 mm gerideBelirli düz plaka ve sensör yerleşimine aittir
Sensör geometrisi20 mm yüzer eleman, 70 mm toplam yükseklik, 0,3 mm aralıkAralıklı ve ideal aralıksız modeller karşılaştırıldıDiğer aralık ve eleman boyutlarına doğrudan genellenemez
Akış zarfı20–40 km, Mach 5–7, toplam 25 temel koşulİrtifa ve Mach sayısına bağlı hata yüzeyi oluşturulduBu zarfın dışı test edilmedi
Akış modeliKararlı RANS, k–ω SST, Sutherland viskozite yasasıHipersonik türbülanslı sınır tabakası ve WSS hesaplandıGeçiş ve zamana bağlı etkiler ayrıca modellenmedi
Duvar ve çıkışKaymama, 300 K izotermal duvar; basınç çıkışıStandartlaştırılmış sayısal rüzgâr tüneli kurulduAerodinamik ısınan duvar davranışı temsil edilmedi
Ağ çözünürlüğüİlk tabaka 2 × 10−5 m, \(y^+\approx1\), yerel inceltmeAğ bağımsızlığı sonrasında kararlı WSS değerleriSayısal model-form hatasını ortadan kaldırmaz
Temel doğrulama25 km, Mach 662,94 Pa CFD; 64,93 Pa referans yöntem; %3,05 farkAynı sensörde fiziksel deney değildir
Aralık etkisi25 km, Mach 6,5; 0,3 mm aralıkÖlçüm noktasında +40,38 Pa sistematik farkBu sayı incelenen koşula özeldir
Düzeltme modeliİkili dördüncü dereceden polinomR² 0,9971; düzeltilmiş R² 0,9931; RMSE 0,9212 PaPolinom katsayıları metinde verilmemiştir
Montaj toleransıΔX, ΔY ve ΔH için ±50 µmTek yönlü hata bütün zarfta en fazla %8,83Düzgün dağılım ve özel geometri varsayımına bağlıdır
Orta koşul belirsizliği30 km, Mach 651,32 ± 4,94 Pa; göreli ±%9,63Yalnızca ele alınan montaj bileşenlerini kapsar
En yüksek belirsizlik20 km, Mach 7204,67 ± 21,41 Pa; göreli ±%10,46Düşük irtifa ve yüksek hız koşuluna aittir
Hata önceliğiÜç uçuş koşulunda yön katkılarıΔH katkısı %94,44–%98,46Başka sensör geometrilerinde değişebilir

Temel mühendislik çıkarımları

  1. Yüzer elemanın çevresindeki aralık, yalnızca mekanik bir zorunluluk değil, ölçüme doğrudan müdahale eden bir akış geometrisidir.
  2. Aralık hatası, tek bir sabit yüzdeyle değil, irtifa ve Mach sayısına bağlı bir düzeltme yüzeyiyle ele alınmalıdır.
  3. Montaj sırasında en sıkı kontrol, yüzer elemanın duvarla aynı düzlemde kalmasını sağlayan ΔH yönünde uygulanmalıdır.
  4. Sayısal düzeltme modeli, gerçek rüzgâr tüneli kalibrasyonu ve açık polinom katsayıları olmadan doğrudan üretim ölçüm sistemine aktarılmamalıdır.
  5. Toplam ölçüm bütçesine elektronik, termal, kalibrasyon ve tekrarlanabilirlik belirsizlikleri de eklenmelidir.

Kaynak ve Yöntem Notu

Çalışmanın özgün adı: An Integrated Dual-Error Control Method for Hypersonic Skin Friction Balance: Gap-Induced Systematic Correction and Installation Uncertainty Evaluation

Yazarlar ve doğru sıralama: Xiaohan Liao, Hongxiao Chao, Xuefeng Chen, Diao Yang, Xiaosong Wang, Xiaohu Zhang, Yong Zhou, Yong He, Yongcai Guo.

Eş birinci yazar: Çalışmada eş birinci yazarlık veya ortak ilk yazarlık belirtilmemiştir. Yazar katkı tablosundaki “equal” ifadeleri belirli katkı rollerini göstermekte, eş birinci yazarlık beyanı oluşturmamaktadır.

Sorumlu yazarlar: Yong He ve Yongcai Guo.

Kurumsal bağlantılar:

  • Key Laboratory of Optoelectronic Technology and Systems of the Education Ministry of China, College of Optoelectronic Engineering, Chongqing University, Chongqing, Çin.
  • Northwest Institute of Mechanical and Electrical Engineering, Çin.

DOI:10.2139/ssrn.6952639

Yayın veya preprint platformu: SSRN.

Resmî kaynak:SSRN çalışma sayfası

Dergi: Hakemli bir dergi yayını doğrulanamamıştır.

Özgün yayınevi: Ayrı bir hakemli dergi yayınevi bulunmamaktadır. Çalışma SSRN preprint platformunda yer almaktadır; nihai yayımlanmış sürüm doğrulanamamıştır.

Yayın yılı: 2026.

Hakemlik durumu: Bu çalışma bir preprinttir ve hakem değerlendirmesinden geçmemiştir.

Kaynak türü: Üç boyutlu sayısal akış modellemesi, sistematik hata düzeltmesi ve Tip B belirsizlik değerlendirmesi içeren mühendislik preprinti.

Finansman: National Natural Science Foundation of China tarafından 92471207 numaralı destek ve Fundamental Research Funds for Central Universities tarafından 2024CDJGF020 numaralı destek kapsamında finanse edilmiştir.

Çıkar çatışması: Yazarlar mali çıkar çatışması bulunmadığını bildirmiştir.

Veri durumu: Yazarlar çalışmayı destekleyen verilerin makale ve ek dosyalarda bulunduğunu belirtmiştir. Bununla birlikte yüklenen sürümde seçilen dördüncü dereceden düzeltme polinomunun katsayıları ve ayrıntılı çözücü girdi dosyaları yer almamaktadır.

Bu Türkçe açıklama, yüklenen 33 sayfalık çalışma; metin, denklemler, şekiller, tablolar, ağ bağımsızlığı incelemesi, doğrulama karşılaştırmaları ve belirsizlik sonuçları temel alınarak hazırlanmıştır. Çalışmada bulunmayan bilimsel bulgular eklenmemiştir. Dış kontrol yalnızca başlık, DOI, SSRN kaydı ve yayın kimliğinin doğrulanması amacıyla kullanılmıştır.

Çalışmanın başlıca sınırları; hakem değerlendirmesinden geçmemiş olması, gerçek sensör üzerinde hipersonik deneysel doğrulama içermemesi, polinom katsayılarının verilmemesi, Denklem 7’deki dinamik basınç tutarsızlığı, etkileşim katsayılarının birimlerinin eksik açıklanması ve belirsizlik hesabının yalnızca seçilen montaj sapmalarını kapsamasıdır.

Sonuçlar, aynı sensör geometrisi ve aynı hesaplama zarfı için yöntemsel değer taşımaktadır. Önerilen düzeltme modeli ve ±50 µm tolerans, farklı sürtünme terazilerine veya gerçek uçuş sistemlerine deneysel doğrulama olmadan doğrudan aktarılmamalıdır.


Paylaş:

Yorumlar incelendikten sonra yayımlanır.Gönderdiğiniz yorum onay sürecine alınır ve uygun bulunduğunda görünür hâle gelir.

Bir yorum bırakın

E-posta adresiniz yayınlanmayacaktır. Gerekli alanlar * ile işaretlenmiştir

Your experience on this site will be improved by allowing cookies Cookie Policy