Evaluation and Evaluation and Design of Water- - Design of Water - - PowerPoint PPT Presentation

evaluation and evaluation and design of water design of
SMART_READER_LITE
LIVE PREVIEW

Evaluation and Evaluation and Design of Water- - Design of Water - - PowerPoint PPT Presentation

Evaluation and Evaluation and Design of Water- - Design of Water Splitting Cycles Splitting Cycles Scott Mullin, Jacob Tarver, Uchenna Uchenna Odi Odi Scott Mullin, Jacob Tarver, University of Oklahoma May 2006 Overview Overview


slide-1
SLIDE 1

Evaluation and Evaluation and Design of Water Design of Water-

  • Splitting Cycles

Splitting Cycles

Scott Mullin, Jacob Tarver, Scott Mullin, Jacob Tarver, Uchenna Uchenna Odi Odi

University of Oklahoma May 2006

slide-2
SLIDE 2

Overview Overview

  • Need for hydrogen

Need for hydrogen

  • Water

Water-

  • splitting cycles as solution

splitting cycles as solution

  • Current evaluation methods

Current evaluation methods

  • Efficiency defined

Efficiency defined

  • Our methodology as improvement

Our methodology as improvement

  • Results of our analysis

Results of our analysis

  • Economics

Economics

  • Conclusions

Conclusions

slide-3
SLIDE 3

Accomplishments Accomplishments

  • Novel methodology

Novel methodology

  • Rapidly screen cycles without detailed

Rapidly screen cycles without detailed process process flowsheets flowsheets

  • Optimize T, P and excess reactants for non

Optimize T, P and excess reactants for non-

  • spontaneous reactions

spontaneous reactions

  • Scoping algorithm

Scoping algorithm

  • Calculations refined for best cycles

Calculations refined for best cycles

  • Found better cycles than currently favored

Found better cycles than currently favored Sulfur Sulfur-

  • Iodine and UT

Iodine and UT-

  • 3

3

slide-4
SLIDE 4

Hydrogen Economy Hydrogen Economy

  • Currently 11 million tons/year

Currently 11 million tons/year

  • In H

In H2

2 economy

economy†

†:

:

  • 200 million tons/year for transportation

200 million tons/year for transportation

  • 450 million tons/year for all non

450 million tons/year for all non-

  • electric

electric

  • H

H2

2 is not a natural resource

is not a natural resource

  • Must be produced

Must be produced

  • Steam reformation of methane

Steam reformation of methane

  • CO

CO2

2 output

  • utput
  • Rising fuel prices

Rising fuel prices

† K. R. Schultz 2003, General Atomics, DOE grant

slide-5
SLIDE 5

Alternative H Alternative H2

2 Production

Production

Petroleum Petroleum

  • CO

CO2

2, expensive

, expensive

Electrolysis, high T electrolysis Electrolysis, high T electrolysis

  • Premature, inefficient

Premature, inefficient

Photocatalytic Photocatalytic reactors reactors

  • Premature

Premature

Thermochemical Thermochemical cycles cycles

  • Efficient, established processing techniques

Efficient, established processing techniques

slide-6
SLIDE 6

2 2 2

Abundant heat, electricity

slide-7
SLIDE 7

Water Water-

  • Splitting Cycles

Splitting Cycles

“New New” ” technology, chosen by DOE technology, chosen by DOE through Nuclear Hydrogen Initiative through Nuclear Hydrogen Initiative

  • Efficient hydrogen production

Efficient hydrogen production

  • 50

50-

  • 60% currently, 80

60% currently, 80-

  • 90%+ possible

90%+ possible

  • Use 950

Use 950º ºC or cooler process heat C or cooler process heat

  • 202 cycles known, but few researched

202 cycles known, but few researched

  • Others can be found, as described by

Others can be found, as described by Holiastos Holiastos and and Manousiouthakis Manousiouthakis 1998 1998

slide-8
SLIDE 8

Economics Economics

  • $1 billion for water

$1 billion for water-

  • splitting facility

splitting facility

  • $100 million range annual energy costs

$100 million range annual energy costs

  • Which cycle is best?

Which cycle is best?

  • Few cycles researched in detail

Few cycles researched in detail

  • Process design too complex

Process design too complex

  • Efficient cycles desirable

Efficient cycles desirable

  • Justify increased equipment costs

Justify increased equipment costs Bottom line: saving few % efficiency has huge Bottom line: saving few % efficiency has huge savings over plant lifetime savings over plant lifetime

slide-9
SLIDE 9

Cycles Cycles

  • Most are

Most are thermochemical thermochemical, some hybrid electric , some hybrid electric

  • Any number of reactions, species

Any number of reactions, species

  • Named after institutions or chemicals

Named after institutions or chemicals

  • Steady

Steady-

  • state operation

state operation

T1 O2 H2O A B, C T2 H2

Sample 2-step cycle

2

A B + C + O ⎯⎯ →

2 2

B + C + H O A + H ⎯⎯ →

T2 T1 T2 T1

slide-10
SLIDE 10

Efficiency Efficiency

  • Theoretical, 1 mol basis for cycle comparison

Theoretical, 1 mol basis for cycle comparison

  • Minimum reversible energy (heating and work)

Minimum reversible energy (heating and work) requirement requirement

  • Performance limit

Performance limit

  • Thermodynamics: JANAF tables for state

Thermodynamics: JANAF tables for state functions, pure component averages functions, pure component averages

f 2

H (H O) Q W η Δ = +

  • Q is total heat requirement

W is separation, electric and shaft work†

†Shaft work (pumping, compression) small compared to other terms

slide-11
SLIDE 11

Previous Surveys Previous Surveys

  • Brown et al 2000 scored cycles based on

Brown et al 2000 scored cycles based on known characteristics known characteristics

  • Good starting point, but not reproducible

Good starting point, but not reproducible

  • Arbitrary criteria, no emphasis on efficiency

Arbitrary criteria, no emphasis on efficiency

  • Elemental abundance,

Elemental abundance, “ “corrosivity corrosivity” ”, # elements , # elements

  • Rejects cycles with

Rejects cycles with “ “too positive too positive” ” free energies free energies

  • Favors well

Favors well-

  • researched cycles

researched cycles

slide-12
SLIDE 12

Score Score†

1 2 3 # reactions 6

  • 5

# separations 10 9 8 7 # elements 7

  • 6
  • Least

abundant element Ir Rh, Tc, Os, Ru, Re, Au Pt, Bi, Pt, Bi, Pd, Hg, Pd, Hg, Se Se Ag, In, Ag, In, Cd Cd, , Sb Sb, , Tm, Tm, Tl Tl, Lu , Lu

† †Adapted from Brown et al 2000

Adapted from Brown et al 2000

Brown’s method is good at identifying cycles based on estimated process complexities, but is not quantitative or reproducible. What happens if you change the weights, or add further scoring criteria?

slide-13
SLIDE 13

Previous Surveys cont Previous Surveys cont’ ’d d

  • Cycles are complex, so Lewis et al 2005

Cycles are complex, so Lewis et al 2005 developed systematic approach developed systematic approach

  • Scoping method based on efficiency

Scoping method based on efficiency

  • Quantitative, standard basis

Quantitative, standard basis

  • Oversimplifications

Oversimplifications

  • Requires detailed

Requires detailed flowsheets flowsheets

  • Not truly scoping

Not truly scoping

  • Assumes 50% loss of all work energy

Assumes 50% loss of all work energy

  • Does not estimate real separation energy

Does not estimate real separation energy

Our method is truly scoping, based on Our method is truly scoping, based on theoretical requirements theoretical requirements

slide-14
SLIDE 14

General Methodology General Methodology

  • Cyclic nature couples all calculations

Cyclic nature couples all calculations

  • Decouple the problem

Decouple the problem

  • Find realistic estimates for Q, W

Find realistic estimates for Q, W

  • Refine calculations for best cycles

Refine calculations for best cycles

  • Account for additional energy requirements

Account for additional energy requirements

  • Economic analysis of best cycles

Economic analysis of best cycles

  • Apply methodology to all cycles

Apply methodology to all cycles

  • Evaluate the 202 from literature

Evaluate the 202 from literature

  • Find unknown cycles

Find unknown cycles

slide-15
SLIDE 15

Equilibrium Equilibrium

  • Excess reactants added to shift reactions to

Excess reactants added to shift reactions to the right the right

  • How do we handle excess after the

How do we handle excess after the reaction? reaction?

  • Requires optimization, coupled equations

Requires optimization, coupled equations

products reactants

i i

v i i eq x v i i

n n K K P P n n

ν ν

⎛ ⎞ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ ∑ ∑ = = ⎛ ⎞ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠

∏ ∑ ∏ ∑

slide-16
SLIDE 16

Excess Reactant Handling Excess Reactant Handling

No recycle: saves separation energy, but negatively shifts equilibrium in most cases and increases heat cascade requirement Immediate recycle: full separation energy costs

We optimize T, P, # excess mols and their handling

2

A B + C + O ⎯⎯ →

2 2

B + C + H O A + H ⎯⎯ →

T2 T1 T2 T1

slide-17
SLIDE 17

Cycles cont Cycles cont’ ’d d

  • Methodology

Methodology accounts for arbitrarily accounts for arbitrarily complex cycles complex cycles

T1 T4 T2 T3 H2O O2 H2 A B C D, H2O E, F

T4 T3 T2 T1 T4 T3 T2 T1

2

A + B C + O ⎯⎯ →

2 2

D + H O E + F + H ⎯⎯ →

2

C + H O B + F ⎯⎯ →

2

E + F A + D + H O ⎯⎯ →

Conditions optimized for each reactor

slide-18
SLIDE 18

Heat Requirements Heat Requirements

  • Maximize heat recovery from exothermic reactions and cooling

Maximize heat recovery from exothermic reactions and cooling streams streams

  • Pinch occurs when there is not enough heat to power reactions or

Pinch occurs when there is not enough heat to power reactions or heat streams, requiring input from the hot utility heat streams, requiring input from the hot utility

slide-19
SLIDE 19

Generic Heat Integration Generic Heat Integration

Hhot is total enthalpy of cooling streams Hcold is total enthalpy of heating streams

slide-20
SLIDE 20

Pinch Point and Approach Temp. Pinch Point and Approach Temp.

Heat is added above the pinch. Heat transfer over the pinch (greater than the minimum heat requirement) goes to cold utility and is wasted. ΔTmin is closest feasible temperature, since complete heat transfer requires infinite exchanger area.

slide-21
SLIDE 21

Heat Integration Method Heat Integration Method†

  • Zonal analysis

Zonal analysis

  • Approach temperature

Approach temperature

  • Simplifying algorithm

Simplifying algorithm

  • Keep track of total

Keep track of total heat usage, advancing heat usage, advancing to successive zones to successive zones and reactors and reactors

  • Cold utility ignored

Cold utility ignored

  • Leftover heat

Leftover heat sometimes useful for sometimes useful for electricity generation electricity generation

† PT&W Plant Design and Economics for Chemical Engineers

slide-22
SLIDE 22

Electrical Work Electrical Work

  • Nernst equation for electrolytic cells

Nernst equation for electrolytic cells

  • Assume steady

Assume steady-

  • state operation of electrolytic cells

state operation of electrolytic cells

  • New electrolysis methods efficient compared to batch process

New electrolysis methods efficient compared to batch process†

  • Hybrid cycles treated same in heat integration

Hybrid cycles treated same in heat integration

elec

W nFE = −

  • T

(298) 298

d(E (T)) E E dT = +∫

  • †Motupally et al 1998
slide-23
SLIDE 23

Separation Work Separation Work

  • Minimum separation estimate

Minimum separation estimate

  • Phases self

Phases self-

  • separate

separate

  • We don

We don’ ’t pay for it t pay for it

  • Estimate separation efficiencies

Estimate separation efficiencies

This provides us with a minimum requirement. Chemical mixing and individual processes will increase W. Assign efficiencies to each process: e.g. assume distillation columns 50% efficient

, sep ideal separation

W W η =

ln Assuming isothermal separation ln ln

sep i i i i i i sep i i i i i i

  • ut

in

W G n RT n x W RT n x n x μ

= = = =

= −Δ = −Δ = − Δ ⎡ ⎤ ⎛ ⎞ ⎛ ⎞ = − ⎢ ⎥ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ ⎝ ⎠ ⎣ ⎦

∑ ∑ ∑ ∑

slide-24
SLIDE 24

Solution Procedure Solution Procedure

  • Most reactions go to completion

Most reactions go to completion

  • No excess reactants to handle

No excess reactants to handle

  • Optimize reactors individually

Optimize reactors individually

  • For other reactions

For other reactions

  • Find equilibrium concentrations

Find equilibrium concentrations

  • Newton method to solve for conversion

Newton method to solve for conversion

  • Know how much product we need from

Know how much product we need from connectivity connectivity

slide-25
SLIDE 25

Solution Procedure cont Solution Procedure cont’ ’d d

  • Computer algorithm finds optimum

Computer algorithm finds optimum efficiency for each T efficiency for each T

  • P easy to find

P easy to find

  • Finds Q and W for each #

Finds Q and W for each # mols mols excess excess

  • Optimize these for each recycle scheme

Optimize these for each recycle scheme

  • Computer crawls through solutions, and

Computer crawls through solutions, and maximizes efficiency maximizes efficiency

slide-26
SLIDE 26

Example Optimization Example Optimization

Cl2 (g) + H2O (g) -> 2HCl (g) + ½O2 (g), ∆Grxn= -17 kJ / cycle mol

Wmin and Excess Cl2 Required for varying excess H2O

35 35.5 36 36.5 37 37.5 38 38.5 39 39.5 0.2 0.3 0.4 0.48227768 0.6 0.7 0.8

Excess H2O (moles) Wmin Separation (kJ/mol)

0.2 0.4 0.6 0.8 1 1.2

Excess Cl2 Required (moles)

Wmin CL2

slide-27
SLIDE 27

Sample Sample Thermochemical Thermochemical Cycles Cycles

  • Julich

Julich

  • Ispra

Ispra Mark 9 Mark 9

T = 1073 K 1 3 4(s) 4(s) 2 3(s) 2(g) 2(g) 2 T = 973 K (s) 2 (g) 3 4(s) 2(g) T = 473 K 2 3(s) 2(g) (s) 4(s)

Fe O + 3FeSO 3Fe O + 3SO + O 3FeO + H O Fe O + H 3Fe O + 3SO 3FeO + 3FeSO ⎯⎯⎯⎯ → ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

T = 923 K 2(s) 2 (g) 3 4(s) (g) 2(g) T = 693 K 3 1 2(g) 3 4(s) (g) 3(l) 2 (g) 2(g) 2 2 T = 423 K 3(s)

3FeCl + 4H O Fe O + 6HCl + H Cl + Fe O + 6HCl 3FeCl + 3H O + O 3FeCl ⎯⎯⎯⎯ → ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

3 2(g) 2(s) 2

Cl + 3FeCl

slide-28
SLIDE 28

Sample Sample Thermochemical Thermochemical Cycles Cycles

  • Sulfur Iodine

Sulfur Iodine

  • US

US-

  • Chlorine

Chlorine

T = 1123 K 1 2 4(g) 2(g) 2 (g) 2(g) 2 T = 573 K (g) 2(g) 2(g) T = 473 K 2(g) 2(g) 2 (l) (g) 2 4(g)

H SO SO + H O + O 2HI I + H I + SO + 2H O 2HI + 1H SO ⎯⎯⎯⎯ → ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

T = 1123 K 1 2(g) 2 (g) (g) 2(g) 2 T = 773 K 2(s) (l) 2(g) T = 473 K (s) (g) 2(s) 2(g)

Cl + H O HCl + O 2CuCl 2CuCl + Cl 2CuCl + 2HCl 2CuCl + H ⎯⎯⎯⎯ → ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

slide-29
SLIDE 29

Sample Sample Thermochemical Thermochemical Cycles Cycles

  • Gaz

Gaz de France de France

  • UT

UT-

  • 3 Tokyo

3 Tokyo

T = 1098 K 2 (s) (g) 2 2(s) T = 998 K (l) (l) 2 (s) 2(g) T = 398 K 1 2 2(s) 2 (l) (s) 2(g) 2

2K O K + K O 2K + 2KOH 2K O + H K O + H O KOH + O ⎯⎯⎯⎯ → ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

T = 1023 K 2(l) 2 (g) (s) (g) T = 873 K 2 (g) 2(s) 2(g) (s) 1 2(s) 2(g) 3 4(s) (g) 2(g) 2 T = 573 K 3 4(s) (g)

CaBr + H O CaO + HBr 4H O + 3FeBr + Br + CaO CaBr + O + Fe O + HBr + H Fe O + 8HBr Br ⎯⎯⎯⎯ → ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

2(g) 2(s) 2 (g)

+ 3FeBr + 4H O

slide-30
SLIDE 30

Sample Sample Thermochemical Thermochemical Cycles Cycles

  • Ispra

Ispra Mark 4 Mark 4

  • Ispra

Ispra Mark 7A Mark 7A

T = 1123 K 1 2(g) 2 (g) (g) 2(s) 2 T = 1073 K 2 (g) (g) 2(g) T = 693 K 3(l) 2(g) 2(s) 2(s)

Cl + H O 2HCl + O H S S + H 2FeCl Cl + 2FeCl 2FeCl ⎯⎯⎯⎯ → ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

T = 473 K (g) (l) 3(s) 2 (g)

+ 2HCl + S 2FeCl + H S ⎯⎯⎯⎯ →

T = 1273 K 3 3 1 2(g) 2 3(s) 3(l) 2(g) 2 2 4 T = 923 K 2(s) 2 (g) 3 4(s) (g) 2(g) T = 693 K 3 3(l) 2(g) 2(s) 2 1 3 4(s) 2(g 4

Cl + Fe O FeCl + O 3FeCl + 4H O Fe O + 6HCl + H 3FeCl Cl + 3FeCl Fe O + O ⎯⎯⎯⎯ → ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

T = 623 K 3 ) 2 3(s) 2 T = 393 K 2 3(s) (g) 3(l) 2 (l)

Fe O Fe O + HCl 2FeCl + 3H O ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

slide-31
SLIDE 31

Sample Sample Thermochemical Thermochemical Cycles Cycles

  • Ispra

Ispra Mark 7B Mark 7B

T = 1273 K 9 3 9 2(g) 2 3(s) 3(l) 2(g) 2 2 4 T = 923 K 2(s) 2 (g) 3 4(s) (g) 2(g) T = 693 K 3 3(l) 2(g) 2(s) 2 3 (g) 2( 2

Cl + Fe O 3FeCl + O 3FeCl + 4H O Fe O + 6HCl + H 3FeCl Cl + 3FeCl 6HCl + O ⎯⎯⎯⎯ → ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

T = 673 K g) 2(g) 2 (g) T = 623 K 3 1 3 4(s) 2(g) 2 3(s) 4 2

3Cl + 3H O Fe O + O Fe O ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

slide-32
SLIDE 32

Sample Hybrid Cycles Sample Hybrid Cycles

  • Westinghouse

Westinghouse

  • Ispra

Ispra Mark 13 Mark 13

  • Hallett

Hallett Air Products Air Products

T = 1123 K 1 2 4(g) 2(g) 2 (g) 2(g) 2 T = 350 K 2(g) 2 (l) 2 4(g) 2(g)

H SO SO + H O + O SO + 2H O H SO + H ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

T = 1123 K 1 2 4(g) 2(g) 2 (g) 2(g) 2 T = 350 K (aq) 2(aq) 2(g) T = 350 K 2(l) 2(g) 2 (l) (g) 2 4(g)

H SO SO + H O + O 2HBr Br + H Br + SO + 2H O HBr + H SO ⎯⎯⎯⎯ → ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

T = 1073 K 1 2(g) 2 (g) (g) 2(g) 2 T = 298 K (g) 2(g) 2(g)

Cl + H O 2HCl + O 2HCl Cl + H ⎯⎯⎯⎯ → ⎯⎯⎯⎯ →

slide-33
SLIDE 33

Results Results

  • Cycle rankings based on Q

Cycle rankings based on QH

H analysis with

analysis with Δ ΔT Tmin

min=0

=0

  • 1. Hallett Air Products

1. US-Chlorine

  • 1. Sulfur Iodine
  • 1. Ispra Mark 13
  • 1. Westinghouse
  • 2. Ispra Mark 9
  • 3. Ispra Mark 4
  • 4. Gaz de France
  • 5. UT-3 Tokyo
  • 6. Julich
  • 7. Ispra Mark 7B
  • 8. Ispra Mark 7A
slide-34
SLIDE 34

Q QH

H analysis with

analysis with Δ ΔT Tmin

min=0

=0

Cycle Efficiencies using Qh for ΔTmin=0

52.3% 52.4% 54.8% 54.9% 75.0% 78.2% 85.7% 100.0% 100.0% 100.0% 100.0% 100.0%

0.0% 10.0% 20.0% 30.0% 40.0% 50.0% 60.0% 70.0% 80.0% 90.0% 100.0% Ispra Mark 7A Ispra Mark 7B Julich UT3 Tokyo Gaz de France Ispra Mark 4 Ispra Mark 9 Westinghouse Ispra Mark 13 Sulfur Iodine US-Chlorine Hallett Air Products

Cycle Efficiency

slide-35
SLIDE 35

Results cont. Results cont.

  • Now we consider

Now we consider W Wsep

sep, , stoich stoich and

and W Welec

elec as well

as well

QH only QH, Wsep, stoich, and Welec Note: arrows indicate only cycles that change 3+ positions

  • 1. Hallett Air Products

1. US-Chlorine

  • 1. Sulfur Iodine
  • 1. Ispra Mark 13
  • 1. Westinghouse
  • 2. Ispra Mark 9
  • 3. Ispra Mark 4
  • 4. Gaz de France
  • 5. UT-3 Tokyo
  • 6. Julich
  • 7. Ispra Mark 7B
  • 8. Ispra Mark 7A

1. US-Chlorine

  • 2. Sulfur Iodine
  • 3. Westinghouse
  • 4. Ispra Mark 9
  • 5. Gaz de France
  • 6. Ispra Mark 4
  • 7. Ispra Mark 13
  • 8. Julich
  • 9. Hallett Air Products
  • 10. UT-3 Tokyo
  • 11. Ispra Mark 7A
  • 12. Ispra Mark 7B
slide-36
SLIDE 36

Q QH

H,

, W Welec

elec, and

, and stoichiometric stoichiometric separation analysis with separation analysis with Δ ΔT Tmin

min=0

=0

47.9% 49.4% 49.6% 51 .1 % 51 .9% 55.7% 73.4% 75.0% 78.6% 85.1 % 88.1 % 96.1 %

0.0% 10.0% 20.0% 30.0% 40.0% 50.0% 60.0% 70.0% 80.0% 90.0% 100.0% Ispra Mark 7B Ispra Mark 7A UT3 Tokyo Hallett Air Products Julich Ispra Mark 13 Ispra Mark 4 Gaz de France Ispra Mark 9 Westinghouse Sulfur Iodine US-Chlorine

Efficiency Cycle Efficiencies using Qh, Welec, and Wsep, stoich for ΔTmin=0

slide-37
SLIDE 37

Results cont Results cont’ ’d d

  • Now we substitute

Now we substitute W Wsep

sep, , stoich stoich with

with W Wsep

sep, excess , excess

QH, Wsep, stoich, and Welec QH, Wsep, excess, and Welec 1. US-Chlorine

  • 2. Sulfur Iodine
  • 3. Westinghouse
  • 4. Ispra Mark 9
  • 5. Gaz de France
  • 6. Ispra Mark 4
  • 7. Ispra Mark 13
  • 8. Julich
  • 9. Hallett Air Products
  • 10. UT-3 Tokyo
  • 11. Ispra Mark 7A
  • 12. Ispra Mark 7B

1. Westinghouse

  • 2. Gaz de France
  • 3. US-Chlorine
  • 4. Sulfur Iodine
  • 5. Ispra Mark 13
  • 6. Ispra Mark 9
  • 7. Julich
  • 8. Hallett Air Products
  • 9. Ispra Mark 7A
  • 10. Ispra Mark 4
  • 11. Ispra Mark 7B
  • 12. UT-3 Tokyo
slide-38
SLIDE 38

Q QH

H,

, W Welec

elec, and

, and excess separation excess separation analysis with analysis with Δ ΔT Tmin

min=0

=0

33.3% 34.0% 38.9% 39.8% 48.9% 49.4% 52.1% 53.0% 55.2% 60.9% 75.0% 85.1%

0.0% 10.0% 20.0% 30.0% 40.0% 50.0% 60.0% 70.0% 80.0% 90.0% 100.0% UT3 Tokyo Ispra Mark 7B Ispra Mark 4 Ispra Mark 7A Hallett Air Products Julich Ispra Mark 9 Ispra Mark 13 Sulfur Iodine US-Chlorine Gaz de France Westinghouse

Efficiency Cycle Efficiencies using Qh, Welec, and Wsep, excess for ΔTmin=0

slide-39
SLIDE 39

Top 6 Top 6 Thermochemical Thermochemical Cycles Cycles

  • Based upon full analysis at

Based upon full analysis at Δ ΔT Tmin

min=0

=0

  • What about

What about Δ ΔT Tmin

min > 0?

> 0?

  • Some cycles more sensitive

Some cycles more sensitive

1. Westinghouse 2. Gaz de France 3. US-Chlorine 4. Sulfur Iodine 5. Ispra Mark 13 6. Ispra Mark 9

slide-40
SLIDE 40

Effect of Effect of Δ ΔT Tmin

min on Q

  • n QH

H

Qh vs ΔTmin for Top 6 Cycles

285 305 325 345 365 385 405 5 10 15 20 25

ΔTmin (K) Qh (kJ / cycle mol)

Ispra Mark 9 Ispra Mark 13 Sulfur Iodine US-Chlorine Gaz de France Westinghouse

slide-41
SLIDE 41

Corresponding Efficiency Corresponding Efficiency

Cycle Efficiencies using Qh vs ΔTmin for Top 6 Cycles

75.0% 80.0% 85.0% 90.0% 95.0% 100.0% 5 10 15 20 25

ΔTmin (K) Cycle Efficiency

Ispra Mark 9 Ispra Mark 13 Sulfur Iodine US-Chlorine Gaz de France Westinghouse

slide-42
SLIDE 42

Effect of Effect of Δ ΔT Tmin

min on

  • n Q

QH

H+W

+Welec

elec+W

+Wsep

sep, , stoich stoich

Qh + Welec + Wsep, stoich vs ΔTmin for Top 6 Cycles

285.0 335.0 385.0 435.0 485.0 535.0 5 10 15 20 25

ΔTmin (K) Qh + Wsep, stoich (kJ / cycle mol)

Ispra Mark 9 Ispra Mark 13 Sulfur Iodine US-Chlorine Gaz de France Westinghouse

slide-43
SLIDE 43

Corresponding Efficiency Corresponding Efficiency

Cycle Efficiencies using Qh + Welec + Wsep, stoich vs ΔTmin for Top 6 Cycles

50.0% 55.0% 60.0% 65.0% 70.0% 75.0% 80.0% 85.0% 90.0% 95.0% 100.0% 5 10 15 20 25

ΔTmin (K) Cycle Efficiency

Julich US-Chlorine Ispra Mark 4 Sulfur Iodine Ispra Mark 7a Gaz de France

slide-44
SLIDE 44

W Wsep

sep, , stoich stoich vs.

  • vs. W

Wsep

sep, excess , excess

Comparison of Wsep, stoich and Wsep, excess for Top 6 Cycles

215.2 43.3 232.0 183.7 0.0 11.7 30.3 17.2 38.7 11.7 0.0 11.7

0.0 50.0 100.0 150.0 200.0 250.0

Ispra Mark 9 Ispra Mark 13 Sulfur Iodine US-Chlorine Gaz de France Westinghouse

Wsep ( kJ / cycle mol)

Wsep (excess) Wsep (stoich)

slide-45
SLIDE 45

Effect of Effect of Δ ΔT Tmin

min on

  • n Q

QH

H+W

+Welec

elec+W

+Wsep

sep, excess , excess

Qh + Welec + Wsep, excess vs ΔTmin for Top 6 Cycles

285.0 335.0 385.0 435.0 485.0 535.0 585.0 5 10 15 20 25

ΔTmin (K) Qh + Wsep, excess (kJ / cycle mol)

Ispra Mark 9 Ispra Mark 13 Sulfur Iodine US-Chlorine Gaz de France Westinghouse

slide-46
SLIDE 46

Corresponding Efficiency Corresponding Efficiency

Cycle Efficiencies using Qh + Welec + Wsep, excess vs ΔTmin for Top 6 Cycles

50.0% 55.0% 60.0% 65.0% 70.0% 75.0% 80.0% 85.0% 90.0% 95.0% 100.0% 5 10 15 20 25

ΔTmin (K)

Qh + Wsep, excess (kJ / cycle mol)

Ispra Mark 9 Ispra Mark 13 Sulfur Iodine US-Chlorine Gaz de France Westinghouse

slide-47
SLIDE 47

Efficiency: Literature Efficiency: Literature Comparison Comparison†

†Brown et al 2000 ‡10% additional efficiency projected with electricity co-generation

Reported Reported (thermal) (thermal) Theoretical Theoretical (thermal) (thermal) Theoretical Theoretical (heat/work) (heat/work)

Sulfur Sulfur-

  • Iodine

Iodine

52% 52%‡

100% 100% 55% 55%

Tokyo UT Tokyo UT-

  • 3

3

49% 49%‡

55% 55% 33% 33%

Westinghouse Westinghouse

50% 50% 100% 100% 85% 85%

slide-48
SLIDE 48

Good Cycle Characteristics Good Cycle Characteristics

  • Hottest reaction exothermic, cascades

Hottest reaction exothermic, cascades heat to power rest of the cycle heat to power rest of the cycle

  • Minimizes Q

Minimizes Q

  • Products phase separate from each other,

Products phase separate from each other, and from reactants and from reactants

  • Minimizes W

Minimizes W

  • No high T, P,

No high T, P, corrosivity corrosivity, etc. as described , etc. as described by Brown et al 2000 by Brown et al 2000

slide-49
SLIDE 49

Economic Methodology Economic Methodology

  • 500 ton/day production target

500 ton/day production target

  • Enough for 0.95 million cars, according to Schultz

Enough for 0.95 million cars, according to Schultz

  • Heat Integration

Heat Integration

  • Temperature intervals

Temperature intervals

  • Cascades

Cascades

  • Heat exchanger network

Heat exchanger network

  • Process Flow Diagrams

Process Flow Diagrams

  • Assumptions

Assumptions

  • Solids handling

Solids handling

  • Capital cost

Capital cost

slide-50
SLIDE 50

Westinghouse Cycle Westinghouse Cycle -

  • Heat Profile

Heat Profile

slide-51
SLIDE 51

Westinghouse Cycle Westinghouse Cycle -

  • Heat Profile

Heat Profile

slide-52
SLIDE 52

Westinghouse Westinghouse -

  • Heat Cascade

Heat Cascade

QH

  • 288.0 kJ

H2O + SO2

  • 110.3 kJ

∆Hrxn 1 184.8 kJ 7 exchangers O2

  • 13.5 kJ

H2

  • 1.5 kJ

∆Hrxn 2 129.5 kJ H2SO4 94.3 kJ H2O 4.7 kJ

  • 1.5 kJ
  • 3.2 kJ
  • 10.3 kJ
  • 84.0 kJ
  • 26.3 kJ
  • 103.2 kJ
  • 184.8 kJ
slide-53
SLIDE 53

Westinghouse Westinghouse -

  • Heat Exchanger

Heat Exchanger Network Network

H2SO4(g) SO2(g) + H2O(g) + 0.5 O2(g) SO2(g) + 2 H2O(l) H2SO4(g)+H2(g) 1073 K 350 K 298 K H2O Zone 2 Zone 1 O2 SO2, H2O H2 ∆Hrxn 1 ∆Hrxn 2 = -110.3 kJ = 94.3 kJ = -28.8 kJ = 184.8 kJ = 129.5 kJ = 2.4 kJ

508 298 HX-2 HX-2

  • 3.2 kJ

318 298 HX-1 318

HX-1

350 298

  • 1.5 kJ

350

P-51

HX-3

H2SO4

HX-3

  • 10.3 kJ

350 455 1173 508 HX-4 HX-4

  • 84.0 kJ

455 1173 1173 500 HX-5

  • 26.3 kJ

500 350 HX-6

  • 184.8 kJ

QH HX-7

  • 103.2 kJ

QH

slide-54
SLIDE 54

Westinghouse Westinghouse -

  • Heat Exchanger

Heat Exchanger Network Network

∆Hrxn 1 O2 H2 ∆Hrxn 2 H2SO4 H2O 350 K 1173 K 298 K 350 K 350 K 1173 K 318 K

HX-1

298 K

HX-2 HX-3

455 K H2O + SO2 1173 K 1173 K 500 K

HX-5

350 K 508 K

HX-2

298 K 350 K

HX-6 HX-7

slide-55
SLIDE 55

Westinghouse Westinghouse -

  • Process Flow

Process Flow Diagram Diagram

Reactor 1 1173 K Separator H2SO4(g) SO2(g), H2O(g), O2(g) O2(g) SO2(g), H2O(g)

HX-4

SO2(g), H2O(l) H2O(l) H2(g) H2SO4(l)

HX-3

O2(g) Electrolyzer 350 K

HX-1

H2(g) H2O(l)

HX-2

O2(g) H2O(l)

HX-5

Electrolyzer Heat Electrolyzer Heat

HX-7

Reactor Heat Reactor Heat

HX-6

slide-56
SLIDE 56

Handling Solids Handling Solids

  • Physical transport of solids difficult

Physical transport of solids difficult

  • Grinders necessary

Grinders necessary

  • Slow heat transfer between solids

Slow heat transfer between solids

  • Use sweep gas as intermediate heat carrier

Use sweep gas as intermediate heat carrier

  • Solid separations

Solid separations

  • Usually oxides and halide salts

Usually oxides and halide salts – – solvent separation solvent separation

slide-57
SLIDE 57

UT UT-

  • 3 University of Tokyo

3 University of Tokyo†

  • Solids do not move – reactors run in parallel batch
  • Preserves efficiency, but increases capital costs and instability
  • Reported thermal efficiency 49%, compared to 55% theoretical

3FeBr2 + 4H2O → Fe3O4 + 6HBr + H2 CaO + Br2 → CaBr2 + 0.5O2 CaBr2 + H2O → CaO + 2HBr

Membrane

Fe3O4 + 8HBr → 3FeBr2 + 4H2O + Br2

Membrane

O2 H2

P-15

H2O

†Adapted from Brown et al 2000

slide-58
SLIDE 58

US Chlorine US Chlorine – – Heat Cascade Heat Cascade

1123 K 773 K 473 K 298 K HCl O2 CuCl Cl2 H2O CuCl2 H2 2CuCl(s) + 2HCl(aq) 2CuCl2(s) + H2(v) Cl2(v) + H2O(v) 2HCl(v) + 1/2O2(v) 2CuCl2(s) 2CuCl(l) + Cl2(v)

slide-59
SLIDE 59

US Chlorine US Chlorine – – Process Flow Process Flow Diagram Diagram

slide-60
SLIDE 60

Gaz Gaz de France de France -

  • Heat Exchange

Heat Exchange Network Network

O2 ∆Hrxn 3 K H2 K2O ∆Hrxn 2 H2O K2O2 KOH ∆Hrxn 1 1098 K 1098 K

HX-1

998 K 1098 K 1040 K 998 K

HX-2

1030 K

HX-3

398 K 621 K 998 K 1098 K

HX-4

998 K 398 K 398 K

HX-5

298 K 398 K

HX-6

998 K

HX-7

298 K 398 K

HX-8

298 K

slide-61
SLIDE 61

Gaz Gaz de France de France

1098K Reactor 1 K2O(s) K(g) K2O2(s)

HX-1 HX-4

998K Reactor 2 KOH(l) H2(g) K2O(S) 998K 398K Reactor 3 K2O2(s) H2O(l) 298K KOH(s) O2(g)

HX-8

O2 (g) 298K K(g)

HX-7

H2 (g) 298K

HX-3

KOH 621K K(g) K(l)

HX-2 Nuclear Reactor HX-6 HX-5

H2O

slide-62
SLIDE 62

Capital Cost Capital Cost

  • New technology

New technology

  • Processes involve highly corrosive

Processes involve highly corrosive materials and high temperatures materials and high temperatures†

  • Resistance to degradation involved within the

Resistance to degradation involved within the cycles cycles

  • High temperature quality material required

High temperature quality material required

  • Research involved for design

Research involved for design

  • Some kinetics are currently unknown

Some kinetics are currently unknown

  • Contract work involved

Contract work involved

†Perret et al 2004

slide-63
SLIDE 63

Capital Cost cont Capital Cost cont’ ’d d

  • 500 tons/day hydrogen production

500 tons/day hydrogen production

  • Equilibrium (complete reaction)

Equilibrium (complete reaction)

  • Maximum heat exchange area possible

Maximum heat exchange area possible

  • Highly corrosive materials

Highly corrosive materials

  • Scale up has never been done

Scale up has never been done

slide-64
SLIDE 64

Capital Cost Results Capital Cost Results

Westinghouse Gaz de France US-Chlorine Efficiency 85% 75% 60% FCI $3,100,000,000 $6,200,000,000 $3,100,000,000 Energy Cost $27,000,000 $39,000,000 $38,000,000 Process Cost, $/lb H2 produced $0.07 $0.11 $0.11

slide-65
SLIDE 65

Conclusions Conclusions

  • Scoping methodology can screen large

Scoping methodology can screen large number of cycles with reasonable number of cycles with reasonable accuracy accuracy

  • Sulfur

Sulfur-

  • Iodine and other popular cycles are

Iodine and other popular cycles are not necessarily best not necessarily best

  • Find cycles with phase separations and

Find cycles with phase separations and good heat cascade good heat cascade

slide-66
SLIDE 66

Questions? Questions?

slide-67
SLIDE 67

References References

Brown, L.C.; Showalter, S.K.; Funk, J.F.; Brown, L.C.; Showalter, S.K.; Funk, J.F.; Nuclear Production of Hydrogen Using Nuclear Production of Hydrogen Using Thermochemical Thermochemical Water Water-

  • Splitting Cycles

Splitting Cycles. 2000. US DOE project under NERI grant . 2000. US DOE project under NERI grant DE DE-

  • FG03

FG03-

  • 99SF21888

99SF21888 Holiastos Holiastos, K.; , K.; Manousiouthakis Manousiouthakis, V.; , V.; Automatic Synthesis of Thermodynamically Feasible Automatic Synthesis of Thermodynamically Feasible Reaction Clusters Reaction Clusters AIChE AIChE Journal Journal, Vol. 44, No. 1, January 1998 pp. 164 , Vol. 44, No. 1, January 1998 pp. 164-

  • 173

173 Lewis, M. A.; Lewis, M. A.; High High-

  • Temperature

Temperature Thermochemical Thermochemical Processes

  • Processes. DOE Hydrogen Program,

. DOE Hydrogen Program, FY 2005 Progress Report. September 2005. Argonne National Labo FY 2005 Progress Report. September 2005. Argonne National Laboratory ratory Motupally Motupally, S.; , S.; Mah Mah, D.T; , D.T; Freire Freire, F.J.; Weidner, J.W.; , F.J.; Weidner, J.W.; Recycling Chlorine from Hydrogen Recycling Chlorine from Hydrogen Chloride Chloride -

  • A new and economical electrolytic process

A new and economical electrolytic process. The Electrochemical Society . The Electrochemical Society Interface, 1998 Interface, 1998 Peters, M.; Peters, M.; Timmerhaus Timmerhaus, K. D.; West, R. E.; , K. D.; West, R. E.; Optimization Application: Pinch Technology Optimization Application: Pinch Technology Analysis Analysis Plant Design and Economics for Chemical Engineers Plant Design and Economics for Chemical Engineers. pp. 414 . pp. 414-

  • 433.

433. McGraw McGraw-

  • Hill, NYC, 2003

Hill, NYC, 2003 Schultz, K. R. Schultz, K. R. Use of the Modular Helium Reactor for Hydrogen Production Use of the Modular Helium Reactor for Hydrogen Production. General . General Atomics Report. September 2003 US DOE Grant No. DE Atomics Report. September 2003 US DOE Grant No. DE-

  • FG03

FG03-

  • 99SF21888

99SF21888