The theoretical background of

Similar documents
Some Different Perspectives on Linear Least Squares

Applying the condition for equilibrium to this equilibrium, we get (1) n i i =, r G and 5 i

A New Method for Solving Fuzzy Linear. Programming by Solving Linear Programming

2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission

Algorithms behind the Correlation Setting Window

Coherent Potential Approximation

Chemistry 163B Introduction to Multicomponent Systems and Partial Molar Quantities

Chemistry 163B Introduction to Multicomponent Systems and Partial Molar Quantities

Stationary states of atoms and molecules

CS 2750 Machine Learning. Lecture 7. Linear regression. CS 2750 Machine Learning. Linear regression. is a linear combination of input components x

Non-degenerate Perturbation Theory

KURODA S METHOD FOR CONSTRUCTING CONSISTENT INPUT-OUTPUT DATA SETS. Peter J. Wilcoxen. Impact Research Centre, University of Melbourne.

We have already referred to a certain reaction, which takes place at high temperature after rich combustion.

7.0 Equality Contraints: Lagrange Multipliers

Department of Mechanical Engineering ME 322 Mechanical Engineering Thermodynamics. Ideal Gas Mixtures. Lecture 31

Some results and conjectures about recurrence relations for certain sequences of binomial sums.

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II

General Method for Calculating Chemical Equilibrium Composition

ROOT-LOCUS ANALYSIS. Lecture 11: Root Locus Plot. Consider a general feedback control system with a variable gain K. Y ( s ) ( ) K

DIFFERENTIAL GEOMETRIC APPROACH TO HAMILTONIAN MECHANICS

A Model Reduction Technique for linear Model Predictive Control for Non-linear Large Scale Distributed Systems

Standard Deviation for PDG Mass Data

The Mathematics of Portfolio Theory

Solutions to problem set ); (, ) (

Capacitated Plant Location Problem:

Functions of Random Variables

A Conventional Approach for the Solution of the Fifth Order Boundary Value Problems Using Sixth Degree Spline Functions

PRACTICAL CONSIDERATIONS IN HUMAN-INDUCED VIBRATION

Introduction. Free Electron Fermi Gas. Energy Levels in One Dimension

D. L. Bricker, 2002 Dept of Mechanical & Industrial Engineering The University of Iowa. CPL/XD 12/10/2003 page 1

Fig (1) A Distillation column

1 Lyapunov Stability Theory

CS 1675 Introduction to Machine Learning Lecture 12 Support vector machines

ENGI 4421 Propagation of Error Page 8-01

3.1 Introduction to Multinomial Logit and Probit

The Mathematical Appendix

DYNAMICS. Systems of Particles VECTOR MECHANICS FOR ENGINEERS: Seventh Edition CHAPTER. Ferdinand P. Beer E. Russell Johnston, Jr.

Support vector machines II

Homework Assignment Number Eight Solutions

1 Onto functions and bijections Applications to Counting

Objectives of Multiple Regression

å 1 13 Practice Final Examination Solutions - = CS109 Dec 5, 2018

DATA DOMAIN DATA DOMAIN

STATISTICAL PROPERTIES OF LEAST SQUARES ESTIMATORS. x, where. = y - ˆ " 1

16 Homework lecture 16

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions.

Parallelized methods for solving polynomial equations

Chapter 9 Jordan Block Matrices

The E vs k diagrams are in general a function of the k -space direction in a crystal

An Introduction to. Support Vector Machine

Relations to Other Statistical Methods Statistical Data Analysis with Positive Definite Kernels

Binary classification: Support Vector Machines

are positive, and the pair A, B is controllable. The uncertainty in is introduced to model control failures.

( ) ( ) ( ( )) ( ) ( ) ( ) ( ) ( ) = ( ) ( ) + ( ) ( ) = ( ( )) ( ) + ( ( )) ( ) Review. Second Derivatives for f : y R. Let A be an m n matrix.

THE ROYAL STATISTICAL SOCIETY 2016 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE MODULE 5

Global Optimization for Solving Linear Non-Quadratic Optimal Control Problems

3D Reconstruction from Image Pairs. Reconstruction from Multiple Views. Computing Scene Point from Two Matching Image Points

MOLECULAR VIBRATIONS

EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM

Third handout: On the Gini Index

Module 7. Lecture 7: Statistical parameter estimation

( ) 2 2. Multi-Layer Refraction Problem Rafael Espericueta, Bakersfield College, November, 2006

2.28 The Wall Street Journal is probably referring to the average number of cubes used per glass measured for some population that they have chosen.

Bounds on the expected entropy and KL-divergence of sampled multinomial distributions. Brandon C. Roy

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

12.2 Estimating Model parameters Assumptions: ox and y are related according to the simple linear regression model

A Helmholtz energy equation of state for calculating the thermodynamic properties of fluid mixtures

Theory study about quarter-wave-stack dielectric mirrors

COV. Violation of constant variance of ε i s but they are still independent. The error term (ε) is said to be heteroscedastic.

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:

Econometric Methods. Review of Estimation

arxiv:math/ v2 [math.co] 2 Feb 1999

SUBCLASS OF HARMONIC UNIVALENT FUNCTIONS ASSOCIATED WITH SALAGEAN DERIVATIVE. Sayali S. Joshi

Unit 9. The Tangent Bundle

Lecture 9: Tolerant Testing

CHAPTER VI Statistical Analysis of Experimental Data

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin

C.11 Bang-bang Control

Lecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

Unsupervised Learning and Other Neural Networks

On the energy of complement of regular line graphs

A Remark on the Uniform Convergence of Some Sequences of Functions

III-16 G. Brief Review of Grand Orthogonality Theorem and impact on Representations (Γ i ) l i = h n = number of irreducible representations.

Multiple Choice Test. Chapter Adequacy of Models for Regression

The number of observed cases The number of parameters. ith case of the dichotomous dependent variable. the ith case of the jth parameter

THE TRUNCATED RANDIĆ-TYPE INDICES

LINEAR REGRESSION ANALYSIS

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations

Fibonacci Identities as Binomial Sums

L5 Polynomial / Spline Curves

Summary of the lecture in Biostatistics

A New Measure of Probabilistic Entropy. and its Properties

7. Equilibrium Statistical Mechanics

Chapter 5 Properties of a Random Sample

Integral Equation Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, Xin Wang and Karen Veroy

Mu Sequences/Series Solutions National Convention 2014

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

Ideal multigrades with trigonometric coefficients

Transcription:

he theoretcal backgroud of

-echologes he theoretcal backgroud of FactSage he followg sldes gve a abrdged overvew of the ajor uderlyg prcples of the calculatoal odules of FactSage.

-echologes he bbs Eergy ree Mawell H, U, F Phase Dagra,c p(), H (),S (),a,v bbs-duhe Legedre rasfor. Equlbra Partal Dervatves wth Respect to, or P Msato bbs Eergy Matheatcal ethods are used to derve ore forato fro the bbs eergy ( of phase(s) or whole systes ) Matheatcal Method Calculatoal result derved fro

-echologes herodyac potetals ad ther atural varables Varables bbs eergy: = (, p,,...) Ethalpy: H = H (S, P,,...) Free eergy: A = A (,V,,...) Iteral eergy: U = U (S, V,,...) Iterrelatoshps: A = U S H = U PV = H S = U PV S

-echologes P µ, V A, SP H, SV U, Mawell-relatos: herodyac potetals ad ther atural varables V P H S P P S U V H A S V ad

-echologes... S,V, cost. for 0, du U...,p, cost. for 0, d herodyac potetals ad ther atural Equlbru codto:... U,, cost. for 0, d A... S,p, cost. for 0, dh H... cost.u, V, for 0, ds S a

-echologes eperature Coposto p p p p p H c S H S, 2 2,,, Use of odel equatos perts to start at ether ed! bbs-duhe tegrato Partal Operator Itegral quatty:, H, S, c p Partal quatty: µ, h, s, c p() herodyac propertes fro the bbs-eergy

-echologes herodyac propertes fro the bbs-eergy J.W. bbs defed the checal potetal of a copoet as:,p Wth ( s a etesve property!) oe obtas

-echologes rasforato to ole fractos : 1 = partal operator herodyac propertes fro the bbs-eergy C p c p C p C p H H H h S s S S

-echologes bbs eergy fucto for a pure substace () (.e. eglectg pressure ters) s calculated fro the ethalpy H() ad the etropy S() usg the well-kow bbs-helholtz relato: H S I ths H() s H H 298 298 cp d ad S() s S S 298 298 c p d hus for a gve -depedece of the c p -polyoal (for eaple after Meyer ad Kelley) oe obtas for (): () A B C l D 2 E 3 F 2

-echologes bbs eergy fucto for a soluto As show above (,) for a soluto cossts of three cotrbutos: the referece ter, the deal ter ad the ecess ter. For a sple substtutoal soluto (oly oe lattce ste wth rado occupato) oe obtas usg the well-kow Redlch-Kster-Muggau polyoal for the ecess ters: (, ) o, R l j 0 j j L ( ) j ( )( j ) j k j k ( L jk ( ) j L jk j ( ) k L jk k ( )) /( j k )

-echologes Equlbru cosderatos a) Stochoetrc reactos Equlbru codto: Reacto : A A + B B +... = S S + +... eerally : or d 0 B 0 For costat ad p,.e. d = 0 ad dp = 0, ad o other work ters: d d

-echologes Equlbru cosderatos a) Stochoetrc reactos For a stochoetrc reacto the chages d are gve by the stochoetrc coeffcets ad the chage eted of reacto d. d d hus the proble becoes oe-desoal. Oe obtas: d d 0 [see the followg graph for a eaple of = () ]

bbs eergy -echologes Equlbru cosderatos a) Stochoetrc reactos = 400K = 500K = 550K 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 Etet of Reacto bbs Eergy as a fucto of etet of the reacto 2NH 3 <=>N 2 + 3H 2 for varous teperatures. It s assued, that the chages of ethalpy ad etropy are costat.

-echologes Equlbru cosderatos a) Stochoetrc reactos Separato of varables results : d dξ µ 0 hus the equlbru codto for a stochoetrc reacto s: µ 0 Itroducto of stadard potetals ad actvtes a yelds: µ µ R l a Oe obtas: µ R l a 0

-echologes Equlbru cosderatos a) Stochoetrc reactos It follows the Law of Mass Acto: µ R l a where the product K a or s the well-kow Equlbru Costat. he REACION odule perts a ulttude of calculatos whch are based o the Law of Mass Acto. K ep R

-echologes Equlbru cosderatos b) Mult-copoet ult-phase approach Cople Equlbra May copoets, ay phases (soluto phases), costat ad p : wth o R l a or p

-echologes Equlbru cosderatos b) Mult-copoet ult-phase approach Massbalace costrat a j b j j = 1,..., of copoets b Lagragea Multplers M j tur out to be the checal potetals of the syste copoets at equlbru: j b j M j

-echologes Equlbru cosderatos b) Mult-copoet ult-phase approach a j j Phase as Slag Lq. Fe C o po et s Syst e C o po et s Fe N O C Ca S Mg Fe 1 0 0 0 0 0 0 N 2 0 2 0 0 0 0 0 O 2 0 0 2 0 0 0 0 C 0 0 0 1 0 0 0 CO 0 0 1 1 0 0 0 CO 2 0 0 2 1 0 0 0 Ca 0 0 0 1 0 0 0 C ao 0 0 1 0 1 0 0 S 0 0 0 0 0 1 0 SO 0 0 1 0 0 1 0 Mg 0 0 0 0 0 0 1 SO 2 0 0 2 0 0 1 0 Fe 2 O 3 2 0 3 0 0 0 0 C ao 0 0 1 0 1 0 0 FeO 1 0 1 0 0 0 0 M go 0 0 1 0 0 0 1 Fe 1 0 0 0 0 0 0 N 0 1 0 0 0 0 0 O 0 0 1 0 0 0 0 C 0 0 0 1 0 0 0 Ca 0 0 0 0 1 0 0 S 0 0 0 0 0 1 0 Mg 0 0 0 0 0 0 1 Eaple of a stochoetrc atr for the gas-etal-slag syste Fe-N-O-C-Ca-S-Mg

-echologes Equlbru cosderatos b) Mult-copoet ult-phase approach Modellg of bbs eergy of (soluto) phases,,p Pure Substace o, o, (, p) (stochoetrc) Soluto phase, ref, d, e S d Choose approprate referece state ad deal ter, the check for devatos fro dealty. See Page 11 for the sple substtutoal case.

-echologes Equlbru cosderatos Mult-copoet ult-phase approach Use the EQUILIB odule to eecute a ulttude of calculatos based o the cople equlbru approach outled above, e.g. for cobusto of carbo or gases, aqueous solutos, etal clusos, gas-etal-slag cases, ad ay others. NOE: he use of costrats such calculatos (such as fed heat balaces, or the occurrece of a predefed phase) akes ths odule eve ore versatle.

-echologes Phase dagras as projectos of bbs eergy plots Hllert has poted out, that what s called a phase dagra s dervable fro a projecto of a so-called property dagra. he bbs eergy as the property s plotted alog the z-as as a fucto of two other varables ad y. Fro the u codto for the equlbru the phase dagra ca be derved as a projecto oto the -y-plae. (See the followg graphs for llustratos of ths prcple.)

-echologes Phase dagras as projectos of bbs eergy plots g a b g b a b a b g a P P Uary syste: projecto fro --p dagra

-echologes Phase dagras as projectos of bbs eergy plots 1.0 0.5 0.0-0.5-1.0 300 400 1.0 N 0.8 0.6 0.4 0.2 N 0.0 Cu 700 600 500 Bary syste: projecto fro -- dagra, p = cost.

-echologes Phase dagras as projectos of bbs eergy plots erary syste: projecto fro - 1-2 dagra, = cost ad p = cost

-echologes Phase dagras geerated wth FactSage Use the PHASE DIARAM odule to geerate a ulttude of phase dagras for uary, bary, terary or eve hgher order systes. NOE: he PHASE DIARAM odule perts the choce of, P, (as R l a), a (as l a), ol () or weght (w) fracto as as varables. Mult-copoet phase dagras requre the use of a approprate uber of costats, e.g. a terary sopleth dagra vs oe olar rato has to be kept costat.

-echologes N-Copoet Syste (A-B-C- -N) Etesve varables q S V A B N U q -P µ A µ B µ N q j Correspodg potetals d U d S P d V d d q bbs-duhe: S d V d P d q d 0

-echologes Choce of Varables whch always gve a rue Phase Dagra N-copoet syste (1) Choose potetals: 1, 2,, (2) Fro the o-correspodg etesve varables (q +1, q +2, ), for (N+1-) depedet ratos (Q +1, Q +2,, Q N+1 ). N 1 Eaple: Q j q N2 J 1 q j 1 N 1 [ 1, 2,, ; Q +1, Q +2,, Q N+1 ] are the the (N+1) varables of whch 2 are chose as aes ad the reader are held costat.

-echologes MgO-CaO Bary Syste Etesve varables ad correspodg potetals S Chose aes varables ad costats 1 = for y-as V -P 2 = -P costat MgO CaO µ MgO µ CaO q q 3 4 MgO CaO Q 3 MgO CaO CaO for -as

-echologes Fe - Cr - S - O Syste S f 1 = (costat) V -P f 2 = -P (costat) O 2 O 2 3 O 2 -as S 2 S 2 4 S 2 -as Fe Cr Fe Cr q q 5 6 Cr Fe Q 5 Cr Fe (costat)

-echologes Fe - Cr - C Syste - proper choce of aes varables S f 1 = (costat) V -P f 2 = -P (costat) C Fe C Fe f 3 = C a C for -as ad Q 4 for y-as Q 4 Cr F e C r C (NO OK) Cr Cr Q 4 F e Cr C r (OK) dq j Requreet: 0 for 3 dq

Mole fracto of Cr -echologes Fe - Cr - C Syste - proper choce of aes varables 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 M 23 C 6 M 7 C 3 bcc fcc -3-2 -1 0 1 2 log(a c ) ceette hs s NO a true phase dagra. Reaso: C ust NO be used forula for ole fracto whe a C s a as varable. NOE: FactSage users are safe sce they are ot gve ths partcular choce of aes varables.