McCabe-Thiele Diagrams for Binary Distillation
|
|
- Hannah McDonald
- 6 years ago
- Views:
Transcription
1 McCabe-Thele Dagrams for Bnary Dstllaton Tore Haug-Warberg Dept. of Chemcal Engneerng August 31st, 2005 F V 1 V 2 L 1 V n L n 1 V n+1 L n V N L N 1 L N L 0 VN+1 Q < 0 D Q > 0 B FIGURE 1: Smplfed pcture of a counter-current dstllaton column. The McCabe Thele dagram for solvng bnary dstllaton problems represents a pllarstone n every chemcal engneerng class. The reason for ths s not ts practcal mportance whch n the era of modern computer technology has become very modest but ts educatonal value based on a very llustratve graphcal pcture of a complcated physcal process. The McCabe Thele dagram n ts most smple form represents a graphcal soluton to the classcal separaton problem of N + 1 deal equlbrum stages connected wth countercurrent vapor and lqud flows as llustrated n Fg.1. The streams are shown as separated flows except at certan dstnct locatons (trays) where complete mxng and thermodynamc equlbrum s assumed. Ths s a useful abstracton but t must be understood that the stuaton n a real plant devates n two respects: Frst of all there s no one-to-one correspondance between a physcal tray (or plate) and a thermodynamc equlbrum stage 1. Secondly, many columns have ther plates replaced by a structured packng whch gves ntmate contact between the vapor and lqud phases at all locatons n the column. The smplfed pcture s stll vald, however, f each tray s consdered to be one theoretcal equlbrum stage regardless of the physcal realzaton of t. 1 It s common practce to talk about tray effcences whch s the fracton of a theoretcal equlbrum stage reached for each physcal plate. Ths effcency may reach values as hgh as but s n many cases much lower. 1
2 2 1. The energy balance To calculate the soluton of the complcated separaton problem n Fg.1 we need 3 sets of mass balances, one for the rectfyng, one for the feed and one for the strppng secton of the column. In addton there s an overall energy conservaton assumpton known as constant molar overflow. The orgn and mportance of ths assumpton s clearfed below. Let n represent any of the nternal stages n the the column except the feed tray. The energy balance for an adabatc stage s Ḣ vap +1 + Ḣlq 1 = Ḣvap + Ḣ lq Snce enthalpy s an Euler homogenous functon of frst order (extensve state varable) the energy balance can also be wrtten h vap +1 V +1 + h lq 1 L 1 = h vap V + h lq L where h, andv and L denote the molar enthalpy and molar flows of the vapor and lqud streams respectvely 2. Combned wth the mass balance V +1 + L 1 = V + L for stage the last equaton becomes ( h lq 1 ) hvap +1 L 1 = ( ) h vap +1 V + ( h lq h+1) vap L Addng and subtractng h vap 1 L 1 and h vap L brngs the equaton on the followng form: ( h lq 1 hvap 1 ) L 1 + ( h vap 1 ) hvap +1 L 1 = ( h vap + ( h lq +1) V ) L + ( h vap +1) L The enthalpy functon has no absolute zero, and t s only the enthalpy dfferences that matter n the energy balance. But, the mportance of the temperature dependency of enthalpy s usually much smaller than the heat of vaporzaton 3. The assumpton that h vap 1 = hvap = h vap +1 s therefore qute legtmate and the energy balance can be smplfed to: ( h lq 1 ) hvap 1 L 1 = ( ) h lq L Fnally, assumng that the heat of vaporzaton s constant for the two-phase mxture makes h lq 1 hvap 1 = hlq and consequently L 1 = L.Thelast argument can easly be extended to yeld L 0 = L 1 = L 2 =...= L n 1 V 1 = V 2 = V 3 =...= V n 2 The dot-notaton has been dropped for the rate symbols. Ths s accordng to common practce n ths feld. 3 Typcal changes per stage are J/mol for sensble heats and J/mol for heats of vaporzaton
3 2. OPERATING LINES 3 whch s known as the constant molar overflow assumpton 4. The same assumpton s vald for the bottom secton as well, but not across the feed tray, where a dscontnuty occurs because of the feed stream enterng from outsde the control volume. 2. Operatng lnes Ths secton (and the rest of the paper) s restrcted to the dstllaton of bnary mxtures. In the top secton of the column an enrchment of the lghter of the two components takes place and n the bottom secton the heaver component s enrched 5. The next task s to express the mass balances of the column and have them ntegrated wth the energy balance from the prevous secton nto somethng called an operatng lne. The composton of the lghter component s desgnated by y and x for the vapor and lqud phases respectvely. By conventon these symbols stand for the lghter component The rectfyng secton. We shall wrte the mass balances as one total balance and one component balance. For the top secton these equatons are V = D + L 1 y V = x 0 D + x 1 L 1 whch means the heaver component s never consdered explctly 6. The equatons are easly solved for the vapor composton: D L 1 y = x 0 + x 1 V V D L 1 = x 0 + x 1 D + L 1 D + L 1 The constant molar overflow assumpton yelds L 0 = L 1 =... = L n 1 L. Furthermore, t s customary to defne the external reflux rato of the column as R L/D [0, : 1 y = x R + x R 1, [1, n] 1 + R Ths s the so-called operatng lne representng the combned mass and energy balances n the rectfyng secton of the column. The equaton s vald for all of the nternal stages = 1, 2,...,n and n partcular for = 1 whch yelds y 1 = x 0 = x D. A bt surprsng maybe, 4 Vald for any number of the chemcal components n the mxture. 5 Ths means that the temperature profle of the column decreases from the bottom towards the top. 6 For bnary mxtures t can easly be elmnated from the gven equatons.
4 4 y 1 y 2 y 3 y 4 slope = R 1+R x 0 1+R x 3 x 2 x 1 x D FIGURE 2: Rectfyng secton of a benzene toluene column (R = 2). but from Fg.1 we see that the condenser has no separatng effect t only condenses the vapor nto a lqud phase 7,8. For = 2 the operatng lne connects the vapor composton from stage 2 wth the lqud composton from stage 1, but here the calculaton stops. We cannot move on before the lqud composton from stage 1 s known! The answer to ths problem les n the assumpton of thermodynamc equlbrum from whch we can derve x 1 = x eq (y 1 ), x 2 = x eq (y 2 ) and so on. Wth ths relatonshp n mnd we got all three ngredents of the McCabe Thele graphcal method: (1) Constant molar overflow assumpton (2) Operatng lne (3) Equlbrum curve The dagram s conventonally drawn wth the heaver component to the left and the lghter component to the rght. That means the equlbrum curve s located above the man 45 dagonal n the dagram. Fg.2 shows a typcal example for an almost deal mxture of benzene (lght) and toluene (heavy). 7 Ths s referred to as a total condenser. In some cases the column s equpped wth a partal condenser whch makes x D x 0. The condenser must then be treated lke an addtonal equlbrum stage, see the dscusson of the strppng secton below. 8 Prove that the heat duty of the condenser s n the range R Q/(D vap h) R + 1for partal versus total operaton.
5 2. OPERATING LINES 5 The starcase (red) shows the nteracton between the equlbrum curve (vertcal step) and the operatng lne (horsontal step). Note that the operatng lne (black) crosses the equlbrum lne (blue) at one pont n the dagram. The column cannot operate to the left of ths cross-over pont because that would be n the one-phase regon of the phase dagram. Hence, the number of theoretcal stages approaches nfnty as the column strves towards a sngular condton The strppng secton. The operatng lne for the bottom secton follows the same recepe as for the top secton. The two mass balances are prncpally the same, V + B = L 1 y V + x B B = x 1 L 1 and agan we solve for the vapor composton on the left hand sde: B L 1 y = x B + x 1 V V B V + B = x B + x 1 V V The constant molar overflow assumpton yelds V n+1 = V n+2 =...= V N+1 V. Furthermore, we may characterze the heat duty 10 of the reboler as S V/(V + B) 0, 1] and wrte 1 S 1 y = x B + x 1, [n + 1, N + 1] S S where S s the fracton of the downcomng lqud whch s beng evaporated. Alternatvely we can defne S as the heat duty of the reboler dvded by the heat whch s needed to evaporate all the downcomng lqud. The operatng lne s vald for all the nternal stages = n + 1, n + 2,...,N and also for the reboler = N + 1. The strppng secton has therefore one external equlbrum stage n addton to the nternal stages of the column (the stuaton s smlar to that of a partal condenser n the rectfyng secton). The cross-over between the operatng lne and the man dagonal s gven by y N+2 = x N+1 = x B. The ndex = N + 2 s strctly speakng outsde the control volume but t does stll serve as an easy-toremember constructon pont of the dagram, see Fg.3. At ths pont t s possble to count the number of theoretcal stages from the specfcatons of the top and bottom products together wth the condenser and reboler dutes. It s more common, however, to specfy 9 Ths s barely vsble n the current dagram but an nfntly large magnfyng glass would reveal an nfnte number of successvely smaller steps. 10 Prove that the heat duty of the reboler s Q = S (1 S ) 1 B vap h.
6 6 y 5 y 6 slope = 1 S y 7 y 8 x B x 7 x 6 x 5 x 4 FIGURE 3: Strppng secton of a benzene toluene column (S = 2/3). Tray numbers are rounded to the nearest nteger. the feed stream rather than the reboler duty. To manage ths we need the operatng lne for the feed tray whch s dscussed n the next secton The feed tray. The operatng lne for the feed tray (q-lne) requres four mass balances and two energy balances (ncludng the constant molar overflow assumpton). The extra mass and energy balances are needed to characterze the physcal state of the feed stream. To begn wth we have the component balances y n V n = x D D + x n 1 L n 1 x n L n = y n+1 V n+1 + x B B for the rectfyng and strppng sectons respectvely. At the feed tray there wll n general be sudden jumps n all the stream varables, but for an optmal placement of the feed tray the ntensve state varables (temperature and composton) are contnuous. Hence, we shall requre that y n = y n+1 y and that x n = x n 1 x n the equatons above. Summaton plus reorganzaton yelds: y (V n+1 V n ) = x (L n L n 1 ) (x D D + x B B) } {{ } x F F
7 2. OPERATING LINES 7 Note that the total mass balance x F F = x D D+ x B B s substtuted on the rght sde. The q-lne s frst wrtten on the prelmnary form: y = x L n L n 1 F x F V n+1 V n V n+1 V n The feed stream s subsequently replaced by the mass balance for the feed tray F = (L n L n 1 ) (V n+1 V n ). A lttle manpulaton yelds y = x L ( ) n L n 1 Ln L n 1 V n+1 V n x F x L V x F ( L V 1 1 V n+1 V n ) where V and L measure the dfferences n the vapor and lqud flows across the feed tray (top-down). Ths could serve as our fnal result, but t s customary to characterze the feed qualty by L q L V where q = 1 V = 0 (saturated lqud feed) and q = 0 L = 0 (saturated vapor feed). Hence, the feed qualty can be nterpreted as the fracton of saturated lqud n the feed. The defnton of the feed qualty s easly nverted to gve L/ V = q/(q 1) whch puts the q-lne on the fnal form: q y = x q 1 x 1 F q 1 Ths s the last of the three operatng lnes used n a McCabe-Thele dagram, but t should be mentoned that the optmal feed tray locaton causes a slght redundancy n the calculaton of the dagram: E.g. n most cases the composton specfcatons x F and x D are gven n addton to the reflux rato R and the feed qualty q. Ths specfes the upper part of the column and the condtons of the strpper wll be fxed once x B or S s known (the other parameter s fxed by the common cross-over pont of the operatng lnes). Anyway, the calculaton s complete n the sense that t reveales all the sgnfcant fgures of a bnary dstllaton column (see also Fg.4): (1) Compostons of the top and bottom products (2) Composton and qualty of the feed stream (3) Reflux rato (and condenser duty) (4) Reboler duty 2.4. Extreme specfcatons. The condtonof mnmum reflux s reached then ) the cross-over between the operatng lnes s located at the phase boundary as llustrated n Fg.2 or ) the upper operatng lne s tangent to the equlbrum curve whle the cross-over pont s nsde the phase dagram
8 slope = q q 1 N top = N bot = FIGURE 4: Complete McCabe-Thele dagram for a benzene toluene column wth parameters x D = 0.9, x F = 0.55, x B = 0.1, R = 2and q = 0.2 (consstent wth S = 2/3). (apples to mxtures wth concave equlbrum curves). In both cases the number of theoretcal stages approaches nfnty. The other extreme s then the reflux s nfnte and the operatng lnes for the column overlap wth the man dagonal. Ths yelds the smallest number of theoretcal stages for a gven separaton. There are also cases when the strppng secton or the rectfyng secton s mssng or neffectve. These cases nclude batch dstllaton (no strpper) and degassng unts (no rectfyer). Mathematcally, there s also the possblty that a secton becomes neffectve because of an extreme feed condton. Ths s readly verfed from the q-lne whch has slope ±1 for q = ±. The frst case (q = ) corresponds to a lqud feed so cold that t condenses all the vapor comng up from the reboler (neffectve rectfyng secton). The second case (q = ) corresponds to a superheated vapor so hot that t evaporates all the lqud comng down from the condenser (neffectve strppng secton).
Supplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationSolution Thermodynamics
Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs
More informationy i x P vap 10 A T SOLUTION TO HOMEWORK #7 #Problem
SOLUTION TO HOMEWORK #7 #roblem 1 10.1-1 a. In order to solve ths problem, we need to know what happens at the bubble pont; at ths pont, the frst bubble s formed, so we can assume that all of the number
More informationExercises of Fundamentals of Chemical Processes
Department of Energ Poltecnco d Mlano a Lambruschn 4 2056 MILANO Exercses of undamentals of Chemcal Processes Prof. Ganpero Gropp Exercse 7 ) Estmaton of the composton of the streams at the ext of an sothermal
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationName: SID: Discussion Session:
Name: SID: Dscusson Sesson: Chemcal Engneerng Thermodynamcs 141 -- Fall 007 Thursday, November 15, 007 Mdterm II SOLUTIONS - 70 mnutes 110 Ponts Total Closed Book and Notes (0 ponts) 1. Evaluate whether
More informationCinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure
nhe roblem-solvng Strategy hapter 4 Transformaton rocess onceptual Model formulaton procedure Mathematcal Model The mathematcal model s an abstracton that represents the engneerng phenomena occurrng n
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationThe ChemSep Book. Harry A. Kooijman Consultant. Ross Taylor Clarkson University, Potsdam, New York University of Twente, Enschede, The Netherlands
The ChemSep Book Harry A. Koojman Consultant Ross Taylor Clarkson Unversty, Potsdam, New York Unversty of Twente, Enschede, The Netherlands Lbr Books on Demand www.bod.de Copyrght c 2000 by H.A. Koojman
More informationTHE SUMMATION NOTATION Ʃ
Sngle Subscrpt otaton THE SUMMATIO OTATIO Ʃ Most of the calculatons we perform n statstcs are repettve operatons on lsts of numbers. For example, we compute the sum of a set of numbers, or the sum of the
More information( ) 1/ 2. ( P SO2 )( P O2 ) 1/ 2.
Chemstry 360 Dr. Jean M. Standard Problem Set 9 Solutons. The followng chemcal reacton converts sulfur doxde to sulfur troxde. SO ( g) + O ( g) SO 3 ( l). (a.) Wrte the expresson for K eq for ths reacton.
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More information10.34 Numerical Methods Applied to Chemical Engineering Fall Homework #3: Systems of Nonlinear Equations and Optimization
10.34 Numercal Methods Appled to Chemcal Engneerng Fall 2015 Homework #3: Systems of Nonlnear Equatons and Optmzaton Problem 1 (30 ponts). A (homogeneous) azeotrope s a composton of a multcomponent mxture
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationEstimation of the composition of the liquid and vapor streams exiting a flash unit with a supercritical component
Department of Energ oltecnco d Mlano Va Lambruschn - 05 MILANO Eercses of Fundamentals of Chemcal rocesses rof. Ganpero Gropp Eercse 8 Estmaton of the composton of the lqud and vapor streams etng a unt
More informationModule 9. Lecture 6. Duality in Assignment Problems
Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationAdiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram
Adabatc Sorpton of Ammona-Water System and Depctng n p-t-x Dagram J. POSPISIL, Z. SKALA Faculty of Mechancal Engneerng Brno Unversty of Technology Techncka 2, Brno 61669 CZECH REPUBLIC Abstract: - Absorpton
More informationIf two volatile and miscible liquids are combined to form a solution, Raoult s law is not obeyed. Use the experimental data in Table 9.
9.9 Real Solutons Exhbt Devatons from Raoult s Law If two volatle and mscble lquds are combned to form a soluton, Raoult s law s not obeyed. Use the expermental data n Table 9.3: Physcal Chemstry 00 Pearson
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
More information9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 - Chapter 9R -Davd Klenfeld - Fall 2005 9 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys a set
More informationCHEMICAL ENGINEERING
Postal Correspondence GATE & PSUs -MT To Buy Postal Correspondence Packages call at 0-9990657855 1 TABLE OF CONTENT S. No. Ttle Page no. 1. Introducton 3 2. Dffuson 10 3. Dryng and Humdfcaton 24 4. Absorpton
More informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
More information(1) The saturation vapor pressure as a function of temperature, often given by the Antoine equation:
CE304, Sprng 2004 Lecture 22 Lecture 22: Topcs n Phase Equlbra, part : For the remander of the course, we wll return to the subject of vapor/lqud equlbrum and ntroduce other phase equlbrum calculatons
More information1 GSW Iterative Techniques for y = Ax
1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn
More informationAffine transformations and convexity
Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/
More informationEnergy, Entropy, and Availability Balances Phase Equilibria. Nonideal Thermodynamic Property Models. Selecting an Appropriate Model
Lecture 4. Thermodynamcs [Ch. 2] Energy, Entropy, and Avalablty Balances Phase Equlbra - Fugactes and actvty coeffcents -K-values Nondeal Thermodynamc Property Models - P-v-T equaton-of-state models -
More informationANALYSIS OF THE ZONE CONNECTING CONSECUTIVE SECTORS IN GENERALIZED DISTILLATION COLUMNS BY USING THE PONCHON-SAVARIT METHOD
Latn Amercan Appled Research 46:127-132 (2016) ANALYSIS OF THE ZONE ONNETING ONSEUTIVE SETORS IN GENERALIZED DISTILLATION OLUMNS BY USING THE PONHON-SAVARIT METHOD A. MARILLA, M.D. SERRANO and J.A. REYES-LABARTA
More informationGeneral Thermodynamics for Process Simulation. Dr. Jungho Cho, Professor Department of Chemical Engineering Dong Yang University
General Thermodynamcs for Process Smulaton Dr. Jungho Cho, Professor Department of Chemcal Engneerng Dong Yang Unversty Four Crtera for Equlbra μ = μ v Stuaton α T = T β α β P = P l μ = μ l1 l 2 Thermal
More informationNon-Ideality Through Fugacity and Activity
Non-Idealty Through Fugacty and Actvty S. Patel Deartment of Chemstry and Bochemstry, Unversty of Delaware, Newark, Delaware 19716, USA Corresondng author. E-mal: saatel@udel.edu 1 I. FUGACITY In ths dscusson,
More informationA Self-Consistent Gibbs Excess Mixing Rule for Cubic Equations of State: derivation and fugacity coefficients
A Self-Consstent Gbbs Excess Mxng Rule for Cubc Equatons of State: dervaton and fugacty coeffcents Paula B. Staudt, Rafael de P. Soares Departamento de Engenhara Químca, Escola de Engenhara, Unversdade
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationElectrical double layer: revisit based on boundary conditions
Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer
More informationModule 3: The Whole-Process Perspective for Thermochemical Hydrogen
"Thermodynamc Analyss of Processes for Hydrogen Generaton by Decomposton of Water" by John P. O'Connell Department of Chemcal Engneerng Unversty of Vrgna Charlottesvlle, VA 2294-4741 A Set of Energy Educaton
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationThe optimal delay of the second test is therefore approximately 210 hours earlier than =2.
THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple
More informationEconomics 101. Lecture 4 - Equilibrium and Efficiency
Economcs 0 Lecture 4 - Equlbrum and Effcency Intro As dscussed n the prevous lecture, we wll now move from an envronment where we looed at consumers mang decsons n solaton to analyzng economes full of
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationLecture 3. Ax x i a i. i i
18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More information...Thermodynamics. If Clausius Clapeyron fails. l T (v 2 v 1 ) = 0/0 Second order phase transition ( S, v = 0)
If Clausus Clapeyron fals ( ) dp dt pb =...Thermodynamcs l T (v 2 v 1 ) = 0/0 Second order phase transton ( S, v = 0) ( ) dp = c P,1 c P,2 dt Tv(β 1 β 2 ) Two phases ntermngled Ferromagnet (Excess spn-up
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationSome modelling aspects for the Matlab implementation of MMA
Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationAPPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS
Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent
More informationChapter 8. Potential Energy and Conservation of Energy
Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal
More informationThe Geometry of Logit and Probit
The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.
More informationWinter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan
Wnter 2008 CS567 Stochastc Lnear/Integer Programmng Guest Lecturer: Xu, Huan Class 2: More Modelng Examples 1 Capacty Expanson Capacty expanson models optmal choces of the tmng and levels of nvestments
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationFoundations of Arithmetic
Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an
More informationAssignment 4. Adsorption Isotherms
Insttute of Process Engneerng Assgnment 4. Adsorpton Isotherms Part A: Compettve adsorpton of methane and ethane In large scale adsorpton processes, more than one compound from a mxture of gases get adsorbed,
More informationHomework Chapter 21 Solutions!!
Homework Chapter 1 Solutons 1.7 1.13 1.17 1.19 1.6 1.33 1.45 1.51 1.71 page 1 Problem 1.7 A mole sample of oxygen gas s confned to a 5 lter vessel at a pressure of 8 atm. Fnd the average translatonal knetc
More information1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More information3. Be able to derive the chemical equilibrium constants from statistical mechanics.
Lecture #17 1 Lecture 17 Objectves: 1. Notaton of chemcal reactons 2. General equlbrum 3. Be able to derve the chemcal equlbrum constants from statstcal mechancs. 4. Identfy how nondeal behavor can be
More informationThermodynamics II. Department of Chemical Engineering. Prof. Kim, Jong Hak
Thermodynamcs II Department of Chemcal Engneerng Prof. Km, Jong Hak Soluton Thermodynamcs : theory Obectve : lay the theoretcal foundaton for applcatons of thermodynamcs to gas mxture and lqud soluton
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationTemperature. Chapter Heat Engine
Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the
More informationThe Expectation-Maximization Algorithm
The Expectaton-Maxmaton Algorthm Charles Elan elan@cs.ucsd.edu November 16, 2007 Ths chapter explans the EM algorthm at multple levels of generalty. Secton 1 gves the standard hgh-level verson of the algorthm.
More informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationDummy variables in multiple variable regression model
WESS Econometrcs (Handout ) Dummy varables n multple varable regresson model. Addtve dummy varables In the prevous handout we consdered the followng regresson model: y x 2x2 k xk,, 2,, n and we nterpreted
More informationOutline. Unit Eight Calculations with Entropy. The Second Law. Second Law Notes. Uses of Entropy. Entropy is a Property.
Unt Eght Calculatons wth Entropy Mechancal Engneerng 370 Thermodynamcs Larry Caretto October 6, 010 Outlne Quz Seven Solutons Second law revew Goals for unt eght Usng entropy to calculate the maxmum work
More informationMeasurement Uncertainties Reference
Measurement Uncertantes Reference Introducton We all ntutvely now that no epermental measurement can be perfect. It s possble to mae ths dea quanttatve. It can be stated ths way: the result of an ndvdual
More information8.6 The Complex Number System
8.6 The Complex Number System Earler n the chapter, we mentoned that we cannot have a negatve under a square root, snce the square of any postve or negatve number s always postve. In ths secton we want
More informationAssortment Optimization under MNL
Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationProf. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model
EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationDynamic Model of Dividing Wall Column for Separation of Ternary System
R. K. Dohare, K. Sngh, R. Kumar, S. Upadhyaya, S. Gupta Dynamc Model of Dvdng Wall Column for Separaton of Ternary System Raeev Kumar Dohare a, Kalash Sngh a, Raesh Kumar a, Sushant Upadhyaya a, Sourabh
More information36.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
ChE Lecture Notes - D. Keer, 5/9/98 Lecture 6,7,8 - Rootndng n systems o equatons (A) Theory (B) Problems (C) MATLAB Applcatons Tet: Supplementary notes rom Instructor 6. Why s t mportant to be able to
More informationCHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR OF SEPARATION PRODUCTS OUTPUT DETERMINATION
Górnctwo Geonżynera Rok 0 Zeszyt / 006 Igor Konstantnovch Mladetskj * Petr Ivanovch Plov * Ekaterna Nkolaevna Kobets * Tasya Igorevna Markova * CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR
More informationNAME and Section No. it is found that 0.6 mol of O
NAME and Secton No. Chemstry 391 Fall 7 Exam III KEY 1. (3 Ponts) ***Do 5 out of 6***(If 6 are done only the frst 5 wll be graded)*** a). In the reacton 3O O3 t s found that.6 mol of O are consumed. Fnd
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationElectrochemical Equilibrium Electromotive Force
CHM465/865, 24-3, Lecture 5-7, 2 th Sep., 24 lectrochemcal qulbrum lectromotve Force Relaton between chemcal and electrc drvng forces lectrochemcal system at constant T and p: consder Gbbs free energy
More informationIntroduction to information theory and data compression
Introducton to nformaton theory and data compresson Adel Magra, Emma Gouné, Irène Woo March 8, 207 Ths s the augmented transcrpt of a lecture gven by Luc Devroye on March 9th 207 for a Data Structures
More information763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.
7636S ADVANCED QUANTUM MECHANICS Soluton Set 1 Sprng 013 1 Warm-up Show that the egenvalues of a Hermtan operator  are real and that the egenkets correspondng to dfferent egenvalues are orthogonal (b)
More informationDigital Signal Processing
Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationLecture. Polymer Thermodynamics 0331 L Chemical Potential
Prof. Dr. rer. nat. habl. S. Enders Faculty III for Process Scence Insttute of Chemcal Engneerng Department of Thermodynamcs Lecture Polymer Thermodynamcs 033 L 337 3. Chemcal Potental Polymer Thermodynamcs
More informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More informationSTATISTICAL MECHANICS
STATISTICAL MECHANICS Thermal Energy Recall that KE can always be separated nto 2 terms: KE system = 1 2 M 2 total v CM KE nternal Rgd-body rotaton and elastc / sound waves Use smplfyng assumptons KE of
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationMath1110 (Spring 2009) Prelim 3 - Solutions
Math 1110 (Sprng 2009) Solutons to Prelm 3 (04/21/2009) 1 Queston 1. (16 ponts) Short answer. Math1110 (Sprng 2009) Prelm 3 - Solutons x a 1 (a) (4 ponts) Please evaluate lm, where a and b are postve numbers.
More informationHMMT February 2016 February 20, 2016
HMMT February 016 February 0, 016 Combnatorcs 1. For postve ntegers n, let S n be the set of ntegers x such that n dstnct lnes, no three concurrent, can dvde a plane nto x regons (for example, S = {3,
More information#64. ΔS for Isothermal Mixing of Ideal Gases
#64 Carnot Heat Engne ΔS for Isothermal Mxng of Ideal Gases ds = S dt + S T V V S = P V T T V PV = nrt, P T ds = v T = nr V dv V nr V V = nrln V V = - nrln V V ΔS ΔS ΔS for Isothermal Mxng for Ideal Gases
More informationFormulas for the Determinant
page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use
More informationLecture 7: Boltzmann distribution & Thermodynamics of mixing
Prof. Tbbtt Lecture 7 etworks & Gels Lecture 7: Boltzmann dstrbuton & Thermodynamcs of mxng 1 Suggested readng Prof. Mark W. Tbbtt ETH Zürch 13 März 018 Molecular Drvng Forces Dll and Bromberg: Chapters
More information