Dynamics of Heating and Cooling Loads: Models

Size: px
Start display at page:

Download "Dynamics of Heating and Cooling Loads: Models"

Transcription

1 Dynamics of Heating and Cooling Loads: Models An Overview of Modelling of Thermostatically Controlled Load Presentation at Smart Grid Journal Club Xinbo Geng Texas A&M University April Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

2 Overview 1 Motivation 2 Models Start Point Deterministic Differential Equation Stochastic Difference Equation Markov Chain Matrix Equation Hybrid Partial Differential Equation 3 Relationships Among Models & Summary Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

3 Motivation Motivation Thermostatically Controlled Load (TCL): air conditioners (ACs), heat pumps, water heaters, and refrigerators, etc. In this presentation, TCL ACs In 2013, 40% of total U.S. energy consumption was consumed in residential and commercial buildings. Residential Thermostatically Controlled Loads (TCLs) such as represent about 20% of the total electricity consumption in the United States. Features of TCL load: installed storage like Energy Storage (batteries) inertia change the temperature set point, the effects will NOT be shown immediately controllable direct/indirect control Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

4 Motivation Need For Models... Models from the Department of Civil Engineering: objective: analyze the indoor environment too many parameters, too complicated too accurate Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

5 Motivation Wanted: Models Need For Models concise (less parameters) output (power) state variables (indoor temperature) and control Need For Models of Aggregated TCLs 1 TCL: 5kW 1 Generator: 1MW only aggregated TCLs matter models for a large population of TCLs Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

6 Models 1 Motivation 2 Models Start Point Deterministic Differential Equation Stochastic Difference Equation Markov Chain Matrix Equation Hybrid Partial Differential Equation 3 Relationships Among Models & Summary Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

7 Control Strategies in Common Dead-band Control (example: heating load) Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

8 Start Point Start Point Newton s Heating/Cooling Law dt dt = k(t T a) (1) reference any high school textbooks t time T temperature of an object T a temperature of its surroundings k proportionality constant. k 1/t. 1/k is a time constant. Verified by numerous experiments. Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

9 Deterministic Differential Equation Model#1: Deterministic Differential Equation dt dt = 1 τ (T T a WT g ) (2) reference Ihara, S., & Schweppe, F. C. (1981). Physically Based Modeling of Cold Load Pickup. IEEE Transactions on Power Apparatus and Systems t time T spatial average temperature of the interior of the house T a T g τ effective thermal time constant ambient temperature temperature gain of heater/cooler, proportional with the AC Power W binary variable. W = 0 AC is OFF; W = 1 AC is ON. Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

10 Stochastic Difference Equation Model#2: Stochastic Difference Equation T [(n+1) t] = e t τ T [n t] (1 e t τ )(Ta W [n t]t g )+V [n t] (3) reference: Mortensen, R. E., & Haggerty, K. P. (1988). A stochastic computer model for heating and cooling loads. IEEE Transactions on Power Systems t time T spatial average temperature of the interior of the house V ( ) error, sequence of independent Gaussian random variables of mean zero and variance τ effective thermal time constant ambient temperature T a T g temperature gain of heater/cooler W binary variable (indicating ON/OFF). n number of the sampling instant Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

11 Markov Chain Matrix Equation Model#3: Markov Chain Matrix Equation Reference: Mortensen, R. E., & Haggerty, K. P. (1988). A stochastic computer model for heating and cooling loads. IEEE Transactions on Power Systems t ON the time of cooling the room from T s + to T s t OFF the time of heating the room from T s to T s + Assumption The t ON satisfies a negative binomial distribution: P(t ON = n) = Cm 1 n 1 (1 p) n m if n m (4) P(t ON = n) = 0 if n < m (5) Estimation Run simulation using Model#2 and get a distribution of t ON, and find the best parameter p ON, m ON to fit the distribution Estimation similarly, we can estimate p OFF, m OFF. Markov Chain then we can formulate a Markov Chain to generate t ON and t OFF Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

12 Markov Chain Matrix Equation Model#3: Markov Chain Matrix Equation Criticism: Estimation good estimation for a long period, but could be very in accurate in a short time. Power States of TCLs Failed to reveal the relationship between P(t) and other variables. Not Very Helpful we cannot formulate control strategies based on this. An Alternative Approach Quantizing the indoor temperature T into discrete level (of Model #2), and use this as the states of a Markov Chain. Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

13 Markov Chain Matrix Equation A very special case: p i,i+1 = 1. Reference: Lu, N., & Chassin, D. P. (2004). A state-queueing model of thermostatically controlled appliances. IEEE Transactions on Power Systems Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19 Model#3.1: A State Queue Model

14 Hybrid Partial Differential Equation Model of Aggregated TCLs 2 Possible Approaches: Track Individual TCL Method have numerous copies of Model#1 or Model#2, add the power together Pros Accurate, (maybe) deterministic Cons privacy issue, (not necessary) to be exactly accurate Model the Distribution Overall behavior ignore many details Estimation Good Estimation when Large Population (Large Number s Law). Control partial differential equation control Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

15 Hybrid Partial Differential Equation Model#4: Hybrid Partial Differential Equation (CPFE Model) Reference: Malhame, R., & Chong, C.-Y. (1985). Electric load model synthesis by diffusion approximation of a high-order hybrid-state stochastic system. IEEE Transactions on Automatic Control Estimation Model TCLs as particles. Estimate by a diffusion process. Fokker-Planck Equation brief summary: describes the time evolution of the probability density function of the velocity of a particle under the influence of drag forces and random forces. check (forward) Chapman-Kolmogorov Equation and Fokker-Planck Equation for more details. 2 Process at the same time 1 process of indoor temperature (continuous) 2 process of state of TCLs (discrete/binary) Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

16 Hybrid Partial Differential Equation Model#4: Hybrid Partial Differential Equation (Details) Coupled Fokker - Planck Equations Model (CFPE Model) f 1 t f 0 t = T [f 1 τ (T (t) T a + T g )] + σ2 2 = T [f 0 τ (T (t) T a)] + σ2 2 2 f 1 T 2 (6) 2 f 0 T 2 (7) f 1 (T, t)dt = P(T < x(t) T + dt AC is ON at time t) (8) f 0 (T, t)dt = P(T < x(t) T + dt AC is OFF at time t) (9) And 6 more equations about boundary conditions: Absorbing Boundaries, Conditions at Infinity, Continuity Conditions, Probability Conservation Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

17 Hybrid Partial Differential Equation Model#4: Criticism not solvable Analytical closed form solutions to the Fokker-Planck partial differential equation are difficult to obtain except in particular cases. discretize to solve these equations numerically, it is necessary first to discretize them both in time and space. the transition to an implementable computer program. drawbacks of discretize whatever virtue the sophisticated continuous-state continuous-time mathematical formalism may possess, some of it is lost in the process of discretizing Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

18 Relationships Among Models & Summary Relationships Among Models & Summary Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

19 Relationships Among Models & Summary Any Questions? Dynamics of Heating and Cooling Loads: Models An Overview of Modelling of Thermostatically Controlled Load Presentation at Smart Grid Journal Club Xinbo Geng Texas A&M University April Xinbo Geng (TAMU) Dynamics of Heating and Cooling Loads: Models April / 19

Set Point Control of a Thermal Load Population

Set Point Control of a Thermal Load Population Set Point Control of a Thermal Load Population 1/15 Set Point Control of a Thermal Load Population CE 291F Final Project Travis Walter Civil Systems Program Civil and Environmental Engineering Department

More information

MODELING AND CONTROL OF THERMOSTATICALLY CONTROLLED LOADS

MODELING AND CONTROL OF THERMOSTATICALLY CONTROLLED LOADS MODELING AND CONTROL OF THERMOSTATICALLY CONTROLLED LOADS Soumya Kundu Nikolai Sinitsyn Scott Backhaus University of Michigan Los Alamos National Laboratory Los Alamos National Laboratory Ann Arbor USA

More information

Some Observations on the Locational Marginal Prices and Loads

Some Observations on the Locational Marginal Prices and Loads on the Locational Marginal Prices and Loads Xinbo Geng Texas A&M University gengxinbo@gmail.com October 15, 2014 Xinbo Geng (TAMU) Relationship between LMP and Load October 15, 2014 1 / 17 Overview 1 Preliminary

More information

Reduced-Order Modeling of Aggregated Thermostatic Loads With Demand Response

Reduced-Order Modeling of Aggregated Thermostatic Loads With Demand Response Reduced-Order Modeling of Aggregated Thermostatic Loads With Demand Response Wei Zhang, Jianming Lian, Chin-Yao Chang, Karanjit Kalsi and Yannan Sun Abstract Demand Response is playing an increasingly

More information

Large population optimal demand response for thermostatically controlled inertial loads

Large population optimal demand response for thermostatically controlled inertial loads Large population optimal demand response for thermostatically controlled inertial loads Gaurav Sharma, Le Xie and P. R. Kumar Department of Electrical and Computer Engineering, Texas A & M University.

More information

2014 IEEE. Digital Object Identifier:

2014 IEEE. Digital Object Identifier: 214 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other users, including reprinting/ republishing this material for advertising or promotional purposes,

More information

2.5 Linear Applications 99 Two-Tank Mixing. Two tanks A and B are assumed to contain A 0 and B 0 liters of brine at t = 0. Let the input for the rst t

2.5 Linear Applications 99 Two-Tank Mixing. Two tanks A and B are assumed to contain A 0 and B 0 liters of brine at t = 0. Let the input for the rst t 98 First Order Dierential Equations 2.5 Linear Applications This collection of applications for the linear equation y 0 + p(x)y = r(x) includes mixing problems, especially brine tanks in single and multiple

More information

Kalman filtering and friends: Inference in time series models. Herke van Hoof slides mostly by Michael Rubinstein

Kalman filtering and friends: Inference in time series models. Herke van Hoof slides mostly by Michael Rubinstein Kalman filtering and friends: Inference in time series models Herke van Hoof slides mostly by Michael Rubinstein Problem overview Goal Estimate most probable state at time k using measurement up to time

More information

2.5 Linear Applications

2.5 Linear Applications 2.5 Linear Applications 111 2.5 Linear Applications This collection of applications for the linear equation y + p(x)y = r(x) includes mixing problems, especially brine tanks in single and multiple cascade,

More information

Lecture 6: Bayesian Inference in SDE Models

Lecture 6: Bayesian Inference in SDE Models Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs

More information

Distributed demand control in power grids

Distributed demand control in power grids Distributed demand control in power grids and ODEs for Markov decision processes Second Conference on the Mathematics of Energy Markets Wolfgang Pauli Institute, Vienna, 4-9 July 217 Ana Bušić Inria Paris

More information

Problems 5: Continuous Markov process and the diffusion equation

Problems 5: Continuous Markov process and the diffusion equation Problems 5: Continuous Markov process and the diffusion equation Roman Belavkin Middlesex University Question Give a definition of Markov stochastic process. What is a continuous Markov process? Answer:

More information

MODELING HETEROGENEOUS POPULATIONS OF THERMOSTATICALLY CONTROLLED LOADS USING DIFFUSION-ADVECTION PDES

MODELING HETEROGENEOUS POPULATIONS OF THERMOSTATICALLY CONTROLLED LOADS USING DIFFUSION-ADVECTION PDES MODELING HETEROGENEOUS POPULATIONS OF THERMOSTATICALLY CONTROLLED LOADS USING DIFFUSION-ADVECTION PDES Scott Moura Victor Ruiz Dept. of Mechanical and Aerospace Engineering University of California, San

More information

Queueing Theory I Summary! Little s Law! Queueing System Notation! Stationary Analysis of Elementary Queueing Systems " M/M/1 " M/M/m " M/M/1/K "

Queueing Theory I Summary! Little s Law! Queueing System Notation! Stationary Analysis of Elementary Queueing Systems  M/M/1  M/M/m  M/M/1/K Queueing Theory I Summary Little s Law Queueing System Notation Stationary Analysis of Elementary Queueing Systems " M/M/1 " M/M/m " M/M/1/K " Little s Law a(t): the process that counts the number of arrivals

More information

Distributed demand side management via smart appliances contributing to frequency control

Distributed demand side management via smart appliances contributing to frequency control Journal of Chongqing University (English Edition) [ISSN 1671-84] Vol. 14 No. 3 September 015 doi:10.11835/j.issn.1671-84.015.03.03 o cite this article: ZHANG Wei-chen. via smart appliances contributing

More information

SECTION 1.2. DYNAMIC MODELS

SECTION 1.2. DYNAMIC MODELS CHAPTER 1 BY RADU MURESAN Page 1 ENGG4420 LECTURE 5 September 16 10 6:47 PM SECTION 1.2. DYNAMIC MODELS A dynamic model is a mathematical description of the process to be controlled. Specifically, a set

More information

Elementary Applications of Probability Theory

Elementary Applications of Probability Theory Elementary Applications of Probability Theory With an introduction to stochastic differential equations Second edition Henry C. Tuckwell Senior Research Fellow Stochastic Analysis Group of the Centre for

More information

ECEN 667 Power System Stability Lecture 25: FFT, Energy Methods

ECEN 667 Power System Stability Lecture 25: FFT, Energy Methods ECEN 667 Power System Stability Lecture 5: FFT, Energy Methods Prof. Tom Overbye Dept. of Electrical and Computer Engineering Texas A&M University, overbye@tamu.edu 1 Announcements Read Chapter 9 Final

More information

State-space Modelling of Hysteresis-based Control Schemes

State-space Modelling of Hysteresis-based Control Schemes European Control Conference (ECC) July 7-9,, Zürich, Switzerland. State-space Modelling of Hysteresis-based Control Schemes Soumya Kundu Ian A. Hiskens Abstract The paper deelops a state-space model for

More information

Markov decision processes with threshold-based piecewise-linear optimal policies

Markov decision processes with threshold-based piecewise-linear optimal policies 1/31 Markov decision processes with threshold-based piecewise-linear optimal policies T. Erseghe, A. Zanella, C. Codemo Dept. of Information Engineering, University of Padova, Italy Padova, June 2, 213

More information

Lecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations

Lecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations Lecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations Department of Biomedical Engineering and Computational Science Aalto University April 28, 2010 Contents 1 Multiple Model

More information

Demand Response Load Following of Source and Load Systems Hu, Jianqiang; Cao, Jinde; Yong, Taiyou; Guerrero, Josep M.; Z. Q. Chen, Michael; Li, Yaping

Demand Response Load Following of Source and Load Systems Hu, Jianqiang; Cao, Jinde; Yong, Taiyou; Guerrero, Josep M.; Z. Q. Chen, Michael; Li, Yaping Aalborg Universitet Demand Response Load Following of Source and Load Systems Hu, Jianqiang; Cao, Jinde; Yong, Taiyou; Guerrero, Josep M.; Z. Q. Chen, Michael; Li, Yaping Published in: I E E E Transactions

More information

Predicting the Electricity Demand Response via Data-driven Inverse Optimization

Predicting the Electricity Demand Response via Data-driven Inverse Optimization Predicting the Electricity Demand Response via Data-driven Inverse Optimization Workshop on Demand Response and Energy Storage Modeling Zagreb, Croatia Juan M. Morales 1 1 Department of Applied Mathematics,

More information

Real-Time Demand Response with Uncertain Renewable Energy in Smart Grid

Real-Time Demand Response with Uncertain Renewable Energy in Smart Grid Forty-Ninth Annual Allerton Conference Allerton House, UIUC, Illinois, USA September 28-3, 211 Real-Time Demand Response with Uncertain Renewable Energy in Smart Grid Libin Jiang and Steven Low Engineering

More information

Derivation of Itô SDE and Relationship to ODE and CTMC Models

Derivation of Itô SDE and Relationship to ODE and CTMC Models Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)

More information

Markovian Models for Electrical Load Prediction in Smart Buildings

Markovian Models for Electrical Load Prediction in Smart Buildings Markovian Models for Electrical Load Prediction in Smart Buildings Muhammad Kumail Haider, Asad Khalid Ismail, and Ihsan Ayyub Qazi LUMS School of Science and Engineering, Lahore, Pakistan {8,,ihsan.qazi}@lums.edu.pk

More information

Modeling Options for Demand Side Participation of Thermostatically Controlled Loads

Modeling Options for Demand Side Participation of Thermostatically Controlled Loads 213 IREP Symposium-Bulk Power System Dynamics and Control -IX (IREP), August 25-3, 213, Rethymnon, Greece Modeling Options for Demand Side Participation of Thermostatically Controlled Loads Maryam Kamgarpour,

More information

Table of Contents [ntc]

Table of Contents [ntc] Table of Contents [ntc] 1. Introduction: Contents and Maps Table of contents [ntc] Equilibrium thermodynamics overview [nln6] Thermal equilibrium and nonequilibrium [nln1] Levels of description in statistical

More information

Power Management of Cooling Systems with Dynamic Pricing

Power Management of Cooling Systems with Dynamic Pricing Power Management of Cooling Systems with Dynamic Pricing Saghar Hosseini, Ran Dai, and Mehran Mesbahi Abstract This paper addresses the optimal power management problems in electric cooling systems based

More information

STA 624 Practice Exam 2 Applied Stochastic Processes Spring, 2008

STA 624 Practice Exam 2 Applied Stochastic Processes Spring, 2008 Name STA 624 Practice Exam 2 Applied Stochastic Processes Spring, 2008 There are five questions on this test. DO use calculators if you need them. And then a miracle occurs is not a valid answer. There

More information

Design of IP networks with Quality of Service

Design of IP networks with Quality of Service Course of Multimedia Internet (Sub-course Reti Internet Multimediali ), AA 2010-2011 Prof. Pag. 1 Design of IP networks with Quality of Service 1 Course of Multimedia Internet (Sub-course Reti Internet

More information

Current Lead Optimization for Cryogenic Operation at Intermediate Temperatures. L. Bromberg, P.C. Michael, J.V. Minervini and C.

Current Lead Optimization for Cryogenic Operation at Intermediate Temperatures. L. Bromberg, P.C. Michael, J.V. Minervini and C. PSFC/JA-09-23 Current Lead Optimization for Cryogenic Operation at Intermediate Temperatures L. Bromberg, P.C. Michael, J.V. Minervini and C. Miles MIT Plasma Science and Fusion Center September 22, 2009

More information

Electrical Energy Modeling In Y2E2 Building Based On Distributed Sensors Information

Electrical Energy Modeling In Y2E2 Building Based On Distributed Sensors Information Electrical Energy Modeling In Y2E2 Building Based On Distributed Sensors Information Mahmoud Saadat Saman Ghili Introduction Close to 40% of the primary energy consumption in the U.S. comes from commercial

More information

Existence, Uniqueness and Stability of Invariant Distributions in Continuous-Time Stochastic Models

Existence, Uniqueness and Stability of Invariant Distributions in Continuous-Time Stochastic Models Existence, Uniqueness and Stability of Invariant Distributions in Continuous-Time Stochastic Models Christian Bayer and Klaus Wälde Weierstrass Institute for Applied Analysis and Stochastics and University

More information

International Conference on Mechanics, Materials and Structural Engineering (ICMMSE 2016)

International Conference on Mechanics, Materials and Structural Engineering (ICMMSE 2016) International Conference on Mechanics, Materials and Structural Engineering (ICMMSE 2016) hermal Model Parameter Estimation for HVAC Facility using Recursive Least Square Method Dae Ki Kim1, a, Kyu Chul

More information

Morgan Jones, Matthew M. Peet

Morgan Jones, Matthew M. Peet Solving Dynamic Programming with Supremum Terms in the Objective and Application to Optimal Battery Scheduling for Electricity Consumers Subject to Demand Charges Morgan Jones, Matthew M. Peet Abstract

More information

40P (2 x 60 x 60) = 2.5 x 10 6 (4200)(5) P = 1.82 x 10 5 W

40P (2 x 60 x 60) = 2.5 x 10 6 (4200)(5) P = 1.82 x 10 5 W NAME : F.3C ( ) Marks: /50 Form 3 Physics Assessment on Heat Time allowed: 45 minutes Section A (34 marks) 1. An indoor swimming pool containing 2.5 x 10 6 kg of water uses 40 identical heaters to maintain

More information

Control of Sampled Switched Systems using Invariance Analysis

Control of Sampled Switched Systems using Invariance Analysis 1st French Singaporean Workshop on Formal Methods and Applications Control of Sampled Switched Systems using Invariance Analysis Laurent Fribourg LSV - ENS Cachan & CNRS Laurent Fribourg Lsv - ENS Cachan

More information

Estimation and Control of Quality of Service in Demand Dispatch

Estimation and Control of Quality of Service in Demand Dispatch Estimation and Control of Quality of Service in Demand Dispatch Yue Chen, Ana Bušić, and Sean Meyn August 2, 2018 arxiv:1609.00051v1 [cs.sy] 31 Aug 2016 Abstract It is now well known that flexibility of

More information

British Physics Olympiad 2006

British Physics Olympiad 2006 1 British Physics Olympiad 2006 Paper 3 Monday 27 th February 2006. Time allowed 3hrs plus 15 minutes reading time. All questions should be attempted. All questions carry 20 marks each. Question 4 of the

More information

Autoregressive Linear Thermal Model of a Residential Forced-Air Heating System with Backpropagation Parameter Estimation Algorithm

Autoregressive Linear Thermal Model of a Residential Forced-Air Heating System with Backpropagation Parameter Estimation Algorithm Autoregressive Linear Thermal Model of a Residential Forced-Air Heating System with Backpropagation Parameter Estimation Algorithm Eric Burger Scott Moura David E. Culler Electrical Engineering and Computer

More information

Coordinated Aggregation of Distributed Resources

Coordinated Aggregation of Distributed Resources Coordinated Aggregation of Distributed Resources UC Berkeley October 13, 2011 Coordinated Aggregation 1 of 39 Co-conspirators Anand Subramanian, Manuel Garcia, Josh Taylor [Berkeley Students] Duncan Callaway,

More information

ECE 422/522 Power System Operations & Planning/ Power Systems Analysis II 3 Load Modeling

ECE 422/522 Power System Operations & Planning/ Power Systems Analysis II 3 Load Modeling ECE 422/522 Power System Operations & Planning/ Power Systems Analysis II 3 Load Modeling Spring 2014 Instructor: Kai Sun 1 References 1. Load Performance for Dynamic Performance Analysis, IEEE Committee

More information

Dynamic Response Analysis of Start-up Transient in Air Conditioning System

Dynamic Response Analysis of Start-up Transient in Air Conditioning System Dynamic Response Analysis of Start-up Transient in Air Conditioning System รองศาสตราจารย ดร.ประกอบ ส รว ฒนาวรรณ 1 และ ส รส ทธ พลพ นธ 2 1รองคณบด คณะว ศวกรรมศาสตร มหาว ทยาล ยเกษตรศาสตร กร งเทพฯ E-mail: fengpsw@ku.ac.th

More information

Hybrid Control and Switched Systems. Lecture #1 Hybrid systems are everywhere: Examples

Hybrid Control and Switched Systems. Lecture #1 Hybrid systems are everywhere: Examples Hybrid Control and Switched Systems Lecture #1 Hybrid systems are everywhere: Examples João P. Hespanha University of California at Santa Barbara Summary Examples of hybrid systems 1. Bouncing ball 2.

More information

and Joachim Peinke ForWind Center for Wind Energy Research Institute of Physics, Carl-von-Ossietzky University Oldenburg Oldenburg, Germany

and Joachim Peinke ForWind Center for Wind Energy Research Institute of Physics, Carl-von-Ossietzky University Oldenburg Oldenburg, Germany Stochastic data analysis for in-situ damage analysis Philip Rinn, 1, a) Hendrik Heißelmann, 1 Matthias Wächter, 1 1, b) and Joachim Peinke ForWind Center for Wind Energy Research Institute of Physics,

More information

Household Energy Disaggregation based on Difference Hidden Markov Model

Household Energy Disaggregation based on Difference Hidden Markov Model 1 Household Energy Disaggregation based on Difference Hidden Markov Model Ali Hemmatifar (ahemmati@stanford.edu) Manohar Mogadali (manoharm@stanford.edu) Abstract We aim to address the problem of energy

More information

M4A42 APPLIED STOCHASTIC PROCESSES

M4A42 APPLIED STOCHASTIC PROCESSES M4A42 APPLIED STOCHASTIC PROCESSES G.A. Pavliotis Department of Mathematics Imperial College London, UK LECTURE 1 12/10/2009 Lectures: Mondays 09:00-11:00, Huxley 139, Tuesdays 09:00-10:00, Huxley 144.

More information

A Learning Agent for Heat-Pump Thermostat Control

A Learning Agent for Heat-Pump Thermostat Control A Learning Agent for Heat-Pump Thermostat Control Daniel Urieli and Peter Stone Department of Computer Science The University of Texas at Austin {urieli,pstone}@cs.utexas.edu Heating, Ventilation, and

More information

arxiv: v1 [cs.sy] 27 Jan 2013

arxiv: v1 [cs.sy] 27 Jan 2013 A Two Level Feedback System Design to Regulation Service Provision Bowen Zhang 1 and John Baillieul 2 arxiv:1316386v1 [cssy] 27 Jan 213 Abstract Demand side management has gained increasing importance

More information

(implicitly assuming time-homogeneity from here on)

(implicitly assuming time-homogeneity from here on) Continuous-Time Markov Chains Models Tuesday, November 15, 2011 2:02 PM The fundamental object describing the dynamics of a CTMC (continuous-time Markov chain) is the probability transition (matrix) function:

More information

Stochastic process. X, a series of random variables indexed by t

Stochastic process. X, a series of random variables indexed by t Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,

More information

Recap. Probability, stochastic processes, Markov chains. ELEC-C7210 Modeling and analysis of communication networks

Recap. Probability, stochastic processes, Markov chains. ELEC-C7210 Modeling and analysis of communication networks Recap Probability, stochastic processes, Markov chains ELEC-C7210 Modeling and analysis of communication networks 1 Recap: Probability theory important distributions Discrete distributions Geometric distribution

More information

Different approaches to model wind speed based on stochastic differential equations

Different approaches to model wind speed based on stochastic differential equations 1 Universidad de Castilla La Mancha Different approaches to model wind speed based on stochastic differential equations Rafael Zárate-Miñano Escuela de Ingeniería Minera e Industrial de Almadén Universidad

More information

Gillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde

Gillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation

More information

Exact Sizing of Battery Capacity for Photovoltaic Systems

Exact Sizing of Battery Capacity for Photovoltaic Systems Exact Sizing of Battery Capacity for Photovoltaic Systems Yu Ru a, Jan Kleissl b, Sonia Martinez b a GE Global Research at Shanghai (email: Yu.Ru@ge.com). b Mechanical and Aerospace Engineering Department,

More information

Systems Driven by Alpha-Stable Noises

Systems Driven by Alpha-Stable Noises Engineering Mechanics:A Force for the 21 st Century Proceedings of the 12 th Engineering Mechanics Conference La Jolla, California, May 17-20, 1998 H. Murakami and J. E. Luco (Editors) @ASCE, Reston, VA,

More information

B6VMAX Series TECHNICAL SPECIFICATIONS. FEATURES and BENEFITS

B6VMAX Series TECHNICAL SPECIFICATIONS. FEATURES and BENEFITS TECHNICAL SPECIFICATIONS B6VMAX Series Air Handler with TXV 13-19 SEER Residential System 18,000-60,000 Btuh (Heat Pump & Air Conditioner) R-410A Refrigerant The B6VMAX Series of air handlers, when combined

More information

Dynamic simulation of DH house stations

Dynamic simulation of DH house stations Article Dynamic simulation of DH house stations Jan Eric Thorsen, Director, DHS Application Centre, Danfoss A/S www.danfoss.com Jan Eric Thorsen, Director, DHS Application Centre, Danfoss A/S Presented

More information

Synchrony in Stochastic Pulse-coupled Neuronal Network Models

Synchrony in Stochastic Pulse-coupled Neuronal Network Models Synchrony in Stochastic Pulse-coupled Neuronal Network Models Katie Newhall Gregor Kovačič and Peter Kramer Aaditya Rangan and David Cai 2 Rensselaer Polytechnic Institute, Troy, New York 2 Courant Institute,

More information

Unit D-2: List of Subjects

Unit D-2: List of Subjects ES312 Energy Transfer Fundamentals Unit D: The Second Law of Thermodynamics ROAD MAP... D-1: From First to Second Law D-2: Second Law of Thermodynamics Unit D-2: List of Subjects Second Law of Thermodynamics

More information

1 Introduction. 2 Diffusion equation and central limit theorem. The content of these notes is also covered by chapter 3 section B of [1].

1 Introduction. 2 Diffusion equation and central limit theorem. The content of these notes is also covered by chapter 3 section B of [1]. 1 Introduction The content of these notes is also covered by chapter 3 section B of [1]. Diffusion equation and central limit theorem Consider a sequence {ξ i } i=1 i.i.d. ξ i = d ξ with ξ : Ω { Dx, 0,

More information

STOCHASTIC PROCESSES Basic notions

STOCHASTIC PROCESSES Basic notions J. Virtamo 38.3143 Queueing Theory / Stochastic processes 1 STOCHASTIC PROCESSES Basic notions Often the systems we consider evolve in time and we are interested in their dynamic behaviour, usually involving

More information

PROBLEM 1.2 ( ) 25 C 15 C dx L 0.30 m Ambient air temperature, T2 (C)

PROBLEM 1.2 ( ) 25 C 15 C dx L 0.30 m Ambient air temperature, T2 (C) PROBLEM 1.2 KNOWN: Inner surface temperature and thermal conductivity of a concrete wall. FIND: Heat loss by conduction through the wall as a function of ambient air temperatures ranging from -15 to 38

More information

IDEA Campus Energy 2016

IDEA Campus Energy 2016 IDEA Campus Energy 2016 Sufficient metering - Best practice for thermal energy metering www.aquametro.com IDEA Campus Energy 2016 Sufficient metering - Best practice for thermal energy metering Overview

More information

Stochastic Modelling Unit 1: Markov chain models

Stochastic Modelling Unit 1: Markov chain models Stochastic Modelling Unit 1: Markov chain models Russell Gerrard and Douglas Wright Cass Business School, City University, London June 2004 Contents of Unit 1 1 Stochastic Processes 2 Markov Chains 3 Poisson

More information

The Wright-Fisher Model and Genetic Drift

The Wright-Fisher Model and Genetic Drift The Wright-Fisher Model and Genetic Drift January 22, 2015 1 1 Hardy-Weinberg Equilibrium Our goal is to understand the dynamics of allele and genotype frequencies in an infinite, randomlymating population

More information

A Unified Stochastic Hybrid System Approach to Aggregate Modeling of Responsive Loads

A Unified Stochastic Hybrid System Approach to Aggregate Modeling of Responsive Loads 1 A Unified Stochastic Hybrid System Approach to Aggregate Modeling of Responsive Loads Lin Zhao and Wei Zhang arxiv:153.6911v3 [cs.sy] 18 Aug 217 Abstract Aggregate load modeling is of fundamental importance

More information

Downloaded from

Downloaded from CHAPTER 12 ELECTRICITY Electricity is a general term that encompasses a variety of phenomena resulting from the presence and flow of electric charge. These include many easily recognizable phenomena such

More information

METHODS OF HEAT COST ALLOCATION

METHODS OF HEAT COST ALLOCATION PERIODICA POLYTECHNICA SER. MECH. ENG. VOL. 44, NO. 2, PP. 227 236 (2000) METHODS OF HEAT COST ALLOCATION István CSOKNYAI Department of Building Service Engineering Budapest University of Technology and

More information

Population Genetics: a tutorial

Population Genetics: a tutorial : a tutorial Institute for Science and Technology Austria ThRaSh 2014 provides the basic mathematical foundation of evolutionary theory allows a better understanding of experiments allows the development

More information

A Bayesian method for the analysis of deterministic and stochastic time series

A Bayesian method for the analysis of deterministic and stochastic time series A Bayesian method for the analysis of deterministic and stochastic time series Coryn Bailer-Jones Max Planck Institute for Astronomy, Heidelberg DPG, Berlin, March 2015 Time series modelling heteroscedastic,

More information

Lecture 2: From Linear Regression to Kalman Filter and Beyond

Lecture 2: From Linear Regression to Kalman Filter and Beyond Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing

More information

Irreducibility. Irreducible. every state can be reached from every other state For any i,j, exist an m 0, such that. Absorbing state: p jj =1

Irreducibility. Irreducible. every state can be reached from every other state For any i,j, exist an m 0, such that. Absorbing state: p jj =1 Irreducibility Irreducible every state can be reached from every other state For any i,j, exist an m 0, such that i,j are communicate, if the above condition is valid Irreducible: all states are communicate

More information

Lecture 1: Brief Review on Stochastic Processes

Lecture 1: Brief Review on Stochastic Processes Lecture 1: Brief Review on Stochastic Processes A stochastic process is a collection of random variables {X t (s) : t T, s S}, where T is some index set and S is the common sample space of the random variables.

More information

16. Working with the Langevin and Fokker-Planck equations

16. Working with the Langevin and Fokker-Planck equations 16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation

More information

NEGST. New generation of solar thermal systems. Advanced applications ENEA. Comparison of solar cooling technologies. Vincenzo Sabatelli

NEGST. New generation of solar thermal systems. Advanced applications ENEA. Comparison of solar cooling technologies. Vincenzo Sabatelli NEGST New generation of solar thermal systems Advanced applications Comparison of solar cooling technologies Vincenzo Sabatelli ENEA vincenzo.sabatelli@trisaia.enea.it NEGST Workshop - Freiburg - June

More information

Nonlinear State Estimation! Extended Kalman Filters!

Nonlinear State Estimation! Extended Kalman Filters! Nonlinear State Estimation! Extended Kalman Filters! Robert Stengel! Optimal Control and Estimation, MAE 546! Princeton University, 2017!! Deformation of the probability distribution!! Neighboring-optimal

More information

LOAD POCKET MODELING. KEY WORDS Load Pocket Modeling, Load Forecasting

LOAD POCKET MODELING. KEY WORDS Load Pocket Modeling, Load Forecasting LOAD POCKET MODELING Eugene A. Feinberg, Dora Genethliou Department of Applied Mathematics and Statistics State University of New York at Stony Brook Stony Brook, NY 794-36, USA Janos T. Hajagos KeySpan

More information

Two mark questions and answers UNIT II SECOND LAW 1. Define Clausius statement. It is impossible for a self-acting machine working in a cyclic process, to transfer heat from a body at lower temperature

More information

Sloppy derivations of Ito s formula and the Fokker-Planck equations

Sloppy derivations of Ito s formula and the Fokker-Planck equations Sloppy derivations of Ito s formula and the Fokker-Planck equations P. G. Harrison Department of Computing, Imperial College London South Kensington Campus, London SW7 AZ, UK email: pgh@doc.ic.ac.uk April

More information

EMPIRICAL MODEL OF A 11 kw (NOMINAL COOLING) R134a WATER-WATER HEAT PUMP

EMPIRICAL MODEL OF A 11 kw (NOMINAL COOLING) R134a WATER-WATER HEAT PUMP EMPIRICAL MODEL OF A 11 kw (NOMINAL COOLING) R134a WATER-WATER HEAT PUMP Sébastien A. Brideau 1, Ian Beausoleil-Morrison 1, and Michaël Kummert 1 Department of Mechanical and Aerospace Engineering, Carleton

More information

Math 166: Topics in Contemporary Mathematics II

Math 166: Topics in Contemporary Mathematics II Math 166: Topics in Contemporary Mathematics II Xin Ma Texas A&M University November 26, 2017 Xin Ma (TAMU) Math 166 November 26, 2017 1 / 14 A Review A Markov process is a finite sequence of experiments

More information

Diffusion in the cell

Diffusion in the cell Diffusion in the cell Single particle (random walk) Microscopic view Macroscopic view Measuring diffusion Diffusion occurs via Brownian motion (passive) Ex.: D = 100 μm 2 /s for typical protein in water

More information

XHDU Cassette type fan-coils for double/four pipe connections 1,98-14 kw HYDRONIC TERMINALS

XHDU Cassette type fan-coils for double/four pipe connections 1,98-14 kw HYDRONIC TERMINALS HYDRONIC TERMINALS XHDU 402 1203 Versions XHDU 2T XHDU-R 2T XHDU 4T XHDE 2T Single coil for double pipe connections Single coil for double pipe connections with heater Double coil for four pipe connections

More information

arxiv: v1 [cs.sy] 27 Sep 2016

arxiv: v1 [cs.sy] 27 Sep 2016 Battery Capacity of Deferrable Energy Demand Daria Madjidian, Mardavij Roozbehani, and Munther A. Dahleh arxiv:6.973v [cs.sy] 7 Sep 6 Abstract We investigate the ability of a homogeneous collection of

More information

Stochastic processes and

Stochastic processes and Stochastic processes and Markov chains (part II) Wessel van Wieringen w.n.van.wieringen@vu.nl wieringen@vu nl Department of Epidemiology and Biostatistics, VUmc & Department of Mathematics, VU University

More information

Session 1: Probability and Markov chains

Session 1: Probability and Markov chains Session 1: Probability and Markov chains 1. Probability distributions and densities. 2. Relevant distributions. 3. Change of variable. 4. Stochastic processes. 5. The Markov property. 6. Markov finite

More information

Q = (c) Assuming that Ricoh has been working continuously for 7 days, what is the probability that it will remain working at least 8 more days?

Q = (c) Assuming that Ricoh has been working continuously for 7 days, what is the probability that it will remain working at least 8 more days? IEOR 4106: Introduction to Operations Research: Stochastic Models Spring 2005, Professor Whitt, Second Midterm Exam Chapters 5-6 in Ross, Thursday, March 31, 11:00am-1:00pm Open Book: but only the Ross

More information

HiLumi LHC FP7 High Luminosity Large Hadron Collider Design Study. Milestone Report. Cryogenic Scenarios for the Cold Powering System

HiLumi LHC FP7 High Luminosity Large Hadron Collider Design Study. Milestone Report. Cryogenic Scenarios for the Cold Powering System CERN-ACC-2014-0065 HiLumi LHC FP7 High Luminosity Large Hadron Collider Design Study Milestone Report Cryogenic Scenarios for the Cold Powering System Ballarino, A (CERN) et al 27 May 2014 The HiLumi LHC

More information

Assessment of Optimal Flexibility in Ensemble of Frequency Responsive Loads

Assessment of Optimal Flexibility in Ensemble of Frequency Responsive Loads Assessment of Optimal Flexibility in Ensemble of Frequency Responsive Loads Soumya Kundu, Jacob Hansen, Jianming Lian and Karan Kalsi Optimization and Control Group acific Northwest National Laboratory,

More information

Investigating the ability of various buildings in handling load shiftings

Investigating the ability of various buildings in handling load shiftings Investigating the ability of various buildings in handling load shiftings Paul Malisani,Bérenger Favre,Stéphane Thiers, Bruno Peuportier, François Chaplais, Nicolas Petit EDF R&D CEP - MINES-ParisTech

More information

EXAMPLES EXAMPLE - Temperature in building

EXAMPLES EXAMPLE - Temperature in building DYNAMICAL SYSTEMS EXAMPLES EXAMPLE - Temperature in building Energy balance: Rate of change = [Inflow of energy] [Outflow of energy] of stored energy where Rate of change of stored energy = cρv dt (c =

More information

CHAPTER 6 THE SECOND LAW OF THERMODYNAMICS

CHAPTER 6 THE SECOND LAW OF THERMODYNAMICS CHAPTER 6 THE SECOND LAW OF THERMODYNAMICS S. I. Abdel-Khalik (2014) 1 CHAPTER 6 -- The Second Law of Thermodynamics OUTCOME: Identify Valid (possible) Processes as those that satisfy both the first and

More information

Integral-based Algorithm for Parameter Identification of the Heat Exchanger

Integral-based Algorithm for Parameter Identification of the Heat Exchanger Proceedings of the International MultiConference of Engineers and Computer Scientists 218 Vol II IMECS 218, March 14-16, 218, Hong Kong Integral-based Algorithm for Parameter Identification of the Heat

More information

Water terminal control solutions

Water terminal control solutions Water Terminal Control Solutions Water terminal control solutions Carrier offers a wide range of fan coil units designed to meet the needs of different systems and applications: from cabinet units in the

More information

p 1 ( Y p dp) 1/p ( X p dp) 1 1 p

p 1 ( Y p dp) 1/p ( X p dp) 1 1 p Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p

More information

Regression Analysis Tutorial 77 LECTURE /DISCUSSION. Specification of the OLS Regression Model

Regression Analysis Tutorial 77 LECTURE /DISCUSSION. Specification of the OLS Regression Model Regression Analysis Tutorial 77 LECTURE /DISCUSSION Specification of the OLS Regression Model Regression Analysis Tutorial 78 The Specification The specification is the selection of explanatory variables

More information

Green Scheduling of Control Systems

Green Scheduling of Control Systems University of Pennsylvania ScholarlyCommons Publicly Accessible Penn Dissertations 1-1-2012 Green Scheduling of Control Systems Truong Xuan Nghiem University of Pennsylvania, truong.nghiem@gmail.com Follow

More information

Brownian Motion and Poisson Process

Brownian Motion and Poisson Process and Poisson Process She: What is white noise? He: It is the best model of a totally unpredictable process. She: Are you implying, I am white noise? He: No, it does not exist. Dialogue of an unknown couple.

More information