Next Article in Journal
Commutation Torque Ripple Suppression Strategy of Brushless DC Motor Considering Back Electromotive Force Variation
Previous Article in Journal
Evaluation and Optimization of the Oil Jet Lubrication Performance for Orthogonal Face Gear Drive: Modelling, Simulation and Experimental Validation
Previous Article in Special Issue
Research on Network Intrusion Detection Based on Incremental Extreme Learning Machine and Adaptive Principal Component Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Automatic Coordination of Internet-Connected Thermostats for Power Balancing and Frequency Control in Smart Microgrids

Department of Mechanical Engineering, San Jose State University, San Jose, CA, 95192, USA
*
Author to whom correspondence should be addressed.
Energies 2019, 12(10), 1936; https://doi.org/10.3390/en12101936
Submission received: 8 February 2019 / Revised: 8 May 2019 / Accepted: 13 May 2019 / Published: 20 May 2019
(This article belongs to the Special Issue Internet of Things and Smart Environments)

Abstract

:
This paper proposes a novel feedback control strategy, so-called clock-like controller (CLC), to balance power supply and demand in smart microgrids by adjusting the setpoint temperatures of air conditioning (AC) loads. In the CLC algorithm, the grid operator communicates with the individual thermostats via the Internet and adjusts their setpoints by discrete temperature intervals (e.g., ±0.5 °C). Numerical simulations indicate that the proposed algorithm is able to deliver a smooth controllability of the aggregate AC power despite discrete temperature offsets. It can also be used for peak load shedding to mitigate the power generation cost. The CLC algorithm is then integrated into the grid frequency control problem, in which both power generators and loads in the network attempt to regulate the frequency of the system despite disturbances from demand, renewable sources, and local weather conditions. An autonomous microgrid model including a steam and a hydro generator, a solar energy source, and a large number of thermostatic loads is developed to evaluate and demonstrate the proposed method. Simulation results indicate that the AC loads with CLC algorithm can help maintain the power system frequency during extreme events when demand exceeds the maximum generation capacity available to the network. Under normal conditions, the contribution of demand-side control is marginalized by the fast responding generators, because of time delays in the frequency measurement and internet communication network.

1. Introduction

This paper addresses the problem of power balancing and frequency control in smart autonomous microgrids using the setpoint adjustment of internet-connected thermostats. Due to the large thermal capacity of residential buildings, and relative flexibility in their temperature setpoints, heating and cooling loads offer a great potential for energy management in smart grids. This paper contributes to the literature in at least three different ways. First, unlike most previous methods which adjust the setpoint temperatures of the loads concurrently and continuously, the proposed control scheme offsets the setpoints in discrete intervals, e.g., ±0.5 °C, to facilitate implementation through the currently available infrastructure, i.e., the Internet. Secondly, the proposed algorithm provides a systematic mechanism to call the loads without any random selection process, to enable fair and equal contribution of all participating loads. Last, this paper develops a grid frequency control system for thermostatic loads in the presence of active automatic generation control, and investigates scenarios in which demand-side frequency control can be effective.
Demand-side energy management has been the subject of many research studies in the past few decades. Since the emergence of smart grid, electric vehicles and thermostatically-controlled thermal loads have been at the center of attention for demand-side energy management [1,2]. Various types of thermostatic loads have been investigated in the literature, including air conditioning systems [3,4], water heaters [5,6], and refrigerators [7]. In most of these studies, setpoint temperature has been selected as the primary control variable, through which the grid operator can access the loads and alter their thermal response based on the grid needs. Renewable energy balancing by power tracking [3,4], frequency regulation [7,8,9], and peak load shedding [10,11] are the most commonly-adopted control objectives. Direct on-off scheduling using a predicted electricity price signal [12,13], energy management based on packet scheduling algorithms [14,15], and privacy preserving-based management strategies [16] have also been studied in the recent literature.
Aggregate power of thermostatic loads has been generally modeled in at least four different ways. The first approach is to simulate a large number of individual load dynamics with randomized parameters and initial conditions, and then aggregate them to replicate a realistic system [3,4,17]. This method, also known as the Monte Carlo approach, is the closest approximation to practice for validating the performance of demand-side controllers. For controller design, however, some studies have looked into developing control-oriented partial differential equation models for the distribution of loads on the thermostatic on and off states. These models are either represented in the form of diffusion [3,17,18] or transport PDE [4], and can be used directly or in the discrete form for stability analysis. The third modeling approach is to approximate the aggregated power dynamics using a linear second order transfer function. These models are developed for both setpoint temperature offset [19] and ambient temperature [20] as the external system inputs. Finally, state queuing [21,22] or stochastic Markov chain [23] models have been developed, which discretize the thermostatic temperature range in a finite number of states and capture load distribution dynamics using matrix-based state transition formula. These models are primarily developed for heterogeneous loads, and tend to provide good approximation of the Monte Carlo model.
Grid frequency control is one of the main operations targeted by demand-side energy management studies. Traditionally, power generators on the grid are responsible for correcting frequency deviations through a process called frequency droop control [24]. In multi-area power networks, automatic generation control (AGC) is used to simultaneously regulate frequency and tie-line power flows between neighboring areas through area control error (ACE) signal [24]. Under AGC, a selected number of generators are equipped with integral controllers, which gradually lift the burden from the initially responding generators, thereby allowing them to return to their scheduled setpoints. The frequency droop controllers stabilize frequency within 3–5 s upon an event. For the demand-side frequency controllers to be effective, a similar or faster response time is required. Given the time delays present in the communication infrastructure and internal load dynamics, competing with the generation-side controllers seems to be a very challenging task.
For microgrids with energy storage recourses such as electric vehicles, different frequency control methods have been proposed, such as model-free sliding mode control [25] and fuzzy logic-based optimal PID controllers [26,27]. Moreover, other load types such as household air conditioning systems and refrigerators have been adopted for grid frequency stabilization [7,9]. These methods are primarily developed around the assumption that the setpoint temperatures of the responding loads can be varied continuously and simultaneously using a centrally broadcast signal or an onboard measurement device for grid frequency. Presently, the market is being dominated by smart internet-connected thermostats under different brand names such as Honeywell, Nest or Ecobee. These thermostats are connected to a central server and can be managed by the user via the thermostat’s controls or a computer application. Moreover, the users may opt in to a demand-response program and allow the central server to temporarily control the system in exchange for an economic incentive.
The setpoint temperatures in the smart thermostats can only be shifted by discrete intervals (e.g., ±0.5 °C or ±1 °F) due to sensor limitations and simplicity of operation. Therefore, many control logics that have been developed based on continuous temperature offset strategy cannot be employed in the current infrastructure. This paper attempts to address this problem by proposing a systematic control framework for internet-connected smart thermostats, building significantly on an earlier publication [28]. The control logic is reformulated for a variety of control problems such as power tracking, peak load shedding, and frequency regulation in the presence of AGC. Several microgrid-level simulations are presented to demonstrate the strengths and deficiencies of the proposed method in different real-word scenarios.
The remainder of the paper is organized as follows: Section 2 presents the model used for simulating the aggregate power of a large number of air conditioning loads. Section 3 formulates the proposed controllers for power tracking and peak-load shedding, along with a mechanism for protecting user comfort. Section 4 formulates the frequency regulation problem in the presence of an automatic generation control system. Numerical simulations are included in all the three sections to evaluate and validate the proposed models and algorithms. Finally, Section 5 summarizes the paper’s key conclusions.

2. Aggregate AC Load Model

In this section, an aggregate AC load model is developed by modeling each AC system with a first-order thermal differential equation as follows [3]:
T ˙ i = 1 C i R i ( T T i s i ( t ) R i Q i ) , i = 1 , 2 ,   ,   N A C
where i is the load index, T and T represent the indoor and ambient temperatures, respectively, C and R are the thermal capacity and thermal resistance of the building, Q is the energy transfer rate by the AC system, and s is the thermostatic switching controller described as:
s i ( t ) = { 0 ,   i f   s i ( t ε ) = 1   a n d   T i ( t ) T m i n , i 1 , i f   s i ( t ε ) = 0   a n d   T i ( t ) T m a x , i s i ( t ε ) , o t h e r w i s e
where ε is an infinitesimal time delay to facilitate the formulation of the thermostatic switching (relay) function in continuous time. T m i n and T m a x represent the lower and upper temperature bounds of the thermostat, respectively. These quantities are related to the setpoint temperature ( T s p ) and the thermostat’s dead-band δ d b .
T m i n , i = T s p , i δ d b 2 , T m a x , i = T s p , i + δ d b 2
The aggregate AC power is the summation of individual AC powers adjusted by their respective coefficient of performance (η):
P A g g ( t ) = i = 1 N A C s i Q i η i
In this study, we assume that the grid operator is able to communicate with the individual thermostats and offset the setpoint temperatures up or down by a limited amount applied separately to each thermostat:
T s p , i ( t ) = T s p 0 , i + u i ( t ) T s t e p , i
where T s p , 0 represents the user-specified setpoint value, T s t e p is the temperature step (e.g., 0.5 °C, 1 °C, or larger), and u is a signed integer determining the number of temperature steps away from the user setpoint.
Simulink provides a simple way to simulate the system in the block diagram form as shown in Figure 1. All the parameters and signals have been vectorized to simplify the Simulink model. The AC temperatures are initialized inside the integrator, and a custom relay function is developed based on Equation (2) to allow arbitrary initial values for the switching state, s. The 1st order Euler ODE solver is used to enable simulation of discrete-time clock-like controller discussed in the next section. The parameters of the system are adopted from [3] and given in Table 1. The table also contains a column on the relative standard deviation of the parameters used to create a realistic population from normal parameter distributions.
To simulate the aggregate system response, a 24-h ambient temperature profile is created based on temperature data for the city of San Jose, California, in a warm summer day. The number of loads (NAC) is set to be 5000 to create a smooth aggregate response. The setpoint temperatures are distributed normally with 21.0 °C (~70 °F) mean value and 10% relative standard deviation. These values are then rounded to the nearest half degree value to create a discrete distribution with 0.5 °C intervals. The system is simulated first with temperatures initially set to be the same as the initial ambient temperature. The simulation is then repeated by setting the initial temperature values same as the final temperature values from the first run. This would represent a more realistic simulation, as the initial temperatures are dependent on the historical ambient temperature trajectory and the user setpoints.
Figure 2 shows the simulation result, with the ambient temperature, aggregate AC power, and five sample indoor temperature trajectories with five different setpoint values plotted. As can be seen, the aggregate power (divided by the number of loads) increases as the ambient temperature rises. However, there is a time delay of around 1–2 h between the ambient temperature peak and the aggregate AC power peak due to the thermal capacitance of the loads. During peak time, the aggregate AC power reaches to around 30% of the maximum power the ACs can possibly demand (1.6 kW vs. 5.6 kW max power per AC unit). This is because of the fact that the loads turn on and off at different times due to different initial conditions and thermal properties. For this simulation, an outdoor-indoor temperature difference of around 28 °C can push the AC to remain always on, according to Equation (1) when T ˙ = 0, and using the nominal values in Table 1.
The simulation in this section shows that the setpoint temperature impacts the AC demand. From the sample temperature trajectories in Figure 2, we can see that the total power-on duration is longer for loads with lower setpoint temperatures. This behavior motivates us to develop an algorithm that provides controllability over the aggregate AC demand through adjustment of the setpoint temperatures via internet-connected smart thermostats.

3. Aggregate Power Control using Clock-Like-Controller

Since temperature intervals are discrete in smart thermostats, the aggregate power controller cannot send temperature step commands to all AC units simultaneously. Otherwise, large power fluctuations will occur with all AC units turning off or on simultaneously. Therefore, the setpoints must be adjusted selectively at different times based on the system needs. In this study, a clock-like control (CLC) framework is introduced for the systematic management of aggregate AC power. The goal of the CLC method is to (i) provide a smooth control over the aggregate AC power, and (ii) avoid randomized selection schemes, and engage all the participating loads in a fair and equal way.
The CLC system includes a clock face (dial) whose periphery is numbered from 1 through the total number of AC units, NAC (see Figure 3). The clock has two hands that can rotate in the clock-wise direction only, and meet all of the AC units in one full turn. When one of the clock hands crosses an AC unit sector, it increases its setpoint temperature by one temperature step (e.g., 0.5 °C). The other clock hands steps down the setpoint temperature of the AC unit it crosses. We can therefore modify the aggregate AC power response by increasing and decreasing the setpoint temperatures of the individual ACs through the rotation of the two clock hands. A feedback control algorithm can provide the required angular positions of the clock hands to deliver the ultimate aggregate power tracking objective.

3.1. Feedback Control Formulation Using the CLC Method

The objective of the CLC scheme is to balance the power supply and demand within the microgrid when AC units are active. We assume that the controller is only provided with the tracking error, defined as:
e ( k ) = { P S u p ( k ) P D e m ( k ) , H e a t i n g P D e m ( k ) P S u p ( k ) , C o o l i n g
The determination of heating or cooling mode is based on the majority of the loads. It is reasonable to assume that the demand-side control is only effective during sufficiently warm or cold days, when the AC units are actively heating or cooling. When there is a mix of heating and cooling loads in moderately warm or cool days, the aggregate power will be insignificant, and therefore not useful for demand-side control.
To develop a control law for the two clock hands, we first define the positive and negative error signal as:
e p ( t ) = { e ( t ) , i f   e ( t ) > 0 0 , i f   e ( t ) 0
e n ( t ) = { e ( t ) , i f   e ( t ) < 0 0 , i f   e ( t ) 0
An integral controller is then used to rotate the clock hands based on the positive or negative error:
θ p ( t ) = K i , p 0 t e p ( τ ) d τ
θ n ( t ) = K i , n 0 t e n ( τ ) d τ
where Ki,p and Ki,n are the integral gains, determining the rotation rate of the clock hands based on the power error. These gains must be positive for the cooling, and negative for the heating loads, for the CLC to function correctly.
To link the clock hand positions to the AC units on the clock periphery, a certain time interval must be dedicated to observe the rotation of the clock hands. Therefore, the CLC scheme has to be implemented in the discrete-time form.
At a given time, t, the index of the AC unit pointed by the positive or the negative clock hand can be calculated from:
z p ( t ) = ( θ p ( t ) 2 π θ p ( t ) 2 π ) N A C
z n ( t ) = ( θ n ( t ) 2 π θ n ( t ) 2 π ) N A C
where · and · indicate the floor and ceiling operations, respectively.
The setpoint temperature change for AC unit i , ( i { 1 , , N A C } ) can be calculated from:
Δ u i ( k ) { 1 i f   z p ( k 1 ) < i z p ( k )   o r i z p ( k ) < z p ( k 1 )   o r z p ( k ) < z p ( k 1 ) < i 1 i f   z n ( k 1 ) < i z n ( k )   o r i z n ( k ) < z n ( k 1 )   o r z n ( k ) < z n ( k 1 ) < i 0 o t h e r w i s e
where Δ u i ( k ) = u i ( k ) u i ( k 1 ) , with k being the discrete time index (k = tt). Note that f(k) is the simplified version of f(kΔt), representing the discrete-time function at time t or time step k.
Equation (13) summarizes the formulation of the CLC algorithm. When the AC unit index is swept by one of the clock hands within a certain time interval Δt, the setpoint temperature for that load is changed by one temperature step up or down. The reason we have three conditions for the setpoint offsets is due to the fact that we wrap θ between 0 to 2π in the computation of AC index, z. Although z is a normally increasing signal, i.e., z ( k ) z ( k 1 ) , there are times when the value of z(k) is reset because of the clock hand crossing 2π, thereby dropping below z ( k 1 ) . In such events, one portion of the covered load sectors lie between z ( k 1 ) to NAC, and another part between 1 to z(k), adding two more conditions per clock hand to the setpoint offset event.

3.2. Accounting for End-User Comfort

To ensure positive end-user response, a certain limit must be applied to the setpoint temperature offset signals to maintain the indoor temperature near the user specified value. This limit can be imposed in the CLC scheme by limiting the difference between clock hand positions. If the clock hands are separated by 2π at a given time, all the homes have experienced one full temperature step deviation from the initial setpoint. Therefore, a limit of 2π seems to be a reasonable limit. To impose a 2π limit on the separation of the clock hands, the integral control laws can be modified to:
θ p ( t ) = K i , p 0 t Proj Δ θ [ e p ( τ ) ] d τ
θ n ( t ) = K i , n 0 t Proj Δ θ [ e n ( τ ) ] d τ
where:
Δ θ ( t ) = θ p ( t ) θ n ( t )
Proj Δ θ [ e p ( t ) ] = { 0 i f   Δ θ ( t ) 2 π   a n d   e p ( t ) > 0 e p ( t ) o t h e r w i s e
Proj Δ θ [ e n ( t ) ] = { 0 i f   Δ θ ( t ) 2 π   a n d   e n ( t ) < 0 e n ( t ) o t h e r w i s e
The role of the projection operation is to simultaneously clamp the separation of the clock hands, and prevent the integrator windup by setting the integrant to zero during the clamp period.

3.3. CLC Simulations

To demonstrate the CLC system response, a Simulink diagram is created first as shown in Figure 4. The upper diagram demonstrates the single-input, single-output feedback loop for the system, and the lower diagram shows the single-input, multi-output CLC subsystem, where the power supply and demand error is received as the input and the control vector u is calculated as the output. Two MATLAB functions are used to program the projected positive and negative errors before integration, i.e., Equations (17) and (18), and the CLC indexing and activation logic, i.e., Equations (11)–(13). Similar to the aggregate AC model, the CLC is implemented in the vector form at and after the second MATLAB function’s output. Once the control vector u is computed, it gets broadcast to all participating thermostats via internet. Each thermostat receives one element of u and responds accordingly. In practice, a time delay in the order of a few seconds may be involved in the internet communication, which is accounted for in the Simulink model.
The first simulation is to track a desired power trajectory comprising both smoothly varying and abruptly changing components as shown in Figure 5. The purpose of this simulation is to evaluate how CLC can respond to different conditions on the grid, such as smooth ramping or sudden changes due to contingencies. The ambient temperature is set to be 30 °C, and kept constant throughout the simulation. The setpoint temperature step is set to be 0.5 °C for all loads, and the communication delay is set to 5 s. To make sure the room temperatures remain near the user setpoint, the average value of desired power is matched with that of the natural aggregate power of the ACs. Mismatch in the desired vs. natural power for a long period of time may result in a continuous drift in one of the clock hands, thereby leading to control saturation.
As can be seen in Figure 5a, the aggregate AC power is able to closely track the desired power trajectory in the smooth (sinusoidal) segment. During abrupt changes, the aggregate power is able reach and stay around the desired level within about 1 min. There is a time period between hours 11:00–12:00 where the aggregate power parts way from the desired power. This is due to the implemented control saturation function, which limits clock hand separation within ±2π to enforce user comfort. As can be seen from Figure 5b, the clock hand separation angle hits the lower limit between hours 11:00–12:00 to prevent setpoint temperature migration beyond the intended ±0.5 °C level. The clock hand separation, however, remains within ±2π for most part, providing control over the aggregate power. Figure 5c shows a sample indoor temperature trajectory under the provided simulation, compared to the natural (uncontrolled) response. In this example, the natural temperature varies between 19.75 °C to 20.25 °C, whereas the controlled temperature varies between 19.25 °C to 20.75 °C, which are ±0.5 °C apart from the limits of the natural response.
To provide a more realistic simulation, the application of the proposed control method in peak load shaving scenario is studied next. We apply the ambient temperature trajectory in Figure 2, resulting in the peak aggregate power of 1.65 kW per AC unit, around 3 pm. The control objective in this simulation is to limit the peak power at or below a particular level, e.g., 1 kW per AC unit, using the CLC method. To enable this process, a modification is needed in the CLC algorithm. We only need to activate the CLC controller when the aggregate power demand goes beyond the desired level, and that happens when the tracking error is positive (in the case of cooling loads). Therefore, we can simply let en = 0 in the controller to allow natural AC operation below the desired power limit. A second modification is to bring the setpoint temperatures back to their initial values by driving the clock hand angle difference to zero once the event is cleared. A simple and creative way to enable this effect is to modify the negative error definition as follows:
e n ( t ) = { α Δ θ ( t ) , i f   e ( t ) < 0 0 , i f   e ( t ) 0   ( P e a k   S h a v i n g   &   C o o l i n g )
where α is a scaling factor, determining the convergence rate of clock hands toward each other. Based on this revised formula, the CLC will correct the separation angle when the power demand falls below the limit. To further analyze the controller dynamics when e(t) < 0, we can replace the revised en in the integral control law (Equation (10)), and take time derivative to get:
θ ˙ n ( t ) + α K i , n θ n ( t ) = α K i , n θ p ( t )
which indicates θ n θ p as t if α K i , n > 0 .
Figure 6 shows the peak load shaving simulations. The aggregate AC power for the cases of natural (uncontrolled) scenario, and controlled scenario, with and without setpoint correction, are shown in Figure 6a. The peak shaving objective is fairly achieved in both controlled cases. However, when the setpoint correction mechanism is implemented (based on Equation (19)), there is a payback period after the natural response falls below the limit. This is due to the fact that, setpoint correction requires the ACs to cool down to the initial temperature setpoints. Figure 6b compares the clock hands separation angle for the two cases, indicating an exponential converge of the angle separation to zero, as anticipated from Equation (20). The convergence of angle separation to zero is necessary for two reasons: (i) returning the occupants’ comfort, and (ii) preparing for the next day operation by bringing the control signal back to the middle position. To show the impact of the control process, a sample load temperature trajectory is shown in Figure 6c. Under the peak load shaving controller, the indoor temperature rises by half a degree above the normal thermostatic range, but is then returned to the range through the setpoint correction mechanism.

4. Microgrid Frequency Control using AGC and CLC

Frequency stability is becoming an increasingly important problem in microgrids due to the growth of disturbances introduced by uncontrollable renewable resources. Traditionally, load frequency control (LFC) and automatic generation control (AGC) methods are used to stabilize grid frequency [24]. With automatic demand-side control, additional stability margin can be introduced to the system. Understanding the dynamic interactions between generation-side and demand-side control has not been investigated in the literature in a good depth. This section attempts to provide some insights into how thermostatic loads can cooperate with AGC to improve grid frequency stability.

4.1. Automatic Generation Control Model

A well-stablished model for AGC is used in this paper to evaluate the simultaneous demand-side and generation-side frequency control problem. This model is based on the frequency droop control in power plants serving the grid in conjunction with integral controllers on selected power plants [24]. Assuming the grid voltage and the frequency dynamics are decoupled, the frequency control problem is reduced to an equivalent single-inertia system, relinquishing the need for network topology. In this paper, a hybrid generation system comprising a steam and a hydro power plant is used for the generation-side microgrid control based on models in [24], as shown in Figure 7. The steam plant has a second order dynamics, and the hydro plant has a non-minimum-phase first order transfer function. These models are the linearized versions of complex physics-based nonlinear models [24], and are suitable for stability analysis and control design. Both power plants are equipped with speed governors (controllers), which respond to changes in frequency immediately after an event (sudden change in load), by adjusting the fuel valve or water valve positions. Besides, the steam power plant is equipped with an integral controller which is responsible for bringing the frequency back to the target level (e.g., 60 Hz) after the initial response from the droop controllers. This way, the steam generator becomes the main power plant in charge of tracking slow changes in the load, and the hydro power plant provides additional bandwidth to suppress the transient disturbances.
The output from the power plants is subtracted from the load, resulting in a net power imbalance in the system. A first order linearized model is traditionally used for analyzing the grid’s inertial dynamics, and is adopted in this study as shown in Figure 7. This model includes two parameters, accounting for the collective grid rotor inertia through time constant M, and the load damping factor, through constant D. The entire generation-side system is represented in the normalized per-unit (pu) system for the simplicity of analysis and consistency with the literature. The parameters of the system used for the simulations of this section are adopted from [24] and listed in Table 2. Most of the parameters are various time constants associated with the response time of the power plants and the controllers. As can be seen, all the time constants are within 10 s, which indicates the frequency can be recovered within a few seconds after an event. In practice, there are limits on the ramping rate of the power plants, which would impact frequency recovery depending on the size of the disturbance. An improved generation-side model is used in the next section to account for these limits.
To simulate the AGC system response, a 5% step increase in load is applied to the system at t = 10 s, and response of the system is recorded for 2 min. Figure 8 shows the frequency deviation in pu, as well as the power plants responses to the step disturbance with and without integral controller on the steam plant.
As can be seen in both cases, the frequency is recovered within a few seconds. The system with integral controller is able to eliminate the steady-state error within a few minutes (with Ki,AGC value of 1). Without the integral controller, both power plants share 50% of the load at steady state due to identical frequency droop coefficients (Rh = Rs). With integral controller, however, the steam plant takes the full burden of the load at steady-state, allowing the hydro plant to go back to its initial state after the preliminary reaction.

4.2. Clock-Like Frequency Controller (CLFC) for Thermostatic Loads

To control grid frequency through internet-connect thermostatic loads, the control error signal is changed from power imbalance in CLC to frequency deviation in CLFC. Therefore, the control strategy must be changed accordingly. Due to the first order dynamic relation between the frequency and power imbalance, a modified error signal is created and used as input to the controller, combining frequency deviation and its time derivative as follows:
ξ ( t ) = { σ Δ f ˙ ( t ) + Δ f ( t ) , H e a t i n g σ Δ f ˙ ( t ) Δ f ( t ) , C o o l i n g
where σ is a positive control parameter, determining the weight of the time derivative of the frequency error in the control process. This modification builds a predictive response in the controller to take action before the frequency deviates significantly. It is also consistent with other methods proposed in the literature for grid frequency control using thermostatic loads [9].
The positive and negative error signals for frequency control are then defined as:
ξ p ( t ) = { ξ ( t ) i f   ξ ( t ) > ε f α Δ θ ( t ) i f   | ξ ( t ) | ε f   a n d   Δ θ ( t ) < 0 0 , o t h e r w i s e
ξ n ( t ) = { ξ ( t ) i f   ξ ( t ) > ε f α Δ θ ( t ) i f   | ξ ( t ) | ε f   a n d   Δ θ ( t ) < 0 0 , o t h e r w i s e
where α Δ θ ( t ) is the correction term to bring clock hands together when the magnitude of the weighted error term, ξ , is smaller than a certain threshold set by ε f . The CLFC gets activated only when the magnitude of ξ exceeds ε f . The positive and negative clock hands are then adjusted based on projected error integrals as before:
θ p ( t ) = K i , p 0 t Proj Δ θ [ ξ p ( τ ) ] d τ
θ n ( t ) = K i , n 0 t Proj Δ θ [ ξ n ( τ ) ] d τ
To demonstrate the operation of CLFC, a diagram is created that separates the three error regions for the controller, and summarizes the clock hand dynamics for each region (see Figure 9). For example, consider a case where the AC loads are in the cooling mode, and Δf and its time derivative are positive. These indicate there is more power generation than demand. In this case, ξ ( t ) is negative, and if it grows to the point where ξ ( t ) < ε f , the negative clock-hand will be activated and the setpoint temperatures will be reduced, increasing the AC power and hence, the total power demand in the network. This will drive ξ ( t ) back to the desirable region.
To demonstrate the performance of the proposed CLFC system, a simulation is created, in which the microgrid frequency is controlled using the aggregate AC power only. The block diagram of this hypothetical system is shown in Figure 10. The CLFC subsystem receives the frequency error, computes the setpoint temperature offsets, and broadcasts them to the AC thermostats with a certain time delay. To simulate the system in the pu scale, the aggregate AC power is shifted down by its initial value to start from zero, and is scaled by a conversion factor (pu_c). This factor is set to be 1/(NACPmax) in this simulation, where Pmax indicates the expected value of the AC electrical power (i.e., 5.6 kW in this paper). The performance of the CLFC is then evaluated by applying a disturbance ΔL and monitoring the frequency error signal.
The simple case of constant ambient temperature (30 °C) is considered for the simulation of this section. The disturbance signal is set to be 5% step increase in load, 1 min after the start of the simulation. For the first simulation, time delay and εf are set zero, the value of σ is set 5, and the two integrator gains are set to be equal (i.e., Ki,p = Ki,n). Figure 11 shows the response of the CLFC system for different values of the integrator gains. As can be seen from Figure 11a, the frequency recovery is accomplished successfully, with the frequency returning to the original level after a temporary departure. The larger the integral gains, the smaller the frequency deviation. The integral gain cannot be increased beyond a certain value due to risk of instability because of time delays. Figure 11b shows that the frequency recovery is achieved by moving the AC power in the opposite direction of the load disturbance, and hence, maintaining the power balance in the system. The clock hands separation angle for this simulation are shown in Figure 11c in the pu scale, with 1 pu indicating 2π radians. It is important to note that if the load disturbance is not corrected externally, the clock hand separation will drift toward 1.0 pu, resulting in control saturation and loss of performance.
To demonstrate the effects of the frequency time derivative, the same simulation is repeated by fixing the integrator gains (Ki,p = Ki,n = 2) and trying different values of σ, ranging from 0 to 20. Time delay and εf are kept at zero. Figure 12 shows the resulting frequency trajectories, underscoring impact of the derivative term in the controller’s transient performance. Although increasing σ enhances frequency error characteristics, it may increase risk of noise amplification. In the simulations of this paper, the value of σ is set to 5 to create a reasonable balance between performance and noise sensitivity in the system.
The next simulation is to evaluate the impacts of time delay. Since internet communication is not delay-free, it is important to evaluate the impacts of different communication delays on the controller performance. Figure 13 shows the performance of the system for different time delay values ranging from 0–10 s. As can be seen, the time delay results in performance drop and possible instability. In practice, an upper bound must be characterized for the time delay, and used to tune the appropriate control gains. The larger the time delay, the smaller must be the CLFC integrator gains. A time delay of 4 s is used in the simulations of this paper, hereafter.
The last simulation of this section investigates the impact of parameter εf, the activation threshold for CLFC. Figure 14 shows the response of the system with εf = 0.01 pu compared to εf = 0. For the nonzero εf, the grid frequency remains within ±εf, except during the disturbance period around t = 20 min. The maximum frequency deviation is, however, similar to the one with εf = 0. It is important to note that in the CLFC formulation, εf represents the upper bound of modified error signal (ξ). Figure 14a shows that this parameter can also serve as an upper bound for the frequency error. The clock hands separation angle in Figure 14b indicates that increasing the value of εf from 0 to 0.01 reduces the sensitivity of the CLFC considerably, thereby enhancing user comfort without compromising the control performance during fault periods.

4.3. Simultaneous Generation and Demand Control

In this section, the demand-side CLFC algorithm is integrated with AGC after some modifications in the AGC system. First, rate limits of 2% and 100% per minute are applied to the steam and the hydro plants, respectively, based on the values suggested in [24]. Moreover, capacity limits of 0.8 pu and 0.2 pu are applied to the steam and the hydro plants to reasonably scale the two power plants and keep their output limited to their maximum value. To further improve the frequency dynamics, a second integrator is added to the steam plant to maintain the integral of the frequency near zero, for maintaining the grid’s clock cycle accuracy. Finally, to avoid integration wind up in the event of rate or capacity saturation, an integrator anti-windup algorithm is applied to both integrators. This logic compares the control signal with the rate and saturation limits, and stops integration during periods when the limits are activated. The full power system diagram with both CLFC and modified AGC system is shown Figure 15. The ambient temperature and the user setpoint inputs are included inside the aggregate AC model subsystem. In addition, a renewable source is included in the simulation.
The initial simulation is for the simple case of constant ambient temperature (30 °C), and a 5% step increase in the demand. The renewable source is set to zero for this simulation. Figure 16a shows the response of the system for four different scenarios: steam plant only, steam and CLFC, steam and hydro, and all three systems together.
When only the steam plant controller is active, the frequency recovery is very slow due to the 2% per min rate limit on the steam plant. In this case, the inclusion of CLFC helps to improve frequency recovery. However, when the hydro plant is included in the control loop, the benefit of CLFC is diminished, since the hydro plant responds quicker than both steam plant and CLFC due to faster ramp rate (100% per minute). Figure 16b shows the outputs of the power plants, the aggregate AC power, as well as the injected load disturbance for the case of having all three systems included in the control loop. The hydro power plant responds to the load immediately, and then drops gradually as the steam plant overtakes the load. The aggregate AC power truncates about a quarter of the load disturbance and maintains its state within the 10 min simulation interval. In the extended simulation, however, the aggregate AC power returns to its initial state through the setpoint correction mechanism built in the CLFC algorithm.
The previous simulation reveals that the CLFC system can be outpaced by the hydro plant under normal conditions. In the next section, a scenario is investigated, in which the power demand exceeds the maximum generation capacity of the microgrid. This can happen on a very hot day, when the AC units work excessively to maintain desirable indoor temperatures. To provide a realistic simulation, the ambient temperature shown in Figure 2 is shifted up by 5 °C to resemble a hot day with 22 °C low and 37 °C peak temperature. The California Independent System Operator (CAISO) data is used to obtain sample grid demand and renewable power profiles for the simulation. The data for the same day as the ambient temperature data is used for consistency. To separate the aggregate AC power from the total grid demand, the true ambient temperature is used first to obtain and subtract the aggregate AC power from the demand. To covert the power data to the pu scale, a margin of 2% is assumed for the peak demand below the generation capacity (combined steam and hydro). The net demand, defined as the total power demand minus renewable power, has a much larger clearance margin from the generation capacity (around 30%).
Figure 17 shows the resulting power trajectories. The peak aggregate AC power is around 30% of the total grid demand for the true ambient temperature. This value would rise to around 59% when the ambient temperature is offset by 5 °C. The base demand (dashed line in Figure 17), defined as the total grid demand minus the aggregate AC power, is used as the demand input in the simulation, with aggregate AC power being calculated via CLFC, and injected to the system separately.
It is important to note that all the power inputs (AC, base demand, and renewable) are offset to start from zero in the simulation to avoid initial imbalance in the system. The steam and hydro plant saturation limits are offset as well to account for these initial conditions. The hydro plant is assumed to be in its mid-point, initially, moving the saturation limits from [0 0.2] to [−0.1 0.1], and the rest of the load minus the renewable generation is initially picked up by the steam plant.
Figure 18 shows the simulation results with and without CLFC, for the case of elevated ambient temperature. As can be seen in Figure 18a, in the absence of CLFC, the frequency deviates significantly during certain time periods. The reason is that the net demand (total power demand minus renewable power) exceeds the maximum generation capacity as shown in Figure 18b, resulting in the saturation of both steam and hydro power plants. The integrator anti-windup plays a crucial role in the timely recovery of frequency after the net demand returns to the trackable range. To maintain the clock cycle accuracy, there is an over-frequency period (between Hours 20–30), during which the integral of frequency is corrected to zero. The integrator gain for this purpose, i.e., Ki,G2, is set to 0.0001. Note that, the simulation time is extended to 36 h (repeating the first 12 h) for a better view of the system responses.
When the CLFC mechanism is active, the frequency deviation is maintained within ±εf (0.005 pu). This is achieved by curtailing the aggregate AC power such that the net grid demand remains below the generation capacity, as seen in Figure 18c. The clock hands separation angle in Figure 18d shows how CLFC responds to the frequency error. When the error is within the specified limit, the frequency control is accomplished by the steam and the hydro plants, and the clock hands are locked unless separation angle is nonzero. When the frequency approaches the threshold limits due to AGC saturation, the CLFC takes action automatically to contain the frequency within the specified range. Clock hand separation angle returns slowly to zero when |ξ| < εf, via the setpoint correction mechanism. The value of α is set to 0.0001/2π in this simulation.
Evaluating the total simulation time for the simultaneous AGC and CLFC system with 5000 AC units indicates that the proposed CLFC method is a computationally-friendly algorithm. It takes around 3 milliseconds per simulation time step to compute both control signals and the evolution of AC temperatures using a computer with 8 GB memory and an Intel Core i7 processor. The proposed control method can therefore accommodate millions of AC systems for large-scale grid-wide updates every 2–4 s. The main limitation remains to be the internet communication bandwidth and delays associated with broadcasting control signals to the connected thermostats.
The simulations of this section underscore the importance of demand-side energy management in resource limited microgrids when the conventional generation system reaches its maximum capacity. This research can be extended to investigate other scenarios such as presence of energy storage, economic load dispatch, and larger power grids with multi-area networks.

5. Conclusions

This paper proposed a novel control method for demand-side energy management in microgrids through internet-connected smart thermostats. The advantage of the proposed clock-like control framework is providing smooth power demand control by offsetting setpoint temperatures of loads at discrete intervals, within the capability of thermostats, in a systematic (non-random) manner. The formulation of the clock-like controller for power tracking and frequency control in microgrids was presented and evaluated through various simulations. The proposed controller is able to achieve tracking of desired power trajectories, as well as limiting the peak demand under a certain threshold, with a reasonable accuracy. The frequency control simulations indicated that proposed controller is able to stabilize the grid frequency by itself and in cooperation with generation control. In the presence of fast power plants, however, the impact of demand-side controller is diminished due to time delays in the communication network as well as the internal dynamics of the loads. The demand-side frequency controller provides an additional security margin in scenarios where the power demand exceeds the maximum generation capacity of the microgrid due to unusual weather conditions and contingences. A robust internet communication network between the grid operator and the thermostats, with a predictable time delay upper bound, is a key to the successful tuning and implementation of the proposed control method.

Author Contributions

Conceptualization, S.B., methodology, S.B., K.L.L.; software, S.B., K.L.L.; validation, S.B., K.L.L.; formal analysis, S.B., K.L.L.; investigation, S.B., K.L.L.; writing—original draft preparation, S.B.; writing—review and editing, S.B.; supervision, S.B.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

SymbolDescriptionSymbolDescription
iAC unit indexTmaxUpper thermostatic temperature limit
kDiscrete time indexTminLower thermostatic temperature limit
nNegative variable indexTspSetpoint temperature
pPositive variable indexTsp0User-specified setpoint temperature
ePower tracking errorTstepSetpoint temperature step
s Thermostatic switching stateTAmbient temperature
tTimeαConvergence rate of clock hands
zAC unit index pointed by clock hand σControl parameter
CThermal capacitanceθClock hand angle
KiIntegral control gainεfFrequency error threshold
NACTotal number of AC unitsξModified frequency error signal
PDemPower demand δ d b Thermostatic deadband temperature
PSupPower supplyηCoefficient of performance
QEnergy transfer rateΔfGrid frequency error
R Thermal resistanceΔtTime step
TIndoor room temperatureΔuSetpoint temperature offset coefficient

References

  1. Wang, Z.; Tang, Y.; Chen, X.; Men, X.; Cao, J.; Wang, H. Optimized daily dispatching strategy of building integrated energy systems considering vehicle to grid technology and room temperature control. Energies 2018, 11, 1287. [Google Scholar] [CrossRef]
  2. Callaway, D.S.; Hiskens, I.A. Achieving controllability of electric loads. Proc. IEEE 2011, 99, 184–199. [Google Scholar] [CrossRef]
  3. Callaway, D.S. Tapping the energy storage potential in electric loads to deliver load following and regulation, with application to wind energy. Energy Convers. Manag. 2009, 50, 1389–1400. [Google Scholar] [CrossRef]
  4. Bashash, S.; Fathy, H.K. Modeling and control of aggregate air conditioning loads for robust renewable power management. IEEE Trans. Control Syst. Technol. 2013, 21, 1318–1327. [Google Scholar] [CrossRef]
  5. Gustafson, M.W.; Baylor, J.S.; Epstein, G. Direct water heater load control. Estimating program effectiveness using an engineering model. IEEE Trans. Power Syst. 1993, 8, 137–143. [Google Scholar] [CrossRef]
  6. Ericson, T. Direct load control of residential water heaters. Energy Policy 2009, 37, 3502–3512. [Google Scholar] [CrossRef] [Green Version]
  7. Short, J.; Infield, D.G.; Freris, L.L. Stabilization of grid frequency through dynamic demand control. IEEE Trans. Power Syst. 2007, 22, 1284–1293. [Google Scholar] [CrossRef]
  8. Olama, M.M.; Kuruganti, T.; Nutaro, J.; Dong, J. Coordination and control of building HVAC systems to provide frequency regulation to the electric grid. Energies 2018, 11, 1852. [Google Scholar] [CrossRef]
  9. Bashash, S.; Fathy, H.K. Power grid stabilization through setpoint temperature control of frequency-responsive air conditioning loads. In Proceedings of the 2012 ASME Dynamic Systems and Control Conference, Fort Lauderdale, FL, USA, 17–19 October 2012. [Google Scholar]
  10. Risbeck, M.J.; Maravelias, C.T.; Rawlings, J.B.; Turney, R.D. Cost optimization of combined building heating/cooling equipment via mixed-integer linear programming. In Proceedings of the 2015 American Control Conference, Chicago, IL, USA, 1–3 July 2015; pp. 1689–1694. [Google Scholar]
  11. Werminski, S.; Jarnut, M.; Benysek, G.; Bojarski, J. Demand side management using DADR automation in the peak load reduction. Renew. Sustain. Energy Rev. 2017, 67, 998–1007. [Google Scholar] [CrossRef]
  12. Chan, K.; Bashash, S. Modeling and energy cost optimization of air conditioning loads in smart grid environments. In Proceedings of the 2017 Dynamic Systems and Control Conference, Tysons, VA, USA, 11–13 October 2017. [Google Scholar]
  13. Bhattacharya, S.; Kar, K.; Chow, J.H. Economic operation of thermostatic loads under time varying prices: An optimal control approach. IEEE Trans. Sustain. Energy 2018. [Google Scholar] [CrossRef]
  14. Almassalkhi, M.; Frolik, J.; Hines, P. Packetized energy management: Asynchronous and anonymous coordination of thermostatically controlled loads. In Proceedings of the 2017 American Control Conference, Seattle, WA, USA, 24–26 May 2017; pp. 1431–1437. [Google Scholar]
  15. Espinosa, L.A.D.; Almassalkhi, M.; Hines, P.; Frolik, J. System properties of packetized energy management for aggregated diverse resources. In Proceedings of the IEEE Power Systems Computation Conference, Dublin, Ireland, 11–15 June 2018. [Google Scholar]
  16. Halder, A.; Geng, X.; Kumar, P.R.; Xie, L. Architecture and algorithms for privacy preserving thermal inertial load management by a load serving entity. IEEE Trans. Power Syst. 2017, 32, 3275–3286. [Google Scholar] [CrossRef]
  17. Malhamé, R.; Chong, C.Y. Electric load model synthesis by diffusion approximation of a high-order hybrid-state stochastic system. IEEE Trans. Autom. Control 1985, 30, 854–860. [Google Scholar] [CrossRef]
  18. Moura, S.; Ruiz, V.; Bendsten, J. Modeling heterogeneous populations of thermostatically controlled loads using diffusion-advection PDEs. In Proceedings of the 2013 ASME Dynamic Systems and Control Conference, Palo Alto, CA, USA, 21–23 October 2013. [Google Scholar]
  19. Perfumo, C.; Kofman, E.; Braslavsky, J.H.; Ward, J.K. Load management: Model-based control of aggregate power for populations of thermostatically controlled loads. Energy Convers. Manag. 2012, 55, 36–48. [Google Scholar] [CrossRef] [Green Version]
  20. Mahdavi, N.; Braslavsky, J.H.; Perfumo, C. Mapping effect of ambient temperature on the power demand of populations of air conditioners. IEEE Trans. Smart Grid 2016, 9, 1540–1550. [Google Scholar] [CrossRef]
  21. Lu, N.; Chassin, D.P.; Widergren, S.E. Modeling uncertainties in aggregated thermostatically controlled loads using a state queuing model. IEEE Trans. Power Syst. 2005, 20, 725–733. [Google Scholar] [CrossRef]
  22. Chen, Y.B.P.; Zhu, X.; Hu, M. The extended 2-dimensional state-queuing model for the thermostatically controlled loads. Electr. Power Energy Syst. 2019, 105, 323–329. [Google Scholar]
  23. Mathieu, J.L.; Koch, S.; Callaway, D.S. State estimation and control of electric loads to manage real-time energy imbalance. IEEE Trans. Power Syst. 2013, 28, 430–440. [Google Scholar] [CrossRef]
  24. Kundur, P. Power System Stability and Control; McGraw-Hill: New York, NY, USA, 1993; pp. 581–587. [Google Scholar]
  25. Khooban, M.H. Secondary load frequency control of time-delay stand-alone microgrids with electric vehicles. IEEE Trans. Ind. Electron. 2018, 65, 7416–7422. [Google Scholar] [CrossRef]
  26. Khooban, M.H.; Niknam, T.; Shasadeghi, M.; Dragicevic, T.; Blaabjerg, F. Load frequency control in microgrids based on a stochastic noninteger controller. IEEE Trans. Sustain. Energy 2018, 9, 853–861. [Google Scholar] [CrossRef]
  27. Gheisarnejad, M.; Khooban, M.H.; Dragicevic, T. The future 5G network based secondary load frequency control in maritime microgrids. IEEE J. Emerg. Sel. Top. Power Electron. 2019. [Google Scholar] [CrossRef]
  28. Chekan, J.A.; Bashash, S. IoT-oriented demand-side energy management of thermostatically-controlled loads. In Proceedings of the 2017 American Control Conference, Seattle, WA, USA, 24–26 May 2017. [Google Scholar]
Figure 1. Simulink diagram of the aggregate AC system, with inputs and outputs highlighted.
Figure 1. Simulink diagram of the aggregate AC system, with inputs and outputs highlighted.
Energies 12 01936 g001
Figure 2. Simulation of 5000 AC loads: Ambient temperature, aggregate AC power, and five sample individual indoor temperatures with setpoint values between 18–22 °C.
Figure 2. Simulation of 5000 AC loads: Ambient temperature, aggregate AC power, and five sample individual indoor temperatures with setpoint values between 18–22 °C.
Energies 12 01936 g002
Figure 3. Schematic representation of the Clock-Like-Controller (CLC) mechanism.
Figure 3. Schematic representation of the Clock-Like-Controller (CLC) mechanism.
Energies 12 01936 g003
Figure 4. The Simulink block diagram of the CLC algorithm.
Figure 4. The Simulink block diagram of the CLC algorithm.
Energies 12 01936 g004
Figure 5. Sample power tracking simulation: (a) Power tracking performance, (b) Clock angle separation, and (c) Sample load temperature.
Figure 5. Sample power tracking simulation: (a) Power tracking performance, (b) Clock angle separation, and (c) Sample load temperature.
Energies 12 01936 g005
Figure 6. Sample load shedding simulation: (a) Aggregate power response, (b) Clock angle separation, and (c) Sample load temperature.
Figure 6. Sample load shedding simulation: (a) Aggregate power response, (b) Clock angle separation, and (c) Sample load temperature.
Energies 12 01936 g006
Figure 7. Block diagram of the microgrid AGC system with a steam and a hydro power plant.
Figure 7. Block diagram of the microgrid AGC system with a steam and a hydro power plant.
Energies 12 01936 g007
Figure 8. AGC’s transient response to 5% step increase in load: (a) Frequency and (b) power plant outputs, with and without integral controller on the steam plant.
Figure 8. AGC’s transient response to 5% step increase in load: (a) Frequency and (b) power plant outputs, with and without integral controller on the steam plant.
Energies 12 01936 g008
Figure 9. Diagram of CLFC dynamics.
Figure 9. Diagram of CLFC dynamics.
Energies 12 01936 g009
Figure 10. Block diagram of the CLFC system.
Figure 10. Block diagram of the CLFC system.
Energies 12 01936 g010
Figure 11. CLFC response to step load: (a) Frequency (b) aggregate AC power, and (c) clock hands separation angle, for different values of integrator gains, with σ = 5, εf = 0, and zero time delay.
Figure 11. CLFC response to step load: (a) Frequency (b) aggregate AC power, and (c) clock hands separation angle, for different values of integrator gains, with σ = 5, εf = 0, and zero time delay.
Energies 12 01936 g011
Figure 12. CLFC response to step load for different values of σ, with Ki,p = Ki,n = 2, εf = 0, and zero time delay: (a) Frequency, and (b) clock hand separation angle.
Figure 12. CLFC response to step load for different values of σ, with Ki,p = Ki,n = 2, εf = 0, and zero time delay: (a) Frequency, and (b) clock hand separation angle.
Energies 12 01936 g012
Figure 13. CLFC response to step load for different time delay values, with σ = 5, Ki,p = Ki,n= 2, εf = 0: (a) Frequency, and (b) clock hand separation angle.
Figure 13. CLFC response to step load for different time delay values, with σ = 5, Ki,p = Ki,n= 2, εf = 0: (a) Frequency, and (b) clock hand separation angle.
Energies 12 01936 g013
Figure 14. CLFC response to step load for εf = 0 and εf = 0.01, with σ = 5, Ki,p = Ki,n= 2, and 4 s time delay: (a) Frequency, and (b) clock hands separation.
Figure 14. CLFC response to step load for εf = 0 and εf = 0.01, with σ = 5, Ki,p = Ki,n= 2, and 4 s time delay: (a) Frequency, and (b) clock hands separation.
Energies 12 01936 g014
Figure 15. Combined generation-side and demand-side power system control diagram.
Figure 15. Combined generation-side and demand-side power system control diagram.
Energies 12 01936 g015
Figure 16. Simultaneous generation-side and demand-side frequency control for step increase in load, and constant ambient temperature: (a) Frequency, and (b) power outputs.
Figure 16. Simultaneous generation-side and demand-side frequency control for step increase in load, and constant ambient temperature: (a) Frequency, and (b) power outputs.
Energies 12 01936 g016
Figure 17. The total demand and renewable power trajectories obtained from CAISO website, as well as aggregate AC power, for computing input trajectories for the simulation.
Figure 17. The total demand and renewable power trajectories obtained from CAISO website, as well as aggregate AC power, for computing input trajectories for the simulation.
Energies 12 01936 g017
Figure 18. Simultaneous generation-side and demand-side control simulation for time-varying demand, renewable power, and ambient temperature inputs: (a) Frequency trajectory with and without CLFC, (b) power trajectories with AGC only, (c) power trajectories with AGC and CLFC, and (d) clock-hand separation angle (with active AGC and CLFC).
Figure 18. Simultaneous generation-side and demand-side control simulation for time-varying demand, renewable power, and ambient temperature inputs: (a) Frequency trajectory with and without CLFC, (b) power trajectories with AGC only, (c) power trajectories with AGC and CLFC, and (d) clock-hand separation angle (with active AGC and CLFC).
Energies 12 01936 g018
Table 1. System parameters adopted from [3].
Table 1. System parameters adopted from [3].
ParameterMean ValueUnitRel. Stand. Deviation
R, Thermal Resistance2°C/kW0.1
C, Thermal Capacitance10kWh/°C0.1
η, Coefficient of Performance 2.5 -0.0
δdb, Thermostat Deadband0.5°C0.0
Q, Energy Transfer Rate14kW0.1
Table 2. Generation-side system parameters adopted from [24].
Table 2. Generation-side system parameters adopted from [24].
ParameterValueParameterValue
M10.0 sD1.0
Rs0.05Rh0.05
TG0.2 sTG,h0.2 s
FHP0.3RT0.38
TRH7.0 sTR5.0 s
TCH0.3 sTW1.0 s

Share and Cite

MDPI and ACS Style

Bashash, S.; Lee, K.L. Automatic Coordination of Internet-Connected Thermostats for Power Balancing and Frequency Control in Smart Microgrids. Energies 2019, 12, 1936. https://doi.org/10.3390/en12101936

AMA Style

Bashash S, Lee KL. Automatic Coordination of Internet-Connected Thermostats for Power Balancing and Frequency Control in Smart Microgrids. Energies. 2019; 12(10):1936. https://doi.org/10.3390/en12101936

Chicago/Turabian Style

Bashash, Saeid, and Kai Lun Lee. 2019. "Automatic Coordination of Internet-Connected Thermostats for Power Balancing and Frequency Control in Smart Microgrids" Energies 12, no. 10: 1936. https://doi.org/10.3390/en12101936

APA Style

Bashash, S., & Lee, K. L. (2019). Automatic Coordination of Internet-Connected Thermostats for Power Balancing and Frequency Control in Smart Microgrids. Energies, 12(10), 1936. https://doi.org/10.3390/en12101936

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop