Presented at the COMSOL Conference 2010 Paris

Similar documents
Understanding of corrosion mechanisms after irradiation : effect of ion irradiation of the oxide layer on the corrosion rate of M5 alloy

Lecture 28: Kinetics of Oxidation of Metals: Part 1: rusting, corrosion, and

Defect Ch em Ch istry 1

Nuclear Data Uncertainty Analysis in Criticality Safety. Oliver Buss, Axel Hoefer, Jens-Christian Neuber AREVA NP GmbH, PEPA-G (Offenbach, Germany)

Basic Concepts of Electrochemistry

Loads on RPV Internals in a PWR due to Loss-of-Coolant Accident considering Fluid-Structure Interaction

EE 143 Microfabrication Technology Fall 2014

Tutorials : Corrosion Part 1: Theory and basics

Lecture 14. Electrolysis.

Massachusetts Institute of Technology Department of Materials Science and Engineering

Chapter 8: Internal Energy and the Laws of Thermodynamics

State-Space Modeling of Electrochemical Processes. Michel Prestat

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

Lecture 29: Kinetics of Oxidation of Metals: Part 2: Wagner Parabolic

5. Defects Thermal defects Planar defects Linear defects Point defects stochiometric compounds Schottky defects Anti-Schottky defects

Diffusion. Diffusion = the spontaneous intermingling of the particles of two or more substances as a result of random thermal motion

Advanced Analytical Chemistry Lecture 12. Chem 4631

Nucleophilic attack on ligand

Electrochemistry. Goal: Understand basic electrochemical reactions. Half Cell Reactions Nernst Equation Pourbaix Diagrams.

surface c, c. Concentrations in bulk s b s b red red ox red

Multidimensional, Non-Isothermal, Dynamic Modelling Of Planar Solid Oxide Fuel Cells

Diffusion in multicomponent solids. Anton Van der Ven Department of Materials Science and Engineering University of Michigan Ann Arbor, MI

Section 7: Diffusion. Jaeger Chapter 4. EE143 Ali Javey

Understanding the enhancement effect of high-temperature steam on the carbonation reaction of CaO with CO 2

Initial position, x p (0)/L

INTRODUCTION TO THE DEFECT STATE IN MATERIALS

FUEL CELLS in energy technology (4)

MATERIALS THERMODYNAMICS KINETICS SUMMARY INTRODUCTION

This unit explores the problem solving process. We will see that there are organized ways to approach and solve problems. Applied mathematics forms

Extrinsic Defect Reactions in

Modeling ion and electron profiles in methane-oxygen counterflow diffusion flames

Separation through Dialysis

The mathematical description of the motion of Atoms, Molecules & Other Particles. University of Rome La Sapienza - SAER - Mauro Valorani (2007)

Thermodynamics of CeO 2 thermochemical fuel production - Supporting Information

Previous class. Today. Porous Electrodes. Impedance of a film on electrode surface. Cylindrical pore

Novel Ceramic Membrane for High Temperature Carbon Dioxide Separation

A First Course on Kinetics and Reaction Engineering Example 1.2

Design constraints Maximum clad temperature, linear power rating

Thermodynamic and Stochiometric Principles in Materials Balance

Characterisation of uranium vacancies in hyper stoichiometric uranium dioxide

Effect of Fission Fragments on the Properties of UO 2 Fuel of Pressurized Water Reactors

Fluid Dynamics and Balance Equations for Reacting Flows

Accurate measurement of elemental impurities in metals and metal alloys using the Thermo Scientific icap TQ ICP-MS

ELECTROCHEMICAL SYSTEMS

Exergy. What s it all about? Thermodynamics and Exergy

Improved PWR Simulations by Monte-Carlo Uncertainty Analysis and Bayesian Inference

On the use of SERPENT code for few-group XS generation for Sodium Fast Reactors

Electrochemical Cell - Basics

Finite Element Modules for Enhancing Undergraduate Transport Courses: Application to Fuel Cell Fundamentals

WDHS Curriculum Map: Created by Erin Pence September 2010

EEC 503 Spring 2009 REVIEW 1

Reactions and Reactors

Oxide growth model. Known as the Deal-Grove or linear-parabolic model

Atomic Transport & Phase Transformations. PD Dr. Nikolay Zotov

Electrons, Holes, and Defect ionization

Lecture notes for /12.586, Modeling Environmental Complexity. D. H. Rothman, MIT October 22, 2014

Quiz, Physics & Chemistry

Definitions and Concepts

Periodic Classification and Properties Page of 6

PROTON CERAMIC ELECTRODICS

Modeling as a tool for understanding the MEA. Henrik Ekström Utö Summer School, June 22 nd 2010

Carbon dioxide removal processes by alkanolamines in aqueous organic solvents Hamborg, Espen Steinseth

Introduction into defect studies. in ceramic materials(iii) Structure, Defects and Defect Chemistry. Z. Wang. January 18, 2002

The Thermodynamic Characterization of Spontaneous Electrochemical Reactions

3.022 Microstructural Evolution in Materials

MSE 360 Exam 1 Spring Points Total

Topics to discuss...

TABLE OF CONTENT. Chapter 4 Multiple Reaction Systems 61 Parallel Reactions 61 Quantitative Treatment of Product Distribution 63 Series Reactions 65

Kinetics. Rate of change in response to thermodynamic forces

SEMICONDUCTOR PHYSICS REVIEW BONDS,

1. Why are chemical reactions important to energy, environmental and process engineering? Name as many reasons as you can think of.

Modeling of Liquid Water Distribution at Cathode Gas Flow Channels in Proton Exchange Membrane Fuel Cell - PEMFC

Table of Contents Preface PART I. MATHEMATICS

Modeling Hydrogen Permeation through a TiO 2 Film and Palladium

9.4 Effusion and Diffusion of Gases

Electrical Double Layer Energy Storage Capacitors Up to 3 V Operating Voltage

EMA4303/5305 Electrochemical Engineering Lecture 03 Electrochemical Kinetics

The moderator temperature coefficient MTC is defined as the change in reactivity per degree change in moderator temperature.

The Mechanism of Electropolishing of Nb in Hydrofluoric-Sulfuric Acid (HF+H 2 SO 4 ) Electrolyte

DATA SHEET. PDTA114E series PNP resistor-equipped transistors; R1 = 10 kω, R2 = 10 kω DISCRETE SEMICONDUCTORS

8 Phenomenological treatment of electron-transfer reactions

Ruggedized Electrical Double Layer Energy Storage Capacitors

Similarities and differences:

KINETICS OF OXIDE GROWTH ON METAL SURFACES

Appendix 4. Appendix 4A Heat Capacity of Ideal Gases

Nicholas Cox, Pawel Drapala, and Bruce F. Finlayson Department of Chemical Engineering, University of Washington, Seattle, WA, USA.

Pre-AP Biology Energy Unit Study Guide Part I

In this chapter, we shall derive the DC and AC conductivities for the Bloch electrons by using a relaxation time approximation.

Carriers Concentration and Current in Semiconductors

Electrolytic processes Notes

Phase-Field Modeling of Technical Alloy Systems Problems and Solutions

Cerium dioxide redox cycle for fuel production

THORIUM SELF-SUFFICIENT FUEL CYCLE OF CANDU POWER REACTOR

Reactivity of the Aluminium Surface in Aqueous Solutions

8.1 Drift diffusion model

Example Resistive Heating

Fundamentals of Nuclear Reactor Physics

3.155J/6.152J Microelectronic Processing Technology Fall Term, 2004

Appendix 1: List of symbols

single single single yellow S HH H Maximum Ratings at T A = 25 C, unless otherwise specified Parameter Symbol Value Unit Diode reverse voltage V R

Transcription:

Presented at the CMSL Conference 2010 Paris

Solid State Transport of Reactive Charged Particles Application to Metal xidation P. BUTTIN, B. MALKI, P. BARBERIS, B. BARUX

Introduction Zirconium corrosion in power plant coolant Fuel assembly life Zr Nuclear Power Plant Picture AREVA NP Fuel assembly Picture AREVA NP Understanding corrosion key factor impact on oxide growth Entreprise - AREVA NP - All rights reserved - p.3

bjectives Develop a numerical model for zirconium oxidation Simulation using CMSL Multiphysics Multicomponent transport coupled with internal chemical reaction Zr 2 Zr V Coolant K + = V e V V µ e Anionic conductivity transport V, V, e Electrons trapping V - + e = V - Constraints on system : Electroneutrality + Stoichiometry + Chemical reaction Entreprise - AREVA NP - All rights reserved - p.4

Conservative Ensembles Considering ensemble transport (J. Maier, 1994) Conservative ensembles xygen ensemble (Vacancies, oxygen interstitial) Cross effect Electron ensemble (Electrons, hole, ) V V e Thermodynamic of irreversible processes J J e = s = s e η η s s e ee η η e e (Lankhorst, 1996) Ensembles are conservatives for the chemical reaction Chemical reaction : Cross effect between ensembles fluxes Entreprise - AREVA NP - All rights reserved - p.5

xygen Flux xygen flux given by: (Wagner, 1977) σ J = µ = D c 4F amb δ 2 Chemical diffusion coefficient σ amb dµ dc δ σ amb dµ D = 2 4F dc : Ambipolar conductivity : Chemical capacity Constraints on system c + 2 c c = 0 Electroneutrality V K e V c + c + c = 0 Stoichiometry V V = c V c V c e Chemical Reaction Transport equation: c t = D δ c All system s information are given by D δ Entreprise - AREVA NP - All rights reserved - p.6

xide Growth xide growth given by: de dt ox 1 = J Boundary condition on internal interface co,int Keeping subdomain constant: x = xreal 1 ; e = e reduced real reduced ox ox Simulation in one dimension: xide growth model formulated as boundary weak form Metal xide Coolant e ox de ox /dt, δ int Zr 2-δ δ ext Entreprise - AREVA NP - All rights reserved - p.7

xide Growth Model parameters Symbol K D V Description Chemical equilibrium between vacancies Diffusion coefficient of vacancies V D V Diffusion coefficient of vacancies V D e δ int δ ext T Diffusion coefficient of electrons Zirconia non stoichiometry at the M/ interface Zirconia non stoichiometry at the surface Temperature Entreprise - AREVA NP - All rights reserved - p.8

Simulations Comparison with experimental data xide thickness (µm) 5 4 3 2 1 Experimental data* xidation model 0 0 100 200 300 Time (days) * J. Vermoyal, A. Hammou, L. Dessemond, A. Frichet, Elec. Acta, 2002 Simulation gives good results according to experimental data Entreprise - AREVA NP - All rights reserved - p.9

Effect on reaction and transport coefficient Simulations xide thickness (µm) 3 2 1 D V D V D V = D V = 0,1D V = 10D V V ionized trapping 0 0 100 200 300 Time (days) Effect of electron trapping on oxide growth kinetic Entreprise - AREVA NP - All rights reserved - p.10

Simulations Effect on zirconia non stoichiometry on inner interface Zr 2-δ xide thickness (µm) 1,8 δ 1 1,2 0,6 δ 2 δ 3 δ 1 > δ 2 > δ 3 0,0 0 100 200 300 Time (days) Effect oxygen chemical potential gradient through the oxide Entreprise - AREVA NP - All rights reserved - p.11

Conclusions Development of an oxide growth model for anionic transport Model is not specific for zirconium oxidation Assessment of the control parameters on growth kinetics Transport coefficient, stoichiometry, Environmental effect could be add on the model 2D model version when symmetry is broken Entreprise - AREVA NP - All rights reserved - p.12

Thank you for your attention! Entreprise - AREVA NP - All rights reserved - p.13

Any reproduction, alteration, transmission to any third party or publication in whole or in part of this document and/or its content is prohibited unless AREVA NP has provided its prior and written consent. This document and any information it contains shall not be used for any other purpose than the one for which they were provided. Legal action may be taken against any infringer and/or any person breaching the aforementioned obligations. Entreprise - AREVA NP - All rights reserved - p.14