燃料電池の実効電圧を計算・可視化できるオンラインツールの紹介

目次

はじめに

燃料電池は、水素と酸素の化学反応を利用して電気を生み出すクリーンな電源で燃料電池は水素と酸素の電気化学反応によって電気を生み出すクリーンな電源です。しかし、理論上の最大電圧(理論電圧、25℃・1気圧で約1.23 V)と、実際に取り出せる電圧(実効電圧)は一致しません

この差を生み出す原因は主に3種類の損失です。活性化過電圧・内部抵抗によるオーム損・濃度過電圧がそれぞれ異なる電流密度領域で電圧を下げます。これらを個別に定量化することで、「どの損失が性能を制限しているか」を把握し、触媒・膜・ガス拡散層(GDL)の改善につなげることができます。

このページでは3種類の過電圧を考慮した実効電圧の計算方法と、パラメータを入力するだけでI-V特性曲線を自動描画できるツールを紹介します。

実効電圧の計算方法

燃料電池の実効電圧は次の式で表せます:Veff=VtheoηactηohmηconcV_{eff} = V_{theo} – \eta_{act} – \eta_{ohm} – \eta_{conc}

記号意味支配的な電流密度領域
VtheoV_{theo}理論電圧(≈1.23 V)
ηact\eta_{act}活性化過電圧低電流密度域
ηohm=IRint\eta_{ohm} = I \cdot R_{int}内部抵抗によるオーム損中電流密度域
ηconc\eta_{conc}濃度過電圧高電流密度域

各過電圧の計算方法

内部抵抗(IRドロップ)

内部抵抗 RintR_{int}(電解質膜・電極・接触抵抗の合計)による電圧降下は次の式です。ηohm=IRint\eta_{ohm} = I \cdot R_{int}

RintR_{int}は交流インピーダンス測定(EIS)または高周波抵抗計(HFR)で実測できます。IR補正後の電圧(VIR-free)は Vcell+IRintV_{cell} + I \cdot R_{int}で求められます。

活性化過電圧(ターフェル式)

低電流密度領域では電極反応の活性化エネルギー障壁による電圧降下が支配的です。ターフェル式で表されます。ηact=a+blogi\eta_{act} = a + b \log i

  • aa:切片(交換電流密度i₀に依存。i₀が大きいほど活性化過電圧が小さい)
  • bb:ターフェル勾配(反応機構・移動係数αcに依存。典型値:約120 mV/decade)

実験データからI-V特性のIR-free電圧を取得し、低電流域を片対数プロット(ターフェルプロット)して線形フィッティングすることでa・bを求めます。

詳細な測定・解析手順は→ 活性化過電圧とは?測定方法・ターフェル式による解析の基本をご覧ください。

濃度過電圧

高電流密度域では電極近傍でのガス拡散律速が生じ、電流密度が限界電流密度 iLi_LiL​ に近づくほど電圧が急落します。ηconc=cln(1iiL)\eta_{conc} = -c \cdot \ln\left(1 – \frac{i}{i_L}\right)

  • iLi_L:限界電流密度(ガス拡散係数・GDL厚さ・ガス濃度で決まる)
  • cc:係数(c=RT/nFc = RT/nF、25℃・n=4のとき約0.0065 V)

実験データの高電流域に対して非線形フィッティングを行い、iLi_LiL​とcを推定します。

詳細な解析手順は→ 濃度過電圧とは?測定方法・限界電流密度の求め方を解説をご覧ください。

ツールの使い方

入力項目の説明

入力項目説明
理論電圧 VtheoV_{theo}Vtheo​ [V]通常1.23 V(温度・圧力で変化)
内部抵抗 RintR_{int}Rint​ [Ω·cm²]HFR等で実測した値
ターフェル切片 a [V]低電流域フィッティングから取得
ターフェル勾配 b [V/decade]同上(典型値:0.06〜0.12)
限界電流密度 iLi_LiL​ [A/cm²]高電流域フィッティングから取得
係数 c [V]通常は理論値(RT/nFRT/nFRT/nF)を使用
電流密度の範囲計算・グラフ表示する範囲を指定

計算結果の見方

ツールでは入力パラメータをもとに以下が自動計算・グラフ表示されます。

  • I-V特性曲線: 実効電圧 VeffV_{eff}vs 電流密度 i のグラフ
  • 各過電圧の内訳: 活性化過電圧・オーム損・濃度過電圧それぞれの電流密度依存性を重ねて表示

パラメータを変えてシミュレーションすることで「iLi_LiL​を大きくするとI-V特性がどう改善するか」「RintR_{int}Rint​を下げるとどこが変わるか」といった設計感度の把握に活用できます。

🔧実効電圧の可視化ツール

計算例

条件: VtheoV_{theo}​=1.23 V、RintR_{int}=0.1 Ω·cm²、a=−0.4 V、b=0.08 V/decade、iLi_L=2.0 A/cm²、c=0.006 V

電流密度 i=1.0 A/cm² での各過電圧と実効電圧:ηact=0.4+0.08×log(1.0)=0.4 V\eta_{act} = -0.4 + 0.08 \times \log(1.0) = -0.4 \text{ V}ηohm=1.0×0.1=0.1 V\eta_{ohm} = 1.0 \times 0.1 = 0.1 \text{ V} ηconc=0.006×ln(11.02.0)0.004 V\eta_{conc} = -0.006 \times \ln\left(1 – \frac{1.0}{2.0}\right) \approx 0.004 \text{ V}Veff=1.230.40.10.004=0.726 VV_{eff} = 1.23 – 0.4 – 0.1 – 0.004 = 0.726 \text{ V}

この電流密度での出力密度は P=0.726×1.0=0.726P = 0.726 \times 1.0 = 0.726 W/cm²となります。

関連記事

よかったらシェアしてね!
  • URLをコピーしました!
  • URLをコピーしました!
目次