General Method for Calculating Chemical Equilibrium Composition

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

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

Idea is to sample from a different distribution that picks points in important regions of the sample space. Want ( ) ( ) ( ) E f X = f x g x dx

Numerical Analysis Formulae Booklet

Mu Sequences/Series Solutions National Convention 2014

FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM

Homework Assignment Number Eight Solutions

UNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS

ENGI 4421 Propagation of Error Page 8-01

G S Power Flow Solution

MOLECULAR VIBRATIONS

Newton s Power Flow algorithm

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

MATH 371 Homework assignment 1 August 29, 2013

2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America

Logistic regression (continued)

Midterm Exam 1, section 1 (Solution) Thursday, February hour, 15 minutes

Midterm Exam 1, section 2 (Solution) Thursday, February hour, 15 minutes

CHEMICAL EQUILIBRIA BETWEEN THE ION EXCHANGER AND GAS PHASE. Vladimir Soldatov and Eugeny Kosandrovich

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

means the first term, a2 means the term, etc. Infinite Sequences: follow the same pattern forever.

THE ROYAL STATISTICAL SOCIETY 2016 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE MODULE 5

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

Physics 114 Exam 2 Fall Name:

Chapter 9 Jordan Block Matrices

Lecture Notes 2. The ability to manipulate matrices is critical in economics.

Chapter 4 (Part 1): Non-Parametric Classification (Sections ) Pattern Classification 4.3) Announcements

Module 1 : The equation of continuity. Lecture 5: Conservation of Mass for each species. & Fick s Law

STK4011 and STK9011 Autumn 2016

( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices.

Simple Linear Regression

Fibonacci Identities as Binomial Sums

Third handout: On the Gini Index

CHAPTER 4 RADICAL EXPRESSIONS

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

ESS Line Fitting

Lecture Notes Types of economic variables

ECONOMETRIC THEORY. MODULE VIII Lecture - 26 Heteroskedasticity

Chapter 5 Properties of a Random Sample

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

1 0, x? x x. 1 Root finding. 1.1 Introduction. Solve[x^2-1 0,x] {{x -1},{x 1}} Plot[x^2-1,{x,-2,2}] 3

Summary of the lecture in Biostatistics

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

NATIONAL SENIOR CERTIFICATE GRADE 11

Lecture 07: Poles and Zeros

X X X E[ ] E X E X. is the ()m n where the ( i,)th. j element is the mean of the ( i,)th., then

Department of Agricultural Economics. PhD Qualifier Examination. August 2011

TESTS BASED ON MAXIMUM LIKELIHOOD

The Mathematical Appendix

ECE 595, Section 10 Numerical Simulations Lecture 19: FEM for Electronic Transport. Prof. Peter Bermel February 22, 2013

Rademacher Complexity. Examples

Numerical Differentiation

Multiple Choice Test. Chapter Adequacy of Models for Regression

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

Chapter 8. Inferences about More Than Two Population Central Values

A stopping criterion for Richardson s extrapolation scheme. under finite digit arithmetic.

Multiple Regression. More than 2 variables! Grade on Final. Multiple Regression 11/21/2012. Exam 2 Grades. Exam 2 Re-grades

Basic Concepts in Numerical Analysis November 6, 2017

Assignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class)

Lecture Note to Rice Chapter 8

Statistics Descriptive and Inferential Statistics. Instructor: Daisuke Nagakura

MA/CSSE 473 Day 27. Dynamic programming

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

The equation is sometimes presented in form Y = a + b x. This is reasonable, but it s not the notation we use.

STA302/1001-Fall 2008 Midterm Test October 21, 2008

Chemistry 163B Introduction to Multicomponent Systems and Partial Molar Quantities

CH E 374 Computational Methods in Engineering Fall 2007

hp calculators HP 30S Statistics Averages and Standard Deviations Average and Standard Deviation Practice Finding Averages and Standard Deviations

L5 Polynomial / Spline Curves

The conformations of linear polymers

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

End of Finite Volume Methods Cartesian grids. Solution of the Navier-Stokes Equations. REVIEW Lecture 17: Higher order (interpolation) schemes

( t) ( t) ( t) ρ ψ ψ. (9.1)

Chapter Statistics Background of Regression Analysis

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

Taylor s Series and Interpolation. Interpolation & Curve-fitting. CIS Interpolation. Basic Scenario. Taylor Series interpolates at a specific

Log1 Contest Round 2 Theta Complex Numbers. 4 points each. 5 points each

GMS Equations From Irreversible Thermodynamics

Chapter 5. Curve fitting

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.

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

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

Chemistry 163B Introduction to Multicomponent Systems and Partial Molar Quantities

Functions of Random Variables

Investigation of Partially Conditional RP Model with Response Error. Ed Stanek

Chapter -2 Simple Random Sampling

Example. Row Hydrogen Carbon

Numerical Differentiation

Reliability evaluation of distribution network based on improved non. sequential Monte Carlo method

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

To use adaptive cluster sampling we must first make some definitions of the sampling universe:

Dimensionality Reduction and Learning

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

ECON 482 / WH Hong The Simple Regression Model 1. Definition of the Simple Regression Model

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

M2S1 - EXERCISES 8: SOLUTIONS

Introduction to Matrices and Matrix Approach to Simple Linear Regression

Lecture 3 Probability review (cont d)

ON THE MOTION OF PLANAR BARS SYSTEMS WITH CLEARANCES IN JOINTS

The Selection Problem - Variable Size Decrease/Conquer (Practice with algorithm analysis)

Transcription:

AE 6766/Setzma Sprg 004 Geeral Metod for Calculatg Cemcal Equlbrum Composto For gve tal codtos (e.g., for gve reactats, coose te speces to be cluded te products. As a example, for combusto of ydroge wt ar we mgt cose te followg products:,,,,,,, ad. I terms of a coverso reacto, we would wrte: φ + ( + 376. + + + + + + + were φ represets a arbtrary umber of moles of te fuel (ad also te way ts equato was cose to be wrtte t s also called te equvalece rato, were a value of φ represets just te rgt umber of moles of fuel to react wt te oxyge to form oly te most stable combusto products, ts case. For a mxture of M speces, tere are geerally M+ ukows: M cocetratos,, ad two tesve propertes, e.g., T ad. For a mxture cosstg of R atoms, e.g. R4 for te C/// system, we ca wrte R atom coservato equatos. You ca tk of eac of te coservato equatos as tal codtos or costrats. (ote, tese R equatos may ot be depedet, wc case we ca oly use as may as are depedet. If we furter specfy two termodyamc propertes, we typcally ave M-R ukows. e could specfy te T ad of te products, or for example adabatc flame temperature calculatos, you would specfy ad of te products. We solve for te M-R ukows usg stocometrc reacto relatosps to gve us eoug depedet p. To come up wt te M-R reactos, oe metod s to wrte formato reactos for eac speces preset, except for te elemet speces. (ote: ts metod s ot so elpful f oe of te elemet speces ot part of te mxture. I our example (M8,R3, we eed 5 formato reactos: + + + + 3 ext, wrte equlbrum relatosps for eac formato reactos usg te pf, for eac. For our ydroge/ar example, we ave:

f, f, f, f, f, 4 ow, we clude te (depedet atom coservato equatos. Aga for our example, we get: atoms atoms atoms φ ( ( 3.76 ( were s te al umber of product moles per φ moles of ad s ukow at ts pot. 5 To remove te depedece, we use atom balace ratos (pyscally, t s tese ratos, ot te al umber of moles, wc are most mportat, ad we add te costrat tat te mole fractos must sum to uty, e.g., atoms atoms atoms atoms φ ( 3.76 ( ( ( 6 ow gve two termodyamc propertes, we ave eoug formato to solve for te ukow. If T ad of te products are kow, te soluto smply cossts of determg te pf, from a source suc as te JAAF tables. If te fal temperature s ukow, for example a adabatc flame temperature calculato, te te soluto s teratve: guessg T, fdg te product composto, te calculatg ts assocated T ad usg t to mprove your guess at T. As a alteratve for calculatg te fal temperature, oe ca realze tat most of te eergy s assocated wt te presece of te major (largest speces. Teore you ca gore all te oter/ mor speces o your frst terato ad get a very close estmate of T usg te major speces oly. Te go back ad reterate, ow cludg te mor speces. For roug estmates of product compostos, you ca smply take te temperature ad te speces mole fractos foud from te major speces product calculatos ad use tem, alog wt stocometrc reactos tat form te mor speces from te major speces (.e., approprate p, to calculate te mor speces cocetratos. Ts approac s kow as te major-mor model or major-mor speces approxmato. I some cases, lke te - -

ydroge/ar example descrbed above, oe ca get a smple algebrac soluto for te mole fractos of te major products usg te major speces model (see below. f course, te easest way to solve te problem s to use a cemcal equlbrum computer code/tool. You stll ave to determe te products to be cluded te calculato, ad te tal codtos, e.g., tal atom ratos, but te te computer ca perform te termodyamc property evaluatos ad te teratos!! 7 Major-Mor Model: To llustrate te use of te major-mor model, let s estmate te flame temperature for te ydroge/ar combusto example. Frst, we coose te major speces; we let te products be, ad eter (for lea mxtures or (for rc mxtures. Wrtg te reactos for te two cases ad deotg te stocometrc coeffcets for te products terms of φ from smple atom balaces, we ave φ for φ< : φ + ( + 376. φ + + 88. for φ> : φ + ( + 376. + ( φ + 88. Major Speces Mole Fractos: From te above reacto equatos ad wt algebra, we get: φ< (lea φ> (rc φ φ+ + 88 φ+ 88. + 88. φ. ( φ+ 88. φ + + 88. 88. ( φ+ 88. ( φ ( φ + 88. Tus smply gve φ (te : rato, we kow te product composto ad ca calculate te fal temperature. For example wt φ.3 (rc combusto, we get, 0.345 0.59 0.0943 Adabatc Flame Temperature: Assumg a tal temperature of -55 C (8 for te reactats, we fd te adabatc flame temperature usg 0 ad R0,.e., roducts ( T ( ad Reactats [( + ] [( + ] Tad T roducts 8 Reactats 8 T - 3 -

mol Wrtg out te summato for eac product ad reactat, ad usg te umber of moles of eac for our φ.3 flame we ave,.3mol [( T + ] +.88mol [( ] [( ] T + + 0.3mol T + ad, ad ad f T [( ] [( ] [( ] 8 + + 0.5mol 8 + +.88mol 8 + Sce te etalpy of formato for elemets s zero, we get mol.3 [( T + ] +.88mol [( ] [( ] T + 0.3mol T ad ad ad T moles [( ] [( ] [( ] 8 + 0.5moles 8 +.88moles 8 Data for f of water ad te sesble etalpy cages for eac speces * ca be foud a umber of sources, e.g., te JAAF tables. Usg ts data, oe fds Tad 300. Mor Speces Mole Fractos: We ca wrte te followg stocometrc relatosps betwee te mor speces ad te major speces of our rc ydroge/ar flame (for empass, te mor speces are wrtte bold letters. Eac reacto represets a metod for producg te mor speces usg oly te major products of our ydroge/ar flame. + + + + + + ow we ca wrte te followg expressos for te mole fractos of te mor speces terms of te major speces : f, f, f, f, f, f, f, f, * T It s ot reasoable to assume tat c p s costat c ( T T T p dt for our large temperature rage. - 4 - T

were we ave used te fact tat p for eac of te stocometrc reactos s smply a fucto of te formato equlbrum costats of te speces te reacto. For example for te reacto + te equlbrum costat s gve by p f, f, f, f, sce te formato equlbrum costat pf of a elemet s uty (by defto. Usg te p f from te JAAF tables at 300, te estmated for te major speces, ad assumg a pressure of bar, we get te values for te mor speces lsted te table below. As a comparso, te table below also cludes results from a complete soluto obtaed wt te STAJA cemcal equlbrum code. Wle ot completely accurate, te major-mor model does a good job of predctg te flame temperature (+5 or ~% relatve error ad te major speces mole fractos (<% relatve error, ad gves a reasoable estmate of te mor speces mole fractos (wt a 40% relatve error, ad certaly better ta oe order-of-magtude. Tus for roug approxmatos, te major-mor model s smple ad reasoably accurate. Speces/Tad Major-Mor Model Full Calculato Error(% Tad 300 75. 59.% 59.0% 0.7 3.5% 3.% 0.96 9.4% 9.4% 0 0.30% (3000 ppm 0.6% (600 ppm 5 0.6% (600 ppm 0.3% (300 ppm 3 0 ppm 80 ppm 49 ppm 36 ppm 39 3.4 ppb.5 ppb 36 ppmparts per mllo ( 0-6, ppbparts per bllo ( 0-9 - 5 -