ME425/525: Advanced Topics in Building Science

Similar documents
Ordinary Differential Equations

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

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

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

Advanced Organic Chemistry

Subway stations energy and air quality management

Numerical Dispersion

Announcements: Warm-up Exercise:

AP Chemistry--Chapter 12: Chemical Kinetics

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...

Vehicle Arrival Models : Headway

System of Linear Differential Equations

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

LAPLACE TRANSFORM AND TRANSFER FUNCTION

Chapter 2. First Order Scalar Equations

Some Basic Information about M-S-D Systems

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

EXERCISES FOR SECTION 1.5

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

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

Economics 8105 Macroeconomic Theory Recitation 6

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

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

Chapter 13 Homework Answers

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

Kinematics and kinematic functions

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

Solutions to Assignment 1

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

EE363 homework 1 solutions

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

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

2.4 Cuk converter example

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

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

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

Notes on Kalman Filtering

Solutions Problem Set 3 Macro II (14.452)

ECE 2100 Circuit Analysis

Sliding Mode Controller for Unstable Systems

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

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

EECE 301 Signals & Systems Prof. Mark Fowler

Problem Set on Differential Equations

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

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

Ordinary dierential equations

KINEMATICS IN ONE DIMENSION

SUPPLEMENTARY INFORMATION

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

Today in Physics 218: radiation reaction

INDEX. Transient analysis 1 Initial Conditions 1

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

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

Traveling Waves. Chapter Introduction

5.1 - Logarithms and Their Properties

3.6 Derivatives as Rates of Change

Chapter 14 Homework Answers

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

CHEAPEST PMT ONLINE TEST SERIES AIIMS/NEET TOPPER PREPARE QUESTIONS

The Arcsine Distribution

Class Meeting # 10: Introduction to the Wave Equation

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

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

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

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

Unsteady Flow Problems

Math 4600: Homework 11 Solutions

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

BU Macro BU Macro Fall 2008, Lecture 4

Initial Value Problems

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

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

Homework 2 Solutions

ENV 6015 Solution to Mixing Problem Set

Chapter 2: Principles of steady-state converter analysis

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

Two Coupled Oscillators / Normal Modes

Ordinary differential equations. Phys 750 Lecture 7

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

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

EECE 301 Signals & Systems Prof. Mark Fowler

Sub Module 2.6. Measurement of transient temperature

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

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

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

ECE 2100 Circuit Analysis

15. Vector Valued Functions

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

Section 4.4 Logarithmic Properties

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

1 Answers to Final Exam, ECN 200E, Spring

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

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

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

Basilio Bona ROBOTICA 03CFIOR 1

6.003 Homework #8 Solutions

Probabilistic Robotics

04. Kinetics of a second order reaction

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

Transcription:

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

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?

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

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)

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?

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?

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

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?

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

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?

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

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

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 = 2.5 5 2+ 2.5 8 2 = 40 m 2 o Assume ceiling is ceiling ile Area = 5 8 = 40 m 2 i v d,i A i V

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 125.0 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=0 17.5 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

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 1 + 340 m3 /h 100m 3 m ppb 5.5 2 + 1.8 10 5.6 ppb + h 40m2 + 5.5 m h 40m2 + 5.5 m h 40m2 h 100m 3 α = 5.4 h -1

Indoor O3 concenraion (ppb) Simulaion resuls 2.5 2 1.5 1 0.5 0 C-O3: AC, Hom, He C-O3: Hom, He C-O3: He 0 1 2 3 4 5 6 7 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?

Inerpreing resuls Compue he conribuion o he overall sink from each erm in he equaion for alpha: Conribuion o indoor sinks 3.53 h -1 0.10 h -1 1.25 h -1 0.50 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).

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

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!

Simplifying assumpions Wha are wo unrealisic assumpions we have made in performing his analysis? Oudoor ozone is consan Island average O3 ( g/m 3 ) 100 80 60 40 20 0 Ozone Temperaure 0 5 10 15 20 25 Hour of day (average from 1/1/2007-31/12/2007) Figure 1. Average O 3 and emp. in Singapore 31 30 29 28 27 26 25 Temperaure (deg C) Oudoor ozone highly variable, diurnal rend ha is relian on solar radiaion o iniiae he chemisry (phoochemisry)

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

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

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

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

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

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

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.

Indoor ozone concenraion (ppb) Discreizing soluions and error 2.5 2 1.5 1 Discree, Δ = 0.05 h Discree, Δ = 0.1 h Discree, Δ = 0.005 h Analyical 0.5 0 0 0.2 0.4 0.6 0.8 1 Time (h)