IV. Transport Phenomena Lecture 35: Porous Electrodes (I. Supercapacitors)

Similar documents
Electrochemical Equilibrium Electromotive Force. Relation between chemical and electric driving forces

The Hyperelastic material is examined in this section.

Grand Canonical Ensemble

Chapter 2 Theoretical Framework of the Electrochemical Model

Physics of Very High Frequency (VHF) Capacitively Coupled Plasma Discharges

A Note on Estimability in Linear Models

CHAPTER 33: PARTICLE PHYSICS

From Structural Analysis to FEM. Dhiman Basu

:2;$-$(01*%<*=,-./-*=0;"%/;"-*

Jones vector & matrices

Review - Probabilistic Classification

ACOUSTIC WAVE EQUATION. Contents INTRODUCTION BULK MODULUS AND LAMÉ S PARAMETERS

Lecture 14. Relic neutrinos Temperature at neutrino decoupling and today Effective degeneracy factor Neutrino mass limits Saha equation

Economics 600: August, 2007 Dynamic Part: Problem Set 5. Problems on Differential Equations and Continuous Time Optimization

External Equivalent. EE 521 Analysis of Power Systems. Chen-Ching Liu, Boeing Distinguished Professor Washington State University

September 27, Introduction to Ordinary Differential Equations. ME 501A Seminar in Engineering Analysis Page 1. Outline

Phys 774: Nonlinear Spectroscopy: SHG and Raman Scattering

Lecture 3: Phasor notation, Transfer Functions. Context

The root mean square value of the distribution is the square root of the mean square value of x: 1/2. Xrms

ECE507 - Plasma Physics and Applications

Ερωτήσεις και ασκησεις Κεφ. 10 (για μόρια) ΠΑΡΑΔΟΣΗ 29/11/2016. (d)

Lecture 23 APPLICATIONS OF FINITE ELEMENT METHOD TO SCALAR TRANSPORT PROBLEMS

1- Summary of Kinetic Theory of Gases

MECH321 Dynamics of Engineering System Week 4 (Chapter 6)

te Finance (4th Edition), July 2017.

Physical and Interfacial Electrochemistry Lecture 4. Electrochemical Thermodynamics. Thermodynamics of electrochemical systems

1) They represent a continuum of energies (there is no energy quantization). where all values of p are allowed so there is a continuum of energies.

Heisenberg Model. Sayed Mohammad Mahdi Sadrnezhaad. Supervisor: Prof. Abdollah Langari

EDGE PEDESTAL STRUCTURE AND TRANSPORT INTERPRETATION (In the absence of or in between ELMs)

CHAPTER 4. The First Law of Thermodynamics for Control Volumes

Soft k-means Clustering. Comp 135 Machine Learning Computer Science Tufts University. Mixture Models. Mixture of Normals in 1D

A NEW GENERALISATION OF SAM-SOLAI S MULTIVARIATE ADDITIVE GAMMA DISTRIBUTION*

signal amplification; design of digital logic; memory circuits

SPECTRUM ESTIMATION (2)

Discrete Shells Simulation

167 T componnt oftforc on atom B can b drvd as: F B =, E =,K (, ) (.2) wr w av usd 2 = ( ) =2 (.3) T scond drvatv: 2 E = K (, ) = K (1, ) + 3 (.4).2.2

Fakultät III Wirtschaftswissenschaften Univ.-Prof. Dr. Jan Franke-Viebach

The Fourier Transform

Fakultät III Univ.-Prof. Dr. Jan Franke-Viebach

Analyzing Frequencies

Phy213: General Physics III 4/10/2008 Chapter 22 Worksheet 1. d = 0.1 m

10/7/14. Mixture Models. Comp 135 Introduction to Machine Learning and Data Mining. Maximum likelihood estimation. Mixture of Normals in 1D

An Improved Mathematical Model of Lead-Acid Batteries for Simulation of VRLA Batteries

Relate p and T at equilibrium between two phases. An open system where a new phase may form or a new component can be added

IV. First Law of Thermodynamics. Cooler. IV. First Law of Thermodynamics

Two-Temperature Plasma Modeling of Argon Gas Tungsten Arcs.

8-node quadrilateral element. Numerical integration

Numerical Study of Two-fluid Flowing Equilibria of Helicity-driven Spherical Torus Plasmas

Basic Electrical Engineering for Welding [ ] --- Introduction ---

A Quick introduction to Quantum Monte Carlo methods

HORIZONTAL IMPEDANCE FUNCTION OF SINGLE PILE IN SOIL LAYER WITH VARIABLE PROPERTIES

Electrochemistry L E O

Variational Approach in FEM Part II

Lecture #15. Bipolar Junction Transistors (BJTs)

COMPLEX NUMBER PAIRWISE COMPARISON AND COMPLEX NUMBER AHP

CHAPTER 7d. DIFFERENTIATION AND INTEGRATION

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

Folding of Regular CW-Complexes

UNIT 8 TWO-WAY ANOVA WITH m OBSERVATIONS PER CELL

JEE-2017 : Advanced Paper 2 Answers and Explanations

Derivation of Electron-Electron Interaction Terms in the Multi-Electron Hamiltonian

Heating of a solid cylinder immersed in an insulated bath. Thermal diffusivity and heat capacity experimental evaluation.

Background: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals.

genius PHYSICS Energy Bands. Types of Solids.

Group Codes Define Over Dihedral Groups of Small Order

Fundamentals of Continuum Mechanics. Seoul National University Graphics & Media Lab

Chapter 2: The Photosphere

Polytropic Process. A polytropic process is a quasiequilibrium process described by

Quasi-Classical States of the Simple Harmonic Oscillator

Lucas Test is based on Euler s theorem which states that if n is any integer and a is coprime to n, then a φ(n) 1modn.

LEBANESE UNIVERSITY FACULTY OF ENGINEERING

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.

Einstein Equations for Tetrad Fields

[ ] 1+ lim G( s) 1+ s + s G s s G s Kacc SYSTEM PERFORMANCE. Since. Lecture 10: Steady-state Errors. Steady-state Errors. Then

A Probabilistic Characterization of Simulation Model Uncertainties

Where k is either given or determined from the data and c is an arbitrary constant.

Modelling of new generation plasma optical devices

Differentiation of Exponential Functions

Lesson 7. Chapter 8. Frequency estimation. Bengt Mandersson LTH. October Nonparametric methods: lesson 6. Parametric methods:

Physics 256: Lecture 2. Physics

Elements of Statistical Thermodynamics

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

??? Dynamic Causal Modelling for M/EEG. Electroencephalography (EEG) Dynamic Causal Modelling. M/EEG analysis at sensor level. time.

Voltage, Current, Power, Series Resistance, Parallel Resistance, and Diodes

EEC 686/785 Modeling & Performance Evaluation of Computer Systems. Lecture 12

Lecture 21. Boltzmann Statistics (Ch. 6)

Search sequence databases 3 10/25/2016

Outlier-tolerant parameter estimation

The pn junction: 2 Current vs Voltage (IV) characteristics

Electrochemical Energy Systems Spring 2014 MIT, M. Z. Bazant. Midterm Exam

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory

Optimal Ordering Policy in a Two-Level Supply Chain with Budget Constraint

P A = (P P + P )A = P (I P T (P P ))A = P (A P T (P P )A) Hence if we let E = P T (P P A), We have that

Chapter 6. pn-junction diode: I-V characteristics

ARC-ELECTRODE INTERACTIONS IN SILICON AND FERROSILICON FURNACES

Linear Algebra Provides a Basis for Elasticity without Stress or Strain

ST 524 NCSU - Fall 2008 One way Analysis of variance Variances not homogeneous

Radial Cataphoresis in Hg-Ar Fluorescent Lamp Discharges at High Power Density

Properties of ferromagnetic materials, magnetic circuits principle of calculation

Vibrational analysis for the nuclear electronic orbital method

Transcription:

IV. Transort Phnomna Lctur 35: Porous Elctrods (I. Surcaactors) MIT Studnt (and MZB) 1. Effctv Equatons for Thn Doubl Layrs For surcaactor lctrods, convcton s usually nglgbl, and w dro out convcton trms hr. Lt s focus on ffctv quatons govrnng th transorts and lctrostatcs n lctrolyt. Fgur 1. Flat Elctrod Surfac Scs consrvaton quatons, Nrnst-Plank flux consttutv quatons, and Posson quaton mak u Posson-Nrnst-Plank (PNP) st of quatons (bold fonts ndcat that th varabls ar n vctor quantty): c t F 0 (1) F Dc 2 zc (2) (3) Elctrostatc constran and flux constrans on th surfac scfy boundary condtons. Elctrostatc constran can b ntrrtd dffrntly, gvn dffrnt scfd varabls. Whn thr s a scfd amount of surfac charg, w can hav Gaussan law satsfyng th lctrostatc constran:

Lctur 35 10.626 Elctrochmcal Enrgy Systms (2011) Bazant nˆ ( ) qs (4) On th othr hand, whn gvn a scfd surfac otntal, w can us th followng aroxmaton, Equaton (5), nstad of Equaton (4), satsfyng th lctrostatc constran. Ths boundary condton could b mor smlfd assumng nglgbl caacty n Strn layr. (x = 0) nˆ (x = 0) (x = 0) s (5) W can us thr Equaton (4) or (5) to satsfy th lctrostatc constran, dndng on scfd varabls at th surfac. In addton, flux constrans of scs scfy th rst of ncssary boundary condtons: nf ˆ R (6) W now aly th abov st of quatons (PNP) as wll as boundary condtons to th orous lctrod wth doubl layr thcknss far thnnr than th or lngth scal. Fgur 2. Thn Doubl Layr n a Por Whn th or lngth scal s far largr than th lngth scal of doubl layr,, w hav saraton of lngth scals. Th mathmatcal structur of thn doubl layr roblm s wll undrstood from th rsctv of sngular rturbaton analyss, n whch ach of two rgons rqurs a dffrnt aroxmaton. Two dffrnt aroxmatons ar constrand by matchd asymtotc xansons. W frst assgn notatons for dffrnt varabls. From now on throughout ths lctur, w us th followng notatons. : Concntraton of Scs n Doubl Layr : Concntraton of Scs n Bulk Elctrolyt

Lctur 35 10.626 Elctrochmcal Enrgy Systms (2011) Bazant : Chmcal Potntal of Scs n Doubl Layr : Chmcal Potntal of Scs n Bulk Elctrolyt Fgur 3. Varabls n Two Dffrnt Rgons In quas-nutral bulk lctrolyt, w can us th quas-nutral aroxmaton, and conclud wth zro dvrgnc of currnt dnsty. z c 0 (7) j z F 0 (8) As statd abov dscusson, asymtotc xansons of aroxmatd varabls n two dffrnt rgons should b matchd wth th corrsondng ar rsctvly. In a or, varabls ar wll aroxmatd by th quas-qulbrum aroxmaton, and chmcal otntals n th two rgons ar aroxmatly constant across th or. lmcˆ lmc xˆ ˆ x0 (9) (10) Now, condtons constranng th varabls n doubl layr can b found by dfnng surfac varabls as shown n Fgur 3 (Strn layr caacty gnord): cˆ lm ˆ 0 c dx (11) x F D S S (12) Thn th ffctv boundary condton usng th varabls abov s (Chu and Bazant, 2007): S S ˆ S F R n F (13) t

Lctur 35 10.626 Elctrochmcal Enrgy Systms (2011) Bazant whr : Excss Surfac Concntraton of Scs r Ara : Surfac Flux of Scs : Dmnsonlss Chmcal Potntal of Scs : Nt racton rat of Scs Usng Bulk Varabls ( Total dffus charg dnsty can b calculatd wth xcss surfac concntratons. q z (14) 2. Porous Elctrods Formal drvaton by volum avragng (homognzaton) gos from th thn doubl layr quatons nsd ors to macroscoc artal dffrntal quatons. An arorat modl would nvolv volum avrags of varabls ovr a volum lmnt small comard to th ovrall dmnsons (L), but larg than th or structur lngth scal (h ). Hnc, to hav ths modl vald, th followng condton s ncssary. L h D (15) Fgur 4. Lngth Scals n Porous Elctrod Not always, but for surcaactors, w do not dlt th salt ons n lctrolyt. Thus t s a vald aroxmaton to assum that concntratons ( and otntals ( and ) ar varyng slowly n th macroscoc vwont. Ths assumton justfs th volum-avragng or homognzaton. In ths macroscoc tratmnt, w do not consdr th actual gomtrc dtals of th ors. Rathr, w dfn macroscoc otntal n lctrolyt, otntal n sold matral, and on concntratons to b contnuous and wll-dfnd functons of sac coordnats. As a rsult, th orous lctrod n ths modl s rrsntd by th suroston of two contnuous mda

Lctur 35 10.626 Elctrochmcal Enrgy Systms (2011) Bazant wthout mcrostructur, on corrsondng to th lctrolyt soluton and th othr corrsondng to th sold matral matrx. In ths modl both mda ar dfnd n th whol doman. Thrfor, n ths macroscoc modl, otntal n lctrolyt ( ) as wll as otntal n conductng sold matral ( ) ar dfnd n th whol doman of sac, whras thy wr only dfnd n ach has formrly. Volum-avragd concntratons n macroscoc vwont ar dffrnt from and rlatd to th concntratons n bulk lctrolyt and n doubl layrs as shown n th followng quaton (n formr lcturs, was usd for th concntratons n rsrvors or n nlt flows of ful. In ths lctur, s usd for volum-avragd concntratons n macroscoc vwont). Th macroscoc on concntratons ar dfnd throughout th whol volum as wll. c c a (16) Th frst trm on rght hand sd corrsonds to th contrbutons from bulk lctrolyt, and th rd scond trm s from xcss concntratons n doubl layrs. Nwman s book (3 Ed, 2004) dos not tak account of th scond trm. Howvr, ths trm may hav sgnfcant ffcts scally whn w consdr surcaactors whos chargs ar mostly stord n th doubl layrs. Wth th varabls dfnd wth th macroscoc modl, w can construct th scs consrvaton quatons, as wll as th flux consttutv quatons: c x F R (17) F Dc (18) W could bttr modl th macroscoc flux quaton, ncludng th transort n doubl layrs as shown n Equaton (19) and Fgur 5. Transort n doubl layr (surfac transort) would b sgnfcant whn th stord dffus charg s larg. F D c D a S (19) Fgur 5. Transort n Doubl Layr (Surfac Transort)

Lctur 35 10.626 Elctrochmcal Enrgy Systms (2011) Bazant For th macroscoc volum-avragd racton trm ( ), w could modl t by usng th surfac racton trm basd on th bulk concntratons and otntals: S R a R ( c,, ) (20) whr : Dmnsonlss Chmcal Potntal of Scs : Effctv Dffusvty n orous mdum (s Lctur 34) ϵ a : Effctv Surfac Dffusvty : Macroscoc Porosty : Surfac Ara Dnsty Wth th macroscoc concntratons and fluxs abov, w can also dfn macroscoc currnt dnsty as wll as charg dnsty. For th charg dnsty, only xcss surfac concntratons mattr, snc bulk lctrolyt ks quas-nutralty. j zf z c a z a q D (21) (22) Du to macroscoc nutralty,. Ths macroscoc charg dnsty s consrvd by th followng quaton. Th ngatv sgn n Equaton (23) s du to th dfnton of currnt dnsty, whch dfns anodc currnts to b ostv. Now, w dfn lctronc currnt dnsty n a orous lctrod. j j F zr (23) t j (24) Du to macroscoc charg consrvaton, th dvrgnc of total currnt s zro. j j 0 (25) Wth th govrnng quatons dscussd abov, w can now undrstand what th solutons look lk for dffrnt alcatons. In ths lctur, w focus on gttng and undrstandng th soluton for surcaactors. Th solutons for battrs and ful clls wll b dscussd n th nxt lctur (Lctur 36). In surcaactor alcatons, w do not hav sgnfcant faradac ractons on th lctrod or surfacs. Thus, w can gt th soluton, sttng th faradac

Lctur 35 10.626 Elctrochmcal Enrgy Systms (2011) Bazant transfrnc currnt to b zro ( ). On th othr hand, for battrs and ful clls, w can nglct caactanc of doubl layrs, and focus on faradac contrbutons. Also, t s rqurd to modl racton roducts, n addton to th ractng ons and lctrons. 3. Surcaactors In a modl for surcaactors, w can nglct faradac transfrnc currnt ( ) and changs n salt concntraton du to caacty chargng of doubl layrs (c.f. salt dlton du to doubl layr chargng s consdrabl n caacty dsalnaton alcatons). Thn, w can dfn ffctv onc conductvty (, aroxmatd to b constant. j (26) Th charg consrvaton quaton, Equaton (23), thn bcoms: s j (27) t t Combnng Equaton (24) and (25), j j 2 2 (28) W nd to mloy a doubl layr caactanc modl (C D ), such as DH modl or GC modl, to rlat th charg dnsts to th otntals. q t D C (,, c,...) t D (29) Pluggng n Equaton (22) and (28) to th abov Equaton (29), w now hav an quaton govrnng th otntals n lctrolyt and sold matral. ac 2 2 D t (30) Equaton (30) contans two govrnng artal dffrntal quatons for otntal n lctrolyt ( ) and otntal for sold matral ( ). To solv th govrnng quatons, w nd boundary condtons and ntal condtons. Th boundary condtons ar shown n Fgur 6. On th surfac contactng wth a sarator, lctron flux s not allowd and th normal comonnt of otntal gradnt n sold matral s zro. In smlar mannr, on th surfac contactng wth a currnt collctor, on flux s not allowd and th normal comonnt of otntal gradnt n lctrolyt s zro. Th otntal n lctrolyt at th surfac contactng a sarator s st to b th zro rfrnc valu. And w st otntal (V), may b a functon of tm, n conductng sold matral at th surfac contactng a currnt collctor.

Lctur 35 10.626 Elctrochmcal Enrgy Systms (2011) Bazant Fgur 6. Boundary Condtons for Macroscoc Potntals For ntal condtons, lt s start from charg-otntal rlatonsh. Th charg dnsty (q) s rlatd to th caactanc and th otntal dffrnc wth th followng quaton: q( x, t) C (31) W consdr startng from th qulbrum voltag (on crcut voltag, V o ). Ths lads us to hav a unform ntal charg dnsty throughout th orous lctrod. In gnral, ths ntal condton can b xrssd usng Equaton (32) and (33). q( x,0) unform CVo ( x,0) ( x,0) V o (32) (33) Th govrnng quatons, Equaton (30), could b solvd for surcaactor alcatons wth th boundary condtons shown n Fgur 6, and th ntal condton, Equaton (33). Th scfc solutons wll b dscussd n th nxt lctur, Lctur 36. 4. Transmsson Ln Modl In fact, th govrnng quatons, Equaton (30), could b quvalntly obtand by usng th transmsson ln modl. Ths can b asly sn by manulatng xrssons of conductvty and caactanc, usng rsstancs and caactanc r lngth.

Lctur 35 10.626 Elctrochmcal Enrgy Systms (2011) Bazant Fgur 7. Macroscoc Cross-sctonal Ara R R 1 A C Aa C 1 A D : Rsstanc to Elctron Currnt r Lngth : Rsstanc to Ion Currnt r Lngth : Total Caactanc of Doubl Layr r Lngth Wth th nwly dfnd aramtrs abov, w can rarrang th govrnng quaton, Equaton (30), and show that ths quaton s quvalnt to th transmsson ln modl. By dong so, w can now undrstand whch hyscal bass th crcut lmnts ar basd on. 1 1 C R R t 2 2 (34) Fgur 8. Transmsson Ln Equvalnt Crcut Modl In nxt lctur, Lctur 36, w wll solv th govrnng quatons, Equaton (30) or quvalntly Equaton (34), usng th boundary condtons shown n Fgur 6, and th ntal condton, Equaton (33), for surcaactor alcatons.

MIT OnCoursWar htt://ocw.mt.du 10.626 / 10.462 Elctrochmcal Enrgy Systms Srng 2011 For nformaton about ctng ths matrals or our Trms of Us, vst: htt://ocw.mt.du/trms.