The Green ’ s type matrices for consumption reduction in a heterogeneous population model of TCLs with diffusion

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.


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 reference [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 a proper choice of 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]: 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 i th load.Except differential equation ( 1), Θ i is subject to a natural constraint Here, Θ min,i and Θ max,i are the minimal and maximal admissible temperatures.Then, the set-point temperature of the i th load, Θ sp,i , is determined in terms of Θ min,i and Θ max,i as follows: On the other hand, introducing the width of the temperature band of the i th load, σ i , we have 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 , Θ min,i , and Θ max,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]: where 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:

Preprints
where k is the performance coefficient.The quantities Θ sp,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].

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 where and β is a diffusivity coefficient with unit (°C) 3 /s.Homogeneous system ( 8) is equipped with the following boundary conditions: The following initial conditions are also given: Remark 2. Note that system (8) is coupled via boundary conditions ( 9)- (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.

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 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 Boundary conditions ( 9), (10)  Following to [14], the general solution of ( 12)-( 11) can be represented as follows: where G ij = G ij (Θ, ϑ, t, τ), i, j = 1, 2, are the elements of the matrix of Green's type.An efficient numerical scheme for computing G ij has been developed by Melnikov in [14].
Remark 3. Note that solution ( 13) is linear with respect to the initial state (u 0 , v 0 ).

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

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 Θ sp and σ accordingly.Such a problem for the simplest model described by (4) has been considered, e.g., in [? ].The problem is to choose appropriate values for Θ sp and σ in order to ensure 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 Θ sp 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 [13], 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 [13].Substituting ( 14) into (15), we make the dependence R = R Θ sp , σ explicit.Remark 6.It should be noted that R depends on Θ sp , σ 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 Θ sp , σ > 0 such that 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 ε, Remark 7. Depending on values of κ 0 , (18) may not be satisfied by any values of Θ sp , σ with a given ε.See also Remark 4. In such cases, inequality (15) must be evaluated.

Simulations
In this section, we provide results of numerical simulation of the diffusive system above.Following to [? ], 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 ij , i, j = 1, 2, developed in [14], 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 TCLs/°C σ (0.5, Numerical evaluation of G ij , 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 Θ sp , σ with a unique global minimum.Moreover, for fixed Θ sp ∈ [17, 21], κ increases when σ increases in [0.5, that for fixed σ ∈ [0.5, 1.5], κ decreases when Θ sp increases in [17, 21].Therefore, the global minimum of κ is achieved at Θ sp , σ = (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  Now let us examine the corresponding solution (u, v).In the contour plots below (Figures 2-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 Θ ∈ [Θ min , Θ max ], 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, Θ max ].Then, the values of u oscillate between lower and higher values until it is stabilized near ũ = 400 for t ≥ 9.A similar behavior is observed also for v with the difference that v is stabilized near ṽ = 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.

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.

Figure 1 .
Figure 1.The minimal aggregate power consumption against time for N = 800, 1000 and 1500 TCLs

Figure 2 .
Figure 2. Contour plots of the distributions u and v for N = 800

Figure 3 .Figure 4 .
Figure 3. Contour plots of the distributions u and v for N = 1000

Table 1 .
Values of inputs used during simulation