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

18. 모델링 사다리

학습 목표
  • 0-D에서 3-D까지의 사다리에서 모델링 질문을 옳은 단에 놓고, 답할 수 있는 가장 작은 모델을 고른다.
  • 수송 시간 척도(종횡비) 논증으로 두께 방향 1-D 모델을 정당화하고, 그 정당화가 깨지는 기하를 짚는다.
  • 오르기 규칙 두 가지(한 단 아래로 검증, 물리보다 파라미터가 먼저 썩음)를 실제 설계 질문을 단으로 분류하는 데 적용한다.
요약
  • 질문에 답할 수 있는 가장 작은 모델: 질문이 어느 단에 사는지 알아보는 것이 모델링 기술의 대부분이다.
  • 비용이 오르는 다섯 단: (0) 분극 방정식, (1) MEA 두께 방향 1-D, (2) 채널 방향 +1-D, (3) 셀 3-D CFD, (4) 스택 & 시스템.
  • 단은 소프트웨어가 아니라 질문이 고르게 하고, 각 단을 한 단 아래로 검증하며, 물리보다 파라미터가 먼저 썩는다는 것을 기억하라.

가장 큰 모델이 필요한 게 아니라, 질문에 답할 수 있는 가장 작은 모델이 필요합니다. 연료전지 모델링은 자연스럽게 단(rung)으로 배열되고, 질문이 어느 단에 사는지 아는 것이 이 기술의 대부분입니다. (배터리 쪽 사다리도 같은 모양입니다: 등가회로 → SPM/SPMe → P2D → 팩 네트워크, 9. 단순화 모델 (SPM / SPMe) 참조 .)

그림 · 모델링 사다리. 0D 럼프드(온도·전압 하나)에서 1D(층 방향 기울기), 2D(채널 방향), 3D(전체 형상)로 올라갈수록 잡아내는 물리와 계산 비용이 함께 커진다. 올바른 모델은 가장 정교한 것이 아니라, 질문에 답하는 가장 낮은 사다리 칸이다.
그림 · 모델링 사다리. 0D 럼프드(온도·전압 하나)에서 1D(층 방향 기울기), 2D(채널 방향), 3D(전체 형상)로 올라갈수록 잡아내는 물리와 계산 비용이 함께 커진다. 올바른 모델은 가장 정교한 것이 아니라, 질문에 답하는 가장 낮은 사다리 칸이다.Gemini로 생성한 그림

0단: 분극 방정식

14. 분극 곡선 읽기의 4-파라미터 V(j)V(j). 답할 수 있는 것: 운전점, 효율-출력 트레이드, 스택 사이징, 셀이 여러 부품 중 하나인 시스템 연구. 비용: 밀리초. 못 하는 것: 무엇도 국소화하지 못하고, 파라미터가 왜 변했는지 말하지 못합니다.

1단: MEA 두께 방향 1-D

15. PEMFC 해부학: MEA의 샌드위치를 두께 방향으로 분해합니다: 촉매층에 분포된 Butler-Volmer 소스, 이오노머 물 수송(끌림 vs 역확산), Bruggeman 보정 기체 확산, 결합 전위장. 배터리 P2D(8. P2D 지배방정식)의 연료전지 형제입니다: 같은 다공성 전극 수학에 기체상 화학종과 수분 의존 전도도가 얹힌 것 . 답: 부하 시 멤브레인 가습, 촉매층 활용도, 음극/양극 손실 분리. 비용: 초. 이 플랫폼에 미래의 "Fuel Cell 1D" 도구가 생긴다면 차지할 단이 바로 여기입니다.

2단: 채널 방향 +1-D

기체 조성이 입구에서 출구로 진화합니다: 산소는 고갈되고 물은 쌓입니다. 1단 슬라이스를 채널을 따라 행진시키는 "1+1-D" 모델은 16. 물 관리의 마른-입구/잠긴-출구 공존을 CFD 대비 하찮은 비용으로 잡습니다. 답: 화학량비 선택, 대향류 vs 병행류 배치, 가습 스케줄.

3단: 셀 3-D CFD

실제 채널/랜드 기하에서 Navier-Stokes 전체 + 2상 물. 답: 유동장 설계, 랜드 아래 전류밀도 지도, 물방울·응축 위치. 비용: 케이스당 시간 단위. 탐색 도구가 아니라 설계 검증 도구입니다.

4단: 스택 & 시스템

셀을 네트워크 노드로(12. 셀에서 팩으로의 논리), 유량을 불균등하게 나누는 매니폴드, 냉각 루프, 압축기 맵, 제어기. 각 노드 안의 전기화학은 의도적으로 0단과 1단으로 강등됩니다. 사다리가 자기 안으로 접힙니다.

언제 올라야 하는가

0단은 각 층을 하나의 균일한 상태로 뭉칩니다; 1단은 그 안의 기울기를 분해합니다. 그 도약은 언제 정말 필요할까요? 답은 봉투 뒷면에서 가늠할 수 있는 시간 척도 경주입니다.

집중(lumped) 모델이 깨지는 때: Damköhler 논증유도 보기
이 유도가 딛는 것

집중(0단) 모델은 층마다 농도가 하나로 균일하다고 가정합니다. 이는 반응이 소비하는 속도보다 화학종이 층을 더 빨리 확산해 건널 때만 성립하니, 두께 LL인 층에서 두 시간 척도를 비교합시다.

층을 건너는 확산 시간:

τdiff∼L2D\tau_\mathrm{diff} \sim \frac{L^2}{D}

소비 시간, 곧 전류밀도 jj에서 반응이 면적당 재고 cLcL로부터 화학종을 면적당 j/nFj/nF로 끌어쓰는 시간은:

τrxn∼cL j/nF =nFcLj\tau_\mathrm{rxn} \sim \frac{cL}{\,j/nF\,} = \frac{nFcL}{j}

둘의 비가 Damköhler 수, 곧 반응 속도 대 수송 속도입니다:

Da=τdiffτrxn=jLnFDc=jjlim,jlim≡nFDcLDa = \frac{\tau_\mathrm{diff}}{\tau_\mathrm{rxn}} = \frac{jL}{nFDc} = \frac{j}{j_\mathrm{lim}}, \qquad j_\mathrm{lim} \equiv \frac{nFDc}{L}

그러니 집중 가정은 사실 Da≪1Da \ll 1이라는 진술입니다: 소비가 확산 한계 전류 jlimj_\mathrm{lim}의 작은 일부에 그친다는 것. DaDa가 jj에 선형으로 커지므로, 저부하에서 충실하던 0-D 모델은 전류가 오르면 막 두께 방향과 전극 두께 방향 기울기를 소리 없이 놓칩니다. 바로 1단으로 올라서야 하는 영역입니다. ∎

뿌리 사슬τdiff∼L2/D\tau_\mathrm{diff}\sim L^2/Dτrxn∼nFcL/j\tau_\mathrm{rxn}\sim nFcL/jDa=τdiff/τrxn=j/jlimDa=\tau_\mathrm{diff}/\tau_\mathrm{rxn}=j/j_\mathrm{lim}Da≪1Da\ll 1일 때만 집중 타당

1-D가 충분한 때: 종횡비 논증

종횡비와 1D 근사얇고 넓은 셀에서 두께 방향 길이가 면 방향보다 훨씬 짧아, 두께 방향 확산이 빨라 1D 모델이 정당화된다.L⊥ ~ 100 µmL∥ ~ 0.1 m (in-plane)τ⊥ / τ∥ ~ (L⊥/L∥)² ~ 10⁻⁶두께 방향이 백만 배 빠름 → 1D면 충분
그림 · 종횡비가 1D 모델을 정당화한다. 확산 시간은 길이의 제곱에 비례하므로, 두께(~100 µm)와 면 방향(~0.1 m)의 시간비는 (L⊥/L∥)² ~ 10⁻⁶이다. 두께 방향 평형이 면 방향보다 백만 배 빨리 잡히므로, 얇고 넓은 셀은 두께 방향만 푸는 1D 모델로 충분하다.

0단에서 1단으로 오르면 두께 방향 기울기를 얻습니다. 1단은 그 기울기를 지키되 면 방향의 모든 위치를 똑같다고 취급합니다. 면 방향을 버리는 것이 언제 정직할까요? 0단을 정리한 바로 그 수송 시계가 이것도 정리하는데, 이번엔 셀의 두 방향 사이의 경주입니다.

두께 방향 1-D가 충분한 때: 종횡비 논증유도 보기
이 유도가 딛는 것

막-전극 접합체는 얇은 판입니다: 두께 방향으로 수백 마이크로미터, 면 방향으로 수십 센티미터. 화학종이 각 방향으로 그것을 확산해 건너는 데 걸리는 시간을 물읍시다. Fick 법칙은 확산 시간을 거리의 제곱을 확산도로 나눈 것으로 놓으니, 같은 DD가 양방향에 작용하면

τ⊥∼L⊥2D(두께 방향, 얇은 쪽),τ∥∼L∥2D(면 방향, 넓은 쪽)\tau_\perp \sim \frac{L_\perp^2}{D} \quad(\text{두께 방향, 얇은 쪽}), \qquad \tau_\parallel \sim \frac{L_\parallel^2}{D} \quad(\text{면 방향, 넓은 쪽})

둘을 나누면 확산도가 상쇄되고 순수 기하만 남습니다:

τ⊥τ∥∼(L⊥L∥)2\frac{\tau_\perp}{\tau_\parallel} \sim \left(\frac{L_\perp}{L_\parallel}\right)^2

두께 대 폭 비가 10−210^{-2}에서 10−310^{-3} 정도이면 그 제곱은 10−410^{-4}에서 10−610^{-6}입니다: 두께 방향이 면 방향보다 만 배에서 백만 배 빨리 평형에 이릅니다. 두께 방향 기울기가 세워져 준정상 상태에 도달하는 시간 척도에서 면 방향 분포는 거의 움직이지 않았고, 사실상 동결되어 있습니다. 그러면 모델은 빠르고 강한 두께 방향 기울기를 온전히 분해하고, 면 방향의 각 위치에 자기만의 독립적인 1-D 해를 쥐여 줄 수 있습니다(균일하게, 또는 2단의 느린 입구-출구 행진으로). 이 축약은 편의가 아니라, 제곱된 종횡비가 허가하는 것입니다. ∎

뿌리 사슬τ∼L2/D\tau\sim L^2/D (Fick)τ⊥/τ∥∼(L⊥/L∥)2\tau_\perp/\tau_\parallel\sim(L_\perp/L_\parallel)^2얇고 넓은 셀 ⇒\Rightarrow 비 ≪1\ll 1면 방향 동결 ⇒\Rightarrow 두께 방향만 분해
실전 예제· 얇고 넓은 셀의 두 확산 시간

작동 셀 하나를 잡습니다: 0.5 mm0.5~\mathrm{mm} 두께의 MEA 샌드위치가 10 cm10~\mathrm{cm} 폭의 활성 면적에 펼쳐져 있고, 대표적인 용존 화학종 확산도는 D=1×10−9 m2/sD = 1\times10^{-9}~\mathrm{m^2/s}입니다. 두 확산 시간과 그 비를 계산합니다.

양기호값
두께 방향 두께L⊥L_\perp5×10−4 m5\times10^{-4}~\mathrm{m}
면 방향 폭L∥L_\parallel0.1 m0.1~\mathrm{m}
확산도(양방향)DD1×10−9 m2/s1\times10^{-9}~\mathrm{m^2/s}
두께 방향 시간τ⊥=L⊥2/D\tau_\perp = L_\perp^2/D2.5×102 s2.5\times10^{2}~\mathrm{s}
면 방향 시간τ∥=L∥2/D\tau_\parallel = L_\parallel^2/D1×107 s1\times10^{7}~\mathrm{s}
비τ⊥/τ∥\tau_\perp/\tau_\parallel2.5×10−52.5\times10^{-5}

두께 방향은 약 4분에 정착하고, 면 방향은 여러 달이 걸립니다. 그 비 2.5×10−52.5\times10^{-5}은 정확히 (L⊥/L∥)2=(5×10−3)2(L_\perp/L_\parallel)^2 = (5\times10^{-3})^2이고, 1보다 다섯 자릿수 작으며, DD가 상쇄되므로 확산도와 무관합니다. 따라서 여기서 두께 방향 1-D 모델은 완전히 허가됩니다: 분해할 값어치가 있는 어떤 과도 응답도 면 방향 분포가 뒤척이기 한참 전에 끝납니다.

오르기 규칙

세 규칙이 오르기를 정직하게 지킵니다 :

  1. 물리보다 파라미터가 먼저 썩는다. GDL 투과도를 찍어 넣은 3단 모델은 측정된 파라미터의 1단 모델만 못합니다(20. 파라미터는 어디서 오는가).
  2. 한 단 아래로 검증하라. 모든 모델은 극한에서 아래 단을 재현해야 합니다. 채널 평균낸 1+1-D는 자기가 자라난 분극 방정식을 되찾아야 합니다.
  3. 단은 질문이 고른다. 가진 소프트웨어가 고르게 두지 마세요.

연습문제

Q1

어떤 소형 시험 셀은 단면이 거의 정사각형입니다: 활성 폭이 L∥=4 mmL_\parallel = 4~\mathrm{mm}, 두께 방향 총 스택이 L⊥=1 mmL_\perp = 1~\mathrm{mm}. 양방향에 같은 확산도를 써서 τ⊥/τ∥\tau_\perp/\tau_\parallel을 계산하고, 두께 방향 1-D 모델이 여전히 안전한지 판단하세요.

풀이 보기

확산도가 상쇄되므로 비는 순수 기하입니다:

τ⊥τ∥∼(L⊥L∥)2=(14)2=0.0625\frac{\tau_\perp}{\tau_\parallel} \sim \left(\frac{L_\perp}{L_\parallel}\right)^2 = \left(\frac{1}{4}\right)^2 = 0.0625

약 6 %입니다. 두께 방향이 면 방향보다 겨우 열여섯 배쯤 빠를 뿐이라, 면 방향 변화는 두께 방향 평형 시간에 걸쳐 더는 무시할 수 없습니다. 여기서 1-D 모델은 아슬아슬합니다: 면 방향 균일성에 민감한 질문(전류 분포, 국소적으로 뜨겁거나 잠긴 자리)은 2단이나 3단이 필요합니다. 폭을 두께 쪽으로 줄이면 1-D의 유일한 허가증인 제곱된 종횡비가 증발합니다.

Q2

각 질문을 충분한 가장 낮은 단에 배정하고, 한 구절로 이유를 대세요: (a) 시스템 연구를 위한 100셀 스택의 정격 출력, (b) 채널의 출구가 잠기는 동안 입구가 마르는지 여부, (c) 리브(랜드) 아래 대 열린 채널 위에서 전류밀도가 몰리는 방식.

풀이 보기

(a) 0단(셀당 집중 V(j)V(j)). 스택을 둘러싼 매니폴드·냉각·제어기 네트워크에 대해서만 4단으로 승격됩니다. (b) 2단, 입구-출구 조성을 분해하는 채널 방향 +1+1-D 행진. (c) 3단, 3-D CFD: 리브 대 채널의 몰림은 밀리미터 규모의 면 방향, 기하가 주도하는 변화이고 리브 폭과 MEA 두께가 서로 비슷합니다. 바로 종횡비 비가 1 규모라 두께 방향 1-D 허가증이 깨지는 영역입니다. 각 답은 질문이 지목하는 물리를 아직 담고 있는 가장 작은 단입니다.

이 장이 다루지 못하는 것

열화 모델(자체의 시간축 사다리)과, 2단과 3단에서 훈련된 뒤 그 단들을 점점 단락시키는 데이터 기반 서로게이트.

참고문헌

[1] S. G. Marquis, V. Sulzer, R. Timms, C. P. Please, S. J. Chapman (2019). An Asymptotic Derivation of a Single Particle Model with Electrolyte J. Electrochem. Soc. 166, A3693. doi:10.1149/2.0341915jes ↗

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

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

함께 보기

8. P2D 지배방정식 · 14. 분극 곡선 읽기 · 20. 파라미터는 어디서 오는가 · O'Hayre 외 6장; Weber & Newman, Chem. Rev. 104 (2004) 4679.

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

jinyongkim@hanyang.ac.kr
https://real.hanyang.ac.kr

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