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

16. 물 관리

학습 목표
  • 멤브레인을 지나는 순 물 유속을 전기삼투 끌림과 Fick 역확산으로 나누고, 양성자 유속과 Fick 제1법칙에서 균형 Nwater=ndj/FDwdcw/dxN_\mathrm{water} = n_d\,j/F - D_w\,dc_w/dx를 조립합니다.
  • 끌림이 역확산을 앞질러 음극 쪽이 마르기 시작하는 문턱 전류를 계산하고, 멤브레인 두께가 그 문턱을 어떻게 옮기는지 예측합니다.
  • 두 고장 모드(건조 vs 플러딩)를 분극 곡선과 전압 서명으로 진단하고, 연료전지 시스템이 능선에 머물기 위해 쓰는 네 제어 레버를 짚어냅니다.
요약
  • PEM 멤브레인은 양성자를 전도하려면 젖어 있어야 하고, 양극 기공은 산소가 지나가도록 열려 있어야 한다. 물 관리는 그 능선 위의 삶이다.
  • 멤브레인을 지나는 순(net) 물 유속은 줄다리기다: 전기삼투 끌림(ndj/Fn_d j/F, 음극→양극)에서 Fick 역확산(양극→음극)을 뺀 것.
  • 능선을 사이에 둔 두 고장 모드는 건조(전도도→발열 나선)와 플러딩(산소 굶주림)이며, 한 셀 안에 공존할 수 있다.

PEM 연료전지는 동시에 젖어 있으면서 말라 있어야 하는 기계입니다: 멤브레인은 양성자를 전도하려면 포화를 요구하고, 양극의 기공은 산소가 지나가도록 충분히 열려 있어야 합니다. 물 관리는 그 능선 위에서 사는 기술입니다.

물은 어디서 와서 어디로 가는가

세 흐름이 겹칩니다:

  • 생성: 모든 암페어가 양극 촉매층에서 물을 만듭니다(면적당 j/2Fj/2F, 부하에 가차 없이 비례).
  • : 멤브레인을 건너는 양성자마다 물 분자 1에서 2.5개를 음극→양극으로 끌고 갑니다 . 그래서 고전류는 생성이 양극을 잠기게 하는 바로 그 순간 멤브레인의 음극 쪽을 말립니다.
  • 역확산: 끌림과 생성이 만든 농도 기울기가 물을 양극→음극으로 되밀어 비대칭을 부분 치유합니다. 얇은 멤브레인이 역확산에 유리한 것도 얇은 쪽이 이기는 또 하나의 이유입니다.

여기에 설계된 흐름(가습된 입구 기체, GDL을 지나 채널로의 액체 제거)을 더하면 균형 잡을 수도꼭지가 다섯입니다.

유속 균형

그림 · PEMFC 물 균형. 양성자가 멤브레인을 건널 때 물을 끌고 가는 전기삼투 끌림(음극→양극, 전류에 비례)과, 양극에서 생성된 물이 다시 확산하는 역확산(양극→음극)이 경쟁한다. 고전류에서 끌림이 이기면 음극이 말라 멤브레인 저항이 치솟는다.
그림 · PEMFC 물 균형. 양성자가 멤브레인을 건널 때 물을 끌고 가는 전기삼투 끌림(음극→양극, 전류에 비례)과, 양극에서 생성된 물이 다시 확산하는 역확산(양극→음극)이 경쟁한다. 고전류에서 끌림이 이기면 음극이 말라 멤브레인 저항이 치솟는다.Gemini로 생성한 그림

앞의 세 흐름은 독립된 장부 항목이 아닙니다: 한 멤브레인을 공유하므로 수화를 지배하는 것은 그 입니다. 멤브레인을 건너는 물을 따라가 봅시다.

멤브레인 물 유속 균형유도 보기
이 유도가 딛는 것

음극(anode)→양극(cathode) 방향을 양(+)으로 잡고, xx를 음극 면(x=0x=0)에서 양극 면(x=Lx=L)까지 두어 물의 몰 유속을 따라갑니다.

끌림. 전류밀도 jj는 양성자 유속 j/Fj/F를 정하는데, 그 양성자 유속은 가정이 아니라 두 공리의 장부 결과이니, 그 위에 쌓기 전에 증명해 둘 값어치가 있습니다.

보조정리양성자 유속은 j/F이다 (Faraday와 전하 보존)보조정리 보기
이 유도가 딛는 것

멤브레인은 전자를 막으므로 정상 상태에서 전하 보존이 이온 전류 전체를 양성자가 나르도록 강제합니다: 양성자 전류밀도는 셀 전류밀도 jj와 같습니다. Faraday 상수 FF는 단위 전하 운반체 1몰당 전하이므로, 전류밀도 jj(면적·시간당 전하)는 양성자 몰 유속

NH+=jFN_{\mathrm{H^+}} = \frac{j}{F}

(면적·시간당 몰 H⁺)으로 전달됩니다. 동역학은 들어오지 않고, 순수한 전하 세기일 뿐입니다. ∎

양성자마다 물 분자 ndn_d개(전기삼투 끌림 계수)를 끌고 가니, +x+x 방향 물 유속은:

Ndrag=ndjFN_\mathrm{drag} = n_d\,\frac{j}{F}

역확산. 물은 양극에서 생성되므로 농도 cwc_w가 거기서 가장 높고 음극 쪽으로 낮아집니다. Fick 법칙이 그 기울기를 따라 확산 유속을 만듭니다 :

Ndiff=DwdcwdxN_\mathrm{diff} = -D_w\,\frac{dc_w}{dx}

cwc_w가 양극 쪽으로 높아지므로(dcw/dx>0dc_w/dx > 0) 이 항은 음수이고, 양극→음극을 가리키며 끌림을 거스릅니다.

합. 둘을 겹치고, 역확산 크기를 NbdDwdcw/dx>0N_\mathrm{bd} \equiv D_w\,dc_w/dx > 0으로 두면:

Nwater=ndFjDwdcwdx=ndFjNbdN_\mathrm{water} = \frac{n_d}{F}\,j - D_w\,\frac{dc_w}{dx} = \frac{n_d}{F}\,j - N_\mathrm{bd}

건조 조건. 음극 면은 순 유속이 양극 쪽을 가리킬 때, 즉 Nwater>0N_\mathrm{water} > 0일 때 물을 잃습니다:

ndFj>Nbdj>jFndNbd\frac{n_d}{F}\,j > N_\mathrm{bd} \quad\Longleftrightarrow\quad j > j^\star \equiv \frac{F}{n_d}\,N_\mathrm{bd}

끌림 항은 전류에 선형으로 커지고, 역확산은 멤브레인이 버틸 수 있는 가장 가파른 농도 기울기에 묶여 있습니다. 그래서 문턱 jj^\star 위에서 음극 쪽이 줄다리기에 지며 마르고, 그 순간 양극은 잠깁니다. 아래 두 고장 모드는 이 한 방정식을 양 끝에서 읽은 것입니다. ∎

뿌리 사슬양성자 유속 =j/F=j/F끌림 Ndrag=ndj/FN_\mathrm{drag}=n_d\,j/FFick 역확산 Ndiff=Dwdcw/dxN_\mathrm{diff}=-D_w\,dc_w/dxNwater=ndFjDwdcwdxN_\mathrm{water}=\tfrac{n_d}{F}\,j-D_w\tfrac{dc_w}{dx}

두 계수에는 되먹임이 숨어 있습니다: ndn_dDwD_w는 멤브레인 국소 수분 함량 λ\lambda(황산기 자리당 물 분자 수)에 따라 커지는데, 그 λ\lambda가 바로 이 방정식이 기술하는 유속으로 정해집니다. 그래서 수화(hydration)는 움직이는 표적입니다. 전류를 바꾸면 λ\lambda 분포 전체가 다시 자리 잡고, 자신을 정한 계수를 이동시킵니다.

실전 예제· 1 A/cm²에서 끌림 vs 역확산

정격 부하의 얇은 멤브레인에서 전기삼투 끌림 유속을 역확산 유속과 비교합니다. nd=1n_d = 1, 멤브레인 두께 L=25 μmL = 25~\mu\mathrm{m}, 물 확산계수 Dw=1×1010 m2/sD_w = 1\times10^{-10}~\mathrm{m^2/s}, 멤브레인을 가로지르는 수화차 Δλ=4\Delta\lambda = 4(양극이 음극보다 젖음)로 잡습니다. Nafion 계열 멤브레인의 고정 전하 농도는 cSO31800 mol/m3c_{\mathrm{SO_3}} \approx 1800~\mathrm{mol/m^3}이므로 Δcw=ΔλcSO3\Delta c_w = \Delta\lambda\,c_{\mathrm{SO_3}}입니다.

물리량
전류밀도j=1 A/cm2j = 1~\mathrm{A/cm^2}1×104 A/m21\times10^4~\mathrm{A/m^2}
끌림 유속Ndrag=ndj/F=1×104/96485N_\mathrm{drag} = n_d\,j/F = 1\times10^4/964850.104 molm2s10.104~\mathrm{mol\,m^{-2}\,s^{-1}}
물 기울기Δcw/L=(4×1800)/(25×106)\Delta c_w/L = (4\times1800)/(25\times10^{-6})2.9×108 molm42.9\times10^{8}~\mathrm{mol\,m^{-4}}
역확산 유속Nbd=DwΔcw/L=1010×2.9×108N_\mathrm{bd} = D_w\,\Delta c_w/L = 10^{-10}\times2.9\times10^{8}0.029 molm2s10.029~\mathrm{mol\,m^{-2}\,s^{-1}}
순 유속Nwater=0.1040.029N_\mathrm{water} = 0.104 - 0.029+0.075 molm2s1\mathbf{+0.075~mol\,m^{-2}\,s^{-1}}

끌림이 역확산을 약 3.6 대 1로 이기고, 순 유속이 음극→양극을 가리키므로 (Nwater>0N_\mathrm{water} > 0) 1 A/cm21~\mathrm{A/cm^2}에서 음극 면은 물을 적극적으로 잃고 있습니다. 이것이 음극 공급을 가습하는 이유이자, Δcw/L\Delta c_w/L을 가파르게 해 NbdN_\mathrm{bd}를 높이는 얇은 멤브레인이 수화 논쟁에서 이기는 이유입니다.

두 가지 고장 모드

그림 · 물 관리의 두 고장. 물이 너무 많으면 기공이 액체로 막혀 반응 기체가 못 들어가고(침수), 너무 적으면 멤브레인이 말라 양성자 저항이 치솟는다(건조). PEMFC는 이 둘 사이의 좁은 습도 창에서만 잘 돌아가며, 물 관리가 곧 그 창을 지키는 일이다.
그림 · 물 관리의 두 고장. 물이 너무 많으면 기공이 액체로 막혀 반응 기체가 못 들어가고(침수), 너무 적으면 멤브레인이 말라 양성자 저항이 치솟는다(건조). PEMFC는 이 둘 사이의 좁은 습도 창에서만 잘 돌아가며, 물 관리가 곧 그 창을 지키는 일이다.Gemini로 생성한 그림

건조. 멤브레인 수분 함량↓ → 양성자 전도도↓ → 저항 손실↑(14. 분극 곡선 읽기의 중간 기울기가 눈에 띄게 가팔라짐) → 국소 발열↑ → 더 마름. 이것이 수 분 안에 나선을 그릴 수 있는 양의 되먹임입니다. 고온·저가습·고화학량비(과잉 건조 공기)가 단골 범인.

플러딩. 액체 물이 GDL 기공과 채널 모서리를 채우고, 산소 길이 좁아지고, 한계 전류가 무너집니다. 서명: 이르고 무른 농도 추락 + 물 슬러그가 간헐적으로 빠지며 요동치는 전압. 저온·저풍량·과가습이 원인이고, 영하 시동에서 생성수가 촉매층 안에서 얼음이 되는 특수 사례도 포함합니다.

두 모드는 한 셀 안에 공존합니다: 공기 입구 근처는 마르고 출구 근처는 잠기고. 그 기울기가 채널 분해능 모델(18. 모델링 사다리의 3단)이 존재하는 이유입니다.

제어 레버

입구 가습, 셀 온도(포화압은 운전 창에서 10 K마다 대략 1.5에서 2배로 가파르게 상승), 공기 화학량비, 채널 압력 강하(물방울을 끌어냄). 운전 전략은 이 넷 위의 스케줄링 문제이고, 연료전지 시스템이 어떤 배터리에도 없는 가습기를 달고 다니는 이유입니다.

연습문제

Q1

워크드 예제의 역확산 유속 Nbd=0.029 molm2s1N_\mathrm{bd} = 0.029~\mathrm{mol\,m^{-2}\,s^{-1}}nd=1n_d = 1을 써서 순 멤브레인 물 유속이 0이 되는(끌림이 역확산과 정확히 균형) 문턱 전류밀도 jj^\star를 구하세요. 정격 1 A/cm21~\mathrm{A/cm^2}은 그보다 위입니까 아래입니까, 그리고 그것이 음극에 무엇을 뜻합니까?

풀이 보기

유속 균형에서 Nwater=0N_\mathrm{water} = 0으로 두면 ndj/F=Nbdn_d\,j^\star/F = N_\mathrm{bd}이므로:

j=FndNbd=96485×0.02912.8×103 A/m2=0.28 A/cm2j^\star = \frac{F}{n_d}\,N_\mathrm{bd} = \frac{96485 \times 0.029}{1} \approx 2.8\times10^{3}~\mathrm{A/m^2} = 0.28~\mathrm{A/cm^2}

0.28 A/cm20.28~\mathrm{A/cm^2} 아래에서는 역확산이 음극을 젖은 채로 유지하고, 그 위에서는 음극 쪽이 마릅니다. 정격 1 A/cm21~\mathrm{A/cm^2}은 약 3.6배 높으니, 정격 부하에서 음극 면은 확실히 건조 영역에 있고 외부 가습이 필요합니다. 워크드 예제에서 구한 순 유속 그대로입니다.

Q2

멤브레인을 같은 수화차 Δcw\Delta c_w와 확산계수 DwD_w에서 절반 두께 (12.5 μm12.5~\mu\mathrm{m})로 교체합니다. 역확산 유속과 Q1의 문턱 전류 jj^\star를 다시 계산하고, 얇은 멤브레인이 음극 건조에 강한 이유를 한 문장으로 설명하세요.

풀이 보기

역확산은 Nbd=DwΔcw/LN_\mathrm{bd} = D_w\,\Delta c_w/L로 스케일합니다. Δcw\Delta c_w를 고정한 채 LL을 반으로 하면 기울기가 두 배가 되므로

Nbd=2×0.029=0.058 molm2s1,j=96485×0.05815.6×103 A/m2=0.56 A/cm2N_\mathrm{bd} = 2 \times 0.029 = 0.058~\mathrm{mol\,m^{-2}\,s^{-1}}, \qquad j^\star = \frac{96485 \times 0.058}{1} \approx 5.6\times10^{3}~\mathrm{A/m^2} = 0.56~\mathrm{A/cm^2}

문턱이 두 배가 됩니다. 더 얇은 멤브레인은 같은 양 끝 수화차에 대해 더 가파른 물 기울기를 지탱하므로 역확산이 두 배 빠르고 건조 개시를 두 배 전류로 밀어냅니다. 멤브레인이 대략 175에서 20 마이크로미터 아래로 얇아진 정량적 이유입니다.

이 장이 다루지 못하는 것

정량적 2상 수송(포화도 의존 투과도, 모세관압 곡선), 동결 시동 공학, 물방울 동역학을 위한 채널 설계: 모델링 사다리 3단 이상.

참고문헌

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

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

함께 보기

15. PEMFC 해부학: MEA · 17. 열 관리 · 14. 분극 곡선 읽기 · 4. 수송: 확산과 이동 · O'Hayre 외 5장; Barbir 7장.

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.