محاسبه‌ تقریبی ضریب شدت تنش ترک‌های سطحی در ذرات کروی الکترودها در اثر جدایش فازی

آزمایش‌ها نشان داده‌اند که جدایش فازی در الکترود باتری‌های یون لیتیومی می‌تواند باعث تشکیل خرابی‌های مکانیکی و کاهش ظرفیت در این باتری‌ها شود. هدف این تحقیق، محاسبه‌ عددی و تحلیلی ضریب شدت تنش ترک‌های سطحی در الکترودهای متشکل از ذرات کروی ‌شکل در حضور جدایش فازی طی فرایند الکتروشیمیایی دشارژ است. به این منظور، توزیع غلظت و تنش‌های ناشی از نفوذ در یک ذره‌ کروی با حل عددی معادلات حاصل بر اساس مدل میدان فازی به‌دست می‌آید. به کمک یک مدل هندسی که در مطالعات پیشین پیشنهاد شده، ضریب شدت تنش برای عمیق‌ترین نقطه از ترک با هندسه‌ نیم‌بیضوی به کمک روش تابع وزن به‌دست می‌آید. به علاوه، به کمک مدل هسته پوسته با مرز فازی تیز و بهره‌گیری از یک روش تحلیلی، ضریب شدت تنش حداکثر که در عمیق‌ترین نقطه‌ ترک‌های سطحی طی فرایند دشارژ کامل حادث می‌شود، تخمین زده می‌شود. نتایج عددیِ توزیع غلظت و توزیع تنش حلقوی ارائه، و همچنین ضرایب شدت تنش به‌دست آمده بر اساس مدل میدان فازی با نتایج مدل هسته پوسته مقایسه می‌شود. نتایج حاصل از مدل عددی نشان می‌دهد که با کاهش اختلاف غلظت در دو فاز، ضرایب شدت تنش حداکثر، وابستگی بیشتری نسبت به تغییر نرخ جریان اعمالی بر سطح ذره داشته و می‌تواند با افزایش نرخ جریان اعمالی در محدوده پارامترهای مورد بررسی در این مطالعه تا بیش از دو برابر افزایش یابد.

Estimation of the Stress Intensity Factors for Surface Cracks in Spherical Electrode Particles Subject to Phase Separation

Experiments have frequently shown that phase separation in lithiumion battery electrodes could lead to the formation of mechanical defects, hence causing capacity fading. The purpose of the present work has been to examine stress intensity factors for preexisting surface cracks in spherical electrode particles during electrochemical deintercalation cycling using both analytical and numerical methods. To this end, we make use of a phase field model to examine the timedependent evolution of the concentration and stress profiles in a phase separating spherical electrode particles. By using a geometrical approximation scheme proposed in the literature, stress intensity factors at the deepest point of the preexisting surface cracks of semielliptical geometry are calculated with the aid of the wellestablished weight function method of fracture mechanics. By taking advantage of a sharpinterphase coreshell model, an analytical solution for the maximum stress intensity factors arising at the deepest point of the surface cracks during a complete deintercalation halfcycle is also developed. Numerical results for evolution of the concentration profile and the distribution of the hoop stresses in the particle are presented; further, the stress intensity factors found numerically based on the phase field model are compared with those predicted by the analytical coreshell model. The results of the numerical model suggest that the maximum stress intensity factor could significantly vary with changes in the surface flux, increasing potentially by a factor of two within the range of parameters considered here, when the concentration difference between the two phases is decreased.