ME425/525: Advanced Topics in Building Science

Size: px
Start display at page:

Download "ME425/525: Advanced Topics in Building Science"

Transcription

1 ME425/525: Advanced Topics in Building Science Indoor environmenal qualiy for susainable buildings: Lecure 6 Dr. Ellio T. Gall, Ph.D.

2 Lecure 6 Today s objecives o Error propagaion Apply o SS soluion (venilaion, emission, air-cleaning) o Chemical reacions Heerogeneous Homogeneous o Compeiion beween sinks Applicaion o indoor spaces, coupled mass balances o Adsorpion/desorpion processes o Quesions?

3 Compeiive phenomena Compeiion amongs sinks We solved for C() under scenario of o Infilraion (Q C o ) o Exfilraion (Q C) o Emission (E) o Air cleaning o Heerogeneous chemisry o Homogeneous chemisry C = C =0 e α + λ α C o + E αv 1 e α Where α = λ + nq c V + σ j k c C j + σ i v d,i A i V C o Q C hom. chem Q c Qc he. chem E Q

4 A realisic example We wan o rack indoor ozone dynamics o Need ime-varian mass balance Le s assume he following scenario: o A 100 m 3, well-mixed indoor space o A ceiling heigh of 2.5 m o An air exchange rae of 0.5 h -1 o No indoor ozone emissions (E = 0)

5 Air cleaning for ozone Le s assume we operae on low speed and ha ozone is removed wih an efficiency of 50% (ypical for a mid-life acivaed carbon filer). Wha is he CADR?

6 Homogeneous chemisry Wha else do we need o model our sysem? Homogeneous chemisry o Rae consan o Concenraion of second gas-phase reacan Le s assume a simple sysem, where ozone is reacing wih one oher compound, so + R hom = k c C j Wha does his C represen he concenraion of?

7 Ozone and limonene Ozone reacs rapidly wih he erpene limonene o Ubiquious indoors o Obained commercially from cirus o Presen in cleaning producs Scening agen (lemony fresh!) i. Consumer/self-care producs ii. Air fresheners Solven properies i. Cleaning producs o 45,000 ons/year produced o Iself may be an irrian o Bu, ozone + limonene yields harmful producs For now, we ll focus on jus he reacans

8 Indoor VOCs Wha are ypical indoor limonene levels? o Sudy of personal 48-h exposures in 201 individuals Mean personal limonene concenraion was found o be 32 μg/m 3 R hom = k c C j Do you hink his indoor limonene level is higher or lower han oudoors?

9 Reacion rae consan R hom = k c C j Noe unis! Rae consan is given in 1/(ppbs), while we ve have ozone, limonene levels in µg/m 3 (in his example) Many resources for hese: NIST Chemisry Webook Je Propulsion Laboraory Review aricles in he lieraure

10 Wha are ypical oudoor ozone levels? Oudoor ozone C = C =0 e λ + λ α C o + E αv 1 e α Box plos of daily 8-h maximum average ozone concenraions in six locaions across Singapore. Le s assume an average of 35 µg/m 3 of ozone. The U.S. EPA NAAQS is 70 ppb, how does our assumpion compare?

11 Heerogeneous chemisry Ozone reacs wih surfaces othree prominen surfaces Drywall (pained), Carpe, Ceiling ile i v d,i A i V

12 Heerogeneous chemisry Maerial surface imporan; T, RH less so Increase wih increase in mixing condiion, why? i v d,i A i V For 3 maerials: 1) carpe, 2) ceiling ile, 3) pained drywall

13 Heerogeneous chemisry Each v d is associaed wih an area o Assume floor is carpe Area = 5 8 = 40 m 2 o Assume walls are pained drywall Area = = 40 m 2 o Assume ceiling is ceiling ile Area = 5 8 = 40 m 2 i v d,i A i V

14 Organize he parameers lambda 0.5 h-1 air exchange rae V 100 m 3 Volume of indoor space E-ozone 0 ppb/h*m^3 ozone emission rae CADR m^3/h air cleaner clean air delivery rae kc-j 5.10E-06 1/(ppb-s) gas-phase ozone-limonene rae consan kc-j 1.84E-02 1/(ppb-h) vd,c 5.5 m/h ozone-carpe deposiion velociy vd,c 2.5 m/h ozone-ceiling ile deposiion velociy vd, pdw 0.5 m/h ozone - pained drywall deposiion velociy A, c 40 m 2 carpe area A, c 40 m 2 ceiling ile area A, pdw 65 m 2 pained drywall ara Co 35 μg/m 3 Co 17.5 ppb oudoor ozone level Clim 32 μg/m 3 Clim 5.6 ppb indoor limonene concenraion C= ppb iniial indoor ozone level C o Q C hom. chem Q c = convered parameer Qc he. chem E Q C = C =0 e λ + λ α C o + E αv 1 e α Where α = λ + nq c V + σ j k c C j + σ i v d,i A i V

15 Solving he inegraed M.B. C = C =0 e α + λ α C o + E αv 1 e α Where α = λ + nq c V + σ j k c C j + σ i v d,i A i V C o Q C hom. chem Qc he. chem E Q Q c Le s firs evaluae α wha should he unis of alpha be? α = 0.5 h m3 /h 100m 3 m ppb ppb + h 40m m h 40m m h 40m2 h 100m 3 α = 5.4 h -1

16 Indoor O3 concenraion (ppb) Simulaion resuls C-O3: AC, Hom, He C-O3: Hom, He C-O3: He Time (h) α = 4.0 h -1 α = 4.1 h -1 α = 5.4 h -1 Lower alpha resuls in higher concenraion (hink of alpha as a firs-order loss erm) Wha is he dominan removal mechanism?

17 Inerpreing resuls Compue he conribuion o he overall sink from each erm in he equaion for alpha: Conribuion o indoor sinks 3.53 h h h h -1 Air exchange Air cleaning Heerogegneous chemisry Homogeneous chemisry Surface reacions conribue mos subsanially o he removal of indoor ozone (in his example, bu also in general).

18 Indoor O 3 concenraion (ppb) Compare o α=λ C-O3: AC, Hom, He C-O3: Hom, He C-O3: He C-O3: No sinks Time (h) Wha is he implicaion of α=λ? Wha is he seady-sae concenraion in such a scenario, assuming E = 0?

19 Simplifying assumpions Wha are wo unrealisic assumpions we have made in performing his analysis? Indoor limonene is consan R hom = k c C limonene C ozone Primary ozonide Chemisry is complex, bu for every 1 mole of ozone removed, 1 mole of limonene is removed -> limonene levels will vary as well!

20 Simplifying assumpions Wha are wo unrealisic assumpions we have made in performing his analysis? Oudoor ozone is consan Island average O3 ( g/m 3 ) Ozone Temperaure Hour of day (average from 1/1/ /12/2007) Figure 1. Average O 3 and emp. in Singapore Temperaure (deg C) Oudoor ozone highly variable, diurnal rend ha is relian on solar radiaion o iniiae he chemisry (phoochemisry)

21 Implicaions of our assumpions Are we over or under-predicing loss of ozone due o limonene reacion? o Bu here are many oher compounds ozone could reac wih Assuming ozone is consan? o Could be over-predicing or under-predicing o Time-averaged soluions are common Appropriae if you ime-averaging is longer han your sysem s dynamics o approaching seadysae

22 How o resolve? We are essenially assuming limonene is a seady-sae dc lim d 0 Wha should we do? o Mass balance on limonene o Assume we have Air exchange An emission source Q Homogeneous chemisry E Q c C o C hom. chem

23 A dynamic mass balance for limonene Coupled mass balances dc lim d = λc o,lim λc lim () + E O3 V k lim O3C lim ()C O3 () A dynamic mass balance for O3 dc O3 d = λc o,o3 λc O3 () + E O3 V k lim O3C lim ()C O3 () i v d,i A i V These equaions are coupled Wih respec o limonene concenraion And ozone concenraion

24 Oudoor levels are also ime-varian A dynamic mass balance for limonene dc lim d = λc o,lim () λc lim () + E O3 V k lim O3C lim ()C O3 () A dynamic mass balance for O3 (noe no air-cleaner) dc O3 d = λc o,o3 () λc O3 () + E O3 V k lim O3C lim ()C O3 () i v d,i A i V Oudoor levels will also vary wih respec o ime

25 Numerical soluions A dynamic mass balance for limonene dc lim d = λc o,lim () λc lim () + E O3 V k lim O3C lim ()C O3 () A dynamic mass balance for O3 dc O3 d = λc o,o3 () λc O3 () + E O3 V k lim O3C lim ()C O3 () i v d,i A i V C O3 Analyical soluions for some simple scenarios may exis Numerical soluions are ofen pursued for models

26 Forward euler mehod Forward euler mehod Break derivaive ino discree erms For limonene: Where h is he ime-sep C + lim C lim = λc o,lim λc lim + E O3 V k lim O3C lim C O3 A dynamic mass balance for O3 C + O3 C O3 = λc o,o3 λc O3 + E O3 V k lim O3C lim C O3 i v d,i A i V C O3

27 Forward Euler Mehod Solve algebraically for C(+Δ) C + lim = λc o,lim λc lim + E O3 V k lim O3C lim C O3 + C lim + C O3 = λc o,o3 λc O3 + E O3 V k lim O3C lim C O3 i v d,i A i V C O3 + C O3 Upsides o Forward Euler: Simple Downsides: Unsable, requires very small ime-sep (herefore compuaion ime) Oher approaches: Modified Euler, Crank-Nicholson, Runge-Kua, ec.

28 Indoor ozone concenraion (ppb) Discreizing soluions and error Discree, Δ = 0.05 h Discree, Δ = 0.1 h Discree, Δ = h Analyical Time (h)

Ordinary Differential Equations

Ordinary Differential Equations Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

CHEMICAL KINETICS: 1. Rate Order Rate law Rate constant Half-life Temperature Dependence

CHEMICAL KINETICS: 1. Rate Order Rate law Rate constant Half-life Temperature Dependence CHEMICL KINETICS: Rae Order Rae law Rae consan Half-life Temperaure Dependence Chemical Reacions Kineics Chemical ineics is he sudy of ime dependence of he change in he concenraion of reacans and producs.

More information

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species

More information

Advanced Organic Chemistry

Advanced Organic Chemistry Lalic, G. Chem 53A Chemisry 53A Advanced Organic Chemisry Lecure noes 1 Kineics: A racical Approach Simple Kineics Scenarios Fiing Experimenal Daa Using Kineics o Deermine he Mechanism Doughery, D. A.,

More information

Subway stations energy and air quality management

Subway stations energy and air quality management Subway saions energy and air qualiy managemen wih sochasic opimizaion Trisan Rigau 1,2,4, Advisors: P. Carpenier 3, J.-Ph. Chancelier 2, M. De Lara 2 EFFICACITY 1 CERMICS, ENPC 2 UMA, ENSTA 3 LISIS, IFSTTAR

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

Announcements: Warm-up Exercise:

Announcements: Warm-up Exercise: Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple

More information

AP Chemistry--Chapter 12: Chemical Kinetics

AP Chemistry--Chapter 12: Chemical Kinetics AP Chemisry--Chaper 12: Chemical Kineics I. Reacion Raes A. The area of chemisry ha deals wih reacion raes, or how fas a reacion occurs, is called chemical kineics. B. The rae of reacion depends on he

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18 A Firs ourse on Kineics and Reacion Engineering lass 19 on Uni 18 Par I - hemical Reacions Par II - hemical Reacion Kineics Where We re Going Par III - hemical Reacion Engineering A. Ideal Reacors B. Perfecly

More information

Economics 8105 Macroeconomic Theory Recitation 6

Economics 8105 Macroeconomic Theory Recitation 6 Economics 8105 Macroeconomic Theory Reciaion 6 Conor Ryan Ocober 11h, 2016 Ouline: Opimal Taxaion wih Governmen Invesmen 1 Governmen Expendiure in Producion In hese noes we will examine a model in which

More information

CHEMICAL KINETICS Rate Order Rate law Rate constant Half-life Molecularity Elementary. Complex Temperature dependence, Steady-state Approximation

CHEMICAL KINETICS Rate Order Rate law Rate constant Half-life Molecularity Elementary. Complex Temperature dependence, Steady-state Approximation CHEMICL KINETICS Rae Order Rae law Rae consan Half-life Moleculariy Elemenary Complex Temperaure dependence, Seady-sae pproximaion Chemical Reacions Kineics Chemical ineics is he sudy of ime dependence

More information

1. Consider a pure-exchange economy with stochastic endowments. The state of the economy

1. Consider a pure-exchange economy with stochastic endowments. The state of the economy Answer 4 of he following 5 quesions. 1. Consider a pure-exchange economy wih sochasic endowmens. The sae of he economy in period, 0,1,..., is he hisory of evens s ( s0, s1,..., s ). The iniial sae is given.

More information

Chapter 13 Homework Answers

Chapter 13 Homework Answers Chaper 3 Homework Answers 3.. The answer is c, doubling he [C] o while keeping he [A] o and [B] o consan. 3.2. a. Since he graph is no linear, here is no way o deermine he reacion order by inspecion. A

More information

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis Chaper EEE83 EEE3 Chaper # EEE83 EEE3 Linear Conroller Design and Sae Space Analysis Ordinary Differenial Equaions.... Inroducion.... Firs Order ODEs... 3. Second Order ODEs... 7 3. General Maerial...

More information

Kinematics and kinematic functions

Kinematics and kinematic functions Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables and vice versa Direc Posiion

More information

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal? EE 35 Noes Gürdal Arslan CLASS (Secions.-.2) Wha is a signal? In his class, a signal is some funcion of ime and i represens how some physical quaniy changes over some window of ime. Examples: velociy of

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

More information

2) Of the following questions, which ones are thermodynamic, rather than kinetic concepts?

2) Of the following questions, which ones are thermodynamic, rather than kinetic concepts? AP Chemisry Tes (Chaper 12) Muliple Choice (40%) 1) Which of he following is a kineic quaniy? A) Enhalpy B) Inernal Energy C) Gibb s free energy D) Enropy E) Rae of reacion 2) Of he following quesions,

More information

EE363 homework 1 solutions

EE363 homework 1 solutions EE363 Prof. S. Boyd EE363 homework 1 soluions 1. LQR for a riple accumulaor. We consider he sysem x +1 = Ax + Bu, y = Cx, wih 1 1 A = 1 1, B =, C = [ 1 ]. 1 1 This sysem has ransfer funcion H(z) = (z 1)

More information

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based

More information

Open loop vs Closed Loop. Example: Open Loop. Example: Feedforward Control. Advanced Control I

Open loop vs Closed Loop. Example: Open Loop. Example: Feedforward Control. Advanced Control I Open loop vs Closed Loop Advanced I Moor Command Movemen Overview Open Loop vs Closed Loop Some examples Useful Open Loop lers Dynamical sysems CPG (biologically inspired ), Force Fields Feedback conrol

More information

2.4 Cuk converter example

2.4 Cuk converter example 2.4 Cuk converer example C 1 Cuk converer, wih ideal swich i 1 i v 1 2 1 2 C 2 v 2 Cuk converer: pracical realizaion using MOSFET and diode C 1 i 1 i v 1 2 Q 1 D 1 C 2 v 2 28 Analysis sraegy This converer

More information

Lecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples

Lecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples EE263 Auumn 27-8 Sephen Boyd Lecure 1 Overview course mechanics ouline & opics wha is a linear dynamical sysem? why sudy linear sysems? some examples 1 1 Course mechanics all class info, lecures, homeworks,

More information

Selective peracetic acid determination in the presence of hydrogen peroxide using the molecular absorption properties of catalase

Selective peracetic acid determination in the presence of hydrogen peroxide using the molecular absorption properties of catalase Analyical and Bioanalyical Chemisry Elecronic Supplemenary Maerial Selecive peraceic acid deerminaion in he presence of hydrogen peroxide using he molecular absorpion properies of caalase Javier Galbán,

More information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3. Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

Solutions Problem Set 3 Macro II (14.452)

Solutions Problem Set 3 Macro II (14.452) Soluions Problem Se 3 Macro II (14.452) Francisco A. Gallego 04/27/2005 1 Q heory of invesmen in coninuous ime and no uncerainy Consider he in nie horizon model of a rm facing adjusmen coss o invesmen.

More information

ECE 2100 Circuit Analysis

ECE 2100 Circuit Analysis ECE 1 Circui Analysis Lesson 35 Chaper 8: Second Order Circuis Daniel M. Liynski, Ph.D. ECE 1 Circui Analysis Lesson 3-34 Chaper 7: Firs Order Circuis (Naural response RC & RL circuis, Singulariy funcions,

More information

Sliding Mode Controller for Unstable Systems

Sliding Mode Controller for Unstable Systems S. SIVARAMAKRISHNAN e al., Sliding Mode Conroller for Unsable Sysems, Chem. Biochem. Eng. Q. 22 (1) 41 47 (28) 41 Sliding Mode Conroller for Unsable Sysems S. Sivaramakrishnan, A. K. Tangirala, and M.

More information

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 3 Signals & Sysems Prof. Mark Fowler Noe Se #2 Wha are Coninuous-Time Signals??? Reading Assignmen: Secion. of Kamen and Heck /22 Course Flow Diagram The arrows here show concepual flow beween ideas.

More information

Problem Set on Differential Equations

Problem Set on Differential Equations Problem Se on Differenial Equaions 1. Solve he following differenial equaions (a) x () = e x (), x () = 3/ 4. (b) x () = e x (), x (1) =. (c) xe () = + (1 x ()) e, x () =.. (An asse marke model). Le p()

More information

Embedded Systems and Software. A Simple Introduction to Embedded Control Systems (PID Control)

Embedded Systems and Software. A Simple Introduction to Embedded Control Systems (PID Control) Embedded Sysems and Sofware A Simple Inroducion o Embedded Conrol Sysems (PID Conrol) Embedded Sysems and Sofware, ECE:3360. The Universiy of Iowa, 2016 Slide 1 Acknowledgemens The maerial in his lecure

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

Ordinary dierential equations

Ordinary dierential equations Chaper 5 Ordinary dierenial equaions Conens 5.1 Iniial value problem........................... 31 5. Forward Euler's mehod......................... 3 5.3 Runge-Kua mehods.......................... 36

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION DOI: 0.038/NCLIMATE893 Temporal resoluion and DICE * Supplemenal Informaion Alex L. Maren and Sephen C. Newbold Naional Cener for Environmenal Economics, US Environmenal Proecion

More information

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),

More information

Today in Physics 218: radiation reaction

Today in Physics 218: radiation reaction Today in Physics 18: radiaion reacion Radiaion reacion The Abraham-Lorenz formula; radiaion reacion force The pah of he elecron in oday s firs example (radial decay grealy exaggeraed) 6 March 004 Physics

More information

INDEX. Transient analysis 1 Initial Conditions 1

INDEX. Transient analysis 1 Initial Conditions 1 INDEX Secion Page Transien analysis 1 Iniial Condiions 1 Please inform me of your opinion of he relaive emphasis of he review maerial by simply making commens on his page and sending i o me a: Frank Mera

More information

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)

More information

Voltage/current relationship Stored Energy. RL / RC circuits Steady State / Transient response Natural / Step response

Voltage/current relationship Stored Energy. RL / RC circuits Steady State / Transient response Natural / Step response Review Capaciors/Inducors Volage/curren relaionship Sored Energy s Order Circuis RL / RC circuis Seady Sae / Transien response Naural / Sep response EE4 Summer 5: Lecure 5 Insrucor: Ocavian Florescu Lecure

More information

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

More information

5.1 - Logarithms and Their Properties

5.1 - Logarithms and Their Properties Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We

More information

3.6 Derivatives as Rates of Change

3.6 Derivatives as Rates of Change 3.6 Derivaives as Raes of Change Problem 1 John is walking along a sraigh pah. His posiion a he ime >0 is given by s = f(). He sars a =0from his house (f(0) = 0) and he graph of f is given below. (a) Describe

More information

Chapter 14 Homework Answers

Chapter 14 Homework Answers 4. Suden responses will vary. (a) combusion of gasoline (b) cooking an egg in boiling waer (c) curing of cemen Chaper 4 Homework Answers 4. A collision beween only wo molecules is much more probable han

More information

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1 SZG Macro 2011 Lecure 3: Dynamic Programming SZG macro 2011 lecure 3 1 Background Our previous discussion of opimal consumpion over ime and of opimal capial accumulaion sugges sudying he general decision

More information

CHEAPEST PMT ONLINE TEST SERIES AIIMS/NEET TOPPER PREPARE QUESTIONS

CHEAPEST PMT ONLINE TEST SERIES AIIMS/NEET TOPPER PREPARE QUESTIONS CHEAPEST PMT ONLINE TEST SERIES AIIMS/NEET TOPPER PREPARE QUESTIONS For more deails see las page or conac @aimaiims.in Physics Mock Tes Paper AIIMS/NEET 07 Physics 06 Saurday Augus 0 Uni es : Moion in

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

CSE 3802 / ECE Numerical Methods in Scientific Computation. Jinbo Bi. Department of Computer Science & Engineering

CSE 3802 / ECE Numerical Methods in Scientific Computation. Jinbo Bi. Department of Computer Science & Engineering CSE 3802 / ECE 3431 Numerical Mehods in Scienific Compuaion Jinbo Bi Deparmen of Compuer Science & Engineering hp://www.engr.uconn.edu/~jinbo 1 Ph.D in Mahemaics The Insrucor Previous professional experience:

More information

2. For a one-point fixed time method, a pseudo-first order reaction obeys the equation 0.309

2. For a one-point fixed time method, a pseudo-first order reaction obeys the equation 0.309 Chaper 3. To derive an appropriae equaion we firs noe he following general relaionship beween he concenraion of A a ime, [A], he iniial concenraion of A, [A], and he concenraion of P a ime, [P] [ P] Subsiuing

More information

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder# .#W.#Erickson# Deparmen#of#Elecrical,#Compuer,#and#Energy#Engineering# Universiy#of#Colorado,#Boulder# Chaper 2 Principles of Seady-Sae Converer Analysis 2.1. Inroducion 2.2. Inducor vol-second balance,

More information

Unsteady Flow Problems

Unsteady Flow Problems School of Mechanical Aerospace and Civil Engineering Unseady Flow Problems T. J. Craf George Begg Building, C41 TPFE MSc CFD-1 Reading: J. Ferziger, M. Peric, Compuaional Mehods for Fluid Dynamics H.K.

More information

Math 4600: Homework 11 Solutions

Math 4600: Homework 11 Solutions Mah 46: Homework Soluions Gregory Handy [.] One of he well-known phenomenological (capuring he phenomena, bu no necessarily he mechanisms) models of cancer is represened by Gomperz equaion dn d = bn ln(n/k)

More information

6.003 Homework 1. Problems. Due at the beginning of recitation on Wednesday, February 10, 2010.

6.003 Homework 1. Problems. Due at the beginning of recitation on Wednesday, February 10, 2010. 6.003 Homework Due a he beginning of reciaion on Wednesday, February 0, 200. Problems. Independen and Dependen Variables Assume ha he heigh of a waer wave is given by g(x v) where x is disance, v is velociy,

More information

BU Macro BU Macro Fall 2008, Lecture 4

BU Macro BU Macro Fall 2008, Lecture 4 Dynamic Programming BU Macro 2008 Lecure 4 1 Ouline 1. Cerainy opimizaion problem used o illusrae: a. Resricions on exogenous variables b. Value funcion c. Policy funcion d. The Bellman equaion and an

More information

Initial Value Problems

Initial Value Problems Iniial Value Problems ChEn 2450 d d f(, ) (0) 0 6 ODE.key - November 26, 2014 Example - Cooking a Lobser Assumpions: The lobser remains a a uniform emperaure. This implies ha he hermal conduciviy of he

More information

ψ ( t) = c n ( t) t n ( )ψ( ) t ku t,t 0 ψ I V kn

ψ ( t) = c n ( t) t n ( )ψ( ) t ku t,t 0 ψ I V kn MIT Deparmen of Chemisry 5.74, Spring 4: Inroducory Quanum Mechanics II p. 33 Insrucor: Prof. Andrei Tokmakoff PERTURBATION THEORY Given a Hamilonian H ( ) = H + V ( ) where we know he eigenkes for H H

More information

המחלקה : ביולוגיה מולקולרית הפקולטה למדעי הטבע

המחלקה : ביולוגיה מולקולרית הפקולטה למדעי הטבע טל : 03-9066345 פקס : 03-9374 ראש המחלקה: ד"ר אלברט פנחסוב המחלקה : ביולוגיה מולקולרית הפקולטה למדעי הטבע Course Name: Physical Chemisry - (for Molecular Biology sudens) כימיה פיזיקאלית )לסטודנטים לביולוגיה

More information

Homework 2 Solutions

Homework 2 Solutions Mah 308 Differenial Equaions Fall 2002 & 2. See he las page. Hoework 2 Soluions 3a). Newon s secon law of oion says ha a = F, an we know a =, so we have = F. One par of he force is graviy, g. However,

More information

ENV 6015 Solution to Mixing Problem Set

ENV 6015 Solution to Mixing Problem Set EN 65 Soluion o ixing Problem Se. A slug of dye ( ) is injeced ino a single ank wih coninuous mixing. The flow in and ou of he ank is.5 gpm. The ank volume is 5 gallons. When will he dye concenraion equal

More information

Chapter 2: Principles of steady-state converter analysis

Chapter 2: Principles of steady-state converter analysis Chaper 2 Principles of Seady-Sae Converer Analysis 2.1. Inroducion 2.2. Inducor vol-second balance, capacior charge balance, and he small ripple approximaion 2.3. Boos converer example 2.4. Cuk converer

More information

PROBLEMS FOR MATH 162 If a problem is starred, all subproblems are due. If only subproblems are starred, only those are due. SLOPES OF TANGENT LINES

PROBLEMS FOR MATH 162 If a problem is starred, all subproblems are due. If only subproblems are starred, only those are due. SLOPES OF TANGENT LINES PROBLEMS FOR MATH 6 If a problem is sarred, all subproblems are due. If onl subproblems are sarred, onl hose are due. 00. Shor answer quesions. SLOPES OF TANGENT LINES (a) A ball is hrown ino he air. Is

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Ordinary differential equations. Phys 750 Lecture 7

Ordinary differential equations. Phys 750 Lecture 7 Ordinary differenial equaions Phys 750 Lecure 7 Ordinary Differenial Equaions Mos physical laws are expressed as differenial equaions These come in hree flavours: iniial-value problems boundary-value problems

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 3 Signals & Sysems Prof. Mark Fowler Noe Se # Wha are Coninuous-Time Signals??? /6 Coninuous-Time Signal Coninuous Time (C-T) Signal: A C-T signal is defined on he coninuum of ime values. Tha is:

More information

Sub Module 2.6. Measurement of transient temperature

Sub Module 2.6. Measurement of transient temperature Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,

More information

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation Hea (iffusion) Equaion erivaion of iffusion Equaion The fundamenal mass balance equaion is I P O L A ( 1 ) where: I inpus P producion O oupus L losses A accumulaion Assume ha no chemical is produced or

More information

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem) Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =

More information

ECE 2100 Circuit Analysis

ECE 2100 Circuit Analysis ECE 1 Circui Analysis Lesson 37 Chaper 8: Second Order Circuis Discuss Exam Daniel M. Liynski, Ph.D. Exam CH 1-4: On Exam 1; Basis for work CH 5: Operaional Amplifiers CH 6: Capaciors and Inducor CH 7-8:

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

EE 435. Lecture 31. Absolute and Relative Accuracy DAC Design. The String DAC

EE 435. Lecture 31. Absolute and Relative Accuracy DAC Design. The String DAC EE 435 Lecure 3 Absolue and Relaive Accuracy DAC Design The Sring DAC . Review from las lecure. DFT Simulaion from Malab Quanizaion Noise DACs and ADCs generally quanize boh ampliude and ime If convering

More information

Section 4.4 Logarithmic Properties

Section 4.4 Logarithmic Properties Secion. Logarihmic Properies 59 Secion. Logarihmic Properies In he previous secion, we derived wo imporan properies of arihms, which allowed us o solve some asic eponenial and arihmic equaions. Properies

More information

Module 4: Time Response of discrete time systems Lecture Note 2

Module 4: Time Response of discrete time systems Lecture Note 2 Module 4: Time Response of discree ime sysems Lecure Noe 2 1 Prooype second order sysem The sudy of a second order sysem is imporan because many higher order sysem can be approimaed by a second order model

More information

1 Answers to Final Exam, ECN 200E, Spring

1 Answers to Final Exam, ECN 200E, Spring 1 Answers o Final Exam, ECN 200E, Spring 2004 1. A good answer would include he following elemens: The equiy premium puzzle demonsraed ha wih sandard (i.e ime separable and consan relaive risk aversion)

More information

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters?

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters? of Topics wih MATLAB Shan He School for Compuaional Science Universi of Birmingham Module 6-3836: Compuaional Modelling wih MATLAB Wha will we learn from his lecure Aim of analsis: Aim of analsis. Some

More information

ADVANCED MATHEMATICS FOR ECONOMICS /2013 Sheet 3: Di erential equations

ADVANCED MATHEMATICS FOR ECONOMICS /2013 Sheet 3: Di erential equations ADVANCED MATHEMATICS FOR ECONOMICS - /3 Shee 3: Di erenial equaions Check ha x() =± p ln(c( + )), where C is a posiive consan, is soluion of he ODE x () = Solve he following di erenial equaions: (a) x

More information

PHYSICS 220 Lecture 02 Motion, Forces, and Newton s Laws Textbook Sections

PHYSICS 220 Lecture 02 Motion, Forces, and Newton s Laws Textbook Sections PHYSICS 220 Lecure 02 Moion, Forces, and Newon s Laws Texbook Secions 2.2-2.4 Lecure 2 Purdue Universiy, Physics 220 1 Overview Las Lecure Unis Scienific Noaion Significan Figures Moion Displacemen: Δx

More information

Basilio Bona ROBOTICA 03CFIOR 1

Basilio Bona ROBOTICA 03CFIOR 1 Indusrial Robos Kinemaics 1 Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables

More information

6.003 Homework #8 Solutions

6.003 Homework #8 Solutions 6.003 Homework #8 Soluions Problems. Fourier Series Deermine he Fourier series coefficiens a k for x () shown below. x ()= x ( + 0) 0 a 0 = 0 a k = e /0 sin(/0) for k 0 a k = π x()e k d = 0 0 π e 0 k d

More information

Probabilistic Robotics

Probabilistic Robotics Probabilisic Roboics Bayes Filer Implemenaions Gaussian filers Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel Gaussians : ~ π e p N p - Univariae / / : ~ μ μ μ e p Ν p d π Mulivariae

More information

04. Kinetics of a second order reaction

04. Kinetics of a second order reaction 4. Kineics of a second order reacion Imporan conceps Reacion rae, reacion exen, reacion rae equaion, order of a reacion, firs-order reacions, second-order reacions, differenial and inegraed rae laws, Arrhenius

More information

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore Soluions of Sample Problems for Third In-Class Exam Mah 6, Spring, Professor David Levermore Compue he Laplace ransform of f e from is definiion Soluion The definiion of he Laplace ransform gives L[f]s

More information