Studies on bi-directional hydrogen isotopes permeation through the first wall of a magnetic fusion power reactor

Similar documents
Isotope Effects on Tritium Mass Transfer in Li 17 Pb 83 Liquid Blanket

Implantation Energy Dependence on Deuterium Retention Behaviors for the Carbon Implanted Tungsten

Tritium Transport Modelling: first achievements on ITER Test Blanket Systems simulation and perspectives for DEMO Breeding Blanket

In-vessel Tritium Inventory in ITER Evaluated by Deuterium Retention of Carbon Dust

1 EX/P4-8. Hydrogen Concentration of Co-deposited Carbon Films Produced in the Vicinity of Local Island Divertor in Large Helical Device

Effects of Water Vapor on Tritium Release Behavior from Solid Breeder Materials

Comparison of deuterium retention for ion-irradiated and neutronirradiated

Tritium Transport and Corrosion Modeling in the Fluoride Salt-Cooled High-Temperature Reactor

SOME CONSIDERATION IN THE TRITIUM CONTROL DESIGN OF THE SOLID BREEDER BLANKET CONCEPTS

Helium effects on Tungsten surface morphology and Deuterium retention

Perspective on Fusion Energy

In-Vessel Tritium Inventory in Fusion DEMO Plant at JAERI

FLIBE ASSESSMENTS ABSTRACT

Tritium Control Using Carbon Outside of Core Stephen T Lam

Fusion: The Ultimate Energy Source for the 21 st Century and Beyond

Comparison of tungsten fuzz growth in Alcator C-Mod and linear plasma devices

Air exposure and sample storage time influence on hydrogen release from tungsten

Impurity accumulation in the main plasma and radiation processes in the divetor plasma of JT-60U

31704 Dynamic Monte Carlo modeling of hydrogen isotope. reactive-diffusive transport in porous graphite

Estimation of Tritium Recovery from a Molten Salt Blanket of FFHR

Steady-state diffusion is diffusion in which the concentration of the diffusing atoms at

Recent Development of LHD Experiment. O.Motojima for the LHD team National Institute for Fusion Science

Role and Challenges of Fusion Nuclear Science and Technology (FNST) toward DEMO

Tritium Fuel Cycle Safety

Elements of Strategy on Modelling Activities in the area of Test Blanket Systems

Molecular Dynamics Study of Plasma Surface Interactions for Mixed Materials

A SUPERCONDUCTING TOKAMAK FUSION TRANSMUTATION OF WASTE REACTOR

Yuntao, SONG ( ) and Satoshi NISHIO ( Japan Atomic Energy Research Institute

Formation of High-b ECH Plasma and Inward Particle Diffusion in RT-1

Measurement and modelling of hydrogen uptake and transport. Alan Turnbull

Gas Chromatography Separation of H 2 -D 2 -Ar Using Pd/K

Comparison of tungsten fuzz growth in Alcator C-Mod and linear plasma devices!

EX/C3-5Rb Relationship between particle and heat transport in JT-60U plasmas with internal transport barrier

Modelling of JT-60U Detached Divertor Plasma using SONIC code

The Spherical Tokamak as a Compact Fusion Reactor Concept

Influence of Amorphous Carbon Deposition on the Probability for Recombination of Neutral Oxygen Atoms on Aluminium Surfaces

PHYSICAL VAPOR DEPOSITION OF THIN FILMS

Numerical Simulation of Detritiation System for NIFS with Commercial Catalyst and Adsorbent

M. Rutigliano a, D. Santoro a, and M. Balat-Pichelin b

Effect of Spiral Microwave Antenna Configuration on the Production of Nano-crystalline Film by Chemical Sputtering in ECR Plasma

DETECTION OF ANOMALOUS ELEMENTS, X-RAY AND EXCESS HEAT INDUCED BY CONTINUOUS DIFFUSION OF DEUTERIUM THROUGH MULTI-LAYER CATHODE (Pd/CaO/Pd)

STEADY-STATE EXHAUST OF HELIUM ASH IN THE W-SHAPED DIVERTOR OF JT-60U

Elastic Recoil Detection Method using DT Neutrons for Hydrogen Isotope Analysis in Fusion Materials. Abstract

Dopant Diffusion. (1) Predeposition dopant gas. (2) Drive-in Turn off dopant gas. dose control. Doped Si region

Lecture 4. Ultrahigh Vacuum Science and Technology

Plasma Wall Interactions in Tokamak

Tritium Safety of Russian Test Blanket Module

Driving Mechanism of SOL Plasma Flow and Effects on the Divertor Performance in JT-60U

Noninductive Formation of Spherical Tokamak at 7 Times the Plasma Cutoff Density by Electron Bernstein Wave Heating and Current Drive on LATE

Comparison of Ion Internal Transport Barrier Formation between Hydrogen and Helium Dominated Plasmas )

Molecular Dynamics Simulation of Chemical Sputtering of Hydrogen Atom on Layer Structured Graphite

Preliminary Safety Analysis of CH HCSB TBM

ITER A/M/PMI Data Requirements and Management Strategy

DEMO Concept Development and Assessment of Relevant Technologies. Physics and Engineering Studies of the Advanced Divertor for a Fusion Reactor

Plasma and Fusion Research: Regular Articles Volume 10, (2015)

Isotope Exchange in High Heat-flux Components of the JET Tritium Neutral Beam Injector

1. Motivation power exhaust in JT-60SA tokamak. 2. Tool COREDIV code. 3. Operational scenarios of JT-60SA. 4. Results. 5.

Physics of the detached radiative divertor regime in DIII-D

Drift-Driven and Transport-Driven Plasma Flow Components in the Alcator C-Mod Boundary Layer

Experimental Facility to Study MHD effects at Very High Hartmann and Interaction parameters related to Indian Test Blanket Module for ITER

Plasma interactions with Be surfaces

The GDT-based fusion neutron source as a driver of subcritical nuclear fuel systems

DIVIMP simulation of W transport in the SOL of JET H-mode plasmas

Design window analysis of LHD-type Heliotron DEMO reactors

Reduction of Turbulence and Transport in the Alcator C-Mod Tokamak by Dilution of Deuterium Ions with Nitrogen and Neon Injection

Predicting the Rotation Profile in ITER

TRITIUM PRODUCTION IN PALLADIUM DEUTERIDE/HYDRIDE IN EVACUATED CHAMBER

Investigation of Water Fragments

Chemical Sputtering of Carbon Materials due to Combined Bombardment by Ions and Atomic Hydrogen

Cesium Dynamics and H - Density in the Extended Boundary Layer of Negative Hydrogen Ion Sources for Fusion

Tritium processing system for the ITER Li/V Blanket Test Module

Dynamic Monte-Carlo modeling of hydrogen isotope reactive diffusive transport in porous graphite

INTRODUCTION TO MAGNETIC NUCLEAR FUSION

Investigation of Intrinsic Rotation Dependencies in Alcator C-Mod

Yasuhiro Iwamura, Takehiko Itoh, Nobuaki Gotoh and Ichiro Toyoda

Preliminary system modeling for the EUROfusion Water Cooled Lithium Lead blanket

Toward the Realization of Fusion Energy

Multiple hydrogen trapping by vacancies: Its impact on defect dynamics and hydrogen retention in tungsten

Lecture 11. Etching Techniques Reading: Chapter 11. ECE Dr. Alan Doolittle

Lower Hybrid Current Drive Experiments on Alcator C-Mod: Comparison with Theory and Simulation

Influence of High Magnetic Field on Fusion Reactor Blanket Processes

Present status of fusion neutronics activity and comments to neutron diagnostics for TBM

Electron Emission from Nano and Micro Structured Materials for Plasma Applications

Thomas Schwarz-Selinger. Max-Planck-Institut for Plasmaphysics, Garching Material Science Division Reactive Plasma Processes

On Dust Particle Dynamics in Tokamak Edge Plasma

Analysis By Time-Of-Flight Secondary Ion Mass Spectroscopy or Nuclear Products In Hydrogen Penetration Through Palladium

ENEA Experience in PbLi Technologies

NanoEngineering of Hybrid Carbon Nanotube Metal Composite Materials for Hydrogen Storage Anders Nilsson

FUSION NEUTRONICS EXPERIMENTS AT FNG: ACHIEVEMENTS IN THE PAST 10 YEARS AND FUTURE PERSPECTIVES

Applicability of Laser-Induced Desorption Quadruple Mass Spectrometry (LID-QMS) for the Determination of Local Deuterium Retention

Divertor Requirements and Performance in ITER

Light-Induced Atom Desorption in Alkali Vapor Cells

Energy fluxes in plasmas for fabrication of nanostructured materials

Vacuum Pumps. Two general classes exist: Gas transfer physical removal of matter. Mechanical, diffusion, turbomolecular

Flux & Energy Dependence of Methane Production From Graphite Due to H + Impact

Vibrationally resolved ion-molecule collisions

Is plasma important? Influence molecule formation?

Permeation Measurement Testing Techniques. Michelle Stevens MOCON, Inc.

Integrated Particle Transport Simulation of NBI Plasmas in LHD )

The role of PMI in MFE/IFE common research

Transcription:

Studies on bi-directional hydrogen isotopes permeation through the first wall of a magnetic fusion power reactor IAEA-CRP Plasma-Wall Interaction with Reduced Activation Steel Surfaces in Fusion Devices Dec. 9-11 th 2015 Yoshi Hirooka National Institute for Fusion Science Graduate University of Advanced Studies Haishan Zhou 周海山

Table of Contents 1. Background and motivation 2. Experimental facility and setup 3. Hydrogen permeation experimental data and modelling 3-1. Implications for reactor fuel loop operations 3-2. Hydrogen transport parameters in the first wall 6. Summary 2

Definition of the first wall: Background (1) All the fusion experimental devices up to ITER: the first wall is a vacuum chamber wall to separate plasma from the environment. Power reactors: the first wall is the plasma-facing surfaces of breeding blanket units. First wall JET LHD The first walls of existing fusion devices. FFHR reactor A. Sagara et al,. research report, NIFS-MEMO-64 (2013). 3

Blanket concepts include: water-cooled solid breeder (e.g. Li 2 TiO 3 ), He-cooled solid breeder (e.g. Li 4 SiO 4 ), self-cooled liquid breeder (e.g. FLiBe), water-cooled liquid breeder (e.g. Li-Pb) Typical of the first wall of a fusion power reactor are: I. Large surface area: ~10 3 m 2 ; II. Made of reduced activation materials; III. High temperature operation, e.g. ~500 o C for ferritic steel alloys; IV. Thin wall design to reduce thermomechanical stress. Background (2) A. Sagara et al., J. Nucl. Mater. 258-263 (1998) 2079-2082. 4

Background (3) In most of the recent reactor studies, the first wall is designed to be 5 mm or even less, although these concepts employ various first wall materials. Thickness optimization [2] [2] A. Sagara et al., Fusion Technol. 39(2001)753-757. Y. Ueda et al. J. Nucl. Mater. 313-316(2003)32-41. 5

Potential Issues with the thin wall design First wall ~5mm Edge plasma Breeder/coolant D + + T + Plasma driven permeation(pdp) T 2 Gas-driven permeation(gdp) Potential issues, but not quite well addressed, with the thin wall design are: Plasma-driven permeation (PDP) of D/T fuel into the blanket side Necessitates isotope separation Gas-driven permeation (GDP) of bred T into the plasma side (i.e. gas puff) Leads to an unwanted increase in edge plasma density Affect core plasma confinement

Tritium equilibrium pressures in breeders Solubility of gases in metals (Sieverts' law): C = ST ( ) p T in wall T 2 Assuming 0.1 ppm, the total tritium inventory is of the order of 10 g in the FLiBe blanket. For blanket employing FLiBe, the tritium thermodynamic equilibrium pressure is ~10 4 Pa at 527 o C at a tritium concentration of ~0.1 ppm. S. Fukada, Y. Edao, Y. Maeda and T. Norimatsu, Fusion Eng. Des. 83, 747 (2008). 7

DIFFUSE-code predictions on FLiBe-blankets (PDP: 100eV, e16 H + /cm 2 5 mm α-fe, GDP: 10 4 Pa H 2 5 mm α-fe at 300 800K) Permeation fluxes Concentration profiles PDP and GDP fluxes (H-atoms/cm 2 /s) 1.0E+16 1.0E+14 1.0E+12 1.0E+10 1.0E+08 1.0E+06 1.0E+04 1.0E+02 1.0E+00 1.0E-02 1.0E-04 1.0E-06 1.0E-08 1.0E-10 1.0E-12 1.0E-14 1.0E-16 1.0E-18 1.0E-20 Temp. PDP flux PDP GDP Temperature(K) GDP~1Torr liter/s/m 2 0 200 400 600 800 1000 Time(s) 1000 900 800 700 600 500 400 300 200 100 0 Temperature (K) Hydrogen concentration (H-atoms fraction) 1.40E-05 1.20E-05 1.00E-05 8.00E-06 6.00E-06 4.00E-06 2.00E-06 GDP t=1.0e-03s GDP t=2.0e+02s GDP t=1.0e+03s 100XPDP t=1.0e- 03s 100XPDP t=2.0e+2s 100XPDP t=1.0e+3s 0.00E+00 0.00 0.10 0.20 0.30 0.40 0.50 Distance from the upstream surface (cm)

F82H: Hydrogen PDP and GDP fluxes Hydrogen PDP and GDP fluxes measured in VEHICLE-1 for a 5 mm thick F82H membrane (under the FLiBe-blanket conditions) Membrane temperature ~500 o C J + = D L J K 0 r Flow from blanket J= DT ( ) dc dx Flow from edge plasma Net Implantation flux for PDP H 2 driving pressure for GDP 2.0 10 16 H cm-2 s-1 (-100 V bias) ~10 4 Pa(Fukada) (Thermodynamic equilibrium data) Hydrogen transport in the first wall is dominated by the flow from the blanket. Using the experimental data, the hydrogen recycling rate has been estimated to be R= 1.025. 9

Evaluation of bi-directional permeation Conditions for DIFFUSE calculation: T 2 pressure: 1 10 4 Pa (Fukada) Implantation flux assumptions: Γ D = 5 10 15 D cm-2 s-1 Γ T = 5 10 15 T cm-2 s-1 Membrane: 5 mm thick α-fe Temperature: 550 o C D/T release fluxes at the plasma-facing surface T GDP flow Boundary conditions: Gas side: Sieverts law T Plasma side: recombination D T Intrinsic trap density: 1% Trapping energy: 0.62 ev The tritium release flux at the plasma-facing surface is a total flux from GDP and re-emission. 10

PDP and GDP under reactor-relevant conditions PDP and GDP through a 5 mm thick F82H membrane (FLiBe-blanket conditions) Hydrogen recycling coefficient: R = Γ reflection = γ Γ plasma Γ +Γ +Γ Membrane temperature Net Implantation total flux for PDP T 2 driving pressure for GDP reflection re emission T GDP Γ plasma (γ is the particle reflection coefficient) 550 o C 1 10 16 (2 10 16 ) DT cm-2 s-1 10 4 Pa (Fukada) (Thermodynamic equilibrium data) Γ plasma re-emission First wall Γ GDP T-GDP Γ PDP D,T-PDP Assuming a particle reflection coefficient of 0.5, the total incident flux is 2 10 16 D&T cm -2 s -1 and the first wall recycling rate has been estimated to be R= 1.018. 11

Motivation Potential issues associated with bi-directional DT permeation: PDP D lowers the recovery efficiency of T from the breeder. Gas-T permeation increases recycling on the first wall side. What is necessary to address the bi-directional hydrogen permeation and its associated reactor operation issues? Understandings of hydrogen PDP and GDP individually in detail; Evaluation of the tritium equilibrium pressure in the fuel recovery loop system; Determination of missing hydrogen isotopes transport parameters (solubility, diffusivity and surface recombination coefficient). No literature data available for the surface recombination coefficient (for F82H), which is important for recycling and retention. 12

Objectives I. To understand the mechanisms driving hydrogen isotopes transport processes. II. To demonstrate experimentally hydrogen transport phenomena that are predicted for the first wall of a fusion power reactor. III. To establish a database on hydrogen transport parameters for designing fusion power reactors. 13

2. Experimental facility and setup

Experimental facility VEHICLE-1 Plasma parameters: n e : ~10 10 cm -3 T e : 3~4 ev Ion flux : ~10 16 cm -2 s -1 Magnet Probe 1 Pneumatic shutter Test chamber 1kW ECR Plasma source Wave guide Plasma column Optical spectrometer IR pyrometer Probe 2 Turbo pump Y. Hirooka et al., J. Nucl. Mater. 337-339(2005)585-589. 15

Plasma characteristics in VEHICLE-1 Species mix modelling Langmuir Probe data ne Te 200W ECR power ne Te H 2 16

Experimental setup Permeation membrane sample 70 mm H-plasma Membrane φ35 mm H 2 gas Heater Bias Thickness: 0.5-5.0 mm Diameter: 35 mm Materials: F82H (Fe-8Cr-2W) SUS304 (Fe-19Cr-11Ni) The plasma- and gas-driven permeation fluxes are measured by two H 2 partial pressure gauges, respectively. Ion bombardment energy is provided by a negative bias (-100 V or -50 V). Temperature: ~200-520 o C 17

3. Bi-directional hydrogen permeation experiments and modelling 18

PDP Bi-directional permeation For the self-cooled breeder blankets, hydrogen isotopes will penetrate through the first wall by plasma-driven permeation (PDP) in one direction and gas-driven permeation (GDP) in the opposite direction. T + T + Edge plasma D + D + T + D + T + D + T + D + T + T + D + Implantation re-emission First wall ~5 mm ~500 o C T-GDP Diffusion D,T-PDP Diffusion Breeding blanket n + Li T + He T 2 T 2 T 2 T 2 GDP: Bi-directional permeation T 2 Dissociation T 2 & Solution T 2 T 2 Important parameters GDP: Solubility Diffusivity External pressure PDP: Surface recombination coefficient Diffusivity Implantation flux Reflection coef. 19

Hydrogen bi-directional permeation experiment H 2 partial pressure Membrane H-plasma Hα Spectrometer H 2 gas H 2 pressure Absolute pressure gauge H plasma: Te: ~10 ev Ne: ~1.0 10 10 cm -3 Bias: -50V H 2 pressure: Membrane: ~7x10 4 Pa 0.6 mm thick F82H H GDP flows in the counter direction to H PDP flow and affects the upstream plasma. The GDP flux has been measured to be 9.9x10 15 H/cm 2 /s. 20

Hydrogen bi-directional permeation modelling GDP-T 2 pressure: ~7x10 4 Pa Ion flux: ~8.5 10 15 H cm-2 s-1 Membrane: 0.6 mm thick α-fe Temperature: Gas-facing side: ~580 o C Plasma-facing side: and ~550 o C Boundary conditions: Gas side: Sieverts law Plasma side: recombination Intrinsic trap density: 1% Trapping energy: 0.62 ev Comparison of permeability 9.9x10 15 H/cm 2 /s The experimental result is in relatively good agreement with the prediction by modelling. 21

F82H: Hydrogen PDP and GDP fluxes Hydrogen PDP and GDP fluxes measured in VEHICLE-1 for a 5 mm thick F82H membrane (under the FLiBe-blanket conditions) Membrane temperature ~500 o C J + = D L J K 0 r Flow from blanket J= DT ( ) dc dx Flow from edge plasma Net Implantation flux for PDP H 2 driving pressure for GDP 2.0 10 16 H cm-2 s-1 (-100 V bias) ~10 4 Pa(Fukada) (Thermodynamic equilibrium data) Hydrogen transport in the first wall is dominated by the flow from the blanket. Using the experimental data, the hydrogen recycling rate has been estimated to be R= 1.025. 22

4. Modelling on the reactor fuel loop operation with hydrogen isotopes bi-directional permeation through the first wall 4-1. Re-evaluation of hydrogen bidirectional permeation fluxes through the first wall for FLiBe blankets 23

Re-evaluation of the tritium pressure for FLiBe blankets Re-evaluation of the tritium flows in a FLiBe loop has been performed, taking into account tritium leakage from the first wall.,,, M T : tritium inventory in the loop; J 1 : tritium production rate; J 2 : tritium PDP flow into blanket; J 3 : tritium GDP flow from blanket and GDP leak from pipes; J 4 : recovered tritium; J 5 : tritium GDP leak from heat exchanger Overall tritium inventory: dm T dt = J + J J J J 1 2 3 4 5 24

Re-evaluation of the tritium pressure for FLiBe blankets Conditions for analysis: *Song et al., Plasma Fusion Res. 7 (2012) 2405016. This work Fusion power 3 GW (Sagara) 1 GW Tritium breeding ratio 1.3 (Sagara) 1.25 *Song et al. Blanket surface area 3000 m 2 (Tanaka) 489 m 2 FLiBe flow rate 2.2 m 3 /s 2.2 m 3 /s Tritium recovery rate 0.99 (Sagara) 0.98 First wall Plasma flux Calculated tritium pressure 5 mm thick F82H 5 10 15 D cm 2 /s 5 10 15 T cm 2 /s 10 mm thick F82H with coatings No PDP assumed 1.1 10 3 Pa 4.3 10 3 Pa 25

4-2. Isotopes effects on hydrogen permeation 26

Modelling of bi-directional permeation involving multiple hydrogen isotopes Bi-directional GDP of the same hydrogen isotope Membrane T T T 2 gas T 2 gas p 1 p 2 Equal pressure: p 1 = p 2 GDP-T 2 pressure: ~10 4 Pa (Fukada) Membrane: 5 mm thick α-fe Temperature: 527 o C Boundary condition: Sieverts law Intrinsic trap density: 1% Trapping energy: 0.62 ev The tritium concentration profiles interact with each other in the two counter flows, finally reaching a flat profile with no net directional flow (quasi-thermodynamic equilibrium). 27

Modelling of bi-directional permeation involving multiple hydrogen isotopes Bi-directional PDP-D/T and GDP-T 2 flows. D+T Plasma T D First wall T T 2 gas T 2 pressure: ~1 Pa Implantation flux: 5 10 15 D cm -2 s-1 5 10 15 T cm-2 s-1 Membrane: 5 mm thick α-fe Temperature: 527 o C Boundary conditions: Gas side: Sieverts law Plasma side: recombination Intrinsic trap density: 1% Trapping energy: 0.62 ev The tritium concentration profiles interact with each other in the two counter flows. Deuterium flow appears to be independent of these tritium flows, driven by its own concentration gradient. 28

4-3. Re-evaluation of hydrogen bidirectional permeation fluxes through the first wall for FLiBe blankets 29

Re-evaluation of the bi-directional permeation process Conditions for DIFFUSE calculation: T 2 pressure: 1.1 10 3 Pa Implantation flux assumptions: Γ D = 5 10 15 D cm-2 s-1 Γ T = 5 10 15 T cm-2 s-1 Membrane: 5 mm thick α-fe Temperature: 550 o C T Boundary conditions: Gas side: Sieverts law D T Plasma side: recombination D/T release fluxes at the plasma-facing surface Intrinsic trap density: 1% Trapping energy: 0.62 ev The tritium release flux at the plasma-facing surface is a total flux from GDP and re-emission. 30

PDP and GDP under reactor-relevant conditions PDP and GDP through a 5 mm thick F82H membrane (FLiBe-blanket conditions) Hydrogen recycling coefficient: R Γ +Γ +Γ reflection re emission T GDP Γ plasma Γ plasma Γ reflection = γ Γ plasma (γ is the particle reflection coefficient) First wall Γ GDP T-GDP Membrane temperature 550 o C re-emission Γ PDP D,T-PDP Net Implantation flux for PDP 1 10 16 (2 10 16 ) DT cm-2 s-1 T 2 driving pressure for GDP 1.1x10 3 Pa (This work) Assuming a particle reflection coefficient of 0.5, the total incident flux is 2 10 16 D&T cm -2 s -1 and the first wall recycling rate has been estimated to be R= 1.006. 31

5. Experimental work on hydrogen isotopes transport parameter 5-1. Gas-driven permeation (GDP) experiments for the evaluation of solubility and diffusivity 32

Gas-driven permeation J Solubility of gases in metals (Sieverts' law): C= ST ( ) p H 2 S(T): Sieverts constant as a function of temperature PH 2 : hydrogen pressure Diffusion limited gas-driven permeation flux: Low pressure = DT ( ) dc dx DT ( ) ST ( ) = ( ph2 upsteam ph2 downsteam ) L L: membrane thickness D(T): diffusion coefficient PH 2 -downstream: hydrogen pressure at the low pressure side, usually PH 2 -downstream 0 ph downstream 2 H 2 Permeation H 2 H 2 H 2 x Membrane H H H H H Diffusion H H L C High pressure ph upstream 2 H 2 H 2 H 2 H 2 H 2 H 2 H 2 H 2 H 2 H 2 H 2 H 2 H 2 Dissociation & Solution H 2 H 2 H 2 Permeability for GDP: PT ( ) D( T) ST ( ) 33

Diffusion coefficient from GDP 3 methods to evaluate the diffusion coefficient, D: 1. Transient curve fitting method [1] Solving the diffusion equation: Cxt t 2 (,) Cxt (,) = DT ( ) 2 DT ( ) ST ( ) p 2 2 H n π J t = + DT t 2 L d 2 n ( ) [1 2 ( 1) exp( ( ) )] x n= 1 Transient GDP flux curve 2. Time lag method [2] : t = L 2 /6D l Jt ( l ) 0.63 Jt ( ) Time-integrated GDP flux 3. Breakthrough time method [2] : t b = L 1 1 ( ) D π L 2 2 2 6 15.3 D [1] E. Serra et al., J. Nucl. Mater. 245 (1997) 108-114. [2] MAV Devanathan et al., Proc. R. Soc. Lond. A270 (1962) 90-102. Time- 34

SUS304: Validation of the experimental setup and data analysis GDP through 0.65 mm thick SUS304 Diffusion coefficient Permeability PT ( ) D( T) ST ( ) = J GDP p H L 2 Diffusivity: 3 0.57 [ev] D = 5.5 10 exp [cm 2 s-1 ] kt Permeability: 10 0.62 [ ev] P = 9.6 10 exp( ) [mol cm -1 s-1 Pa-1/2 ] kt Solubility: 7 0.05 [ ev] S = 1.8 10 exp( ) [mol cm kt -3 Pa-1/2 ] Sun et al., Materials Science and Engineering, Al14 (1989) 179-187. B. L. Doyle and D. K. Brice, J. Nucl. Mater. 122 & 123 (1984) 1523-1530. 35

F82H: pressure and thickness dependence of GDP flux GDP (1 mm thick F82H) GDP at ~500 o C A linear relation between GDP flux and the square-root of upstream pressure has been found at all the temperatures examined in this work. GDP is diffusion-limited model: ( ) ( ) J = DT ( ) dc = D T S T p dx L GDP H 2 Even for a 5 mm thick membrane, reactor-relevant, hydrogen GDP has been found to be diffusion limited. 36

F82H: hydrogen diffusion coefficient evaluation The diffusion coefficient of H through F82H has been evaluated from the transient permeation behavior. A breaking point has been found for the diffusion coefficient data, which we attribute to the trapping effect. In the presence of traps: Measured diffusion coefficient: >~250 4 0.14 [ev] D = 7.5 10 exp kt < ~250 o C 0.50 [ev] D = 1.9 exp kt [cm 2 s -1 ] [cm 2 s -1 ] Cxt 2 (,) (,) t = D ( T) 2 t x t = D ( T) E. Serra et al., J. Nucl. Mater. 245 (1997) 108. V. Shestakov et al., J. Nucl. Mater. 307-311 (2002) 1494. T.V. Kulsartov et al., Fusion Eng. Des. 81 (2006) 701. a a 2 Cxt Cxt (,) 2 x C (,) xt DT ( ) D( T) (with traps) 37

F82H: hydrogen permeability and solubility The permeability and solubility of H in F82H have been evaluated from the steady state temperature dependent GDP data. 10 0.39 [ ev] P = 2.3 10 exp( ) kt [mol cm -1 s-1 Pa-1/2 ] E. Serra et al., J. Nucl. Mater. 245 (1997) 108. A. Pisarev et al., Phys. Scr. T94 (2001) 121. T.V. Kulsartov et al., Fusion Eng. Des. 81 (2006) 701. 7 0.25 [ev] S = 3.1 10 exp kt [mol cm -3 Pa-1/2 ] E. Serra et al., J. Nucl. Mater. 245 (1997) 108. Z. Yao et al., Fusion Sci. Tech. 48 (2005) 1285. K. Forcy et al., J. Nucl. Mater. 160 (1988) 117. 38

5-2. Plasma-driven permeation (PDP) experiments for the evaluation of surface recombination coefficients 39

Steady state plasma-driven permeation models Upstream Edge Plasma d Membrane Downstream J + Permeation Flux conservation: J = J + J 0 + J 0 Implantation c ( J - >> J + ) (DD) Diffusion c x Diffusion J + = d L J 0 (RD) Recombination (RR) J re-emission Recombination K r c d x x L Diffusion Recombination K 2 Recombinative desorption flux: J = Kr CH (C H is hydrogen atomic concentration at the upstream surface) L J + = D L r J K KL J = + K + K B.L. Doyle, J. Nucl. Mater. 111 & 112 (1982) 628-635. L 0 r J 0 40

Estimation of the net implantation flux for PDP experiments The ion species mix has been estimated by modelling. Particle reflection coefficient R: ~0.5 H 3 + is the dominant ion species in the electron temperature range of the experiments. The concentration of H + increases as the increase of electron temperature and becomes the dominant species when the electron temperature is higher than 4 ev. 41

F82H: membrane thickness effects on steady state PDP Steady state PDP data for F82H membranes at ~220 o C and ~500 o C (-100 V bias) Te ~3 ev ne ~1x10 10 cm -3 Net Implantation flux ~5x10 15 cm -2 s -1 Permeation flux ~10 13 cm -2 s -1 PDP flux ratio ~10-3 RD limited flux: J + = D L J K 0 r The steady state permeation flux is inversely proportional to the membrane thickness, meaning that the hydrogen transport process is in the recombination-diffusion limited regime. 42

F82H: surface recombination coefficient measurements The surface recombination coefficient of H on F82H has been measured from the steady state temperature dependent PDP data. The surface recombination coefficient: K r 21 [ ] 4.8 10 exp( 0.48 e V = ) kt [cm 4 s-1 ] The recombination coefficient for F82H has been estimated for the first time from the experimental data. M.I. Baskes, J. Nucl. Mater. 92 (1980) 318. S. Hirata et al., JAEA report, JAERI-M-91-081, pp. 276. 43

5-3. Surface condition effects on PDP 44

Earlier studies related to surface condition effects Surface condition effects interpretations in earlier studies: (1) changes in surface recombination coefficient due to sputtering or deposition of contaminations (Causey et al.); (2) changes in surface recombination coefficient and diffusivity due to ion-induced surface defects (Winter et al.). This work will: (1) propose a new hypothesis to interpret the surface contamination effects and; (2) propose a new model to describe the surface area/roughness effects. 45

Surface contamination effects on PDP (1) (1) Surface contamination effects (a) Clean surface Re-emission Concentration (RD-regime) Membrane Impurity film (b) Thin impurity film Membrane Permeation Permeation Implantation d L Range: d (~nm) (c) Thick impurity film Thickness: L A new hypothesis: Comparing the implantation range and the thickness of the impurity layer. Impurity layer d Membrane L Permeation 46

Surface contamination effects on PDP (2) The contamination effects have been investigated by the surface oxidization method. PDP with clean surface O 2 gas 1x10-2 Pa at 450-500 o C Surface oxidization PDP with oxidized surface Clean surface Chemical sputtering of surface oxides 2H+Fe x O y H 2 O+Fe x O y-1 H-plasma on H-plasma off Impurity layer d Membrane L After oxidization, it takes longer time for the permeation flux to reach steady state, suggesting a lower effective diffusion coefficient for the oxidized samples. 47

Surface contamination effects on PDP (2) Surface composition of the samples are examined by X-ray photoelectron spectroscopy (XPS). Polished surface After 45 min oxidization Impurity layer C O Fe Cr W C O Fe Cr W Membrane Depth (nm) d L After 45 min oxidization and PDP Hydrogen implantation profile C O Fe Cr W The thickness of the impurity layer is larger than the implantation range. 48

Surface area effect on PDP behavior (a) PDP through a flat surface Surface area A 0 (Net implantation) J 0 J 0- (Re-emission) (b) PDP through a modified surface Surface area A 1 (Net implantation) J0 = J0 + J0+ J1 = J1 + J1 θ J 1 + (Re-emission) J 1- J 0+ (Permeation) J 1+ (Permeation) 2 Surface concentration J0 J0 = KC r 0 ratio: 2 C / C J1 J1 = KC r 1 = J / J Surface flux ratio: J 1 /J 0 = A 0 /A 1 J = DT ( ) dc + dx The permeation flux ratio: 1 0 1 0 J 1 / J + 0+ A 1 / A 1 0 49

Surface area effects on PDP (1) To verify surface area effect, PDP experiments have been performed using samples with well controlled surface morphology. Modified surface-1 Three samples: Flat surface Modified surface-1 Modified surface-2 Modified surface-2 1 mm 0.4 mm Surface area ratio: A 1 /A 0 =6.4 Effective thickness: 4 mm< L eff < 5 mm Surface area ratio: A 1 /A 0 =3.2 Effective thickness: 4 mm< L eff < 5 mm 50

Surface area effects on PDP (2) The measured steady state permeation flux has been found to be inversely proportional to the square root of surface area, which is in good agreement with the model prediction. Flat surface A 1 /A 0 =1 Surface area ratio: A 1 /A 0 =6.4 Surface area ratio: A 1 /A 0 =3.2 1 51

Surface area effects on PDP (3) Polished surface Plasma-modified surface Compared with the polished surface, the steady state PDP flux for the plasma modified surface has been found to decrease by a factor of ~1.7, indicating a surface area ratio of A 1 /A 0 =~2.8. The surface area is increased due to the surface modification. 52

Combined two surface condition effects The modified surface 1 mm A combined effect has been observed as follows: Combined surface effect oxidization effect surface area effect ~3.2 ~ 1.4 ~2.5 53

Summary and future plans I. Bi-directional hydrogen permeation has been demonstrated for the first time in a laboratory-scale steady state plasma facility. II. III. IV. Revaluation of the dynamic tritium pressure in a fuel recovery loop has been performed. Hydrogen PDP and GDP through F82H have been investigated under some of the reactor-relevant conditions. Hydrogen transport parameters have been evaluated for F82H, including the recombination coefficient. V. Surface condition effects on hydrogen permeation have been examined and a new model has been proposed to interpret the data. VI. Future plans include: deuterium PDP and GDP experiments to evaluate the isotope effects and also surface coatings effects (i.e. W- coatings). 54