Magnezyum alüminat spinel (MgAl₂O₄), MgO–Al₂O₃ sisteminde oluşan kübik spinel yapılı bir oksit seramiğidir. Atmosferik plazma püskürtme (Atmospheric Plasma Spraying, APS), toz parçacıklarını sıcak bir plazma jetinde ısıtıp hızlandırarak bir altlık üzerinde kaplama oluşturan termal püskürtme yöntemidir. İncelenen araştırma, APS ile MgAl₂O₄ kaplama üretiminde ark akımının 400, 500 ve 600 A olarak değiştirilmesini üç boyutlu modelleme, uçuş hâlindeki parçacık ölçümleri ve kaplama karakterizasyonuyla değerlendirdi. Bu koşullarda kaplama kalınlığı 135 µm’den 175 µm’ye, Vickers sertliği 769’dan 900’e yükseldi. Sonuçlar incelenen gaz debileri, toz beslemesi ve 100 mm püskürtme mesafesi için geçerlidir; bütün APS sistemleri için optimum akım belirlemez.

Akım artışıyla torç gücü 25,6 kW’den 36,6 kW’ye çıktı; plazmanın sıcak ve hızlı bölgeleri genişledi. X = 0,10 m kesitindeki hesaplanan ortalama parçacık sıcaklığı 2766 K’den 2918 K’ye, ortalama parçacık hızı 123 m/s’den 154 m/s’ye yükseldi. Aynı kesite gelmeden hesaplama alanını terk eden parçacıkların kütle kesri 0,639’dan 0,302’ye düşerken erimiş olarak ulaşanların kütle kesri 0,234’ten 0,472’ye çıktı. Deneyle karşılaştırma için ayrıca seçilen dar bölgede modelin ortalama sıcaklık hatası %11,49–15,57, hız hatası %0,90–5,08 aralığındaydı. Bu dar bölgenin ortalamaları, bütün kesitin ortalamalarıyla karıştırılmamalıdır.

Araştırma problemi: Torç ayarından kaplama yapısına uzanan zincir

MgAl₂O₄; yüksek sıcaklık, elektriksel yalıtım, aşınma ve işlevsel yüzey uygulamaları bakımından ilgi gören bir kaplama malzemesidir. Kaynağın girişinde bu ilgi; yüksek ergime sıcaklığı, düşük yoğunluk, sertlik, kimyasal kararlılık ve ısıl şok direnciyle ilişkilendirilir. Yeterince yoğunlaştırılmış spinelin görünür ve kızılötesi bölgede optik geçirgenlik gösterebilmesi de arka plan bilgisidir. Ancak bu araştırmada optik geçirgenlik, dielektrik dayanım, aşınma ömrü veya uzun süreli yüksek sıcaklık dayanımı ölçülmemiştir. Bir malzemenin olası uygulamaları, burada üretilen kaplamaların doğrulanmış performansları değildir.

Kaynak literatür değerlendirmesi; sol-jel, plazma destekli kimyasal buhar biriktirme ve aerosol biriktirme ile oluşturulan spinel kaplamalarını, ayrıca termal püskürtme uygulamalarını özetler. Sayısal APS çalışmalarında FAST-3D, LAVA, COMSOL, CFX ve FLUENT gibi farklı çözücülerle Al₂O₃, W, Fe, Ni, NiCrAlY, ZrO₂ ve itriya ile stabilize zirkonya (YSZ) parçacıkları incelenmiştir. Kaynağın belirlediği boşluk, doğrudan MgAl₂O₄ parçacıklarını hedefleyen ve uçuş davranışını deneyle karşılaştıran modellemelerin sınırlı olmasıdır. Burada özgün araştırma sorusu, diğer işlem değişkenleri sabitken ark akımının plazma alanını, parçacıkların uçuş durumunu ve oluşan kaplamayı birlikte nasıl değiştirdiğidir.

APS sırasında olaylar kısa zaman ölçeklerinde ve çok yüksek sıcaklıklarda gerçekleşir. Yalnız kaplama sonrasında ölçüm yapmak, parçacığın uçuş boyunca yeterince ısınıp ısınmadığını veya sıcak çekirdeği ne kadar süre izlediğini doğrudan göstermez. Bu nedenle araştırma, süreç içi optik ölçümü ve hesaplamayı birleştirir. Deney, modelin belirli çıktılarının sınanmasını sağlar; model ise ölçüm noktasına gelmeden önceki uzaysal davranışı yorumlamaya yardım eder. Bu tamamlayıcılık, modelin her fiziksel ayrıntısının doğrulandığı anlamına gelmez.

Atmosferik plazma püskürtmede ark akımı parçacıkları nasıl etkiler?

Ark akımı, torca verilen elektriksel güç üzerinden plazma jetinin sıcaklık ve hız alanlarını değiştirir. Bu çalışmada 400–600 A arasında artan akım, hem parçacıkların ısınabileceği hem de ivmelenebileceği bölgeleri genişletmiştir. Parçacıklar daha hızlı hareket ettiğinde sıcak bölgede kalış süresi kısalabilse de incelenen koşullarda güç ve sıcak bölge genişlemesinin etkisi baskın gelmiş, ortalama parçacık sıcaklığı ve hızı birlikte yükselmiştir. Sonuç, yalnız akımdan bağımsız olmayan torç gücü ve akış alanı değişimleriyle birlikte yorumlanmalıdır.

Torç, aktarılmayan doğru akım arkıyla çalışır: ark katot ve anot arasında torç içinde oluşur. Plazma jeti dışarı çıktıktan sonra, jet doğrultusuna dik toz enjeksiyonu yapılır. Parçacıklar plazmadan enerji ve momentum alır; yeterince eriyen parçacıklar altlığa çarpıp yayılır. Çarpma sonrasında oluşan yassı birikim birimleri, kaynakta kullanılan ifadeyle splat’lar, üst üste gelerek kaplama tabakasını oluşturur. Parçacığın erime durumu kadar altlığa ulaşma hızı ve yörüngesi de son yapıyla ilişkilidir.

Neden küçük ve büyük parçacıklar aynı davranmaz?

Kaynak modelde parçacık çapları 12–63 µm aralığındadır. Daha küçük parçacıklar jetin sıcak ve hızlı merkez bölgesini izleyerek daha etkin ısınır ve hızlanır. Daha büyük parçacıklar ise dik enjeksiyon sırasında jetin altına doğru nüfuz edebilir; daha soğuk ve yavaş bölgelere ilerlediklerinde yeterince ısınmadan hesaplama alanının alt sınırını terk edebilir. Bu ayrım, kaynak Şekil 8’de sıcaklıkla, Şekil 10’da hızla renklendirilen yörüngelerde görülür. Aynı toz beslemesi içinde bile tek bir sıcaklık veya tek bir hız bulunmamasının nedeni budur.

Bu açıklama, bütün küçük parçacıkların mutlaka kaplamaya katıldığı veya bütün büyük parçacıkların kaybolduğu biçiminde okunmamalıdır. Kaynak yörüngeleri çap ile akış alanının birlikte etkisini gösterir. Ayrıca yöntem bölümünde 12–63 µm aralığı verilirken bazı sonuç paragrafları en büyük parçacığı 64 µm olarak tanımlar. Bu küçük fakat gerçek tutarsızlık, yeni bir boyut dağılımı varsayılarak giderilmemiştir.

DPV sistemi parçacık sıcaklığını ve hızını nasıl ölçer?

DPV, uçuş hâlindeki parçacıkların yaydığı ışığı kullanarak sıcaklık ve hızlarını birlikte belirleyen bir optik tanılama sistemidir. Sıcaklık iki farklı dalga boyundaki ışınım sinyallerinin oranından, hız ise iki optik ölçüm konumu arasındaki geçiş süresinden hesaplanır. Sıcaklık ölçümünde iki dalga boyundaki yayma katsayılarının yaklaşık eşitliği varsayılır; bu nedenle sinyal oranı yöntemi bütün ölçüm belirsizliklerini ortadan kaldırmaz.

Kaynakta ölçüm hacmi yaklaşık 1 × 10⁻⁹ m³, yani 1 mm³; en yüksek veri toplama hızı 800 parçacık/s olarak verilir. Bunlar cihazın özellikleridir; her deneyde 800 parçacık/s hızına ulaşıldığı bildirilmemiştir. Algılama kafasındaki iki yarıklı maske ve optik fiber demeti üzerinden taşınan ışık, dikroik ayna ile ayrılır. Merkez dalga boyları λ₁ = 785 ± 25 nm ve λ₂ = 995 ± 25 nm olan bant geçiren filtrelerden sonra iki silisyum fotodedektör sinyalleri algılar. Kaynak, DPV ailesini tanıtır ve DPV-EVOLUTION kılavuzuna atıf verir; kullanılan cihazın tam modelini yöntem metninde açıkça kesinleştirmez.

Işınım oranından sıcaklığa: Denklem (1)–(4)

Denklem (1)

E(λᵢ), ilgili dalga boyunda algılanan ışınım büyüklüğüdür; ε(λᵢ) boyutsuz yayma katsayısı, dᵢ parçacık çapı, T yüzey sıcaklığıdır. K, ölçüm sistemi ve ışınım bağıntısına bağlı orantı katsayısı; K₂ = hc/k ise ikinci ışınım sabitidir. Buradaki h Planck sabiti, c ışık hızı, k Boltzmann sabitidir. Dalga boyu ve çap metre, sıcaklık kelvin olarak tutarlı kullanılmalıdır; K₂’nin birimi m·K olduğundan üstel ifadenin içi boyutsuzdur. Kaynak E ve K için ayrı bir cihaz sinyal birimi vermez. Bu bağıntı, sıcaklığın ışınımı etkilediğini, ancak mutlak sinyalin çap ve yayma katsayısından da etkilendiğini gösterir.

Denklem (2)

R, iki kanaldaki ışınımın boyutsuz oranıdır. Aynı parçacığın iki dalga boyundaki yayma katsayıları yaklaşık eşit kabul edilince ortak çarpanlar oran içinde sadeleşir. Denklem (1) paydasındaki “−1” terimi Denklem (2)’de Wien yaklaşımıyla ihmal edilmiş biçime geçer; bu nedenle (2), bütün koşullarda (1)’in yaklaşık kabulsüz tam oranı değildir. Kaynağın belirttiği gri cisim ve Wien kabulleri, ölçümden sıcaklık çıkarılmasının koşullarıdır.

Denklem (3)

Bu ifade, Denklem (2)’nin sıcaklık için çözümüdür. Logaritmalar boyutsuz oranlara uygulanır; ön çarpanın birimi kelvindir. Tₚ parçacığın optik yöntemle belirlenen yüzey sıcaklığıdır. Sinyal oranındaki hata sıcaklık hesabına taşınır; aynı mutlak ışık şiddeti farklı çaplardan gelebileceği için sıcaklığın tek kanal şiddetiyle özdeşleştirilmemesi gerekir.

Denklem (4)

A₁ ve A₂, iki dedektör kanalındaki sinyal şiddetleridir. Kaynağın ölçüm tanımında sıcaklığı belirleyen giriş, bu sinyallerin oranıdır. Aₚ ile gösterilen parçacık yüzey alanı, ilerideki ısı denklemine aittir; A₁ ve A₂ alan değildir. Kanal kalibrasyon katsayıları veya ham sinyaller kaynakta verilmediği için burada sayısal bir sıcaklık hesabı yeniden yapılmamıştır.

Denklem (5)

Vₚ parçacık hızı, d_slit iki yarık arasındaki mesafe, TOF_f iki ölçüm konumu arasındaki uçuş süresidir. Mesafe metre ve süre saniye ise sonuç m/s olur. Aynı optik mesafede daha kısa geçiş süresi daha yüksek hız anlamına gelir. Kaynak, cihaz kılavuzunda ölçüm hacmi içindeki yarık görüntüleri arasındaki uzaklığa değinir; fiziksel maske aralığının ölçüm uzayına dönüşümünü ve sayısal aralığı raporlamaz. Bu nedenle formül, kalibrasyonu bilinmeyen bir maskeden bağımsız bir hız hesabı yapmak için yeterli değildir.

İki bölge, iki çözücü ve iki yönlü parçacık–plazma etkileşimi

Torç bölgesi, OpenFOAM v2112 içinde rhoPimpleFoam tabanlı üç boyutlu çözücüyle hesaplanır. Yerel termodinamik denge (Local Thermodynamic Equilibrium, LTE) kabulü altında elektromanyetik alan ve akış birlikte ele alınır; Lorentz kuvveti momentum denklemine, Joule ısınması enerji denklemine katılır. Torç çıkışındaki sıcaklık ve hız dağılımları, ANSYS FLUENT 2021 R1 ile çözülen kaplama bölgesinin giriş koşullarıdır. Bu bölgesel aktarım, iki yazılımın tek bir ortak zaman çözümünde bütünüyle eşzamanlı çalıştığının kanıtı değildir; kaynak sınır koşulu aktarımı tanımlar.

Kaplama bölgesindeki plazma, kararlı ve sıkıştırılabilir sürekli ortam olarak modellenir. Ayrık faz modeli (Discrete Phase Model, DPM), parçacıkları Lagrange yaklaşımıyla yörüngeleri boyunca izleyen bir modelleme yöntemidir. İki yönlü etkileşim, plazmanın parçacıkları ısıtıp hızlandırmasının yanında parçacıkların da plazmanın enerji ve momentum alanını değiştirmesidir. UDF olarak anılan kullanıcı tanımlı fonksiyonlar uygulanmıştır; ayrıntılı kaynak kod verilmez. Ar–H₂ karışımının termodinamik ve taşınım özellikleri sıcaklık ve gaz bileşimine bağlıdır. Torç dışındaki elektrik alan etkisi ihmal edilir; bu bölgede parçacık davranışının ısıl ve aerodinamik etkileşimlerle belirlenmesi esas alınır.

Plazma alanının korunma denklemleri

Denklem (6)

Kararlı akışın kütle korunumu denklemidir. ρ plazma yoğunluğu (kg/m³), v hız vektörüdür (m/s). Akı ρv’nin birimi kg/(m²·s), diverjansının birimi kg/(m³·s) olur. Sıfır sağ taraf, yazılan kararlı denklemde yerel net kütle birikimi bulunmadığını ifade eder. Sıkıştırılabilir akışta yoğunluk değişebileceğinden bu ifade doğrudan ∇·v = 0 ile aynı değildir.

Denklem (7)

Momentum taşınımını basınç, viskoz gerilmeler ve yerçekimiyle ilişkilendirir. p basınç (Pa), τ gerilme tensörü (Pa), g yerçekimi ivmesidir (m/s²). Her terim hacim başına kuvvet boyutundadır. Hızların yan yana yazılması tensörel momentum akısını temsil eder. Kaynak bu denklemi kaplama bölgesi bağlamında verir; torç içindeki Lorentz terimi ayrıca metinde anlatılır. İki yönlü etkileşimden söz edilmesine karşın parçacık kaynak terimleri basılı denklemde açıkça yer almaz; burada denkleme raporlanmamış bir terim eklenmemiştir.

Denklem (8)

Enerji denkleminde E özgül toplam enerji olarak okunur (J/kg); böylece ρE ile basınç aynı boyuttadır. Kaynak E’yi kısaca toplam enerji diye adlandırır. k_eff etkin ısıl iletkenlik (W/(m·K)), τ_eff etkin gerilme tensörü (Pa), Q_rad hacimsel ışınım kaybıdır (W/m³). Negatif ışınım terimi enerji kaybını gösterir. Isıl iletim ve gerilme işi akılarının diverjansı da W/m³ boyutundadır. Türbülans standart k–ε modeliyle hesaplanır. Bu modeldeki k türbülans kinetik enerjisidir; optik bağıntıdaki Boltzmann sabiti k ile aynı büyüklük değildir.

Parçacığın hızlanması ve kaynak Reynolds sayısı sorunu

Denklem (9)

mₚ parçacık kütlesi (kg), vₚ parçacık hızı (m/s), ρₚ parçacık yoğunluğu (kg/m³), τᵣ gevşeme süresidir (s). İlk sağ taraf terimi parçacık ile plazma arasındaki hız farkının sürükleme yoluyla azalmasını, ikinci terim yoğunluk farkıyla düzeltilmiş yerçekimi etkisini temsil eder. Bütün terimlerin birimi newtondur. Parçacık plazmadan yavaşsa sürükleme onu hızlandırabilir; plazma akışı aşağı akımda yeterince yavaşladığında aynı etkileşim parçacığın yavaşlamasına katkı verebilir.

Denklem (10)

dₚ parçacık çapı (m), μ plazmanın dinamik viskozitesi (Pa·s), C_d sürükleme katsayısı ve Reₚ parçacık Reynolds sayısıdır. C_d ve Reₚ boyutsuz olduğunda ifade saniye verir. Diğer çarpanlar sabit tutulduğunda çapın karesi arttıkça gevşeme süresi büyür; parçacığın akış hızına uyumu daha yavaş olur. Ancak Reynolds sayısı ve sürükleme katsayısı da koşullarla değiştiğinden bütün parçacıklar için yalnız dₚ² ile belirlenen tek bir hızlanma yasası çıkarılamaz.

Denklem (11)

Kaynakta basıldığı biçimiyle aktarılmıştır. PDF’nin 13. sayfasında çap dₚ² olarak yazılmıştır. Verilen SI birimleriyle sağ tarafın boyutu metre olur; dolayısıyla Reynolds sayısı için beklenen boyutsuzluk sağlanmaz. Bu durum metin çıkarma hatası değildir; sayfa görüntüsünde de kare vardır. Boyutsal tutarlılık için çapın birinci kuvvetiyle yazılan ρdₚ|vₚ−v|/μ ifadesi boyutsuz olur; bu açıklama Verianla’nın boyut kontrolüdür, yazarların doğrulanmış düzeltmesi veya kullandıkları kodun tespiti değildir. Sayısal sonuçlar yeniden hesaplanmamış ve basılı formül sessizce değiştirilmemiştir.

Denklem (12)

Kaynak, sürükleme katsayısını beş Reynolds aralığı için parça parça tanımlar. Buradaki Re, parçacık Reynolds sayısına karşılık gelir. C_d boyutsuzdur; aynı bağıntının farklı akış rejimlerine gelişigüzel uygulanmaması için aralık sınırları korunmuştur. Son satırda karekök √Re bulunur. İlk rejimde C_dRe = 24 olduğundan Denklem (10)’daki son çarpan bire eşittir; bu, denklemler arasındaki öğretici bir cebirsel bağlantıdır. Denklem (11)’deki basılı tutarsızlık nedeniyle bu bağıntılar yeniden uygulamaya alınmadan önce yazar kodu veya düzeltilmiş yöntem bilgisiyle teyit edilmelidir.

Parçacık ısınması ve faz geçişlerinin enerji hesabı

Denklem (13)

Bu denklem, taşınımla alınan ısının parçacık sıcaklığına etkisini hesaplar. cₚ özgül ısı (J/(kg·K)), Aₚ parçacık yüzey alanı (m²), h taşınım ısı aktarım katsayısıdır (W/(m²·K)). Buradaki h Planck sabiti değildir. Her iki taraf watt birimindedir. T > Tₚ olduğunda sağ taraf pozitiftir ve parçacık ısınır; T < Tₚ olduğunda soğuma yönünde etki oluşur. Model yığılmış ısıl kapasite kabulü kullanır: tek parçacığın içi her an tek sıcaklıkla temsil edilir. Bu nedenle hesaplanan Tₚ ile optik olarak ölçülen yüzey sıcaklığının eşdeğerliği bir model kabulüne bağlıdır; parçacık içi sıcaklık gradyanları çözülmez.

Erime ve buharlaşmaya karşılık gelen enerji gereksinimleri, sıcaklığa bağlı özgül ısı eğrisine dar aralıklarda dağıtılır. Böylece eklenen enerji, yalnız sıcaklığı yükseltmek için değil faz değişimi enerjisini karşılamak için de kullanılır. Kaynak Şekil 3(b)’deki yüksek ve dar özgül ısı bölgeleri bu sayısal temsilin parçasıdır; bağımsız olarak ölçülmüş yeni bir faz geçiş eğrisi değildir. Özgül ısı fonksiyonunun bütün sayısal noktaları, taşınım katsayısının korelasyonu ve tam UDF kodu verilmediğinden eğri yeniden uydurulmamıştır.

Tablodaki molar entalpiler, kaynak birimleriyle korunmuştur. Bunların J/(kg·K) cinsinden etkin özgül ısıya dönüştürülmesinde kullanılan bileşen ve normalizasyon ayrıntıları verilmediği için yeni cₚ değerleri hesaplanmamıştır. Kaynak, 3250 K civarındaki aşamayı bir yerde MgAl₂O₄’ün Al₂O₃’e ayrışması, başka yerde Al₂O₃ buharlaşmasıyla ilişkili anlatır. Bu dil, ayrıntılı kimyasal kinetik veya çok bileşenli buhar bileşiminin doğrulanmış hesabı olarak genişletilmemelidir.

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

Deney koşulları ve sabit tutulan değişkenler

Bağımsız işlem ayarı ark akımıdır; gerilim tamamen sabit değildir ve üç koşulda 64, 62 ve 61 V olarak raporlanır. Güç değerleri akım ile gerilimin çarpımıyla tutarlıdır. Diğer koşulların sabit tutulması, bu düzenekte akımın etkisini karşılaştırmayı kolaylaştırır. Yine de yalnız üç akım düzeyi incelendiğinden 600 A üzerinde performansın artacağı, güç tüketimi başına optimum verimin bulunduğu veya daha düşük akımların bütün uygulamalarda yetersiz olduğu söylenemez.

Hesaplama geometrisi ve örnekleme

Şekil 2, yan ve aşağı sınırlarda çıkış koşulları bulunan bir uçuş alanını gösterir. Dolayısıyla X = 0,10 m’ye ulaşma, 100 mm altlık konumuna erişimin sayısal göstergesidir; duvara çarpma, yayılma, yapışma veya geri sıçramanın ayrı birikme modeliyle çözüldüğü anlamına gelmez. X = 0,12 m, hesaplama alanının daha aşağı akımdaki değerlendirme kesitidir; 120 mm’de ayrı kaplama deneyleri yapıldığı bildirilmemiştir.

Geometri açıklığı: Bölüm 2.2.2, 4,15 mm’yi iki giriş ekseni arasındaki düşey uzaklık olarak anlatır; sonuç bölümünde enjeksiyon merkezi Z = 10 mm ve jet merkezi Z = 0 olarak verilir. Şekil 2’deki 4,15 mm ölçülendirmesinin referansı, bu eksen ifadesini doğrudan açıklığa kavuşturmaz. Bu nedenle tüm ölçülerden kesin bir üretim çizimi veya tam yeniden çalıştırılabilir geometri türetilmemiştir. Kaynağın yöntem ayrıntıları aktarılırken çelişen ölçüler tek bir varsayımsal değere dönüştürülmemiştir.

Şekil 3(a), çap arttıkça kümülatif kütle kesrinin yükseldiği ölçüm ve model noktalarını, Rosin–Rammler uyum eğrisini ve besleme tozunun SEM görüntüsünü birlikte sunar. Yatay eksen çap (µm), düşey eksen boyutsuz kümülatif kütle kesridir. Görüntü ve eğri, çok boyutlu toz beslemesinin kullanıldığını gösterir; sayısal dağılım parametreleri verilmediğinden grafikten kesin dağılım sabiti türetilmemiştir.

Plazma jetinin sıcaklık ve hız alanları

Bunlar plazmaya ait maksimumlardır; MgAl₂O₄ parçacık sıcaklığı veya hızı değildir. Kaynağın karşılaştırmasında 400 A’ya göre maksimum hız 500 A’da %20, 600 A’da %40 artar; maksimum sıcaklık artışları yaklaşık %2,1 ve %4,2 olarak verilir. Akımın incelenen aralıktaki göreli etkisi, torç çıkışı hızında sıcaklıktan daha büyüktür.

Şekil 4, Y = 0 düzlemindeki sıcaklık alanını; Şekil 6 aynı düzlemdeki hız alanını gösterir. Her iki şekilde a–c panelleri tozsuz, d–f panelleri tozlu akışı, g–i panelleri “tozsuz − tozlu” farkını karşılaştırır. X aşağı akım yönü, Z düşey yöndür; koordinatlar metre cinsindedir. Toz verilmesiyle sıcaklık ve hız alanları asimetrikleşir. Bu farklılık, parçacıkların akıştan enerji ve momentum aldığını gösterir; fark haritası doğrudan parçacık sıcaklığı veya parçacık hızı haritası değildir.

Sıcaklık haritalarında 1000, 2500, 5000, 7500, 10 000 ve 12 500 K izotermleri akım arttıkça genişler. Tozlu ve tozsuz akış arasındaki Δ1000 K sınırı için ısınma etkisi yaklaşık X = 0,004–0,100 m aralığındadır ve akımlar arasındaki değişim küçüktür. Δ2000 K sınırında başlangıç yaklaşık 0,004 m kalırken bitiş 400, 500 ve 600 A’da sırasıyla 0,031, 0,033 ve 0,035 m’ye ilerler. Bir etkileşim bölgesinin uzunluğu, hangi fark eşiğiyle tanımlandığına bağlıdır.

Şekil 5’te merkez çizgisi sıcaklığının tozlu akışta daha hızlı düştüğü görülür. Kaynağın yaklaşık toplu anlatımında tozsuz akış, X = 0,004–0,038 m arasında yaklaşık 14 000 K’den 5160 K’ye; sonra 0,120 m’ye kadar 2300 K’ye iner. Tozlu akışta benzer ilk düşüş 0,004–0,025 m arasında yaklaşık 5140 K’ye ulaşır ve alan sonunda yaklaşık 1500 K görülür. Bu yuvarlatılmış açıklama, üç koşulun bütün eğri noktalarının aynı olduğu şeklinde kullanılmamalıdır.

Sıcaklık farkının maksimum büyüklüğü azalırken konumu aşağı akıma kayar. Bu yüzden “en büyük sıcaklık farkı küçüldü, parçacık ısınması mutlaka kötüleşti” sonucu çıkarılamaz: alanın uzaysal genişliği, parçacık yörüngesi ve plazma düzeyi birlikte önemlidir. Kaynak, başlıca ısınma bölgesini 400 A’da X = 0,004–0,021 m, 500 A’da 0,004–0,023 m, 600 A’da 0,004–0,025 m olarak yorumlar. Başlıca ivmelenme bölgeleri ise aynı sırayla 0,004–0,020; 0,004–0,021 ve 0,004–0,022 m’dir.

Şekil 6’daki 50, 250, 500, 750, 1000 ve 1250 m/s hız sınırları akımla genişler; 600 A’da ayrıca 1500 m/s sınırı görülür. Δ25 m/s fark sınırıyla tanımlanan etkileşim yaklaşık X = 0,005 m’de başlar; 400 A’da 0,115 m, 500 A’da 0,118 m ve 600 A’da 0,120 m’nin ötesine ulaşır. Δ100 m/s üzerindeki sınırlar da genişler. Bu eşik temelli alan, merkez çizgisi farkının maksimumuna göre tanımlanan başlıca ivmelenme bölgesiyle aynı ölçüt değildir.

Toz verildiğinde jetin momentumunun bir bölümü parçacıkları ivmelendirmeye gider ve merkez çizgisi daha erken yavaşlar. Şekil 7’nin yatay ekseni X konumu, düşey ekseni hız veya hız farkıdır; eğriler zaman serisi değildir. Isı ve momentum aktarımının etkili olduğu bölgeler aynı uzunlukta olmadığından parçacık sıcaklığı ile hızının aşağı akımdaki eğilimlerinin farklı olması mümkündür.

Parçacık sıcaklıkları, hızları ve uzaysal dağılımlar

Şekil 8 ve 10’da, X–Z yörüngelerine ek olarak X = 0,10 m ve 0,12 m’deki Z–Y dağılımları gösterilir. Parçacıkların çoğu jet merkezinin altında, Z < 0 bölgesindedir. Küçük ve daha sıcak/hızlı parçacıklar üst kısma, daha büyük ve daha soğuk/yavaş parçacıklar alt kısma yönelir. Akım arttığında parçacık taşınımı jet merkezine doğru daha fazla yoğunlaşır. Bunlar sayısal dağılımlardır; kaynakta her parçacığın ham çap ve yörünge dosyası bulunmaz.

Şekil 9(b) ve 11(b)’nin kutuları %25–75 aralığını, bıyıkları minimum ve maksimumu, küçük beyaz işareti ortalamayı, yatay çizgiyi medyanı gösterir. Bu bıyıklar güven aralığı değildir. Çeyrek değerler ve medyanların kesin sayıları verilmediği için tablolara tahminle eklenmedi. Kesit ortalamaları incelendiğinde sıcaklık 100 mm’den 120 mm’ye bütün koşullarda artmaz: 400 A’da 2766 K’den 2765 K’ye çok küçük bir düşüş vardır. Hız ise üç koşulda da azalır.

Kaçan parçacıklar ortalama sıcaklığın yorumunu neden değiştirir?

Hesaplama alanını terk eden soğuk parçacıklar sonraki kesitlerin ortalamasına katılmaz; geride kalan parçacıkların ortalaması bu nedenle daha yüksek görünebilir. Dolayısıyla aşağı akımda artan ortalama sıcaklık, enjekte edilen bütün tozun sürekli ısındığını tek başına göstermez. Bu çalışmada sıcaklık ve hız ortalamaları, kesitte kalan parçacık nüfusu ve kaçan kütle kesriyle birlikte okunmalıdır.

Kaynak, Z < −0,045 m’den aşağı kaçışı özellikle vurgular. İlk kaçış X = 0,041 m’de 400 A, 0,048 m’de 500 A ve 0,058 m’de 600 A için ortaya çıkar. Kaçıştan önceki bölgede raporlanan en yüksek ortalama sıcaklıklar sırasıyla yaklaşık X = 0,040 m’de 2625 K; 0,047 m’de 2727 K; 0,0474 m’de 2865 K’dir. Bu karşılaştırma, bütün kesitlerde kalan parçacıkların giderek seçilmesinden kaynaklanan yanlılığa karşı kaynağın yaptığı kontrolü temsil eder.

Ortalama hız, enjeksiyonun hemen ardından hızla yükselir ve kaynak anlatımında yaklaşık X = 0,025 m sonrasında genel olarak azalma eğilimine geçer. Kaçış etkisi dışlanan aralıktaki maksimum ortalama hızlar 400 A’da X ≈ 0,040 m’de 123 m/s, 500 A’da X ≈ 0,043 m’de 140 m/s ve 600 A’da X ≈ 0,043 m’de 161 m/s olarak verilir. Hız bölümünde 400 A ilk kaçışı 0,040 m diye yuvarlanırken sıcaklık bölümünde 0,041 m yazılır. Grafiklerin yaklaşık anlatımı ile tekil maksimum konumları arasında daha yüksek hassasiyet varsayılmamıştır.

Kesite ulaşan kütlenin erime sınıfları

Bu sınıflandırma sıcaklık eşiklerine dayanır; her parçacığın eriyen hacim yüzdesinin doğrudan ölçümü değildir. Sayısal enerji modelinde erime 2381–2390 K aralığına yayılmış olmasına rağmen sınıflandırmada 2390 K’ye kadarki grup “erimemiş” adını alır. Ayrıca eşikler kaynakta tam kelvin aralıklarıyla yazılmıştır; sürekli sıcaklıklar için aradaki sınır değerlerinin yazılımda nasıl sınıflandırıldığı açıklanmaz. Kaynağın sınıfları yeni sınırlar eklenmeden aktarılmıştır.

Paylar yalnız kesite varan parçacıklar içinde normalize edilmiş erime yüzdeleri değildir; kaynak kaçan kütleyi aynı değerlendirmeye dahil eder. Akım artarken erimiş payın yanında erimemiş payın da artması bu nedenle çelişki değildir: daha önce kaybolabilecek parçacıkların daha fazlası değerlendirme kesitine ulaşır. “Erimemiş pay arttı, ısınma kötüleşti” yorumu toplam taşınım değişimini gözden kaçırır. Aynı biçimde erimiş kütle kesri, doğrudan tartılarak ölçülmüş kaplama birikme verimi değildir.

Kütle bilançosu notu: Tablodaki dört sınıfın toplamı X = 0,10 m’de sırasıyla 0,998; 0,999; 1,000, X = 0,12 m’de 0,993; 0,995; 0,994 eder. Bunlar Verianla’nın verilen sayılara uyguladığı aritmetik kontroldür, yeni araştırma sonucu değildir. Özellikle 0,12 m toplamlarının eksikliği yalnız standart üç ondalıklı yuvarlamayla açıklanacak büyüklükte değildir; kaynak sebebini belirtmez. Eksik pay başka bir sınıfa eklenmemiş ve veriler toplamı bire eşitleyecek biçimde yeniden normalize edilmemiştir.

Deneyle karşılaştırmada ölçüm alanı neden önemlidir?

DPV ile, optik sinyalin en yoğun parçacık geçişini gösterdiği 1 mm² alanda 2000 parçacık ölçülür. Modelde 1300 yörüngeden 1 mm²’den küçük yoğun bölgeye yalnız 14–18 parçacık düştüğü için yazarlar sayısal karşılaştırmayı 5 mm² alana genişletir. Bu alana 400, 500 ve 600 A’da sırasıyla 142, 116 ve 146 parçacık girer. Kaynak son sayının ardından yanlışlıkla A birimini kullanır; burada akım değil parçacık sayısı olduğu bağlamıyla belirtilmiştir. Alan genişletme örnek sayısını artırır fakat deneyle model arasında tam mekânsal eşleşmeyi bozar.

Karşılaştırma Bölüm 3.3’te X = 0,10 m için anlatılır. Şekil 13 açıklamasında hem 0,10 hem 0,12 m yazmasına rağmen 0,12 m için ayrı deneysel doğrulama serisi verilmez. Bu nedenle aşağıdaki tablolar yalnız yöntem metninde tanımlanan X = 0,10 m karşılaştırmasına aittir. Deneyin 1 mm² analiz alanı, cihazın 1 mm³ ölçüm hacmi ve simülasyonun 5 mm² örnekleme alanı farklı niceliklerdir.

Bağıl hata değerleri, ortalamalar arasındaki mutlak farkın deney ortalamasına bölünmesinin yüzde karşılığıyla tutarlıdır. Model sıcaklığı her koşulda yüksek tahmin eder ve bağıl hata akımla artar. Hızda ise ortalama uyumu daha iyidir; 400 A’da model daha düşük, 500 ve 600 A’da daha yüksek değer verir. Fakat ortalamaların yakın olması dağılımların bütününün aynı olduğu anlamına gelmez. Örneğin 600 A’da deneysel en yüksek hız 379 m/s, modelde 271 m/s’dir. Kaynak kutu grafiklerindeki dağılım farkları da bu nedenle önemlidir.

Yazarlar sıcaklık farkını en azından kısmen iki renkli pirometrinin belirsizliğiyle ilişkilendirir. Bu açıklama bütün farkın yalnız cihazdan kaynaklandığını kanıtlamaz; farklı örnekleme alanı ve model kabulleri de yorum sınırıdır. Hız için kaynak, uçuş süresi yönteminin genel belirsizliğinin %2’den küçük olduğuna değinir; bu ifade burada yapılan her ölçüm için çıkarılmış özel bir belirsizlik bütçesi değildir. Ayrıca 400 ve 500 A’daki model–deney ortalama hız farkları bu %2 düzeyinden yüksektir.

Dikkat edilmesi gereken başka bir ayrım, deneydeki bazı maksimum sıcaklıkların modelin yaklaşık 3774 K üst düzeyini aşmasıdır. 400 A’da 4302 K, 600 A’da 4386 K maksimumları raporlanmıştır. Bu tekil optik değerler, modelin “buharlaşmış parçacık öngörülmedi” sonucunun deneysel olarak tam doğrulandığı biçiminde kullanılamaz. Model sınıflandırması ile ölçüm sinyalinden çıkarılan sıcaklık ayrı kanıt türleridir.

Kaplama mikroyapısı, fazı, kalınlığı ve sertliği

Şekil 14, üç akımda üretilen kaplamaların kesit taramalı elektron mikroskobu (SEM) görüntüleridir. Görüntülerde altlık, MgAl₂O₄ tabakası ve numune hazırlama reçinesi ayırt edilir; gözenekler ve kısmen erimiş parçacıklar işaretlenmiştir. Ölçek çubukları 50 µm’dir. Kaynak görüntüleri, kaplamanın akım arttıkça kalınlaşmasını ve mikroyapı farklarını gösterir; tek başına görüntüye bakılarak yeni bir gözeneklilik yüzdesi veya tane boyutu çıkarılmamıştır.

Kalınlık ve gözeneklilik, her akım koşulu için alınan on kesit SEM görüntüsü analiz edilerek nicelleştirilmiştir. On görüntü, on bağımsız kaplama üretim tekrarı anlamına gelmez; bağımsız üretim sayısı belirtilmez. Şekil 15’in sertlik ekseni HV₀,₁ olarak etiketlidir. Ancak baskı sayısı, yük uygulama süresi ve cihaz modeli gibi sertlik yöntem ayrıntıları verilmez. Grafikte hata çubukları vardır fakat bunların sayısal değerleri ve standart sapma, standart hata veya başka bir yayılım ölçütü olup olmadığı açıkça tanımlanmaz. Bu nedenle tabloya hata çubuğu veya güven aralığı uydurulmamıştır.

Ark akımını artırmak kaplama gözenekliliğini her adımda azaltır mı?

Hayır. Bu çalışmada gözeneklilik 400 A’da %5,09 iken 500 A’da %5,50’ye yükselmiş, 600 A’da %4,02’ye düşmüştür. En düşük raporlanan değer 600 A’ya aittir; ancak 400–600 A boyunca tekdüze bir azalma yoktur. Buna karşılık kalınlık ve Vickers sertliği üç koşul boyunca artmıştır. İstatistiksel anlamlılık testi verilmediği için küçük farkların anlamlılığı ayrıca ileri sürülemez.

Yazarlar daha yüksek parçacık hızı ve erimiş ulaşım payını, altlığa daha yüksek kinetik enerjiyle çarpma ve splat’lar arası bağın iyileşmesiyle ilişkilendirir. Bu mekanizma yorumu, simülasyon ile kaplama gözlemlerinin birlikte değerlendirilmesine dayanır. Çarpma enerjisinin bağımsız ölçümü veya doğrudan bağ dayanımı deneyi sunulmaz. Bu yüzden artan sertlik, ölçülmemiş yapışma dayanımının nicel kanıtı olarak kullanılamaz.

X-ışını kırınımı (XRD), kristal fazları kırınım desenleri üzerinden inceleyen yöntemdir. Şekil 15(a)’da yatay eksen 2θ açısı (derece), düşey eksen göreli şiddettir; başlangıç tozu ile üç kaplamanın desenleri karşılaştırılır. Kaynak, belirgin yeni faz veya desen değişimi gözlenmemesini MgAl₂O₄ fazının korunduğu biçiminde yorumlar. Bu gözlem, ölçümün duyarlılığı içindeki faz kararlılığını destekler; iz düzeyinde bileşim değişimini veya tam stokiyometrik eşitliği bağımsız kimyasal analiz gibi kanıtlamaz. Kaynağın sonuç bölümündeki “tespit edilebilir değişim olmadan” ifadesi bu sınırı daha doğru yansıtır.

Şekil etiketleriyle ilgili not: Şekil 14 ve 15’in açıklamalarında 300, 400 ve 500 A yazılır; yöntem, Tablo 1, sonuç metni ve görüntü/grafik içindeki etiketler 400, 500 ve 600 A’yı gösterir. Burada sonuçlar bu birbirini destekleyen metin ve panel etiketlerine göre 400/500/600 A olarak aktarılmış, açıklama çelişkisi gizlenmemiştir. Ayrıca Şekil 15’te gerçek panel (c) sertlik, panel (d) gözenekliliktir; açıklama ve bazı metin atıfları ters sırayı kullanır. Sayısal değerler panel harflerine körü körüne eşlenmemiştir.

Bulguların desteklediği sonuçlar ve uygulama sınırı

Çalışma, sabit gaz debileri ve toz beslemesi altında 400–600 A akım değişiminin plazma alanını ve parçacık taşınımını birlikte etkilediğini destekler. Sayısal sonuçlarda sıcak ve hızlı bölgeler genişler, hesaplama alanından kaçış azalır ve X = 0,10 m’ye erimiş ulaşan kütle payı artar. Deneyde sıcaklık ve hız ortalamalarının akımla yükselmesi aynı yönlüdür; kaplamalar daha kalın ve daha sert olur. 600 A koşulu bu üçlü karşılaştırmada en düşük gözenekliliği verir.

Çalışma; 600 A’nın evrensel optimum olduğunu, bütün seramik tozlara aynı sayısal sonuçların uygulanacağını, endüstriyel üretimde aynı performansın garanti edileceğini veya daha yüksek akımın sınırsız yarar sağlayacağını desteklemez. Elektriksel yalıtım, uzun süreli aşınma, termal çevrim ömrü, korozyon, yapışma dayanımı ve enerji maliyeti için yeni deneyler raporlanmamıştır. Sıcaklık/hız modeliyle kaplama yapısının birlikte değerlendirilmesi bir optimizasyon zemini sunar; tam çok amaçlı optimizasyon ya da saha yeterlilik deneyi değildir.

Kaynak ve Yöntem Notu

Özgün araştırma:Numerical and Experimental Study of In-flight MgAl2O4 Particles in Atmospheric Plasma Spraying: Effect of Arc Current Variation.

Yazarlar: Byeong-Il Min, Byeongryun Jeon, Hansol Kwon, Yeonwoo Yoo, Do Hyun Kim, Youngjin Park, Yong-Jin Kang ve Hunkwan Park. Yazar sırası PDF ile korunmuştur. Eş birinci yazarlık veya eş katkı beyanı raporlanmamıştır. Sorumlu yazar: Hunkwan Park; kaynak iletişim adresi hkpark@kims.re.kr.

Kurum: Extreme Materials Research Institute, Korea Institute of Materials Science, Changwon, 51508, Republic of Korea. Kaynak türü: Hakem değerlendirmesinden geçmemiş SSRN ön baskısı. Kayıt: 6907963. DOI:10.2139/ssrn.6907963. Resmî bağlantı:SSRN araştırma kaydı.

DOI, resmî SSRN kaydının arama dizinindeki bibliyografik sonucuyla eşleştirilmiştir. Kayıt sayfasının tamamı açılamadığı için yüklenme tarihi, güncel sürüm tarihi, açık lisans veya daha sonra çıkmış bir hakemli sürüm doğrulanamamıştır. PDF’de kesin yayın tarihi, dergi adı, cilt ve sayı bilgisi yer almaz. Bu bilgiler tahmin edilmemiştir. Bilimsel anlatımın tamamı yüklenen PDF’ye dayanır; dış kaynaklardan yeni araştırma sonucu eklenmemiştir.

Finansman: Ministry of Science and ICT tarafından finanse edilen National Research Foundation of Korea Nano & Material Technology Development Program, RS-2022-NR068229; ayrıca Korea Institute of Materials Science Fundamental Research Program, PNKB070. Yazarlar, plazma termofiziksel özellik verileri için Avustralya CSIRO’dan Dr. Murphy’ye teşekkür eder.

Çıkar çatışması: Yazarlar, çalışmayı etkilemiş görünebilecek bilinen mali çıkar çatışması veya kişisel ilişki bulunmadığını beyan eder. Veri erişimi: İncelenen PDF’de bağımsız bir veri erişilebilirliği beyanı, ham veri deposu, çözücü dosyası veya UDF kod bağlantısı görülmemiştir. Lisans: Açık yeniden kullanım lisansı doğrulanmamıştır; özgün kaynak figürleri bu dosyaya kopyalanmamıştır.

Yeniden üretilebilirlik ve raporlama sınırları

Kaynakta altlığın malzemesi ve boyutları, yüzey hazırlığı, püskürtme geçiş sayısı ve tarama hızı, bağımsız kaplama tekrarlarının sayısı, SEM/XRD cihaz modelleri ve ayrıntılı ölçüm ayarları tam olarak raporlanmaz. XRD tarama adımı ve ışınım kaynağı, sertlik baskı sayısı ve bekletme süresi de belirtilmez. Sayısal tarafta ağ bağımsızlığı incelemesi, yakınsama eşikleri, zaman adımı veya örnekleme zamanının ayrıntıları, Rosin–Rammler parametreleri, bütün özgül ısı verileri ve tam UDF uygulaması bulunmaz. Örgü hücre sayısının verilmesi tek başına ağdan bağımsız sonuç garantisi değildir.

Açık hipotez testi, p değeri, güven aralığı veya bağımsız deney tekrarlarına dayalı istatistiksel anlamlılık analizi raporlanmaz. Grafikteki kutu/bıyık gösterimleri dağılım bilgisidir; farklı akımlarda raporlanan ortalamalar arasındaki farkların istatistiksel anlamlılığının yerini tutmaz. Sayısal iki yönlü etkileşimin kaynak terimleri basılı denklemlerde ayrıntılandırılmadığı için yalnız makaledeki denklem setiyle birebir yeniden üretim iddiasında bulunulamaz.

Verianla anlatımı; özgün araştırmanın yerine yeni bir deney veya model sunmaz. Denklem (11)’deki boyut sorunu, Tablo 2 yoğunluk birimi, çap aralığındaki 63/64 µm farkı, geometri mesafesinin referansı, Şekil 13 örnekleme konumu, Şekil 14–15 akım etiketleri, Şekil 15 panel sırası ve kütle kesri toplamları açıkça belirtilmiştir. Bunlar sonuçları otomatik olarak geçersiz ilan etmek için değil, kaynakta neyin raporlandığını ve neyin ayrıca doğrulanması gerektiğini ayırmak için kaydedilmiştir.

Tüm tablolar, yüklenen çalışmada açıkça yazılan sayıları korur; yaklaşık belirtilen kaynak değerleri yaklaşık olarak işaretlenmiştir. Öğretici boyut ve toplam kontrolleri Verianla açıklaması olarak ayrılmıştır. Kaynak grafiklerindeki ham eğriler, hata çubukları, SEM mikroyapıları veya XRD desenleri yapay biçimde yeniden üretilmemiştir. Formüller MathJax uyumlu LaTeX olarak ayrı düğümlerdedir; çalıştırılabilir kod içermez. Canlı site editöründe biçim koruma testi yapılmadığından site düzeyindeki MathJax ve biçimlendirme davranışı doğrulanmış sayılmaz.