Solution in semi infinite diffusion couples (error function analysis)

Similar documents
J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

Scattering at an Interface: Oblique Incidence

( ) () we define the interaction representation by the unitary transformation () = ()

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

Notes on the stability of dynamic systems and the use of Eigen Values.

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

Mechanics Physics 151

Linear Response Theory: The connection between QFT and experiments

Mechanics Physics 151

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

Comb Filters. Comb Filters

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

P R = P 0. The system is shown on the next figure:

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

National Exams December 2015 NOTES: 04-BS-13, Biology. 3 hours duration

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

Advanced Machine Learning & Perception

CHAPTER 5: MULTIVARIATE METHODS

TSS = SST + SSE An orthogonal partition of the total SS

Lecture 11 SVM cont

How about the more general "linear" scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )?

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :

Transcription: Messenger RNA, mrna, is produced and transported to Ribosomes

Mechanics Physics 151

On One Analytic Method of. Constructing Program Controls

Chapter Lagrangian Interpolation

Polymerization Technology Laboratory Course

II. Light is a Ray (Geometrical Optics)

PHYS 705: Classical Mechanics. Canonical Transformation

3. OVERVIEW OF NUMERICAL METHODS

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading

Volatility Interpolation

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data

Chapters 2 Kinematics. Position, Distance, Displacement

Robust and Accurate Cancer Classification with Gene Expression Profiling

Chapter 6: AC Circuits

Displacement, Velocity, and Acceleration. (WHERE and WHEN?)

FI 3103 Quantum Physics

Graduate Macroeconomics 2 Problem set 5. - Solutions

Part II CONTINUOUS TIME STOCHASTIC PROCESSES

Chapter 2 Linear dynamic analysis of a structural system

CHAPTER 10: LINEAR DISCRIMINATION

Sampling Procedure of the Sum of two Binary Markov Process Realizations

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

PHYS 1443 Section 001 Lecture #4

Variants of Pegasos. December 11, 2009

Let s treat the problem of the response of a system to an applied external force. Again,

Lecture 2 M/G/1 queues. M/G/1-queue

Relative controllability of nonlinear systems with delays in control

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)

On elements with index of the form 2 a 3 b in a parametric family of biquadratic elds

, t 1. Transitions - this one was easy, but in general the hardest part is choosing the which variables are state and control variables

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

Example: MOSFET Amplifier Distortion

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model

Density Matrix Description of NMR BCMB/CHEM 8190

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005

Robustness Experiments with Two Variance Components

Diffusion of Heptane in Polyethylene Vinyl Acetate: Modelisation and Experimentation

Department of Economics University of Toronto

. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.

Pendulum Dynamics. = Ft tangential direction (2) radial direction (1)

Math 128b Project. Jude Yuen

5-1. We apply Newton s second law (specifically, Eq. 5-2). F = ma = ma sin 20.0 = 1.0 kg 2.00 m/s sin 20.0 = 0.684N. ( ) ( )

[Link to MIT-Lab 6P.1 goes here.] After completing the lab, fill in the following blanks: Numerical. Simulation s Calculations

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

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

Born Oppenheimer Approximation and Beyond

Advanced time-series analysis (University of Lund, Economic History Department)

Testing a new idea to solve the P = NP problem with mathematical induction

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts

Nonequilibrium models for a multi component reactive distillation column

Bernoulli process with 282 ky periodicity is detected in the R-N reversals of the earth s magnetic field

. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue.

CS286.2 Lecture 14: Quantum de Finetti Theorems II

Density Matrix Description of NMR BCMB/CHEM 8190

2/20/2013. EE 101 Midterm 2 Review

Comparison of Differences between Power Means 1

Molecular Dynamics Simulation Study forgtransport Properties of Diatomic Liquids

COMPUTER SCIENCE 349A SAMPLE EXAM QUESTIONS WITH SOLUTIONS PARTS 1, 2

SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β

Homework 8: Rigid Body Dynamics Due Friday April 21, 2017

Cubic Bezier Homotopy Function for Solving Exponential Equations

A Competitive Test for Uniformity of Monotone Distributions

Introduction to Boosting

January Examinations 2012

FTCS Solution to the Heat Equation

Computing Relevance, Similarity: The Vector Space Model

SELFSIMILAR PROCESSES WITH STATIONARY INCREMENTS IN THE SECOND WIENER CHAOS

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Institute for Mathematical Research, Universiti Putra Malaysia, UPM Serdang, Selangor, Malaysia

Lecture Notes 4. Univariate Forecasting and the Time Series Properties of Dynamic Economic Models

Summary: Solving a Homogeneous System of Two Linear First Order Equations in Two Unknowns

Transcription:

Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of pure A and an alloy of A. Unlke he prevous case, here, because of he dfference n he composon, dffuson wll be drven by he chemcal poenal graden. However, we shall show laer ha he soluon we are gong o derve, can be used only n he case where he concenraon and he chemcal poenal dfference of he end members s no much. Tha means dffuson coeffcens do no vary sgnfcanly wh he composon.

Sem nfne means ha, we anneal for a ceran annealng me such ha he end of he nal maerals are no affeced by he dffuson of elemens. Alhough he meanng of he sem nfne ndcaes ha a good par of he end of he dffuson couple should no be affeced bu acually even wh one unaffeced aomc layer n he end s suffcen o consder he sysem as sem nfne dffuson couple! Ths s mporan snce oherwse we can apply he relaon derved here o deermne he dffuson parameer or o calculae he concenraon profle from he known dffuson parameer. For he sake of solvng he Fck s second law, we dvde he alloy wh he composon no n numbers of very hn slces. We consder he dsrbuon of elemen from one slce. Suppose from he slce,, whch s a he dsance of from he nal conac plane,. Then he dsrbuon of elemen followng he soluon for hn flm, we can wre ( ) ep 4 where M ( mole / m )

d 4 ep ) ( Smlarly, f we consder he dsrbuon of elemens from all he slces and supermpose, we shall ge he oal dsrbuon of elemen a a parcular poson as n 4 ep Furher consder Snce we have consdered very hn slces, we can wre d d ep ep d d Tha means y replacng When and when

Snce hs negraon s no very sraghforward and akes very long me o deermne one sngle value, s epressed n erms of error funcon. Advanage of usng hs error funcon s ha one can deermne he values from a able. Error funcon of Z s defned as z ( z) ep( ) d Noe here ha z s acually he upper lm of he negraon. Lower lm of he negraon should be zero. So he epresson n he prevous slde should be rewren such ha can be replaced by he error funcon ( ) ep( ) d ep( ) d Followng, we can wre n erms of error funcon as ( ) snce (-z) (z)

So from he prevous relaon, s apparen ha he sgn of wll depend on whch sde of he, we are neresed o calculae. Furher, f he composon profle s jus he oppose compared o he frs eample, can be shown ha he relaon becomes Noe: We always use a superscrp for he concenraon o denoe rgh hand sde of he couple and - for he lef hand sde of he couple.

Tabulaed error funcon values z ( z) ep( ) d 3 5 7 z z z z... 3! 5! 73! ( z) ( z)

In he prevous case, we consdered he dffuson beween he dffuson couple of wo dfferen blocks wh composon and or oherwse - and. Tha means n boh he cases, was a couple beween one pure elemen and an alloy of A. Error funcon soluons gven prevously can be rewren nerms of normalzed concenraon profles as In some cases, s possble ha wo alloys of A are couples. Now f we consder he dffuson couple beween - and, where > -, he relaon can be wren as shown n Fg. a he relaon can be wren as shown n Fg. a or

Furher, he relaon n he case of dffuson couple, where - and and - > as shown n Fg. b (n he prevous slde), can be epressed as or In he case of carburzng, we can see ha he relaon wll be he same, rrespecve of he sde on whch carburzng s done, snce n he case of Fg. a, s posve bu n Fg. b, s negave o S o S S o S o S For o

Few mporan noes So from he prevous relaons, s apparen ha he sgn of (wheher posve or negave) wll depend on n whch sde of he, we are neresed o calculae. Noe agan ha wll be more or less he consan a any poson s calculaed. I was also one of he assumpons for hs dervaon. From he error funcon analyss above, we can see ha, a from he relaon, comes ( ) from he relaon, S ( S ) comes ( ) S ( ) ( ) I ndcaes ha n a sysem, where wo blocks wh dfferen concenraons are coupled, a, whch s bascally he nal conac plane, he concenraon wll be always average of he concenraon of he nal maerals. Noe ha we need o locae anyway he nal conac plane afer he dffuson process, snce n he equaon s acually measured from he locaon of he nal conac plane. Ths also ndcaes ha hs equaon can only be used when end pars of he couple s no affeced by he dffuson process, snce oherwse a, wll have anoher average value dependng on how much of he end members are affeced.

P P Q Q In he case of carburzng, he surface concenraon should always be he same. Noe agan ha wll be more or less he consan a any poson s calculaed. I was also one of he assumpons for hs dervaon. efore proceedng furher for calculaon usng error funcon, one needs o check f ndeed you can locae he poson of by akng average of he end member concenraons. Ths poson should be possble o locae even by equalzng he areas P and Q (snce he loss from one sde of he couple should be equal o he gan n anoher sde). If he composon profle s no symmerc, can be found mahemacally by d ( )d When for a parcular he negrals have he same values, should be fed as poson.

Alhough looks very rval, fndng hs plane may no be ha easy n mos of he cases. Snce nal conac plane s an unque plane (only one value s possble), should be possble o fnd from any of he concenraon profles of elemens A or, as shown n he prevous slde. However, from our eperence, we have seen ha s possble o locae or fnd a sngle value from any of he profles, only when molar volume changes deally n he sysem, as shown n he fgure above. Noe ha he molar volume daa are requred o calculae he concenraon. When devaes (posvely or negavely) from he dealy, one fnds wo dfferen values from wo dfferen composon profles. fference depends on he devaon from he dealy. I feels as f he sysem los s nal conac plane. I wll be dscussed furher n deals a laer sage, wh eamples. So he lmaon of he error funcon analyss s no only ha dffuson coeffcen s more or less he consan bu also he molar volume does no change much wh he composon or changes deally.