Single chamber fuel cell models
1 / 32

Single-Chamber Fuel Cell Models - PowerPoint PPT Presentation

  • Uploaded on

Single-Chamber Fuel Cell Models. D. G. Goodwin, Caltech. Develop validated physics-based models of SCFC operation Use models along with test results to develop understanding of factors determining performance Use to aid in design optimization. Multiple models. Computational expense.

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
Download Presentation

PowerPoint Slideshow about 'Single-Chamber Fuel Cell Models' - bridie

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
Single chamber fuel cell models
Single-Chamber Fuel Cell Models

D. G. Goodwin, Caltech

  • Develop validated physics-based models of SCFC operation

  • Use models along with test results to develop understanding of factors determining performance

  • Use to aid in design optimization

Multiple models
Multiple models

Computational expense

  • Model 1: a simple model for qualitative parametric studies

    • Allows rapid exploration of the effects of various parameters on performance

  • Model 2: Solves 2D channel flow assuming fully developed flow. Computes

    • Species concentration profiles

    • Current density profiles

    • Power output vs. load

  • Model 3: Solves 2D reacting channel flow accurately (in development, Yong Hao)

Seconds on a laptop PC

Minutes on a linux workstation

Minutes to hours

Model 1 a zero dimensional fuel cell model

Can be used to model single- or dual-chamber designs

No consideration of gas flow

Approximate equilibrium treatment of hydrocarbon oxidation

Includes diffusion through electrodes, activation polarizations, ohmic losses

Can compute current-voltage curves

Written in a simple scripting language (Python)

Uses the Cantera software package to evaluate thermodynamic and transport properties, and compute chemical equilibrium (

Good for semi-quantitative parametric studies

Model 1: a “zero-dimensional” fuel cell model

Idealized geometry

Uniform cathode-side gas

Porous Cathode

Porous Anode

Uniform anode-side gas

Idealized Geometry

  • Each side exposed to uniform gas with specified composition

    • No depletion in gas

    • Corresponds to limit of fast transport

    • Compositions can be set equal (single-chamber) or each independently specified (dual-chamber)

Electrochemical reactions
Electrochemical Reactions

  • Anode reactions

    • H2 + O2- = H2O

    • CO + ½ O2- = CO2

  • Cathode reaction

    • O2 = 2O2-

  • Catalyst selectivity

    • Reactions allowed to occur at opposite electrode with relative rate 0 < Fc < 1

    • Fc > 0 lowers OCV

    • At Fc = 0, OCV = 0

Gas composition
Gas Composition





  • Approximate treatment of partial oxidation

  • Assume gas is a mixture of the input gas composition + equilibrium composition

  • No selectivity assumed – CO, CO2, H2, and H2O all present


1 - Feq

Input: 1:3:12 C3H8/O2/He


gas mixture

600 C

Transport through electrodes

concentration gradients in electrode drive diffusion


Assumed uniform gas composition



Transport through electrodes

  • Gas composition at electrode/electrolyte interface determined by diffusion through porous electrode

  • Effective diffusion coefficients account for pore size, porosity, and tortuosity of electrode microstructure

  • Concentrations at electrode/electrolyte interface used to calculate Nernst potential

Cathode activation polarization
Cathode Activation Polarization

  • Represents largest loss

  • Dependence on oxygen partial pressure assumed first-order

Range considered

Anode activation polarization
Anode Activation Polarization

  • Assumed not to be rate-limiting

  • Anode exchange current density set to a large multiple of cathode exchange current density (100 – 1000)

Electrolyte ohmic loss
Electrolyte Ohmic Loss

Value for GDC used

Current density computation
Current Density Computation

  • Nernst potential calculated using concentrations at electrode/electrolyte interfaces, and includes effects of back reaction

  • Given Eload, this equation is solved for the current density

Simulation of test results with ni sdc sdc ssc pt sdc at 600 c

Simulation of Test Results with Ni-SDCSDCSSC-Pt-SDC at 600 C

Test results can be accounted for with physically reasonable parameters
Test results can be accounted for with physically-reasonable parameters

  • Experimental Ni-SDCSDCSSC-Pt-SDC results at 600 C best fit by

    • I0,c = 70 mA/cm2

    • 80% electrode selectivity

    • 50% conversion to equilibrium products

Accurate modeling of transport limit requires more accurate treatment of transport processes – see Model 2 results

Gas composition effects
Gas Composition Effects parameters

  • Increasing percent conversion to equilibrium products moves the transport limit to higher current densities

  • For fuel-rich input mixtures, equilibrium composition contains significant CO and H2, in addition to CO2 and H2O

  • Therefore, non-electrochemical oxidation of CO and H2 not likely to be a problem as long as a fuel-rich mixture is used



Single chamber vs dual chamber
Single Chamber vs. Dual Chamber parameters

  • Dual chamber calculation sets cathode gas composition to air, and eliminates the back reactions at the electrodes



Catalyst selectivity effects
Catalyst selectivity effects parameters

Catalysts must have reasonable selectivity for electrochemical reactions in order for SCFC to function

More selective

Less selective

Scfc loss mechanisms
SCFC Loss Mechanisms parameters

  • Dominated by losses due to

    • Low cathode activity

    • Incomplete cathode and anode selectivity

Model overview
Model Overview parameters

  • Inputs

    • Inlet gas composition, temperature, pressure

    • Load potential

    • Parameters characterizing kinetics, electrode transport, geometry, etc.

  • Outputs

    • 2D spatial distributions of C3H8, CH4, CO, H2, CO2, and H2O in channel

    • Current density profile J(x)

  • Assumes isothermal, isobaric conditions

  • Includes an unsealed, non-catalytic plate (interconnect) separating anode and cathode gas streams

Model geometry

Anode parameters

Model Geometry



Cathode-side flow channel


Fuel / air mixture

Non-catalytic partition

Anode-side flow channel

Mathematical model

Species equations finite-differenced and integrated in time to steady state.

Porous electrodes handled by locally modifying diffusion coefficients

Species equations solved simultaneously with equation for current density

Mathematical Model

Model problem
Model Problem to steady state.

  • Channel height = 700 mm, length = 10 mm

  • 200 mm anode, 50 mm cathode

  • Electrode porosity 0.4, pore size 0.1 mm

  • 15 mm GDC electrolyte

  • T = 600 C, P = 1 atm

  • Premixed 1:3 C3H8 / air

  • Partial oxidation rate at anode set to give nearly complete consumption of propane

  • Other parameters same as in zero-D model

Porous electrode transport

reactants to steady state.



Porous Electrode Transport

  • Gas must diffuse through porous electrodes to reach electrochemically-active triple-phase boundary

  • Process modeled with effective diffusion coefficients for each species that interpolate between Knudsen and ideal gas limits

  • Effective diffusion coefficient close to the Knudsen limit

Partial oxidation
Partial Oxidation to steady state.

  • Global partial oxidation reaction C3H8 + 3/2 O2 => CO + 4H2

    • Produces electrochemically-active species

    • assumed to occur throughout the anode

    • May occur on the cathode also

  • Rate modeled as first-order in C3H8 and O2

  • Magnitude set to lead to nearly complete conversion in the anode-side exhaust

    • ample residence time for complete conversion (50-100 ms vs. 1 ms)

    • Degree of conversion can be tuned experimentally by material choice, and anode fabrication methods

Velocity profile
Velocity Profile to steady state.

This velocity profile is imposed, based on known solution for viscous fully-developed flow

Porous cathode

Porous anode

Species distributions at max power
Species Distributions at Max Power to steady state.

Cathode on right

Anode on left


Current density distribution
Current Density Distribution to steady state.

  • Movie shows steady-state J(x) for load potentials ranging from zero to 0.9 V

Predicted performance at 600 c
Predicted Performance at 600 C to steady state.

  • Predicted OCV = 0.9 V, peak power density = 85 mW/cm2

  • Easily meets target SCFC performance of 50 – 100 mW/cm2.

Conclusions to steady state.

  • Performance targets appear to be easily achievable

  • Largest potential gains in performance:

    • improved cathode catalytic activity

    • improved electrode selectivity

  • Separator plate may not be necessary

  • As long as gas composition is fuel rich, non-electrochemical oxidation of CO and H2 will not go to completion, and therefore nonselective catalyst for partial oxidation is acceptable.

Future work
Future Work to steady state.

  • Validation against all available test data – single-chamber, dual-chamber, etc.

  • Prediction of coking behavior

  • Prediction of low-temperature performance

  • Integration with Swiss Roll heat exchanger model to predict operating temperature

Summary to steady state.

  • Two numerical models have been developed to predict SCFC performance.

    • A simple model useful for interpreting test data

    • A channel flow model useful for predicting micropower generator performance

  • Test results can be accounted for with physically-reasonable kinetic parameters

  • Using these parameters in the channel-flow model leads to performance at 600 C that meets our targets

  • Both models are suitable for use in design and optimization studies, including system studies with the Swiss Roll heat exchanger.