Büyük silindirik lityum iyon bataryaların aşırı soğukta çalıştırılmadan veya şarj edilmeden önce ısıtılması gerekir. Ancak bataryanın yalnızca ortalama sıcaklığını yükseltmek yeterli değildir. Hücre içinde sıcak ve soğuk bölgeler oluşması; eşit olmayan elektrokimyasal davranışa, yerel yaşlanmaya, ısıl gerilmelere ve güvenlik sorunlarına yol açabilir.

Bu çalışmada araştırmacılar, silindirik bir 4680 lityum iyon bataryanın iki iç nikel ısıtıcı ve dış kenara yerleştirilmiş bir üçüncü nikel ısıtıcıyla ön ısıtılmasını sayısal olarak incelemiştir. Isıtıcıların konumları doğrudan milimetreyle değil, aralarında bulunan jelly-roll katmanlarının sayısıyla tanımlanmıştır. Böylece ısıtıcı yerleşimi, hücrenin katmanlı üretim yapısıyla ilişkilendirilen ayrık bir tasarım problemine dönüştürülmüştür.

Birinci ve ikinci katman sayıları 8 ile 24 arasında değiştirilmiş ve toplam 289 düzen COMSOL tabanlı parametrik simülasyonlarla karşılaştırılmıştır. Bütün düzenlerde toplam ısıtma gücü 522,5 W, etkin ısıtma süresi 60 saniye ve enerji girdisi 31,35 kJ olarak sabit tutulmuştur. Isıtıcılar kapatıldıktan sonra sıcaklığın hücre içinde yeniden dağılması için 30 saniyelik dinlenme dönemi modellenmiştir.

Referans yerleşimde üç bölgenin her birinde 17 jelly-roll katmanı bulunurken optimum olarak bildirilen yerleşim 14, 19 ve 18 katmandan oluşmuştur. Bu düzenin 90 saniyelik sürecin sonundaki maksimum sıcaklığı 24,82 °C, minimum sıcaklığı 19,16 °C ve hücre içi sıcaklık farkı 5,662 °C’dir. Referans çalışmada kullanılan 12,2 °C değerine göre sıcaklık düzensizliği %53,6 azalmıştır.

Bununla birlikte bu sonuç, fiziksel bir 4680 hücre üzerinde doğrulanmamıştır. Model sıcaklığa bağlı malzeme özelliklerini, katmanlar arası temas direncini, elektrokimyasal ısı üretimini, şarj durumunu ve lityum kaplanmasını içermemektedir. Doğrulama da gerçek deney verisine karşı değil, daha önce yayımlanmış başka bir sayısal modelin sonuçlarına karşı yapılmıştır.

Ayrıca 15–35 °C koşulu yalnızca 30 saniyelik dinlenme döneminin sonuna uygulanmıştır. Optimum düzenin sıcaklık-zaman grafiğinde maksimum hücre sıcaklığı, aktif ısıtmanın sona erdiği 60. saniyede yaklaşık 41–42 °C’ye çıkmakta ve daha sonra 24,8 °C’ye düşmektedir. Dolayısıyla hücrenin bütün ön ısıtma süresi boyunca 35 °C’nin altında kaldığı söylenemez.

Aşırı soğuk lityum iyon bataryaları neden zorlar?

Lityum iyon bataryaların performansı sıcaklığa güçlü biçimde bağlıdır. Düşük sıcaklıkta elektrolitin viskozitesi artar, lityum iyonlarının taşınması yavaşlar, elektrotlardaki difüzyon sınırlanır ve yük aktarım direnci yükselir. Bunun sonucunda kullanılabilir kapasite ve güç azalabilir, şarj süresi uzayabilir ve bataryanın gerilim tepkisi zayıflayabilir.

Düşük sıcaklıkta şarj sırasında önem kazanan başka bir risk, lityum iyonlarının grafit anot içine yeterince hızlı yerleşememesi ve anot yüzeyinde metalik lityum olarak birikmesidir. Lityum kaplanması olarak adlandırılan bu süreç geri döndürülemez kapasite kaybına yol açabilir ve bazı koşullarda iç kısa devre olasılığını artırabilir.

Çalışmada lityum iyon bataryalar için istenen çalışma sıcaklığı 15–35 °C aralığı olarak kabul edilmiştir. Araştırmanın hedefi, hücreyi yalnızca bu aralığa taşımak değil, aynı zamanda hücre içindeki maksimum ve minimum sıcaklık arasındaki farkı küçültmektir.

Dıştan ve içten ısıtma arasındaki fark nedir?

Dış ısıtmada ısı; hücre kabuğuna temas eden plaka, hava, sıvı veya başka bir termal sistem üzerinden bataryaya aktarılır. Bu yaklaşım üretim açısından görece basit olabilir. Ancak ısı, büyük silindirik bir hücrenin dış yüzeyinden merkezine iletimle ilerlediği için merkez ve dış bölgeler arasında yüksek sıcaklık farkları oluşabilir.

İç ısıtmada ise ince metalik elemanlar bataryanın katmanlı yapısına yerleştirilir. Elektrik akımı bu elemanlardan geçtiğinde Joule ısısı oluşur. Isı doğrudan aktif katmanlara yakın bölgelerde üretildiği için ısı aktarım yolu kısalabilir ve soğuk başlangıç süresi azaltılabilir.

İç ısıtmanın kendi tasarım sorunları vardır. Isıtıcılar birbirine çok yakın yerleştirilirse termal alanlar aşırı örtüşebilir ve sıcak noktalar oluşabilir. Çok uzak yerleştirilirse aralarında yeterince ısıtılmayan bölgeler kalabilir. Isıtıcı sayısı ve gücü değiştirilmese bile yalnızca konumlarının değiştirilmesi, hücre içindeki sıcaklık dağılımını önemli ölçüde etkileyebilir.

Çalışmanın araştırma boşluğu nedir?

Önceki çalışmaların çoğu şu değişkenlere odaklanmıştır:

• İç ısıtıcıların sayısı,
• Uygulanan elektrik gücü,
• Ortam sıcaklığı,
• Taşınımla ısı aktarım katsayısı,
• Isıtma ve bekleme süreleri,
• Darbeli veya değişken görev çevrimli ısıtma stratejileri.

Bu çalışmanın ayırıcı yönü, iki iç ısıtıcının radyal konumlarını hücredeki jelly-roll katmanlarının sayısıyla tanımlaması ve bu konumlar için geniş bir ayrık tarama gerçekleştirmesidir.

Araştırmanın temel sorusu şöyledir:

Toplam güç, ısıtma süresi, dış boyutlar ve sınır koşulları değiştirilmeden, ısıtıcıların jelly-roll içindeki yerleri değiştirilerek hücre sıcaklığı daha eşit dağıtılabilir mi?

4680 hücrenin katmanlı yapısı nasıl modellenmiştir?

Hücrenin aktif bölgesi, tekrar eden sekiz malzeme katmanından oluşan jelly-roll birimleriyle temsil edilmiştir:

Bu sekiz katmanın toplam kalınlığı 424 µm’dir. Üç temel jelly-roll bölgesinde toplam 51 tekrar birimi bulunmaktadır.

Şekil 1’de hücre şu bölgelere ayrılmıştır:

• Polimer merkez ile birinci iç ısıtıcı arasındaki Pack 1,
• Birinci ve ikinci iç ısıtıcı arasındaki Pack 2,
• İkinci iç ısıtıcı ile dış ısıtıcı arasındaki Pack 3.

Şeklin kendi açıklama kutusunda model “2D axisymmetric model” olarak tanımlanmıştır. Başka bir ifadeyle geometri, hücrenin ekseni etrafındaki simetri kullanılarak iki boyutlu bir radyal-eksenel kesit üzerinden çözülmektedir. Buna karşın yöntem metninin ilerleyen bölümünde “üç boyutlu sayısal model” ifadesi kullanılmıştır. Bu iki tanım birbiriyle uyumlu değildir.

Isıtıcılar nasıl modellenmiştir?

İç ısıtıcılar 30 µm kalınlığında nikel folyolar olarak modellenmiştir. Her nikel katmanın iki tarafında, elektriksel teması önlemek için 28 µm kalınlığında PET yalıtım kaplaması bulunmaktadır. Isının jelly-roll yapısına aktarılmasını desteklemek amacıyla 9 µm kalınlığında bir bakır ısı dağıtım katmanı da modele eklenmiştir.

Birinci iç ısıtıcıya 118,9 W, ikinci iç ısıtıcıya 224,2 W ve dış ısıtıcıya 179,4 W güç uygulanmıştır. Dış ısıtıcının gücü, dış kabuğa yakın bölgelerde aşırı ısı birikimini azaltmak amacıyla ikinci iç ısıtıcının gücünün %80’i olarak belirlenmiştir.

Katman tabanlı tasarım değişkenleri nasıl tanımlanmıştır?

Üç bölgedeki jelly-roll katmanı sayıları \(n_1\), \(n_2\) ve \(n_3\) ile gösterilmiştir:

• \(n_1\): Polimer merkez ile birinci iç ısıtıcı arasındaki katman sayısı,
• \(n_2\): İki iç ısıtıcı arasındaki katman sayısı,
• \(n_3\): İkinci iç ısıtıcı ile dış ısıtıcı arasındaki katman sayısı.

Toplam aktif katman sayısı sabit tutulmuştur:

\[ n_1+n_2+n_3=N_{JR} \]

\[ N_{JR}=51 \]

Bu nedenle üçüncü bölgenin katman sayısı şu şekilde hesaplanmıştır:

\[ n_3=51-n_1-n_2 \]

\(n_1\) ve \(n_2\), 8 ile 24 arasında 17 farklı tamsayı değer almıştır. Böylece:

\[ 17\times17=289 \]

farklı ısıtıcı yerleşimi oluşturulmuştur.

Arama uzayının küçük bir asimetrisi vardır. \(n_1\) ve \(n_2\) 8–24 aralığıyla sınırlandırılırken \(n_3\) için aynı sınır konulmamıştır. Denklem gereği üçüncü bölge teorik olarak 3 ile 35 katman arasında değişebilmektedir. En iyi sonuçların bulunduğu bölgede \(n_3\) 18 veya 19’dur; ancak bütün 289 konfigürasyonun aynı üretim kısıtları altında değerlendirilip değerlendirilmediği ayrıca açıklanmamıştır.

Isı transferi hangi denklemle çözülmüştür?

Her katı malzeme bölgesi için geçici ısı iletim denklemi kullanılmıştır:

\[ \rho_i C_{p,i}\frac{\partial T_i(\mathbf{x},t)}{\partial t} =\nabla\cdot[k_i\nabla T_i(\mathbf{x},t)]+\dot{q}_i \]

• \(\rho_i\): i numaralı malzemenin yoğunluğudur.
• \(C_{p,i}\): Özgül ısı kapasitesidir.
• \(T_i\): Yerel sıcaklıktır.
• \(k_i\): Isıl iletkenliktir.
• \(\dot{q}_i\): Birim hacimde üretilen ısıdır.

Isı üretimi yalnızca nikel ısıtıcılara uygulanmıştır:

\[ \dot{q}_{h,j}=\frac{P_{h,j}}{V_{h,j}} \]

Isıtıcı hacmi şu şekilde tanımlanmıştır:

\[ V_{h,j}=A_{h,j}\delta_h \]

• \(P_{h,j}\): j numaralı ısıtıcının gücüdür.
• \(A_{h,j}\): Etkin ısıtıcı alanıdır.
• \(\delta_h\): Nikel folyonun kalınlığıdır.

Dış yüzeyde taşınımla ısı kaybı uygulanmıştır:

\[ -k_i(\nabla T_i\cdot\mathbf{n}) =h(T_s-T_\infty) \]

Burada ortam sıcaklığı \(T_\infty=-50\ ^\circ C\), taşınım katsayısı ise \(h=5\ \text{W/(m²·K)}\) olarak verilmiştir.

Simetri ekseninde sıfır ısı akısı kabul edilmiştir:

\[ \nabla T\cdot\mathbf{n}=0 \]

Modelin temel varsayımları nelerdir?

• Yoğunluk, özgül ısı ve ısıl iletkenlik sıcaklıktan bağımsız ve sabittir.
• Katmanlar arasında kusursuz termal temas vardır.
• Katmanlar arası temas direnci sıfır kabul edilmiştir.
• Bataryanın elektrokimyasal reaksiyonları modellenmemiştir.
• Şarj durumu ve akım üretimi dikkate alınmamıştır.
• Lityum kaplanması kinetiği hesaplanmamıştır.
• Isıtıcı gücü hacimsel ısı kaynağı olarak uygulanmıştır.
• Malzeme deformasyonu ve termal gerilmeler çözülmemiştir.

Bu varsayımlar, yerleşimin saf geometrik etkisini incelemeyi kolaylaştırmaktadır. Ancak sonuçların gerçek hücreye aktarılmasını sınırlandırmaktadır.

Başlangıç sıcaklığı nedir?

Çalışmanın yöntem metninde hesaplama alanının 20 °C’lik düzgün bir başlangıç değeriyle başlatıldığı yazılmıştır. Tablo 4 de “geçici ısı aktarımı çözümü için başlangıç değeri”ni 20 °C olarak vermektedir. Buna karşılık Şekil 3, 4 ve 7–9’daki maksimum, ortalama ve minimum sıcaklık eğrilerinin tümü yaklaşık −50 °C’den başlamaktadır.

Bu iki bilgi aynı fiziksel başlangıç koşulunu temsil edemez. Olası açıklamalardan biri, 20 °C’nin çözücü için sayısal başlangıç tahmini; −50 °C’nin ise gerçek hücre başlangıç sıcaklığı olmasıdır. Ancak çalışma böyle bir ayrım yapmamaktadır. Bu nedenle simülasyonun fiziksel başlangıç sıcaklığı metinde tutarlı biçimde raporlanmış değildir.

Isıtma çevrimi nasıl uygulanmıştır?

İlk 60 saniyede ısıtıcılar aktiftir. Daha sonraki 30 saniyede güç kesilmekte ve ısının hücre içinde yeniden dağılması beklenmektedir. Optimizasyonda kullanılan son sıcaklık değerleri, bu 90 saniyelik sürecin sonundaki değerlerdir.

Sıcaklık eşitliği nasıl ölçülmüştür?

Hücrenin hacim ortalamalı sıcaklığı şu şekilde tanımlanmıştır:

\[ \overline{T}(t) =\frac{1}{V_{cell}} \int_{\Omega_{cell}}T(\mathbf{x},t)\,dV \]

Aynı hesaplama alanındaki en yüksek ve en düşük sıcaklıklar:

\[ T_{max}(t)=\max_{\mathbf{x}\in\Omega_{cell}}T(\mathbf{x},t) \]

\[ T_{min}(t)=\min_{\mathbf{x}\in\Omega_{cell}}T(\mathbf{x},t) \]

olarak belirlenmiştir. Ana termal eşitlik ölçütü ise:

\[ \Delta T(t)=T_{max}(t)-T_{min}(t) \]

biçimindedir.

Küçük \(\Delta T\), hücre içindeki sıcaklıkların birbirine daha yakın olduğunu gösterir. Ancak yalnızca sıcaklık farkını küçültmek yeterli değildir. Örneğin bütün hücrenin −20 ile −15 °C arasında bulunması yalnızca 5 °C fark oluşturur; buna rağmen hücre yeterince ısıtılmış olmaz.

Optimizasyon problemi nasıl yazılmıştır?

Çalışmanın matematiksel amacı, 90 saniyenin sonundaki sıcaklık farkını en aza indirmektir:

\[ \min_{n_1,n_2}\Delta T(t_f) \]

Kısıtlar şöyledir:

\[ n_1,n_2\in\{8,9,\ldots,24\} \]

\[ n_3=51-n_1-n_2 \]

\[ 15\ ^\circ C\leq T_{min}(t_f) \leq T_{max}(t_f)\leq35\ ^\circ C \]

Bu kısıtlar yalnızca son zaman olan \(t_f\) için geçerlidir. Aktif ısıtma sırasında oluşan en yüksek sıcaklık için ayrı bir üst sınır yazılmamıştır.

Ağ bağımsızlığı ne kadar güçlüdür?

Şekil 2’de yaklaşık 600 bin, 760 bin ve 1,3 milyon elemanlı ağlar karşılaştırılmıştır. Eleman sayısı artırıldığında:

• Maksimum sıcaklık yaklaşık 27,48 °C’den 26,36 °C’ye,
• Ortalama sıcaklık yaklaşık 24,30 °C’den 23,41 °C’ye,
• Minimum sıcaklık yaklaşık 16,69 °C’den 14,33 °C’ye

değişmektedir.

Şekilde minimum sıcaklığın değişimi yaklaşık −%14,5 olarak gösterilmektedir. Bu büyüklük, sıcaklık farkının optimize edildiği bir çalışmada önemsiz kabul edilebilecek kadar küçük değildir. Ayrıca sonraki analizlerde hangi ağ yoğunluğunun seçildiği ve seçim ölçütünün ne olduğu açıkça belirtilmemiştir.

Bu nedenle çalışma bir ağ duyarlılığı incelemesi sunmaktadır; ancak “sonuçlar ağ boyutundan bağımsızdır” sonucu özellikle minimum sıcaklık bakımından yeterince güçlü değildir.

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

Araştırmacılar eşit aralıklı referans düzeni kullanmıştır:

\[ n_1=n_2=n_3=17 \]

Bu geometri, önceki bir 4680 iç ısıtma çalışmasının sınır koşulları ve güçleriyle çözülmüş; mevcut modelin sıcaklık eğrileri önceki sayısal sonuçlarla karşılaştırılmıştır.

Maksimum ve minimum değerler arasındaki farkın %1’den küçük olduğu belirtilmiştir. Sıcaklık konturlarının genel biçimi de önceki simülasyona benzemektedir.

Bu işlem deneysel doğrulama değildir. Mevcut model, başka bir sayısal modelin sonuçlarını yeniden üretmektedir. Dolayısıyla iki modelin aynı kabulleri veya aynı eksiklikleri paylaşması mümkündür.

289 yerleşim ne göstermiştir?

Şekil 5’te \(n_1\) ve \(n_2\) değerlerinin her birleşimi için son sıcaklık farkı bir ısı haritasıyla gösterilmiştir. Haritadaki değerler yaklaşık 5,7 °C’den 133,6 °C’ye kadar uzanmaktadır.

Isıtıcıların çok dengesiz konumlandırıldığı bazı yerleşimlerde sıcaklık farkı aşırı yükselmiştir. Isıtıcılar belirli orta aralıklarda dağıtıldığında ise termal alanlar birbirini daha dengeli tamamlamış ve sıcaklık farkı küçülmüştür.

Şekil 6, sıcaklık farkı ile maksimum sıcaklık arasındaki ilişkiyi göstermektedir. Bazı konfigürasyonlar düşük \(\Delta T\) üretse de sıcaklıkları istenen 15–35 °C aralığında değildir. Bu nedenle araştırmacılar yalnızca ısı haritasının en küçük değerini değil, son maksimum ve minimum sıcaklıkları da değerlendirmiştir.

Optimum olarak hangi yerleşim seçilmiştir?

Çalışmada optimum olarak bildirilen yapı şöyledir:

Bu yerleşimin son sıcaklıkları şöyledir:

• Maksimum sıcaklık: 24,82 °C,
• Minimum sıcaklık: 19,16 °C,
• Sıcaklık farkı: 5,662 °C.

Şekil 7(a), dış ve iç ısıtıcıların termal alanlarının hücre boyunca daha dengeli dağıldığını göstermektedir. En sıcak bölgeler ısıtıcılara yakın iç alanlarda bulunurken dış ve merkez bölgeler arasındaki fark referans yerleşime göre azalmıştır.

Optimum çözüm benzersiz midir?

Hayır. Tablo 5’te şu iki yerleşim aynı sıcaklık farkını üretmektedir:

Her iki yerleşim de matematiksel kısıtlardaki 15–35 °C aralığını sağlamaktadır. Dolayısıyla Denklem 12 ve 13’te yazılan amaç fonksiyonu açısından iki düzen de eşit optimumdur.

Araştırmacılar 14–19–18 düzenini, aynı güç altında biraz daha yüksek minimum ve maksimum sıcaklık üretmesi ve sıcaklık dağılımını daha uygun görmeleri nedeniyle tercih etmiştir. Bu, anlaşılabilir bir mühendislik seçimi olabilir; ancak “daha güçlü genel ısıtma” veya belirli bir kontur niteliği matematiksel amaç fonksiyonuna ikincil ölçüt olarak eklenmemiştir.

Bu nedenle çalışma tek bir matematiksel optimumdan çok, eşit sıcaklık farkına sahip iki aday arasından ek bir nitel değerlendirmeyle seçilmiş tercih edilen düzen sunmaktadır.

Diğer düşük sıcaklık farkına sahip düzenler nasıldır?

Şekil 9’daki 12–21–18 düzeninde sıcaklık farkı yalnızca 6,06 °C olmasına rağmen ısıtıcı çevrelerinde daha belirgin sıcaklık yığılması görülmektedir. Bu karşılaştırma, aynı skaler sıcaklık farkına yakın iki alanın uzamsal dağılımlarının farklı olabileceğini göstermektedir.

Hücre gerçekten 15–35 °C aralığında mı kalmıştır?

Yalnızca son durumda. Optimum yerleşimin Şekil 7(b)’deki sıcaklık eğrisinde maksimum sıcaklık, aktif ısıtmanın bittiği 60. saniyede yaklaşık 41–42 °C’ye ulaşmaktadır. Isıtıcılar kapatıldıktan sonra sıcak bölgedeki ısı daha soğuk bölgelere dağılmakta ve maksimum sıcaklık 90. saniyede 24,82 °C’ye inmektedir.

Bu nedenle çalışmanın 15–35 °C kısıtı:

• Isıtmanın her anı için uygulanmamıştır.
• Isıtıcı yüzeyindeki geçici tepe sıcaklığı sınırlanmamıştır.
• Yalnızca 30 saniyelik dinlenme sonrasındaki son sıcaklıklara uygulanmıştır.

Gerçek bir hücrede yalıtım tabakalarının, elektrolitin veya yerel elektrot bölgelerinin kısa süreli yüksek sıcaklıklara nasıl tepki vereceği ayrıca değerlendirilmelidir.

%53,6 iyileşme nasıl hesaplanmıştır?

Referans çalışmanın sıcaklık farkı 12,2 °C, optimum düzenin sıcaklık farkı 5,662 °C olarak alınmıştır:

\[ \frac{12.2-5.662}{12.2}\times100 \approx53.6\% \]

Ancak mevcut modelin kendi referans simülasyonunda maksimum sıcaklık 26,9 °C ve minimum sıcaklık 14,9 °C olduğundan fark 12,0 °C’dir. Optimum değer bu 12,0 °C ile karşılaştırılırsa azalma yaklaşık %52,8 olur.

Aradaki fark küçüktür; ancak çalışmanın %53,6 değeri kendi referans simülasyonuna değil, önceki çalışmanın 12,2 °C değerine dayanmaktadır.

Enerji değerlendirmesi ne göstermektedir?

Bütün yerleşimlerde toplam güç ve ısıtma süresi aynıdır:

\[ E_{in}=P_{total}t_{heating} \]

\[ E_{in}=522.5\ \text{W}\times60\ \text{s}=31.35\ \text{kJ} \]

Çalışmada ayrıca enerjiye göre normalize edilmiş bir sıcaklık eşitliği geliştirme göstergesi tanımlanmıştır:

\[ \eta_{\Delta T} =\frac{\Delta T_{ref}-\Delta T_{opt}}{E_{in}} \]

\[ \eta_{\Delta T} =\frac{12.2-5.662}{31.35} =0.208\ ^\circ\text{C/kJ} \]

Bu değer klasik anlamda enerji verimliliği değildir. Hücrenin aldığı yararlı enerjiyi, elektriksel kayıpları veya ortam kaybını hesaplamaz. Yalnızca aynı enerji verilen iki yerleşim arasındaki sıcaklık farkı iyileşmesini birim enerjiyle ifade eden, çalışmaya özgü karşılaştırma göstergesidir.

Isıtıcı yerleşiminin üretim açısından anlamı nedir?

Isıtıcı konumlarının katman sayısıyla tanımlanması, tasarımın jelly-roll sarım sürecine bağlanması açısından yararlıdır. Optimum düzen:

• Isıtıcı sayısını artırmamaktadır.
• Toplam ısıtma gücünü yükseltmemektedir.
• Hücrenin dış boyutlarını değiştirmemektedir.
• Isı üretim noktalarını aktif yapının içinde yeniden dağıtmaktadır.

Bununla birlikte çalışma bazı bölümlerde iyileştirmenin “hücrede yapısal değişiklik yapılmadan” elde edildiğini söylemektedir. İç ısıtıcıların jelly-roll içindeki yerinin değiştirilmesi, dış boyutları değiştirmese de hücrenin iç katman düzeninde bir tasarım değişikliğidir. Daha doğru ifade, “dış boyutlar, ısıtıcı sayısı ve toplam güç değiştirilmeden” olmalıdır.

Pratik uygulamada hangi güvenlik sorunları vardır?

• Nikel folyonun aktif katmanlardan güvenilir biçimde yalıtılması,
• PET kaplamanın sarım, titreşim ve çevrim yaşlanmasına dayanması,
• Isıtıcı veya bağlantı arızasının iç kısa devre oluşturmaması,
• Yerel sıcak noktaların algılanması ve gücün zamanında kesilmesi,
• Isıtıcının hücre kapasitesi ve enerji yoğunluğundan ne kadar hacim aldığı,
• Ek metal ve polimerlerin geri dönüşüm sürecine etkisi,
• Isıtıcının termal kaçak sırasında nasıl davrandığı,
• Üretim toleranslarının optimum katman aralıklarını nasıl değiştirdiği

çalışmanın sayısal optimizasyonunda çözülmemiştir.

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

• Isıtıcıların radyal konumu, aynı güç altında hücre içi sıcaklık dağılımını önemli ölçüde değiştirebilir.
• Katman sayısı, silindirik hücrelerde üretimle ilişkili bir ayrık yerleşim değişkeni olarak kullanılabilir.
• Isıtıcıların aşırı yakın veya aşırı uzak olması sıcak noktaları ya da soğuk bölgeleri artırabilir.
• İncelenen model ve sınır koşullarında 14–19–18 ile 13–20–18 düzenleri en düşük son sıcaklık farkını vermiştir.
• 14–19–18 düzeni, referans sayısal sonuca göre sıcaklık farkını yaklaşık yarı yarıya azaltmıştır.
• Yerleşim optimizasyonuyla ek elektrik enerjisi verilmeden daha eşit son sıcaklık alanı elde edilebilir.

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

• 14–19–18 yerleşiminin bütün 4680 hücre kimyaları için optimum olduğunu kanıtlamamaktadır.
• Gerçek bir bataryada %53,6 iyileşme elde edileceğini göstermemektedir.
• Hücrenin bütün ısıtma süresince 15–35 °C arasında kaldığını göstermemektedir.
• İç ısıtıcıların çevrim ömrünü uzattığını deneysel olarak kanıtlamamaktadır.
• Lityum kaplanmasının önlendiğini göstermemektedir.
• Isıtıcıların elektriksel yalıtım ve kısa devre güvenliğini doğrulamamaktadır.
• Modelin ağ boyutundan tamamen bağımsız olduğunu güçlü biçimde göstermemektedir.
• Önerilen yerleşimin seri üretimde sarılabilir ve güvenilir olduğunu kanıtlamamaktadır.
• İlave nikel, bakır ve yalıtım katmanlarının enerji yoğunluğuna etkisini hesaplamamaktadır.
• Diğer ortam sıcaklıkları, güç seviyeleri ve taşınım koşullarında aynı optimumun korunacağını göstermemektedir.

Türkiye’deki elektrikli araç ve batarya çalışmaları açısından anlamı nedir?

Araştırma belirli bir ülkenin kullanım verisine bağlı değildir. Büyük silindirik bataryaların soğuk ortamda ön ısıtılmasıyla ilgili genel bir termal tasarım problemi ele alınmaktadır.

Türkiye’de benzer bir tasarımın değerlendirilmesi için şu çalışmalar önemlidir:

• Farklı hücre kimyalarıyla gerçek 4680 numunelerinin üretilmesi,
• −10, −20, −30 ve daha düşük sıcaklıklarda deneysel doğrulama,
• Isıtıcıların gerçek jelly-roll sarımına mekanik olarak eklenmesi,
• Yalıtım delinmesi ve iç kısa devre testleri,
• Isıtma sırasında hücrenin farklı noktalarında termokupl ölçümleri,
• Aktif ısıtma sırasındaki geçici sıcaklık tepesinin sınırlandırılması,
• Şarj başlangıcıyla lityum kaplanması riskinin birlikte modellenmesi,
• Yaşlandırma ve çevrim ömrü deneyleri,
• Batarya paketindeki komşu hücreler ve soğutma sistemiyle etkileşimin incelenmesi,
• Ek ısıtıcı malzemelerinin kapasite, maliyet ve geri dönüşüme etkisinin hesaplanması

gereklidir. Çalışmadaki katman aralıkları, deneysel doğrulama olmadan doğrudan bir üretim reçetesi veya güvenlik standardı olarak kullanılmamalıdır.

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

Teknik yöntem özeti

Referans ve optimize edilmiş yerleşimin karşılaştırması

Raporlanan temel kazanımlar

Metin ve görsellerdeki kritik tutarsızlıklar

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

• Isıtıcı konumunu doğrudan jelly-roll katman sayısıyla ilişkilendirmiştir.
• 289 farklı ayrık yerleşimi aynı güç ve sınır koşulları altında karşılaştırmıştır.
• Maksimum, minimum, ortalama sıcaklık ve sıcaklık farkını birlikte değerlendirmiştir.
• Yalnızca ortalama sıcaklığı değil, uzamsal sıcaklık konturlarını da incelemiştir.
• Isıtma ve ısı dağılımı için ayrı aktif ve dinlenme dönemleri modellemiştir.
• Malzeme katmanlarını tek bir eşdeğer ortam yerine ayrı termofiziksel özelliklerle temsil etmiştir.
• Önerilen düzenin dış boyut veya toplam enerji artışı gerektirmediğini göstermiştir.
• Üretim, yalıtım, güvenlik ve geri dönüşüm sorunlarını sonuç tartışmasına dâhil etmiştir.
• Finansman, veri erişimi, çıkar çatışması ve yapay zekâ kullanım beyanları verilmiştir.

Çalışmanın temel sınırlılıkları

• Fiziksel hücre veya prototip üzerinde deney yapılmamıştır.
• Doğrulama bağımsız deney verisi yerine başka bir sayısal modele dayanmaktadır.
• Başlangıç sıcaklığı ve model boyutu tutarlı raporlanmamıştır.
• Ağ bağımsızlığı özellikle minimum sıcaklık bakımından ikna edici değildir.
• Malzeme özellikleri sıcaklıktan bağımsız kabul edilmiştir.
• Katmanlar arası termal temas direnci ihmal edilmiştir.
• Elektrokimyasal ısı üretimi ve şarj durumu bulunmamaktadır.
• Lityum kaplanması yalnızca motivasyon olarak anılmış, modellenmemiştir.
• Isıtıcıların aktif olduğu dönemde 35 °C üzerindeki geçici sıcaklık sınırlandırılmamıştır.
• Yalıtım delinmesi, iç kısa devre ve termal kaçak değerlendirilmemiştir.
• Isıtıcıların aktif malzeme hacmi ve enerji yoğunluğu üzerindeki etkisi hesaplanmamıştır.
• Optimum yalnızca tek ortam sıcaklığı ve tek güç dağılımı için belirlenmiştir.
• \(n_3\) değişkenine \(n_1\) ve \(n_2\) ile aynı alt ve üst sınırlar uygulanmamıştır.
• İki eşit amaç değerinden birinin seçilmesini sağlayan ikincil ölçüt matematiksel olarak tanımlanmamıştır.
• Simülasyon dosyaları ve tam veri seti açık bir depoda yayımlanmamıştır.

Sonuçları güçlendirecek ek çalışmalar

• Fiziksel 4680 prototip üzerinde termokupl ve termal kamera doğrulaması,
• İç ısıtıcıların jelly-roll üretimine gerçek olarak eklenmesi,
• Sıcaklığa bağlı termofiziksel özelliklerin kullanılması,
• Katmanlar arası temas dirençlerinin ölçülmesi ve modele eklenmesi,
• Termal-elektrokimyasal eşlenik model geliştirilmesi,
• Isıtma sırasında lityum kaplanması riskinin hesaplanması,
• Aktif ısıtma tepe sıcaklığına ayrı bir kısıt konulması,
• Ortam sıcaklığı ve güç için dayanıklılık optimizasyonu yapılması,
• Isıtıcı konumu ile güç dağılımının birlikte optimize edilmesi,
• Farklı hücre kimyaları ve yaşlanma düzeylerinin sınanması,
• Gerçek batarya paketi içinde komşu hücrelerin etkisinin incelenmesi,
• Yalıtım ve iç kısa devre güvenlik testleri,
• Enerji yoğunluğu, üretim maliyeti ve geri dönüşüm etkisinin hesaplanması,
• Daha güçlü ağ yakınsama çalışması ve optimumun farklı ağlarda yeniden hesaplanması,
• Çok ölçütlü optimizasyonda \(\Delta T\), ortalama sıcaklık, tepe sıcaklığı ve enerji tüketiminin birlikte kullanılması.

Kaynak ve Yöntem Notu

• Çalışmanın özgün adı: Discrete Optimization of Internal Heater Arrangement Based on Thermal Uniformity for Hybrid Preheating of 4680 Lithium-Ion Batteries in Extremely Cold Environments
• Yazarlar: Hadi Fallah Ardeshir; Alireza Alizadeh Milani
• Yazar sıralaması: Hadi Fallah Ardeshir birinci, Alireza Alizadeh Milani ikinci yazardır.
• Eş katkı veya eş birinci yazarlık: Beyan bulunmamaktadır.
• Sorumlu yazar: Hadi Fallah Ardeshir
• Sorumlu yazar e-postası: fallah.hadi@sut.ac.ir
• Kurum: Faculty of Mechanical Engineering, Sahand University of Technology, Tabriz, Iran
• Çalışma yılı: 2026
• Kesin SSRN yükleme tarihi: Erişilebilen çalışma metninde açıkça belirtilmemiş ve resmî kayıttan kesin gün olarak doğrulanamamıştır.
• Kaynak türü: Hakem değerlendirmesinden geçmemiş preprint ve sayısal modelleme çalışması
• Yayın platformu: SSRN
• Preprint DOI: 10.2139/ssrn.7196306
• Resmî kayıt:SSRN resmî çalışma sayfası
• Kalıcı DOI bağlantısı:10.2139/ssrn.7196306
• Dergi: Hakemli bir dergide yayımlandığı doğrulanmamıştır.
• Özgün yayınevi: Hakemli dergi yayınevi bulunmamaktadır; çalışma SSRN preprint platformunda sunulmuştur.
• Hakemlik durumu: Çalışmanın her sayfasında hakem değerlendirmesinden geçmediği belirtilmektedir.
• Finansman: Araştırmacılar kamu, ticari veya kâr amacı gütmeyen bir kuruluştan özel finansman alınmadığını bildirmiştir.
• Çıkar çatışması: Bilinen finansal çıkar çatışması veya kişisel ilişki bulunmadığı beyan edilmiştir.
• Etik onay: İnsan veya hayvan çalışması içermediği için uygulanabilir olmadığı belirtilmiştir.
• Veri erişimi: Sayısal verilerin sorumlu yazardan makul talep üzerine edinilebileceği belirtilmiştir.
• Açık veri deposu: Tam simülasyon veri seti kamuya açık bir depoya yüklenmemiştir.
• Kod ve COMSOL dosyaları: Açık erişim bağlantısı verilmemiştir.
• Üretken yapay zekâ beyanı: ChatGPT’nin dil düzenleme, biçimlendirme ve okunabilirliği geliştirme amacıyla kullanıldığı; yazarların içeriği gözden geçirerek sorumluluğu üstlendiği belirtilmiştir.

Hadi Fallah Ardeshir’in CRediT katkıları; kavramsallaştırma, yöntem, yazılım, doğrulama, analiz, araştırma, kaynaklar, veri düzenleme, görselleştirme, ilk taslak, gözden geçirme, denetim ve proje yönetimini kapsamaktadır. Alireza Alizadeh Milani’nin katkıları yazılım, doğrulama, araştırma, veri düzenleme, görselleştirme ve yazım olarak bildirilmiştir.

Bu Verianla makalesinin bilimsel içeriği yalnızca yüklenen çalışmanın metni, denklemleri, malzeme tabloları, parametrik geometri şeması, ağ duyarlılığı grafiği, sıcaklık konturları, 289 düzenin ısı haritası, optimum karşılaştırmaları ve sonuç tartışması esas alınarak hazırlanmıştır. Dış kaynaklar yalnızca bibliyografik kimliğin, kurumun, DOI’nin ve resmî SSRN kaydının doğrulanması amacıyla kullanılmış; dış kaynaklardan yeni bilimsel bulgu eklenmemiştir.

Çalışmanın Şekil 1’inde model iki boyutlu eksenel simetrik olarak gösterilirken yöntem metninde üç boyutlu olarak adlandırılmıştır. Başlangıç sıcaklığı metinde 20 °C verilmiş, ancak grafiklerde sıcaklıklar yaklaşık −50 °C’den başlamıştır. Bu uyumsuzluklar sessizce düzeltilmemiştir.

Ağ duyarlılığı incelemesinde minimum sıcaklığın farklı ağlar arasında yaklaşık %14,5 değişmesi nedeniyle modelin ağ bağımsızlığı tam olarak gösterilmiş kabul edilmemelidir. Seçilen ağın eleman sayısı da açıkça belirtilmemiştir.

14–19–18 ve 13–20–18 düzenleri aynı 5,662 °C sıcaklık farkını üretmekte ve yazılan son sıcaklık kısıtlarını sağlamaktadır. Dolayısıyla bildirilen matematiksel optimizasyon problemi tek başına 14–19–18 düzenini benzersiz optimum olarak seçmemektedir. Araştırmacılar bu düzeni daha yüksek son sıcaklıklar ve görsel olarak daha uygun sıcaklık dağılımı gerekçesiyle tercih etmiştir.

%53,6 azalma, önceki referans çalışmanın 12,2 °C değerine göre hesaplanmıştır. Mevcut modelin kendi referans sonucu 12,0 °C olduğundan aynı hesaplama yaklaşık %52,8 verir. Ayrıca 15–35 °C kısıtı yalnızca 90. saniyedeki son duruma uygulanmış; aktif ısıtma sırasında görülen yaklaşık 41–42 °C’lik geçici maksimumu kapsamamıştır.

Çalışma hakem değerlendirmesinden geçmemiştir. Sonuçlar gerçek hücre deneyiyle doğrulanmadığından önerilen katman düzeni üretim reçetesi, güvenlik sınırı veya bütün 4680 hücreler için evrensel optimum olarak değerlendirilmemelidir.