
Bu araştırma, elektro-hidrolik frenlerin bırakma sürecini gerçekleştiren elektro-hidrolik itici (electro-hydraulic thruster, EHT) içindeki akışın nasıl girdaplara dönüştüğünü ve bu girdapların akışkan viskozitesi ile motor hızından nasıl etkilendiğini geçici hesaplamalı akışkanlar dinamiği (CFD) ile incelemektedir. Araştırmacılar EHT içinde α, β ve γ olarak adlandırılan üç temel dolaşım yapısı tanımlamıştır. Model sonuçlarına göre α ve β girdapları sıvının piston altına etkili biçimde taşınmasına engel olan ve ek enerji kaybına yol açan dolaşımlar; γ girdabı ise pistonun uzamasını destekleyen temel dolaşım alanıdır. Viskozitenin 54,5 Pa·s'den 101,5 Pa·s'ye yükselmesiyle γ girdabının tepe hızı 0,39 m/s'den 0,21 m/s'ye ve piston altı ortalama basınç 0,178 MPa'dan 0,096 MPa'ya düşmüştür. Motor hızının 200 rad/s'den 500 rad/s'ye yükselmesi ise γ girdabının tepe hızını 0,09 m/s'den 0,58 m/s'ye, piston altı ortalama basıncı 0,0263 MPa'dan 0,435 MPa'ya çıkarmıştır. Bununla birlikte bütün bu sonuçlar geçici CFD simülasyonuna dayanmaktadır; çalışma fiziksel prototip deneyi veya saha doğrulaması sunmamaktadır.
Elektro-hidrolik frenin başlangıç durumu frenleme durumudur. EHT motoru çalışmadığında piston geri çekilmiş konumdadır ve yay sistemi fren kollarını tambura doğru kuvvetlendirir. Fren bırakılacağı zaman EHT motoru pervaneyi döndürür; sıvı pervane ve kılavuz plaka çevresinde dolaştırılarak piston altına basınç oluşturulur. Piston yukarı doğru uzadığında yay kuvveti yenilir ve fren kolları açılarak fren tamburu serbest bırakılır.
Çalışmanın temel bulgusu, EHT içindeki bütün girdapların aynı işlevi taşımadığıdır. α girdabı pervane üstü ile kılavuz plaka altı arasında; β girdabı kılavuz plaka ve piston altı bölgesinde akışın bir bölümünü geri dolaştırarak ana yukarı yönlü aktarımı zorlaştırmaktadır. Buna karşılık γ girdabı dış bölgelerden piston altına yükselen ve merkezden tekrar pervaneye dönen ana çevrimdir. Araştırmacılar bu nedenle γ girdabının korunmasını EHT'nin fren bırakma performansı açısından temel bir gereksinim olarak yorumlamaktadır.
Viskozite yükseldiğinde bütün girdapların mutlak tepe hızları azalmaktadır; fakat faydalı γ girdabının zayıflaması daha güçlüdür. 54,5 Pa·s'de α/γ hız oranı %94,9 ve β/γ oranı %64,1 iken 101,5 Pa·s'de bu oranlar sırasıyla %166,7 ve %104,8'e çıkmıştır. Bu artış, zararlı girdapların mutlak olarak güçlendiğini değil, γ girdabının daha hızlı zayıflaması nedeniyle göreli hâkimiyetlerinin arttığını göstermektedir.
Motor hızının artırılması ise üç girdabı da güçlendirmiştir. Ancak 200–500 rad/s aralığında γ girdabının tepe hızındaki artış kaynakta %544,4 olarak hesaplanırken α ve β için artışlar sırasıyla %169,6 ve %173,3'tür. Bu nedenle yüksek hız ana tahrik dolaşımını kuvvetlendirebilmekte, fakat aynı zamanda enerji kaybı oluşturan α ve β girdaplarını da büyütmektedir. Çalışma bu yüzden en yüksek motor hızının otomatik olarak en iyi işletme noktası olduğunu savunmamaktadır.
Türkiye açısından çalışma; vinç, kaldırma ve hidrolik fren kullanılan endüstriyel mekanizmalarda düşük sıcaklıkta yağ viskozitesi ile motor hızının birlikte dikkate alınması gerektiğine ilişkin bir sayısal tasarım çerçevesi sunmaktadır. Ancak bildirilen viskozite, hız ve basınç değerleri tek bir EHT geometrisi ve belirli CFD varsayımlarına aittir. Farklı bir fren veya itici tasarımında bu değerlerin doğrudan ayar noktası olarak kullanılması yerine ilgili sistemin kendi geometrisi ve akışkanı üzerinde deneysel veya doğrulanmış sayısal analiz yapılması gerekir.
Elektro-hidrolik fren nasıl çalışıyor?
Çalışmada incelenen elektro-hidrolik fren (electro-hydraulic brake, EHB), sabit vinç tipi bir kapak kaldırma sisteminin frenleme ve fren bırakma mekanizmasına uygulanmıştır. Sistem; elektro-hidrolik itici, yay grubu, ayarlanabilir bağlantı çubuğu, aktarım plakası, fren tamburu, destek elemanları ve iki fren kolundan oluşmaktadır.
Frenleme durumunda EHT motoru durur ve piston geri çekilir. Gergin yay, aktarım mekanizmasını aşağı çekerek iki fren kolunu tambura doğru kapatır. Böylece tambur kaldırma milini kavrayarak mekanik frenleme oluşturur.
Fren bırakma sırasında EHT motoru çalışır. Pervane sıvıya kinetik enerji verir ve sıvı pistonun altına taşınarak yukarı yönlü hidrolik kuvvet meydana getirir. Bu kuvvet yeterli olduğunda piston yay kuvvetini yenerek uzar ve fren kolları iki yana açılır.
EHT'nin içinde sıvı hangi yolu izliyor?
EHT'nin üç boyutlu modeli üst gövde, piston çubuğu, piston, kılavuz plaka, yağ haznesi, pervane, alt gövde ve motor bileşenlerinden oluşmaktadır.
Pervane döndüğünde alt bölgedeki sıvı hem çevresel olarak dönmekte hem radyal yönde dışarı doğru hareket etmektedir. Kılavuz plaka kanatları arasından yukarı çıkan sıvı piston altına ulaştığında piston üzerinde hidrolik basınç oluşturur. Akış daha sonra merkezi açıklıktan aşağıya, pervane üstüne dönerek kapalı bir dolaşım meydana getirir.
Bu akış yalnız tek yönlü bir pompalama hareketi değildir. Çevresel dönme, radyal hareket ve eksenel akış aynı anda bulunduğundan kılavuz plaka ve pervane çevresinde yerel geri dolaşımlar oluşmaktadır.
α, β ve γ girdaplarının farkı nedir?
Araştırmada üç farklı girdap yapısı tanımlanmıştır:
| Girdap | Temel bölge | Akış üzerindeki rolü |
|---|---|---|
| α | Pervane üstü ile kılavuz plaka altı | Akışın bir bölümünü geri dolaştırır; piston yönündeki aktarımı sınırlar ve enerji kaybıyla ilişkilendirilir. |
| β | Kılavuz plaka boyunca ve piston altı bölgesine doğru | Kılavuz kanatların akışı bloke etmesi ve ayrılmasıyla ilişkilendirilen yerel geri dolaşımdır; ana akışı engeller. |
| γ | Pervane–kılavuz plaka–piston altı ana dolaşım hattı | Piston altına akışı ve basıncı taşıyan yararlı ana çevrim olarak tanımlanır. |
Bu sınıflandırma, çalışmanın en önemli mekanik sonucudur. Girdap kavramı burada otomatik olarak “istenmeyen türbülans” anlamına gelmemektedir. Aynı sistem içinde bazı dolaşımlar hidrolik enerji kaybına, başka bir dolaşım ise işlevsel piston basıncının kurulmasına katkıda bulunmaktadır.
Üç girdap nasıl tespit edildi?
Araştırmacılar kılavuz plaka çevresindeki basınç ve hız dağılımını incelemek için üç eksenel izleme hattı tanımlamıştır.
- İzleme hattı 1: kılavuz plaka kanatları arasında yüksek basınçlı bölgeden geçen eksenel çizgi,
- İzleme hattı 2: EHT'nin merkezi ekseni,
- İzleme hattı 3: kanatlar arasında düşük basınçlı bölgeden geçen eksenel çizgi.
İzleme koordinatı \(Y=4\)–79 mm aralığındadır ve üç bölgeye ayrılmıştır:
- a bölgesi: \(Y=4\)–25 mm; pervane üstü ile kılavuz plaka altı,
- b bölgesi: \(Y=25\)–45 mm; kılavuz plakanın alt ve üst yüzeyi arası,
- c bölgesi: \(Y=45\)–79 mm; kılavuz plaka üstü ile piston altı.
Kılavuz plaka kanatlarının izleme hattını doğrudan kesmesini önlemek amacıyla eksenel çizgilerden ikisi üst görünüşte radyal doğrultudan 5° ötelenmiştir.
Basınç dağılımı girdapları nasıl gösteriyor?
İzleme hattı 1 ile izleme hattı 3 arasındaki basınç eğrileri a, b ve c bölgelerinde yaklaşık karşılıklı/simetrik değişim göstermektedir. Hat 1'deki basınç genel olarak hat 3'tekinden daha yüksektir. Bu fark, iki nokta arasında çevrimsel bir akışın bulunduğuna işaret etmektedir.
Merkezi izleme hattı 2'nin basıncının da hat 1'den düşük olması, dış bölgeden merkez doğrultusuna geri dolaşımı destekleyen basınç gradyanını göstermektedir.
Piston tabanına yaklaşıldığında, \(Y=79\) mm civarında üç hattın basıncı birbirine yaklaşmaktadır. Aynı noktada eksenel hızlar yaklaşık sıfıra iner. Araştırmacılar bunu akışın dinamik basıncının piston altında daha çok statik basınca dönüşmesi ve pistonun uzamasını sağlamasıyla ilişkilendirmektedir.
Hız yönleri α girdabını nasıl tanımlıyor?
a bölgesinde izleme hattı 1'deki akış negatif Y yönünde, izleme hattı 3'teki akış pozitif Y yönündedir. İki taraftaki ters yönlü eksenel hız, pervane üstü ile kılavuz plaka altında kapalı bir dolaşım meydana getirir. Çalışma bu yapıyı α girdabı olarak tanımlamaktadır.
α girdabı pervane dönmesinin oluşturduğu akış ayrışmasıyla ilişkilendirilmekte ve sıvının doğrudan piston altına taşınmasına katkı sağlamadığı için olumsuz bir yapı olarak değerlendirilmektedir.
β girdabı nerede oluşuyor?
b ve c bölgelerinde yönler tersine dönmektedir. İzleme hattı 1'de sıvı pozitif Y yönünde ilerlerken hat 3'te negatif Y yönlü bir geri akış oluşmaktadır.
Araştırmacılar bu dolaşımı β girdabı olarak adlandırmakta ve kılavuz plaka kanatlarının akışı engellemesiyle ilişkilendirmektedir. Bu dolaşım da ana piston besleme akışına karşı bir yerel enerji kaybı olarak yorumlanmaktadır.
γ girdabı neden yararlı kabul ediliyor?
Merkezi izleme hattı 2'de akış a, b ve c bölgelerinde negatif Y yönünde, yani piston bölgesinden tekrar pervane doğrultusuna geri dönmektedir. Aynı anda dış bölgelerde yukarı yönlü akış bulunması, pervane–kılavuz plaka–piston arasında büyük bir kapalı çevrim oluşturur.
Çalışmada bu dolaşım γ girdabı olarak tanımlanmıştır. γ girdabı piston altında sürekli sıvı beslemesi ve basınç oluşturduğu için EHT'nin fren bırakma işlevini destekleyen temel akış yapısıdır.
CFD modeli hangi denklemlere dayanıyor?
Akış alanı için süreklilik ve momentum korunumu denklemleri kullanılmıştır. Sıkıştırılamaz akış için süreklilik denklemi kaynakta:
\[ \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} =0 \]
şeklinde verilmiştir. Burada \(u\), \(v\) ve \(w\) sırasıyla Kartezyen koordinatlardaki hız bileşenleridir.
Momentum denklemleri her hız bileşeni için genel olarak atalet, hacim kuvveti ve viskoz yayınım terimlerini birleştirmektedir. Örneğin x yönündeki denklem:
\[ \frac{\partial u}{\partial t} + u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} + w\frac{\partial u}{\partial z} = F_x+ \frac{\mu}{\rho} \left( \frac{\partial^2u}{\partial x^2} + \frac{\partial^2u}{\partial y^2} + \frac{\partial^2u}{\partial z^2} \right) \]
biçimindedir. \(F_x\) x yönündeki hacim kuvveti bileşenini, \(\mu\) dinamik viskoziteyi ve \(\rho\) akışkan yoğunluğunu ifade eder. Aynı yapı y ve z hız bileşenleri için de kullanılmıştır.
Neden RNG k-ε türbülans modeli seçildi?
Çalışmada türbülans için RNG \(k-\varepsilon\) modeli kullanılmıştır. Araştırmacılar bu modelin dönel akış, güçlü akış çizgisi eğriliği, akış ayrılması ve yüksek gerinim oranlı bölgeleri standart \(k-\varepsilon\) modeline göre daha uygun biçimde temsil edebildiğini gerekçe göstermektedir.
Türbülans kinetik enerjisinin taşıma denklemi:
\[ \rho_1\frac{\partial k}{\partial t} + \rho_1\frac{\partial}{\partial x_i}(ku_i) = \frac{\partial}{\partial x_j} \left[ \alpha_k(\mu_1+\mu_t) \frac{\partial k}{\partial x_j} \right] + G_k-\rho_1\varepsilon \]
ve türbülans disipasyon oranı denklemi:
\[ \rho_1\frac{\partial \varepsilon}{\partial t} + \rho_1\frac{\partial}{\partial x_i}(\varepsilon u_i) = \frac{\partial}{\partial x_j} \left[ \alpha_\varepsilon(\mu_1+\mu_t) \frac{\partial \varepsilon}{\partial x_j} \right] + C_{1\varepsilon}\frac{\varepsilon}{k}G_k - C_{2\varepsilon}\rho_1\frac{\varepsilon^2}{k} \]
olarak verilmiştir.
Kaynakta model katsayıları:
- \(C_{1\varepsilon}=1,42\)
- \(C_{2\varepsilon}=1,68\)
olarak kullanılmıştır.
Sayısal çözüm hangi koşullarda yapılmıştır?
Geçici CFD analizleri ANSYS Fluent 2021 R1 ortamında gerçekleştirilmiştir. Basınç tabanlı örtük çözücü ve basınç-hız eşleştirmesi için SIMPLE algoritması kullanılmıştır. Türbülans kinetik enerjisi ve disipasyon oranı taşınım terimleri ikinci mertebe upwind şemasıyla ayrıklaştırılmıştır.
| Simülasyon parametresi | Kaynakta bildirilen değer |
|---|---|
| CFD yazılımı | ANSYS Fluent 2021 R1 |
| Çözüm tipi | Geçici / transient |
| Basınç-hız eşleştirmesi | SIMPLE |
| Ayrıklaştırma | İkinci mertebe upwind |
| Artık yakınsama kriteri | \(10^{-6}\) |
| Referans yoğunluk | 866 kg/m³ |
| Referans dinamik viskozite | 76,5 Pa·s |
| Referans pervane hızı | 300 rad/s |
| Duvar koşulu | No-slip |
| Duvar y+ | < 5 |
| Zaman adımı | \(1,0\times10^{-3}\) s |
| Toplam simülasyon süresi | 0,5 s |
Sıcaklık için hangi akışkan değerleri kullanılmıştır?
| Sıcaklık | Yoğunluk (kg/m³) | Dinamik viskozite (Pa·s) |
|---|---|---|
| 0 °C | 872 | 101,5 |
| 10 °C | 866 | 76,5 |
| 20 °C | 869 | 54,5 |
Bu değerler kaynakta verilen birimler ve sayılarla korunmuştur. Çalışma düşük sıcaklık bağlamını bu üç durumla temsil etmektedir; parametrik CFD tablosu 0 °C'nin altındaki bir sıcaklık noktası içermemektedir.
Ağ bağımsızlığı nasıl kontrol edilmiştir?
EHT akış hacmi polihedral elemanlarla ağlanmıştır. Pervane, pervane alt bağlantı bölgesi ve kılavuz plaka kanatları yerel olarak inceltilmiştir.
| Eleman sayısı | Piston altı basınç (MPa) | Kaynakta verilen simülasyon hatası (%) |
|---|---|---|
| 10.043.090 | 0,0962 | 0,31 |
| 6.212.859 | 0,0968 | 0,00 |
| 5.270.888 | 0,0965 | 0,21 |
Araştırmacılar sonuç farklarının %0,35'in altında olduğunu belirterek hesaplama yükünü azaltmak amacıyla 5.270.888 elemanlı ağı nihai analizlerde kullanmıştır.
Viskozite yükseldiğinde ne oluyor?
Viskozite 54,5 Pa·s'den 101,5 Pa·s'ye yükseldiğinde üç girdabın da tepe hızı azalmaktadır. Ancak bu azalma eşit değildir.
Verianla Live: Viskozitenin üç girdaba ve piston altı basınca etkisi
Aşağıdaki değerler çalışmanın Tablo 3'ü ve Bölüm 4.1'de doğrudan bildirilen iki uç viskozite durumunu karşılaştırmaktadır. Ara değer üretilmemiştir.
| Dinamik viskozite (Pa·s) | α tepe hızı (m/s) | β tepe hızı (m/s) | γ tepe hızı (m/s) | α/γ (%) | β/γ (%) | Piston altı ortalama basınç (MPa) |
|---|---|---|---|---|---|---|
| 54,5 | 0,37 | 0,25 | 0,39 | 94,9 | 64,1 | 0,178 |
| 101,5 | 0,35 | 0,22 | 0,21 | 166,7 | 104,8 | 0,096 |
Verianla Live: Grafik yalnız yukarıdaki görünür bilimsel tablodan oluşturulur. Değerler kaynakta bildirilen iki viskozite uç durumudur.
Kaynakta verilen azalma oranları α için %5,4, β için %12,0 ve γ için %46,2'dir. Dolayısıyla yüksek viskozite zararlı kabul edilen α ve β girdaplarını da zayıflatmaktadır; fakat piston tahrikine katkı veren γ girdabını çok daha güçlü biçimde baskılamaktadır.
Bu yüzden α/γ ve β/γ oranları yükselmektedir. Göreli oranların yükselmesi α ve β'nin mutlak hızlarının arttığı anlamına gelmez; γ'nın daha hızlı zayıflaması sonucu oluşmaktadır.
Kaynağın viskozite konusunda kendi içinde çeliştiği nokta nedir?
Makalenin özetinde daha yüksek viskozitenin girdap gelişimini ve basınç dalgalanmasını artırdığı yazmaktadır. Buna karşılık Bölüm 4.1'de basınç dalgalanma genliğinin viskoziteyle azaldığı, üç girdabın tepe hızlarının da düştüğü gösterilmiştir. Tartışma bölümü de yüksek viskozitenin bütün girdapları bastırdığını belirtmektedir.
Bu nedenle sayısal tablolar ve ayrıntılı sonuç bölümü dikkate alındığında desteklenen ifade şudur: viskozite artışı mutlak girdap hızlarını azaltır, ancak faydalı γ girdabını oransal olarak daha fazla zayıflatır ve piston altı basıncı düşürür. Özet bölümündeki aksi yöndeki ifade kaynak içi tutarsızlık olarak korunmalıdır.
Motor hızının artırılması ne yapıyor?
Motor hızı 200, 300, 400 ve 500 rad/s seviyelerinde analiz edilmiştir. Hız arttıkça üç izleme hattındaki basınçlar ve basınç farkları büyümüş, bütün girdapların tepe hızları yükselmiştir.
Verianla Live: Motor hızının üç girdaba ve piston altı basınca etkisi
Tablo çalışmanın Tablo 4'ünde verilen 200 ve 500 rad/s uç değerleri ile Bölüm 4.2'deki piston altı basınç sonuçlarını bir araya getirir.
| Motor hızı (rad/s) | α tepe hızı (m/s) | β tepe hızı (m/s) | γ tepe hızı (m/s) | α/γ (%) | β/γ (%) | Piston altı ortalama basınç (MPa) |
|---|---|---|---|---|---|---|
| 200 | 0,23 | 0,15 | 0,09 | 255,6 | 166,7 | 0,0263 |
| 500 | 0,62 | 0,41 | 0,58 | 106,9 | 70,7 | 0,435 |
Verianla Live: Görselleştirme yalnız kaynakta açıkça verilen uç değerleri kullanır; 300 ve 400 rad/s için grafikten yaklaşık sayı üretilmez.
Kaynak, 200–500 rad/s arasında α girdabı tepe hızındaki artışı %169,6, β için %173,3 ve γ için %544,4 olarak bildirmektedir. Bu nedenle γ'nın göreli ağırlığı hız yükseldikçe artmaktadır: α/γ oranı %255,6'dan %106,9'a, β/γ oranı %166,7'den %70,7'ye düşmektedir.
Ancak motor hızındaki artış α ve β'yi ortadan kaldırmamaktadır. Tam tersine bu iki zararlı girdabın mutlak hızları da büyümektedir. Çalışmanın mühendislik yorumu bu nedenle “motor hızını mümkün olduğunca yükseltmek” değil, fren bırakma kuvveti ile hidrolik enerji kaybı arasında uygun bir hız aralığı belirlemektir.
Viskozite ile motor hızı birlikte nasıl yorumlanmalı?
Çalışmanın sayısal sonuçları iki farklı düzenleme mekanizması göstermektedir. Viskozite yükseldiğinde viskoz sürtünme ve akış direnci artmakta, tüm dolaşımlar zayıflamaktadır. Fakat γ daha hızlı zayıfladığı için piston tahrik kapasitesi olumsuz etkilenmektedir.
Motor hızı yükseldiğinde sisteme verilen kinetik enerji ve basınç gradyanı artmakta, üç girdap da güçlenmektedir. γ girdabının büyüme oranı daha yüksek olduğundan piston basıncı iyileşmekte; fakat α ve β kaynaklı enerji kaybı da yükselmektedir.
Kaynak bu iki parametrenin sıcaklık, yük ve tepki süresiyle uyumlu seçilmesini önermektedir. Bununla birlikte makalede viskozite ve hızın bütün kombinasyonlarını aynı deney matrisinde tarayan ayrı bir iki faktörlü etkileşim analizi verilmemiştir.
Çalışmanın desteklediği sonuçlar
- İncelenen EHT geometrisinin geçici CFD modelinde α, β ve γ olarak üç farklı dolaşım yapısı tanımlanmıştır.
- α ve β girdapları ana piston besleme akışını engelleyen ve ek hidrolik enerji kaybıyla ilişkilendirilen yapılar olarak değerlendirilmiştir.
- γ girdabı piston altındaki basıncın oluşumuna katkı sağlayan ana dolaşım olarak belirlenmiştir.
- 54,5–101,5 Pa·s viskozite aralığında üç girdabın tepe hızları azalmış, γ girdabındaki azalma kaynakta %46,2 olarak bildirilmiştir.
- Aynı viskozite değişiminde piston altı ortalama basınç 0,178 MPa'dan 0,096 MPa'ya düşmüştür.
- 200–500 rad/s motor hızında üç girdap da güçlenmiş; γ girdabının tepe hızındaki artış %544,4 olarak raporlanmıştır.
- Aynı hız aralığında piston altı ortalama basınç 0,0263 MPa'dan 0,435 MPa'ya yükselmiştir.
- Ağ bağımsızlığı kontrolünde seçilen 5.270.888 elemanlı ağın piston altı basınç sonucu daha yoğun ağlarla yakın bulunmuştur.
Çalışmanın desteklemediği veya test etmediği sonuçlar
- CFD sonuçları fiziksel EHT prototipiyle deneysel olarak doğrulanmamıştır.
- Çalışma gerçek saha koşullarında fren bırakma süresi veya güvenilirlik testi sunmamaktadır.
- RNG k-ε modelinin başka türbülans modellerine göre sonuç duyarlılığı test edilmemiştir.
- Kavitasyon etkisi modele dahil edilmemiştir.
- Akışkan içindeki sıcaklık gradyanları çözülmemiştir.
- İmalat toleransları, montaj boşlukları ve bileşen aşınması modele dahil edilmemiştir.
- Üç viskozite noktası ve 200–500 rad/s hız aralığı bütün EHT çalışma koşullarını temsil etmemektedir.
- Kaynakta tek bir evrensel optimum motor hızı verilmemiştir.
- Kaynakta bütün EHT tasarımlarına uygulanabilecek evrensel optimum viskozite belirlenmemiştir.
- Viskozite ve motor hızının tüm kombinasyonlarını içeren tam faktöriyel bir etkileşim analizi sunulmamıştır.
Çalışmanın Yöntemi ve Bulguları
EHT geometrisi ve akış yapısı
Çalışmada kullanılan EHT'nin ana bileşenleri üst gövde, piston çubuğu, piston, kılavuz plaka, yağ haznesi, pervane, alt gövde ve elektrik motorudur. Pervane sıvıyı çevresel ve radyal yönde hızlandırmakta, kılavuz plaka akışı piston altına yönlendirmekte ve merkez bölgesi dönüş hattını oluşturmaktadır.
CFD modelinin ana amacı pistonun mekanik hareketini tam bir çoklu fizik modeli olarak çözmekten çok, fren bırakma sürecindeki iç akış alanının basınç ve hız karakterini ve bunun girdaplara dönüşümünü incelemektir.
Akışkan sınır koşulları
| Çalışma sıcaklığı | Yoğunluk (kg/m³) | Dinamik viskozite (Pa·s) |
|---|---|---|
| 0 °C | 872 | 101,5 |
| 10 °C | 866 | 76,5 |
| 20 °C | 869 | 54,5 |
Kaynağın ana referans çözümünde 866 kg/m³ yoğunluk, 76,5 Pa·s viskozite ve 300 rad/s pervane hızı kullanılmaktadır.
Sayısal çözücü ayarları
| Teknik ayar | Değer/yöntem |
|---|---|
| Yazılım | ANSYS Fluent 2021 R1 |
| Akış çözümü | Geçici |
| Türbülans modeli | RNG k-ε |
| Basınç-hız bağlantısı | SIMPLE |
| Konvektif ayrıklaştırma | İkinci mertebe upwind |
| Artık yakınsama kriteri | \(10^{-6}\) |
| Pervane sınır koşulu | Dönen duvar |
| Diğer duvarlar | No-slip |
| Duvar y+ | < 5 |
| Zaman adımı | \(1,0\times10^{-3}\) s |
| Toplam süre | 0,5 s |
Ağ bağımsızlığı testi
| Karakteristik ağ boyutu | Eleman sayısı | Piston altı basınç (MPa) | Simülasyon hatası (%) |
|---|---|---|---|
| 1,2 mm | 10.043.090 | 0,0962 | 0,31 |
| 1,6 mm | 6.212.859 | 0,0968 | 0,00 |
| 2,0 mm | 5.270.888 | 0,0965 | 0,21 |
Çalışma nihai hesaplamalarda 2,0 mm karakteristik ağ boyutuna karşılık gelen 5.270.888 elemanlı modeli seçmiştir.
Viskoziteye bağlı girdap sonuçları
| Dinamik viskozite (Pa·s) | α tepe hızı (m/s) | β tepe hızı (m/s) | γ tepe hızı (m/s) | α/γ | β/γ |
|---|---|---|---|---|---|
| 54,5 | 0,37 | 0,25 | 0,39 | %94,9 | %64,1 |
| 101,5 | 0,35 | 0,22 | 0,21 | %166,7 | %104,8 |
| Kaynakta verilen değişim | −%5,4 | −%12,0 | −%46,2 | — | — |
Bu tablo, yüksek viskozitenin zararlı girdapları tamamen bastırmadığını değil, üç girdabı da zayıflattığını ve γ'nın en fazla zayıflayan dolaşım olduğunu göstermektedir.
Motor hızına bağlı girdap sonuçları
| Motor hızı (rad/s) | α tepe hızı (m/s) | β tepe hızı (m/s) | γ tepe hızı (m/s) | α/γ | β/γ |
|---|---|---|---|---|---|
| 200 | 0,23 | 0,15 | 0,09 | %255,6 | %166,7 |
| 500 | 0,62 | 0,41 | 0,58 | %106,9 | %70,7 |
| Kaynakta verilen değişim | +%169,6 | +%173,3 | +%544,4 | — | — |
Motor hızının artırılması faydalı γ dolaşımını göreli olarak daha hızlı güçlendirmiştir; ancak α ve β girdaplarının mutlak hızları da yükselmiştir. Bu nedenle sonuç yalnız yüksek hızın avantajını değil, tahrik ile enerji kaybı arasındaki bir optimizasyon problemini göstermektedir.
Piston altı basıncın iki temel parametreye tepkisi
| Değişken | Başlangıç koşulu | Başlangıç piston altı ortalama basıncı | Son koşul | Son piston altı ortalama basıncı |
|---|---|---|---|---|
| Dinamik viskozite | 54,5 Pa·s | 0,178 MPa | 101,5 Pa·s | 0,096 MPa |
| Motor hızı | 200 rad/s | 0,0263 MPa | 500 rad/s | 0,435 MPa |
Sayısal modelde yüksek viskozite piston altı basıncını düşürürken yüksek motor hızı piston altı basıncını yükseltmiştir. Bu sonuçlar yalnız incelenen EHT geometrisi ve tanımlanan CFD koşulları için geçerlidir.
Çalışmanın ana model varsayımları ve sınırları
- Geometri sabittir; imalat ve montaj toleransları hesaba katılmamıştır.
- Bileşen aşınması modellenmemiştir.
- Kavitasyon modele dahil edilmemiştir.
- Akışkan içi sıcaklık gradyanları çözülmemiştir.
- Sınır bozucularının tümü modellenmemiştir.
- Yalnız RNG k-ε türbülans modeli kullanılmıştır.
- Deneysel PIV, basınç sensörü veya prototip testiyle doğrulama yapılmamıştır.
- Duvar ve girdap çekirdeği çözünürlüğü hesaplama kaynaklarıyla sınırlıdır.
Gelecek çalışma olarak ne öneriliyor?
Araştırmacılar sonraki çalışmalarda prototip deneyleri, farklı türbülans modellerinin karşılaştırılması, uyarlanabilir ağ inceltme ve daha geniş viskozite/hız parametre aralıklarının incelenmesini önermektedir. Bu unsurlar mevcut araştırmanın gerçekleştirilmiş sonuçları değil, gelecekteki doğrulama adımlarıdır.
Kaynak ve Yöntem Notu
Tam özgün çalışma adı: Research on the Formation Mechanism of Vortices and Key Parameter Regulation in the Electro-Hydraulic Thruster
Yazarlar: Yanan Sun, Zezheng Tian, Na Li, Haiyong Jiang, Chao Yang, Chongchong Chen, Lei Yang, Lei Xing, Lijie Zhang.
Yazar sırası: Kaynaktaki sıra aynen korunmuştur.
Sorumlu yazar: Lijie Zhang.
Eş katkı/eş birinci yazar: Kaynakta böyle bir beyan bulunmamaktadır.
Kurumlar: College of Mechatronical & Electrical Engineering, Hebei Agricultural University, Baoding, China; Hebei Province Intelligent Agricultural Equipment Technology Innovation Center, Baoding, China.
Dergi: Machines.
Yayınevi: MDPI.
Bibliyografik kayıt: Machines 2026, 14, 669.
DOI: 10.3390/machines14060669.
Resmî yayın bağlantısı:https://doi.org/10.3390/machines14060669
Yayın tarihi: 8 Haziran 2026.
Kaynak türü ve hakemlik durumu: Hakemli özgün araştırma makalesidir. Bilimsel bulgular geçici CFD simülasyonlarına dayanmaktadır; fiziksel prototip doğrulaması yapılmamıştır.
Lisans: Creative Commons Attribution (CC BY).
Finansman: Çalışma; Hebei Province Innovation Team of Modern Agricultural Industry Technology System (HBCT2024050206), Hebei Province Major Science and Technology Support Program Project (242N1901Z), Hebei Province Agricultural Science and Technology Achievement Transformation Fund Project (2025JNZ-S19) ve Hebei Province üniversitelerinin temel bilimsel araştırma giderleri projesi (KY2025044) tarafından kısmen desteklenmiştir.
Veri erişilebilirliği: Özgün simülasyon verileri, ağ dosyaları ve model parametre dosyalarının makul talep üzerine sorumlu yazardan sağlanabileceği belirtilmiştir.
Çıkar çatışması: Yanan Sun, Na Li, Haiyong Jiang, Chao Yang, Chongchong Chen, Lei Yang ve Lijie Zhang'ın Hebei Province Intelligent Agricultural Equipment Technology Innovation Center bünyesinde görev yaptığı belirtilmiştir. Yazarlar çalışmanın ticari veya finansal bir ilişki olmadan yürütüldüğünü bildirmektedir.
Yazar katkıları: Yanan Sun yöntem, araştırma ve biçimsel analiz; Zezheng Tian ANSYS Fluent 2021 R1 yazılımı ve biçimsel analiz; Na Li inceleme/düzenleme ve proje yönetimi; Haiyong Jiang yöntem ve doğrulama; Chao Yang görselleştirme ve veri düzenleme; Chongchong Chen biçimsel analiz; Lei Yang doğrulama; Lei Xing kaynak sağlama; Lijie Zhang gözetim ve kavramsallaştırma görevlerini üstlenmiştir.
Bu Verianla makalesindeki EHB/EHT çalışma prensibi, akış alanı açıklamaları, denklemler, CFD koşulları, ağ bağımsızlığı verileri, girdap sınıflandırması, viskozite ve motor hızı sonuçları ile sınırlılıklar incelenen kaynak çalışmaya dayanmaktadır. Bibliyografik kimlik doğrulaması dışında dış kaynaklardan yeni bilimsel deney, mekanizma veya sayısal sonuç eklenmemiştir.
Temel yöntemsel sınır: Çalışmada fiziksel EHT prototip testi, saha ölçümü veya yayımlanmış deney verisiyle doğrulama yapılmamıştır. Yalnız RNG k-ε türbülans modeli kullanılmış ve türbülans modeli duyarlılığı incelenmemiştir. Kavitasyon, akışkan sıcaklık gradyanı, üretim toleransları, montaj boşlukları, aşınma ve bazı geçici sınır bozucuları modele dahil edilmemiştir.
Kaynak içi tutarsızlık notu: Makalenin özetinde yüksek viskozitenin girdap gelişimini ve basınç dalgalanmasını artırdığı belirtilmektedir. Buna karşılık ayrıntılı Bölüm 4.1, Tablo 3 ve Tartışma bölümünde viskozite yükseldikçe üç girdabın da tepe hızının azaldığı ve basınç dalgalanma genliğinin küçüldüğü belirtilmektedir. Verianla makalesinde bu iki ifade sessizce birleştirilmemiş; sayısal sonuçlar ile kaynak içi sözlü uyuşmazlık ayrı tutulmuştur.
Parametrik etkileşim sınırı: Kaynak özetinde viskozite ve motor hızının bağımsız ve birleşik etkilerinden söz edilmekle birlikte, sunulan temel parametrik sonuçlar iki değişkeni ayrı ayrı değiştiren analizlerden oluşmaktadır. Tam iki faktörlü bir etkileşim matrisi veya istatistiksel coupling testi sunulmamıştır.
Birim notu: Dinamik viskozite kaynakta 54,5; 76,5 ve 101,5 Pa·s olarak verilmiştir. Bu değerler kaynak sadakati gereği aynı birimle korunmuştur.

Bir yorum bırakın
E-posta adresiniz yayınlanmayacaktır. Gerekli alanlar * ile işaretlenmiştir