Bu çalışma, gelişmiş üç aşamalı hafif gaz fırlatıcısında ikincil uçucu plakanın hem çok yüksek hıza ulaştırılması hem de hedefe mümkün olduğunca düz biçimde çarpması arasındaki mühendislik çatışmasını incelemiştir. Araştırmacılar, kenar profili bir kuvvet fonksiyonuyla biçimlendirilen ve merkez tabanına yüksek yoğunluklu platin katman eklenen çok katmanlı heterojen kademeli yoğunluklu çarpıcıyı geliştirmiştir. ANSYS Autodyn 2022 R1 ile kurulan iki boyutlu eksenel simülasyon modeli, daha önce yayımlanmış üç aşamalı gaz tabancası deneyleriyle karşılaştırılmış; ikincil uçucu plakanın terminal hızındaki hesaplama hatası yüzde 2,0 ile yüzde 4,0 arasında kalmıştır.

En iyi tasarımın taban yarıçapı R = 7 mm, kenar profil üssü a = 0,9 ve merkezdeki platin katman kalınlığı d = 0,82 mm olarak belirlenmiştir. Bu yapı, ikincil uçucu plakanın merkez terminal hızını 11,92 km/s’den 13,35 km/s’ye çıkararak yüzde 12 artış sağlamıştır. Merkezden 1 mm yarıçaptaki hedef sinyali ile merkez sinyali arasındaki varış zamanı farkı 2,71 nanosaniyeden 0,35 nanosaniyeye düşmüş; böylece 2 mm çaplı ölçüm bölgesindeki zaman farkı yüzde 88 azalmıştır.

İyileştirme, bütün uçucu plaka noktalarını aynı hızda hareket ettirmekten kaynaklanmamıştır. Platin katman merkez basıncını ve yükleme süresini artırarak merkez bölgeyi kenarlardan biraz daha hızlı hâle getirmiştir. Bu kontrollü radyal hız farkı, başlangıçta içbükey biçimde deforme olan uçucu plakayı hedefe ulaşıncaya kadar daha düz bir geometriye taşıyan telafi mekanizması olarak çalışmıştır. Bununla birlikte aşırı kalın platin katmanı merkezi gereğinden fazla hızlandırmış ve dışbükey deformasyona yol açmıştır. Sonuçlar fiziksel MHGDI deneyi değil, doğrulanmış fakat idealize edilmiş hidrodinamik simülasyon sonuçlarıdır.

Araştırmanın temel mühendislik sorunu nedir?

Üç aşamalı hafif gaz fırlatıcıları, aşırı basınç ve sıcaklık altında malzemelerin davranışını incelemek amacıyla ince uçucu plakaları 10 km/s’nin üzerine çıkarabilmektedir. Ancak yüksek terminal hız tek başına yeterli değildir. Uçucu plakanın merkezi ve kenarları hedefe farklı zamanlarda ulaşırsa çarpma yüzeyi eğrilir, hedefin kullanılabilir ölçüm alanı küçülür ve farklı noktalardan alınan basınç veya parçacık hızı sinyalleri eş zamanlı olmaktan çıkar.

Gelişmiş hiperhızlı fırlatıcıda büyük çaplı bir kademeli yoğunluklu çarpıcı, kendisinden çok daha küçük çaplı ikincil uçucu plakaya çarpmaktadır. Büyük çap farkı, fırlatıcı namlusunun yüksek empedanslı kenarlarıyla etkileşen basınç dalgalarının çevreden merkeze doğru yakınsamasını sağlar. Bu yakınsama merkez bölgesini ek olarak hızlandırarak yüksek terminal hıza katkı verir. Aynı mekanizma, yüklemenin farklı yarıçaplarda farklı zamanlarda gerçekleşmesine ve uçucu plakanın ciddi biçimde eğrilmesine de neden olmaktadır.

Araştırma bu nedenle iki hedefi birlikte ele almıştır:

• İkincil uçucu plakanın terminal hızını azaltmadan, mümkünse artırarak 10 km/s’nin belirgin biçimde üzerine çıkarmak.
• Uçucu plakanın farklı yarıçaplarının hedefe ulaşma zamanları arasındaki farkı azaltarak daha düz bir çarpma yüzeyi elde etmek.

İki aşamalı ve üç aşamalı hafif gaz fırlatıcıları arasındaki fark nedir?

İki aşamalı hafif gaz tabancasında ilk güç kaynağı bir pistonu hızlandırır. Piston, helyum gibi hafif gazı sıkıştırır; yüksek basınçlı gaz da sabot içindeki mermiyi veya uçucu plakayı hedef odasına doğru fırlatır. Çalışmada geleneksel iki aşamalı sistemlerin yaklaşık 8 km/s düzeyine kadar kullanıldığı belirtilmiştir.

Üç aşamalı sistemde ise ikinci pompa tüpünden çıkan ilk uçucu eleman doğrudan nihai hedefe gitmez. Bu eleman, farklı yoğunluk ve dalga empedansına sahip çok sayıda katmandan oluşan bir kademeli yoğunluklu çarpıcıdır. Çarpıcı, üçüncü aşamada bulunan ince ikincil uçucu plakaya yük uygular ve onu daha yüksek hıza çıkarır.

Kademeli yoğunluk yapısı, düşük empedanslı polimerlerden yüksek empedanslı metallere doğru düzenlenerek basınç dalgasının zamana yayılmasını ve enerji aktarımının kontrol edilmesini amaçlamaktadır. Gelişmiş fırlatıcı düzeninde çarpıcının yarıçapı ikincil uçucu plakanın yarıçapından büyük tutulur. Böylece namlu kenarlarından merkeze yönelen yakınsayan basınç dalgaları ek hız kazancı sağlar.

GDI, HGDI ve MHGDI ne anlama gelmektedir?

HGDI kenar profili, namlu kenarında oluşan yansıma ve seyrelme dalgalarının zamanını ve şiddetini değiştirmektedir. MHGDI’de eklenen platin katman ise merkezdeki yansımış basıncı ve yükleme süresini artırmaktadır. Araştırmanın temel tasarım yaklaşımı, kenar geometrisi ile merkezdeki yüksek yoğunluklu katmanı birbirinin etkisini dengeleyecek biçimde birlikte ayarlamaktır.

Geleneksel çarpıcı hangi katmanlardan oluşmaktadır?

Başlangıç yapısındaki geleneksel GDI, yükleme yönü boyunca sırasıyla tungsten, bakır, Ti-6Al-4V titanyum alaşımı, alüminyum, magnezyum ve polimetil metakrilat katmanlarından oluşmaktadır.

Bu katmanların başlangıç yarıçapı 12,5 mm’dir. Çarpıcının önünde 1,2 mm kalınlığında ve 14 mm yarıçapında TPX tampon katmanı bulunmaktadır. İkincil uçucu plaka, 0,5 mm kalınlık ve 5 mm yarıçapa sahip tantaldan oluşmaktadır. İlk aşama çarpıcısının başlangıç hızı 6,9 km/s olarak tanımlanmıştır.

Kenar geometrisi nasıl matematiksel olarak tanımlanmıştır?

MHGDI’nin kenar profili şu kuvvet fonksiyonuyla oluşturulmuştur:

\[ Y(x)=h\left[1-\left(\frac{x}{r}\right)^a\right] \]

Burada Y(x), merkezden x kadar uzaklıktaki kalan profil yüksekliğini; h çarpıcının toplam kalınlığını; r profilin kesildiği en büyük radyal uzaklığı; a ise profil kuvvet üssünü göstermektedir. Uzunluk değişkenleri aynı birimle kullanıldığında parantez içindeki oran boyutsuzdur.

• a ≤ 1: Profil dışbükey veya yaklaşık doğrusal biçimdedir.
• a > 1: Profil daha içbükey hâle gelmektedir.
• R: Profilin altında kalan düz merkez tabanının yarıçapıdır.
• d: Merkez tabanına eklenen platin katmanın kalınlığıdır.

Çalışmanın denklem açıklamasında yoğunluk sembolü ρ da tanımlanmıştır; ancak ρ bu profil denkleminde yer almamaktadır. Bu nedenle söz konusu tanımın başka bir ilişkiden kalmış olması veya denklem açıklamasına yanlışlıkla eklenmiş olması mümkündür. Metin bu durumu açıklamamaktadır.

Sayısal model nasıl kurulmuştur?

Hesaplamalar ANSYS Autodyn 2022 R1 yazılımının açık zaman integrasyonlu doğrusal olmayan hidrodinamik çözücüsüyle gerçekleştirilmiştir. Fırlatıcı düzeneği iki boyutlu eksenel simetrik olarak modellenmiştir. Bu yaklaşım, çevresel yönde bütün yapının kusursuz simetrik olduğunu varsaymaktadır.

Büyük şekil değiştirmeleri ve farklı malzemelerin yüksek gerinim hızlı etkileşimlerini çözmek için Euler tipi çok malzemeli hesaplama alanı kullanılmıştır. Malzemeler arası arayüzler hacim kesri taşınımıyla izlenmiş ve ayrı temas tanımları oluşturulmamıştır. Hesaplama alanının dış sınırlarında akış çıkış koşulu uygulanmıştır.

Şok cephesindeki fiziksel olmayan sayısal salınımları bastırmak amacıyla doğrusal yapay viskozite katsayısı 0,2, karesel yapay viskozite katsayısı 1,0 olarak belirlenmiştir. İkincil uçucu plakanın ve hedefin serbest yüzeylerine Lagrange izleyici parçacıkları yerleştirilmiştir. Bu izleyiciler malzemeyle birlikte hareket ederek belirli malzeme noktalarındaki parçacık hızının zamana bağlı geçmişini kaydetmiştir.

İki boyutlu eksenel model neyi dışarıda bırakmaktadır?

Eksenel simetri, çarpıcı ve uçucu plakanın merkez eksenine göre kusursuz biçimde hizalandığını varsaymaktadır. Gerçek deneylerde oluşabilecek küçük açısal eğiklikler, eksen kaçıklığı, katman kalınlığı farklılıkları, bağlama kusurları, çevresel asimetri ve üç boyutlu dalga düzensizlikleri modele dahil edilmemiştir.

Araştırmacılar bu seçimi, kenar profili ile radyal dalga yayılımının ideal yapısal etkisini rastgele montaj hatalarından ayırmak amacıyla yapmıştır. Dolayısıyla model optimum tasarımın teorik üst sınırını gösterebilir; gerçek deneyde aynı düzlemselliğin otomatik olarak elde edileceğini kanıtlamaz.

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

Çalışmada 0,005, 0,01, 0,02, 0,04 ve 0,08 mm hücre boyutlarının karşılaştırıldığı görülmektedir. İkincil uçucu plakanın terminal hızı yakınsama ölçütü olarak kullanılmıştır. 0,01 mm ile 0,02 mm ağlar arasındaki terminal hız farkının yüzde 3,5’ten küçük olduğu bildirilmiştir.

Hesaplama süresi ile doğruluk arasında denge kurmak amacıyla sonraki analizlerde 0,02 mm sabit hücre boyutu seçilmiştir. Metin beş ayrı ağ boyutunun incelendiğini belirtmesine rağmen aynı cümlede yalnızca 0,08, 0,02 ve 0,005 mm değerlerini “kaba, orta ve ince” olarak sıralamaktadır; 0,01 ve 0,04 mm değerleri ilgili grafikte görülmektedir.

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

Yeni MHGDI tasarımı doğrudan deneysel olarak denenmemiştir. Bunun yerine temel hidrodinamik model, daha önce gerçekleştirilen üç aşamalı gaz tabancası deneylerindeki tantal uçucu plaka hızlarıyla karşılaştırılmıştır.

Hesaplanan değerler üç koşulda da deneysel hızlardan biraz düşük kalmıştır. Araştırmacılar bu farkı malzeme özellikleri ve ölçüm belirsizlikleriyle ilişkilendirmiştir. Uçucu plakanın 1,0 ve 1,4 mikrosaniyedeki yoğunluk dağılımı ile deformasyon biçimi de önceki sayısal çalışmadaki konturlarla nitel olarak karşılaştırılmış ve benzer davranış gösterdiği belirtilmiştir.

Bu doğrulama, çözücünün temel GDI düzenindeki terminal hızı ve genel deformasyon eğilimini makul doğrulukla temsil edebildiğini göstermektedir. Ancak optimum MHGDI geometrisinin performansını doğrudan doğrulayan fiziksel deney bulunmamaktadır.

Uçucu plakanın düzlemselliği nasıl ölçülmüştür?

İkincil uçucu plaka, lityum florür hedefe doğru hareket ettirilmiştir. LiF hedefinin merkezinde ve merkezden 1, 1,5 ile 2,5 mm uzaklıktaki radyal noktalarda parçacık hızının yükselmeye başladığı varış zamanları izlenmiştir.

Belirli bir radyal konum ile hedef merkezi arasındaki zaman farkı şu genel büyüklükle ifade edilmiştir:

\[ \Delta t_{Tx}=t(r_{\mathrm{target}}=x)-t(r_{\mathrm{target}}=0) \]

Çalışmada bu ilişki numaralı bir denklem olarak ayrıca verilmemiştir; sembol tanımının anlamı bu şekildedir. Δt değeri sıfıra yaklaştıkça merkez ve ilgili radyal nokta hedefe daha eş zamanlı ulaşmakta ve çarpma yüzeyi daha düz kabul edilmektedir.

• Pozitif Δt: Radyal nokta merkezden daha geç ulaşmaktadır.
• Negatif Δt: Radyal nokta merkezden daha erken ulaşmaktadır.
• Δt ≈ 0: İlgili hedef bölgesinde yüksek varış eş zamanlılığı bulunmaktadır.

En iyi tasarımın seçiminde merkez terminal hızı V_F0 ile merkezden 1 mm uzaklıktaki Δt_T1 temel performans göstergeleri olarak kullanılmıştır. Merkezden ±1 mm’lik bölge toplam 2 mm çap oluşturduğu için özetteki “2 mm hedef bölgesi” ifadesi bu ölçüme karşılık gelmektedir.

Geleneksel GDI ikincil uçucu plakayı nasıl hızlandırmaktadır?

Simülasyonlar geleneksel GDI yüklemesini iki aşamaya ayırmıştır.

İlk aşama: yaklaşık 0–0,6 mikrosaniye

Çarpıcının düz merkez yüzeyi ikincil uçucu plakaya doğrudan çarpmakta ve yükleme yaklaşık tek boyutlu gerçekleşmektedir. Uçucu plakanın merkezinden kenarına kadar izlenen parçacık hızları bu dönemde birbirine yakın kalmaktadır.

İkinci aşama: 0,6 mikrosaniyeden sonra

GDI’nin dış bölgesi yüksek empedanslı tungsten namlu uzantısıyla etkileşmektedir. Empedans farkı çevrede yüksek basınçlı bir yansıyan şok oluşturmakta ve bu dalga merkez eksenine doğru yakınsamaktadır.

Basınç ile parçacık hızındaki değişim arasındaki temel ilişki çalışma içinde şu biçimde ifade edilmiştir:

\[ dP=\rho\,du \]

Burada dP basınç değişimini, ρ malzeme yoğunluğunu ve du parçacık hızındaki değişimi ifade etmektedir. Bu kısa ilişki tam bir durum denklemi değildir; empedans eşleştirme mekanizmasını nitel olarak açıklamak amacıyla kullanılmıştır.

Yakınsayan dalga önce uçucu plakanın kenar bölgelerini hızlandırmakta, ardından merkeze ulaşıp yaklaşık 0,85 mikrosaniyede basıncı merkezde yoğunlaştırmaktadır. Merkez hızı zamanla kenar hızını yakalamakta ve aşmaktadır. Bu merkez aşırı yüklemesi yüksek terminal hızın temel kaynağıdır.

Fakat çevreden merkeze taşınan basınç gecikmeli olduğu için plakanın bütün yarıçapları aynı hız geçmişini yaşamamaktadır. Sonuçta ince plaka eğilmekte, orta bölüm ile kenarlar hedefe farklı zamanlarda ulaşmaktadır.

Kenar profil üssü a terminal hızı nasıl etkilemiştir?

İlk parametrik çalışmada taban yarıçapı R = 7 mm’de sabit tutulmuş, profil üssü a değiştirilmiştir.

Dışbükey veya yaklaşık doğrusal profillerde, yani a ≤ 1 koşullarında terminal hız değişimi yüzde 1’den küçük kalmıştır. İçbükey profillerde a > 1 olduğunda ise terminal hız hızla azalmıştır. Merkez terminal hızı a arttıkça 11,92 km/s’den 10,79 km/s’ye düşmüştür.

İçbükey kenar, şok dalgasının profillendirilmiş serbest yüzeye daha erken ulaşmasına neden olmuştur. Buradan yansıyan seyrelme dalgaları, GDI ile namlu arasındaki yüksek basıncı erken boşaltmış ve merkeze doğru yakınsayan basınç dalgasını zayıflatmıştır.

Radyal hız farkı azalırken düzlemsellik neden kötüleşmiştir?

İçbükey profil büyüdükçe merkez ile kenar arasındaki terminal hız farkı azalmıştır. Yalnızca son hız dağılımına bakıldığında bunun daha düz bir plaka oluşturması beklenebilir. Buna rağmen merkezden 2 mm uzaklıktaki hedef sinyalinin zaman farkı 2,71 nanosaniyeden 4,80 nanosaniyeye yükselmiş ve düzlemsellik kötüleşmiştir.

Araştırmacılar bu görünüşteki çelişkiyi şok cephesinin geometrisiyle açıklamıştır. Profillendirilmiş kenarın oluşturduğu erken seyrelme dalgası yalnızca son hız değerlerini değil, yüklemenin ne zaman ve hangi yüzey biçimiyle gerçekleştiğini de değiştirmiştir. Uçucu plaka daha küçük son hız farkına sahip olsa bile yükleme sırasında oluşan eğrilik hedefe kadar korunabilmiştir.

Bu sonuç, düzlemselliğin yalnızca merkez ve kenarın son hızlarını eşitleyerek optimize edilemeyeceğini göstermektedir. Bütün ivmelenme geçmişi, basınç dalgasının şekli ve yüklemenin radyal zamanlaması birlikte değerlendirilmelidir.

Taban yarıçapı R neden önemlidir?

Taban yarıçapı, çarpıcının namlu ve ikincil uçucu plaka üzerinde etkin biçimde yük oluşturduğu merkez alanın büyüklüğünü belirlemektedir. R küçüldüğünde profillendirilmiş serbest kenar merkeze yaklaşmakta ve seyrelme dalgaları çekirdek bölgeye daha erken ulaşmaktadır.

Dışbükey a = 0,5 profilinde R’nin küçülmesi namlu kenarındaki en yüksek basıncı 293 GPa’dan 266 GPa’ya, yaklaşık yüzde 9 azaltmıştır. İçbükey a = 5 profilinde aynı küçülme basıncı 293 GPa’dan 102 GPa’ya, yaklaşık yüzde 65 azaltmıştır. Büyük profil üssü ile küçük taban yarıçapının birleşmesi erken basınç boşalmasını belirgin biçimde güçlendirmiştir.

Etkin yükleme alanı nasıl tanımlanmıştır?

Araştırmacılar, azaltılmış kütleli bir çarpıcının geleneksel yapıyla yaklaşık aynı terminal hızı üretip üretemeyeceğini belirlemek için göreli hız sapmasını kullanmıştır:

\[ \delta=\frac{|V_{F0}-V_0|}{V_0}\times100\% \]

Burada V_F0 incelenen yapının merkez terminal hızını, V₀ = 11,92 km/s geleneksel yapının referans hızını ve δ yüzde cinsinden göreli sapmayı göstermektedir. δ ≤ %3 koşulu kabul edilebilir hız korunumu sınırı olarak seçilmiştir.

Her profil üssü için bu sınırı sağlayan en küçük taban yarıçapı R_min belirlenmiştir:

Profil içbükeyleştikçe seyrelme dalgalarının etkilediği alan genişlemiş ve aynı terminal hızın korunması için daha büyük düz taban gerekmiştir. Bu harita, gereksiz dış kütleyi azaltırken merkezde yeterli yükleme alanını korumaya yönelik tasarım aralığı sunmaktadır.

Platin katman terminal hızı nasıl artırmıştır?

MHGDI tasarımında yoğunluğu ve şok empedansı yüksek platin, çarpıcının merkez tabanına yerleştirilmiştir. Platin katmanın kalınlığı d arttıkça merkezde yansıyan basınç yükselmiş ve basınç yüklemesi daha uzun süre korunmuştur.

R = 5 mm ve a = 0,5 koşulunda merkez terminal hızı yaklaşık 11,60 km/s’den 12,73 km/s’ye yükselmiştir. R = 11 mm’lik geniş tabanda ise platin kalınlığı 1,6 mm’ye çıkarıldığında merkez terminal hızı 11,92 km/s’den 14,23 km/s’ye kadar ulaşmıştır.

R = 5 mm, a = 5 ve d = 1,6 mm örneğinde 0,5 mikrosaniyedeki taban yansıma basıncı, platin içermeyen HGDI’de 225 GPa iken platinli yapıda 297 GPa’ya yükselmiştir. Platin katman özellikle merkezde basıncı uzatırken, profillendirilmiş kenardan gelen seyrelme dalgaları çevresel basıncı azaltmıştır.

Merkezin kenarlardan hızlı olması neden yararlı olmuştur?

Geleneksel yapıda uçucu plakanın yükleme geçmişi, hedefe ulaşırken içbükey bir yüzey oluşturmuştur. Merkezin terminal hızının kenarlardan kontrollü biçimde yüksek tutulması, uçuş sırasında merkez bölümün bu gecikmeyi kapatmasını sağlamıştır.

Bu nedenle radyal hız gradyanı bu çalışmada yalnızca bir düzensizlik değildir. Doğru büyüklükte seçildiğinde başlangıç deformasyonuna zıt yönde çalışan bir şekil düzeltme mekanizmasıdır. Merkez ve kenar hızlarını tamamen eşitlemek yerine, hedefteki varış zamanlarını eşitleyecek bir hız farkı oluşturulmuştur.

Platin katman neden sınırsız biçimde kalınlaştırılamaz?

Platin kalınlığı arttıkça hedef varış zamanı farkı önce sıfıra yaklaşmış, optimum noktadan sonra yeniden büyümüştür. Bunun nedeni merkez bölgenin kenarlardan aşırı hızlı hâle gelmesidir.

• Yetersiz platin kalınlığı: Merkez gecikmesi tamamen giderilemez ve plaka içbükey kalır.
• Uygun platin kalınlığı: Kontrollü merkez hızlanması ilk eğriliği telafi eder ve hedef varış zamanları yaklaşır.
• Aşırı platin kalınlığı: Merkez kenarları geçer, telafi aşırıya kaçar ve plaka dışbükey hâle gelir.

Dolayısıyla terminal hızı en yüksek yapan yapı ile en iyi düzlemselliği sağlayan yapı aynı olmak zorunda değildir. Tasarım çok amaçlı bir hız-düzlemsellik dengesi gerektirmektedir.

En iyi MHGDI yapısının parametreleri nelerdir?

Bu yapı hesaplamalı tarama içinden, yüksek merkez terminal hızı ile düşük hedef varış zamanı farkının birlikte sağlandığı çözüm olarak seçilmiştir.

En iyi tasarımın nicel performansı nedir?

Zaman farkındaki düşüş, plakanın bütün yüzeyinin matematiksel olarak kusursuz düz olduğu anlamına gelmemektedir. Sonuç, seçilen LiF hedef noktalarında parçacık hızı sinyallerinin başlangıç zamanlarının birbirine önemli ölçüde yaklaştığını göstermektedir.

Basınç alanı optimum yapıda nasıl değişmiştir?

Platin katman 0,7 mikrosaniye civarında hem ikincil uçucu plakanın merkezindeki iç basıncı hem de namlu sınırlarındaki arayüz basıncını güçlendirmiştir. Daha sonraki basınç boşalma aşamasında kenar profili çevresel yükün azalmasını sağlarken radyal basınç merkeze doğru yakınsamaya devam etmiştir.

Bu iki mekanizma birlikte:

• Mutlak terminal hızı yükseltmiştir.
• Merkez ile çevre arasındaki terminal hız farkını kontrollü biçimde artırmıştır.
• Başlangıçtaki içbükey deformasyonu hedefe ulaşmadan önce azaltmıştır.
• Aşırı merkez basıncının dışbükey deformasyon oluşturmasını kenar profiliyle sınırlandırmıştır.

Şekillerin sunduğu temel bilgiler nelerdir?

• Şekil 1: İki aşamalı ve üç aşamalı hafif gaz tabancalarının düzenek farklarını göstermektedir.
• Şekil 2: Model doğrulama, profil üssü, taban yarıçapı, platin kalınlığı ve optimum tasarım adımlarını özetlemektedir.
• Şekil 3: Geleneksel GDI katmanlarını, MHGDI geometrisini ve LiF hedefteki ölçüm noktalarını göstermektedir.
• Şekil 4: Ağ yakınsamasını, deney-simülasyon hız karşılaştırmasını ve deformasyon konturu doğrulamasını sunmaktadır.
• Şekil 5: Geleneksel GDI’de çevreden merkeze ilerleyen basınç dalgasının hız farkı ve ciddi eğrilik oluşturmasını göstermektedir.
• Şekil 6–8: Kenar profil üssünün terminal hız, zaman farkı ve seyrelme dalgası üzerindeki etkisini açıklamaktadır.
• Şekil 9–12: Taban yarıçapının etkin yükleme alanını ve gereken minimum yarıçapı nasıl belirlediğini göstermektedir.
• Şekil 13–16: Platin kalınlığının terminal hız, radyal hız farkı, hedef zaman farkı ve basınç alanı üzerindeki etkisini sunmaktadır.
• Şekil 17: Farklı R, a ve d birleşimlerinde plakanın içbükey ya da dışbükey biçime dönüşebildiğini göstermektedir.
• Şekil 18: Optimum R = 7 mm, a = 0,9 ve d = 0,82 mm yapısını ve 13,35 km/s ile 0,35 ns sonuçlarını göstermektedir.
• Şekil 19: Geleneksel ve optimum yapıdaki basınç dağılımını ve deformasyon gelişimini doğrudan karşılaştırmaktadır.

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

• Kademeli yoğunluklu çarpıcının bütün yarıçapı terminal hıza aynı ölçüde katkı sağlamamaktadır; kritik bir etkin yükleme alanı bulunmaktadır.
• Kenar profilinin etkisi taban yarıçapı küçüldükçe güçlenmektedir.
• Güçlü içbükey profil ve küçük taban yarıçapı erken seyrelme oluşturarak merkezdeki yakınsayan basıncı ciddi biçimde azaltmaktadır.
• Merkezdeki platin katman yansıyan basıncı ve yükleme süresini artırarak terminal hızı yükseltmektedir.
• Kontrollü bir merkez-kenar hız farkı, başlangıç deformasyonunu telafi ederek hedefte daha eş zamanlı çarpma sağlayabilmektedir.
• Platin katmanın aşırı kalınlaştırılması merkez bölgeyi fazla hızlandırarak düzlemselliği yeniden bozabilmektedir.
• R = 7 mm, a = 0,9 ve d = 0,82 mm hesaplanan tasarımlar içinde en uygun hız-düzlemsellik dengesini sağlamıştır.
• Optimum yapı terminal hızı yüzde 12 artırmış, seçilen hedef bölgesindeki varış zamanı farkını yüzde 88 azaltmıştır.

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

• Optimum MHGDI fiziksel olarak üretilmemiş veya üç aşamalı gaz tabancasında denenmemiştir.
• Gerçek deneyde 13,35 km/s hız ve 0,35 ns zaman farkının kesin olarak elde edileceği gösterilmemiştir.
• Katmanların üretim toleransları, arayüz bağlanması, boşluklar veya eksen kaçıklıkları değerlendirilmemiştir.
• İki boyutlu eksenel model üç boyutlu eğiklik, yalpalama ve çevresel asimetriyi temsil etmemektedir.
• Platin katmanın üretim maliyeti, işlenebilirliği ve yeniden kullanılabilirliği incelenmemiştir.
• Uçucu plakanın kırılma, parçalanma veya kopma olasılığı ayrıntılı bir hasar ölçütüyle doğrulanmamıştır.
• Yalnızca hedef varış zamanı farkının azalması, bütün hedef yüzeyinde kusursuz düzlemsellik bulunduğunu kanıtlamamaktadır.
• Sonuçlar gerçek uzay aracı koruma sisteminin, zırhın veya meteorit çarpmasının performans değerlendirmesi değildir.

Türkiye açısından ne ifade etmektedir?

Çalışmanın Türkiye’ye aktarılabilir temel yönü, çok katmanlı malzemelerde yoğunluk, empedans ve geometrinin yalnızca toplam enerji aktarımını değil, basınç dalgasının uzay ve zaman içindeki dağılımını da belirlediğini göstermesidir. Bu yaklaşım; yüksek hızlı çarpma fiziği, havacılık ve uzay malzemeleri, parçacık çarpma deneyleri, şok dayanımı ve aşırı koşullardaki malzeme karakterizasyonu açısından kullanılabilir.

Araştırma ayrıca mühendislik optimizasyonunda yalnızca en yüksek hız veya en yüksek basınç gibi tek bir ölçütün yeterli olmadığını göstermektedir. Daha yüksek merkez basıncı hızı artırırken plaka biçimini bozabilir; kenar boşalması deformasyonu azaltırken terminal hızı düşürebilir. Başarılı tasarım, bu etkileri birlikte ayarlayan çok amaçlı bir optimizasyon gerektirmektedir.

Türkiye’de uygulama için yerel üretim yöntemleriyle hazırlanmış kademeli yoğunluklu numuneler, yüksek hızlı görüntüleme veya hız interferometrisi, üç boyutlu modelleme, katman arayüzü dayanımı, üretim toleransları ve gerçek fırlatıcı deneyleriyle ayrıca doğrulama yapılması gerekir.

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

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

Parametrik tasarım değişkenleri

Model doğrulama sonuçları

• Beş ağ boyutu karşılaştırılmış ve 0,02 mm hücre boyutu seçilmiştir.
• 0,01 mm ile 0,02 mm ağların terminal hızları arasındaki fark yüzde 3,5’ten düşük kalmıştır.
• Üç deneysel hız koşulunda hesaplama hataları yüzde 2,0, yüzde 3,9 ve yüzde 4,0 olmuştur.
• Hesaplanan terminal hızlar deneysel değerlerden sistematik biçimde biraz düşük çıkmıştır.
• 1,0 ve 1,4 mikrosaniyedeki deformasyon konturları önceki sayısal sonuçlarla nitel benzerlik göstermiştir.

Geleneksel GDI’nin mekanik yanıtı

Kenar profilinin nicel etkileri

• R = 7 mm ve a ≤ 1 koşullarında terminal hız değişimi yüzde 1’den küçük kalmıştır.
• a arttıkça merkez terminal hızı 11,92 km/s’den 10,79 km/s’ye düşmüştür.
• İçbükey profillerde merkezden 2 mm uzaklıktaki zaman farkı 2,71 ns’den 4,80 ns’ye yükselmiştir.
• Son hız farkı azalmasına rağmen şok cephesi geometrisi nedeniyle düzlemsellik kötüleşmiştir.

Taban yarıçapının nicel etkileri

• R = 11 mm’de profil değişikliği merkez hızını yalnızca 11,92’den 11,91 km/s’ye değiştirmiştir.
• R = 5 mm’de merkez hızı 8,83 km/s’ye kadar düşmüştür.
• R = 5 mm’de zaman farkının işareti +2,79 ns’den −6,09 ns’ye dönmüş ve yüzey biçimi içbükeyden dışbükeye geçmiştir.
• a = 0,5 profilinde namlu kenarı tepe basıncı yüzde 9; a = 5 profilinde yüzde 65 azalmıştır.
• Yüzde 3 hız sapması ölçütüne göre gerekli minimum yarıçap 5,0–8,9 mm arasında değişmiştir.

Platin katmanın nicel etkileri

• R = 5 mm ve a = 0,5 yapısında merkez hızı 11,60’tan 12,73 km/s’ye yükselmiştir.
• Bir örnekte Δt_T1 3,64 ns’den 0,74 ns’ye düşmüştür.
• R = 5 mm ve a = 5 yapısında platin katman 0,5 µs’deki taban yansıma basıncını 225 GPa’dan 297 GPa’ya yükseltmiştir.
• R = 11 mm’de 1,6 mm platin kalınlığı merkez terminal hızını 14,23 km/s’ye kadar çıkarmıştır.
• Zaman farkı platin kalınlığıyla önce sıfıra yaklaşmış, ardından aşırı merkez hızlanması nedeniyle yeniden büyümüştür.

Optimum yapının karşılaştırmalı sonucu

Belirsizlikler ve yöntemsel sınırlılıklar

• Optimum yapı için bağımsız fiziksel deney yapılmamıştır.
• Model iki boyutlu ve eksenel simetriktir; gerçek üç boyutlu kusurlar dışarıda bırakılmıştır.
• Ağ farkı ve geçmiş deney karşılaştırması dışında kapsamlı bir belirsizlik yayılımı sunulmamıştır.
• R, a ve d parametrelerinin üretim toleranslarına duyarlılığı hesaplanmamıştır.
• Katman arayüzlerinin bağ dayanımı ve ayrılma olasılığı açıklanmamıştır.
• Malzeme dayanım ve hasar modelleri, parçalanma ölçütleri ve sıcaklığa bağlı özellikler ayrıntılı biçimde raporlanmamıştır.
• Düzlemsellik yalnızca seçilmiş LiF hedef noktalarındaki zaman farkıyla değerlendirilmiştir.
• Optimizasyon yönteminin tarama çözünürlüğü ve bütün olası yapıların aranıp aranmadığı ayrıntılı olarak verilmemiştir.
• Veriler herkese açık bir depoda paylaşılmamış, talep üzerine sağlanacağı belirtilmiştir.

Metin ve şekillerde belirlenen teknik tutarsızlıklar

• Tablo 1’de yoğunluk birimi kg/m³ olarak verilmiştir; 0,83–21,44 değerleri bu birimle uyumlu değildir. Değerler muhtemelen g/cm³ veya 10³ kg/m³ ölçeğindedir, ancak çalışma bunu açıklamamıştır.
• Tablo 1’deki empedans değerleri için yazılan birim ve sayısal ölçek arasında da açıklanmayan bir uyuşmazlık bulunmaktadır.
• Denklem 1’in açıklamasında ρ yoğunluğu tanımlanmış, ancak ρ denklemde kullanılmamıştır.
• Ağ bağımsızlığı metni beş ağ boyutundan söz etmekte ancak aynı cümlede yalnızca üç değer sıralamaktadır; grafik beş değeri göstermektedir.
• Şekil 10 açıklamasındaki son HGDI yapı parametresi görsel üzerindeki a = 5 etiketiyle uyuşmamaktadır.
• Platinli ve platinsiz basınç karşılaştırmasının “Şekil 8”de olduğu yazılmıştır; ilgili karşılaştırma Şekil 14’te gösterilmektedir.
• Şekil 16 açıklamasında 0,3, 0,4, 0,5 ve 0,7 µs yazılırken görsel panellerinde 0,4, 0,5, 0,7 ve 1,0 µs değerleri bulunmaktadır.

Kaynak ve Yöntem Notu

Çalışmanın tam özgün adı: Design of Multilayer Heterogeneous Graded Density Impactor for Planarity-improved in the Enhanced Hypervelocity Launcher

Yazarlar: Xiaolong Xin; Jian Zhang; Chengcheng Guo; Ruizhi Zhang; Zhiguo Li; Qiang Shen; Guoqiang Luo.

Yazar sıralaması: Yukarıdaki liste çalışmanın özgün yazar sırasını korumaktadır.

Eş birinci yazar: Eş birinci yazarlık veya eş katkı beyanı bulunmamaktadır.

Sorumlu yazarlar: Jian Zhang ve Ruizhi Zhang.

Sorumlu yazar e-posta adresleri: zhangjian178@whut.edu.cn; zhangrz027@whut.edu.cn.

Kurumlar:

• State Key Laboratory of Advanced Technology for Materials Synthesis and Processing, Wuhan University of Technology, Wuhan 430070, China.
• School of Materials Science and Engineering, Wuhan University of Technology, Wuhan 430070, China.
• National Key Laboratory of Shock Wave and Detonation Physics, Institute of Fluid Physics, China Academy of Engineering Physics, Mianyang 621900, China.

Kurum-yazar ilişkisi: Xiaolong Xin, Jian Zhang, Chengcheng Guo, Ruizhi Zhang, Qiang Shen ve Guoqiang Luo Wuhan University of Technology bünyesindeki iki kuruma; Zhiguo Li ise China Academy of Engineering Physics bünyesindeki kuruma bağlıdır.

DOI:10.2139/ssrn.6952640

Dergi: Nihai hakemli dergi yayını doğrulanmamıştır.

Yayın platformu: SSRN.

Özgün yayınevi: Nihai bir dergi yayınevi bulunmamaktadır; çalışma SSRN üzerinde preprint olarak yayımlanmıştır.

Yayın yılı: 2026.

Kaynak türü: Deneysel karşılaştırmalarla doğrulanmış hidrodinamik simülasyon ve parametrik mühendislik optimizasyonu içeren preprint araştırma çalışması.

Hakemlik durumu: Çalışma hakem değerlendirmesinden geçmemiştir. Bulgular bu sınırlılık dikkate alınarak değerlendirilmelidir.

Resmî kaynak bağlantısı:https://ssrn.com/abstract=6952640

Yazar katkıları: Xiaolong Xin özgün taslak, doğrulama, biçimsel analiz ve veri düzenleme görevlerini üstlenmiştir. Jian Zhang gözden geçirme ve düzenleme, gözetim, kaynaklar ve yöntem çalışmalarına katkı sağlamıştır. Chengcheng Guo doğrulama ile gözden geçirme ve düzenlemeye katkıda bulunmuştur. Ruizhi Zhang gözden geçirme ve düzenleme, kaynaklar ve yöntem çalışmalarını yürütmüştür. Zhiguo Li, Qiang Shen ve Guoqiang Luo gözden geçirme ve düzenleme görevlerine katkı sağlamıştır.

Finansman: Araştırma, Hubei Province Central Government-Guided Local Science and Technology Development Special Fund tarafından 2025CSA004 numaralı destek ve National Natural Science Foundation of China tarafından 52494933 numaralı destek kapsamında finanse edilmiştir.

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

Veri erişimi: Verilerin talep üzerine paylaşılacağı belirtilmiştir.

Bu Türkçe bilimsel açıklama, yüklenen çalışmanın metni, denklemi, tabloları, grafik ve basınç konturlarının baştan sona incelenmesiyle hazırlanmıştır. Bilimsel içerik yalnızca çalışmada sunulan hesaplamalara ve açıklamalara dayanmaktadır. Çalışmada bulunmayan deneysel sonuçlar veya dış literatür bulguları eklenmemiştir. Dış kaynaklar yalnızca özgün başlık, yazar listesi, DOI ve SSRN kaydının bibliyografik doğrulanması amacıyla kullanılmıştır.

Çalışmanın en önemli sınırlılığı, önerilen optimum MHGDI’nin fiziksel deneyle doğrulanmamış olmasıdır. İki boyutlu eksenel model gerçek montaj hatalarını, üç boyutlu düzensizlikleri, üretim toleranslarını ve katman arayüzü kusurlarını dışarıda bırakmaktadır. Yüzde 12 hız artışı ve yüzde 88 zaman farkı azalması, doğrulanmış temel çözücüyle elde edilen teorik tasarım sonuçlarıdır; gerçek sistem performans garantisi olarak yorumlanmamalıdır.