seca core program – recent development of … library/research/coal/energy systems...seca core...
TRANSCRIPT
1/30/2003
1
SECA Core Program – Recent Development of Modeling Activities at PNNL
MA KhaleelEmail: [email protected] Phone (509) 375-2438
KP Recknagle, DR Rector, J Deibler, RT Pagh,RE Williford, LA Chick
Pacific Northwest National LaboratoryRichland, WA 99352
June 19, 2002
1/30/2003
2
Technical Issues and R& D ObjectivesTechnical Issues• Concurrent management and control of thermal, physical, chemical
and electrochemical processes over various SOFC operational parameters.
ObjectivesDevelop modeling and simulation tools to be used by the SECA
vertical teams as an integral part of the design process. Tools to be used for:
– System design requirements roll-out– Sub-system component design– Microstructural and material design/optimization– Control design – Life prediction
1/30/2003
3
Technical ApproachFlow Analysis
Thermal Analysis
Stress Analysis
Flow, thermal & Electrochemistry analysis
Thermal Shock
Electrical Power System
Thermal system
Stack
Cell
Modeling levels
Cont
inuu
m-le
v el e
le ctro
chem
is try
Too l
s an d
Met
hodo
logi
esfo
r Rap
id st
art-u
p
Validation and Property measurement
Methane distributionCell-
level
Microstructure-level
Power density profile
H2O concentration
%H2 97.0% 0.970 HEAT GENERATION
200 sccm CO 0.0% -0.2088 W/cm2
1.49E-04 mol/s H2O 3.0% (Minus sign denotes exotherm)CO2 0.0% i = 0.5 A/cm2
N2 0.0% Vi = 0.869 voltsTotal 100.0% P= 0.00 kW 1.65 W
P= 0.43 W/cm2Fraction of CO that is water-gas shifted = 0.5 Fuel Utilization = 6.82%
mol/sec Molar % moles/sec Molar %H2 1.44E-04 97.00% Active cell area= 3.8 cm2 H2 1.34E-04 90.38%CO 9.64E-10 0.00% CO 8.85E-10 0.00%
H2O 4.46E-06 3.00% Thickness ,µm %Porosity Tortuosity H2O 1.43E-05 9.62%CO2 3.64E-11 0.00% Electrolyte= 10 na na CO2 1.15E-10 0.00%
N2 4.13E-21 0.00% Anode= 600 30 2.50 N2 4.13E-21 0.00%CH4 0.00E+00 0.00% Interconnect= 0 na na CH4 0.00E+00 0.00%Po2 ------ 2.67E-23 Cathode= 50 30 2.50 O2 ----- 3.16E-22
Total 1.49E-04 1.00E+00 Contact Resistance= 0 Ohm-cm2 Total 1.49E-04 1.00E+00Ave. Stack Temp= 750 °C
750 °C Inlet Temp 1023 K 750 °C Outlet Temp1023 K B-V Parameters: α = 0.518 1023 K
750 °C Inlet Temp Px = 127 750 °C Outlet Temp1023 K QH2= 0.425 Eact = 110 1023 K
Dsurf-H2 0.0107mol/sec P, atm mol/sec P, atm
O2 4.50E-05 0.21 O2 4.01E-05 0.19N2 1.69E-04 0.79 Offset voltage due to leak= -0.06 N2 1.69E-04 0.81
Total 0.0002 1.00 Total 2.09E-04 1.00
ELECTRIC WORK
Anode Exhaust
Cathode Inlet Gas Stream Cathode Outlet Gas Stream
Fuel Stream
Fuel Gas Input Parameters
e-
Excel SystemModel
1/30/2003
4
%H2 97.0% 0.970 HEAT GENERATION
200 sccm CO 0.0% -0.2088 W/cm2
1.49E-04 mol/s H2O 3.0% (Minus sign denotes exotherm)CO2 0.0% i = 0.5 A/cm2
N2 0.0% Vi = 0.869 voltsTotal 100.0% P= 0.00 kW 1.65 W
P= 0.43 W/cm2Fraction of CO that is water-gas shifted = 0.5 Fuel Utilization = 6.82%
mol/sec Molar % moles/sec Molar %H2 1.44E-04 97.00% Active cell area= 3.8 cm2 H2 1.34E-04 90.38%CO 9.64E-10 0.00% CO 8.85E-10 0.00%
H2O 4.46E-06 3.00% Thickness ,µm %Porosity Tortuosity H2O 1.43E-05 9.62%CO2 3.64E-11 0.00% Electrolyte= 10 na na CO2 1.15E-10 0.00%
N2 4.13E-21 0.00% Anode= 600 30 2.50 N2 4.13E-21 0.00%CH4 0.00E+00 0.00% Interconnect= 0 na na CH4 0.00E+00 0.00%Po2 ------ 2.67E-23 Cathode= 50 30 2.50 O2 ----- 3.16E-22
Total 1.49E-04 1.00E+00 Contact Resistance= 0 Ohm-cm2 Total 1.49E-04 1.00E+00Ave. Stack Temp= 750 °C
750 °C Inlet Temp 1023 K 750 °C Outlet Temp1023 K B-V Parameters: α = 0.518 1023 K
750 °C Inlet Temp Px = 127 750 °C Outlet Temp1023 K QH2= 0.425 Eact = 110 1023 K
Dsurf-H2 0.0107mol/sec P, atm mol/sec P, atm
O2 4.50E-05 0.21 O2 4.01E-05 0.19N2 1.69E-04 0.79 Offset voltage due to leak= -0.06 N2 1.69E-04 0.81
Total 0.0002 1.00 Total 2.09E-04 1.00
ELECTRIC WORK
Anode Exhaust
Cathode Inlet Gas Stream Cathode Outlet Gas Stream
Fuel Stream
Fuel Gas Input Parameters
e-
Input Parameters:•Fuel composition and flow rate•Temperature•Adjustable cell parameters
Output Data:•I-V curve and power•Heat generation•Fuel utilization and exhaust composition
Cell Electrochemistry Performance Model
1/30/2003
5
Single Ce ll Da ta
0 .4
0 .5
0 .6
0 .7
0 .8
0 .9
1 .0
1 .1
0 .0 0 .2 0 .4 0 .6 0 .8 1 .0 1 .2 1 .4 1 .6 1 .8 2 .0Curre nt De ns ity, A/c m 2
Cel
l Vol
tage
, V
97% H2, 750°C75% H2, 750°C50% H2, 750°C25% H2, 750°C10% H2, 750°CM odel, M ole F rac . H2= 0.1M odel, M ole F rac . H2= 0.25M odel, M ole F rac . H2= 0.49M odel, M ole F rac . H2= 0.74M odel, M ole F rac . H2= 0.97
Single Cell Data
0.40
0.50
0.60
0.70
0.80
0.90
1.00
1.10
0.00 0.50 1.00 1.50 2.00
Current Density, A/cm2
Cel
l Vol
tage
, V
800°C, 97%H2750°C, 97%H2700°C, 97%H2650°C, 97%H2Model, T=800Model, T=750Model, T=700Model, T=675Model, T=650
•Cell parameters can be adjusted so that one set of cell parameters provide excellent fit of a family of IV curves for a “unit” cell operated over a range of temperatures and a range of hydrogen concentrations. • The “calibrated” model can then be used to predict large stack performance by applying it within a CFD code to computational unit cells.
Temperaturevariations
H2 concentrationvariations
Model Output
1/30/2003
6
Topics for Continuum Electrochemistry Modeling
Calculation Flow Diagram For 1-Cell StackGeneric Cross-Flow CaseAlternate Flow ConfigurationsTransition from Transient Heating to Steady StateCalculation Flow Diagram For Multiple-Cell StacksMultiple Cell Stack Modeling Results
1/30/2003
7
STAR-CD/EC (1-Cell Stack) Flow Diagram
STAR-CD
POSDAT
Temperature, Gas Chemistry (mass conc.) Model Geometry
ECT, pp(i)
Current, Scalar (molar), and Enthalpy source terms
Scalar
Mass
Enthalpy
SORSCA
FLUINJ
SORENT
Source Terms
(Source Terms)
Cell Voltage
1/30/2003
8
MARC/EC (1-Cell Stack) Flow Diagram
MARCFlow – Thermal - Electrical -
Mechanical
STRESSES
FLUX
Heat Flux PLOTV
ELEVAR
EC
T, pp(i)Cell Voltage
Gas, Chemistry,
Heat Generation
UpdateState
Variables
Model Geometry, Temperature, Current Distribution
1/30/2003
9
Cross-Flow Case: ConditionsInflow Air & Fuel Temperature = 625°CAir delivery rate = 15 gm/sec/60cellsFuel delivery rate = 1.08 gm/sec/60cells (9.5x10-4 moles/sec)• Composition: shifted to equilibrium at 625 °C
– [H2] = 0.37443– [H2O] = 0.03449– [CO] = 0.33662– [CO2] = 0.06759– [N2] = 0.18687
Cell Voltage = 0.7 (as in all other cases)Cyclic boundary conditions at top and bottom of unit cell.
1/30/2003
10
Cross-Flow Case: ResultsCurrent Density & Heat Generation
62% Fuel Utilization Case:Current Density = 0.300-1.46 (0.687 A/cm2)Heat Generation = 0.21 - 0.99 (0.477 W/cm2)PEN Temperature = 643 - 912 (769 °C)
Air Temperature, Kelvin
Fuel
Air
H2, Mass Fraction
O2, Mass Fraction PEN Temperature, K
in, moles/s out, moles/sh2 3.5480E-04 1.3963E-04h2o 3.2680E-05 2.4776E-04co 3.1898E-04 1.1871E-04co2 6.4050E-05 2.6423E-04n2 1.7706E-04 1.7704E-04moles/s 9.4756E-04 9.4737E-04
Fuel Utilization= 61.7%
1/30/2003
11
Cross-Flow Case: Results (Continued)CO, Mass Fraction
CO2, Mass Fraction
H2O, Mass Fraction
Heat Generation, W/cm2
1/30/2003
12
Alternate Flow Configurations – Steady StateFlow configuration
Conditions Results
Air delivery, gm/s @ °C
Fuel delivery, gm/s @ °C
∆TPEN, °C TPEN, °C Iave, A/cm2 Fuel Utilization
Cross-flow 0.25 @ 625
0.018 @ 625
269 769 0.69 62%
Co-flow 0.25 @ 625
0.018 @ 625
184 763 0.71 64%
Counter-flow
0.25 @ 595
0.018 @ 595
267 758 0.73 63%
1/30/2003
13
Alternate Flow Configurations – Steady State
Temperature,Degrees C(∆T=180 co-∆T=270 cross-)
04/09/02
Current Density,A/cm2
H2 MassConc.,kg/kg
CROSS-FLOW CO-FLOW COUNTER-FLOW
1/30/2003
14
Modeling of Transition from Transient to Steady State – Cross Flow Case
2% of transition time 10% 20% 30%
40% 50% 80% 100% Steady State Solution
Start with Transient Temperatures 04/09/02(temperature range ~900-1075K)
(temperature range = 931-1185K)Is the heating method sound?Will start of reaction cause instabilities?Can we actually shorten the “startup” time?
1/30/2003
15
STAR-CD/EC (Multiple-Cell) Flow Diagram
Sum over cell to get total cell currentCell Current = Stack Current?
yes no
AdjustCell
Voltage
STAR-CD
POSDAT
Temperature, Gas Chemistry (mass conc.) Model Geometry
ECT, pp(i)
Current, Scalar (molar), and Enthalpy source terms
Scalar
Mass
Enthalpy
SORSCA
FLUINJ
SORENT
Source Terms
(Source Terms)
Cell Voltage
1-cell stack
Multiple-Cell StackIs an extension ofThe 1-cell stack
calculation
1/30/2003
16
Steady State: 16-Cell Stack ModelFuel Delivery: 8E-6 kg/s/cell @ 944KAir Delivery: 0.25 kg/s/cell @ 944K
Output: 245 mW/cm2
Tcell(ave) = 751C
Full 3-D Temperature dataset available for computing thermal stress
1/30/2003
17
Topics for Microstructural ElectrochemistryMethod.Advantages.Sample simulation results.• Reaction zones in the anode• Internal reformation
1/30/2003
18
Effective Property MethodDiscretize gas channels, electrodes and electrolyte into nodes with 5-25 µm thicknessEach node has effective transport and reaction properties (gas diffusion, surface diffusion, surface area density, TPB density, etc.)Solve flow and transport equations using lattice Boltzmann to obtain three-dimensional distributions for• Gas velocity, density• Gas species concentrations• Adsorbed species, solid diffusing species (oxygen)• Energy, temperature
1/30/2003
19
Microstructural Electrochemistry MethodGeometry may be generated using statistical data taken from digitized pictures of the porous materialProperties include effective gas diffusion, surface diffusion, solid diffusion, triple phase boundary density, etc.Effective properties may be determined using lattice Boltzmannsimulations of the discrete microstructure
1/30/2003
20
Advantages
Spatially varying electrode propertiesParallel transport paths for oxygen (gas, surface, solid) Distributed reaction zones (TPB, internal reformation)Link cell performance to electrode microstructure
1/30/2003
21
Simulation Results from Effective Property Model
Fuel Channel
Anode
Cathode
Air Channel
800 nodes
∆ X = 10µm5 10 -6
5.2 10 -6
5.4 10 -6
5.6 10 -6
5.8 10 -6
6 10 -6
0.04 0.06 0.08 0.1 0.12 0.14
Hyd
roge
n co
ncen
trat
ion
(mol
e/cm
3)
Location (cm.)
4 10-7
6 10-7
8 10-7
1 10-6
1.2 10-6
1.4 10-6
1.6 10-6
1.8 10-6
2 10-6
0.04 0.06 0.08 0.1 0.12 0.14
Wat
er C
once
ntra
tion
(mol
e/cm
3)
Location (cm.)
2.2 10-6
2.25 10-6
2.3 10-6
2.35 10-6
2.4 10-6
2.45 10-6
0.12 0.125 0.13 0.135 0.14
Oxy
gen
conc
entr
atio
n (m
ole/
cm3)
Location (cm.)
280nodes
1/30/2003
22
Internal ReformationMethane distribution
Methane is reformed upon contact withAnode surface
Methane is continually depletedGradient is driving it to surface
Hydrogen distribution
High diffusion into the channelAnd low diffusion into the anode
1/30/2003
23
Topics for Start-upStart-up issues and challenges.Computational tools for start-up simulations and stack design.Thermal controller for rapid heating of stacks.Structural parametric results.Optimization studiesExperimental validation of structural models.
1/30/2003
24
Rapid Start-up Issues for SOFC
• Flow through stack must be “uniformly” distributed.• Maintain thermal stresses within material set strength.• Stack pressure drop to be small.• Minimize time to heat stack to initiation temperature of 700 oC (within a few
minutes ultimately)• Issues necessitate survey of designs, with given material set, to discover
working options …Stack Geometry, Gas channel and manifold dimensions, flow configurations.
1/30/2003
25
Transient simulationsTarget Structure – Basic Planar Stack Design
Air inletmanifold
Air outletmanifold
Design Variables• Gas flow distribution channel
dimensions/types• Ceramic or metallic interconnect• Material thickness• PEN strength• Rigid or compression seals• Manifold dimensions• Flow configuration
Square, cross flow
1/30/2003
26
Transient - Rapid Start-upThermal-Fluids Stack Model
Air - Inflow manifold
Air - (porous media) 2.5mm channel ht
Fuel - Outflow manifoldMetallic or CeramicInterconnectPositive-Electrolyte-Negative(PEN) Composite in Model
Fluid cells removed
470,000 computational cells36 hours (x8 CPU) for 5 min
1/30/2003
27
Stack Model with Temperature Control, Modified Geometry and Boundary Conditions
Non-Uniform stack flow and heating…..
Would also mean non-uniform reactions and heating during steady operation.
Fix: increased outlet manifold dimension
• Flow channel heightshortened to (1.5mm)
• User routine defined free convection BC at walls ( Te)
• User routine defined temperature control
1/30/2003
28
Prediction of Temperature Distribution and Subsequent Thermal Stresses
Temperatures from CFD model Stresses in various ComponentsFrom FEA Models
Guidance for Modifications:• Heating Method• Flow Channels • Manifold Dimensions
Updated Geometry or Boundary Conditions.New CFD Prediction of Flow and Temperature
1/30/2003
29
Temperature ControllerTemperature Controller
0200400600800
100012001400
0 10 20 30
Time (min.)
Tem
pera
ture
(K)
OptimizedStandard
245724380640504523
Pen ∆TAve Pen
CFD model is run with standard controller to generate temperature distribution
FEA modeling performed using temperature distribution to determine maximum PEN ∆T allowed based upon strength of material
CFD model is re-run using optimized temperature controller to meet the maximum ∆T at each average PEN temperature
1/30/2003
30
Optimization studiesController for heat-up time optimizationGeometrical optimizationOptimization of mesh stiffness
Bounding (minimum) heat-up is the order of 10 minutes(heat rate of 19.6 KJ/sec and total heat input of 6.81 MJ)
Stamped 2-piece design
fold
1/30/2003
31
Effect of Temperature Profile and Seal Compliance on Stresses in the PEN
0
20
40
60
80
100
120
1 100 10000 1000000 1E+08 1E+10
k
1 M
pa
Linear - PENParabolic - PENLinear - stackParabolic - stack
1:1 100/10 Stack (18.4 MPa)
1:1 Linear Temp Stack(15.7 MPa)
RigidFree
1/30/2003
32
Effect of Temperature Profile and Seal Compliance on Stresses in the PEN
0102030405060708090
100
1 100 10000 1000000 1E+08 1E+10
k
σ1
Mpa Counter flow
Co flowCross flow
CROSS-FLOW CO-FLOW COUNTER-FLOW
Fuel
Air
FuelFuelA
ir
Air
Tave = 771°C∆T = 250°C
Tave = 761°C∆T = 183°C
Tave = 756°C∆T = 265°C
1/30/2003
33
Layered model results
0.06326.3Counter flow
0.03510.3Co - flow
0.03127.3Cross flow
VerticalDeflection
(mm)
Anode stress(MPa)
Design
PEN out-of-plane deflection
Anode principal stress (Pa)
1/30/2003
34
Results 3-Stack Simulations
Anode σ1 MPa
SS σeqv MPa
1% (A) Bottom 24 370Top 22 547
0.001% (D) Bottom 49.6 313Top 32.4 609
0.001% (E) Bottom 56.4 410Simple BC Top 33.8 585
0.001% (F) Bottom 52.5 32610% glass Top 28.9 617
Will the stack survive thermal stresses? (based on stress/strength failure criteria)
What is the effect of out of plane stiffness?
Will softer glass reduce stresses?How do the B.C.’s change stress
profiles?
1/30/2003
35
Experimental Validation of Structural Modeling
Rapid (<30 sec.) heating of ceramic PEN to 700°C with 20 KW infrared heaters. Temperature profile controlled with parabolic shaped mask
Infrared image of temperature profile
Finite element modeling of test
1/30/2003
36
Applicability to SOFC CommercializationModeling tools developed by PNNL for design, optimization and operation of SOFC materials, stacks and systems.Engineering insights and guidance regarding SOFC materials, stacks and systems.
1/30/2003
37
Future ActivitiesEnhancement of continuum level electrochemistry models for full stacks and steady state parametric studies.Micro-structural level electrochemistry for microstructural optimization and simulating internal reformationPredictive models for strength and lifeMaterial properties and model correlation/validation.