## How Can Smart Buildings Be Price-Responsive?

Ricardo Fern´ andez-Blanco Group OASYS University of Malaga Malaga, Spain Ricardo.FCarramolino@uma.es

## Juan Miguel Morales

Salvador Pineda

Dept. Applied Mathematics University of Malaga Malaga, Spain Juan.Morales@uma.es

Abstract -The prospective participation of smart buildings in the electricity system is strongly related to the increasing active role of demand-side resources in the electrical grid. In addition, the growing penetration of smart meters and recent advances on home automation technologies will spur the development of new mathematical tools to help optimize the local resources of these buildings. Within this context, this paper first provides a comprehensive model to determine the electrical consumption of a single-zone household based on economic model predictive control. The goal of this problem is to minimize the electricity consumption cost while accounting for the heating dynamics of the building, smart home appliances, and comfort constraints. This paper then identifies and analyzes the key parameters responsible for the price-responsive behaviour of smart households.

Index Terms -Model predictive control, price-responsive loads, smart appliances, smart buildings.

## NOMENCLATURE

A bar above ( x ) or below ( x ) any variable denotes its maximum and minimum bounds, respectively.

## A. Sets and Indices

- C i Set of power phases of load i , indexed by c .
- C Set of user-defined comfort constraints.
- I Set of uninterruptible loads, indexed by i .
- T Set of time periods, indexed by t .
- T i Set of desirable operation times of load i .

T occ Set of time periods during which the household is occupied, indexed by t .

- T Set of technical constraints pertaining to control actions.

## B. Parameters

- CT i Cycle duty of uninterruptible load i [# time periods]. N T Number of time periods.

Dept. Electrical Engineering University of Malaga Malaga, Spain spinedamorente@gmail.com

z amb t

Outdoor ambient temperature in period t [ ◦ C].

z

hwd

t

Hot water demand in period t [l].

z

occ

t

Occupancy schedule in period t .

φ

t

Luminance from natural light [lumen].

η

le

Indoor luminous efficacy [lumen/kW].

λ

t

Electricity price in period t [ e /kWh].

ρ

l

Penalty term for the light level bounds [ e /lumen].

ρ temp

Penalty term for temperature bounds [ e / ◦ C].

## C. Variables

l

t

Light level in period t [lumen].

u b t [kW].

Power consumption of the building in period t

u bl t Window blind position in period t [p.u.]. u ı t

Power consumption in period t . Superscript ı = { heat , cool , hp , wh , rf , al } and they denote heating and cooling system, heat pump for the floor heating system, water heater, refrigerator, and artificial lighting, respectively [kW].

u ul it

Power consumption of uninterruptible load i in period t [kW].

v l t

Slack variable for the light level bounds in period t [lumen].

v r/wh/rf t

Slack variable for the temperature bounds of the room/water heater/refrigerator in period t [ ◦ C].

x β t

Temperature in period t . Superscrit β = { r, f, w, wh, rf } and they denote temperature of room, floor, water in floor heating pipes, water heater, refrigeration chamber, respectively [ ◦ C].

δ a cit

Auxiliary variable to relate power phase c with the scheduling of uninterruptible load i in period t . Binary variable indicating whether uninterruptible load i is scheduled on (1) or off (0) in period t .

δ ul it

P

ul

Power of phase c of uninterruptible load i [kW].

Standby power consumption in period t [kW].

ci P sb t

TR it

Binary-valued parameter specifying the periods of operation of the uninterruptable load i .

r,a

- UA the ambient [W/ C].

Heat transfer coefficient between the room air and ◦

This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No 755705). This project has also been supported by Fundaci´ on Iberdrola Espa˜a 2018. The authors thankfully acknowledge the n computer resources, technical expertise and assistance provided by the SCBI (Supercomputing and Bioinformatics) center of the University of Malaga.

## D. Vectors

- u Vector of control variables.
- v Vector of slack variables.
- x Vector of state variables.
- y Vector of measured signals.
- z Vector of external disturbances.

## E. Matrices

- A Matrix of heat transfer coefficients (dynamics matrix). B Matrix of coefficients relating the state and control variables (control matrix).
- C Matrix of coefficients relating the state variables and the measured signals (sensor matrix).
- D Matrix of coefficients relating the measured signals directly with the control variables (direct term).
- E Matrix of coefficients relating the state variables and external disturbances (disturbances matrix).

## I. INTRODUCTION

Demand-side resources are expected to play an important role in future electricity systems [1]. Their active participation in providing flexibility to the grid will also help integrate more and more renewable energy, thus paving the way towards the decarbonization of the electricity sector. Prospective demandside resources are storage devices and smart buildings. Largescale integration of the former is expected to increase provided capital cost reductions and increased charging/discharging efficiencies [2], whereas the latter is gaining more attention due to recent advances on home automation and the growing penetration of smart meters in the distribution network [3]. This paper is focused on the smart buildings' ability to provide flexibility, which may be of utmost interest to incentivize the active participation of demand-side resources, thus leading to new ways of electricity trading [4].

There is a vast literature on building simulation by: (i) individually modelling thermal loads such as refrigerators [5], water heaters [6], or smart solar tanks [7]; or (ii) modelling the building heating dynamics [8]. Halvgaard et al. [8] analyzed a single-zone building with a water-based floor heater (FH) via economic model predictive control (MPC) in order to shift power to periods with lower prices. Note that an economic MPC is a model that represents complex system dynamics while minimizing an economic objective function.

The active participation of smart buildings in the electricity system has been widely investigated in the literature [9]-[15]. Most of the works consider thermostatically-controlled loads (TCLs) only [9]-[12]. Li et al. [9] devised a market-based coordination of a pool of TCLs which were modelled by a second-order equivalent thermal parameter (ETP) model. Authors of [10] investigated the effects of dynamic-price retail contracts on integrated retail and power market operations wherein the households are represented by a simplified ETP model of their air conditioning systems. Authors of [11] simulated a pool of price-responsive households equipped with heating pumps to apply inverse optimization in order to forecast their consumption. Zhao and Zhang [12] developed a framework to aggregate price-responsive loads represented by a second-order ETP model. However, little attention has been paid to an integrated formulation of the smart building with smart appliances [13]-[15], which is crucial to accurately capture the potential flexibility of their buildings.

Anjos et al. [14] proposed a Dantzig-Wolfe decomposition approach to deal with numerous and heterogeneous buildings managed by a central aggregator. Although, they model one of each load type that a building may have, the building heating dynamics were disregarded, thus ignoring the effect of the comfort constraints. Contreras-Oca˜a n et al. [13] presented one of the most completed building models and analyzed possible interactions between commercial buildings and an aggregator of electric vehicles. Gonzalez et al. [15] presented a residential load simulator with appliances. However, their focus was on providing load profiles rather than on how the household could provide flexibility. Recently, Junker et al. [16] devised a flexibility metric to characterize generic buildings (or districts) that can measure their reactions to different penalty signals.

The major contributions of this paper are thus twofold:

- 1) We propose a compact formulation for modelling smart buildings with smart appliances, and a comprehensive formulation for a single-zone household. It includes a state-space model with five states in order to capture the building heat dynamics, comfort constraints with user-defined parameters, and technical constraints. As major novelties, the model includes the implementation of two different space heaters (water-based floor heater or HVAC systems) along with other thermal loads (e.g. refrigerator and/or water heater); and the uninterruptible loads are modelled with discrete variables to accurately account for their variable cycle power.
- 2) We analyze the key drivers leading to price-responsive households based on their comfort settings and structural characteristics.

## II. COMPACT PROBLEM FORMULATION

Managing buildings with smart appliances call for an integrated operational model that takes into account the heating dynamics of the building. Economic MPC becomes thus a solid option to control the building and it has been applied to building dynamic simulation in the last decade (see [17] and references therein). This section provides a compact mathematical formulation for modelling a smart building by using economic MPC. The model can be mathematically expressed as the following optimization problem:

$$\min _ { u, v, x, y } \quad h ( u, v )$$

$$\text{subject to} \colon \ \dot { x } = A x + B u + E z \quad \quad ( 1 b )$$

$$y = C x + D u \text{ \quad \ \ } ( 1 c )$$

$$u, v, y \in \mathbb { C }$$

$$u \in \mathbb { T }, \quad \ \ \ \ ( 1 e )$$

where h ( ) · comprises the electricity cost incurred by the smart building and the penalty cost due to discomfort, which is minimized in constraint (1a).

Constraints (1b)-(1c) define the state-space model of the building heat dynamics (including thermal loads), which are represented by a linear system of first-order differential equations. External disturbances can be accounted for in this system. Specifically, constraints (1b) model the heat balance in each building zone and thermal load (e.g. refrigerator or water heater). Constraints (1c) account for the relationship among the measured signals with the state and control variables. Equation (1d) models the set of user-defined comfort constraints, which include control over the indoor air temperature and lighting, among others. Finally, equation (1e) represents the set of technical constraints including the scheduling of uninterruptible loads such as washing machine, tumble dryer, dishwasher, smart oven, and so on.

The system of differential equations can be discretized using Euler's method and thus the problem can be recast as:

$$\min _ { u, v, x, y } \ h ( u, v ) \quad \ ^ { - } & ( 2$$

subject to:

$$x _ { t } = A _ { d } x _ { t - 1 } + B _ { d } u _ { t -$$

$$y _ { t } = C _ { d } x _ { t - 1 } + D _ { d } u _ { t -$$

$$u, v, y \in \mathcal { C } & & ( 2 d ) & & \quad \quad \quad$$

$$u \in \mathbb { T }, & & ( 2 e ) = \stackrel { \mu \nu } { \lambda$$

where expressions (2a), (2d)-(2e) are identical to (1a), (1d)(1e), in that order, but expressed in discrete time. Similarly, constraints (2b)-(2c) comprise the discrete state-space model described in (1b)-(1c). Subscript d denotes that the matrices are in discrete time.

## III. PRICE-RESPONSIVE HOUSEHOLD FORMULATION

Given a single-zone household, we explore two options for space heating: 1) a water-based FH, and 2) an HVAC 1 system. Apart from these thermal loads, we model a residential refrigerator and a water heater, as well as the indoor lighting. Uninterruptible loads are also represented in the model. Assuming that the state variables are directly observed (i.e., C d and D d are, respectively, equal to the identity and null matrices of appropriate dimensions), a price-responsive household can be formulated as:

$$\min _ { \Xi } \sum _ { t \in \mathcal { T } } \lambda _ { t } u _ { t } ^ { b } + \sum _ { t \in \mathcal { T } } \left [ \rho ^ { t e m p } \left ( v _ { t } ^ { r } + v _ { t } ^ { w h } + v _ { t } ^ { r f } \right ) + \rho ^ { l } v _ { t } ^ { l } \right ] \ ( 3 ) \quad \text{the insi} \\ \dots$$

subject to:

$$& \text{subject to:} \\ & \text{subject to:} \\ & \underbrace { \left ( \begin{smallmatrix} x _ { t } ^ { r } \\ x _ { t } ^ { f } \\ x _ { t } ^ { f } \\ x _ { t } ^ { r f } \\ x _ { t } ^ { w h } \end{smallmatrix} } ^ { \infty } & = A _ { d } \left ( \begin{smallmatrix} x _ { t - 1 } ^ { r } \\ x _ { t - 1 } ^ { f } \\ x _ { t - 1 } ^ { f } \\ x _ { t - 1 } ^ { r f } \\ x _ { t - 1 } ^ { w h } \end{smallmatrix} } ^ { \infty } + B _ { d } \left ( \begin{smallmatrix} u _ { t - 1 } ^ { h p } \\ u _ { t - 1 } ^ { h e a t } \\ u _ { t - 1 } ^ { t } \\ u _ { t - 1 } ^ { b l } \\ u _ { t - 1 } ^ { r f } \\ u _ { t - 1 } ^ { w h - 1 } \end{smallmatrix} } ^ { \infty } + E _ { d } \left ( \begin{smallmatrix} z _ { t - 1 } ^ { a m b } \\ z _ { t - 1 } ^ { e d c - 1 } \\ z _ { t - 1 } ^ { w h d } \end{smallmatrix} \right ) } ^ { \text{In a} } \\ & \text{refr} \\ & \text{$\forall t \in \mathcal{T}$ (4)} \\ & \text{$\text{the} \,d$} } ^ { \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$
$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$
$\text$\text{$\text{$\text{$\text{$\text{$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\end{smallmatrix} } } } } } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a}\ } ^ { \text{In a} } } ^ { \text{In a} } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} } } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} } } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a}\ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} { \text{In a} \ } ^ ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \} ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a}" \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In ao} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{In a} \ } ^ { \text{Ina} \ } ^ { \text{In a} \ } ^ { \text{In a} \ }$$

∀

t

∈ T

(4)

$$\underline { x } _ { t } ^ { r } - v _ { t } ^ { r } \leq x _ { t } ^ { r } \leq \overline { x } _ { t } ^ { r } + v _ { t } ^ { r } ; \quad \forall t \in \mathcal { T } & & ( 5 ) \quad \underset { v _ { t } ^ { r }, \, v _ { t } ^ { r f } } { \sum } \text{and}$$

x

wh

t

-

v

wh

t

≤

x

wh

t

≤

x

wh

t

+

v

wh

t

;

∀

t

∈ T

(6)

$$\underline { x } _ { t } ^ { r f } - v _ { t } ^ { r f } \leq x _ { t } ^ { r f } \leq \overline { x } _ { t } ^ { r f } + v _ { t } ^ { r f } ; \quad \forall t \in \mathcal { T } \quad \quad \ \ ( 7 ) \quad \text{defined}$$

$$\underline { l } - v _ { t } ^ { l } \leq l _ { t } \leq \overline { l } + v _ { t } ^ { l } ; \quad \forall t \in \mathcal { T } ^ { o c c } & & ( 8 ) \quad \text{empera}$$

$$l _ { t } = \phi _ { t } u _ { t } ^ { b l } + \eta ^ { l e } u _ { t } ^ { a l } ; \quad \forall t \in \mathcal { T } ^ { o c c } \quad \quad ( 9 ) \quad \text{each} \, \text{ir}$$

$$u _ { t } ^ { b l } & \geq \underline { u } ^ { b l } ; \quad \forall t \in \mathcal { T } ^ { o c c } & ( 1 0 ) \quad \text{tempera} _ { \alpha ^ { r } \, - \, \tau }$$

$$\stackrel { \mu } { u } _ { t } ^ { b } & = \stackrel { \nu } { \sum _ { i } } ^ { \nu } u _ { t } ^ { i } + \sum _ { i \in \mathcal { I } } u _ { i t } ^ { u l } + P _ { t } ^ { s b } ; \quad \forall t \in \mathcal { T } & \stackrel { \nu \mu \nu } { ( 1 1 ) } & \stackrel { x _ { t } ^ { r } \, = \, x } { \text{tempera} }$$

1 FH stands for Floor Heater, while HVAC means Heating, Ventilation, and Air Conditioning system.

$$\text{public} _ { \text{asher}, \quad } 0 \leq u _ { t } ^ { \imath } \leq \overline { u } ^ { \imath } ; \forall \imath = \{ h e a t, \text{cool}, h p, w h, r f, a l \}, t \in \mathcal { T } \ \ ( 1 2 )$$

$$\delta _ { i t } ^ { u l } \leq T R _ { i t } ; \quad \forall i \in \mathcal { I }, t \in \mathcal { T }$$

$$\sin _ { \substack { \mu u - \mu \mu u, \\ \sum _ { t \in \mathcal { T } \\ \dots \sum _ { t } } } \delta _ { i t } ^ { u l } = C T _ { i } ; \quad \forall i \in \mathcal { I } \quad \ \ & \quad ( 1 4 )$$

$$( 2 \mathfrak { a } ) \quad \stackrel { \Phi } { t + C T _ { i } - 1 } { \sum _ { k = t } ^ { \Phi } \delta _ { i k } ^ { u l } \geq C T _ { i } ( \delta _ { i t } ^ { u l } - \delta _ { i, t - 1 } ^ { u l } ) ; \forall t \leq N _ { T } - C T _ { i } + 1 \, ( 1 5 ) } { \Phi } \quad \stackrel { \Phi } { \sim } \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \ \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \PhPh \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \PhI \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Php \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \PhP \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phs \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phis \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phj \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Ph } \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phip \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phim \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \Phi \quad \Phi \quad \Phi \quad \Phi \quad \Phi \Phi \quad \Phi \Phi \Phi \quad \Phi \Phi \quad \Phi \Phi \quad \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Psi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi\Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phs \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \ Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi$$

$$( 2 \mathfrak { m } _ { \mathcal { A } _ { \cdot } } ) \quad \delta _ { 1 i t } ^ { a } \geq \delta _ { i t } ^ { u l } - \delta _ { i, t - 1 } ^ { u l } ; \quad \forall i \in \mathcal { I }, t \in \mathcal { T }$$

$$( 2 c ) \quad \delta ^ { a } _ { c i t } \geq \delta ^ { a } _ { c - 1, i, t - 1 } ; \quad \forall c > 1, i \in \mathcal { I }, t \in \mathcal { T } \quad \quad ( 1 7 )$$

$$( 2 d ) \int _ { \substack { c \in \mathcal { C } \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots 0 \end{array}$$

$$( 1 \mathrm d ) = \quad \delta ^ { a } _ { c i t } \in [ 0, 1 ] ; \quad \forall c \in \mathcal { C }, i \in \mathcal { I }, t \in \mathcal { T } \quad \quad ( 1 9 )$$

$$\text{ilarly}, \quad \delta _ { i t } ^ { u l } \in \{ 0, 1 \} ; \quad \forall i \in \mathcal { I }, t \in \mathcal { T }$$

$$\underset { \text{rices} } { \text{node} } \quad u _ { t } ^ { b l } \in [ 0, 1 ] ; \quad \forall t \in \mathcal { T } \quad \text{(2)}$$

$$& \text{rices} ^ { - \mu _ { t } - \mu _ { t }, - \mu _ { t } } \cdots - \mu _ { t } \\ & \underset { v _ { t } } { \underbrace { \left ( v _ { t } ^ { r }, v _ { t } ^ { r f }, v _ { t } ^ { w h }, v _ { t } ^ { l } \right ) } ^ { T } } _ { v _ { t } } \geq 0 ; \quad \forall t \in \mathcal { T }, \quad \underset { v _ { t } ^ { b } } { \forall t } \in \mathcal { T }, \quad \underset { \lambda ^ { a } \cdot \lambda ^ { u l } } { \lambda ^ { b } \cdot \lambda ^ { u l } } \text{ is the set of variables} }$$

The goal of this optimization problem, given in (3), is the minimization of the electricity costs and penalty costs on violations of user-defined comfort constraints. We assume that a set of prices is known a priori by the household's energy management system. The state-space model of the household heat dynamics is given by constraints (4), whose matrices are described in [18]. This model includes a threestate model for the water-based FH, which is represented by the inside air temperature, the floor temperature, and the water temperature in the floor heating pipes, as similarly done in [8]. In addition, the state-space model (4) takes into account two additional states, i.e., the thermal dynamics of a residential refrigerator and a water heater located outside the household. The interested reader is referred to [5], [6] for complex models on these two appliances. Moreover, the HVAC system is also represented in the matrix equations given in (4) by a singlestate model (only the indoor air temperature). The external disturbances are the ambient temperature, room occupancy, and hot water demand.

where Ξ = ( l t , x t , u t , u b t , v t , δ a cit , δ ul it ) is the set of variables.

Comfort constraints are set in (5)-(10), expressions (11)(18) represent technical constraints, whereas (19)-(22) declare the character of variables δ a cit , δ ul it , u bl t , and the slack variables v r t , v rf t , v wh t , v l t .

Regarding the comfort constraints, (5)-(7) set the userdefined minimum and maximum bounds on the indoor air temperature, water temperature in the water heater, and air temperature in the refrigerator chamber, in that order, for each time period t . The comfort bounds for the indoor air temperature can be chosen as: x r t = x r,set t + αw t and x r t = x r,set t -αw t , where x r,set t is the set-point indoor air temperature, α is the maximum temperature difference with respect to the set-point that the user is able to withstand, and w t is a vector of continuous parameters varying between 0 and 1. Thus, comfort bounds were set up in two ways:

- · Price-independent comfort bounds (PI-CB): w t is a vector of ones.
- · Price-dependent comfort bounds (PD-CB): w t is a vector of normalized prices over a given day.

Constraints (8)-(10) model the household light levels when occupants are in the household, as described in [13]. Constraints (8) enforce recommended minimum and maximum light levels. The household light level given in (9) can be computed as the natural light coming through the windows, which depends on the blind positions (i.e., φ u t bl t ), and the indoor artificial lighting (i.e., η le u al t ). Constraints (10) may enforce a lower bound for the position of the blinds.

Regarding the technical constraints, expression (11) sets the building electricity consumption equal to the contributions from the heat pump of the FH, HVAC system, indoor artificial lighting, water heater, refrigerator, uninterruptible loads, and stand-by power. Constraints (12) impose the bounds on the power for the heating pump of the FH, HVAC system, artificial lighting, water heater, and refrigerator. The set of constraints (13)-(18) defines the operation of the uninterruptible loads. We assume that each uninterruptible load is characterized by its cycle time CT i and the cycle power at each phase c , i.e., P ul ci . In addition, the user may pre-define the time interval T i at which the load could be scheduled on so that the binary-valued parameter TR it = 1 if the uninterruptible load i could be on in period t , and TR it = 0 otherwise. Expressions (13) enforce the uninterruptible loads to be scheduled off outside the predefined time interval, whereas constraints (14) enforce binary variables δ ul it to be 1 during the cycle time within the predefined time interval. Constraints (15) impose that the scheduling during the cycle time must be consecutive. Constraints (16)-(17) set the relationship between the scheduling variable δ ul it and the activation variable δ a cit . Finally, expressions (18) define the power of the uninterruptible load i in period t as the power in the corresponding cycle phase.

## IV. CASE STUDY

The single-zone household presented in [8] is used to analyze what comfort settings and structural attributes will help make it price-responsive. The floor and window areas are respectively 30 and 1 m . For the water-based FH, 2 the data are based on [8]. The mass of water in the FH is 400 kg. The compressor has a coefficient of performance (COP) of 3 and a nominal power of 1 kW. When the household is equipped with an HVAC system, the COP for the heating and cooling system is respectively 1.67 and 3.67 and its electrical capacity is 1 kW. The household is equipped with a 30-l water heater located outside and a residential refrigerator. Data can be found in Table I. Inlet water temperature and the hot water daily consumption profile is given in [18].

Artificial lighting capacity is 60 W with an indoor luminous efficacy equal to 90 lumen/W. We assume that the luminous efficacy of daylight is 105 lumen/W in order to compute the outdoor illuminance in lux and the comfort light levels are set to 100 and 10000 lux. Note that the outdoor illuminance

TABLE I

DATA FOR THE REFRIGERATOR [5], [15] AND WATER HEATER [6]

|                                    |   Refrigerator |   Water heater |
|------------------------------------|----------------|----------------|
| Nominal power capacity [kW]        |          0.35  |           1.26 |
| COP                                |          0.76  |           0.92 |
| Thermal capacity [Wh/ ◦ C]         |          6.65  |          34.85 |
| Heat transfer coefficient [W/ ◦ C] |          0.678 |           0.5  |

TABLE II DATA FOR THE COMFORT CONSTRAINTS

|                                                                          | noflex                                          | flex        | extraflex   |
|--------------------------------------------------------------------------|-------------------------------------------------|-------------|-------------|
| x r,set t [ ◦ C] α [ ◦ C]                                                | 20 0                                            | 20 2        | 20 5        |
| { x wh , x wh } [ ◦ C]                                                   | { 54, 56 }                                      | { 50, 60 }  | { 45, 65 }  |
| { x rf , x rf } [ ◦ C]                                                   | { 4.9, 5.1 }                                    | { 4, 5 }    | { 3, 6 }    |
| Washing machine Dishwasher (first) Dishwasher (second) Tumble dryer Oven | 06:00-14:00 16:00-00:00 15:00-00:00 10:00-15:00 | 06:00-14:00 |             |

and the comfort light levels are respectively multiplied by the window and floor area to convert them to lumen.

We consider 4 uninterruptible loads: oven, washing machine, tumble dryer, and dishwasher. Their schedules are described in [18], their desirable operating hours are given in Table II, and their power on each cycle phase can be found in [19] for year 2020.

The discretization step ∆ t is assumed to be 15 min to properly capture the building dynamics and the penalty terms ρ temp and ρ l are set to 1000. We run daily simulations with 15-min time steps for one year. For each simulation, we use a look-ahead window of one day to account for the future impacts of thermal dynamics when control actions are taken in a given period of time. Ambient temperature, solar radiation, and electricity prices are given in [18] for the sake of reproducibility. This reference also includes the occupancy schedules. Note that the standby power is neglected.

Three different cases regarding the degree of household flexibility are analyzed: noflex , flex , and extraflex cases. Table II also provides the data on setting comfort bounds for the indoor air, refrigerator, and water heater temperature.

The simulations have been performed on a Linux-based server with one CPU clocking at 2.6 GHz and 2 GB of RAM using CPLEX 12.6.3 [20] under Pyomo 5.2 [21]. Optimality gap is set to 0%.

## A. Effect of Comfort Settings

Table III provides results for a single-zone household under two types of space heating (FH and HVAC), three cases for the degree of flexibility ( noflex , flex , and extraflex ), and the strategy followed to set the comfort bounds (PI-CB and PD-CB). This table shows the annual electricity cost, the total violations (i.e., number of Celsius degrees out of the corresponding comfort bounds over the 35040 time periods of the year), the percentage of 15-minute time intervals in which the indoor air temperature lies within three different ranges,

TABLE III RESULTS FOR A SINGLE-ZONE HOUSEHOLD

| Type   | Comfort bounds   | Case                  | Annual cost [ e ]   | Violations [ ◦ C · h]   | Freq. at 20 ◦ C [%] a   | Freq. [18,20)- (20,22] ◦ C [%] a   | Freq. [15,18)- (22,25] ◦ C [%] a   | Building cons. [kWh]   | u b [%] b      | u hp/heat [%] b   |
|--------|------------------|-----------------------|---------------------|-------------------------|-------------------------|------------------------------------|------------------------------------|------------------------|----------------|-------------------|
| FH     | PI-CB            | noflex flex extraflex | 103.5 78.9 68.4     | 862.5 0.0 0.0           | 37.5 1.5 0.0            | 62.5 98.5 0.3                      | 0.0 0.0 99.7                       | 1944.0 1569.1 1372.7   | 54.5 62.7 63.5 | 51.3 86.6 85.8    |
| FH     | PD-CB            | noflex flex extraflex | 103.5 93.2 87.0     | 862.5 55.9 35.2         | 37.5 22.2 21.7          | 62.5 77.8 78.3                     | 0.0 0.0 0.0                        | 1944.0 1843.9 1751.3   | 54.5 61.9 64.1 | 51.3 70.0 73.1    |
| HVAC   | PI-CB            | noflex flex extraflex | 107.0 89.9 76.1     | 0.0 0.0 0.0             | 100.0 1.1 0.5           | 0.0 98.9 19.6                      | 0.0 0.0 79.9                       | 2044.1 1793.0 1538.8   | 57.2 63.5 64.7 | 60.2 77.2 80.1    |
| FH     | PD-CB            | noflex flex extraflex | 107.0 94.5 84.2     | 0.0 0.0 0.0             | 100.0 11.1 8.0          | 0.0 88.9 61.6                      | 0.0 0.0 30.4                       | 2044.1 1867.7 1691.2   | 57.2 62.3 64.4 | 60.2 70.6 73.9    |

a Columns 6-8 represent the percentage of 15-minute time intervals lying within the given temperature intervals throughout the year.

b Columns 10-11 represent the share of power lying within low-price periods.

the building energy consumption, and the share of building and heating power lying in low-price periods (i.e., those periods where the annual normalized prices are lower than 0.5). Note that we show the frequency of time intervals in which the indoor air temperature lies within: (i) 20.0 ◦ C (column 6), (ii) [18, 20)-(20, 22] ◦ C (column 7), and (iii) [15, 18)-(22, 25] ◦ C (column 8). For all cases, computing times for running annual simulations were in the range of 960-3600 s.

The HVAC system leads to more expensive solutions than the FH, but provides in general less discomfort assuming that the reference set-point of 20 C is the temperature at which ◦ the occupants experience the highest comfort. Note that for the noflex case, the HVAC is able to keep the air temperature within comfort bounds (i.e., at the reference set-point), unlike a household with FH (this is also due to the absence of air conditioning when considering FH). We can notice that the annual cost decreases when increasing the degree of flexibility at the expense of a higher degree of discomfort regardless of the comfort settings. Also, the cost reduction is higher under the PI-CB strategy when increasing the degree of flexibility (33.9% vs. 15.9% for the extraflex case with respect to the noflex case with the FH). Similar conclusions can be drawn with the HVAC system. However the PD-CB strategy leads to higher costs than under the PI-CB when the occupants are flexible regardless of the type of space heater.

The PD-CB strategy results in time-varying comfort bounds throughout the day, which is the reason why it leads to higher costs than the PI-CB. However, we can clearly observe that the PD-CB leads to less discomfort (assuming that the reference set-point of 20 ◦ C is the temperature at which the occupants experience the highest comfort). Note that the percentage of 15-minute time intervals lying within [15, 18)-(22, 25] ◦ C is significantly reduced under the PD-CB compared to that percentage under the PI-CB.

Columns 10-11 in Table III provide the share of power (total in the building and that of the space heating system) lying in low-price periods. This would be an indicator whether a household is more or less price-responsive. Percentage of total power (column 10) in low-price periods increases as the

TABLE IV UA r,a

EFFECT OF ON ANNUAL COST AND DISCOMFORT METRICS

| Type   | Cost and metrics                                                           | Factor       | Factor       | Factor               | Factor         |
|--------|----------------------------------------------------------------------------|--------------|--------------|----------------------|----------------|
| Type   | Cost and metrics                                                           | 0.5          | 1            | 2                    | 4              |
| FH     | Cost [ e ] Violations [ ◦ C · h] Freq. at 20 ◦ C [%] [18,20)-(20,22] ◦ C ◦ | 76.6 0 0 100 | 86 0 0 100 0 | 105.8 1.1 2.9 97 0.1 | 149.8 93.4 4.2 |
| FH     | Freq. [%]                                                                  |              |              |                      | 95.4           |
| FH     | Freq. [15,18)-(22,25] C [%]                                                | 0            |              |                      | 0.2            |
| HVAC   | Cost [ e ]                                                                 | 87           | 106.5        | 146                  | 225            |
| HVAC   | Violations [ ◦ C · h]                                                      | 0            | 0            | 0                    | 0              |
| HVAC   | Freq. at 20 ◦ C [%]                                                        | 0.9          | 0.6          | 0.3                  | 0.2            |
| HVAC   | Freq. [18,20)-(20,22] ◦ C [%]                                              | 99.1         | 99.4         | 99.7                 | 99.8           |
| HVAC   | Freq. [15,18)-(22,25] ◦ C [%]                                              | 0            | 0            | 0                    | 0              |

occupants are more flexible. This trend becomes even more noticeable when it comes to the heating power since space heaters have a higher degree of control than other loads such as refrigerators, dishwashers, etc. Both space heating systems allow for shifting energy to low-price periods when increasing the flexibility (around 7.2-7.5% of difference between the share of building power for the extraflex and noflex cases with HVAC versus 9.0-9.6% with FH). Although there are some differences, both systems can be price-responsive when using a look-ahead window.

## B. Effect of Structural Parameters

Structural parameters such as the heat transfer coefficients 2 could be also responsible for the household priceresponsiveness. Table IV shows the effect of the heat transfer coefficient between the room air and the ambient, UA r,a , for the flex case on the annual cost, total number of violations of comfort bounds, and the percentage of 15-minute time intervals in which the room temperature lies within three different ranges, whereas Table V shows the share of building power 3 lying within low-price periods for different flexibility cases. In order to isolate the effect of the heat transfer coefficient, we

2 The heat transfer coefficient UA is the product of the heat conductivity U and the surface area A where the heat transfer takes place.

3 Building power refers to the total consumption of the smart household.

TABLE V EFFECT OF UA r,a ON THE SHARE OF BUILDING POWER CONSUMPTION LYING WITHIN LOW-PRICE PERIODS [%]

|        | FH     | FH   | FH        | HVAC   | HVAC   | HVAC      |
|--------|--------|------|-----------|--------|--------|-----------|
| Factor | noflex | flex | extraflex | noflex | flex   | extraflex |
| 0.5    | 54.9   | 64.2 | 65.5      | 56.1   | 63.0   | 65.2      |
| 1.0    | 54.0   | 67.2 | 68.1      | 55.9   | 62.3   | 64.6      |
| 2.0    | 53.1   | 71.0 | 71.7      | 55.7   | 60.2   | 62.2      |
| 4.0    | 51.7   | 72.5 | 73.4      | 55.6   | 57.8   | 58.6      |

assume that the household does not have any window and that the lighting constraints are ignored. The value of UA r,a has been multiplied by factors of 0.5, 1, 2, and 4, respectively. The greater the value of UA r,a , the less insulated the household is.

It can be observed than when the household is less insulated, the annual costs increase because there are more thermal exchanges through the walls at the expense of slightly increasing the occupants' discomfort (see Table IV). Since this table shows the results from the flex case, most of the time periods are above 18 ◦ C and below 22 C. ◦

The household with an FH is more price-responsive when increasing the value of UA r,a (see Table V). We can observe that the differences between the extraflex and noflex cases vary between 10.6-21.7% for the factors 0.5-4, respectively. Conversely, the household with an HVAC system is more price-responsive when decreasing the value of UA r,a . This is due to the fact that the HVAC system is an appliance with fast dynamics whereas the water-based FH is characterized by slow dynamics.

## V. CONCLUSIONS

We propose a compact and integrated formulation for a smart building with smart appliances via economic model predictive control. We model five states to capture thermal dynamics of space heaters, refrigerator and water heater; and we accurately model variable power cycles for uninterruptible loads.

Substantial cost savings can be achieved when increasing the comfort bounds under HVAC systems (21-29%) and waterbased floor heating systems (16-34%). Price-responsiveness of smart households is quite similar between space heaters of slow (water-based FH) or fast impact (HVAC system) if a lookahead window is used. The amount of building consumption lying within low-price periods increases from 54.5-57.2% to 64.1-64.7% when using wider comfort bounds for the smart appliances. In general, the more price-responsive the household is, the higher discomfort the occupants experience. However, price-dependent comfort bounds could help reduce the occupants' discomfort, thus reducing the number of time intervals at temperatures far from the reference set-point. Insulated households with HVAC systems are more priceresponsive and they may lead to cost savings up to 50% approximately, whereas price-responsiveness in households with water-based FHs is higher for less insulated ones at the expense of increasing the annual costs.

Further research will be devoted to exploring the use of battery energy storage systems and local renewable generation to increase the flexibility degree of smart buildings.

## REFERENCES

- [1] D. S. Kirschen, 'Demand-side view of electricity markets,' IEEE Transactions on Power Systems , vol. 18, no. 2, pp. 520-527, May 2003.
- [2] Y. Dvorkin, R. Fern´ andez-Blanco, D. Kirschen, H. Pandzic, J.-P. Watson, and C. Silva-Monroy, 'Ensuring profitability of energy storage,' IEEE Transactions on Power Systems , vol. 32, no. 1, pp. 611-623, Jan. 2017.
- [3] 'Benchmarking smart metering deployment in the EU-27 with a focus on electricity COM(2014) 356 final,' tech. rep., European Commission, Brussels, 2014.
- [4] Y. Parag and B. K. Sovacool, 'Electricity market design for the prosumer era,' Nature Energy , vol. 1, no. 4, p. 16032, 2016.
- [5] G. T. Costanzo, F. Sossan, M. Marinelli, P. Bacher, and H. Madsen, 'Grey-box modeling for system identification of household refrigerators: A step toward smart appliances,' IYCE 2013 - 4th International Youth Conference on Energy , 2013.
- [6] F. Sossan, A. M. Kosek, S. Martinenas, M. Marinelli, and H. Bindner, 'Scheduling of domestic water heater power demand for maximizing PV self-consumption using model predictive control,' 2013 4th IEEE/PES Innovative Smart Grid Technologies Europe, ISGT Europe 2013 , 2013.
- [7] R. Halvgaard et al. , 'Model predictive control for a smart solar tank based on weather and consumption forecasts,' Energy Procedia , vol. 30, pp. 270-278, 2012.
- [8] R. Halvgaard, N. K. Poulsen, H. Madsen, and J. B. Jorgensen, 'Economic model predictive control for building climate control in a smart grid,' 2012 IEEE PES Innovative Smart Grid Technologies (ISGT) , pp. 1-6, 2012.
- [9] S. Li, W. Zhang, J. Lian, and K. Kalsi, 'Market-based coordination of thermostatically controlled loads -Part I: A mechanism design formulation,' IEEE Transactions on Power Systems , vol. 31, no. 2, pp. 1170-1178, 2016.
- [10] A. G. Thomas and L. Tesfatsion, 'Braided cobwebs: Cautionary tales for dynamic pricing in retail electric power markets,' IEEE Transactions on Power Systems , vol. 33, no. 6, pp. 6870-6882, Nov. 2018.
- [11] J. Saez-Gallego and J. M. Morales, 'Short-term forecasting of priceresponsive loads using inverse optimization,' IEEE Transactions on Smart Grid , vol. 9, no. 5, pp. 4805-4814, Sep. 2018.
- [12] L. Zhao and W. Zhang, 'A unified stochastic hybrid system approach to aggregate modeling of responsive loads,' IEEE Transactions on Automatic Control , vol. 63, pp. 4250-4263, Dec. 2018.
- [13] J. E. Contreras-Oca˜ na, M. R. Sarker, and M. A. Ortega-Vazquez, 'Decentralized coordination of a building manager and an electric vehicle aggregator,' IEEE Transactions on Smart Grid , vol. 9, no. 4, pp. 26252637, Jul. 2018.
- [14] M. F. Anjos, A. Lodi, and M. Tanneau, 'A decentralized framework for the optimal coordination of distributed energy resources,' IEEE Transactions on Power Systems , vol. 34, no. 1, pp. 349-359, Jan. 2019.
- [15] J. M. Gonzalez et al. , 'Smart residential load simulator for energy management in smart grids,' IEEE Transactions on Industrial Electronics , vol. 66, no. 2, pp. 1443-1452, Feb. 2019.
- [16] R. G. Junker et al. , 'Characterizing the energy flexibility of buildings and districts,' Applied Energy , vol. 225, pp. 175-182, Sep. 2018.
- [17] R. Halvgaard, J. B. Jørgensen, N. Kjølstad, and H. Madsen, Model Predictive Control for Smart Energy Systems . Thesis, Kgs. Lyngby: Technical University of Denmark (DTU), 2014.
- [18] R. Fern´ andez-Blanco, J. M. Morales, and S. Pineda, 'How can smart buildings be price-responsive? - Auxiliary data.' [Online]. Available: https://goo.gl/WqrhPp, 2019.
- [19] R. Stamminger, G. Broil, C. Pakula, H. Jungbecker, M. Braun, I. R¨ udenauer, and C. Wendker, 'Synergy potential of smart appliances. Report of the Smart-A project,' tech. rep., 2008.
- [20] 'IBM ILOG CPLEX Optimisation Studio.' [Online]. Available: https://www.ibm.com/analytics/cplex-optimizer, 2019.
- [21] W. E. Hart, J. P. Watson, and D. L. Woodruff, 'Pyomo: Modeling and solving mathematical programs in Python,' Mathematical Programming Computation , vol. 3, no. 3, pp. 219-260, 2011.
Ricardo Fernández-Blanco Juan Miguel Morales Salvador Pineda
How Can Smart Buildings Be Price-Responsive?
The prospective participation of smart buildings in the electricity system is
strongly related to the increasing active role of demand-side resources in the
electrical grid. In addition, the growing penetration of smart meters and
recent advances on home automation technologies will spur the development of
new mathematical tools to help optimize the local resources of these buildings.
Within this context, this paper first provides a comprehensive model to
determine the electrical consumption of a single-zone household based on
economic model predictive control. The goal of this problem is to minimize the
electricity consumption cost while accounting for the heating dynamics of the
building, smart home appliances, and comfort constraints. This paper then
identifies and analyzes the key parameters responsible for the price-responsive
behaviour of smart households.