Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Similar documents
Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

Basic Circuit Elements Professor J R Lucas November 2001

Vehicle Arrival Models : Headway

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

IB Physics Kinematics Worksheet

Some Basic Information about M-S-D Systems

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KINEMATICS IN ONE DIMENSION

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

Math 333 Problem Set #2 Solution 14 February 2003

Sub Module 2.6. Measurement of transient temperature

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits

Final Spring 2007

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

4.6 One Dimensional Kinematics and Integration

EE650R: Reliability Physics of Nanoelectronic Devices Lecture 9:

UCLA: Math 3B Problem set 3 (solutions) Fall, 2018

Note: For all questions, answer (E) NOTA means none of the above answers is correct.

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

! ln 2xdx = (x ln 2x - x) 3 1 = (3 ln 6-3) - (ln 2-1)

AP Chemistry--Chapter 12: Chemical Kinetics

Solutions to Assignment 1

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du.

Chapter 2. First Order Scalar Equations

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

t 2 B F x,t n dsdt t u x,t dxdt

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

Non-uniform circular motion *

CHAPTER 12 DIRECT CURRENT CIRCUITS

Linear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation

2. Nonlinear Conservation Law Equations

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

Unsteady Mass- Transfer Models

The Quantum Theory of Atoms and Molecules: The Schrodinger equation. Hilary Term 2008 Dr Grant Ritchie

= ( ) ) or a system of differential equations with continuous parametrization (T = R

Chapter 7 Response of First-order RL and RC Circuits

Motion along a Straight Line

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

Module 4: Time Response of discrete time systems Lecture Note 2

6.003 Homework 1. Problems. Due at the beginning of recitation on Wednesday, February 10, 2010.

Matlab and Python programming: how to get started

Math 2214 Solution Test 1B Fall 2017

ES.1803 Topic 22 Notes Jeremy Orloff

Sterilization D Values

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation

4.1 - Logarithms and Their Properties

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit

5.1 - Logarithms and Their Properties

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC

Section 4.4 Logarithmic Properties

Hall effect. Formulae :- 1) Hall coefficient RH = cm / Coulumb. 2) Magnetic induction BY 2

AP Calculus BC Chapter 10 Part 1 AP Exam Problems

15. Vector Valued Functions

10. State Space Methods

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

Predator - Prey Model Trajectories and the nonlinear conservation law

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Linear Response Theory: The connection between QFT and experiments

20. Applications of the Genetic-Drift Model

4.5 Constant Acceleration

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

) were both constant and we brought them from under the integral.

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18

Chapter 14 Chemical Kinetics

INDEX. Transient analysis 1 Initial Conditions 1

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

6.2 Transforms of Derivatives and Integrals.

ENV 6015 Solution to Mixing Problem Set

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

ψ(t) = V x (0)V x (t)

6. Solve by applying the quadratic formula.

CHEAPEST PMT ONLINE TEST SERIES AIIMS/NEET TOPPER PREPARE QUESTIONS

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

Chapter 14 Chemical Kinetics

Section 7.4 Modeling Changing Amplitude and Midline

Structural Dynamics and Earthquake Engineering

Math 10B: Mock Mid II. April 13, 2016

EXERCISES FOR SECTION 1.5

CHEMICAL KINETICS: 1. Rate Order Rate law Rate constant Half-life Temperature Dependence

PHYS 1401 General Physics I Test 3 Review Questions

Solution: b All the terms must have the dimension of acceleration. We see that, indeed, each term has the units of acceleration

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

HOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.

MOMENTUM CONSERVATION LAW

A. Using Newton s second law in one dimension, F net. , write down the differential equation that governs the motion of the block.

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Heat Transfer. Revision Examples

Transcription:

Module Fick s laws of diffusion

Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms ec. (m) posiion parameer dc/d concenraion gradien Fick s firs law saes ha he flu is proporional o he negaive concenraion gradien., diffusion coefficien is he proporionaliy consan. Minus sign comes from he fac ha maer flows down he concenraion gradien. I is no necessarily rue in all he cases. Maer may also diffuse up he concenraion gradien, which is called uphill diffusion. As can be seen in he ne module ha no he concenraion bu he chemical poenial gradien is imporan o know he direcion of flu.

Fick s firs law can direcly be applied in he case of seady sae, as shown in he eample below. Seady sae means ha here will no be any change in he composiion profile wih ime. If we pass carburizing gas hrough he pipe as shown in he figure and decarburizing gas ouside he pipe, a seady composiion profile may develop. oncenraion gradien can be calculaed following: dc d i d o o d i From his, one can calculae he flu from he known diffusion coefficiens of carbon or he diffusion coefficien from he flu deermined.

However, jus he Fick s firs law may no be useful o relae he flu and he concenraion profile, since, in mos of he cases, sysems develop non seady sae concenraion profile. I means ha concenraion a a paricular posiion changes wih ime. For eample, if we keep concenraion of carbon on one of he Fe surfaces anneal, composiion profile will change wih ime, as shown in he figure. We can apply Fick s firs law direcly o evaluae he concenraion profile ha develops during diffusion in his case, since, as shown in he figure, composiion changes wih annealing ime a a paricular posiion. So we need o have a relaion, which can relae ime along wih he concenraion and he posiion. For ha Fick s second law is derived. I is derived using conservaion of mass and Fick s firs law.

Le us consider a very hin slab in he maerial. J 1 is he incoming flu, J is he ougoing flu So he amoun of carbon coming in shor ime δ is J 1 δ (mole/m ) and going ou is J δ (here J 1 >J ) If he slab hickness is, hen δ ( J1 J ) δ Eq. 1 Furher, flu change in he hin slab can be considered linear and we can wrie J J J1 J1 J Eq. From Eq. 1 and Eq. J Using Fick s firs law ( J ):

Soluion o he Fick s second law Soluion of he Fick s second law depends on iniial and boundary condiions. Furher,, in some cases, may be considered consan. Tha means does no change wih concenraion or posiion. However, his condiion mees only in very limied cases. In mos of he pracical eamples, is a funcion of concenraion. In his case, soluion o he Fick s second law is complicaed. So in he beginning, we shall solve he Fick s second law for consan. Soluions are mainly for wo differen ypes of condiions, small and large ime values. When diffusion annealing ime is small, soluion is epressed wih inegrals or error funcions. For long annealing ime, i is epressed in erms of rigonomerical series. Noe ha he long or shor annealing ime is raher relaive. y saying long annealing ime, we mean ha he complee sample is affeced by he diffusion process and may lead o homogenizaion. y saying shor annealing ime, we mean ha eperimens are conduced such ha whole maerial is no affeced by he diffusion process.

Soluion for a hin film source Le us consider ha is consan. Tha means does no change wih he composiion. Since a a paricular locaion changes coninuously wih annealing ime, or a paricular changes is locaion coninuously. From he assumpion, we can sae ha is he same a any locaion. The meaning of he above saemen will be more clear, when we shall consider he change of wih he change in concenraion. Le us consider he siuaion, when a very hin film of elemen is sandwiched beween maerial A and hen annealed a elevaed emperaure.

Noe: One migh ask, how hin i is? y saying hin we mean ha he amoun of maerial is very low and even afer oal miing (full homogenizaion) elemen can be considered as impuriies in he maerial A. Tha means, afer deposiion of on A, he chemical poenial gradien is negligible. In oher sense, we can consider his as an eample of diffusion in he absence of any driving force (alhough small driving force will be here from enropy of miing) Acual meaning will be clear, when we discuss abou he aomic mechanism of diffusion! We can consider ha is consan, since he maerial in which i diffuses has fied composiion. For consan, he Fick s nd law can be wrien as As shown in he previous slide, i is seen ha he elemen disribuion afer he diffusion can be epressed by eponenial relaion. Following he Fick s second law he disribuion of elemen in A can be epressed as ( ) 0 ep 1/ where 4 o is a consan This relaion is developed based on he fac ha composiion profile shows eponenial decay from he hin film source.

The correcness of he soluion can be checked by differeniaion of he equaion wih respec o and and hen using hem in he Fick s second law o find he equal values on boh he sides. Furher, he boundary condiion ha 0, a a 0 and, a 0 a 0 also mee. Now one migh ge confused, when we say, since concenraion of elemen (X /V m 1/V m ) can be infinie anyime and will have some definie value even a compleely pure sae. Here is noional and means ha he elemen is no dilued by any oher elemen, which furher means ha i is in pure sae and for sysem i means ha i has infinie source of elemen before saring he eperimen. We shall show ha he absolue values of () or o are no imporan bu he raio ()/ o is imporan, where his soluion is acually used o calculae he racer or impuriy diffusion coefficien of he species.

Toal maerial (mole/m ) ha was sandwiched before annealing can be found following M + d So we can wrie M + 0 ep 1/ d 4 Furher, we consider λ means ( ) λ d d Since inegraion gives, we ge π M ( λ ) dλ 0 π + 0 ep M ep Replacing his for M, we find ( ) π 4

vs. describes he disribuion of elemens. d /d describes (following he Fick s law) he change in flu wih respec o disance a a paricular annealing ime. d d d d vs. I eplains he rae of change of elemen The negaive values indicae he region, where i loses he elemen and posiive value indicaes ha he region where i gains he elemen. Noe ha he region where i is losing or gaining he elemen depends on he ime of annealing.

The change in disribuion of elemens wih he increase in ime is shown in he figure. M ( ) ep π 4 Facor comes from he fac ha elemens diffuse boh he sides from he source

If he elemen diffuses o one direcion hen he facor should no be considered. M 4 ep ) ( π M 4 ln ) ( ln π From he calculaed slope, one can deermine he diffusion coefficien from he known annealing ime.