Next Article in Journal
An Empirical Test on Harrod’s Open Economy Dynamics
Previous Article in Journal
Approximation of Fixed Points for Suzuki’s Generalized Non-Expansive Mappings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Controlling Power Consumption in a Heterogeneous Population Model of TCLs with Diffusion: The Green’s Function Approach

by
Md Musabbir Hossain
1,* and
Asatur Zh. Khurshudyan
2,3
1
School of Mechatronic Engineering and Automation, Shanghai University, Shanghai 200444, China
2
Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, China
3
Department on Dynamics of Deformable Systems and Coupled Fields, Institute of Mechanics, National Academy of Sciences of Armenia, Yerevan 0019, Armenia
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(6), 523; https://doi.org/10.3390/math7060523
Submission received: 29 April 2019 / Revised: 30 May 2019 / Accepted: 5 June 2019 / Published: 8 June 2019

Abstract

:
We consider a control problem for a diffusive PDE model of heterogeneous population of thermostatically controlled loads (TCLs) aiming to balance the aggregate power consumption within a given amount of time. Using the Green’s function approach, the problem is formulated as an approximate controllability problem for a residue depending on control parameters nonlinearly. A sufficient condition for approximate controllability is derived in terms of initial temperature distribution, operation time of TCLs and threshold value of the aggregate power consumption. Case studies allow to reveal the advantages of the proposed solution from numerical calculations point of view.

1. Introduction

The usage of renewable energy sources is becoming more efficient as an alternative to power sources. Such advancement is exciting but performs significant challenges to the power industry. One of the main challenges is inherent variability in achieving significant penetration of renewable energy in modern power supply system [1]. Hence, a necessity of efficient power supply strategies, as well as balancing schemes in presence of intermittent power source occurs. One of best solutions is the elaboration of programmable calculating devices and power control means. Usually, thermostatically controlled loads (TCLs) are the typical control means in such systems and it may be a better option to provide necessary generation-balancing services. Because it is feasible for loads to compensate for power imbalances more rapidly than thermal generators [2]. By controlling aggregated TCLs, a promising opportunity can be obtained to mitigate the mismatch between power generation and demand, thus enabling renewable energy penetration and improving grid reliability [3].
TCLs account for more than one-third of electricity consumption in the United States [4]. Nowadays, the model of TCLs has gained extensive attention. Several methods have been used to model TCL populations. In [5], multi-agent reinforcement learning was developed for modeling of TCLs. Model accuracy scales linearly with the number of agents and gives evidence for the increased agency to further sensing, domain knowledge or data gathering time. In [6], multi-state operating reserve model of aggregate TCLs has been presented for power system short-term reliability evaluation. This analytical model for characterizing dynamics of the operating reserve is firstly developed based on the migration of TCLs’ room temperature inside temperature hysteresis band. A frequency control technique is used in the literature [7] based on decision tree concept by conducting TCLs at smart grids. A new controlling action is introduced here to limit the probabilities of separating loaded feeders from smart grids due to maloperation of under-frequency load shedding relays. In [8], authors propose a two-dimensional state-queuing modeling approach for the TCLs that increases the state vector to a two-dimensional state matrix. This method improves the accuracy of the state-queuing model by adding temperature-varying-rate information.
Loads are usually regulated by switches between on and off regimes, as it is the case for example air conditioners. Switches are used when the thermostat temperature approaches the minimal or maximal admissible values. Thus, large populations of TCLs can be manipulated by small deviations of set-point temperature without causing any considerable inconvenience to consumers. In particular, a common temperature setpoint offset needs to be manipulated for controlling the aggregate power requirement of a population of TCLs [9]. This opportunity allows to develop efficient elaboration schemes aiming to reduce power consumption by altering the set-point temperature and the minimal and maximal admissible values of the temperature.
Currently, there exist multiple mathematical models allowing to describe the basic principle of TCL units. The simplest models are the hybrid ODE models that have been investigated in many previous works including [9,10,11]. In general statement, the hybrid ODE model provides a relationship for the conditioned mass temperature, Θ i as follows [11]:
Θ ˙ i = 1 R i C i [ Θ , i t Θ i t S i t R i P i + w t ] , i = 1 , 2 , , N , t > 0 ,
where Θ , i t is the external temperature, C i R + is the thermal capacitance (kWh/ C), R i R + is the thermal resistance ( C/kW), and P i R (kW) is the cooling power supplied by the TCLs when switched on, the discrete function S i { 0 , 1 } represents the hardware on/off state, and w represents unpredictable heat gains or losses due to occupancy: all the above quantities with a subscript i are for the ith load.
Except differential Equation (1), Θ i is subject to a natural constraint
Θ m i n , i Θ i t Θ m a x , i , i = 1 , 2 , , N , t > 0 .
Here, Θ m i n , i and Θ m a x , i are the minimal and maximal admissible temperatures. Then, the set-point temperature of the ith load, Θ s p , i , is determined in terms of Θ m i n , i and Θ m a x , i as follows:
Θ s p , i = Θ m a x , i + Θ m i n , i 2 .
On the other hand, introducing the width of the temperature band of the ith load, σ i , we have
Θ m i n , i = Θ s p , i σ i 2 , Θ m a x , i = Θ s p , i + σ i 2 .
More sophisticated, but more realistic models are described by PDEs. Let the states u and v denote the density of TCL units per temperature Θ and time t in the on and off states, respectively. Assume a homogeneous population of TCLs, i.e., assume that the parameters R i , C i , P i , Θ m i n , i , and Θ m a x , i are equal for all TCLs omitting the subscript i from the corresponding symbols. Then, the couple u , v satisfies the following system of first order PDEs [3,12]:
t α λ Θ , t Θ α u Θ , t = 0 , t + α μ Θ , t Θ α v Θ , t = 0
λ Θ m a x , t u Θ m a x , t = μ Θ m a x , t v Θ m a x , t
μ Θ m i n , t v Θ m i n , t = λ Θ m i n , t u Θ m i n , t ,
where
α = 1 R C > 0 , λ Θ , t = μ Θ , t R P > 0 , μ Θ , t = Θ t Θ > 0 .
Remark 1.
Despite the fact that system (4) is coupled via boundary conditions (5), (6), there exist various methods for exact and numerical determination of u , v .
At this, some important characteristics as numbers of loads in the on and off states and aggregate power consumption are evaluated in terms of u , v as follows:
n o n t = Θ m i n Θ m a x u Θ , t d Θ , n o f f t = Θ m i n Θ m a x v Θ , t d Θ ,
κ t = P k Θ m i n Θ m a x u Θ , t d Θ ,
where k is the performance coefficient.
The quantities Θ s p , i and σ i are usually considered as control parameters in order to control the aggregate power consumption (7). For a relevant study, see, e.g., [4].

2. A PDE Model Considering Diffusion Phenomenon

Numerical simulation of various TCL populations reveals a diffusive phenomenon. Therefore, in [3] a diffusion equation has been proposed for a proper analysis of heterogeneous population models. The model equations read as
D λ u = 0 , D μ v = 0 , t > 0 , Θ Θ m i n , Θ m a x ,
where
D λ t β Θ Θ α λ Θ , t Θ α , D μ t β Θ Θ + α μ Θ , t Θ α ,
and β is a diffusivity coefficient with unit ( C) 3 /s.
Homogeneous system (8) is equipped with the following boundary conditions:
λ u μ v = 0 , Θ = Θ m i n , Θ m a x , t > 0 ,
Θ u + Θ v = 0 , Θ = Θ m i n , Θ m a x , t > 0 .
The following initial conditions are also given:
u Θ , 0 = u 0 Θ , v Θ , 0 = v 0 Θ .
Remark 2.
Note that system (8) is coupled via boundary conditions (9) and (10), ensuring that the total number of TCLs is conserved. Unlike the model without a diffusive term (see Remark 1), analysis of (8)–(11) is sophisticated.

2.1. Solution via Matrices of Green’s Type

Let us derive the closed form solution of (8)–(11) by means of the concept of Green’s function. To this end, we introduce new functions
u ¯ t = θ t 0 u t , v ¯ t = θ t 0 v t ,
where θ is the Heaviside function, and t 0 in its argument indicates that the point t = 0 is approached from the right. Then, (8) is reduced to
D λ u ¯ = u 0 Θ δ t , D μ v ¯ = v 0 Θ δ t , t > 0 , Θ Θ m i n , Θ m a x .
Boundary conditions (9), (10) refrain their forms.
Following to [13], the general solution of (11) and (12) can be represented as follows:
u ¯ v ¯ = 0 t Θ m i n Θ m a x G 11 G 12 G 21 G 22 u 0 ϑ δ τ v 0 ϑ δ τ d ϑ d τ = = Θ m i n Θ m a x G 11 G 12 G 21 G 22 | τ = 0 u 0 ϑ v 0 ϑ d ϑ ,
where G i j = G i j Θ , ϑ , t , τ , i , j = 1 , 2 , are the elements of the matrix of Green’s type. An efficient numerical scheme for computing G i j has been developed by Melnikov in [13].
Remark 3.
Note that solution (13) is linear with respect to the initial state u 0 , v 0 .

2.2. Aggregate Power Consumption

The closed-form solution (13) allows to express the aggregated power consumption in terms of system parameters. Substituting the firs raw of (13) into (7), we will obtain
κ t = P k Θ m i n Θ m a x Θ m i n Θ m a x u 0 ϑ G 11 Θ , ϑ , t , 0 + v 0 ϑ G 12 Θ , ϑ , t , 0 d ϑ d Θ .

3. Control of the Diffusive PDE Model

In this section, we consider the problem of controlling the amount of the aggregate power consumption (7) by choosing the quantities Θ s p and σ accordingly. Such a problem for the simplest model described by (4) has been considered, e.g., in [4]. The problem is to choose appropriate values for Θ s p and σ in order to ensure
R Θ s p , σ = max 0 t T κ t κ 0 ε ,
with a desired accuracy ε > 0 . Here, κ 0 0 is a given threshold value of admissible consumption, T is the operation time of the TCLs. In other words, we choose Θ s p and σ in such a way that the total power consumption during the operation interval 0 , T does not exceed a desired amount.
Remark 4.
Even though we call ε an accuracy, it does not imply that it should be a small quantity. Indeed, as soon as κ 0 = 0 , ε can not be smaller than the minimal possible value of κ in 0 , T . On the other hand, when κ 0 > 0 , meaningful values of ε 1 .
Remark 5.
Note that in the terminology of [14], this is an approximate controllability problem.
Since the mathematical model is linear in u , v , for our purpose, we involve the Green’s function approach developed in [14]. Substituting (14) into (15), we make the dependence R = R Θ s p , σ explicit.
Remark 6.
It should be noted that R depends on Θ s p , σ through the limits of integration, as well as through G 11 and G 12 . Thus, the minimization of R is not that straightforward.
The main result of the paper can be stated in the following theorem.
Theorem 1.
If for given u 0 , v 0 , T, κ 0 and ε, there exist a bounded constant C Θ s p , σ > 0 such that
u 0 ϑ G 11 Θ , ϑ , t , 0 + v 0 ϑ G 12 Θ , ϑ , t , 0 C Θ s p , σ ,
uniformly for 0 t T , Θ m i n Θ , ϑ Θ m a x , then
R Θ s p , σ P k σ 2 C Θ s p , σ κ 0 .
From Theorem 1 it directly follows
Corollary 1.
For approximate controllability of (8)–(11) it is sufficient that for given u 0 , v 0 , T, κ 0 and ε,
P k σ 2 C Θ s p , σ κ 0 ε .
Remark 7.
Depending on values of κ 0 , (18) may not be satisfied by any values of Θ s p , σ with a given ε. See also Remark 4. In such cases, inequality (15) must be evaluated.
Proof. 
When (16) holds, then the estimate
P k Θ m i n Θ m a x Θ m i n Θ m a x u 0 ϑ G 11 Θ , ϑ , t , 0 + v 0 ϑ G 12 Θ , ϑ , t , 0 d Θ d ϑ P k Θ m i n Θ m a x Θ m i n Θ m a x u 0 ϑ G 11 Θ , ϑ , t , 0 + v 0 ϑ G 12 Θ , ϑ , t , 0 d Θ d ϑ P k C Θ s p , σ Θ m i n Θ m a x Θ m i n Θ m a x d Θ d ϑ P k C Θ s p , σ Θ m a x Θ m i n 2
is obtained immediately. Taking into account (3) in (15), we arrive at (17). □

4. Simulations

In this section, we provide results of numerical simulation of the diffusive system above. Following to [3], we choose the system parameters according to Table 1. Involving the numerical scheme for computing the elements of the matrix of Green’s type G i j , i , j = 1 , 2 , developed in [13], we compute u , v according to (13). Then, κ is computed using (14).
Recalling Remark 3, note that in this case the solution of the diffusion model u , v is linear in N.
Numerical evaluation of G i j , i , j = 1 , 2 shows that for the values of parameters represented in Table 1, there always exists a constant C such that (16) holds. Therefore, Theorem 1 is valid and (17) holds. Moreover, it has been observed that for these values, sufficient condition (18) holds implying that the diffusive model above is approximate controllable.
Numerical analysis shows that for any fixed t, κ is a convex function of Θ s p , σ with a unique global minimum. Moreover, for fixed Θ s p 17 , 21 , κ increases when σ increases in 0.5 , 1.5 , and that for fixed σ 0.5 , 1.5 , κ decreases when Θ s p increases in 17 , 21 . Therefore, the global minimum of κ is achieved at Θ s p , σ = 21 , 0.5 . Figure 1 expresses the behavior of the minimal aggregate consumption against the number of TCLs when κ 0 = 1 2 max 0 t T κ t . As it is expected, κ increases with increase of N. It is also seen from Figure 1 that as time increases, the minimal aggregate consumption stabilizes around corresponding κ 0 .
Now let us examine the corresponding solution u , v . In the contour plots below (Figure 2, Figure 3 and Figure 4), dark blue regions correspond to near-zero values and the lighter the color of the domain, the higher the corresponding value. For N = 800 , as Figure 2 shows, for t < 1 , in the whole range of Θ Θ m i n , Θ m a x , u is in the dark blue region meaning that the density of TCLs in the on state is low. A similar region occurs near t = 2 . In between these regions, u approaches its maximal value in the range of Θ 20.2 , Θ m a x . Then, the values of u oscillate between lower and higher values until it is stabilized near u ˜ = 400 for t 9 . A similar behavior is observed also for v with the difference that v is stabilized near v ˜ = 500 TCLs/ C much earlier, for t 5 .
A similar behavior is observed also in the case of N = 1000 , but unlike the previous case, in this case, the stabilization of both u and v occurs much later. However, when N = 1500 , v is stabilized surprisingly fast, for t 6.5 .

5. Conclusions

A diffusive behavior in modelling of smart grids of TCLs has been observed recently. The diffusion naturally affects the aggregate power consumption which is a very important criterion in estimating the efficiency of any smart grid of TCLs. In this paper, we consider the controllability property of the diffusive model of TCLs described by a one-dimensional diffusion equation with variable coefficients. The analysis of the governing system is sophisticated by the fact that the system is coupled through the boundary conditions. Involving the Melnikov method for numerical computation of the corresponding Green’s matrix, we derive a sufficient condition for approximate controllability of the grid providing minimal aggregate power consumption. The set-point value and the bandwidth of the temperature act as control parameters.
Numerical analysis shows that in the considered range of the parameter values, the population model of TCLs with diffusion is approximately controllable. An extensive numerical analysis reveals that the power consumption is a convex function of control parameters having a unique global minimum. The global minimum is found at the point with maximal set-point and minimal bandwidth values. Corresponding densities of TCLs in on and off states are examined numerically.

Author Contributions

Both authors contributed equally to this manuscript.

Funding

This research received no external funding.

Acknowledgments

M.M.H. would like to thank School of Mechatronic Engineering and Automation, Shanghai University, Shanghai 200444, China. The work of A.Z.K. is supported by State Administration of Foreign Experts Affairs of China.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Moura, S.; Ruiz, V.; Bendsten, J. Modeling heterogeneous populations of thermostatically controlled loads using diffusion-advection PDEs. In Proceedings of the ASME Dynamic Systems and Control Conference, ASME, Palo Alto, CA, USA, 21–23 October 2013. [Google Scholar] [CrossRef]
  2. Sinitsyn, N.A.; Kundu, S.; Backhaus, S. Safe protocols for generating power pulses with heterogeneous populations of thermostatically controlled loads. Energy Convers. Manag. 2013, 67, 297–308. [Google Scholar] [CrossRef] [Green Version]
  3. Moura, S.; Bendtsen, J.; Ruiz, V. Parameter identification of aggregated thermostatically controlled loads for smart grids using PDE techniques. Int. J. Control 2014, 87, 1373–1386. [Google Scholar] [CrossRef]
  4. Ghaffari, A.; Moura, S.; Krstić, M. Modeling, Control, and Stability Analysis of Heterogeneous Thermostatically Controlled Load Populations Using Partial Differential Equations. J. Dyn. Syst. Meas. Control 2015, 137, 101009. [Google Scholar] [CrossRef]
  5. Kazmi, H.; Suykens, J.; Balint, A.; Driesen, J. Multi-agent reinforcement learning for modeling and control of thermostatically controlled loads. Appl. Energy 2019, 238, 1022–1035. [Google Scholar] [CrossRef]
  6. Ding, Y.; Cui, W.; Zhang, S.; Hui, H.; Qiu, Y.; Song, Y. Multi-state operating reserve model of aggregate thermostatically-controlled-loads for power system short-term reliability evaluation. Appl. Energy 2019, 241, 46–58. [Google Scholar] [CrossRef] [Green Version]
  7. Eissa, M.M.; Ali, A.A.; Abdel-Latif, K.M.; Al-Kady, A.F. A frequency control technique based on decision tree concept by managing thermostatically controllable loads at smart grids. Int. J. Electr. Power Energy Syst. 2019, 108, 40–51. [Google Scholar] [CrossRef]
  8. Bao, Y.-Q.; Chen, P.-P.; Zhu, X.-M.; Hu, M.-Q. The extended 2-dimensional state-queuing model for the thermostatically controlled loads. Int. J. Electr. Power Energy Syst. 2019, 105, 323–329. [Google Scholar] [CrossRef]
  9. 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]
  10. Malhame, R.; Chong, C. 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]
  11. 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] [Green Version]
  12. 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]
  13. Melnikov, Y.A.; Melnikov, M.Y. Green’s Function. Construction and Application; Walter de Gruyter: Berlin, Germany, 2012. [Google Scholar]
  14. Avetisyan, A.S.; Khurshudyan, A.Z. Controllability of Dynamic Systems: The Green’s Function Approach; Cambridge Scholars Publishers: Cambridge, UK, 2018. [Google Scholar]
Figure 1. The minimal aggregate power consumption against time for N = 800 , 1000 and 1500 TCLs.
Figure 1. The minimal aggregate power consumption against time for N = 800 , 1000 and 1500 TCLs.
Mathematics 07 00523 g001
Figure 2. Contour plots of the distributions u and v for N = 800 .
Figure 2. Contour plots of the distributions u and v for N = 800 .
Mathematics 07 00523 g002
Figure 3. Contour plots of the distributions u and v for N = 1000 .
Figure 3. Contour plots of the distributions u and v for N = 1000 .
Mathematics 07 00523 g003
Figure 4. Contour plots of the distributions u and v for N = 1500 .
Figure 4. Contour plots of the distributions u and v for N = 1500 .
Mathematics 07 00523 g004
Table 1. Values of inputs used during simulation.
Table 1. Values of inputs used during simulation.
QuantityValueDimension
α 0.051/h
P14kW
R2 C/kW
k2.5
Θ 32 C
β 0.01( C) 3 /s
N{800, 1000, 1500}TCLs
u 0 Θ 0TCLs/ C
v 0 Θ 30 N 16 σ 1 4 σ 2 Θ Θ s p 2 2 TCLs/ C
σ 0.5 , 1.5 C
Θ s p 17 , 21 C

Share and Cite

MDPI and ACS Style

Hossain, M.M.; Zh. Khurshudyan, A. Controlling Power Consumption in a Heterogeneous Population Model of TCLs with Diffusion: The Green’s Function Approach. Mathematics 2019, 7, 523. https://doi.org/10.3390/math7060523

AMA Style

Hossain MM, Zh. Khurshudyan A. Controlling Power Consumption in a Heterogeneous Population Model of TCLs with Diffusion: The Green’s Function Approach. Mathematics. 2019; 7(6):523. https://doi.org/10.3390/math7060523

Chicago/Turabian Style

Hossain, Md Musabbir, and Asatur Zh. Khurshudyan. 2019. "Controlling Power Consumption in a Heterogeneous Population Model of TCLs with Diffusion: The Green’s Function Approach" Mathematics 7, no. 6: 523. https://doi.org/10.3390/math7060523

APA Style

Hossain, M. M., & Zh. Khurshudyan, A. (2019). Controlling Power Consumption in a Heterogeneous Population Model of TCLs with Diffusion: The Green’s Function Approach. Mathematics, 7(6), 523. https://doi.org/10.3390/math7060523

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