PEMFC 水管理シミュレーター 戻る
燃料電池・水管理

PEMFC 水管理シミュレーター

固体高分子形燃料電池(PEMFC)の生命線である「水バランス」を可視化するツールです。電流密度・温度・圧力・膜厚・加湿条件を変えると、反応水生成・電気浸透抗力・空気持出量から膜含水率λとプロトン伝導度σが Springer モデルで即時計算され、フラッディング/ドライアウト境界が見えます。

パラメータ設定
セル有効面積 A
cm²
電流密度 i
A/cm²
セル温度 T
°C
セル圧力 P
bar
膜厚(Nafion グレード)
膜厚と乾燥時の λ を選択
空気量論比 λ_air
必要 O₂ に対する空気供給倍率
アノード加湿 RH_an
%
カソード加湿 RH_ca
%
計算結果
反応水生成 (g/s)
電気浸透抗力 λ_drag
カソード λ
プロトン伝導度 σ (S/m)
膜抵抗 R (Ω·cm²)
水バランス状態
PEMFC セル断面 — 水の流れアニメーション

アノード(左)→ 膜 → カソード(右)の断面。緑:H⁺、青:H₂O。電気浸透抗力で水がカソードへ運ばれ、生成水と合わせて空気で排出されます。膜の色が含水率λを表します。

プロトン伝導度 σ vs 含水率 λ
水フラックス収支(mol/s ×10⁻³)
理論・主要公式

$$\sigma(\lambda) = (0.005139\lambda - 0.00326)\exp\left[1268\left(\frac{1}{303}-\frac{1}{T}\right)\right]$$

λ は含水率(mol H₂O / mol SO₃⁻、範囲 0–22)。Springer (1991) の Nafion 標準モデルで、σ の単位は S/cm。

$$\lambda(a) = 0.043 + 17.81\,a - 39.85\,a^{2} + 36\,a^{3}$$

水蒸気活量 a(相対湿度/100)から膜の含水率 λ を求める Springer 等温線。a=1 で λ≈14。

$$\dot{n}_{\text{drag}} = \frac{I\,\lambda_{\text{drag}}}{F}, \qquad \lambda_{\text{drag}} = \frac{2.5\,\lambda_{\text{dry}}}{22}$$

電気浸透抗力によるアノード→カソードへの水フラックス。F は Faraday 定数、I は総電流。

PEMFC 水管理 — 加湿・ドライアウト・フラッディング

🙋
PEMFC(固体高分子形燃料電池)って、燃料の H₂ と空気の O₂ から水と電気を作るんですよね?水ができるなら、わざわざ「加湿」って必要なんですか?
🎓
いい質問だね。鍵は真ん中にある「Nafion」っていう高分子膜にあるんだ。この膜は中に水分子があってはじめてプロトン(H⁺)を運べる。乾いた膜は電気を通さない絶縁体に近くて、抵抗が10倍にも跳ね上がる。だから反応で水ができる前の起動直後や、低負荷で水が足りない時は、入口の空気・水素を加湿器でわざわざ湿らせて膜を「水浸し」に保つんだ。Toyota Mirai みたいな実車では、この加湿器のコンパクト化が長年の開発テーマだった。
🙋
なるほど、膜は濡らしておかないとダメと。じゃあ電流密度を上げて反応水がドバドバ出てくれば、加湿は要らなくなる?
🎓
理屈ではそうだけど、今度は逆の問題が出てくる。電流密度が 1 A/cm² を超えると、反応水+電気浸透で運ばれた水が、カソードのガス流路や触媒層に液水で溜まり始める。これが「フラッディング」。空気が触媒に届かなくなって、電圧がストンと落ちる。左で電流密度を上げてみて。水バランスがプラスに振り切れて「フラッディング」判定になるはず。
🙋
本当だ、1.0 A/cm² でもう「フラッディング」って出てます。ということは、空気量論比 λ_air を上げて乾かせばいいんですよね?
🎓
そう、それが基本対策。λ_air を 2.0 から 3.0 に上げると空気モル流量が 1.5 倍になって、出口で持ち去る水蒸気が増える。ただし、ブロワーの消費電力が λ_air に対してほぼ3乗で効いてくるから、ずっと λ_air=3 で回すとシステム効率が落ちる。実機では「普段は λ_air=1.8〜2.0、フラッディングを感知したら一時的に 3 まで上げてパージする」という運転制御をするんだ。Plug Power や Ballard の据置電源、Daimler GLC F-CELL あたりはみんなこの考え方。
🙋
逆にドライアウト側だと、何が起こるんですか?膜抵抗が増えるだけ?
🎓
「だけ」じゃ済まない。膜抵抗が増えると IR 損失で発熱が増え、もっと乾く負のフィードバックが起きる。さらに乾いた状態で機械応力が掛かると、Nafion にピンホールが空いて H₂ がクロスオーバーする — これは膜の不可逆損傷で、フラッディングと違って戻らない。だから設計では「やや湿った側」にマージンを取るのが鉄則。膜厚も N212(50μm)みたいに薄い膜は応答が速い反面、ドライアウトに弱い。Mirai 2nd gen が N115 系の比較的厚い膜を使うのもこの保険の意味があるんだ。
🙋
膜の含水率 λ ってカソードで 3.5 とか出てますけど、これって絶対値で何を意味してるんですか?
🎓
λ は「SO₃⁻ 1モルあたりに H₂O が何モルくっついているか」で、Nafion では 0〜22 の範囲を取る。λ=14 が水蒸気飽和(RH 100%)、λ=22 が液水接触の上限。λ<3 だとプロトンが動けなくなって伝導度がほぼゼロ、λ=10〜14 が実用域。今のカソード λ=3.5 はかなり乾いている側で、もしアノードを 100%、カソードを 70% くらいまで湿らせれば λ_avg は 8〜10 に上がって伝導度も2倍以上になる。スライダーで試してみて。

よくある質問

セル内の水収支は「生成量+アノードから運ばれる量−カソード出口で空気が持ち出す量」で計算します。反応水は n_H2O = I/(2F)、電気浸透抗力は n_drag = I·λ_drag/F、空気持出は飽和水蒸気圧と空気モル流量から決まります。本ツールではこの三者を足し合わせて正なら「フラッディング寄り」、負なら「ドライアウト寄り」と判定します。
Springer (1991) が提案した Nafion 膜の経験モデルで、膜の含水率 λ(SO3⁻ 1モルあたりの H2O モル数、0〜22)から水蒸気活量とプロトン伝導度を関係付ける標準的なモデルです。σ(λ,T) = (0.005139λ − 0.00326)·exp[1268·(1/303 − 1/T)] で、λ が高いほど、また T が高いほど σ が大きくなります。本ツールはこの式で σ と膜抵抗 R を計算します。
どちらも性能を落としますが性質が違います。フラッディングはガス流路や触媒層に液水が溜まり、酸素到達が阻害されて高電流で急激に電圧が落ちます。ドライアウトは膜抵抗が跳ね上がり、全電流域で IR 損失が増え、ひどい場合は膜のピンホール/劣化に直結します。一般的にフラッディングは負荷を下げれば回復しますが、ドライアウト由来の膜損傷は不可逆なので、設計では「やや湿った側」に寄せることが多いです。
λ_air を上げるとカソード側の空気モル流量が増え、出口で持ち去る水蒸気も比例して増えます。つまり「乾かす」方向に働きます。λ_air=2 付近はバランス点としてよく使われ、λ_air を 3〜4 に上げるとフラッディング解消には効きますが、ブロワー消費電力が増えてシステム効率が落ちます。低負荷で詰まりやすい運転点では、一時的に λ_air を上げるパージ動作で対応するのが定石です。

実世界での応用

燃料電池自動車(FCV):Toyota Mirai(2nd gen, 134 kW)、Honda Clarity、Hyundai Nexo、Daimler GLC F-CELL などで、加湿器と水管理制御は車載パッケージの心臓部です。Mirai 2nd gen ではセルの3D流路と自己加湿性の高い膜の組合せで、外部加湿器の小型化/一部削減を実現しました。低温起動時の凍結対策(マイナス30℃始動)も水管理の延長線上にあります。

定置型コージェネ・バックアップ電源:Toshiba H2One、Bloom Energy(SOFC ですが類似課題)、Plug Power(フォークリフト電源)、Ballard(バス・物流車)、Cummins-Hydrogenics などが PEMFC を商品化しています。定置用は車載に比べ熱・水を回収しやすく、自己加湿に振った設計で外部加湿器を撤廃したシステムが主流になりつつあります。

自己加湿膜・先進材料:外部加湿器のコストとスペースをなくすため、Nafion 3M PFIA、Aquivion(Solvay)、Pt 微粒子分散による自己加湿膜などが実用化されています。空気中の H₂ クロスオーバーを利用してカソード側で水を生成させる仕組みで、システム simplification と低温運転性能の両方に効きます。

CAE による水分布シミュレーション:市販 CAE(COMSOL Multiphysics、ANSYS Fluent、STAR-CCM+)には PEMFC モジュールが用意され、ガス拡散層(GDL)内の液水分布や、流路内の水滴生成を3次元で解きます。本ツールのような0次元バランス計算は、その3次元 CAE モデルの境界条件決定や妥当性検証の起点として使われます。

よくある誤解と注意点

まず最大の誤解が、「100% RH で加湿すれば水管理は完璧」という思い込み。実は入口を 100% RH にすると、わずかな冷えで液水が凝縮してフラッディングの種になります。実機では「アノード 100%、カソード 50〜70%」のような非対称加湿(humid-dry 構成)が多用され、水蒸気活量勾配で意図的にバックディフュージョンを起こします。本ツールのデフォルトもこの考え方に従っています。

次に、「Springer モデルの λ_drag は定数 2.5」と教科書通りに使ってしまうこと。実際の電気浸透抗力係数は膜の含水状態に依存し、乾いた膜では 1 程度、飽和液水中では 3 を超えることもあります。本ツールは Springer 原典に従って λ_drag = 2.5·λ_dry/22 の簡易式を使っていますが、低 RH 運転や薄膜(N212)では実測値とずれることがあります。研究用途では Zawodzinski の温度依存モデルなど、より精緻な相関を用いてください。

最後に、「膜抵抗 R が小さければ性能は上がる」という単純化。R を下げるには薄膜化(N212→N115p の逆)が効きますが、薄膜は H₂ クロスオーバーが増え、ピンホール/劣化リスクも増します。さらに薄いほどドライアウト時の含水勾配が急になって、機械的疲労が早く現れます。膜抵抗、クロスオーバー、機械強度の三角形で最適点を探す必要があり、単純な「薄ければよい」は危険です。

使い方ガイド

  1. セル面積(cm²)と電流密度(A/cm²)を入力:例えば面積50cm²、電流密度0.8A/cm²の場合、総電流は40Aになります
  2. セル温度(℃)と加湿圧力(bar)を設定:PEMFC標準運転は60~80℃で、加湿圧力は3~4barが一般的です
  3. 計算ボタンを押すと、反応水生成量(g/s)、電気浸透抗力λ_drag、カソード含水率λ、プロトン伝導度σ(S/m)、膜抵抗R(Ω·cm²)、水バランス状態をリアルタイム出力します

具体的な計算例

セル面積100cm²、電流密度1.0A/cm²、温度70℃、加湿圧力3.5barの条件で:反応水生成は約0.0373g/s、電気浸透抗力λ_dragは2.8、カソード含水率λは14.2です。プロトン伝導度σは温度依存性を考慮して約0.085S/m、膜抵抗Rは約0.118Ω·cm²となり、水バランス状態は「適正」と判定されます。Nafion117膜(厚さ183μm)の場合、膜内抵抗は約1.47mΩ·cm²です。

実務での注意点