Bu çalışma, lityum-iyon batarya elektrotlarındaki etkin malzeme parçacıkları, iletken bağlayıcı alanlar ve elektrolit dolu gözeneklerden oluşan karmaşık üç boyutlu mikroyapıyı büyük ölçekte simüle etmek için GPU ile hızlandırılmış bir yumuşatılmış sınır modeli geliştirmiştir. Yöntem, fazlar arasındaki düzensiz arayüzleri geometriye tam olarak uyan karmaşık bir sonlu eleman ağıyla izlemek yerine, düzenli Kartezyen ağ üzerinde sürekli geçiş gösteren faz fonksiyonlarıyla temsil etmektedir. Böylece etkin malzemede lityum difüzyonu, elektrolitte iyon taşınımı, yük korunumu ve Butler–Volmer elektrokimyasal tepkimeleri aynı hesaplama alanında çözülebilmektedir.

Araştırmacılar yumuşatılmış sınır modelini, geometrik çoklu ağ önkoşullandırmalı eşlenik gradyan çözücüsü ve kırmızı–siyah Gauss–Seidel yumuşatıcısıyla birleştirmiştir. CUDA C++ ile geliştirilen sistemde düzenli şablon işlemleri, birleştirilmiş GPU çekirdekleri, önceki zaman adımının çözümünü başlangıç değeri olarak kullanma ve donanıma göre ayarlanan çoklu ağ seviyeleri uygulanmıştır. Hesaplama bakımından baskın elektrolit potansiyeli denkleminin karşılaştırmasında özel çözücü, 64³ düğümlü ağda AMGX’ten 169 kat, 128³ düğümlü ağda 248 kat daha hızlı çalışmıştır. AMGX 256³ ağda bellek yetersizliği nedeniyle çalışamazken önerilen yöntem aynı ağı V100 GPU üzerinde zaman adımı başına 234 milisaniyede çözmüştür.

MetaX C550’nin 64 GB belleği kullanıldığında 512 × 512 × 512, başka bir ifadeyle yaklaşık 134 milyon düğümlü problem zaman adımı başına 1.314,5 milisaniyede çözülmüştür. Ancak bu süre bütün elektrokimyasal modelin eksiksiz bir deşarj simülasyonuna değil, performans karşılaştırması için ayrılan elektrolit potansiyeli denklemine aittir. Çalışmanın “yaklaşık 200 kat hızlanma” özeti, farklı ağlarda ölçülen 169 ve 248 katlık iki karşılaştırmanın genel ifadesidir.

Model önce 10 mikrometre çapında küresel bir etkin malzeme parçacığının bulunduğu basit bir geometride COMSOL sonlu eleman çözümüyle karşılaştırılmıştır. Etkin malzemedeki lityum yoğunluğu, elektrolitteki iyon yoğunluğu ve elektrolit potansiyeli, fazların iç bölgelerinde referans çözüme yakın bulunmuştur. Yumuşatılmış arayüz çevresinde gözle görülür farklar oluşmuş, ancak ağ inceltildikçe faz içlerindeki farkların azaldığı bildirilmiştir. Çalışma bu karşılaştırma için sayısal bir L₂, L_∞, ortalama hata veya güven aralığı vermemiştir.

Gerçekçi elektrot uygulamasında NMC622 katot; yüzde 57,4 etkin malzeme, yüzde 11 fiziksel CBD ve yüzde 31 toplam gözeneklilikle modellenmiştir. CBD’nin kendi iç gözenekliliği yüzde 50 kabul edildiğinden hesaplama geometrisinde yüzde 22’lik eşdeğer gözenekli CBD bölgesi oluşturulmuştur. Deneysel deşarj noktalarıyla en iyi uyum CBD’nin açıkça temsil edildiği modelde elde edilmiştir. CBD çıkarıldığında model kapasiteyi olduğundan yüksek ve kutuplaşmayı daha düşük tahmin etmiştir.

Alan dağılımları, CBD’nin elektrolit içindeki lityum iyonlarının akım toplayıcı yönüne taşınmasını zorlaştırdığını göstermiştir. 6 mA/cm² akım yoğunluğunda ve 1,25 mAh/cm² kapasitede, etkin malzemenin lityumlanması ayırıcıya yakın bölgede yoğunlaşmış; akım toplayıcı tarafındaki elektrolitte lityum iyonu tükenmesi ve yaklaşık 0,42 voltluk potansiyel farkı oluşmuştur. Bu sonuçlar, CBD’nin yalnızca elektron iletimi sağlayan yararlı bir faz olmadığını; iyon taşınımını engelleyerek yoğunluk ve tepkime kutuplaşmasını artırabileceğini göstermektedir.

Araştırmacılar bu sınırlamayı azaltmak için elektrot boyunca uzanan, 2 mikrometre çapında dört simetrik kanal içeren bir delme tasarımı sınamıştır. 3,2 volt kesme gerilimine göre hesaplanan kapasite artışı 3, 6, 8, 10 ve 12 mA/cm² akım yoğunluklarında sırasıyla yaklaşık yüzde 20, yüzde 58, yüzde 76, yüzde 98 ve yüzde 117 olmuştur. 1 mA/cm² düzeyinde ise delmenin etkisi ihmal edilebilir kalmıştır. Bu değerler sayısal simülasyon sonuçlarıdır; delikli elektrot fiziksel olarak üretilmemiş ve elektrokimyasal hücrede sınanmamıştır.

Türkiye açısından: Geliştirilen yaklaşım, Türkiye’de lityum-iyon hücre, katot, enerji depolama sistemi ve elektrikli araç bataryası geliştiren araştırma merkezleri ile yüksek başarımlı hesaplama gruplarında gerçek elektrot görüntülerinden sayısal tasarım yapılması için değerlendirilebilir. Yöntem; parçacık büyüklüğü, CBD dağılımı, gözeneklilik ve lazerle oluşturulacak iyon taşıma kanallarının üretimden önce taranmasına yardımcı olabilir. Ancak uygulamadan önce Türkiye’de kullanılan etkin malzeme kimyaları, bağlayıcı ve iletken katkı oranları, kalınlıklar, presleme koşulları ve gerçek üç boyutlu tomografi yapılarıyla yeniden kurulmalıdır. Model çıktıları bağımsız hücre deneyleriyle doğrulanmalı; delme işleminin etkin malzeme kaybı, hacimsel enerji yoğunluğu, mekanik dayanım, üretim maliyeti ve çevrim ömrü üzerindeki etkileri birlikte ölçülmelidir. Çalışmadan Türkiye’de üretilecek belirli bir elektrotun kapasitesinin yüzde 117 artacağı, belirli bir GPU’nun aynı hızlanmayı sağlayacağı veya yöntemin batarya güvenliğini doğrudan iyileştireceği sonucu çıkarılamaz.

Araştırmanın çözmeye çalıştığı problem nedir?

Lityum-iyon batarya elektrotları makroskopik ölçekte düzgün bir tabaka gibi görünse de mikrometre ölçeğinde üç temel bileşenden oluşan oldukça düzensiz yapılardır:

• Elektrik enerjisinin depolandığı etkin malzeme parçacıkları
• Elektron iletimini ve parçacıkların mekanik bağlanmasını sağlayan iletken bağlayıcı alanlar
• Elektrolitin doldurduğu ve lityum iyonlarının taşındığı gözenek ağı

Parçacık büyüklüğü dağılımı, parçacıkların birbirine temas biçimi, CBD kümelenmesi, gözeneklerin bağlantısı ve elektrot kalınlığı; iyonların hangi yolları izleyeceğini ve hangi bölgelerde tepkime gerçekleşeceğini belirler. Aynı ortalama gözenekliliğe sahip iki elektrot, gözeneklerin ve CBD’nin uzaysal dağılımı farklı olduğunda aynı elektrokimyasal performansı göstermeyebilir.

P2D gibi homojenleştirilmiş modeller elektrodu ortalama büyüklüklerle temsil eder. Bu modeller hücre düzeyindeki gerilim ve kapasiteyi düşük maliyetle hesaplayabilir; ancak yerel iyon tükenmesini, parçacıklar arasındaki farklı tepkime hızlarını, CBD kümelerinin oluşturduğu engelleri ve karmaşık gözenek yollarını açıkça çözemez.

Tam alanlı heterojen modeller ise gerçek üç boyutlu yapıyı çözebilir. Bunun karşılığında etkin malzeme, CBD ve gözeneklerin düzensiz sınırlarına uygun bir ağ oluşturulması gerekir. On milyonlarca düğüm içeren bir elektrot yapısında geometriye uyumlu sonlu eleman ağı hazırlamak zaman alabilir, başarısız elemanlar üretebilir ve çözümün hesaplama maliyetini büyütebilir.

Yumuşatılmış sınır yöntemi hangi yaklaşımı kullanmaktadır?

Yumuşatılmış sınır yöntemi, fazlar arasındaki keskin sınırları ağın yüzeyleriyle bire bir eşleştirmez. Bunun yerine faz kimliğini 0 ile 1 arasında değişen sürekli fonksiyonlarla tanımlar. Çalışmada iki temel gösterge kullanılmıştır:

• \(\psi(\mathbf{r})\): etkin malzeme ile etkin malzeme dışındaki bölgeleri ayırır.
• \(\gamma(\mathbf{r})\): CBD ile CBD dışındaki bölgeleri ayırır.

Bir düğüm fazın içindeyse gösterge yaklaşık 1, dışındaysa yaklaşık 0, arayüz bölgesindeyse 0 ile 1 arasında değer alır. Malzeme özellikleri bu değerler kullanılarak düğüm düzeyinde ağırlıklandırılır.

Bu yapı sayesinde etkin malzeme, CBD ve elektrolit için ayrı ağlar hazırlanması gerekmez. Aynı düzenli Kartezyen ağ bütün hesaplama alanını kapsar. Geometrik karmaşıklık ağ bağlantılarında değil, \(\psi\) ve \(\gamma\) alanlarında saklanır.

Şekil 1’deki iki elektrokimyasal arayüz neyi temsil etmektedir?

Şekil 1, bir etkin malzeme parçacığını elektrolit ve CBD ile birlikte göstermektedir. Araştırmacılar iki farklı etkin malzeme–elektrolit arayüzü tanımlamıştır:

• Arayüz I: Etkin malzemenin doğrudan saf elektrolit gözenek alanıyla temas ettiği yüzey.
• Arayüz II: Etkin malzemenin, elektrolit içeren gözenekli CBD üzerinden temas ettiği yüzey.

CBD tamamen katı bir karbon tabakası olarak değil, kendi içinde elektrolit bulunan gözenekli bir bölge olarak modellenmiştir. Buna rağmen CBD içindeki iyon taşınımı saf elektrolit bölgesine göre daha yavaştır ve etkin malzeme yüzeyinin yalnızca belirli bir bölümü elektrokimyasal olarak erişilebilir kabul edilmiştir.

Bu nedenle iki arayüz için tepkime sabitleri farklılaştırılmıştır. Doğrudan elektrolit temasında \(k_I=k_0\), CBD üzerinden temas eden bölgede ise \(k_{II}=(1-\gamma)k_0\) ilişkisi kullanılmıştır.

Etkin malzemedeki lityum taşınımı nasıl modellenmiştir?

Etkin malzeme içindeki lityum yoğunluğu Fick’in ikinci yasasıyla ifade edilmiştir:

\[ \frac{\partial c_s}{\partial t} = \nabla\cdot(D_s\nabla c_s) \]

• \(c_s\), etkin malzeme içindeki lityum yoğunluğudur.
• \(D_s\), etkin malzemedeki lityum difüzyon katsayısıdır.
• \(t\), zamandır.

Etkin malzeme–elektrolit sınırındaki akı, elektrokimyasal tepkime hızıyla ilişkilendirilmiştir:

\[ -D_s\nabla c_s\cdot\mathbf{n} = \frac{j}{F} \]

• \(\mathbf{n}\), etkin malzeme dışına yönelen birim normaldir.
• \(j\), yerel elektrokimyasal tepkime akı yoğunluğudur.
• \(F\), Faraday sabitidir.

Yumuşatılmış sınır dönüşümü sonrasında hacim denklemi ve arayüz sınır koşulu tek bir eşitlikte birleştirilmiştir:

\[ \psi\frac{\partial c_s}{\partial t} = \nabla\cdot \left( \psi D_s^{*}\nabla c_s \right) - |\nabla\psi|\frac{j}{F} \]

Bu denklemde \(|\nabla\psi|\) yalnızca yumuşatılmış faz sınırında anlamlı büyüklüğe ulaştığı için yüzey tepkimesi ayrı bir geometrik sınır listesi tanımlanmadan hesaplama alanına eklenmektedir.

Elektrolitteki iyon taşınımı nasıl ele alınmıştır?

Elektrolit iki farklı bölgeye ayrılmıştır. Saf gözeneklerde elektrolit hacim oranı 1, CBD içindeki gözeneklerde ise \(\varepsilon_{\mathrm{CBD}}\) olarak alınmıştır.

Etkin iyon difüzyon katsayısı:

\[ D_{\mathrm{eff}}=D_e\varepsilon^b \]

biçimindedir. Burada \(D_e\) saf elektrolitteki difüzyon katsayısı, \(\varepsilon\) yerel gözeneklilik ve \(b\) Bruggeman katsayısıdır. Çalışmada \(b=1{,}5\) kullanılmıştır.

CBD bölgesinde gözenekliliğin azalması, etkin difüzyon ve iyon iletkenliğinin düşmesine yol açar. Yumuşatılmış modelde yerel gözeneklilik:

\[ \varepsilon^{*} = 1-\gamma+\gamma\varepsilon_{\mathrm{CBD}} \]

şeklinde ağırlıklandırılmıştır. CBD dışında \(\gamma=0\) olduğunda değer 1’e, CBD içinde \(\gamma=1\) olduğunda \(\varepsilon_{\mathrm{CBD}}\)’ye yaklaşır.

Elektrokimyasal tepkime nasıl hesaplanmıştır?

Etkin malzeme yüzeyindeki lityum yerleştirme tepkimesi Butler–Volmer denklemiyle modellenmiştir:

\[ j = j_0 \left[ \exp\left( \frac{\alpha_aF\eta}{RT} \right) - \exp\left( -\frac{\alpha_cF\eta}{RT} \right) \right] \]

Aşırı potansiyel:

\[ \eta = \phi_s-\phi_e-U(c_s) \]

ve değişim akım yoğunluğu:

\[ j_0 = Fk^{*} c_e^{\alpha_c} c_s^{\alpha_c} (c_{s,\max}-c_s)^{\alpha_c} \]

olarak tanımlanmıştır.

• \(\phi_s\), katı faz potansiyelidir.
• \(\phi_e\), elektrolit potansiyelidir.
• \(U(c_s)\), yoğunluğa bağlı açık devre potansiyelidir.
• \(c_{s,\max}\), etkin malzemenin en yüksek lityum yoğunluğudur.
• \(\alpha_a\) ve \(\alpha_c\), anoduk ve katoduk aktarım katsayılarıdır.

Elektron taşınımının iyon taşınımından çok hızlı olduğu varsayılarak etkin malzeme parçacıkları eş potansiyelli kabul edilmiştir. Katı faz içinde ayrı bir uzaysal potansiyel denklemi çözülmemiş; bütün yüzeylerdeki tepkime akılarının uygulanan toplam akıma eşit olması şartı kullanılmıştır:

\[ \int_{\Gamma}j\,ds=I_{\mathrm{app}} \]

Elektrolit potansiyeli neden hesaplama bakımından önemlidir?

Elektrolit içindeki iyonik akım, potansiyel gradyanı ile yoğunluk gradyanının birlikte etkisiyle hesaplanmıştır:

\[ \nabla\cdot (\kappa_{\mathrm{eff}}\nabla\phi_e) + \nabla\cdot (\kappa_D\nabla\ln c_e) = 0 \]

Yumuşatılmış sınır dönüşümünde yüzey akısı denkleme gömülmüş ve aşağıdaki tam alan biçimi elde edilmiştir:

\[ \nabla\cdot \left[ (1-\psi)\kappa_e^{*}\nabla\phi_e \right] + \nabla\cdot \left[ (1-\psi)\kappa_D^{*}\nabla\ln c_e \right] = -j \]

Bu Poisson benzeri denklem, her zaman adımında büyük ve seyrek bir doğrusal denklem sistemi oluşturur. Çalışmanın AMGX karşılaştırması, bütün elektrokimyasal denklemler yerine özellikle bu hesaplama bakımından baskın denklem üzerinde yapılmıştır.

Zaman ve uzay ayrıklaştırması nasıl yapılmıştır?

Araştırmacılar zaman için birinci dereceden geri fark, uzay için düzenli ağ üzerinde ikinci dereceden merkezî fark şeması kullanmıştır.

Her yeni zaman adımında işlem sırası şöyledir:

1. Butler–Volmer denklemi ve galvanostatik toplam akım şartı kullanılarak yeni tepkime akısı \(j^{n+1}\) hesaplanır.
2. Etkin malzeme yoğunluğu \(c_s^{n+1}\) çözülür.
3. Elektrolit yoğunluğu \(c_e^{n+1}\) çözülür.
4. Güncellenmiş yoğunluk alanı kullanılarak elektrolit potansiyeli \(\phi_e^{n+1}\) çözülür.

Denklemler tamamen eş zamanlı ve monolitik olarak çözülmemiş, ardışık biçimde ayrıştırılmıştır. Bu seçim hesaplama maliyetini azaltabilir; ancak çok güçlü eşleşmelerde zaman adımı hassasiyetinin ayrıca incelenmesi gerekir. Çalışma farklı zaman adımlarına ilişkin kapsamlı bir duyarlılık tablosu sunmamıştır.

GPU çözücüsünün temel mimarisi nedir?

Şekil 2, geometrik çoklu ağ ile eşlenik gradyan yönteminin birlikte kullanımını göstermektedir. Çoklu ağ yöntemi farklı uzaysal frekanstaki hataları farklı ağ seviyelerinde azaltır:

• İnce ağlar yüksek frekanslı yerel hataları işler.
• Kaba ağlar yavaş değişen düşük frekanslı hataları azaltır.
• Artık değer kaba ağa aktarılır.
• Kaba ağ düzeltmesi yeniden ince ağa enterpole edilir.

Tek başına çoklu ağ yöntemi, fazlar arasındaki yüksek özellik farklarında kararsızlaşabileceği için ana çözücü olarak eşlenik gradyan yöntemi kullanılmıştır. Geometrik çoklu ağ, eşlenik gradyan sisteminin önkoşullandırıcısıdır.

Her ağ düzeyindeki yumuşatma işlemi kırmızı–siyah Gauss–Seidel yöntemiyle gerçekleştirilmiştir. Düğümler satranç tahtasındaki gibi iki gruba ayrılmıştır. Aynı renkteki düğümler birbirine doğrudan bağımlı olmadığı için paralel olarak güncellenebilmektedir.

GPU’ya özgü hangi iyileştirmeler uygulanmıştır?

• Şablon tabanlı matris-vektör çarpımı: Yapılandırılmış ağda komşu düğüm indeksleri çalışma sırasında hesaplanmış, seyrek matris sütun indekslerinin bellekte saklanması önlenmiştir.
• GPU çekirdeği birleştirme: Artık hesabı, vektör güncellemesi ve bazı iç çarpımlar tek çekirdek çağrısında birleştirilmiştir.
• Bellek geçişlerinin azaltılması: Aynı alanların küresel GPU belleğinden tekrar tekrar okunması sınırlandırılmıştır.
• Renk temelli erken dönüş: Her düğüm için bir iş parçacığı başlatılmış, yalnızca işlenen renge ait düğümler güncellenmiştir. Bu seçim düzenli bellek erişimini korumuştur.
• Sabit komşu sayısı: İki boyutta dört, üç boyutta altı komşu derleme zamanında sabitlenerek gereksiz dallanmalar kaldırılmıştır.
• Gereksiz eş zamanlama çağrılarının kaldırılması: CUDA’nın varsayılan akış sıralaması kullanılarak her tarama sonundaki fazladan beklemeler azaltılmıştır.
• Uyarlamalı çoklu ağ derinliği: Ağ boyutuna göre farklı sayıda kaba seviye kullanılmıştır.
• Sıcak başlangıç: Önceki zaman adımının çözümü bir sonraki adımın başlangıç tahmini yapılmıştır.

Geliştirilen çözüm CUDA C++ ile uygulanmış ve araştırmacıların Advanced Battery Designer, kısaca ABD adını verdikleri yazılıma eklenmiştir. Çalışmada bu yazılımın veya kaynak kodun herkese açık bir deposu verilmemiştir.

Basit geometrili doğrulama nasıl hazırlanmıştır?

Doğrulama geometrisi, 50 × 50 × 50 mikrometrelik küp içinde merkezlenmiş 10 mikrometre çapında küresel etkin malzeme parçacığından oluşmaktadır.

Yumuşatılmış sınır çözümü için üç ağ kullanılmıştır:

• 50 × 50 × 50 düğüm
• 100 × 100 × 100 düğüm
• 200 × 200 × 200 düğüm

Sistem 3C deşarj oranında 300 saniye çalıştırılmıştır. Kürenin yüzeyinde düzgün tepkime akısı, x = 50 mikrometre sınırında uygulanan akıma eşit iyon akısı ve sıfır volt referans potansiyeli uygulanmıştır. Diğer dış sınırlar için sıfır akılı ayna koşulları kullanılmıştır.

Şekil 3’te kırmızı referans çizgisi boyunca üç büyüklük karşılaştırılmıştır:

• Etkin malzemedeki lityum yoğunluğu
• Elektrolitteki lityum iyonu yoğunluğu
• Elektrolit potansiyeli

COMSOL eğrileriyle yumuşatılmış sınır noktaları faz içlerinde büyük ölçüde örtüşmektedir. Arayüzde ise yumuşatılmış model alanı birkaç ağ hücresine yaydığı için keskin arayüzlü COMSOL çözümünden ayrılmaktadır. Ağ inceltildikçe faz içindeki sonuçların referansa yaklaştığı görülmektedir.

Çalışma “düşük” veya “ihmal edilebilir” hata ifadelerini kullanmasına rağmen hata normlarını tablo hâlinde vermemiştir. Bu nedenle doğrulamanın nicel hata sınırı yeniden hesaplanamamaktadır.

AMGX karşılaştırması nasıl yapılmıştır?

Performans deneyi, 512 × 512 × 512 mikrometrelik bir küp ve merkezinde 256 mikrometre çapında küresel bölge kullanılarak hazırlanmıştır. Küre yüzeyine 10 mA/cm² sabit akım yoğunluğu uygulanmıştır.

Her test 100 kez tekrarlanmış ve ortalama zaman adımı süresi raporlanmıştır. Yakınsama ölçütü:

\[ \frac{\|r\|}{\|b\|}<10^{-8} \]

olarak belirlenmiştir.

AMGX genel amaçlı bir doğrusal cebir kütüphanesidir. Önerilen çözücü ise sabit komşuluğa sahip düzenli SBM ağı için özel olarak tasarlanmıştır. Bu nedenle karşılaştırma iki genel amaçlı eşdeğer uygulamadan çok, probleme özel çözücü ile genel amaçlı çözücü arasındaki farkı göstermektedir.

Farklı GPU’larda hangi sonuçlar elde edilmiştir?

A100, belleğe sığan üç ağda en düşük süreleri sağlamıştır. MetaX C550 küçük problemlerde daha yavaş olmakla birlikte 256³ ağda A100’e yaklaşmıştır. En büyük problem yalnızca 64 GB belleğe sahip C550’ye sığmıştır.

512³ ağ yaklaşık 134.217.728 düğüm içerir. Zaman adımı başına 1,3145 saniyelik sonuç, büyük problemin tek GPU belleğinde çözülebildiğini göstermektedir. Ancak tam deşarj analizi binlerce veya daha fazla zaman adımı gerektirebilir; dolayısıyla bu değer tek başına bütün hücre simülasyonunun gerçek zamanlı tamamlandığını göstermez.

Gerçekçi elektrot geometrisi nasıl oluşturulmuştur?

Elektrot yapısı daha önce yayımlanan deneysel morfoloji bilgileri kullanılarak yeniden kurulmuştur. Etkin malzeme parçacıkları ayrık eleman yöntemiyle üretilmiş, CBD geometrisi MatBox aracıyla eklenmiştir.

Parçacık büyüklükleri şekil parametresi 2,62 ve ölçek parametresi 2,17 olan Gamma dağılımına dayandırılmıştır. Bildirilen ortalama çap 5,69 mikrometre, standart sapma 3,51 mikrometredir.

Şekil 4’te etkin malzeme yeşil, CBD kırmızı ve saf elektrolit gözenekleri mavi gösterilmiştir. CBD’nin parçacıklar arasında düzensiz ve bağlantılı bir ağ oluşturduğu görülmektedir.

Gerçekçi elektrot modelinde hangi malzeme parametreleri kullanılmıştır?

Tabloda elektromigrasyon katsayısı \(t_+\) için sayısal değer verilmemiştir. Kinetik parametrelerin deneysel gerilim sonuçlarına uyum sağlamak amacıyla kalibre edildiği belirtilmiştir.

Deneysel deşarj eğrileriyle karşılaştırma nasıl yapılmıştır?

Model 1–12 mA/cm² aralığında galvanostatik olarak çalıştırılmıştır. Çalışmada yaklaşık 3 mA/cm² değerinin 1C’ye karşılık geldiği belirtilmiştir.

Şekil 5’te:

• Noktalar daha önce yayımlanmış deneysel ölçümleri,
• Sürekli çizgiler yumuşatılmış sınır modeli sonuçlarını göstermektedir.

CBD’nin açık biçimde temsil edildiği model, 1, 3, 6, 8, 10 ve 12 mA/cm² akım yoğunluklarında deneysel gerilim–kapasite eğrilerinin genel biçimini izlemiştir. Akım yükseldikçe kapasite ve çalışma geriliminin azalması model tarafından yeniden üretilmiştir.

Bununla birlikte çalışma:

• Ortalama mutlak hata, RMSE veya başka bir sayısal uyum ölçüsü vermemiştir.
• Deneysel tekrarların dağılımını veya hata çubuklarını göstermemiştir.
• Kinetik parametreleri aynı eğrilere göre kalibre etmiştir.
• Ayrı bir kalibrasyon ve bağımsız doğrulama veri kümesi tanımlamamıştır.

Bu nedenle karşılaştırma modelin deney eğrilerini temsil edebildiğini göstermektedir; bütünüyle kör ve bağımsız bir öngörü doğrulaması değildir.

CBD modelden çıkarıldığında ne olmuştur?

CBD’siz karşılaştırma, gerçek bir üretilebilir elektrodu temsil etmek amacıyla yapılmamıştır. Modelde \(\gamma(\mathbf{r})\) her yerde sıfıra ayarlanmış ve CBD alanı saf elektrolit olarak değerlendirilmiştir.

Bu değişiklik sonucunda:

• İyonların elektrot kalınlığı boyunca ilerlemesi kolaylaşmıştır.
• Elektrolit yoğunluğu gradyanı azalmıştır.
• Tepkime için gereken aşırı potansiyel düşmüştür.
• Model deneysel ölçümlerden daha yüksek kapasite tahmin etmiştir.
• Gerilim kutuplaşması olduğundan düşük hesaplanmıştır.

CBD iki mekanizmayla kutuplaşmayı artırmaktadır. Birincisi, parçacıklar arasındaki gözenek hacmini işgal ederek iyonların yolunu uzatması ve etkin difüzyonu azaltmasıdır. İkincisi, etkin malzeme yüzeyinin bir bölümünü kapatarak doğrudan elektrolit temasını ve tepkimeye açık yüzey alanını azaltmasıdır.

CBD’nin nano-gözenekli yapısı doğrudan çözülmüş müdür?

Hayır. CBD içindeki gerçek gözenekler onlarca veya yüzlerce nanometre ölçeğinde olabilir. Bütün elektrot mikroyapısıyla birlikte bu gözenekleri doğrudan çözmek, ağ sayısını çok büyük ölçüde artıracaktır.

Bu nedenle CBD, eşdeğer gözenekliliğe ve etkin taşıma özelliklerine sahip homojen bir gözenekli ortam olarak temsil edilmiştir. Çalışmada kullanılan Bruggeman katsayısı ve yüzde 50 gözeneklilik yaklaşık 1,4’lük düşük bir kıvrımlılık etkisi oluştururken, kaynak alınan deneysel çalışmada CBD’nin etkin iyon iletkenliğinin saf elektrolitin yalnızca yüzde 1,78’i olduğu bildirilmiştir.

Araştırmacılar bu farkı, normalden düşük seçilen 2,1 × 10⁻¹¹ m²/s elektrolit difüzyon katsayısıyla dolaylı biçimde temsil etmiştir. Dolayısıyla model CBD’nin gerçek gözenek geometrisini çözmemekte; onun ortalama engelleme etkisini kalibre edilmiş bir katsayıyla içermektedir.

6 mA/cm² koşulundaki alan dağılımları ne göstermektedir?

Şekil 6, 6 mA/cm² akım yoğunluğunda 750 saniye deşarj sonrasında, 1,25 mAh/cm² kapasitedeki alanları göstermektedir.

CBD içeren durumda:

• Ayırıcıya yakın etkin malzeme parçacıkları daha fazla lityumlanmıştır.
• Akım toplayıcı tarafındaki parçacıkların önemli bir bölümü daha az etkin kalmıştır.
• Akım toplayıcı yakınındaki elektrolit lityumu neredeyse tükenmiştir.
• Ayırıcı tarafında iyon birikimi oluşmuştur.
• Elektrolit boyunca yaklaşık 0,42 V potansiyel farkı hesaplanmıştır.

Lityum iyonları ayırıcı tarafından elektrota girmekte ve ilerlerken etkin malzeme tarafından tüketilmektedir. Akım toplayıcıya yakın bölgeye ulaşabilmek için bütün elektrot kalınlığını geçmeleri gerekir. CBD bu yolu daha da zorlaştırdığında uzak bölgelerin iyonla beslenmesi yetersiz kalmaktadır.

Aynı toplam akımın sürdürülebilmesi için sistem daha yüksek aşırı potansiyel üretmek zorunda kalmakta ve hücre gerilimi daha hızlı düşmektedir.

Delikli elektrot tasarımı nasıl oluşturulmuştur?

Araştırmacılar elektrot boyunca x yönünde uzanan dört simetrik kanal tanımlamıştır. Her kanalın çapı 2 mikrometredir. Kanalların bulunduğu bölgelerde hem etkin malzeme hem CBD kaldırılmış ve alan elektrolitle dolu taşıma yolu hâline getirilmiştir.

Şekil 7’de kanallar, farklı kesitlerde devam eden dört mavi silindir olarak görülmektedir. Tasarımın amacı iyonların kıvrımlı CBD ve gözenek ağı içinden geçmek yerine ayırıcıdan akım toplayıcı yönüne daha doğrudan ilerlemesini sağlamaktır.

Delme tasarımı kapasiteyi nasıl değiştirmiştir?

Şekil 8’de delikli ve deliksiz elektrotların simülasyon eğrileri karşılaştırılmıştır. Kesme gerilimi 3,2 V olarak alındığında:

Düşük akımda iyon taşınımı dengeye yakın kaldığı için doğrudan kanalların ek yararı sınırlıdır. Akım yükseldikçe taşıma sınırlaması büyümekte ve kanalların sağladığı kısa yol daha önemli hâle gelmektedir.

Yüzde 117 artış, deliksiz yüksek hızlı elektrotun düşük kalan kapasitesine göre hesaplanan göreli bir artıştır. Elektrotun teorik kapasitesinin yüzde 117 arttığı veya enerji yoğunluğunun iki katına çıktığı anlamına gelmez.

Delme önerisinin temel ödünleşimi nedir?

Kanallar iyon taşınımını iyileştirirken etkin malzeme ve CBD hacmini azaltmaktadır. Delik sayısı veya çapı aşırı artırılırsa:

• Alan ve hacim başına etkin malzeme yüklemesi düşebilir.
• Teorik kapasite azalabilir.
• Elektron iletim ağı kesintiye uğrayabilir.
• Elektrotun mekanik dayanımı değişebilir.
• Üretim süresi ve maliyeti artabilir.

Çalışmada yalnızca dört adet 2 mikrometrelik kanal değerlendirilmiştir. Delik çapı, aralığı, deseni ve toplam çıkarılan hacim optimize edilmemiştir. Mekanik, termal veya çevrim ömrü analizi yapılmamıştır.

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

• Etkin malzeme, CBD ve elektroliti ayrı fiziksel özelliklere sahip üç faz olarak çözmektedir.
• İki farklı etkin malzeme–elektrolit arayüzünü tepkime kinetiğine dâhil etmektedir.
• Geometriye uyumlu karmaşık ağ gereksinimini azaltan düzenli bir hesaplama alanı kullanmaktadır.
• Yöntem basit geometrili keskin arayüzlü FEM çözümüyle karşılaştırılmıştır.
• Farklı ağ çözünürlüklerinde yakınsama davranışı görsel olarak incelenmiştir.
• AMGX ile aynı yakınsama ölçütünde performans karşılaştırması yapılmıştır.
• V100, A100 ve MetaX C550 olmak üzere üç GPU platformu değerlendirilmiştir.
• Yaklaşık 134 milyon düğümlü bir problem tek GPU üzerinde çalıştırılmıştır.
• Gerçek elektrot özelliklerine dayalı üç boyutlu AM–CBD–gözenek geometrisi kullanılmıştır.
• Farklı akım yoğunluklarındaki deneysel deşarj eğrileriyle karşılaştırma yapılmıştır.
• CBD’nin makroskopik gerilim eğrisine etkisi yerel yoğunluk ve potansiyel alanlarıyla açıklanmıştır.
• Delikli elektrot tasarımı yüksek hızlı deşarj koşullarında karşılaştırmalı olarak sınanmıştır.

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

• Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.
• Doğrulama grafiklerinde nicel hata normları raporlanmamıştır.
• Yumuşatılmış arayüz çevresindeki farklar “ihmal edilebilir” değildir ve açıkça görülmektedir.
• Performans karşılaştırması bütün çok alanlı zaman adımından çok elektrolit potansiyeli denklemine odaklanmıştır.
• AMGX karşılaştırması probleme özel çözücü ile genel amaçlı çözücüyü karşılaştırmaktadır.
• CPU tabanlı FEM veya başka bir tam alanlı elektrot yazılımıyla uçtan uca süre karşılaştırması verilmemiştir.
• Elektrik enerjisi tüketimi veya GPU maliyet analizi yapılmamıştır.
• Gerçekçi elektrot yapısı doğrudan tomografi hacminden değil, dağılımlara göre üretilen sentetik geometriden oluşturulmuştur.
• CBD’nin nano-gözenekli yapısı doğrudan çözülmemiş, homojen gözenekli ortam kabul edilmiştir.
• CBD kıvrımlılığı kalibre edilmiş düşük elektrolit difüzyon katsayısı içinde dolaylı biçimde temsil edilmiştir.
• Kinetik parametreler karşılaştırılan deneysel eğrilere göre ayarlanmıştır.
• Bağımsız kalibrasyon ve test veri kümeleri tanımlanmamıştır.
• Deneysel eğriler için hata çubukları ve sayısal uyum ölçütleri verilmemiştir.
• Delikli elektrot tasarımı deneysel olarak üretilip sınanmamıştır.
• Delik sayısı, çapı ve konumu optimize edilmemiştir.
• Delmenin mekanik dayanım ve hacimsel enerji yoğunluğu üzerindeki etkisi hesaplanmamıştır.
• Etkin malzeme genleşmesi, parçacık çatlaması ve hareketi modele dâhil değildir.
• SEI büyümesi, lityum kaplanması ve kapasite kaybı modellenmemiştir.
• Isı üretimi ve sıcaklık dağılımı çözülmemiştir.
• Çevrim yaşlanması veya güvenlik olayı simüle edilmemiştir.
• Tek GPU kullanılmış; çoklu GPU ölçeklenmesi gelecek çalışmaya bırakılmıştır.
• Kaynak kodu ve tam veri kümesi açık bir depoda yayımlanmamıştır.

Çalışma neyi desteklemektedir?

• Yumuşatılmış sınır yöntemi, düzensiz çok fazlı elektrot geometrilerini düzenli ağ üzerinde temsil edebilir.
• Basit küresel testte faz içindeki yoğunluk ve potansiyel sonuçları COMSOL çözümüne yakın kalmıştır.
• Probleme özel çoklu ağ–eşlenik gradyan GPU çözücüsü, incelenen potansiyel denkleminde AMGX’ten belirgin biçimde hızlıdır.
• 134 milyon düğümlü potansiyel problemi yeterli GPU belleği bulunduğunda yaklaşık 1,3 saniyede çözülebilmiştir.
• CBD’nin açıkça temsil edilmesi, incelenen NMC622 elektrotun deneysel deşarj eğrilerinin yeniden üretilmesini geliştirmiştir.
• CBD, modelde iyon taşınımını sınırlandırmış ve elektrokimyasal kutuplaşmayı artırmıştır.
• Dört adet doğrudan kanal, yüksek akım yoğunluğundaki sayısal kapasiteyi deliksiz duruma göre artırmıştır.
• Yöntem geniş parametrik tarama ve elektrot geometrisi karşılaştırmaları için kullanılabilecek hesaplama kapasitesi sunmaktadır.

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

• Bütün lityum-iyon elektrotların 200 kat daha hızlı simüle edileceğini kanıtlamaz.
• Bütün elektrokimyasal zaman adımının 134 milyon düğümde 1,3 saniyede tamamlandığını göstermez.
• Yumuşatılmış sınır yönteminin her geometri ve çözünürlükte FEM kadar doğru olduğunu kanıtlamaz.
• Arayüz yakınındaki hataların bütün uygulamalar için önemsiz olduğunu göstermez.
• CBD’nin her elektrot kimyasında aynı ölçüde kutuplaşma oluşturacağını göstermez.
• Kalibre edilen parametrelerin başka bir NMC622 elektrotuna doğrudan aktarılabileceğini kanıtlamaz.
• Dört delikli tasarımın gerçek hücrede yüzde 20–117 kapasite artışı sağlayacağını göstermez.
• Delikli elektrotun hacimsel enerji yoğunluğunu veya çevrim ömrünü geliştirdiğini kanıtlamaz.
• Yöntemin termal kaçak, yangın veya başka batarya güvenliği sorunlarını tahmin ettiğini göstermez.
• Advanced Battery Designer yazılımının bağımsız kullanıcılar tarafından yeniden üretildiğini göstermez.
• Tek bir tasarımın üretim için en uygun elektrot mimarisi olduğunu belirlemez.

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

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

Basit geometri doğrulamasının özeti

GPU performans deneyinin koşulları

Performans sonuçlarının doğru yorumlanması

CBD etkisinin mekanizma özeti

Delikli elektrot sonuçlarının sınırları

Şekillerin temel mesajları

• Şekil 1: Etkin malzeme, CBD ve elektrolit fazları ile doğrudan ve CBD üzerinden oluşan iki arayüzü göstermektedir.
• Şekil 2: Çoklu ağ V çevrimini ve kırmızı–siyah düğüm gruplandırmasını kullanan GPU çözücü mimarisini göstermektedir.
• Şekil 3: Küresel etkin malzeme testinde SBM sonuçlarının üç ağ çözünürlüğünde COMSOL eğrilerine yaklaşmasını göstermektedir.
• Şekil 4: Ayrık eleman yöntemiyle oluşturulan etkin malzeme parçacıklarını ve AM–CBD–gözenek kesitini göstermektedir.
• Şekil 5: CBD içeren modelin deneysel deşarj eğrilerini CBD’siz modelden daha iyi izlediğini göstermektedir.
• Şekil 6: CBD’nin etkin malzeme lityumlanmasını, elektrolit yoğunluğu gradyanını ve potansiyel farkını artırmasını göstermektedir.
• Şekil 7: Elektrot kalınlığı boyunca oluşturulan dört simetrik iyon taşıma kanalını göstermektedir.
• Şekil 8: Delmenin yüksek akım yoğunluklarında kutuplaşmayı azaltıp hesaplanan kapasiteyi artırmasını göstermektedir.

Kanıt düzeyine göre bulguların durumu

Kaynak ve Yöntem Notu

Çalışmanın tam özgün adı: A GPU-Accelerated Smoothed Boundary Model for Large-Scale Heterogeneous Multiphase Electrochemical Simulations of Lithium-Ion Electrodes

Yazarlar: Yunkai Wu, Wangwang Kuang, Ling Xiao, Jin Dai, Kui Gong, Xushan Zhao ve Yin Wang.

Yazar sıralaması: Yunkai Wu birinci, Wangwang Kuang ikinci, Ling Xiao üçüncü, Jin Dai dördüncü, Kui Gong beşinci, Xushan Zhao altıncı ve Yin Wang yedinci yazardır.

Eşit katkı: Yunkai Wu ve Wangwang Kuang eşit katkıda bulunmuştur.

Sorumlu yazarlar:

• Wangwang Kuang: wwkuang@foxmail.com
• Kui Gong: gongkui@hzwtech.com
• Xushan Zhao: ZhaoXS@catl-21c.com

Kurumsal bağlantılar:

1. Hongzhiwei Technology (Shanghai) CO.LTD., Shanghai 200120, Çin.
2. HZW-CATL Lab, Shanghai 200120, Çin.
3. 21C Innovation Laboratory, Contemporary Amperex Technology Ltd. (CATL), Ningde, Çin.

Yazar–kurum eşleştirmesi:

• Yunkai Wu: Hongzhiwei Technology ve HZW-CATL Lab.
• Wangwang Kuang: Hongzhiwei Technology ve HZW-CATL Lab.
• Ling Xiao: Hongzhiwei Technology ve HZW-CATL Lab.
• Jin Dai: CATL 21C Innovation Laboratory.
• Kui Gong: Hongzhiwei Technology ve HZW-CATL Lab.
• Xushan Zhao: HZW-CATL Lab ve CATL 21C Innovation Laboratory.
• Yin Wang: Hongzhiwei Technology ve HZW-CATL Lab.

DOI:10.2139/ssrn.7197194

Yayın platformu: SSRN.

Resmî kaynak bağlantısı:SSRN çalışma kaydı

SSRN yükleme tarihi: 28 Temmuz 2026.

Sayfa sayısı: 34.

Kaynak türü: Çok fazlı elektrokimyasal model geliştirme, GPU hızlandırma, sayısal doğrulama ve elektrot tasarımı uygulaması içeren preprint.

Hakemlik durumu: Çalışma hakem değerlendirmesinden geçmemiştir. Her sayfada “This preprint research paper has not been peer reviewed” uyarısı bulunmaktadır.

Dergi: Hakemli bir dergi adı, kabul tarihi, cilt, sayı veya nihai makale numarası yer almamaktadır.

Özgün yayınevi: Nihai hakemli yayın için yayınevi bilgisi bulunmamaktadır. İncelenen sürüm SSRN preprintidir.

Finansman: Çalışma, Çin Halk Cumhuriyeti Sanayi ve Bilgi Teknolojileri Bakanlığının 2024ZD0607100 numaralı fonuyla desteklenmiştir.

Yazar katkıları:

• Yunkai Wu: Yöntem, biçimsel analiz, yazılım, araştırma ve veri düzenleme.
• Wangwang Kuang: Kavramsallaştırma, yöntem, ilk taslak, gözden geçirme ve düzenleme.
• Ling Xiao: Veri düzenleme, yazılım, görselleştirme ve doğrulama.
• Jin Dai: Biçimsel analiz, veri düzenleme ve doğrulama.
• Kui Gong: Kaynaklar, finansman edinimi ve proje yönetimi.
• Xushan Zhao: Danışmanlık ve finansman edinimi.
• Yin Wang: Danışmanlık, gözden geçirme ve düzenleme.

Veri erişimi: Verilerin talep üzerine paylaşılacağı belirtilmiştir. Açık veri deposu bağlantısı verilmemiştir.

Kod erişimi: CUDA C++ çözücüsünün veya Advanced Battery Designer yazılımının açık kaynak deposu belirtilmemiştir.

Çıkar çatışması: Yüklenen sürümde açık bir çıkar çatışması beyanı yer almamaktadır.

Etik kurul: Çalışma insan veya hayvan katılımcı içermeyen sayısal mühendislik araştırmasıdır. Etik kurul bilgisi verilmemiştir.

Bu Türkçe açıklama, yüklenen 34 sayfalık çalışmanın matematiksel modeli, sayısal yöntemleri, GPU uygulaması, tabloları ve Şekil 1–8’deki sonuçları temel alınarak hazırlanmıştır. Çalışmada bulunmayan fiziksel hücre deneyi, güvenlik sonucu, çevrim ömrü, üretim maliyeti veya ticari başarı iddiası eklenmemiştir. Dış kaynaklar yalnızca başlık, DOI, SSRN kayıt tarihi, yayın durumu ve kurumsal kimliğin bibliyografik olarak doğrulanması amacıyla kullanılmıştır.

Çalışmanın “yaklaşık 200 kat hızlanma” sonucu, bütün çok fazlı elektrokimyasal sistemin baştan sona süresi olarak yorumlanmamalıdır. Tablo 3’te karşılaştırılan bölüm, elektrolit potansiyeli için oluşturulan doğrusal denklem sistemidir. Özel çözücü 64³ ve 128³ ağlarda sırasıyla 169 ve 248 kat hızlanma sağlamıştır.

134 milyon düğüm için bildirilen yaklaşık 1,3 saniyelik süre de tek bir potansiyel çözüm adımına aittir. Tam deşarj simülasyonu, lityum yoğunluğu ve tepkime hesaplarıyla birlikte çok sayıda zaman adımı gerektirir.

CBD içeren modelin deneysel gerilim eğrilerine yakınlığı önemli olmakla birlikte bazı kinetik parametreler aynı deneysel verilere göre kalibre edilmiştir. Modelin bağımsız bir elektrot, farklı bir üretim partisi veya başka bir kimya üzerinde dış doğrulaması yapılmamıştır.

Delme tasarımındaki yüzde 20–117 kapasite artışları yalnızca seçilen sayısal geometri, malzeme parametreleri, kesme gerilimi ve akım yoğunlukları için geçerlidir. Kanalların kaldırdığı etkin malzeme, mekanik bütünlük, hacimsel enerji yoğunluğu ve üretim uygulanabilirliği birlikte değerlendirilmeden bu oranlar gerçek hücre performansına aktarılamaz.

Araştırmacıların GPU hızlandırmalı SBM’nin bu alandaki ilk uygulamalardan biri olduğu yönündeki yenilik iddiası, çalışmanın kendi literatür incelemesine dayanmaktadır. Bu açıklama kapsamında bağımsız ve eksiksiz bir öncelik veya patent araştırması yapılmamıştır.