리튬 이온 배터리 실리콘 나노입자 전극의 빠른 변형을 위한 개선된 화학 포텐셜 모델

요약

저변형률 가정을 제거한 정제된 화학 포텐셜 표현 제시, 실리콘 나노입자 전극 내 확산 유도 응력과 리튬 분포 정확 모델링 가능. 정전위 및 정전류 사이클링 조건에서 새 표현과 전통 표현 비교, 정전위 동작 시 현저한 차이 확인, 정전류 결과는 크게 영향 없음. 임계 반경(리튬 플럭스와 정수압 응력이 두 모델서 본질적으로 동일해지는 지점) 최초 규명.

서론

리튬화 시 부피 팽창 400% 실리콘 등 고용량 전극 소재는 휴대 전자기기, 전기차, 그리드 저장 핵심. 하지만 균질 부피 팽창 유발 확산 유도 응력 취성 파단 유발, 사이클 수명 단축, 용량 감소. 복합 전극 복잡 형상 순수 이론 기술 난제, 대부분 모델 물질 특성 공간적 변화 고려하나 계면 효과 무시. 이론 고려 전극 형상: 구형, 원통, 평판, 1-2차원 모델.

실리콘 나노입자 전극 최근 연구: 리튬화 유도 상변이, 형태 진화, 응력 생성 상호작용 조명. 고체 열역학 평형 기반, 원자 확산 응력 영향, 그 반대도 성립. Larche와 Cahn 프레임워크: 작고 등방적 변형 가정 응력 의존 화학 포텐셜 도입. Wu 공식: Eshelby 운동량 텐서 포함. Cui 등 유한 변형 확장, 그러나 낮은 변형률(확산 대비) 가정 시 정확. 실리콘 급속 충전 시 큰 부피 팽창 고려 시 이 가정 큰 오류 유발 가능.

본 연구: 저변형률 가정 없는 빠른 변형 유효 화학 포텐셜 표현 유도, 전극 형상 독립적. 실리콘 나노입자 전극 정전위/정전류 제어 조건 적용, 응력과 리튬 분포 영향 정량화. 두 모델 수렴 임계 반경 규명, 응력 영향 공간적 한계 통찰 제공.

방법

역학 방정식

리튬 삽입 부피 변화 유발. 라그랑주와 오일러 좌표로 변형 기술:
$$\mathbf{U}=\mathbf{x}-\mathbf{X}\quad (1)$$ U 변위장. 변형 기울기 텐서:
$$\mathbf{F}=\frac{\partial \mathbf{x}}{\partial \mathbf{X}}=\mathbf{I}+\mathrm{Grad}\,\mathbf{U}\quad (2)$$

구형 입자 좌표 (R,Θ,Φ)(라그랑주), (r,θ,φ)(오일러)서 기울기 텐서 단순화:
$$\mathbf{F}=\begin{bmatrix}{F}_R&0&0\0&{F}_{\Theta}&0\0&0&{F}_{\Phi}\end{bmatrix}=\begin{bmatrix}1+\partial u/\partial R&0&0\0&1+u/R&0\0&0&1+u/R\end{bmatrix}\quad (3)$$

변형 리튬화 유도 형상 변화와 가역 탄성 변형 분해:
$$\mathbf{F}=\mathbf{F}^e\mathbf{F}^c\quad (4)$$ F^c 등방:
$$\mathbf{F}^c=(1+\Omega C)^{1/3}\mathbf{I}\quad (5)$$ Ω 편몰 부피. 결과:
$$\mathbf{F}^e=(1+\Omega(R)C)^{-1/3}\begin{bmatrix}1+\partial u/\partial R&0&0\0&1+u/R&0\0&0&1+u/R\end{bmatrix}\quad (6)$$

그린-라그랑주 응변 텐서 E:
$$\mathbf{E}=\frac{1}{2}\left(\mathbf{F}^T\mathbf{F}-\mathbf{I}\right)\quad (7)$$ 반경, 접선 성분:
$$E_R^e=\tfrac{1}{2}\Big[\tfrac{(1+\partial u/\partial R)^2}{(1+\Omega(R)C)^{2/3}}-1\Big]\quad (9)$$
$$E_{\Theta}^e=E_{\Phi}^e=\tfrac{1}{2}\Big[\tfrac{(1+u/R)^2}{(1+\Omega(R)C)^{2/3}}-1\Big]\quad (10)$$

선형 탄성 가정, 응변 에너지 밀도:
$$W=\tfrac{1}{2}\mathbf{E}^e:\mathbf{C}:\mathbf{E}^e=\det(\mathbf{F}^c)\frac{E_h}{2(1+\nu)}\Big[\tfrac{\nu}{1-2\nu}\big(tr\,\mathbf{E}^e\big)^2+tr\,\mathbf{E}^e\mathbf{E}^e\Big]\quad (12)$$ Eₕ, ν 영률, 푸아송비. 제1 피올라-커르호프 응력:
$$\mathbf{P}=\det(\mathbf{F}^c)\frac{E_h}{1+\nu}\Big[\tfrac{\nu}{1-2\nu}tr\,\mathbf{E}^e+\mathbf{E}^e\Big]\mathbf{F}^e(\mathbf{F}^c)^{-1}\quad (13)$$

성분별:
$$P_R=\left(1+\Omega C\right)^{1/3}\frac{E_h}{2(1+\nu)(1-2\nu)}\left(1+\tfrac{\partial u}{\partial R}\right)\Big[(1-\nu)E_R^e+2\nu E_{\Theta}^e\Big]\quad (14)$$
$$P_{\Theta}=P_{\Phi}=\left(1+\Omega C\right)^{1/3}\frac{E_h}{2(1+\nu)(1-2\nu)}\left(1+\tfrac{u}{R}\right)\Big[\nu E_R^e+E_{\Theta}^e\Big]\quad (15)$$

체적 힘 없는 평형:
$$\frac{\partial P_R}{\partial R}+\frac{2(P_R-P_{\Phi})}{R}=0\quad (16)$$ 경계 조건:
$$u(0,t)=0,\;P_R(R_0,t)=0\quad (17)$$

질량 수송 방정식

리튬 농도 C(X,t), 플럭스 J(X,t):
$$\frac{\partial C}{\partial t}+\mathrm{Div}\,\mathbf{J}=0\quad (18)$$. 구대칭 시:
$$\frac{\partial C(R,t)}{\partial t}+\frac{1}{R^2}\frac{\partial\big(R^2J(R,t)\big)}{\partial R}=0\quad (19)$$. 플럭스 화학 포텐셜 구동력 피크 법칙:
$$J=-\frac{CD}{R_gTF_{11}F_{11}}\frac{\partial \mu}{\partial R}\quad (20)$$.

화학 포텐셜 총 내부 에너지 밀도
$$\Pi=\varphi(C)+W(C,\mathbf{E}^e)\quad (22)$$ 유도:
$$\mu=\frac{\partial\Pi}{\partial C}=\mu_0(C)+\tau(\mathbf{E}^e,C)\quad (23)$$ τ 응력 의존 부분. 희석 용액:
$$\mu_0(C)=\mu_0+R_gT\ln(\gamma C)\quad (24)$$ γ=1.

전통적 헬름홀츠 자유 에너지 접근:
$$\tau_H=\frac{\partial W}{\partial C}= -\det(\mathbf{F}^e)\sigma_m\Omega\quad (25)$$ 결과:
$$\mu={\mu}_0+R_gT\ln C-\det(\mathbf{F}^e)\Omega\sigma_m\quad (26)$$. σₘ 코시 정수압 응력:
$$\sigma_m=\tfrac{1}{3}tr\,(\det^{-1}\mathbf{F}\mathbf{PF}^T)\quad (27)$$. 식 (26) 전통 화학 포텐셜 모델.

대안 깁스 자유 에너지 처리:
$$\tau_G=\frac{\partial W}{\partial C}\quad (28)$$ 결과:
$$\mu={\mu}_0+R_gT\ln C-\frac{\partial W}{\partial C}\quad (29)$$, 본 연구 개발 모델. 질량 수송 시스템 (19),(20),(26),(29) 결합. 정전위(C(R₀,t)=Cₘₐₓ), 정전류(J(R₀,t)=j₀(1+u/R)²) 경계 조건:
$$C(R_0,t)=C_{\max},\;J(R_0,t)=j_0\big(1+u/R\big)^2\quad (30-31)$$ 초기:
$$C(R,0)=0\quad (32)$$

수치 구현

해석적 해 난해, COMSOL Multiphysics로 결합 PDE 풀이, 표 1 물성치 사용. 무차원 변수 명료성 위해 사용. 충전 상태(SOC):
$$SOC=\frac{\int_0^{R_0}C(R,t)R^2dR}{\int_0^{R_0}C_{\max}R^2dR}\quad (33)$$ 두 화학 포텐셜 모델 유도 플럭스:
$$J_H=\frac{\partial\tau_H}{\partial R},\;J_G=\frac{\partial\tau_G}{\partial R}\quad (34)$$

결과 및 토의

그림 1: 구형 Si 나노입자 정전류 사이클링 시 리튬 농도, 반경 응력, 후프 응력 프로파일. 개발 모델(실선), 전통 모델(삼각형) 거의 동일 농도 분포 모든 SOC. 응력 프로파일 미세 차이, 46.7%, 65.5% SOC서 중심부 약간 상승.

그림 2: 정전위 사이클링 유사 데이터. 두 모델 간 차이 더 현저, 특히 중간 SOC 중심 응력. 정전위 동작 표면 농도 고정 → 충전률, 결과적 변형률 정전류보다 높음. 주어진 SOC 전체 변형 유사 → 핵심 요인 절대 변형 아닌 변형률임.

그림 3(참조 [17]): 620 nm Si 나노입자 2 V vs Li서 1분 내 완전 팽창 가능 → 빠른 변형 모델링 중요성 입증. 이런 영역서 두 모델 현저히 다른 응력 상태 예측, 실험 응력 데이터 부재로 정량 검증 불가.

그림 4: 정전류 사이클링 다양한 충전 전류서 플럭스 비율 Jₕ/JG. 입자 전체서 1 근처, 중심부 근접 약간 증가, 표면 근처 감소. 전류 증가 시 벌림 확대 → 변형률 민감도 강함. 모든 곡선 임계 반경 근처 교차 → 화학 포텐셜 독립 영역(CIR) 명명, 플럭스 모델 무관 수렴. CIR 전류 증가 시 표면 쪽으로 이동.

그림 5,6: 두 모델 코시 정수압 응력 분포 비교. CIR 내 σₘ 리튬화 대부분 기간 ~0 → 모델 일치 이유. CIR 밖 σₘ 발산 → 화학 포텐셜, 플럭스 차이 증폭.

결론

빠른 변형 유효, 저변형률 가정 제거 화학 포텐셜 표현 도입. 전통 공식 비교: 정전류 사이클링 무시할 차이, 정전위 제어 현저 편차 → 변형률 주도. 임계 반경(CIR) 등장, 플럭스·응력 사실상 동일, 고전류서 표면 쪽으로 이동. 실리콘 나노입자 전극 급속 충전 시 응력 진화 정확 예측 위해 개발 모델 필요성 부각.