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

14. 분극 곡선 읽기

학습 목표
  • 측정한 V(j)V(j) 곡선을 세 손실(활성화, 저항, 농도)로 분해하고, 각 손실이 차지하는 영역과 물리적 범인을 짚어냅니다.
  • 각 손실을 뿌리 법칙에서 유도합니다: 활성화 절벽은 Butler-Volmer의 Tafel 극한, 저항 기울기는 옴의 법칙, 끝의 추락은 Fick 제1법칙과 한계 전류에서.
  • ENernstE_\mathrm{Nernst}, Tafel 기울기 bb, 면적당 저항, 한계 전류로부터 주어진 전류밀도에서의 셀 전압을 계산하고, 출력과 효율을 읽어냅니다.
요약
  • 분극 곡선 V(j)V(j)는 셀의 이력서입니다: 가역 전압에서 세 손실(활성화, 저항, 농도)을 왼쪽부터 차례로 뺀 것.
  • 손실마다 영역과 범인이 있습니다: 낮은 jj의 ORR 양극 절벽, 중간의 가습이 좌우하는 저항 기울기, 한계 전류 jLj_L의 물질전달 추락.
  • 출력 P=VjP = Vj는 떨어지는 전압과 오르는 전류가 만나는 곳에서 정점을 찍고, 정격점은 그 정점보다 충분히 왼쪽에서 출력을 조금 내주고 효율을 얻습니다.

분극 곡선(셀 전압 vs 전류밀도)은 연료전지의 이력서입니다. 한 번의 스윕으로 촉매·멤브레인·기체 공급의 건강을 알 수 있습니다. 어느 영역이 어느 손실의 책임인지 안다면.

세 손실의 해부학

분극 곡선의 세 영역전류밀도에 따른 셀 전압 곡선: 초기 활성화 급강하, 중간 옴 선형 감소, 말단 농도 붕괴.전류밀도 j셀 전압활성화옴 (선형)농도 (붕괴)OCV
그림 · 분극 곡선의 세 영역. 낮은 전류에서 느린 산소 반응이 활성화 급강하를 만들고(빨강), 중간에서는 멤브레인 옴 저항이 선형 감소를 주며(금색), 한계 전류 근처에서 반응물 고갈이 전압을 붕괴시킨다(파랑). 각 영역이 서로 다른 물리적 손실을 드러낸다.

Nernst 전압에서 출발해 전류축을 따라 오른쪽으로 걸으면:

V(j)=ENernstblog10(j/j0)활성화ASRj저항cln ⁣(1jjL)농도V(j) = E_\mathrm{Nernst} - \underbrace{b\log_{10}(j/j_0)}_{\text{활성화}} - \underbrace{\mathrm{ASR}\cdot j}_{\text{저항}} - \underbrace{\left|c\ln\!\left(1-\tfrac{j}{j_L}\right)\right|}_{\text{농도}}

  • 활성화 (시작의 절벽). 산소 환원 반응의 작디작은 j0j_0(3. 반응 동역학: Butler-Volmer)가 쓸모 있는 전류가 흐르기도 전에 약 0.25에서 0.4 V를 받아갑니다. Pt 위 ORR의 b60b \approx 60에서 70 mV/decade. 가팔라지면 촉매가 피독되었거나 늙은 것.
  • 저항 (곧은 중간). 멤브레인 + 접촉을 (ASR, Ω·cm²) 하나로 요약. 가습이 실시간으로 움직입니다: 마르는 멤브레인은 수 분 안에 중간 기울기가 자라는 것으로 나타납니다(16. 물 관리).
  • 농도 (끝의 추락). 산소는 jLj_L보다 빨리 촉매에 도달할 수 없고, 거기 다가가면 Nernst 항이 무너집니다. 플러딩은 jLj_L을 동적으로 깎습니다. 익사하는 양극의 서명은 이르고 무릎이 무른 추락입니다.

세 영역마다 짧은 유도가 뒤에 있습니다. 전류축을 따라 오르며 만나는 순서로.

양극이 활성화 손실을 떠안는 이유 (ORR vs HOR)유도 보기
이 유도가 딛는 것

두 전극은 직렬이라 같은 전류밀도 jj를 흘립니다. Butler-Volmer의 Tafel(고과전압) 극한(3. 반응 동역학: Butler-Volmer)에서 각 전극의 활성화 과전압은

ηact=blog10jj0\eta_\mathrm{act} = b\log_{10}\frac{j}{j_0}

입니다. 이 단일 전극 Tafel 형태는 가정이 아니라, 한 가지가 평형에서 멀리 구동될 때 Butler-Volmer 법칙이 되는 모습이니, 여기에 기대기 전에 증명해 둘 값어치가 있습니다.

보조정리단일 전극 Tafel 형태, Butler-Volmer에서보조정리 보기
이 유도가 딛는 것

Butler-Volmer를 따르는 전극 하나를 잡습니다. j=j0[exp(αaFη/RT)exp(αcFη/RT)]j = j_0[\exp(\alpha_a F\eta/RT) - \exp(-\alpha_c F\eta/RT)]. 산화 방향으로 세게 구동하면(ηRT/F\eta \gg RT/F, 실온에서 약 26 mV) 역방향 지수항이 0으로 무너지고 정방향 가지만 살아남습니다.

jj0exp ⁣(αaFηRT)j \approx j_0\exp\!\left(\frac{\alpha_a F\eta}{RT}\right)

자연로그를 취해 η\eta에 대해 풀고 밑 10으로 바꾸면:

η=RTαaFlnjj0=2.303RTαaFlog10jj0blog10jj0\eta = \frac{RT}{\alpha_a F}\ln\frac{j}{j_0} = \frac{2.303\,RT}{\alpha_a F}\log_{10}\frac{j}{j_0} \equiv b\log_{10}\frac{j}{j_0}

이라 b=2.303RT/(αaF)b = 2.303\,RT/(\alpha_a F)이고, 하나로 정해지는 decade당 고정된 mV 수입니다 . ∎

ηact\eta_\mathrm{act}에서 전극마다 다른 값은 교환 전류밀도 j0j_0 하나뿐이고, 백금 위에서 두 반응은 극과 극입니다:

j0,HOR103 to 101 A/cm2,j0,ORR108 to 106 A/cm2j_{0,\mathrm{HOR}} \sim 10^{-3}\text{ to }10^{-1}\ \mathrm{A/cm^2}, \qquad j_{0,\mathrm{ORR}} \sim 10^{-8}\text{ to }10^{-6}\ \mathrm{A/cm^2}

수소 산화는 너무 빨라 교환 전류가 산소 환원보다 5에서 6자릿수 큽니다. 같은 jj를 양쪽에 흘리고 과전압을 빼면:

ηORRηHOR=blog10j0,HORj0,ORRb×(5 to 6)0.3 V\eta_\mathrm{ORR} - \eta_\mathrm{HOR} = b\log_{10}\frac{j_{0,\mathrm{HOR}}}{j_{0,\mathrm{ORR}}} \approx b\times(5\text{ to }6) \approx 0.3\ \mathrm{V}

b60b \approx 60 mV/decade로 잡으면 그렇습니다. 음극은 몇 mV를 돌려주고 양극은 몇백 mV를 받아갑니다. 공학적 정확도에서 곡선의 활성화 항 전체가 곧 ORR 과전압이고, 그래서 연료전지 촉매 연구는 압도적으로 더 나은 산소 양극을 찾는 사냥입니다 . ∎

뿌리 사슬같은 jj, 직렬 전극BV의 Tafel 극한ηact=blog10(j/j0)\eta_\mathrm{act}=b\log_{10}(j/j_0)ηORRηHOR0.3\eta_\mathrm{ORR}-\eta_\mathrm{HOR}\approx 0.3 V
저항 손실: 멤브레인과 접촉에 걸친 옴의 법칙유도 보기
이 유도가 딛는 것

저항 항은 전기화학이 필요 없는 유일한 손실, 그냥 옴의 법칙입니다. 전류밀도 jj(A/cm²)를 셀의 직렬 저항에 흘리면 전압이 떨어지고, 그 저항을 단위 전극 면적으로 정규화한 것이 면적당 저항 RΩR_\Omega(Ω·cm²)이라, 강하는

ηohm=RΩj\eta_\mathrm{ohm} = R_\Omega\, j

전류에 선형입니다. 곡선의 곧은 중간 구간. 지배적 기여는 멤브레인을 건너는 양성자 전류입니다: 두께 LL, 이온 전도도 σ\sigma이면

RΩLσ+RcontactR_\Omega \approx \frac{L}{\sigma} + R_\mathrm{contact}

Nafion 계열 멤브레인의 σ\sigma는 건조에서 포화로 가며 대략 한 자릿수 오르므로, RΩR_\Omega(따라서 이 구간의 기울기)는 고정된 물질 상수가 아니라 가습의 실시간 지표입니다. ∎

뿌리 사슬jj가 직렬 RR을 지남단위 면적으로 정규화ηohm=RΩj\eta_\mathrm{ohm}=R_\Omega\,jRΩ=L/σ+RcontactR_\Omega=L/\sigma+R_\mathrm{contact}
농도 손실: 한계 전류와 ln(1 − j/j_L) 형태유도 보기
이 유도가 딛는 것

산소는 유로와 반응 자리 사이의 얇은 정체막(두께 δ\delta의 Nernst 확산층)을 확산으로 건너 촉매에 닿습니다. Fick 제1법칙이 그 막을 가로지르는 농도차에서 몰 플럭스를 정합니다:

N=DcbulkcsurfδN = D\,\frac{c_\mathrm{bulk} - c_\mathrm{surf}}{\delta}

O₂ 1몰이 소비될 때마다 전하 nFnF가 지나가므로 j=nFNj = nFN. 플럭스는(따라서 jj는) 표면이 완전히 굶주릴 때(csurf0c_\mathrm{surf}\to 0) 최대입니다. 그 천장이 곧 한계 전류밀도 jLj_L입니다:

jL=nFDcbulkδj_L = \frac{nFD\,c_\mathrm{bulk}}{\delta}

두 식을 나누어 DD, δ\delta, cbulkc_\mathrm{bulk}를 소거하면:

jjL=cbulkcsurfcbulk=1csurfcbulk        csurfcbulk=1jjL\frac{j}{j_L} = \frac{c_\mathrm{bulk}-c_\mathrm{surf}}{c_\mathrm{bulk}} = 1 - \frac{c_\mathrm{surf}}{c_\mathrm{bulk}} \;\;\Rightarrow\;\; \frac{c_\mathrm{surf}}{c_\mathrm{bulk}} = 1 - \frac{j}{j_L}

표면 농도가 낮아지면 전극의 평형 전위가 Nernst 항(2. 열역학: Nernst와 OCV)만큼 이동합니다. 반응물 고갈이 물리는 전압은

ηconc=RTnFlncbulkcsurf=RTnFln ⁣(1jjL)\eta_\mathrm{conc} = \frac{RT}{nF}\ln\frac{c_\mathrm{bulk}}{c_\mathrm{surf}} = -\frac{RT}{nF}\ln\!\left(1 - \frac{j}{j_L}\right)

jjLj\to j_L에서 발산합니다. 곡선 오른쪽 끝의 추락입니다. 아래 탐색기는 바로 이 ηconc\eta_\mathrm{conc}를 빼는데, 셀 전압 식에서는 +cln(1j/jL)+\,c\ln(1 - j/j_L) 항으로 나타납니다(로그 안이 분수라 음수). 여기서 경험적 c0.05c \approx 0.05 V가 RT/nFRT/nF의 역할을 맡습니다. 이 값이 순수 Nernst 값보다 조금 큰 것은 굶주린 표면이 교환 전류 j0j_0까지 떨어뜨리기 때문이라, 같은 고갈이 Nernst와 동역학에서 두 번 청구되는 셈입니다. ∎

뿌리 사슬Fick 제1법칙j=nFNj=nFNjL=nFDcbulk/δj_L=nFDc_\mathrm{bulk}/\deltaηconc=RTnFln(1j/jL)\eta_\mathrm{conc}=-\tfrac{RT}{nF}\ln(1-j/j_L)
분극 곡선 탐색기 (실시간) — 회색 = 활성화 손실만, 점선 = +저항 손실, 파랑 = 세 손실 모두(농도 손실 포함). 주황 = 면적당 출력. 슬라이더로 각 손실 영역을 분리해 보세요.
최대 출력 ≈ 1.08 W/cm² @ 1.47 A/cm²
0.00.40.81.21.6전류밀도 j (A/cm²)00.350.71.051.3V (V)00.30.60.91.2P (W/cm²)

V(j) = 1.23 − b·log₁₀(j/j₀) − ASR·j + 0.05·ln(1 − j/j_L). 교육용 0-D 모델 — 검증된 스택 모델이 아닙니다.

위 탐색기에서 슬라이더를 하나씩만 움직여 그 영역만 변형되는 걸 보세요. 그 분리가 진단 엔지니어가 측정 곡선을 기준선과 대조해 읽는 방법 그대로입니다.

실전 예제· 세 손실 법칙으로 동작점에서의 셀 전압 구하기

어떤 PEMFC 단셀이 ENernst=1.18 VE_\mathrm{Nernst} = 1.18~\mathrm{V}, j0=1×106 A/cm2j_0 = 1\times10^{-6}~\mathrm{A/cm^2}, Tafel 기울기 b=60 mV/decadeb = 60~\mathrm{mV/decade}, 면적당 저항 ASR=0.15 Ωcm2\mathrm{ASR} = 0.15~\Omega\cdot\mathrm{cm^2}, 한계 전류 jL=1.5 A/cm2j_L = 1.5~\mathrm{A/cm^2}, c=0.05 Vc = 0.05~\mathrm{V}를 가진다고 합시다. j=0.6 A/cm2j = 0.6~\mathrm{A/cm^2}에서의 셀 전압을 구합니다.

손실
활성화blog10(j/j0)=0.060log10(6×105)b\log_{10}(j/j_0) = 0.060\log_{10}(6\times10^{5})0.060×5.78=0.347 V0.060 \times 5.78 = 0.347~\mathrm{V}
저항ASRj=0.15×0.6\mathrm{ASR}\cdot j = 0.15 \times 0.60.090 V0.090~\mathrm{V}
농도cln(1j/jL)=0.05ln(0.6)\lvert c\ln(1 - j/j_L)\rvert = 0.05\,\lvert\ln(0.6)\rvert0.05×0.511=0.026 V0.05 \times 0.511 = 0.026~\mathrm{V}
셀 전압1.180.3470.0900.0261.18 - 0.347 - 0.090 - 0.0260.717 V\mathbf{0.717~V}

이 곡선 중간 지점에서 활성화가 여전히 지배적(0.347 V)이고, 저항 기울기는 적당한 0.090 V, 농도는 jjjLj_L의 40 %에 불과해 겨우 깨어난 정도(0.026 V)입니다. 숫자에서 영역 순서를 바로 읽어냅니다: 추락은 아직 오른쪽으로 한참 멀리 있습니다.

출력의 관점

출력밀도 P=VjP = V \cdot j는 오르고, 정점을 찍고, 농도 손실이 접수하면 떨어집니다. 설계점 둘이 중요합니다: 정격점(정점 출력보다 충분히 왼쪽, 효율 ≈ V/1.25V/1.25가 아직 괜찮은 곳)과 최대 출력점(짧은 버스트 전용, 그 너머는 열적으로 자해입니다).

배터리에도 있다

배터리의 V-I 관계도 비슷해 보이지만(같은 물리니까) SOC가 떨어지며 표류합니다. 배터리 시험이 정상 분극점 유지 대신 C-rate 스윕을 도는 이유 중 하나입니다. 연료전지 곡선은 정상(stationary)입니다: 공급을 고정하면 곡선이 가만히 서 있습니다. 그 정상성이 이 곡선을 그렇게 깨끗한 진단 도구로 만들고, 18장 모델링 사다리가 착취하는 성질이기도 합니다.

연습문제

Q1

워크드 예제의 동작점(j=0.6 A/cm2j = 0.6~\mathrm{A/cm^2}에서 V=0.717 VV = 0.717~\mathrm{V})을 써서 출력밀도 P=VjP = Vj를 계산하고 전압 효율 V/1.25\approx V/1.25를 추정하세요. 이 점은 최대 출력점의 왼쪽입니까 오른쪽입니까?

풀이 보기

출력밀도는

P=Vj=0.717×0.6=0.43 W/cm2P = Vj = 0.717 \times 0.6 = 0.43~\mathrm{W/cm^2}

이고 전압 효율은 0.717/1.250.570.717/1.25 \approx 0.57, 약 57 %입니다. 최대 출력점은 더 오른쪽, 아직 오르는 jj가 농도 추락에 아직 압도되지 않은 곳에 있습니다. j=0.6 A/cm2j = 0.6~\mathrm{A/cm^2}(jLj_L의 40 %에 불과)에서 곡선은 그 정점의 넉넉히 왼쪽에 있고, 바로 그곳이 정격점이 있어야 할 자리입니다: 더 높은 전압의 더 나은 효율을 위해 출력을 조금 내주는 것.

Q2

진단. 셀의 멤브레인이 마릅니다: 이온 전도도 σ\sigma가 젖었을 때 값의 절반으로 떨어지고 접촉 저항은 그대로입니다. 0.15 Ωcm20.15~\Omega\cdot\mathrm{cm^2}의 ASR 중 멤브레인 항을 0.12 Ωcm20.12~\Omega\cdot\mathrm{cm^2}, 접촉을 0.03 Ωcm20.03~\Omega\cdot\mathrm{cm^2}로 잡으세요. 새 ASR와 j=0.6 A/cm2j = 0.6~\mathrm{A/cm^2}에서의 새 저항 손실을 구하고, 곡선의 어느 영역이 움직이는지 말하세요.

풀이 보기

멤브레인 항은 RΩL/σR_\Omega \approx L/\sigma이라, σ\sigma를 반으로 하면 두 배가 됩니다: 0.120.24 Ωcm20.12 \to 0.24~\Omega\cdot\mathrm{cm^2}. 접촉은 0.030.03으로 그대로이므로

ASRnew=0.24+0.03=0.27 Ωcm2,ηohm=0.27×0.6=0.162 V\mathrm{ASR}_\mathrm{new} = 0.24 + 0.03 = 0.27~\Omega\cdot\mathrm{cm^2}, \qquad \eta_\mathrm{ohm} = 0.27 \times 0.6 = 0.162~\mathrm{V}

0.090 V0.090~\mathrm{V}에서 올라갑니다. 곧은 중간만 가팔라지고, 활성화 절벽과 한계 전류 jLj_L은 건드려지지 않습니다. 절벽과 추락은 그대로인데 중간 기울기가 수 분에 걸쳐 자라는 것이 마르는 멤브레인의 지문(16. 물 관리)이고, 그래서 가습을 이 구간에서 바로 읽어냅니다 .

이 모델이 다루지 못하는 것

위 0-D 식은 모든 층을 숫자 넷으로 뭉칩니다. 채널 방향 국소화(입구는 마르고 출구는 잠기는 동시 발생이 흔함), 과도 응답(시동, 부하 스텝), 음극/양극 손실 분리는 못 합니다. 18. 모델링 사다리의 위 칸들이 필요합니다.

참고문헌

[1] A. J. Bard, L. R. Faulkner (2001). Electrochemical Methods: Fundamentals and Applications, 2nd ed. Wiley. 원문 링크 ↗

[2] R. O'Hayre, S.-W. Cha, W. Colella, F. B. Prinz (2016). Fuel Cell Fundamentals, 3rd ed. Wiley. 원문 링크 ↗

함께 보기

13. 연료전지 기초 · 16. 물 관리 · 5. 임피던스(EIS)의 언어 · O'Hayre 외 3에서 5장.

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.