SUMMARY OF STOICHIOMETRIC RELATIONS AND MEASURE OF REACTIONS' PROGRESS AND COMPOSITION FOR MULTIPLE REACTIONS

Similar documents
between standard Gibbs free energies of formation for products and reactants, ΔG! R = ν i ΔG f,i, we

2.3 Nilpotent endomorphisms

CinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure

MEM 255 Introduction to Control Systems Review: Basics of Linear Algebra

G = G 1 + G 2 + G 3 G 2 +G 3 G1 G2 G3. Network (a) Network (b) Network (c) Network (d)

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

Formulas for the Determinant

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,

Quantum Mechanics for Scientists and Engineers. David Miller

DEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA. 1. Matrices in Mathematica

APPENDIX A Some Linear Algebra

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

y i x P vap 10 A T SOLUTION TO HOMEWORK #7 #Problem

Selected Student Solutions for Chapter 2

3. Be able to derive the chemical equilibrium constants from statistical mechanics.

Section 3.6 Complex Zeros

A new Approach for Solving Linear Ordinary Differential Equations

Difference Equations

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

= z 20 z n. (k 20) + 4 z k = 4

( ) 1/ 2. ( P SO2 )( P O2 ) 1/ 2.

1 Matrix representations of canonical matrices

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

1 Generating functions, continued

Lecture 12: Discrete Laplacian

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

1 GSW Iterative Techniques for y = Ax

Energy Storage Elements: Capacitors and Inductors

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

6.01: Introduction to EECS I Lecture 7 March 15, 2011

PRIMES 2015 reading project: Problem set #3

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to

However, since P is a symmetric idempotent matrix, of P are either 0 or 1 [Eigen-values

where I = (n x n) diagonal identity matrix with diagonal elements = 1 and off-diagonal elements = 0; and σ 2 e = variance of (Y X).

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Homework Assignment 3 Due in class, Thursday October 15

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Lecture 5 Decoding Binary BCH Codes

LECTURE 9 CANONICAL CORRELATION ANALYSIS

Chapter 7 Generalized and Weighted Least Squares Estimation. In this method, the deviation between the observed and expected values of

Canonical transformations

Errors for Linear Systems

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Expected Value and Variance

Estimation of the composition of the liquid and vapor streams exiting a flash unit with a supercritical component

Non-Ideality Through Fugacity and Activity

CHEMICAL REACTIONS AND DIFFUSION

EE215 FUNDAMENTALS OF ELECTRICAL ENGINEERING

8.6 The Complex Number System

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

5 The Rational Canonical Form

Foundations of Arithmetic

ME 300 Exam 2 November 18, :30 p.m. to 7:30 p.m.

TEST 5 (phy 240) 2. Show that the volume coefficient of thermal expansion for an ideal gas at constant pressure is temperature dependent and given by

Implicit Integration Henyey Method

Be true to your work, your word, and your friend.

NAME and Section No. it is found that 0.6 mol of O

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

Homework Notes Week 7

The KMO Method for Solving Non-homogenous, m th Order Differential Equations

Kernel Methods and SVMs Extension

Randić Energy and Randić Estrada Index of a Graph

Linear Feature Engineering 11

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Solution Thermodynamics

Linear Approximation with Regularization and Moving Least Squares

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014

Poisson brackets and canonical transformations

Chemical Equilibrium. Chapter 6 Spontaneity of Reactive Mixtures (gases) Taking into account there are many types of work that a sysem can perform

Srednicki Chapter 34

Chapter 3 Differentiation and Integration

18.781: Solution to Practice Questions for Final Exam

From Biot-Savart Law to Divergence of B (1)

Effects of Ignoring Correlations When Computing Sample Chi-Square. John W. Fowler February 26, 2012

MAE140 - Linear Circuits - Fall 13 Midterm, October 31

Outlet temperature of a WGS reactor (Stage I) for the conversion of CO, applied for the abatement of CO to a fixed value.

Cyclic Codes BCH Codes

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Norms, Condition Numbers, Eigenvalues and Eigenvectors

The Fundamental Theorem of Algebra. Objective To use the Fundamental Theorem of Algebra to solve polynomial equations with complex solutions

GAUTENG DEPARTMENT OF EDUCATION SENIOR SECONDARY INTERVENTION PROGRAMME PHYSICAL SCIENCES GRADE 12 SESSION 1 (LEARNER NOTES)

Lecture 16 Statistical Analysis in Biomaterials Research (Part II)

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b

More metrics on cartesian products

Problem Set 9 Solutions

FINITELY-GENERATED MODULES OVER A PRINCIPAL IDEAL DOMAIN

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

International Mathematical Olympiad. Preliminary Selection Contest 2012 Hong Kong. Outline of Solutions

Publication 2006/01. Transport Equations in Incompressible. Lars Davidson

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Transcription:

UMMAY OF TOICHIOMETIC ELATION AND MEAUE OF EACTION' POGE AND COMPOITION FO MULTIPLE EACTION UPDATED 0/4/03 - AW APPENDIX A. In case of multple reactons t s mportant to fnd the number of ndependent reactons. Ths ermnes for how many speces the change n the number of moles may occur ndependently,.e for how many speces we must montor the change n the number of moles caused by reactons n order to compute the composton of the whole system (.e to compute the changes n the number of moles of all other speces). Usually the problem s stated as reactons between speces: A 0;,.. where s the stochometrc coeffcent of speces n reacton (ecall > 0 products, < 0 reactants). Then: a) Fnd the number of ndependent reactons. b) Fnd the molar etents of the ndependent reactons ( ś X,,...) n terms of the measured changes n the number of moles,, of key speces. ( F F o,,,...) c) Epress the changes n moles of the remanng - speces, (,,...) n terms of changes n moles of the selected key speces. d) Determne the etents of the remanng (dependent) reactons n terms of the calculated etents (f so desred). Frst we wll ntroduce the classcal procedure and then show a short-cut method whch accomplshes all of the aboe n a sngle step.

Classcal Method: Consder the matr of stochometrc coeffcents of all orgnal reacton equatons and speces (.e. the (' X ) matr): orgnal ( X )...... 5... 5 5 and fnd ts rank,.e, the largest nonzero ermnant. Use for eample Gaussan elmnaton (see for eample E. Kreyszg "Adanced Engneerng Mathematcs" Fourth Edton, Wley, N.Y. 979, Chapter 7-6). Ths procedure wll generate - rows of zeros when s the rank of the matr and the number of ndependent reactons. All reactons for whch rows of zeros are generated are the dependent ones. The rest s pcked as a set of ndependent reactons for whch the matr of stochometrc coeffcents s: new ( X)............ Please understand that reactons hae now been renumbered, snce some were deleted as beng dependent, and the elements of the new matr of stochometrc coeffcents do not necessarly correspond n order one to one to the elements of the old matr. Now usng the law of stochometry we can wrte: T ( ) X ś ( X)X ( ) () where T s the transpose of the matr.e., s matr wth columns and rows nterchanged.

Epanded ew of eq () s: th row.................................... X& X&...... X&...... X& (a) Consderng now ust the -th row of the ector, and rememberng that matr multplcaton s rows nto columns, we get: X ś X ś... ś F { number of moles of produced n the reacton system 3 Σ oer all ndependent reactons X ν ś 3 X number of moles of produced by reacton (b) We can sole equaton () and we wll show a formal soluton n terms of Cramer's rule. (Ths of course s not practcal and actual soluton should proceed by Gaussan elmnaton or Gauss Jordan technque). The soluton of () s: where T T ( ) X& () () s the nerse matr of ( T ) how to obtan t by Gauss Jordan technque]. [see Kreyszg or some other lnear algebra tet to fnd nce the alue of the ermnant s unchanged f rows and columns are nterchanged, the soluton n terms of Cramer's rule s: ( ) () X & F,,... (3) 3

where ()......... ν (4)......... and () F s the ( ) where the -th row s replaced by F,, etc. ( )......... ν... (5)...... Ths soles parts (a), of fndng the number of ndependent reactons and (b), of relatng etents of ndependent reactons to measured changes n moles of key speces. We hae related etents, X, to the changes n moles of speces, F,,.... If < we need to epress the changes n moles of the remanng - speces n terms of the changes n moles of the frst speces (part c). tochometry dctates that eq (b) s ald for all speces, not ust for the frst speces for whch we hae used t so far. Hence, we can wrte k k X& ; k,,..., (6) Ths ges the change n moles of the remanng - speces n terms of the now known etents. We can proceed a step further and substtute equaton (3) nto equaton (6) to obtan: k k ( ) () F (7) 4

After rearrangement we recognze the followng: 0 k () k () () F ( ) k (8) The left hand sde of the aboe equaton s ust an epanson n mnors of the ermnant on the rght hand sde when the epanson s done about the elements of the last ( ) column. If we deelop ths ermnant wth respect to the frst row (the alue of the ermnant s properly zero ndcatng that F k, k,,..., s lnearly dependent on the frst s ) we get: ' k () ( ) ) 0 (9) k whch fnally yelds k ( ) ) k () n terms of the frst k ; k,..., where ( ) ) s ( ) k F s : k, k wth the -th column replaced by,.... k (0) throw ( ) ) k... k k k......... Note that n prong eq (8) and (9) one uses the property that the alue of a ermnant must be multpled by (-) each tme two rows or two columns are nterchanged. Our results can be summarzed as follows: The etents are obtaned from eq (3) whch now becomes () () () ( ) () F X& ( ) ;,... (3a) (3a) s () wth row and column deleted. The change n moles for the non-key speces s gen by: 5

( ) k () k ; k,... (0) hortcut Method: Ths method ges all the nformaton n one sngle sweep. Consder the augmented matr T F orgnal ( )........................ Use Gauss-Jordan technque on ths augmented matr to generate 's on the dagonal and zeros elsewhere by normalzng the pot element and by generatng. The procedure stops when we cannot generate any more 's on the dagonal and the rows below the last are flled wth zeros. emember, we can always nterchange rows, keepng track to renumber our reactons accordngly, to get the rows wth zeros on the bottom. When we are fnshed wth the Gauss-Jordan procedure the augmented matr has the form: k 6

rows wth ' s on dagonal dagonal untary matr ermnes No. of ndependent reactons 64748 0 0...0 00...0 0 0...0... 0 0 0... Coeffcents for eactons 64748 λ... λ λ λ λ 3... λ... λ...... λ 3... α α α etents rows of zeros 0 0 0...0 0 0 0...0... 0 0 0...0 0...0 0...0... 0...0... β β β narants In the upper left hand corner the ( X ) matr wth 's on the dagonal and zeros eerywhere else ermnes the rank and the number of ndependent reactons. The frst rows on the top n the rght hand corner, to the rght of the partton, ge the etents: X& α ;,... () where α hae been ermned by the procedure. By comparson wth eq (3a) we see that the α generated by the procedure must be ( ) () α ( ) ;,,..., ;,,..., 7

The lower rght corner,.e the last - rows to the rght of the partton, defnes the ( - ) narants of the system: k β k 0; k,,... () These lnear combnatons of ' s (narances) must be zero. Otherwse the last - rows would not hae a 0 0 relaton and the system would be mpossble and hae no soluton, meanng that a mass balance was olated wth a wrong stochometrc coeffcent someplace. Comparson of equaton () and equaton (0) ges: ( ) k β k ; k,,..., ;,,..., (3) () Note that β k hae been calculated by the Gauss-Jordan procedure. Fnally, the ( X( ) matrof λ' s, squeezed between the dagonal ( X ) dentty matr and the partton, yelds the coeffcents for epressng the etents of dependent reactons n terms of the ndependent ones: X& k λ X& k k, ;,..., (4) Let us document ths by an eample: PECIE 3 4 C CO C O O O CO CO CO C O CO CO ( ) ( ) ( 3) ( 4) () 0 ( ) 0 () 3 0 orgnal 0 0 0 8

Use short-cut method. Wrte the transpose matr: T orgnal 0 0 0 T Then proceed wth the augmented matr ( orgnal ) as follows: 0 0 3 0 4 0 0 0 3 0 4 0 0 0 0 0 3 0 0 0 4 etents narants There are (two) ndependent reactons. Ther etents n terms of molar changes of the frst two speces are:?x ( F F 0 ) ( F 0 F )?X ( F 0 F ) ( F 0 F ) There are ndependent reactons and ther etents are gen aboe. There are - 4 - two narants. The changes n the number of moles of speces 3 and 4 can be epressed n terms of the change n moles of speces and usng the two narants as shown below: 3 ( ) ( F F 0 F F 0 ) 4 ( F 0 F ) F 0 F ( ) 9

The etent of the thrd reacton s: etent of the dependent reacton { X ś 3 X ś X ś.e the thrd reacton s the sum of the frst two dded by two. Ths n ths smple case s obous by nspecton. Check the classcal method results. By ths method we would hae frst ermned and would hae taken the frst two reactons as ndependent. The X X matr of stochometrc coeffcents usng the frst two components as the key ones, s: 0 ; () 0 () () ( ) ; ( ) ; ( ) 0 () ( ) 0 ( ) Then α ( ) α ( ) ; a ( ) 3 0 ; α ( ) 4 0 Hence:?X?X 0 QED checks QED Now: β 3 0 β 3 4 4 0 ; β 4 0 0 ; β 4 3 0; 4 0 checks Once we know the number of ndependent reactons and ther etents we can wrte for all speces,,..., 0

F F o ν? X ;,,... (5) F tot F F tot,o ν? X ; ν ν (6) Another hortcut Method: Another approach to the aboe problem s to count all the chemcal elements present n the reacton system under consderaton (say total number of elements s E) and form the matr of coeffcents E so that ν s the number of atomc speces present n the chemcal speces. The matr E s a (EX) matr. One should fnd the rank of matr E, say E. Then the number of ndependent reactons s equal to the total number of chemcal speces mnus the rank of the matr E,.e. E (7)