Application of Predictive Control in Scheduling of Domestic Appliances

: In this work, an algorithm for the scheduling of household appliances to reduce the energy cost and the peak-power consumption is proposed. The system architecture of a home energy management system (HEMS) is presented to operate the appliances. The dynamics of thermal and non-thermal appliances is represented into state-space model to formulate the scheduling task into a mixed-integer-linear-programming (MILP) optimization problem. Model predictive control (MPC) strategy is used to operate the appliances in real-time. The HEMS schedules the appliances in dynamic manner without any a priori knowledge of the load-consumption pattern. At the same time, the HEMS responds to the real-time electricity market and the external environmental conditions (solar radiation, ambient temperature, etc.). Simulation results exhibit the beneﬁts of the proposed HEMS by showing the reduction of up to 70% in electricity cost and up to 57% in peak power consumption.


Introduction
Building energy consumption contributes approximately 20 to 40% to the total energy consumption in developed countries and has exceeded other major sectors like industrial and transportation [1]. Smart grid is a major step forward to an energy efficient future of the humankind which allows integration of advanced sensing and communication technologies and various control methodologies in order to achieve optimum energy flow [2]. One of the important features of smart grid is the development and incorporation of demand-response strategies [3].
Demand-response [4,5] refers to the change in power consumption of an electric utility end-use customer to better match the load-demand with the power supply for reliable function of the power grid. The customer may adjust the load-demand by postponing some tasks that require larger amount of electric power or may decide to pay higher prices for the electricity in order to complete the tasks. The utility companies may provide incentives to the customer for demand-response actions. Obviously, it is difficult for the consumer to manually participate in demand-response programs. This explicates the need of an automated system to communicate with smart grid and make optimal decisions on behalf of the customer to achieve the goal of demand-response strategy.
This automated system usually called home energy management system (HEMS) is connected to smart-grid by means of bi-directional communications. Advanced sensors and smart-meters are employed to receive/send data and control signals between smart-home and the power utility [6]. Advanced control algorithms are incorporated into HEMS to schedule or operate the appliances in the desired way. The development of appropriate scheduling algorithms has been isolated as one of the crucial challenges for the next generation of real-time systems [7]. Extensive research has been dedicated to the problem of electric load scheduling with focus on different objectives, such as customer comfort [8], minimization of electricity cost [9], reduction of energy usage [10], shifting the electric load [11], etc. Moreover, various methods have been used to investigate this problem, which include, mixed-integer programming [12][13][14][15], stochastic programming [16,17], evolutionary algorithms [18,19], heuristic-based algorithms [20][21][22][23][24], game-theory [25], and learning-based algorithms [26,27].
Setlhaolo et al. [28] presented a study on optimal scheduling of typical household appliances. The scheduling model has been formulated as a non-linear integer program by considering the electricity cost, participation incentive and scheduling-inconvenience in the objective function. The approach showed a 25% reduction in the electricity cost for the particular case. Shirazi et al. [13] proposed a home energy management system for a smart home equipped with distributed energy resources (DERs) and thermal storage facility. The proposed scheduling technique considered the energy cost and the peak-load demand in the multi-objective function and analysed the results under different weather season scenarios. Sou et al. [14] investigated deterministic problem of day ahead scheduling of appliances by modeling the decision problem in more realistic way and show benefit of the proposed approach by looking into two case studies based on different tariffs in Sweden and NYC. However, the cost calculation in the approach is not based on real-time prices but on tariff which is known 24 h in advance. A few studies on the appliance-scheduling with different pricing schemes, such as time-of-use (TOU) pricing [28], inclining block rate (IBR) pricing [10], and their combinations have also been reported in the literature. In Reference [29], a new pricing scheme active consumption level pricing scheme (ACLP), based upon the consumption level (CL) of consumers has been proposed. The scheme encourages consumers to keep their energy consumption within a price-invariant-band (PIB). The proposed scheme is able to reduce the electricity cost by up to 53% and peak load by up to 35%. In addition to that, some researchers have focused on the load-scheduling based upon the priorities of the appliance operation [30,31]. The authors in Reference [32] design a price-based HEM framework where priorities of operating different appliances are interpreted as the value of loss load (VOLL). The reliability cost which is a function of VOLL is incorporated in the objective function. The results demonstrated 7.5% and 12% reduction in electricity cost with TOU and IBR pricing schemes, respectively.
The present work focuses on real-time appliance scheduling, i.e., to meet the immediate load demand without a priory knowledge of future load profile. In other words, instead of modeling the future load demand and optimizing, the optimization problem is solved at each time step using MPC (receding horizon control). MPC is an advanced method of control that emerges from application in process industry in late 70s and early 80s [33]. MPC represents a class of advanced control methods in which the model of the process is considered explicitly to predict the future evolution of the process to optimize the control input while respecting certain constraints. So, in this case, the optimization problem is solved at each time step with updated values of the real-time electricity prices and other relevant information like external weather conditions. In literature, many studies have investigated the use of MPC strategy for demand response [34][35][36] to achieve peak load reduction. In these studies, the focus is on the control of HVAC system in the building [37][38][39]. Another popular application of MPC is in residential microgrid management [40][41][42][43]. In these studies, the residential house is considered to be equipped with various distributed energy resource and MPC is used to optimally control them. However, the scope of the present work is different in the manner that it investigates use of predictive control with receding horizon strategy to schedule thermal and non-thermal appliances both. For the thermal appliances, at each time step, the future external disturbance and constraints are updated; whereas, for non-thermal appliances the future deadline of the appliance is updated at each time step. The user has the flexibility to change the deadline as desired while appliance is still in operation, which consequently change the optimization problem.
In literature, there are some studies which investigate the problem of scheduling of thermal and deferrable appliances using MPC. The authors in Reference [44] investigated the problem of load-scheduling of thermal and non-thermal appliances using MPC under different price schemes and achieves the total electricity cost saving of up to 20%. But the previous study do not focus on the peak-load reduction by employing constraints on the total maximum power consumption. In Reference [45], the author addressed the similar problem of deferrable appliance-scheduling considering distributed generation (DG) using a multi-time stochastic MPC to consider the uncertainty emerged due to the intermittent DG power. In the best case scenario, the total electricity cost reduction of 53% is achieved in this study. The limitation of the previous study is that some deferrable appliances are operated during a fixed optimal time interval only which is obtained using genetic algorithm. This takes away the flexibility for the customer to intervene for using the appliances at desired time intervals.

Research Contribution
In the present study, an architecture of HEMS is presented for automated scheduling of appliances with an objective of reducing the peak-power consumption and the total electricity cost. It is assumed that the bi-directional communication between the house and the power grid is present, which enables the HEMS to receive/send data and control signal from the power utility. HEMS is employed in real-time electricity pricing environment.
The HEMS is designed to operate in two modes of operation (MOO) based upon the load category (thermal/non-thermal) (see Figure 1). Mode 1 remains in operation if any of the non-thermal appliances is active. On the other hand, in Mode 2, only the thermal appliances are activated. At each time step, an optimization problem is formulated into a linear-program (LP) with appropriate objective function and constraints. The HEMS functions in a receding horizon framework, i.e., at each time step, updated system states and predicted information is provided to the controller and the optimization process is repeated. In both the modes, soft constraints on the maximum power consumption of the appliances is also imposed. This way, the total peak power consumption is maintained under the specified peak-consumption capacity. The constraint values are decided by the utility company or the aggregator based upon the demand of electricity and operational cost of the grid. These values are communicated to the smart-meter in homes for a fixed duration in advance and updated at each time step. The detailed functionalities of HEMS are described in later sections. The implemented scheduling algorithm does not require a priori knowledge of the load-demand. The non-thermal appliance can be activated at any time by the customer. The customer is required to set a deadline at the time of activation by which the task assigned to the appliance should be finished. The HEMS outputs a warning if the deadline set by the customer is shorter than the running time of that appliance.

Load Categorization
Before going into a detailed description and working methodology of the HEMS, the classification of electrical load is detailed in this section. The load is classified based upon the dynamics of the appliances. The classification is important for the HEMS, based upon this classification it makes the choice that in which MOO it should run. So, the overall electrical load is divided into two categories.

Thermal Load
This category includes the electrical load consumed by the appliances which are thermal in nature, i.e., their system dynamics include temperature as a system-state, e.g., heat pump, room heater, water heater, refrigerator, etc.

Non-Thermal Load
In this category, electrical load consumed by the appliances which are non-thermal in nature is included. These non-thermal appliances are assumed to have following characteristics.

•
Preemptive-The appliance can be paused or interrupted whenever needed.

•
Fixed duration-The appliance run for a fixed duration (Ω) to finish the task assigned to it. • Deadline-The appliance has a deadline associated by which it needs to finish the task assigned to the appliance. The deadline is chosen by the user and cannot be smaller than the running duration of the appliance which is an operational constraint [44]. Examples include washing machine, clothes dryer, dishwasher, etc.

System Layout
This section provides a detailed description of the HEMS. The system architecture of the HEMS is presented in the Figure 1. The main components of the HEMS are described below.

•
The first component of the HEMS is the forecaster. It is deputized to make predictions about the external weather and the heat-gains due to solar irradiance in the house. Data-driven prediction techniques [46,47] are commonly used to make these predictions based on available historical data about the variable. The second component of the HEMS uses these predictions to optimally schedule the thermal appliances.

•
The second and the key component of the HEMS is the Predictive MILP controller. This component operates in two MOO, which is determined automatically within the component based upon the electrical load categories (see Section 3.1). In each mode, a different mathematical optimization problem is solved by exploiting the predictions made by the forecaster about the weather conditions. The electricity prices and the consumption-capacity profile is provided to HEMS by the smart-grid operator using smart-meter communication.

State-Space Model
In order to apply the proposed scheduling algorithm, the dynamics of all the appliances are modeled mathematically into a discrete linear time-invariant state-space representation. The discrete state-space model of a linear system with control inputs (m), outputs (p), state variables (n), and external disturbances (r) is written as in Equation (1) [48]: (1) x ∈ R n is the state vector, u ∈ R m is the control vector, y ∈ R p is called the output vector, and d ∈ R r is called the disturbance vector. Similarly, A ∈ R n×n is the state matrix, B ∈ R n×m is the control matrix, C ∈ R p×n is the output matrix, and E ∈ R r×n is the disturbance matrix. k denotes the time instant.

Thermal Load Appliances Modeling
For the purpose of heating the indoor environment, the house is considered to be equipped with a water source heat pump through floor heating pipes [49]. For hot water requirements, building is connected to a solar water heater system. The dynamics of both the systems is influenced by the solar irradiance and the ambient temperature.
Heat pump: A reduced order model of a residential building equipped with a heat-pump is developed in Reference [49,50] as follows using the heat balance equations: The system of Equation (2) is discretized with appropriate sample time using zero-order-hold model and written into state-space representation as in Equation (1). The first sub-equation in Equation (2) represent the thermal dynamics of the indoor environment of the building which depends upon the floor temperature T f , ambient temperature T a and solar gains φ s . The second sub-equation captures the heat transfer between floor and the underfloor heating system. The last sub-equation relates the temperature of the floor T f to the work done W by the compressor.
Subsequently, for heat pump state-space model, we get Solar water tank: The heat dynamics of the solar water tank is described using simple first-order differential equations as follows [51] Equation (3) captures the thermal dynamics of the tank as a function of inlet water temperature T i , power consumption of the heating element P h , and solar gains φ s . After discretization and rewriting Equation (3) in state-space representation, we get x = y = T t ,

Non-Thermal Load Appliances Modeling
In the present work, a washing machine and a dishwasher are considered in the non-thermal load appliances category. Any non-thermal load appliance can be modeled in the same way since they have the same characteristics as described in Section 3.1.
The dynamics of a non-thermal load appliance is modeled in the following way: where ζ is the variable defined to keep track of the duration up to which the appliance has run, and ψ is a binary integer variable which can have the values from the set {0,1}. It represents that if the appliance is running at time instant k then the value of ψ(k) will be 1 and if the appliance is halted, ψ(k) will be 0.
In both the cases, the value of ψ is added to ζ. Therefore, for non-thermal appliances, Equation (4) can be represented in form of Equation (1), which subsequently gives x = y = ζ and u = ψ. It should be noted that, whenever ψ take the value of 1, a certain amount of power is consumed by the appliance.

Problem Formulation and Scheduling Algorithm
In this section, the task of the scheduling of appliances is formulated as an MILP. The appropriate constraints associated with the appliances are also described. The objective of the program is to reduce the total electricity cost and peak power consumption at the same time. As mentioned previously, an MPC framework is applied to operate the appliances. At each time step an MILP is formulated with a prediction horizon of N time steps ahead from the current time step.

Energy Cost of Thermal Appliances
The electricity cost occurred by only the energy consumption of thermal appliances is described in Equation (5) below: where N ∈ {1, 2, . . . , N}, N is the prediction horizon. J P is the energy cost for next N time steps for running the thermal appliances. c k is the real-time electricity price at kth time step. u k,i represents the power consumption by the ith thermal appliance at the time step k.

Comfort Zone Constraints of Thermal Appliances
The user would prefer to keep the temperature of the building into a prescribed comfort zone. Similarly, temperature of water in the water tank also should be in a certain temperature range according to users' preferences. These constraints are imposed in Equation (6) as where y max,i and y min,i are upper and lower bound on the temperature of the ith thermal appliance.

Maximum Power Consumption Constraint of Thermal Appliances
Each thermal appliance is assumed to consume up to a certain amount of power at each time step and cannot exceed that value, i.e., where u min,i and u max,i represent the upper and lower bound on the power consumption of ith thermal appliance.

Energy Cost of Non-Thermal Appliances
Same as the thermal appliances, the electricity cost induced by the non-thermal appliances energy consumption can be considered as J Q is the total energy cost for operating the activated non-thermal appliances. P k,j is the power consumed by jth non-thermal appliance at kth time step. ψ k,j is the binary variable associated to the jth non-thermal appliance at kth time step, which indicates that the appliance is ON or OFF if the value of ψ is 1 or 0, respectively.

Deadline Constraint of Non-Thermal Appliances
Each non-thermal appliance run for a fixed duration as described in the Section 3.1. It has an associated deadline, by which it should finish the assigned task. This deadline is set by the user while activating the appliance and cannot be less than the running time of appliance, i.e., if user set a smaller deadline then HEMS will give a warning to reset the deadline. This can be assured by imposing the constraint in Equation (9) Ω j ≤ (e j − b j ), j ∈ Q, (9) where Ω j is total running time of jth appliance for which the appliance has to run in order to finish the task. e j and b j is the deadline (end time) and activation time (beginning time) of the jth appliance. It is assumed, when the appliance is turned on HEMS does not schedule it immediately rather it waits for current time slot to complete or next nearest time step to start. It should be noted that, in the present work, this constraint is not implemented in the optimization problem explicitly, rather the starting time and deadlines for non-thermal appliances are determined offline for entire simulation period as per this constraint. However, in the real-life implementation of home energy management system this constraint can be implemented directly in the optimization problem in real-time.

The Non-Thermal Appliance Start-Up Cost
In order to prevent the non-thermal appliances from being interrupted very often, there is a start-up cost associated with the each appliance, which can be represented as Equation (10) γ j refers to the start-up cost of jth non-thermal appliance.

Operation Time Constraint of Non-Thermal Appliances
To finish the task assigned, the non-thermal appliances need to run for a fixed duration Ω. This condition can be imposed as the following constraint: The constraint in Equation (11) ensures that appliance run for the its total running time between the activation time and the deadline of the appliance.

Total Capacity Constraint
This constraint put a restrain to the maximum total power consumption by all the appliances at each time step. It can be represented as following C k is the total maximum power available for consumption at the time step k. The capacity constraint value is decided by the utility company. In reality, the capacity constraint value can be time varying based upon the demand of electricity and operational cost of the grid. This constraint is a hard constraint, i.e., if this constraint is not satisfied the optimization problem will halt which is not desired.
To relax this constraint on total capacity, a soft constraint is introduced in next subsection.

Soft Constraints
Slack variables are introduced to allow the violation of some of the relevant constraints in Equation (13). For instance, if the desired temperature is not achieved or total power consumption exceeds the maximum available capacity, optimization routine automatically compromises between the optimal scheduling and minimization of constraint violation. y k,min,i ≤ y k,i + v k,i k ∈ N , i ∈ P, (13a) where v and w are the introduced slack variable to allow violation of constraints.

Scheduling Algorithm
As the HEMS is invoked, first, it checks if any non-thermal appliance is activated from the value of activation flag (1 or 0) associated with each non-thermal appliance. If the value of activation flag is 1, then that particular non-thermal appliance is running otherwise not. Based on the value of the activation flag, the HEMS decides to operate in the following two modes:

Mode 1
The HEMS operates in this mode if any one or more than one non-thermal appliances are activated. In this case, the optimization problem is formulated as an MILP as follows: subject to Here, α v (v) ≥ 0 and β w (w) ≥ 0 are convex penalty cost functions associated with the slack variables v and w. It should be noted that the value of β w (w) is time dependent and is specified by the power utility depending upon the power-grid operation and stability. This information is shared with HEMS dynamically on a real-time basis.

Mode 2
In Mode 2, since only the thermal appliances are running, the optimization problem is formulated just as a linear program (LP) rather than MILP min J P , (16) subject to The above constraint considers the total power consumption by only the thermal appliances, since no non-thermal appliance is running.
In both the modes, at each time step, a control decision policy {u * k } N k=1 is obtained for next N time steps by solving the formulated optimization. MILP and only the first control input {u * 1 } is applied to the system and rest of the control inputs {u * k } N k=2 are discarded. The optimization problem is formulated and solved again in the similar manner at the next time step as per the receding horizon strategy.
For the present study, it is assumed that perfect electricity prices forecast and weather forecast are available as this is not the focus of the paper.

Simulation Setup
As mentioned previously, we consider heat pump and solar water heater as thermal appliances for current simulation studies. Table 1 contains the values of the parameters used in modeling of thermal appliances as described in Reference [49,51]. To put a constraint on the indoor air temperature, (18 • C, 22 • C) is specified as the desired comfort range. These constraints can be chosen to be time varying also. The heat pump is considered to consume maximum of 1 kW power at each time step. Similarly, the water temperature in. the solar water heater tank is kept bounded between (50 • C, 70 • C) and the maximum power consumption of the water heater is considered to be 2 kW. Everyday water withdrawal demand profile is obtained from the DHWcalc toolbox [52], which generates realistic domestic hot water demand depending upon various characteristics like number of households, total mean daily draw off volume, etc.
Without any loss of generality, a washing machine and dishwasher are considered as the non-thermal household appliances. Table 2 details the parameter values for the washing machine and the dishwasher. Every day, non-thermal appliance requests are activated at random time steps with a random deadline which is greater than the total running time of the appliance. Table 1. Thermal appliance parameters. The entire simulation setup is implemented in MATLAB. For simulation purpose, ambient temperature T a and solar radiation φ s are obtained from ASHRAE IWEC (International Weather for Energy Calculations) weather data files for Dublin, Ireland [53]. The electricity prices are taken from wholesale Danish Energy Market data [54] and employed for real time implementation of the proposed algorithm. A maximum capacity constraint of 4 kW is imposed on the total power consumption of all the appliances. In reality, the capacity constraints can be time varying. The setup is simulated over a 7-day period with a prediction horizon of 24 h (Mode 1). To solve the MILP, Gurobi solver is used in integration with YALMIP toolbox [55] and MATLAB. To summarize, Table 3 shows the number of constraints and optimization methods associated with each mode.  Figure 2 depicts the variation of the indoor air temperature and the power consumed by the heat pump to maintain that temperature. During the entire simulation, the indoor temperature stays within the subjected comfort zone boundaries and the constraint on the power consumption of the heat pump is also satisfied.

Results and Discussion
Similarly, Figure 3 presents the temperature variation of water in the water tank, and the power consumed by the water heater, respectively. The initial water temperature is considered to be at 20 • C (see Figure 3a). As the time passes, the temperature gradually increases over next few hours and then it is maintained inside the specified temperature range. This also explains the fact that the power consumed by the water heater is high in beginning of the simulation period in Figure 3b. Figure 4a,b shows the power consumed by the washing machine and the dishwasher, respectively. It can be seen that the appliances operation is interrupted occasionally in order to minimize the electricity cost and satisfy the capacity constraints. It should be noted that, the HEMS has no knowledge about the arrival request of washing-machine and dishwasher. The arrival requests are generated randomly to investigate the capability of the HEMS to handle immediate request situations. In order to assess the performance of the proposed method, the following different scenarios were simulated.

Without HEMS and Linear Quadratic Regulation (LQR) Control for Thermal Appliances
In this case, it was considered that there is no HEMS to schedule the appliances, however, the thermal appliances were operated using LQR feedback control. This control strategy is a state-feedback control method in which the control input u is given by the feedback control law: u = −Kx, where x is the state of the system and K is the control gain which is obtained by solving Riccati differential equation. Figure 5a shows the total power consumption of all the appliances in this case. It can be seen in the figure that initially the power consumption violates the capacity constraints because the water heater runs at full power due to the low initial temperature of water. The maximum power consumption is at 11 kW. The operation schedule of appliances in this case is given in Figure 6. In this case, thermal appliances are almost always running, since the controller tries to track the set point at each time step. However, this could be changed by providing time-varying constraints on the set points. Clearly, we see that the non-thermal appliances are never interrupted and run continuously once started. In Figure 6, on the noon of the fourth day of the week, the water heater, washing machine, and dishwasher are operating all together at their respective maximum rated-power, which brings the total power consumption to 11 kW.

Without HEMS and MPC Control for Thermal Appliances
In this case, there also is no HEMS to schedule the appliances in this case also. But the thermal appliances are operated using MPC strategy, which is an advanced control technique in comparison to LQR control. Figure 5b shows the total power consumption in this case. It was observed that the peak power consumption in this case is reduced to around 7.5 kW which occurs at around 12:30 PM on the fourth day. It can seen in the Figure 7 that all the appliances (thermal and non-thermal appliances) are in operation at this time which implies that only 0.5 kW power is consumed by the thermal appliances. This happen because the MPC has prior knowledge of the external disturbances (MPC), so it operates the appliances to comply with the subjected constraint of maximum peak-power consumption. Clearly, using MPC instead of LQR controller is useful to reduce the peak-power consumption.

With HEMS
In this scenario, there is HEMS employed in the house to schedule the appliances. All the appliances are scheduled using MPC strategy in this case. Figure 5c clearly depicts the benefits of using HEMS and presents the total power consumed by all the household appliances. In Figure 5c, it can be seen that for initial time steps, total power consumption goes beyond the capacity limit of 4 kW, which happens due to the high power consumption of water heater since initial water temperature in the tank is low. At later time steps, it can be seen that the maximum total power consumed by all the appliances stays around the capacity limit. The constraints on the capacity limit are allowed to violate to ensure the feasibility of the optimization problem (soft constraints). Figure 8 shows the operation schedule of the appliances in this case. It can be seen that, the non-thermal appliances are interrupted occasionally to satisfy the capacity constraints and minimize the electricity cost. For example, at the noon of the fourth day of the week, washing-machine and the dishwasher is not running rather they are shifted to run later in the day.
The total electricity cost as per the time varying prices used in this study, is detailed in the Table 4. It should be noted that this includes the cost of running both thermal and non-thermal appliances for entire simulation period. In comparison to the case without HEMS and LQR controller for the thermal appliances, a reduction of 70% is achieved in total electricity cost by using the HEMS. Similarly, in comparison to the case without HEMS and MPC controller for the thermal appliances, the HEMS achieves a reduction of 45%. To see the effectiveness of HEMS in reducing the peak-load, the total violation of the capacity constraint over the simulation period and the maximum peak-load consumed by the appliances for all the cases was calculated, and is presented in Table 5. The total violation for the case when appliances are scheduled with HEMS is about 10 times lower than the case of no HEMS with MPC control for thermal appliances. Similarly, by using HEMS, the peak power consumption is reduced by 57% and 35% in comparison with the case of no HEMS with LQR control and MPC control for thermal appliances, respectively. The PAR for all the cases is detailed in the Table 6. It is observed that by using the HEMS, PAR is reduced by 35% and 24% in comparison to other two cases.

Computation Time of HEMS
The simulations were carried out in MATLAB on a computer with Intel(R) Core(TM) i7-5600 CPU with 16 GB RAM. Since, the non-thermal appliances can start at any time of the day, 10 different simulation with random starting time of non-thermal appliances are run. The computation time for each case is recorded and presented in Figure 9. The average computation time for each case came out to be 1243.6 s, i.e., 1.8 s per time step.

Conclusions
In this work, the problem of scheduling of appliances to achieve reduction in the peak-load and total electricity cost has been investigated using model predictive control (MPC). An architecture of a home energy management system (HEMS) was presented which operate in two modes of operations (MOOs) based upon the categorization of the appliances (thermal and non-thermal). The system dynamics of the appliances was modeled into a state-space formulation.
The proposed approach was able to schedule the appliances dynamically as the optimization problem is formulated and solved at each time step. The proposed framework was shown to provide flexibility to the user to turn the non-thermal appliances at any time step while achieving the objective of cost minimization. The simulation results showed the reductions of up to 70% in the total electricity cost and up to 57% in the peak power consumption, as well. The operation schedule of appliances with and without HEMS was also compared. The non-thermal appliances were able to interrupt and shift the operation to the low electricity price periods while finishing the task by the deadline.
The future work will include the consideration of energy storage facility (thermal and electric) [56,57], renewable energy generation, and integration of plug-in electric vehicle (PEVs) into the framework. In addition, the framework will be scaled up to a small community or a group of houses to investigate the demand-response from the perspective of power utility or aggregator companies.