JP2010236989A - 粒径分布モデル作成方法、当該粒径分布モデル作成方法を用いた燃料電池触媒の劣化予測方法、及び当該燃料電池触媒の劣化予測方法を用いた燃料電池の制御方法 - Google Patents
粒径分布モデル作成方法、当該粒径分布モデル作成方法を用いた燃料電池触媒の劣化予測方法、及び当該燃料電池触媒の劣化予測方法を用いた燃料電池の制御方法 Download PDFInfo
- Publication number
- JP2010236989A JP2010236989A JP2009084428A JP2009084428A JP2010236989A JP 2010236989 A JP2010236989 A JP 2010236989A JP 2009084428 A JP2009084428 A JP 2009084428A JP 2009084428 A JP2009084428 A JP 2009084428A JP 2010236989 A JP2010236989 A JP 2010236989A
- Authority
- JP
- Japan
- Prior art keywords
- particle size
- fuel cell
- size distribution
- distribution model
- catalyst
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000002245 particle Substances 0.000 title claims abstract description 364
- 238000009826 distribution Methods 0.000 title claims abstract description 193
- 239000003054 catalyst Substances 0.000 title claims abstract description 178
- 239000000446 fuel Substances 0.000 title claims abstract description 136
- 238000000034 method Methods 0.000 title claims abstract description 111
- 230000015556 catabolic process Effects 0.000 title claims abstract description 15
- 238000006731 degradation reaction Methods 0.000 title claims abstract description 15
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 claims description 173
- 230000006866 deterioration Effects 0.000 claims description 85
- 229910052697 platinum Inorganic materials 0.000 claims description 85
- 230000007423 decrease Effects 0.000 claims description 29
- 230000010354 integration Effects 0.000 claims description 29
- 238000013178 mathematical model Methods 0.000 claims description 23
- 238000005259 measurement Methods 0.000 claims description 23
- 238000006243 chemical reaction Methods 0.000 claims description 17
- 238000004090 dissolution Methods 0.000 claims description 14
- 238000007254 oxidation reaction Methods 0.000 claims description 9
- MUMZUERVLWJKNR-UHFFFAOYSA-N oxoplatinum Chemical compound [Pt]=O MUMZUERVLWJKNR-UHFFFAOYSA-N 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 description 25
- 239000012528 membrane Substances 0.000 description 21
- 239000003792 electrolyte Substances 0.000 description 14
- 239000007789 gas Substances 0.000 description 14
- 230000008569 process Effects 0.000 description 13
- 238000010248 power generation Methods 0.000 description 12
- 230000001186 cumulative effect Effects 0.000 description 11
- 238000005315 distribution function Methods 0.000 description 11
- 238000012423 maintenance Methods 0.000 description 11
- 238000002474 experimental method Methods 0.000 description 9
- 239000005518 polymer electrolyte Substances 0.000 description 9
- 239000007787 solid Substances 0.000 description 9
- 238000004088 simulation Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000009792 diffusion process Methods 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 230000000052 comparative effect Effects 0.000 description 4
- 239000012495 reaction gas Substances 0.000 description 4
- 230000003247 decreasing effect Effects 0.000 description 3
- 238000010828 elution Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 239000007800 oxidant agent Substances 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 230000007704 transition Effects 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 230000003197 catalytic effect Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 239000002826 coolant Substances 0.000 description 2
- 239000000498 cooling water Substances 0.000 description 2
- 238000010030 laminating Methods 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 230000003647 oxidation Effects 0.000 description 2
- 230000001590 oxidative effect Effects 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000006555 catalytic reaction Methods 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 238000003411 electrode reaction Methods 0.000 description 1
- 238000005370 electroosmosis Methods 0.000 description 1
- 239000010419 fine particle Substances 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- GPRLSGONYQIRFK-UHFFFAOYSA-N hydron Chemical compound [H+] GPRLSGONYQIRFK-UHFFFAOYSA-N 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003475 lamination Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000004043 responsiveness Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000012798 spherical particle Substances 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04992—Processes for controlling fuel cells or fuel cell systems characterised by the implementation of mathematical or computational algorithms, e.g. feedback control loops, fuzzy logic, neural networks or artificial intelligence
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/02—Investigating particle size or size distribution
- G01N15/0266—Investigating particle size or size distribution with electrical classification
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F15/00—Digital computers in general; Data processing equipment in general
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C99/00—Subject matter not provided for in other groups of this subclass
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M4/00—Electrodes
- H01M4/86—Inert electrodes with catalytic activity, e.g. for fuel cells
- H01M4/90—Selection of catalytic material
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M4/00—Electrodes
- H01M4/86—Inert electrodes with catalytic activity, e.g. for fuel cells
- H01M4/90—Selection of catalytic material
- H01M4/9041—Metals or alloys
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04082—Arrangements for control of reactant parameters, e.g. pressure or concentration
- H01M8/04089—Arrangements for control of reactant parameters, e.g. pressure or concentration of gaseous reactants
- H01M8/04119—Arrangements for control of reactant parameters, e.g. pressure or concentration of gaseous reactants with simultaneous supply or evacuation of electrolyte; Humidifying or dehumidifying
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/10—Fuel cells with solid electrolytes
- H01M2008/1095—Fuel cells with polymeric electrolytes
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/30—Hydrogen technology
- Y02E60/50—Fuel cells
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Chemical & Material Sciences (AREA)
- Electrochemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Materials Engineering (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- Manufacturing & Machinery (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Computer Hardware Design (AREA)
- Immunology (AREA)
- Dispersion Chemistry (AREA)
- Pathology (AREA)
- Fuel Cell (AREA)
- Mathematical Optimization (AREA)
- Medical Informatics (AREA)
- Databases & Information Systems (AREA)
- Computational Mathematics (AREA)
- Artificial Intelligence (AREA)
- Automation & Control Theory (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Fuzzy Systems (AREA)
- Algebra (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
Abstract
【解決手段】粒径範囲決定工程、粒径範囲決定工程で決定された粒径範囲における粒子の出現頻度を積分する積分工程、積分工程で得られた出現頻度の積分値を用いて、区分点となった粒径を決定する区分点決定工程、並びに、最小粒径、最大粒径及び各区分点の粒径を代表点として決定する代表点決定工程、を含み、各代表点の粒径を有する粒子を含有し、且つ、各代表点の粒径を有する粒子の出現頻度が当該代表点で区分された個々の領域の積分値と等しい粒径分布を仮想し、当該仮想された粒径分布を、粒径分布モデルとして得ることを特徴とする、粒径分布モデル作成方法。
【選択図】図6
Description
H2 → 2H+ + 2e− (I)
(I)式で生じる電子は、外部回路を経由し、外部の負荷で仕事をした後、カソード(酸化剤極)に到達する。そして、(I)式で生じたプロトンは、水和した状態で、固体高分子電解質膜内をアノード側からカソード側に、電気浸透により移動する。
2H+ + (1/2)O2 + 2e− → H2O (II)
カソードで生成した水は、主としてガス拡散層を通り、外部へと排出される。このように、燃料電池では、水以外の排出物がなく、クリーンな発電装置である。
本発明は、上記実情を鑑みて成し遂げられたものであり、より計算時間が短く、且つ、より精確な粒径分布モデル作成方法、当該粒径分布モデル作成方法を用いた燃料電池触媒の劣化予測方法、及び当該燃料電池触媒の劣化予測方法を用いた燃料電池の制御方法を提供することを目的とする。
論文1中には、膜・電極接合体中のカソード電極中の白金触媒の実測粒径分布を示したグラフ(Figure7)、及び、白金触媒の表面積を実測した結果を示したグラフ(Figure4)が開示されている。
論文1中のFigure7は、横軸に触媒粒径を、縦軸に粒子数を取り、初期状態の白金触媒粒子と、発電後の白金触媒粒子について、それぞれ200個の球状の粒子を選び、その粒径をすべて透過型電子顕微鏡(TEM)観察によって測定した結果得られたグラフである。発電後の白金触媒粒子(平均粒径=5.9nm)の粒径分布は、初期状態の白金触媒粒子(平均粒径=2.8nm)の粒径分布と比較して、分布の山が低く、且つ幅が広い、いわゆるブロードの状態となっている。これは、運転時間を経るにつれて白金触媒の溶解が進行した結果、粒径が比較的小さい粒子は溶解によって粒径がさらに小さくなり、且つ、溶け出した白金が、粒径が比較的大きい粒子と凝集することによって、粒径がさらに大きくなるために生じる現象である。
論文1中のFigure4は、横軸に発電サイクル数を、縦軸に触媒粒子の総表面積を取ったグラフである。上述したFigure7からも分かるように、白金触媒の溶解が生じることによって、溶け出した白金と、粒径が比較的大きい粒子との凝集が生じることから、発電サイクル数を重ねるごとに触媒粒子の総表面積が減少する。触媒粒子の総表面積が減少すると、燃料電池の電極反応を司る触媒反応場が減少することから、燃料電池触媒の劣化につながる。
図7は、論文2における白金触媒粒径分布を模式的に示したグラフであり、横軸に粒径を、縦軸に白金原子数を取ったグラフである。図から分かるように、論文2における白金触媒粒径分布モデルは、粒子数の変化を考慮しない2種類の異なる粒径の白金粒子で代表させている。
論文2中には、横軸に運転時間を、縦軸に触媒粒子の総表面積を取ったグラフ(Figure3)が開示されている。論文2中のFigure3から分かるように、上記図7に示したモデルを用いたシミュレーション結果を示した2つのグラフ(論文2中のFigure3中の実線で示されたグラフ)は、実験から得られた結果を示したグラフ(論文2中のFigure3中の破線で示されたグラフ)とはいずれも一致していない。これは、粒子数の変化を考慮しない2種類の異なる粒径の白金粒子で代表させた粒径分布モデルを用いているため、白金粒子表面積の変化が不連続となるからである。
図8(b)は、横軸に運転時間を、縦軸に触媒粒子の総表面積を取ったグラフである。図8(b)から分かるように、上記図8(a)に示したモデルを用いたシミュレーション結果を示したグラフ(実線)は、実験から得られた結果を示したグラフ(破線)とはいずれも一致していない。したがって、従来技術の単なる応用に過ぎない粒径分布モデルでは、実験で得られた結果をシミュレーションすることは困難であることが分かる。
このように、従来の方法においては、燃料電池の運転時間が長時間にわたる際の、実験で観察される白金触媒表面積の指数関数的な低下を再現することができなかった。
本発明は、従来技術のような粒径分布モデル作成方法とは異なり、異なる粒子の粒径の差を変数として用いることによって、より計算時間が短く、且つ、精確な粒径分布モデル作成方法を提供することを目的とする。
出現頻度を粒子数で表現した場合には、上述した論文1中のFigure7に示されるように、粒径分布は左右対称の分布(好ましくは正規分布)となる。
一方、出現頻度を質量又は体積で表現した場合には、粒径が大きくなるほど一般的に質量又は体積も大きくなることから、粒径分布は左右対称の正規分布ではなく、非対称の分布となってしまう。しかし、このような粒径分布は、運転時間初期においては実験値と若干の誤差を生じるものの、長時間の運転をシミュレーションする場合には、まったく問題とならない。
まず初めに行うべき粒径範囲決定工程とは、模擬すべき粒子群に含まれる粒子の最小粒径及び最大粒径を決定する工程のことである。この工程においては、最小粒径及び最大粒径と共に、平均粒径を決定してもよい。
図1(a)は、本発明の区分決定工程を具体的に示した模式図である。図1(a)は、縦軸に出現頻度f(x)を、横軸に粒径xをとったグラフであり、曲線20は、実際の粒径分布(正規分布)を示した曲線である。上述した粒径範囲決定工程において、模擬すべき粒子群に含まれる粒子の最小粒径xmin及び最大粒径xmaxは既に決定している。また、上述した積分工程において、粒子の出現頻度を、最小粒径xmin及び最大粒径xmaxを起点及び終点として積分しており、既に積分値(曲線20と横軸に挟まれた面積)は決定している。
図1(a)においては、前記積分値を用いて、図中の面積a=面積b=面積c=…=面積i=面積jとなるように、第1の区分を行うことによって、積分領域を10の領域に区分した様子を示している。これにより、区分された個々の領域の積分値が均等になるような区分点x1、x2、…x9が決定する。
最後の代表点決定工程においては、これら区分点x1、x2、…x9の粒径、並びに、最小粒径xmin及び最大粒径xmaxを代表点として決定する。
以下、正規分布の例について示す。正規分布は、下記式(1)に示すような確率密度関数を持つ。なお、下記式(1)中、xは確率変数を、μは平均値を、σは分散を、それぞれ示す。
Y1、Y2、…Ynの値は、例えば、以下のように決定することができる。図2は、累積分布関数F(X)を模式的に示したグラフである。Y1、Y2、…Ynを、F(X)軸上に等間隔に取ることによって、これらY1、Y2、…Ynに対応する、X1、X2、…Xnを得ることができる。
このような操作で得られたX1、X2、…Xnを、区分点として採用することができる。
特に、長時間にわたる粒径分布の推移を調べるための粒径分布モデルを作成する方法として有効であるという観点から、本発明の粒径分布モデル作成方法の一形態としては、第2の区分を、第1の区分によって区分された所定数の領域の内、最大粒径を含む領域について行うという構成をとることができる。
その理由は以下の通りである。上述した論文1中のFigure7において説明したように、粒径が比較的小さい粒子は、ごく短時間で溶け出し、他の粒子と凝集してしまうことから、非常に長い時間にわたる粒径分布の変動を予測する際には、それらの挙動を特に細かく追跡する必要はない。これに対し、粒径が比較的大きい粒子は、非常に長い時間にわたる粒径分布の変動を予測する際においては、その挙動の追跡が特に重要になる。しかし、上述した図1(b)に示したように、特に最大粒径を含む領域については、代表点が1点のみしか存在せず、細かい粒径分布変動の追跡には必ずしも十分とはいえない。したがって、最大粒径を含む領域(図1(a)中においては面積jを有する領域)について、区分された個々の領域の積分値が均等になるように所定数の領域にさらに区分することによって、長時間にわたる粒径分布の推移を調べるための粒径分布モデルを得ることができる。
このように、最小粒径、最大粒径及び出現頻度の積分値を均等に区分する各点を代表点として用いるため、本発明に係る粒径分布モデル作成方法により作成された粒径分布モデルは、従来技術の粒径分布モデルよりも、より少ない変数で、実験で得られる粒径分布を精確にシミュレーションすることができる。
なお、粒径分布は、精確に実際の粒径分布を再現することができれば、正規分布に限定されることはなく、様々な確率分布を用いることができる。確率分布には、離散型、絶対離散型、連続型等の種類があるが、本発明は、燃料電池触媒の劣化予測を主な目的の1つとし、現実の触媒粒子の分布は連続的な分布であることから、連続型の確率分布を用いることが好ましい。
連続型の確率分布としては、対数正規分布、指数分布、t分布、カイ2乗分布、ガンマ分布、ベータ分布、F分布、コーシー分布、アーラン分布、三角分布、ラプラス分布、レイリー分布、ロジスティック分布、パレート分布、ワイブル分布、実測した白金粒径分布を参考にした関数等を例示することができる。
まず、粒径分布モデル作成方法を用いて、燃料電池触媒の粒径分布モデルを作成し、さらに、粒径分布モデル及び上記数式モデルを用いて、初期触媒粒径、触媒層構造、これまでの運転パターン及び運転時間等の物理的・化学的データから、現在の白金触媒粒子の表面積を推定する。
図3に示された運転時間初期においては、白金触媒粒子の表面積は急激に低下する。この段階では、燃費向上よりも、電池性能の寿命を延ばすことを優先し、上限電位を下げる制御、又は、上限電位保持時間を短くする制御を行う。具体的な上限電位としては、1.0V以下、好ましくは0.9V以下、特に好ましくは0.85V以下とすることができる。
図3に示された運転時間後期においては、白金触媒粒子の表面積の低下度合いは緩やかになる。この段階では、燃費を向上させるため、上限電位を上げる。具体的な上限電位としては、0.5V以下、好ましくは0.7V以下、特に好ましくは0.85V以上とすることができる。
まず、粒径分布モデル作成方法を用いて、燃料電池触媒の粒径分布モデルを作成し、さらに、粒径分布モデル及び上記数式モデルを用いて、初期触媒粒径、触媒層構造、これまでの運転パターン及び運転時間等の物理的・化学的データから、現在の白金触媒粒子の表面積を推定する。
本第2の応用例においては、白金触媒粒子の表面積が低下し、当該表面積が一定の基準に達したら、ユーザーにメンテナンス時期であることを通知することができる。また、本第2の応用例においては、白金触媒粒子の表面積が一定の基準に達していなくても、推定した白金触媒粒子の表面積から、あとどのくらいの使用時間でメンテナンス時期を迎えるかをユーザーに通知することができる。
(1)排水性低下により、過剰水が反応場を減少させる(いわゆるフラッディング)。
(2)白金触媒表面積低下による触媒活性低下。
(3)白金消失層による、プロトン移動抵抗の増大。
なお、原因(3)における「白金消失層」とは、燃料電池の運転時間の増加に伴って、触媒層の電解質膜近傍に形成される層のことをいう。白金触媒は、燃料電池の負荷変動(例えば、当該燃料電池が移動体に用いられる場合には、当該移動体の加減速等)による電位変動の影響を受け、触媒の溶解及び析出を繰り返す。その結果、電解質膜近傍の触媒量減少(消失)という劣化現象が生じる。白金消失層は、当該劣化現象に伴って形成される層のことである。
白金消失層は、触媒層の電解質膜近傍に形成されるため、触媒層のプロトン移動距離が延長され、その結果プロトン移動抵抗が増大し、性能低下を招く原因となる。
性能低下原因が(2)の場合は、加湿に関する制御を行う必要はない。
性能低下原因が(3)の場合は、セル電圧低下を防止するため、燃料電池の加湿を上げるような制御(例えば、加湿モジュールで加湿する制御、ガス流入圧を下げる制御、反応ガス供給量を下げる制御等)が有効である。
このように、性能低下原因が(1)、(2)、(3)のそれぞれの場合について、異なる加湿の必要性の判定及び加湿の制御が必要である。
図4から分かるように、運転時間初期(黒菱形のプロット)においては触媒層において均等に高い値を保っていた白金触媒表面積も、使用時間を経るにつれて(白四角→黒三角→白丸→黒四角のプロット)、特に電解質膜界面側を中心に、急激に減少し始める。そして、白金触媒表面積減少と共に、白金消失層の厚さも、電解質膜界面側を中心に急激に厚くなる。
触媒層厚さ方向の白金触媒表面積分布の推定は、一定の運転時間ごと(例えば、100時間ごと)に計算してもよい。すなわち、運転の都度計算する必要は必ずしもない。
白金消失層の厚さに応じて、抵抗過電圧が生じる。この抵抗過電圧を、例えば下記式(8)のような計算式を用いて事前に計算しておき、制御フィードバックに用いることもできる。
ΔW = I2 × RPt消 × L 式(8)
(ΔW:白金消失層により生じる出力低下(W)、I:電流(A)、RPt消:白金消失層の単位面積当たり抵抗(Ω/cm2)、L:白金消失層の厚さ(cm))
これらの工程のうち、燃料電池触媒劣化予測工程については、本発明に係る燃料電池触媒の劣化予測方法の説明においてすでに述べた。また、セル電圧測定工程、セル抵抗測定工程における、燃料電池のセル電圧を測定する方法、燃料電池のセル抵抗を測定する方法については、それぞれ従来から知られている方法を用いることができる。
以下、加湿制御実行判定工程及び加湿制御実行工程について、詳細に説明する。
例えば、測定したセル電圧Vと、予め設定した閾値V1とを比較することによって、仮にV<V1であれば、何らかの異常が燃料電池内に発生した結果、燃料電池の性能低下が生じていることが分かる。
また、測定したセル抵抗Rと、予め設定したR1(許容可能なセル抵抗の最小値)とを比較することによって、仮にR<R1であれば、フラッディングが発生(上記性能低下原因(1))していると判定することができ、加湿を下げる制御を実行する必要があると判定できる。
さらに、測定したセル抵抗Rと、予め設定したR2(許容可能なセル抵抗の最大値)とを比較することによって、仮にR>R2であれば、乾燥による抵抗値上昇が生じていると判定することができ、加湿を上げる制御を実行する必要があると判定できる。
また、燃料電池触媒劣化予測工程から得られた劣化予測結果を考察することにより、局所的な燃料電池触媒表面積の低下の発生を予測することができれば、白金消失層による電圧低下(上記性能低下原因(3))と判定することができ、加湿を上げる制御を実行する必要があると判定できる。このような場合とは対照的に、局所的な燃料電池触媒表面積の低下が発生していないことを予測することができれば、白金触媒表面積低下による電圧低下(上記性能低下原因(2))と判定することができ、加湿制御は特に必要ではないと判定できる。
以上のように、燃料電池触媒劣化予測工程から得られた劣化予測結果、セル電圧測定工程から得られたセル電圧の値、及び、セル抵抗測定工程から得られたセル抵抗の値を個別に判定することにより、加湿制御の要否を判定することもできるが、これら予測結果及び測定値を総合的に判定することにより、加湿制御の要否を判定することもできる。
具体的な加湿の制御方法としては、加湿モジュールを用いる制御、ガス流入圧制御、反応ガス供給量制御及び冷却水流量制御を挙げることができる。特に、冷却水流量制御は、セル温度を制御することによって、間接的に加湿を制御することが可能である。
第1の加湿制御モードにおいては、ガス流入圧を上げる制御、反応ガス供給量を上げる制御及び冷却水流量を下げる制御(すなわち、セル温度を上げる制御)の少なくともいずれか1つの制御を実行することによって加湿を下げる制御を行う。
第2及び第3の加湿制御モードにおいては、ガス流入圧を下げる制御、反応ガス供給量を下げる制御及び冷却水流量を上げる制御(すなわち、セル温度を下げる制御)の少なくともいずれか1つの制御を実行することによって加湿を上げる制御を行う。
まず、図5中に示したステップS11において、セル電圧Vを測定し、次にステップS12において、測定したセル電圧Vと、予め設定した閾値V1とを比較する。V≧V1であれば、制御を終了し、V<V1であれば、ステップS13へと進む。
ステップS13においては、セル抵抗Rを測定し、さらにステップS14において、測定したセル抵抗Rと、予め設定した閾値R1(許容可能なセル抵抗の最小値)とを比較する。R<R1であれば、フラッディング(上記性能低下原因(1))と判定し、第1の加湿制御モードを選択し、加湿を下げる制御を実行する(ステップS15)。加湿を下げる制御を行った後は、ステップS13へと戻り、以下、抵抗値RがR≧R1となるまで第1の加湿制御モードを選択し、加湿を下げる制御を実行し続ける。ステップS14において、R≧R1であれば、ステップS16へと進む。
ステップS16において、測定したセル抵抗Rと、予め設定した閾値R2(許容可能なセル抵抗の最大値)とを比較する。R>R2であれば、第2の加湿制御モードを選択し、加湿を上げる制御を実行する(ステップS17)。加湿を上げる制御を行った後は、ステップS13へと戻り、以下、抵抗値RがR1≦R≦R2となるまで第2の加湿制御モードを選択し、加湿を上げる制御を実行し続ける。ステップS16において、R≦R2であれば、ステップS18へと進む。
ステップS18において、燃料電池触媒劣化予測工程から得られた劣化予測結果を用いて、電解質膜近傍の白金触媒表面積が低下しているか否かを予測する。当該予測においては、図4に示したような計算結果例のグラフを用いることができる。このようなグラフにより、局所的な燃料電池触媒表面積の低下、具体的には、電解質膜近傍の白金触媒表面積の低下が生じていると予測された場合は、白金消失層による電圧低下(上記性能低下原因(3))と判定し、加湿を上げる制御を実行する第3の加湿制御モードを選択し、加湿を上げる制御を実行する(ステップS21)。なお、局所的な燃料電池触媒表面積の低下、具体的には、電解質膜近傍の白金触媒表面積の低下が生じていないと予測された場合には、白金触媒表面積低下による電圧低下(上記性能低下原因(2))と判定し、特に加湿に関する制御は実行しない(ステップS20)。ステップS20又はS21において処理を終えた後、システム制御を終了させる。
本発明の粒径分布モデル作成方法を用いて、粒径分布モデルを作成した。当該モデルの作成にあたり、表計算ソフトとしてはMicrosoft Office Excel(商品名。Microsoft製)を、関数としては指定した平均と標準偏差に対する正規分布の累積分布関数の逆関数の値を返すNORMINV関数を用いた。なお、NORMINV関数は、NORMINV(確率x,平均μ,標準偏差σ)の3つの引数を有する。
まず、最小粒径として1.9(単位[nm]。以下同じ)、最大粒径として3.1の粒径範囲を決定した。なお、平均粒径を2.5とする。
次に、決定した粒径範囲を10の領域に区分した。具体的には、x=0.1、0.2、0.3、…0.9、μ=2.5、標準偏差σ=0.2をそれぞれNORMINV関数に代入し、各確率に対するF(X)−1の値、すなわち、累積分布関数の逆関数の値を算出し、区分点9個を決定した。算出結果を、下記表2に示す。
さらに、実際の触媒層内の厚さを考慮した、粒径分布モデル1を作成した。いま、上記11の代表点の粒径を有する粒径分布モデルを、触媒層厚さ方向に10層重ねたモデルを作成した(表3)。すなわち、このモデルは、触媒層の厚さ方向に均等な粒径分布を有することになる。
まず、最小粒径として1.4、最大粒径として2.6の粒径範囲を決定する。なお、平均粒径を2とする。
次に、決定した粒径範囲を10の領域に区分する。具体的には、x=0.1、0.2、0.3、…0.9、μ=2(単位[nm]。触媒粒径として好ましい値)、標準偏差σ=0.2をそれぞれNORMINV関数に代入し、各確率に対するF(X)−1の値、すなわち、累積分布関数の逆関数の値を算出し、区分点9個を決定した。算出結果を、下記表4に示す。
[実施例]
本発明の燃料電池触媒の劣化予測方法を用いて、運転時間(発電サイクル数)に対する白金表面積維持率を算出した。なお、粒径分布モデルとして、表3に示した粒径分布モデルを、初期状態の白金触媒粒径分布モデルとして使用した。図6(a)は、当該分布モデルの模式図である。
白金表面積維持率の算出については、数値計算ソフト(gPROMS)を用いた。方程式・境界条件・初期触媒粒径・触媒層の厚さ等の触媒層構造・運転パターンなどの各種パラメータを入力することで、時間変化の数値解が自動的に得られた。
このように、図6(b)に示すような、運転時間(発電サイクル数)に対する白金表面積維持率のグラフが得られた。なお、計算時間は20分程度であった。
実施例の粒径分布モデルとは別に、粒径の差を等間隔に取り、粒子数を変数とした粒径分布モデルを作成した。下記表7に、本比較例の粒径分布モデルの粒径(nm)、及び当該粒径における粒子数の、全ての粒子の合計の数に対する割合(ratio)を示した。図8(a)は、当該分布モデルの模式図である。このような粒径分布モデルを用いて、燃料電池触媒の劣化を予測し、運転時間(発電サイクル数)に対する白金表面積維持率を算出した。
このように、図8(b)に示すような、運転時間(発電サイクル数)に対する白金表面積維持率のグラフが得られた。なお、計算時間は20分程度であった。
これに対し、図8(b)から分かるように、比較例の計算結果(実線のグラフ)は、実施例の計算結果と比較すると、実験結果(グラフ上のプロット)のシミュレーションにおいて精確さに欠けることが分かる。
2…アノード触媒層
3…カソード触媒層
4,5…ガス拡散層
6…アノード電極
7…カソード電極
8…膜・電極接合体
9,10…セパレータ
11,12…ガス流路
20…実際の粒径分布を示した曲線
100…燃料電池の単セル
Claims (9)
- 粒径が異なる複数の粒子が集合してなる粒子群の粒径分布を模擬する粒径分布モデルを作成する方法であって、
模擬すべき粒子群に含まれる粒子の最小粒径及び最大粒径を決定する粒径範囲決定工程、
前記粒径範囲決定工程で決定された粒径範囲における粒子の出現頻度を、前記最小粒径及び前記最大粒径を起点及び終点として積分する積分工程、
前記積分工程で得られた出現頻度の積分値を用いて、前記積分工程で用いた積分領域を、第1の区分を行うことによって、区分された個々の領域の積分値が均等になるように所定数の領域に区分し、区分点となった粒径を決定する区分点決定工程、並びに、
前記最小粒径、前記最大粒径及び前記各区分点の粒径を代表点として決定する代表点決定工程、を含み、
各代表点の粒径を有する粒子を含有し、且つ、各代表点の粒径を有する粒子の出現頻度が当該代表点で区分された個々の領域の積分値と等しい粒径分布を仮想し、当該仮想された粒径分布を、粒径分布モデルとして得ることを特徴とする、粒径分布モデル作成方法。 - 前記区分点決定工程において、前記第1の区分によって区分された所定数の領域の内の一部又は全部について、第2の区分を行うことによって、区分された個々の領域の積分値が均等になるように所定数の領域にさらに区分し、前記第1の区分及び前記第2の区分において区分点となった粒径を決定する、請求項1に記載の粒径分布モデル作成方法。
- 前記第2の区分を、前記第1の区分によって区分された所定数の領域の内、前記最大粒径を含む領域について行う、請求項2に記載の粒径分布モデル作成方法。
- 前記粒径分布モデルにおいて、前記出現頻度が、粒子数、質量又は体積で表される、請求項1乃至3のいずれか一項に記載の粒径分布モデル作成方法。
- 前記粒径分布が正規分布である、請求項1乃至4のいずれか一項に記載の粒径分布モデル作成方法。
- 前記請求項1乃至5のいずれか一項に記載の粒径分布モデル作成方法を用いて、燃料電池触媒の粒径分布モデルを作成する工程、及び、
前記燃料電池触媒の粒径分布モデルを用いて燃料電池触媒の劣化を予測する工程を有することを特徴とする、燃料電池触媒の劣化予測方法。 - 前記燃料電池触媒劣化予測工程において、溶解反応速度の数式モデル、酸化反応速度の数式モデル及び物質収支の数式モデルから選ばれる少なくとも1つ以上の数式モデルを用いる、請求項6に記載の燃料電池触媒の劣化予測方法。
- 前記燃料電池触媒劣化予測工程において、白金溶解反応速度の数式モデル、白金酸化反応速度の数式モデル、酸化白金(II)の溶解反応速度の数式モデル及び物質収支の数式モデルを全て用いる、請求項6又は7に記載の燃料電池触媒の劣化予測方法。
- 前記請求項6乃至8のいずれか一項に記載の燃料電池触媒の劣化予測方法を用いて燃料電池触媒の劣化を予測する工程、
燃料電池のセル電圧を測定する工程、
燃料電池のセル抵抗を測定する工程、
前記燃料電池触媒劣化予測工程から得られた劣化予測結果、前記セル電圧測定工程から得られたセル電圧の値、及び、前記セル抵抗測定工程から得られたセル抵抗の値に基づいて、加湿制御を実行するか否かを判定する工程、並びに、
前記加湿制御判定工程の判定結果において、
前記セル抵抗測定工程から得られたセル抵抗の値が所定の範囲の値未満であれば、加湿を下げる制御を実行する第1の加湿制御モードを、
前記セル抵抗測定工程から得られたセル抵抗の値が所定の範囲の値を超える値であれば、加湿を上げる制御を実行する第2の加湿制御モードを、
前記燃料電池触媒劣化予測工程から得られた劣化予測結果より、局所的な燃料電池触媒表面積の低下が生じていると予測された場合は、加湿を上げる制御を実行する第3の加湿制御モードを、それぞれ選択し実行する工程、を有することを特徴とする、燃料電池の制御方法。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009084428A JP5212224B2 (ja) | 2009-03-31 | 2009-03-31 | 粒径分布モデル作成方法、当該粒径分布モデル作成方法を用いた燃料電池触媒の劣化予測方法、及び当該燃料電池触媒の劣化予測方法を用いた燃料電池の制御方法 |
US12/751,318 US8660827B2 (en) | 2009-03-31 | 2010-03-31 | Method of creating particle size distribution model, method of predicting degradation of fuel cell catalyst using the method of creating particle size distribution model, and method of controlling fuel cell using the method of predicting degradation of fuel cell catalyst |
US14/149,540 US9002689B2 (en) | 2009-03-31 | 2014-01-07 | Method of predicting degradation of fuel cell catalyst using the method of creating particle size distribution model |
US14/149,592 US8843355B2 (en) | 2009-03-31 | 2014-01-07 | Method of controlling fuel cell using the method of predicting degradation of fuel cell catalyst |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009084428A JP5212224B2 (ja) | 2009-03-31 | 2009-03-31 | 粒径分布モデル作成方法、当該粒径分布モデル作成方法を用いた燃料電池触媒の劣化予測方法、及び当該燃料電池触媒の劣化予測方法を用いた燃料電池の制御方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2010236989A true JP2010236989A (ja) | 2010-10-21 |
JP5212224B2 JP5212224B2 (ja) | 2013-06-19 |
Family
ID=42784659
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009084428A Active JP5212224B2 (ja) | 2009-03-31 | 2009-03-31 | 粒径分布モデル作成方法、当該粒径分布モデル作成方法を用いた燃料電池触媒の劣化予測方法、及び当該燃料電池触媒の劣化予測方法を用いた燃料電池の制御方法 |
Country Status (2)
Country | Link |
---|---|
US (3) | US8660827B2 (ja) |
JP (1) | JP5212224B2 (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013171727A (ja) * | 2012-02-21 | 2013-09-02 | Toyota Motor Corp | 燃料電池システムとその運転方法 |
JP5915758B2 (ja) * | 2012-09-12 | 2016-05-11 | 株式会社島津製作所 | 粒度分布データ処理装置及びこれを備えた粒度分布測定装置、並びに、粒度分布データ処理方法及び粒度分布データ処理プログラム |
JP2017224575A (ja) * | 2016-06-17 | 2017-12-21 | 三菱自動車工業株式会社 | 燃料電池システム |
KR20180114971A (ko) * | 2017-04-11 | 2018-10-22 | 현대자동차주식회사 | 연료전지 정지 모드 제어 방법 및 시스템 |
CN108918356A (zh) * | 2018-05-18 | 2018-11-30 | 中建材(合肥)粉体科技装备有限公司 | 一种混合物料挤压粉碎后的粒度预测方法及系统 |
WO2023058394A1 (ja) * | 2021-10-06 | 2023-04-13 | 栗田工業株式会社 | 微粒子測定方法 |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11909077B2 (en) | 2008-06-12 | 2024-02-20 | Massachusetts Institute Of Technology | High energy density redox flow device |
US9786944B2 (en) | 2008-06-12 | 2017-10-10 | Massachusetts Institute Of Technology | High energy density redox flow device |
US8722226B2 (en) * | 2008-06-12 | 2014-05-13 | 24M Technologies, Inc. | High energy density redox flow device |
US9484569B2 (en) | 2012-06-13 | 2016-11-01 | 24M Technologies, Inc. | Electrochemical slurry compositions and methods for preparing the same |
US9362583B2 (en) | 2012-12-13 | 2016-06-07 | 24M Technologies, Inc. | Semi-solid electrodes having high rate capability |
US8993159B2 (en) | 2012-12-13 | 2015-03-31 | 24M Technologies, Inc. | Semi-solid electrodes having high rate capability |
US9310424B2 (en) | 2013-02-25 | 2016-04-12 | International Business Machines Corporation | Monitoring aging of silicon in an integrated circuit device |
US8713490B1 (en) * | 2013-02-25 | 2014-04-29 | International Business Machines Corporation | Managing aging of silicon in an integrated circuit device |
CN106299401B (zh) | 2015-05-20 | 2019-07-16 | 通用电气公司 | 燃料电池系统及其控制方法 |
WO2017124098A1 (en) | 2016-01-15 | 2017-07-20 | 24M Technologies, Inc. | Systems and methods for infusion mixing a slurry-based electrode |
US10615438B2 (en) | 2018-02-23 | 2020-04-07 | Cummins Enterprise Llc | Degradation detecting device for fuel cell stack, fuel cell system and managing method thereof |
CN111630365B (zh) * | 2018-02-27 | 2023-02-24 | 松下知识产权经营株式会社 | 粒子检测传感器 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08162133A (ja) * | 1994-12-05 | 1996-06-21 | Agency Of Ind Science & Technol | 白金系触媒の製造方法 |
JP2003197205A (ja) * | 2001-12-26 | 2003-07-11 | Asahi Glass Co Ltd | 固体高分子型燃料電池の製造方法及びそのためのガス拡散電極の製造方法 |
JP2007190454A (ja) * | 2006-01-17 | 2007-08-02 | Japan Advanced Institute Of Science & Technology Hokuriku | メタノール酸化用PtRu系触媒及びその製造方法 |
JP2009063565A (ja) * | 2007-08-09 | 2009-03-26 | Sharp Corp | 画像処理方法、およびそれを用いた粒度分布測定方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01122570A (ja) | 1987-11-07 | 1989-05-15 | Fuji Electric Co Ltd | 燃料電池の寿命予測方法 |
US6136452A (en) * | 1998-02-27 | 2000-10-24 | The Regents Of The University Of California | Centrifugal synthesis and processing of functionally graded materials |
US6119510A (en) * | 1998-08-31 | 2000-09-19 | Lucent Technologies Inc. | Process for determining characteristics of suspended particles |
US6481268B1 (en) * | 1999-10-12 | 2002-11-19 | Baker Hughes, Inc. | Particle measurement by acoustic speckle |
GB2364052B (en) * | 2000-06-27 | 2003-05-14 | Univ London | A method and apparatus for producing a biomaterial product |
JP2006139935A (ja) | 2004-11-10 | 2006-06-01 | Nissan Motor Co Ltd | 燃料電池システムの制御装置 |
JP2007103115A (ja) | 2005-10-03 | 2007-04-19 | Nissan Motor Co Ltd | 燃料電池システム |
US20100169038A1 (en) * | 2007-04-23 | 2010-07-01 | Tufts University | Methods for calculating particle size distribution |
US8057775B2 (en) * | 2008-04-22 | 2011-11-15 | Us Synthetic Corporation | Polycrystalline diamond materials, methods of fabricating same, and applications using same |
-
2009
- 2009-03-31 JP JP2009084428A patent/JP5212224B2/ja active Active
-
2010
- 2010-03-31 US US12/751,318 patent/US8660827B2/en active Active
-
2014
- 2014-01-07 US US14/149,592 patent/US8843355B2/en active Active
- 2014-01-07 US US14/149,540 patent/US9002689B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08162133A (ja) * | 1994-12-05 | 1996-06-21 | Agency Of Ind Science & Technol | 白金系触媒の製造方法 |
JP2003197205A (ja) * | 2001-12-26 | 2003-07-11 | Asahi Glass Co Ltd | 固体高分子型燃料電池の製造方法及びそのためのガス拡散電極の製造方法 |
JP2007190454A (ja) * | 2006-01-17 | 2007-08-02 | Japan Advanced Institute Of Science & Technology Hokuriku | メタノール酸化用PtRu系触媒及びその製造方法 |
JP2009063565A (ja) * | 2007-08-09 | 2009-03-26 | Sharp Corp | 画像処理方法、およびそれを用いた粒度分布測定方法 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013171727A (ja) * | 2012-02-21 | 2013-09-02 | Toyota Motor Corp | 燃料電池システムとその運転方法 |
JP5915758B2 (ja) * | 2012-09-12 | 2016-05-11 | 株式会社島津製作所 | 粒度分布データ処理装置及びこれを備えた粒度分布測定装置、並びに、粒度分布データ処理方法及び粒度分布データ処理プログラム |
JP2017224575A (ja) * | 2016-06-17 | 2017-12-21 | 三菱自動車工業株式会社 | 燃料電池システム |
KR20180114971A (ko) * | 2017-04-11 | 2018-10-22 | 현대자동차주식회사 | 연료전지 정지 모드 제어 방법 및 시스템 |
KR102336395B1 (ko) * | 2017-04-11 | 2021-12-08 | 현대자동차주식회사 | 연료전지 정지 모드 제어 방법 및 시스템 |
CN108918356A (zh) * | 2018-05-18 | 2018-11-30 | 中建材(合肥)粉体科技装备有限公司 | 一种混合物料挤压粉碎后的粒度预测方法及系统 |
CN108918356B (zh) * | 2018-05-18 | 2021-07-09 | 中建材(合肥)粉体科技装备有限公司 | 一种混合物料挤压粉碎后的粒度预测方法及系统 |
WO2023058394A1 (ja) * | 2021-10-06 | 2023-04-13 | 栗田工業株式会社 | 微粒子測定方法 |
JP2023055462A (ja) * | 2021-10-06 | 2023-04-18 | 栗田工業株式会社 | 微粒子測定方法 |
JP7268710B2 (ja) | 2021-10-06 | 2023-05-08 | 栗田工業株式会社 | 微粒子測定方法 |
Also Published As
Publication number | Publication date |
---|---|
US8843355B2 (en) | 2014-09-23 |
US8660827B2 (en) | 2014-02-25 |
US20140120435A1 (en) | 2014-05-01 |
JP5212224B2 (ja) | 2013-06-19 |
US9002689B2 (en) | 2015-04-07 |
US20100248046A1 (en) | 2010-09-30 |
US20140114628A1 (en) | 2014-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5212224B2 (ja) | 粒径分布モデル作成方法、当該粒径分布モデル作成方法を用いた燃料電池触媒の劣化予測方法、及び当該燃料電池触媒の劣化予測方法を用いた燃料電池の制御方法 | |
JP5044428B2 (ja) | 燃料電池スタックの分極曲線を計算する方法及びそれを用いた燃料電池システム | |
Loo et al. | Characterization of the dynamic response of proton exchange membrane fuel cells–a numerical study | |
JP4273781B2 (ja) | 燃料電池の運転状態判定装置及びその方法 | |
Hu et al. | Modelling and simulations of carbon corrosion during operation of a polymer electrolyte membrane fuel cell | |
CN102412405B (zh) | 预测燃料电池堆中的燃料电池的最小电池电压趋势的方法 | |
US20130189596A1 (en) | Fuel cell system, method and program of determining cause of negative voltage, and storage medium storing program | |
US20140147759A1 (en) | Vehicle application for air storage cathode carbon loss estimation | |
Aindow et al. | Platinum catalyst degradation in phosphoric acid fuel cells for stationary applications | |
JP2009032418A (ja) | 燃料電池の運転方法 | |
Bernhard et al. | Model-assisted analysis and prediction of activity degradation in PEM-fuel cell cathodes | |
US11196064B2 (en) | Fuel cell controller, fuel cell system and method of operation | |
JP2007179749A (ja) | 燃料電池の制御方法及びその制御装置 | |
CN107017425B (zh) | 燃料电池系统和燃料电池系统的性能改进方法 | |
KR101349022B1 (ko) | 연료전지 스택의 촉매 열화 진단 방법 | |
JP2010061887A (ja) | 燃料電池用の電極触媒の性能評価方法および評価装置 | |
US7608351B2 (en) | System and method for controlling cathode stoichiometry to minimize RH excursions during transients | |
JP5428496B2 (ja) | 燃料電池システムの触媒被毒状態の判定 | |
JP2005108673A (ja) | 燃料電池システム | |
JP7248004B2 (ja) | 燃料電池制御指令装置 | |
JP2007042477A (ja) | 燃料電池システム | |
US20230411654A1 (en) | Fuel cell purge system based on tilt location | |
JP2010267516A (ja) | 燃料電池システム | |
TOUIL et al. | A System-Level Modeling Framework for Predicting Pt Catalyst Degradation in Proton Exchange Membrane Fuel Cells | |
JP2007157554A (ja) | 固体高分子型燃料電池の運転方法およびその運転装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20111114 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20120703 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130121 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20130129 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20130211 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 5212224 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20160308 Year of fee payment: 3 |