PEMFC 燃料電池 分極曲線シミュレーター 戻る
エネルギー・水素

PEMFC 燃料電池 分極曲線シミュレーター

固体高分子形燃料電池(PEMFC)の V-i 分極曲線を、活性化・オーミック・濃度の 3 つの過電圧に分けて可視化します。セル温度・動作圧力・交換電流密度・ターフェル勾配・内部抵抗・限界電流密度を変えると、動作点のセル電圧/出力密度/効率/水素消費量がリアルタイムに変化し、Mirai や NEXO のような実用 PEM スタックの設計トレードオフを体験できます。

パラメータ設定
セル温度 T
°C
Nernst 電位と i₀ に影響。70〜85°C が実機の常用域
動作圧力 P
atm
加圧で Nernst 項が上昇するが、コンプレッサ動力で相殺される
交換電流密度 i₀
A/cm²
Pt 触媒の活性指標。Pt 量や合金化で改善
ターフェル勾配 b
V/dec
ORR 反応機構を表す勾配。Pt/C で 60〜70 mV/dec
内部抵抗 R_ohm
Ω·cm²
Nafion 膜厚・含水率・接触抵抗の合計
限界電流密度 i_L
A/cm²
GDL の O₂ 拡散律速で決まる上限
セル有効面積 A
cm²
単セルあたりの MEA 面積(Mirai は約 273 cm²)
計算結果
動作電流密度 (A/cm²)
セル電圧 (V)
出力密度 (W/cm²)
セル出力 (kW)
セル効率 (%)
水素消費 (g/h)
PEMFC 単セル断面 — H₂/O₂ 流入と H⁺・e⁻ の移動

左:アノード(H₂ → 2H⁺ + 2e⁻)、中央:Nafion 膜(H⁺ のみ通過)、右:カソード(½O₂ + 2H⁺ + 2e⁻ → H₂O)。e⁻ は外部回路を流れて負荷で仕事をします。

V-i 分極曲線(3 損失分離)
出力密度 vs 電流密度(最大出力点)
理論・主要公式

$$V = E_{rev} - b\,\log\!\left(\frac{i}{i_0}\right) - i\,R_{ohm} + b\,\log\!\left(1 - \frac{i}{i_L}\right)$$

PEMFC の分極式。b=ターフェル勾配、i₀=交換電流密度、R_ohm=膜・触媒の合計抵抗、i_L=限界電流密度。第 2 項が活性化、第 3 項がオーミック、第 4 項が濃度分極。

$$E_{rev} = 1.229 - 8.5\!\times\!10^{-4}(T-298) + \frac{RT}{2F}\cdot 0.5\,\ln P$$

Nernst 電位の温度・圧力補正(H₂/O₂ 反応の簡易形)。R=8.314 J/(mol·K)、F=96485 C/mol、T は絶対温度。

$$\eta_{cell} = \frac{V_{cell}}{1.481},\qquad \dot m_{H_2} = \frac{i\cdot A}{2F}\cdot M_{H_2}$$

セル効率(HHV 基準、1.481V=H₂ の発熱量電圧)と水素消費速度。M_{H₂}=2.016 g/mol。化学量論からセル電流に対する H₂ 流量が直接出る。

PEMFC 燃料電池の分極曲線設計

🙋
「燃料電池って水素を燃やしてるんですよね?普通のエンジンと何が違うんですか?」
🎓
いや、燃やしてないんだ。これがミソでね、PEMFC は H₂ と O₂ を電気化学的に反応させる装置で、熱を介さないからカルノー効率の壁を越えられる。アノードで H₂ → 2H⁺ + 2e⁻ に分かれて、H⁺ は Nafion 膜の中をスルッと抜け、e⁻ は外の回路をぐるっと回って負荷を駆動して、カソードで O₂ と再会して水になる。理論電圧は 1.229V、実用は 0.6〜0.7V/セル。スタックで何百枚も直列にして 300V くらい作って、Mirai みたいな FCV を走らせるわけ。
🙋
なるほど…でも左で電流密度を 0.5·i_L にしてるグラフを見ると、電圧がストンと 0.78V くらいまで落ちますよね。1.229V から 0.45V も損してる。これって何がそんなに食ってるんですか?
🎓
いいところに気づいた。分極曲線で電圧が落ちる原因は 3 つに分けられる。1 つ目が活性化分極 η_act = b·log(i/i₀) = 0.06·log(1/1e-4) ≒ 0.24V。これが一番大きい。カソードの ORR(酸素還元反応)はそもそも遅くて、Pt 触媒で叩いてもこれだけ食う。2 つ目がオーミック分極 η_ohm = i·R = 1·0.15 = 0.15V。膜と触媒層の電気抵抗で、Nafion 膜を薄くしたり加湿で含水率を上げると下がる。3 つ目が濃度分極 η_conc = 0.018V。i = i_L/2 ならまだ小さいけど、i_L に近づくと一気に効いてくる。
🙋
じゃあ Pt をたくさん使って i₀ を上げれば電圧が上がるってことですか?トヨタが Pt 量を減らすのに苦労してるって聞いた気がするんですけど…
🎓
そう、まさにそれが業界の大テーマ。Pt は 1g 3,000〜4,000 円するから、1 台 30g 積むと触媒だけで 10 万円飛ぶ。i₀ を 10 倍にできれば η_act は 0.06·log(10) = 60 mV 改善するんだけど、Pt を 10 倍にしても i₀ は √比例くらいでしか上がらない。だから Pt-Co、Pt-Ni 合金や、Pt 粒子のナノ構造制御、コアシェル構造で「Pt 量はそのまま、表面積と活性だけ上げる」のが王道。Mirai 2 では初代の半分以下の Pt 量で同じ出力を達成している。
🙋
出力密度のグラフを見ると、最大出力って i_L のだいぶ手前にありますね。なんで i_L まで電流を引かないんですか?
🎓
P = V·i だから、i を増やすと V が下がる戦いになる。i_L 直前は V がほぼゼロまで落ちるから、出力も落ちて下に凸の山になる。だいたい i = 0.5〜0.7·i_L 付近に最大出力点 P_max が来る。実用上はもっと大事な視点があって、効率が良いのは低電流側(V が高い)。FCV では加速時は最大出力点(V=0.6V 付近)、巡航時は高効率点(V=0.7V 付近)と動作点を切り替える。下のグラフで動作点が i_L の半分にあるのは「効率と出力のバランス点」を狙ってるんだ。
🙋
最後に、温度を上げると効率が上がるって本当ですか?でも 95°C を超えたらどうなるんでしょう…
🎓
「ある温度までは」上がる。i₀ がアレニウス的に増えて η_act が下がる効果と、膜のプロトン伝導率が上がって R_ohm が下がる効果のダブルで、60→80°C で出力 20% 改善とかザラ。でも 90°C を超えると Nafion 膜の含水率が落ちて急にカラカラになり、逆に R_ohm が爆増する。だから 80°C 前後を選ぶのが定石で、加湿器とラジエータの設計が燃料電池車の心臓部になる。SOFC(固体酸化物)みたいに 700°C で動かす別系統もあるけど、それは定置用で別物だね。

よくある質問

PEM 燃料電池の電圧降下は (1) 活性化分極(電極反応の活性化エネルギーに起因し、低電流密度で支配的)、(2) オーミック分極(電解質膜・触媒層・接触抵抗による IR 降下で、中電流密度域で線形に効く)、(3) 濃度分極(拡散律速で限界電流密度 i_L に近づくと急峻に立ち上がる)の 3 つです。本ツールは V = E_rev − η_act − η_ohm − η_conc で各損失を分離表示し、どの設計改善(触媒・膜・GDL)がどこを下げるかを直感的に確認できます。
理論起電力(HHV 基準)は 1.481V、Nernst 電位(LHV 基準)は標準条件で 1.229V ですが、活性化・オーミック・濃度損失で実用電圧は 0.6〜0.7V/セルに落ちます。電圧を高く保つ(=効率を高くする)ほど取り出せる電流が小さくなり、逆に低電圧で運転するほど大電流で高出力が得られます。Toyota Mirai 等の FCV では、加速時に 0.6V 付近の高出力点、巡航時に 0.7V 付近の高効率点を使い分け、平均 50〜55% のスタック効率を実現しています。
i_L はガス拡散層(GDL)と触媒層を通って反応サイトに O₂/H₂ が供給できる最大流束で決まります。i_L に近づくと濃度分極 η_conc = −b·log(1 − i/i_L) が対数発散して電圧が急落し、i = i_L では理論上 0V になります。実機では生成水のフラッディングで GDL が詰まり i_L が低下する現象も発生します。本ツールは i_op = 0.5·i_L を動作点としているため、設計余裕として i_op が i_L の半分以下に収まるかを確認できます。
温度上昇は Nernst 電位を −0.85 mV/K で下げる一方、交換電流密度 i₀ とプロトン伝導度を指数関数的に上げるため、トータルでは 60→80°C で出力が大きく改善します。ただし 90°C を超えると Nafion 膜が乾燥し、加湿器設計が難しくなります。圧力上昇は Nernst 項を (RT/4F)·ln(P_O₂·P_H₂²) で押し上げ、1→3 atm で約 30 mV/セル改善しますが、コンプレッサ動力で効率が打ち消されるため、車載 PEMFC では 1.5〜2.5 atm が最適です。

実世界での応用

燃料電池車(FCV):Toyota Mirai(2 代目)は 330 セル、約 128 kW、Hyundai NEXO は 440 セル・95 kW のスタックを搭載。3 分の水素充填で約 700 km 走行可能で、BEV の充電時間問題と航続距離問題を一気に解決します。本ツールで A=250、i_L=2.0 にして電圧 0.65V、出力密度 0.8〜1.0 W/cm² を狙うのが、量産 FCV の典型的な設計範囲です。

商用バス・トラック・物流:長距離トラックは BEV だと電池重量が積載量を圧迫するため、Nikola・Hyundai XCIENT・Toyota Hino などが PEMFC を採用。総出力 150〜300 kW、複数スタック並列構成が標準で、ディーゼル並みの稼働時間と急速充填を両立します。フォークリフトでも Plug Power が累計 7 万台以上の PEM 燃料電池を出荷しています。

定置式コージェネ(エネファーム):家庭用 PEMFC(パナソニック・東芝)は 700W 級で、都市ガスを改質して H₂ を作り、発電+給湯で総合効率 95% を実現。30 万台以上が日本で稼働中。本ツールで A を 100〜150 cm² 程度に下げ、低電流密度で高効率点を狙う設計が定置用の典型です。

ドローン・無人航空機(UAV)・宇宙:PEMFC はリチウムイオン電池の 3〜5 倍のエネルギー密度を持ち、長時間飛行用ドローン(IFC・H2Fly)に採用されています。Apollo・Space Shuttle の電源も PEMFC で、生成水を飲料水として再利用していました。本ツールで A を小さく、効率を最大化する低 i_op の設計が宇宙用途のセオリーです。

よくある誤解と注意点

最初の落とし穴は、「燃料電池の効率=セル電圧 ÷ 1.229V」と思い込むことです。1.229V は Nernst 電位(ΔG / 2F、LHV 基準)で、実際の燃料効率を測るときは HHV 基準の 1.481V で割ります。本ツールは ηcell = V/1.481 を採用しており、これは「H₂ の発熱量のうち何%が電気になったか」を表します。さらに実機ではコンプレッサ・ポンプ・冷却ファンの BOP(Balance of Plant)動力が 10〜15% 取られるため、本ツールの systemEfficiency は cellEfficiency × 0.85 で近似しています。カタログの「効率 60%」が LHV か HHV か、スタック単体か BOP 込みかは必ず確認してください。

2 つ目は、「Pt をたくさん使えば性能が上がる」という誤解。i₀ ∝ Pt 表面積で、Pt 重量を 2 倍にしても表面積は √2 倍程度しか増えないため、活性化分極の改善は b·log(√2) ≈ 9 mV と微々たるものです。一方で Pt は 1g 3,000〜4,000 円と高価で、Mirai 1 世代の Pt 量約 30g を半減させることが商用化の絶対条件でした。現代の PEMFC 設計は「Pt 量を下げる」「Pt-Co/Pt-Ni 合金で活性を上げる」「カーボン担体を最適化して Pt 粒子をナノ分散させる」の 3 軸で攻めます。本ツールの i₀ を 1e-5 → 1e-4 と動かすと、η_act が 0.30 → 0.24V と 60 mV しか変わらないことが体感できます。

3 つ目は、「水管理を無視した設計」。PEMFC は反応で水を作るため、高電流密度ではカソード GDL が水で詰まる「フラッディング」が起き、見かけの i_L が大きく下がります。逆に低負荷・高温では Nafion 膜が乾いて R_ohm が爆増。本ツールの i_L は設計値ですが、実機ではフラッディング監視のためにセル電圧をリアルタイムで分布測定し、加湿度・流量・温度を多重制御しています。「分極曲線がきれいに出る = 実機が動く」ではなく、過渡応答・低温始動(−30°C)・凍結融解耐久まで設計検証して初めて FCV になります。

使い方ガイド

  1. セル温度(60~80°C)と作動圧力(1~3 atm)を設定し、PEMFC の熱力学的動作点を決定
  2. 交換電流密度(0.1~10 mA/cm²)とターフェル勾配(0.04~0.12 V/decade)を調整して、触媒活性と電極反応速度を模擬
  3. 内部抵抗と限界電流密度を入力し、分極曲線を計算。セル電圧・出力密度・効率・水素消費量をリアルタイム表示

具体的な計算例

Toyota Mirai 仕様を想定:セル温度80°C、圧力2.5 atm、交換電流密度2.0 mA/cm²、ターフェル勾配0.065 V/decade、内部抵抗0.25 Ω・cm²、限界電流密度1500 mA/cm²の場合、動作電流密度500 mA/cm²でセル電圧0.68 V、出力密度340 W/cm²、セル効率52%、水素消費2.8 g/h となります。

実務での注意点