5.68J/10.652J Feb 13, 2003 Numerical Solution to Kinetic Equations Lecture 1
|
|
- Harriet Lester
- 5 years ago
- Views:
Transcription
1 5.68J/10.652J Feb 13, 2003 Numercal Soluton to Knetc Equatons Lecture 1 The Problem We want to solve 1 dn l H = ν r ( T l V dt n, ) Transport of ' ' out of V Eq.(1) Where n s the number of moles of speces, ν s the matrx of stochometrc coeffcents, and = Rateof Step 'l' (unts: moles/volume-second) For a typcal bmolecular elementary step reacton A+B = C+D: = k forward ( T )[ A][ B ] k ( T )[ C][ D] For a typcal gas-phase dssocaton reacton A= C + D: = k forward ( T, P )[ A ] k ( T, P )[ C][ D] For a typcal gas-phase somerzaton reacton A = C: = k forward ( T, P )[ A ] k ( T, P )[ C] The rate constants also usually depend (weakly) on the chemcal envronment (e.g. delectrc constant), so they strctly are not really constant but (very weak) functons of all the n s. Note that the forward and rates are normally related by the equlbrum constant, related to the free energy of reacton (thermochemstry). It s mportant to apprecate that except n the smplest cases, there s no analytcal soluton to ths system of dfferental equatons, so we MUST use the computer to solve them numercally. Approxmatons lke quas-steady-state extend the range of systems where analytcal (approxmate) solutons are possble, but you don t have to get very complcated before even the steady-state-approxmaton solutons are ntractable.
2 More About the Fundamental Equaton Chemstry textbooks usually gve nstead of Eq. (1) the equaton: d [ ] l H = ν r ( l, dt n T ) Eq.(2) But Eq.(2) s true only f dv/dt = 0 ( sochorc ) and there s no transport (.e. perfectly homogeneous, no convectve flow). These textbook approxmatons are good for dlute reactons n soluton or for constant volume apparatus (e.g. bomb calormeter, shock tube). In many mportant cases the system s sobarc,.e. dp/dt = 0, so f the temperature ncreases or there s a change n the number of moles of gas V wll change. It s convenent to express the equaton n terms of ntensve quanttes rather than the extensve varables n and V. The ntensve varables that work best are mass fractons yι: n W = [] W y = M total ρ dy W 1dn = Transport of ' ' out of dt ρ V dt V You can plug n the same expresson as above for the term n parentheses. Note that mole fractons x are not such a good choce for knetcs as mass fractons, snce the number of moles can change n chemcal reactons, but the mass does not. However one wrtes the chemcal reacton equatons, one must be explct about the model for dt/dt, snce the knetcs nvarably depends strongly on T. The common textbook assumpton s sothermal (dt/dt = 0) whch s okay f the reacton s very slow and the reagents dlute. In many cases of nterest the system s closer to beng adabatc, so dt/dt s gven by energy conservaton. (For the energy conservaton equaton around a small control volume wth unform T, P, and composton, and an elementary dscusson of how ODE s are solved numercally, see: ). Note that to accurately compute dt/dt one must have a way of computng the enthalpes and heat capactes of all the speces at any T. Thermochemcal Data Because one needs the thermochemcal data both to compute dt/dt and also to compute the reacton rate constants, t s usual to buld up a lbrary fle contanng all ths data for all the speces of nterest. For (stupd) hstorcal reasons the data are usually represented as the coeffcents n a pecewse fttng form known as the NASA polynomals. For some of the arcane detals about NASA polynomals, see the CHEMKIN manual and Beware that there s more than one format for the NASA polynomals! The NIST Webbook gves thermochemcal data n a dfferent Shomate format whch s better
3 behaved than the NASA form. The other popular form s the dscrete values of Benson (used also by THERM and the NIST S&P program) dscussed later n the course. All forms take advantage of the fact that the T-dependence of the enthalpy and the entropy s easly computed (often analytcally) from the heat capacty. So the man ssues are whether one knows H and S exactly at any T for the speces of nterest and then how well the T-dependence of the heat capacty (usually Cp, measured at constant pressure) can be determned and represented. Very often S and Cp are calculated usng statstcal mechancs rather than measured, you need to read carefully to determne whether the numbers reported n a table are measured or just computed, and f they are computed what assumptons were made. The reported S and H are usually for an deal gas state at 1 atm at 298 K, even f the molecules are actually lquds or solds at 298 K and 1 atm. So the reported H(298) s often calculated, perhaps from an expermentally measured H at some dfferent temperature and the computed heat capactes. In many real-world applcatons of knetcs, the estmated barrer heghts of the slow steps are strongly correlated wth the estmates of ther endothermctes, and many of the fast reactons are nearly equlbrated. Hence one s more senstve to uncertantes n the thermochemcal data than to the (usually much larger) uncertantes n the ndvdual rate constants. So a knetcst MUST be concerned about the qualty of the thermochemcal data fle employed n the smulatons. There s plenty to worry about. For example, n 2001 t was dscovered that the heat of formaton of OH was computed ncorrectly many years ago, and that therefore the many speces whose thermochemstry was nferred relatve to OH were also ncorrect. Presumably many smlar errors exst n the databases, watng for some brght graduate student to dentfy them. Note that you wll get observably dfferent results f you use the default therm.dat fle shpped wth CHEMKIN or the GRI-MECH thermo fle. Beware: Inhomogenety can develop Spontaneously If dt/dt s nonzero, often gradents wll develop n the system even f t was ntally homogeneous due to heat losses at the boundares. In some cases these temperature gradents can become extreme (e.g. n a flame). Pressure equlbrates at about the speed of sound, much faster than speces concentraton or temperature gradents can equlbrate, so very often t s a good approxmaton that P s unform n a system. However, n some combuston systems the reacton s so fast that sgnfcant pressure gradents develop; ths can cause loud noses or n an extreme cases detonatons and other types of explosons. These nhomogenetes are reduced by usng a smaller control volume; for reactng flow smulatons people often use control volumes less than a mllmeter across so that ther approxmaton that the system has a sngle temperature and composton s accurate. Requred Inputs What nputs do we need for a Knetc Smulaton? 1) A Thermo Fle (wth data for all the speces n the smulaton) 2) Reacton Mechansm: A lst of the reactons and forward rate constants 3) Reacton condtons (e.g. ntal concentratons, T, P, tmescale, adabatc or sothermal, flow stuaton, etc.) to fully specfy physcal stuaton.
4 4) Smulaton detals: Specfy outputs requred, tolerances, numercal methods, ntal guesses, etc. In CHEMKIN, nput 1 s called somethng lke therm.dat and usually resdes n the /data/ drectory, nput 2 s called somethng lke chem.np (and f there s surface chemstry an addtonal fle surf.np), and nputs 3 and 4 are combned n a sngle fle, e.g. aurora.np. How are ODE s solved numercally? (the short short verson) If we defne Y={y, T} and there s no transport, then Eq.(1) s of the form: dy/dt = F(Y) and we usually know the ntal condtons: Y(t o )=Y o The general procedure s to step forward n tme wth some formula Y(t+ t) = Y(t) + G(Y) t Eq.(3) where G s our best estmate of the average of F(Y(t )) over the trajectory from Y(t =t) to Y(t =t+ t), usng lots of lttle tmesteps t untl we reach t fnal. In the smplest approxmaton called Forward (or Explct) Euler G=F(Y(t)). Ths turns out to be pretty naccurate (just lke the rectangle rule s not a very accurate way to compute numercal ntegrals) and also numercally unstable unless t s very small. The fundamental problem s that we don t know the trajectory Y(t ) all we can do s extrapolate to estmate Y(t ) at future tmes. The extrapolaton at each tmestep ntroduces a lttle bt of error, and then the extrapolaton at the next tmestep s bult on an ncorrect startng pont and an ncorrect estmate of the slope G(Y). If one s not careful, the errors can cumulate n a very unfavorable way, makng the whole procedure numercally unstable (even f the ODEs and the physcal stuaton they descrbe s perfectly stable). There are two general ways of copng wth ths. One s to use relatvely smple estmates of G and just use tny tmesteps t much smaller than any of the physcal tmescales n the system. Ths general approach s called explct, the most famous algorthms that work ths way are called Runge-Kutta methods. The alternatve approach s to use very complcated methods for estmatng G that are guaranteed to be numercally stable, these are called mplct methods because n these methods G s a functon of Y(t+ t), so Eq.(3) becomes an mplct (usually) nonlnear system of equatons, very dffcult to solve. The hope s that by mprovng the numercal stablty, one could use large t s and so reduce the number of tmesteps requred to reach t fnal. The dffcult aspect of chemcal knetcs s that one often has reactve ntermedates wth lfetmes tmes shorter than the relatvely nert startng materals. The short lfetmes of the reactve speces force t n the explct methods to be extremely small, so then a huge number of tmesteps (and a correspondngly huge amount of CPU tme) are requred before the startng materals are consumed sgnfcantly. When the range of tmescales s very large, the ODE system s called stff. Stffness ntroduces many numercal problems, even for mplct ODE solvers; several algorthmc trcks must
5 be used smultaneously n order to solve stff systems wthout ntroducng huge numercal errors. The frst algorthm that could correctly handle stff ODE systems was dscovered by Wllam Gear n the 1970 s. The best programs now avalable for solvng large stff ODE systems are VODE, DASSL/DASPK, and DAEPACK. I beleve the current verson of CHEMKIN calls VODE and DASSL.
Introduction 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 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 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 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 informationNAME and Section No.
Chemstry 391 Fall 2007 Exam I KEY (Monday September 17) 1. (25 Ponts) ***Do 5 out of 6***(If 6 are done only the frst 5 wll be graded)*** a). Defne the terms: open system, closed system and solated system
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 informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationSupplementary 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 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 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 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 informationThermodynamics General
Thermodynamcs General Lecture 1 Lecture 1 s devoted to establshng buldng blocks for dscussng thermodynamcs. In addton, the equaton of state wll be establshed. I. Buldng blocks for thermodynamcs A. Dmensons,
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 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 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 information4.2 Chemical Driving Force
4.2. CHEMICL DRIVING FORCE 103 4.2 Chemcal Drvng Force second effect of a chemcal concentraton gradent on dffuson s to change the nature of the drvng force. Ths s because dffuson changes the bondng n a
More informationGasometric Determination of NaHCO 3 in a Mixture
60 50 40 0 0 5 15 25 35 40 Temperature ( o C) 9/28/16 Gasometrc Determnaton of NaHCO 3 n a Mxture apor Pressure (mm Hg) apor Pressure of Water 1 NaHCO 3 (s) + H + (aq) Na + (aq) + H 2 O (l) + CO 2 (g)
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 informationA quote of the week (or camel of the week): There is no expedience to which a man will not go to avoid the labor of thinking. Thomas A.
A quote of the week (or camel of the week): here s no expedence to whch a man wll not go to avod the labor of thnkng. homas A. Edson Hess law. Algorthm S Select a reacton, possbly contanng specfc compounds
More informationUniversity of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014
Lecture 16 8/4/14 Unversty o Washngton Department o Chemstry Chemstry 452/456 Summer Quarter 214. Real Vapors and Fugacty Henry s Law accounts or the propertes o extremely dlute soluton. s shown n Fgure
More informationConsistency & Convergence
/9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationPh.D. Qualifying Examination in Kinetics and Reactor Design
Knetcs and Reactor Desgn Ph.D.Qualfyng Examnaton January 2006 Instructons Ph.D. Qualfyng Examnaton n Knetcs and Reactor Desgn January 2006 Unversty of Texas at Austn Department of Chemcal Engneerng 1.
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 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 informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
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 informationCHEMISTRY Midterm #2 answer key October 25, 2005
CHEMISTRY 123-01 Mdterm #2 answer key October 25, 2005 Statstcs: Average: 70 pts (70%); Hghest: 97 pts (97%); Lowest: 33 pts (33%) Number of students performng at or above average: 62 (63%) Number of students
More informationLecture 5. Stoichiometry Dimensionless Parameters. Moshe Matalon N X. 00 i. i M i. i=1. i=1
Summer 23 Lecture 5 Stochometry Dmensonless Parameters Stochometry Gven the reacton M dn or n terms of the partal masses dm ( )W N X dn j j M we have seen that there s a relaton between the change n the
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
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 informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
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 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 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 informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationME 501A Seminar in Engineering Analysis Page 1
umercal Solutons of oundary-value Problems n Os ovember 7, 7 umercal Solutons of oundary- Value Problems n Os Larry aretto Mechancal ngneerng 5 Semnar n ngneerng nalyss ovember 7, 7 Outlne Revew stff equaton
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationCurve Fitting with the Least Square Method
WIKI Document Number 5 Interpolaton wth Least Squares Curve Fttng wth the Least Square Method Mattheu Bultelle Department of Bo-Engneerng Imperal College, London Context We wsh to model the postve feedback
More informationOutlet temperature of a WGS reactor (Stage I) for the conversion of CO, applied for the abatement of CO to a fixed value.
Department of Energy Poltecnco d Mlano Va Lambruschn 56 MILAN Exercses of Fundamentals of Chemcal Processes Prof. Ganpero Gropp Exercse utlet temperature of a WGS reactor (Stage I for the converson of
More informationProcess Modeling. Improving or understanding chemical process operation is a major objective for developing a dynamic process model
Process Modelng Improvng or understandng chemcal process operaton s a major objectve for developng a dynamc process model Balance equatons Steady-state balance equatons mass or energy mass or energy enterng
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 informationProblem adapted reduced models based on Reaction-Diffusion Manifolds (REDIMs)
Problem adapted reduced models based on Reacton-Dffuson Manfolds (REDIMs) V Bykov, U Maas Thrty-Second Internatonal Symposum on ombuston, Montreal, anada, 3-8 August, 8 Problem Statement: Smulaton of reactng
More informationbetween standard Gibbs free energies of formation for products and reactants, ΔG! R = ν i ΔG f,i, we
hermodynamcs, Statstcal hermodynamcs, and Knetcs 4 th Edton,. Engel & P. ed Ch. 6 Part Answers to Selected Problems Q6.. Q6.4. If ξ =0. mole at equlbrum, the reacton s not ery far along. hus, there would
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
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 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 informationChemical Equilibrium. Chapter 6 Spontaneity of Reactive Mixtures (gases) Taking into account there are many types of work that a sysem can perform
Ths chapter deals wth chemcal reactons (system) wth lttle or no consderaton on the surroundngs. Chemcal Equlbrum Chapter 6 Spontanety of eactve Mxtures (gases) eactants generatng products would proceed
More informationElectrostatic Potential from Transmembrane Currents
Electrostatc Potental from Transmembrane Currents Let s assume that the current densty j(r, t) s ohmc;.e., lnearly proportonal to the electrc feld E(r, t): j = σ c (r)e (1) wth conductvty σ c = σ c (r).
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationV. Electrostatics. Lecture 25: Diffuse double layer structure
V. Electrostatcs Lecture 5: Dffuse double layer structure MIT Student Last tme we showed that whenever λ D L the electrolyte has a quas-neutral bulk (or outer ) regon at the geometrcal scale L, where there
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 informationV T for n & P = constant
Pchem 365: hermodynamcs -SUMMARY- Uwe Burghaus, Fargo, 5 9 Mnmum requrements for underneath of your pllow. However, wrte your own summary! You need to know the story behnd the equatons : Pressure : olume
More informationCHEMICAL REACTIONS AND DIFFUSION
CHEMICAL REACTIONS AND DIFFUSION A.K.A. NETWORK THERMODYNAMICS BACKGROUND Classcal thermodynamcs descrbes equlbrum states. Non-equlbrum thermodynamcs descrbes steady states. Network thermodynamcs descrbes
More informationIntroduction to Statistical Methods
Introducton to Statstcal Methods Physcs 4362, Lecture #3 hermodynamcs Classcal Statstcal Knetc heory Classcal hermodynamcs Macroscopc approach General propertes of the system Macroscopc varables 1 hermodynamc
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 informationPART I: MULTIPLE CHOICE (32 questions, each multiple choice question has a 2-point value, 64 points total).
CHEMISTRY 123-07 Mdterm #2 answer key November 04, 2010 Statstcs: Average: 68 p (68%); Hghest: 91 p (91%); Lowest: 37 p (37%) Number of students performng at or above average: 58 (53%) Number of students
More informationNice plotting of proteins II
Nce plottng of protens II Fnal remark regardng effcency: It s possble to wrte the Newton representaton n a way that can be computed effcently, usng smlar bracketng that we made for the frst representaton
More informationAnalytical Chemistry Calibration Curve Handout
I. Quck-and Drty Excel Tutoral Analytcal Chemstry Calbraton Curve Handout For those of you wth lttle experence wth Excel, I ve provded some key technques that should help you use the program both for problem
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationTHE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens
THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of
More informationG4023 Mid-Term Exam #1 Solutions
Exam1Solutons.nb 1 G03 Md-Term Exam #1 Solutons 1-Oct-0, 1:10 p.m to :5 p.m n 1 Pupn Ths exam s open-book, open-notes. You may also use prnt-outs of the homework solutons and a calculator. 1 (30 ponts,
More informationFor now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.
Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson
More informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
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 informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
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 informationName ID # For relatively dilute aqueous solutions the molality and molarity are approximately equal.
Name ID # 1 CHEMISTRY 212, Lect. Sect. 002 Dr. G. L. Roberts Exam #1/Sprng 2000 Thursday, February 24, 2000 CLOSED BOOK EXM No notes or books allowed. Calculators may be used. tomc masses of nterest are
More informationParametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010
Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton
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 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 informationAppendix II Summary of Important Equations
W. M. Whte Geochemstry Equatons of State: Ideal GasLaw: Coeffcent of Thermal Expanson: Compressblty: Van der Waals Equaton: The Laws of Thermdynamcs: Frst Law: Appendx II Summary of Important Equatons
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 informationReview of Classical Thermodynamics
Revew of Classcal hermodynamcs Physcs 4362, Lecture #1, 2 Syllabus What s hermodynamcs? 1 [A law] s more mpressve the greater the smplcty of ts premses, the more dfferent are the knds of thngs t relates,
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 informationTurbulent Nonpremixed Flames
School of Aerospace Engneerng Turbulent Nonpremxed Flames Jerry Setzman. 5 Mole Fracton.15.1.5 CH4 HO HCO x 1 Temperature Methane Flame.1..3 Dstance (cm) 15 1 5 Temperature (K) TurbulentNonpremxed -1 School
More informationProblem Set 9 Solutions
Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem
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 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 informationLab session: numerical simulations of sponateous polarization
Lab sesson: numercal smulatons of sponateous polarzaton Emerc Boun & Vncent Calvez CNRS, ENS Lyon, France CIMPA, Hammamet, March 2012 Spontaneous cell polarzaton: the 1D case The Hawkns-Voturez model for
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
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 informationCase A. P k = Ni ( 2L i k 1 ) + (# big cells) 10d 2 P k.
THE CELLULAR METHOD In ths lecture, we ntroduce the cellular method as an approach to ncdence geometry theorems lke the Szemeréd-Trotter theorem. The method was ntroduced n the paper Combnatoral complexty
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationLimited Dependent Variables
Lmted Dependent Varables. What f the left-hand sde varable s not a contnuous thng spread from mnus nfnty to plus nfnty? That s, gven a model = f (, β, ε, where a. s bounded below at zero, such as wages
More informationElectrochemistry Thermodynamics
CHEM 51 Analytcal Electrochemstry Chapter Oct 5, 016 Electrochemstry Thermodynamcs Bo Zhang Department of Chemstry Unversty of Washngton Seattle, WA 98195 Former SEAC presdent Andy Ewng sellng T-shrts
More informationEntropy generation in a chemical reaction
Entropy generaton n a chemcal reacton E Mranda Área de Cencas Exactas COICET CCT Mendoza 5500 Mendoza, rgentna and Departamento de Físca Unversdad aconal de San Lus 5700 San Lus, rgentna bstract: Entropy
More informationE91: Dynamics. E91: Dynamics. Numerical Integration & State Space Representation
E91: Dnamcs Numercal Integraton & State Space Representaton The Algorthm In steps of Δ f ( new ) f ( old ) df ( d old ) Δ Numercal Integraton of ODEs d d f() h Δ Intal value problem: Gven the ntal state
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 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 informationChapter - 2. Distribution System Power Flow Analysis
Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load
More informationPolynomial Regression Models
LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance
More informationCalifornia State Science Fair
Calforna State Scence Far Mathematcal Modelng of Real World Systems Part 1 -- Explorng on Chaos on Your Computer Edward Ruth drruth@x.netcom.com 1) Introducton Mathematcs s the language that scentsts use
More informationUniversity of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015
Lecture 2. 1/07/15-1/09/15 Unversty of Washngton Department of Chemstry Chemstry 453 Wnter Quarter 2015 We are not talkng about truth. We are talkng about somethng that seems lke truth. The truth we want
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 informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More information