REAL 에너지 디바이스 가이드북

4. 수송: 확산과 이동

학습 목표
  • 화학종의 움직임을 확산과 이동으로 나누고, 하나의 전기화학 퍼텐셜 기울기에서 둘을 합쳐 Nernst-Planck 플럭스로 만든다.
  • 고체 입자와 전해질의 두 확산 시간 상수 τ=L2/D\tau = L^2/D를 어림하고, 각 시계가 몇 C에서 추월당하는지 예측한다.
  • Bruggeman 보정으로 벌크 확산계수를 다공성 전극의 유효 확산계수로 바꾸고, 전해질 전위 법칙 i=κϕκDlnci = -\kappa\nabla\phi - \kappa_D\nabla\ln c를 읽어낸다.
요약
  • 동역학은 계면에서 일어나고, 수송이 그 계면을 계속 공급합니다. 배터리·연료전지의 속도 한계 대부분이 여기서 나옵니다.
  • 두 확산 시계 τL2/D\tau \sim L^2/D가 방전 곡선을 지배하고, 둘 중 하나라도 앞지르면 농도 기울기, 곧 회복 가능한(열화가 아닌) 손실이 생깁니다.
  • 다공성 전극은 Bruggeman 세금을 냅니다: 기공이 좁고 구불거려 유효 수송이 ε1.5\varepsilon^{1.5}로 줄어듭니다.

동역학은 계면에서 일어나고, 수송은 그 계면을 계속 공급하는 투쟁입니다. 배터리·연료전지의 속도 한계 대부분은 수송 이야기입니다.

화학종이 움직이는 세 가지 방식

  • 확산: 농도 기울기를 따라, 플럭스 =Dc= -D\,\nabla c. 어디에나 있는 주력: 고체 입자 속 리튬, 전해질 속 염, 연료전지 층 속 용존 기체.
  • 이동(migration): 전기장이 미는 하전종. 이성분 배터리 전해질에서 양이온은 전류의 t+t_+(~0.36)만 나릅니다. 이것이 부하가 걸리면 염이 한쪽 전극에 쌓이고 반대쪽에서 얇아지는 이유입니다 .
  • 대류: 벌크 흐름. 배터리 셀 내부에선 무시 가능; 연료전지 기체 유로에서는 지배적, 설계된 공급 시스템입니다(16. 물 관리).

Nernst-Planck 플럭스

확산과 이동은 억지로 붙여 놓은 두 수송 법칙이 아니라, 한 이온이 느끼는 두 기울기를 한 줄로 쓴 것입니다. 둘을 통합하는 대상이 전기화학 퍼텐셜이고, 8장의 전해질 전위 법칙 전체가 여기서 흘러나옵니다.

Nernst-Planck 플럭스와 전해질 전위유도 보기
이 유도가 딛는 것

전해질 속 이온은 화학적 정체성과 전하를 동시에 지니므로, 실제로 그것을 미는 퍼텐셜은 전기화학 퍼텐셜입니다.

μˉi=μi+RTlnai+ziFϕ\bar\mu_i = \mu_i^\circ + RT\ln a_i + z_i F\phi

여기서 aia_i는 활동도(위 칩이 그 정의를 제공)이고 ϕ\phi는 국소 정전 퍼텐셜입니다. 선형 비가역 열역학은 화학종이 자기 전기화학 퍼텐셜 기울기를 따라, 자신의 확산계수가 정하는 속도로 흐른다고 말합니다.

Ni=DiciRTμˉiN_i = -\frac{D_i c_i}{RT}\,\nabla\bar\mu_i

계수 Di/RTD_i/RT는 확산계수와 이동도를 잇는 Nernst-Einstein 관계입니다. μˉi\bar\mu_i를 넣고 이상용액 활동도 ai=ci/ca_i = c_i/c^\circ를 쓰면 RTlnai=(RT/ci)ciRT\,\nabla\ln a_i = (RT/c_i)\,\nabla c_i이므로,

Ni=Dici    ziDiFRTciϕN_i = -D_i\,\nabla c_i \;-\; \frac{z_i D_i F}{RT}\,c_i\,\nabla\phi

첫 항이 바로 Fick 법칙(ax:fick 칩)입니다: 농도 기울기를 따르는 플럭스. 둘째 항은 이동, 즉 장 E=ϕ\mathbf{E} = -\nabla\phi가 끄는 전하입니다. 둘을 합친 것이 Nernst-Planck 플럭스이며, 모든 다공성 전극 모델의 근간입니다 .

이제 플럭스를 전해질의 전위 법칙으로 바꿉니다. 이온 전류는 화학종 플럭스의 전하 가중합입니다.

i=FiziNii = F\sum_i z_i N_i

보조정리순수 이동은 곧 Ohm 법칙보조정리 보기
이 유도가 딛는 것

농도가 균일한 영역을 택하면 ci=0\nabla c_i = 0이라 각 NiN_i의 이동 항만 살아남습니다. 그러면

i=Fizi(ziDiFRTciϕ)=F2RTizi2Diciκϕi = F\sum_i z_i\left(-\frac{z_i D_i F}{RT}c_i\nabla\phi\right) = -\underbrace{\frac{F^2}{RT}\sum_i z_i^2 D_i c_i}_{\displaystyle \kappa}\,\nabla\phi

이것이 Ohm 법칙 i=κϕi = -\kappa\nabla\phi이며, 이온 전도도 κ\kappa는 괄호에서 바로 읽힙니다. 전하 보존 i=0\nabla\cdot i = 0이 벌크 전해질에서 성립하면 κϕ\kappa\nabla\phi는 국소 우연이 아니라 전류 연속성 진술로 승격됩니다. 그러니 Ohm 법칙은 별도의 가정이 아니라, 농도 기울기를 끈 Nernst-Planck 플럭스입니다. ∎

농도 기울기를 되살리면 합에서 확산 항이 더는 상쇄되지 않습니다. 이성분 전해질에서 남는 항은 lnc\nabla\ln c에 비례하고, 그 앞계수를 유효 전도도 κD\kappa_D로 묶으면 (이 κD\kappa_D는 수송수 t+t_+와 활동도 계수 보정 1+dlnf/dlnc1 + \mathrm{d}\ln f/\mathrm{d}\ln c를 품습니다) 농축용액 전해질 전위 법칙이 나옵니다.

i=κϕ    κDlnci = -\kappa\,\nabla\phi \;-\; \kappa_D\,\nabla\ln c

첫 항은 옴 항이고, 둘째 항은 부하 아래 염이 쌓일 때마다 나타나는 확산 전위입니다. 이 한 줄이 P2D 모델의 전해질상 방정식입니다 . ∎

뿌리 사슬μˉi=μi+RTlnai+ziFϕ\bar\mu_i = \mu_i^\circ + RT\ln a_i + z_iF\phiNi=DiciRTμˉiN_i = -\tfrac{D_ic_i}{RT}\nabla\bar\mu_iNi=DiciziDiFRTciϕN_i = -D_i\nabla c_i - \tfrac{z_iD_iF}{RT}c_i\nabla\phii=κϕκDlnci = -\kappa\nabla\phi - \kappa_D\nabla\ln c

배터리의 두 시계

그림 · 방전 중 입자 내부. 리튬이 표면에서 빠져나가지만 고체 확산이 유한해 중심은 아직 차 있다. 표면 고갈이 클수록 전압이 일찍 컷오프에 닿는다. 확산 시계 τ = r_s²/D_s가 이 기울기를 정한다.
그림 · 방전 중 입자 내부. 리튬이 표면에서 빠져나가지만 고체 확산이 유한해 중심은 아직 차 있다. 표면 고갈이 클수록 전압이 일찍 컷오프에 닿는다. 확산 시계 τ = r_s²/D_s가 이 기울기를 정한다.Gemini로 생성한 그림

방전 곡선에서 보는 거의 모든 것을 두 확산 시간이 지배합니다:

τsolid=rs2Ds,τelyteL2De\tau_\mathrm{solid} = \frac{r_s^2}{D_s}, \qquad \tau_\mathrm{elyte} \sim \frac{L^2}{D_e}

왜 τ ∼ L²/D 인가 (Fick 제2법칙에서 나오는 확산 시간)유도 보기
이 유도가 딛는 것

재료는 Fick 제2법칙 하나뿐입니다:

ct=D2cx2\frac{\partial c}{\partial t} = D\frac{\partial^2 c}{\partial x^2}

여기서 단위를 벗겨냅니다. 길이는 층 두께 LL을, 시간은 아직 모르는 특성 시간 τ\tau를 단위로 재면, 즉 x=Lx^x = L\hat{x}, t=τt^t = \tau\hat{t}이고 x^,t^\hat{x},\hat{t}는 1 정도 크기입니다. 미분은 t=τ1t^\partial_t = \tau^{-1}\partial_{\hat t}, x2=L2x^2\partial_x^2 = L^{-2}\partial_{\hat x}^2로 스케일되므로

1τct^=DL22cx^2\frac{1}{\tau}\frac{\partial c}{\partial \hat{t}} = \frac{D}{L^2}\frac{\partial^2 c}{\partial \hat{x}^2}

양변에 τ\tau를 곱하면:

ct^=DτL22cx^2\frac{\partial c}{\partial \hat{t}} = \frac{D\tau}{L^2}\frac{\partial^2 c}{\partial \hat{x}^2}

무차원 방정식은 그 유일한 계수가 1 정도일 때만 균형이 맞으므로, Dτ/L21D\tau/L^2 \sim 1이 확산 시간을 고정합니다:

τL2D\tau \sim \frac{L^2}{D}

L2L^2의 제곱이 곧, 입자 반경을 반으로 줄이면 확산 시간이 1/4로 주는 이유 전부입니다. 급속충전 전극 설계에서 가장 강력한 지렛대 하나입니다. ∎

뿌리 사슬Fick 제2법칙무차원화 x=Lx^,  t=τt^x=L\hat x,\;t=\tau\hat tDτ/L21D\tau/L^2 \sim 1τL2/D\tau \sim L^2/D

입자 내 고체 확산(Ds1014m2/sD_s \sim 10^{-14}\,\mathrm{m^2/s}, 반경 수 µm)은 분 단위, 셀을 가로지르는 전해질 확산은 수십 초 단위입니다. 시계보다 빠르게 전류를 밀면 기울기가 벌어집니다: 표면 화학량론이 벌크보다 앞서 달리거나 (스튜디오 결과의 표면-vs-벌크 θ 패널이 정확히 이 신호), 전극 안쪽 끝의 전해질 농도가 무너집니다. 둘 다 고율에서 잃는 용량으로 나타나고, 진짜 열화와 달리 휴지 후 회복됩니다.

실전 예제· NMC811/Gr 21700 셀의 두 확산 시계

고체와 전해질 확산 시간을 어림하고, 1C 방전(3600 s)과 5C 방전(720 s)에 견주어 봅니다.

rsr_s (양극 입자 반경)5 μm5~\mu\mathrm{m}
DsD_s (고체 확산계수)1×1014 m2/s1\times10^{-14}~\mathrm{m^2/s}
τsolid=rs2/Ds\tau_\mathrm{solid} = r_s^2/D_s2500 s42 min2500~\mathrm{s} \approx 42~\mathrm{min}
LL (양극 + 분리막 구간)100 μm100~\mu\mathrm{m}
DeD_e (전해질 확산계수)3×1010 m2/s3\times10^{-10}~\mathrm{m^2/s}
τelyte=L2/De\tau_\mathrm{elyte} = L^2/D_e33 s33~\mathrm{s}

1C 방전(3600 s)은 고체 시계(2500 s)보다 느려서 고체 확산이 거의 따라옵니다. 5C 방전(720 s)은 2500 s 안에 한참 들어가므로, 입자 표면이 중심보다 훨씬 먼저 비고 표면-vs-벌크 θ 간극이 벌어집니다. 전해질 시계(33 s)는 둘 다에 비해 짧아, 전해질 고갈은 과도적으로 또는 아주 높은 율에서만 물립니다.

입자 확산 탐색기 (실시간). D_s와 C-rate를 끌면 확산 시계 τ=r²/D와 반경 농도 프로파일이 갱신됩니다 (r_s = 5 µm 고정).
τ ≈ 42 min · Fo ≈ 1.44
확산 영역기울기 형성
0.250.51 (surface)반경 r / r_s0.00.30.61.1c / c_max

확산 시계 τ = r_s²/D_s를 방전 시간과 비교: Fo = t_방전/τ 가 크면 표면과 중심이 같이 차고(잘 섞임), 작으면 표면만 차고 중심은 얼어붙는다. c(r) 형상은 정성적 스케일링이며, 실제 프로파일은 스튜디오 P2D 솔버가 적분한다.

다공성 매질: Bruggeman 세금

그림 · 굴곡도와 Bruggeman 보정. 이온은 빈 통로처럼 곧게 못 가고 입자 사이를 굽이굽이 돌아야 하므로, 유효 전도도와 확산이 벌크값보다 낮다. Bruggeman은 이를 부피분율의 1.5제곱으로 근사한다: κ_eff = ε^1.5·κ. 전극이 촘촘할수록(ε 작을수록) 벌칙이 커진다.
그림 · 굴곡도와 Bruggeman 보정. 이온은 빈 통로처럼 곧게 못 가고 입자 사이를 굽이굽이 돌아야 하므로, 유효 전도도와 확산이 벌크값보다 낮다. Bruggeman은 이를 부피분율의 1.5제곱으로 근사한다: κ_eff = ε^1.5·κ. 전극이 촘촘할수록(ε 작을수록) 벌칙이 커진다.Gemini로 생성한 그림

전극은 스펀지입니다. 각 상은 부피의 ε\varepsilon만 차지하고 길은 구불구불합니다. 유효 물성은 그만큼 줄어드는데, 고전적으로 Deff=Dε1.5D_\mathrm{eff} = D\,\varepsilon^{1.5} (스튜디오 스키마의 brug 지수). 같은 아이디어가 연료전지 GDL의 기체 확산 가격도 매깁니다. 다만 액체 물이 기공을 차지해 열린 분율을 동적으로 깎는다는 반전이 추가됩니다.

Bruggeman 지수 1.5는 어디서 오는가유도 보기
이 유도가 딛는 것

이건 정리가 아니라 닫음(closure)입니다. 무엇을 믿고 있는지 알려면 한 번 봐 둘 가치가 있습니다. 두 가지 기하학적 벌점이 벌크 확산계수 DD를 유효 확산계수로 바꿉니다. 첫째, 어떤 단면이든 상 분율 ε\varepsilon만 수송에 열려 있습니다. 둘째, 열린 길이 구불거려 두께 LL을 건너는 화학종은 실제로 더 먼 거리를 가고, 그 비가 굴곡도 τt1\tau_t \ge 1입니다. 합치면,

Deff=ετtDD_\mathrm{eff} = \frac{\varepsilon}{\tau_t}\,D

빠진 조각은 τt(ε)\tau_t(\varepsilon)인데, 실제 무작위 미세구조에는 정확한 꼴이 없습니다. Bruggeman의 1935년 유효매질 해석(구의 무작위 현탁)이 고전적 경험 닫음을 내놓습니다:

τt=ε1/2\tau_t = \varepsilon^{-1/2}

대입하면 두 벌점이 지수 하나로 뭉칩니다:

Deff=εε1/2D=ε3/2DD_\mathrm{eff} = \frac{\varepsilon}{\varepsilon^{-1/2}}\,D = \varepsilon^{3/2}D

그래서 Deff=Dε1.5D_\mathrm{eff} = D\,\varepsilon^{1.5}입니다. 지수는 법칙이 아니라 피팅 파라미터로 보세요: 실측 전극은 1.5~4를 요구할 때가 많고, 그래서 스튜디오가 brug를 상수가 아니라 손잡이로 노출합니다. ∎

뿌리 사슬Deff=ετtDD_\mathrm{eff} = \tfrac{\varepsilon}{\tau_t}DBruggeman 닫음 τt=ε1/2\tau_t = \varepsilon^{-1/2}Deff=Dε3/2D_\mathrm{eff} = D\varepsilon^{3/2}

기체상 수송 (연료전지)

기체는 액상 이온보다 45자릿수 빨리 확산하지만(D105D \sim 10^{-5} vs 1010m2/s10^{-10}\,\mathrm{m^2/s}), 연료전지는 배터리보다 12자릿수 높은 전류밀도로(정격 ~1 A/cm² vs 1C에서 수 mA/cm²), 역류하는 생성수와 경쟁하는 층을 통과해 돌아갑니다. 14장의 한계 전류 jLj_L은 산소가 소비되는 딱 그 속도로만 도착하는 지점입니다. 그 너머엔 아무것도 없습니다.

고체 시계를 추월해 θ 간극이 벌어지는 것 보기

NMC811/Gr 21700 프리셋을 5C, 25 °C로 방전하고 표면-vs-벌크 θ 패널을 여세요.

  • 표면 화학량론이 벌크보다 앞서 달립니다: 입자가 바깥에서 안으로 비는, 워크드 예제의 τsolid\tau_\mathrm{solid} 이야기 그대로입니다.
  • 방전 뒤 셀을 휴지시키면 간극이 풀립니다. 이 잃은 용량은 회복 가능한 농도 분극이지 열화가 아닙니다.
스튜디오 열기 · 설정 자동 입력

연습문제

Q1

워크드 예제는 rs=5 μmr_s = 5~\mu\mathrm{m}를 썼습니다. 어느 공급사가 같은 화학을 rs=2.5 μmr_s = 2.5~\mu\mathrm{m}로 분쇄해 내놓습니다. 새 고체 확산 시간은 얼마이고, 몇 배로 바뀌며, τ=L2/D\tau = L^2/D의 어느 항이 그 원인입니까?

풀이 보기

τsolid=rs2/Ds\tau_\mathrm{solid} = r_s^2/D_s는 반경의 제곱에 비례합니다. rsr_s를 반으로 줄이면 τ\tau22=42^2 = 4로 나뉩니다: 2500625 s10 min2500 \to 625~\mathrm{s} \approx 10~\mathrm{min}. 길이의 제곱 의존(L2L^2, 여기서는 rs2r_s^2)이 원인이며, 그래서 입자 미세화가 급속충전 전극 설계에서 가장 강력한 단일 지렛대입니다.

Q2

스튜디오 실습. 5C에서 NMC811/Gr 셀은 C/2보다 적은 용량을 냅니다. 그 사라진 용량이 회복 가능한 농도 분극인지 진짜 열화인지 판정하는, 휴지 단계를 쓰는 절차를 설계하라.

풀이 보기

5C로 컷오프까지 방전해 낸 용량을 기록한 뒤, 개방 회로로 휴지시키고(수 분에서 한 시간) C/20 같은 낮은 율로 계속 방전합니다. 휴지 동안 농도 기울기가 풀리므로(전류가 없으면 τ\tau 시계가 거꾸로 돕니다) 회복 가능한 손실은 휴지 후 낮은 율의 추가 용량으로 돌아옵니다. 진짜 열화(활물질 또는 리튬 재고 상실)는 돌아오지 않습니다. 휴지 뒤 회복된 분율이 농도 분극 몫이고, 영구히 남은 나머지가 열화입니다.

이 장이 다루지 못하는 것

결합 다상 유동(한 기공망 속 액체+증기+기체), 멤브레인의 전기삼투 끌림(16장), 실리콘의 응력 결합 확산(11. 역학 (응력-전위) 모델).

참고문헌

[1] J. Newman, K. E. Thomas-Alyea (2004). Electrochemical Systems, 3rd ed. Wiley. 원문 링크 ↗

[2] M. Doyle, T. F. Fuller, J. Newman (1993). Modeling of Galvanostatic Charge and Discharge of the Lithium/Polymer/Insertion Cell J. Electrochem. Soc. 140, 1526. doi:10.1149/1.2221597 ↗

함께 보기

8. P2D 지배방정식 · 14. 분극 곡선 읽기 · 16. 물 관리 · Newman 11~12장; Bruggeman, Ann. Phys. 416, 636 (1935).

REAL

Renewable Engine Analysis Laboratory

Hanyang University ERICA

경기도 안산시 상록구 한양대학로 55 공학관V 502호

Rm #502, Engineering Bldg. V, 55 Hanyangdaehak-ro, Sangnok-gu, Ansan, Gyeonggi-do, Korea

© 2026 REAL — Department of Mechanical Engineering, Hanyang University ERICA. All rights reserved.