Autonomous Controller for Flexible Operation of Heat Pumps in Low-Voltage Distribution Network

: This paper aims to unleash the potential of a heat pump (HP) and its storage system, as a ﬂexible consumer load, in the low-voltage (LV) distribution network by introducing an autonomous controller. Steady-state analysis using DigSILENT Power Factory, a power system analysis tool, is performed to verify the proposed hypothesis. The proposed controller manages to operate the individual HP and storage within the recommended operating limits of grid voltage, by sharing ﬂexibility within the speciﬁc network. It has the capability of sensing local key control parameters for scheduling, re-scheduling, and decision-making on the operation of individual HPs. It also takes the thermal energy comfort of individual consumers into consideration. Measurement of local parameters such as grid voltage, supply temperature and level of cold water in the storage tank deﬁnes the priority for operation of HPs based on operating delays for turning it on and off. This enhances the sharing of ﬂexibility for proper coordination, control, and management of HP systems in LV distribution networks with mutual technical beneﬁts. From the results, the application of the proposed controller is found to be effective to manage grid congestions and local voltage regulation, satisfying the thermal energy requirements of the customer.


Introduction
Hot water is used for space heating and domestic applications. With the concept of power-to-heat (P2H) in smart energy systems to incorporate 100% renewables, heat pumps (HPs) are expected to play an important role for thermal heating and cooling due to their high output-to-input ratio. HPs are also considered to be a key technology supporting the European strategy for thermal systems decarbonization [1]. Furthermore, low-temperature district heating, to reduce transmission loss of thermal energy, requires boosting of hot water for end users [2,3]. The combination of distributed HPs and district heating systems, with high penetration of renewables in power system networks, have a significant potential in low-temperature district heating systems [4]. In [5] the authors investigate the life-cycle cost for use of district heating or distributed HPs and find it depends on the energy system economy. Use of electric boilers (EBs) has been unfavorable due to high electricity prices and taxes in the Swedish system whereas, HPs are providing high flexibility in power-to-heat supply [6] and a high coefficient of performance (COP) when compared to EBs. Individual HPs seem to be the best alternative to district heating to replace individual boilers (based on oil, natural gas or biomass) in rural or suburban areas, where heat densities are lower [7]. This indicates the installation of many HPs in the near future as well as a significant increment in residential electricity demand. Control of HPs plays an immense role in the flexible operation for demand response. HPs can be controlled based on a simple start and stop operation. A wide variety of HPs is available on the market environment to meet its design objectives. The main advantage of this proposed control system is that there is no need for extra investment in communication infrastructure for the management of flexible loads. This article focuses on an on/off control strategy (to simplify the problem) to effectively introduce high shares of individual HPs, in low-voltage (LV) residential networks, as a sustainable solution towards P2H systems. The effect on the control strategy based on the selection of different initial conditions and operating power factor of HPs are well illustrated through various case studies.
This paper is organized as follows. A short overview of the working principle of HP and its modeling with respect to power-to-heat (P2H) is presented in Section 2 along with the concept of a model of a simple stratified hot water storage tank. Section 3 presents an overview of the LV residential network and the grid voltage situation with only residential load (without HP connected). The details of the logical concept of the proposed local controller and its significance for the flexible operation of HPs is discussed in Section 4. Simulation environment with various case studies, based on initial conditions, are discussed in Section 5, followed by result and discussion in Section 6. Finally, the outcome of the research is concluded in Section 7.

Heat Pump and Storage
A HP basically consists of three major components: evaporator, compressor, and condenser. Liquid refrigerant is exposed to the source, through the HP's evaporator coil, at a low temperature and pressure. As this refrigerant picks up energy from the source, it changes its state from liquid to gas within the evaporator coil. Then the gaseous fluid is passed through the compressor pump where it is pressurized. As a consequence, the temperature of this pressurized refrigerant increases. This heated refrigerant is passed through the condenser (a coil wrapped around the tank, or a pumped heat exchanger to feed a separate water storage tank) to discharge its heat energy. After passing through the condenser, the heated refrigerant cools and condenses back into liquid form. It is then passed through an expansion valve, where its pressure is reduced, and the cycle starts over. Thus, HP does not produce hot water instantaneously when turned on. It takes around 10-15 min to achieve a steady state condition [23]. Figure 1 shows the step response of the HP during turn-on and turn-off. The HP is turned on at time T 1 , but it only starts delivering hot water after time T 2 . HPs can operate down to the voltage level of 0.7 pu or below. Different types of HPs are available in the market. They can be generally divided into two categories: inverter-based and non-inverter-based units. Non-inverter-based directly connected HPs exhibits low harmonic content and high inrush current compared to inverter-based [24]. The power factor of most inverter-based HPs, operating within the voltage range of 0.9-1 pu, is between 1 to 0.95 leading, whereas few HPs operates within 1 to 0.95 lagging [25]. Non-inverter-based HP draw inrush current nearly 5 times its nominal value, whereas, inverter driven HP units draw inrush current of the same order of magnitude as per their nominal rating. However, in some inverter drive models of HP, there exist spikes due to the charging of inverter's Direct Current (DC) bus capacitors [25].
The thermal power delivered by the HP is given by Equation (1). Here, P HP is the electrical power rating of the HP (W),Q Heat is the heat delivered by HP (W). The step response of HP is modeled as shown in Figure 1 with a time delay of 10 min. The reactive power delivered or consumed by the HP based on leading or lagging pf respectively is given by Equation (2). The relevant parameters of HP such as flow temperature of hot water T hot = 70 • C and COP = 3 during cold winter season (constant value of COP is selected for simplicity) are selected based on the manufacturer's data-sheet.
The simplified hot water storage tank, based on a single thermocline layer of varying height as shown in Figure 2, separating the storage tank into two layers (upper hot and bottom cold) is modeled as presented in [11,26] with appropriate assumptions and theoretical validation considering charging and discharging of storages tank. It is assumed that

•
The heating element (condenser coil of HP) is attached to the bottom cold layer of the tank.

•
There is an exchange of heat between the condenser, at the bottom of the storage tank, and the water inside the storage during operation of the HP.

•
The thermocline inside the tank is horizontal and independent of radial distance from the side wall of the storage tank.

•
The thermocline remains only when the HP is turned off and there is no exchange of heat between the condenser and the storage tank • When HP is turned on and there is heat transfer to the storage tank, the thermal stratification is destroyed by water turbulence and there is mixing of hot and cold liquid creating uniform temperature inside the storage tank [27].
The stratified layer provides crucial information on the energy status of the storage tank based on the level of cold water filled inside the tank and temperature of the hot water inside the tank. These parameters are vital for the control logic to flexibly operate HPs as distributed local units and this is discussed in Section 4. The temperature of hot water (T hot ) is calculated as per Equation (3) during the presence of the thermocline layer.
Here, T a = 15 • C is considered, and the surface area of storage tank is calculated considering ratio of height to diameter as 2.5.
When the HP is turned on and there is an exchange of heat between the condenser and the storage tank, the thermocline layer gets disturbed by water turbulence and mixing of water inside the storage tank takes place, which gives an average temperature (T avg ) of water inside the storage tank as per Equation (4) [11]. During this period T hot = T avg .
Here, Q demand = heat demand transfer rate due to flow of water (Thermal Demand) [J/s] Q heat = heat flow rate of heating element [J/s]

LV Residential Grid Network
To maintain grid voltage within its operating limit (±10%) in the LV distribution network, when integrated with HPs, good knowledge of voltage profile is necessary. It not only helps to define the control parameters for HP and storage but also provide the environment to evaluate the performance of the proposed controller. A 10/0.4 kV LV transformer substation with 6 feeders and 164 individual households located at Northern Jutland in Denmark is considered and is shown in Figure 3. The detailed grid network has been discussed in [29]. It is a real LV residential grid in North Jutland in Denmark, all the relevant electric grid data were supplied by the distribution grid company that owns the network. Thermal and electrical load profile parameters are discussed in detail in [12] and their daily demand are illustrated in Figure 4a,b, respectively. HPs are connected to each individual house. Steady-state power system analysis is performed using the DigSILENT Power Factory simulation tool to investigate the performance of the proposed controller. Table 1 summarizes the LV grid network.   Figure 5a,b shows the minimum terminal voltages of all feeders with only residential loads (i.e., without any HPs connected). It is clearly seen that feeder F3 and F4 is the strongest feeder of all where minimum voltage attained by terminal T19 and T99 is 0.98 pu, followed by F2-T47 (Feeder 2, terminal 47). Terminals T19, T125, and T162 attain the minimum voltage of 0.95 pu during peak hour in feeders F1, F5, and F6, respectively. F6 is heavily loaded with 421 kWh/day and 39 houses.

Control of HP
The complete system structure is illustrated in Figure 6. The proposed strategy of individual control unit (as discussed later in this section with Figure 7) integrated with HP allows proper coordination of HPs in the respective feeder to operate within the grid voltage limit, satisfying thermal demand, without the need for centralized communication for the operational schedule of HPs. This is possible with the local information such as status information from the HP-storage (T hot -Temperature of hot water, X c -the level of cold water in the storage tank), and information from the grid (V t -terminal voltage at the point of coupling of HP). This information, along with on/off delay parameters for operation of HP based on the terminal voltage and thermal storage status, is the key control strategy. Hereby, HPs with storages have the potential to contribute as a flexible controllable load to support grid performance (grid voltage, line loading, and transformer capacities), by shifting its operating time. The proposed controller technique is also able to shift the flexibility request among the consumers in the same feeder based on their storage and grid status.
T return Figure 6. Block diagram of complete system structure.
OR H1 To determine the control limit to maintain minimum and maximum working temperature of hot water inside the thermal storage tank, it is necessary to know the purpose of thermal demand as well as working limits of commercially available products for the heat exchanger to fulfill user demand. Domestic use of hot water generally includes hot water for the kitchen and bathroom sinks, tubs, and other appliances. It requires an instantaneous water heater with a heat exchanger to ensure that hot water is produced immediately when tapping starts. There are products commercially available that operates with low district heating supply temperature in the range 50-70 • C with return temperature between 13-18 • C. On the other hand, space heating uses hot water with various techniques such as a radiator or floor heating with a thermostat to control room temperature between 8-28 • C. Thus, in this work the return temperature (T return ) of hot water to the system is selected to be 30 • C considering that the thermal demand used for both space heating and domestic uses and is supplied by the same source.
The overall control structure of the HP consists of five different control strategies to ensure proper flexibility in the operation of the HP maintaining grid voltage above its lower operating limit of −10% as well as fulfilling consumer demand. The five control strategies are performed using a level controller, temperature controller, voltage controller, critical controller, and a delay controller.
The voltage controller supports to maintains the terminal voltage (V t ), at the point of coupling of the HP to the electrical network, within the lower operating voltage of −10%. HPs in off state at the terminals where V t < 0.94 pu are not allowed to operate until V t ≥ 0.97 pu to avoid hunting effect and is realized using RS-flipflop logic (without considering its inherent time delays) as shown in Figure 7a. However, HPs which are in on state continues to operate until V t ≥ 0.92 pu. This set-up allows the smooth operation of the HPs which are in on state and prevents the operation of HPs that are in off state, during the period when 0.92 ≤ V t < 0.94.
The level controller operates to maintains the level of cold water inside the thermal storage tank (X c ) below a minimum level (X min = 25%) so that the temperature of hot water (T hot ) in the storage tank is maintained between 60-70 • C in the storage tank. The level controller sends operating signal (H), to operate the HP, when X c ≥ X min until T hot ≥ T max (T max -maximum allowable temperature in storage tank). This hysteresis in level control is realized using the RS-flipflop logic as shown in Figure 7b.
The temperature controller operates to maintains the temperature of hot water (T hot ) between 60-70 • C in the storage tank by sending an operating signal (T), to operate the HP, when the minimum temperature (T min ) inside the storage tank is below 60 • C. This situation arrives when there is no thermal demand for a long time and due to heat loss from the storage tank, the temperature of hot water drops over time. The temperature control is realized using RS-flipflop logic, as seen in Figure 7c, to maintain hysteresis in temperature control i.e., temperature control operates when T hot < T min until T hot ≥ T max .
The operation of the HP is mainly dominated by the voltage controller until a critical situation arrives, which is defined by the condition T hot < 50 • C or X c > 40%. The critical controller operates during these critical conditions of the thermal storage tank to maintain T hot between 50-70 • C despite terminal voltage situation until T hot ≥ 70 • C. This ensures a continuous supply of thermal demand for end-user applications such as space heating or domestic use. Figure 7d shows the schematic of a critical controller realized using RS-flipflop logic.
The delay controller determines the delay for the operation of HPs based on V t . The delay logic based on V t plays a vital role to allow operation of HPs in the weak feeder to meet the thermal demand. Delays with Lower pick-up and higher drop-off time corresponding to lower terminal voltage (V t ) prioritize the operation of HPs at the lower end of the distribution feeder, where the V t is lower than the ones closer to the feeder terminal. The controller selects its delay value corresponding to V t based on information provided in Table 2. If V t < 0.94 pu, HPs are not allowed to connect to the grid until storage is in critical conditions as discussed above. Thus, the critical controller has the highest priority followed by the voltage controller. The temperature and level controller has the least priority in control of HP operation as seen from Figure 7d.   Figure 8.
In case the V t is lower than 0.92 pu for 60 s (two samples to avoid transients), the drop-off time is 0 s and system reconnects only after 30 min (1800 s) of time delay. This is to ensure that the grid does not get stressed frequently due to connection and disconnection of HPs. During the event of critical condition, as defined in Figure 7d, HP connects instantaneously without any delay, irrespective of V t , until the temperature of hot water in the storage tank reaches its maximum value T max . Table 3 shows the value of control parameters associated with the control variables shown in Figure 7 for decision-making to operate HPs. Figure 9 shows the validation of the proposed logic for operation of HPs. Around 07:00 hrs, X c > X min and the logic signal H is set to 1 (Figure 9d). HP delivers heat output 15 min after the HP is turned on (Figure 9f). Around 13:00 hrs, HP again tries to connect as X c > X min (Figure 9d). During this time when HP connect to the grid, the grid voltage goes below 0.92 pu for more than 60 s (Figure 9a). As a consequence, voltage control turns off the HP. Once the HP is disconnected V t > 0.97 pu but HP gain turns on only after 30 min (Figure 9f). Again the V t < 0.92 pu and HP turns off. This time V t does not recover its recovery limit of 0.97 pu (Figure 9c) as a result status of the voltage controller is low, meaning unfavorable condition to connect HP at this terminal. Around 19:00 hrs X c is greater than the critical level of cold water in the storage tank (X crit ), i.e., X c > X crit (Figure 9d) and HP connects instantaneously without any delay and delivers thermal power until storage is fully charged (Figure 9f).  The decentralized control system based on the control logic with local information presented in this section supports the management of HPs in LV-weak grids and is further discussed in the result and discussion section.

Simulation Environment
Details on the residential electrical and thermal load demand profile for peak winter is well analyzed in [12]. The houses with higher electricity demand are allocated with higher thermal demand. Table 4 shows the allocation of HP and storage size based on thermal demand to ensure the continuous supply of hot water for at least 8 h during normal operation when fully charged until X c > 25%. Figure 10 shows the rating of HPs allocated to the individual houses in the proposed LV network.  Performance of the proposed decentralized controller discussed in Section 4 for integrating HPs in the LV electrical grid is investigated with various case studies. In this work, steady-state power system analysis in the time domain is performed using DigSILENT Power Factory simulation tool. It has the flexibility for developing models and control algorithms of flexible loads and energy storages such as HPs, and further test distribution grid integration and impact studies [30,31]. The model of the heat pump and storage is discussed in Section 2 and the LV residential network is discussed in Section 3.
Various case studies are performed based on different initial conditions for the level of cold water in the storage tank (X c ), the temperature of hot water from the storage (T Ctrl ), status of the heating element( C a ), and the operating power factor of the HP as shown in Table 5. In case I, all the storage tanks are fully charged. In case II, all the storage tanks are fully discharged. This situation closely corresponds to the scenario of a power failure due to various causes, such as natural calamities or fault in the power system, for a longer time period. The random situation is selected for case III. The random numbers are generated from MATLAB. Inverter base HPs with inrush current equal to nominal rating, and coefficient of performance (COP) of 3 for winter condition are considered. Since inverter-based HPs operates either at leading or lagging power factor as discussed in Section 2, it is interesting to observe their synergy with the electricity network for both conditions. Therefore, case IV is set up with lagging power factor with initial conditions as in the case I.

Results and Analysis
Figures 11 and 12 shows the detailed evaluation and comparison of controller performance with respect to transformer loading, line loading, maximum voltage drop, and performance of the hot water storage tank for the case with leading power factor of HP (case I, II, and III) as mentioned in Table 5. Sub-figures in column 1, 2, and 3 are for case I, II, and III, respectively. The summary of results for day 1, 2, and 3 are tabulated in Tables 6 and 7, whereas figures for day 2 are not presented here. Minimum terminal voltages at the different feeders for all case studies are summarized in Table 7.    Figure 11a shows the transformer loading, and total percentage of HPs operated during the 24 h period starting from midnight. Figure 11b shows the maximum cable loading attained. Figure 11c illustrates minimum terminal voltages attained in feeder F1, F2, F3 and Figure 11d for F4, F5, F6. The corresponding terminal number with minimum voltage in respective feeders are presented in legend. In order to verify that all the thermal demand are fulfilled without any interruption, the temperature profile of hot water and the corresponding level of cold water in a storage tank, which attains maximum temperature drop in each feeder is shown in Figure 11e,f respectively. The associated terminal with the maximum temperature drop in each feeder is illustrated as a legend. The feeder with a maximum voltage drop is farther from the feeder terminal (Figure 11c,d) whereas, the maximum temperature drop is closer to the feeder terminal (Figure 11e). This suggests that the thermal unit at the far end of the feeder has higher priority to operate, based on the delay setting to turn on and off (as detailed in Table 2) compared to ones closer to the feeder terminal. Similar observations and graphs are presented for Day 3 in Figure 12.

Level of Cold Water in Storage Tank
Day 1 analysis as seen in Table 6 clearly shows that for case II (where the storage tank is fully discharged), all HPs operates together, and the transformer is overloaded 22% more than its rated capacity. Also, the maximum line loading is 135% for cable C45 Feeder F6 and a few others are heavily loaded (Figure 11(b2)). However, the line loading for all other feeders except F6 can sustain 100% penetration of HPs in operating condition. On the other hand, cable C45 in feeder F6 along with a few other cables are also overloaded for case III (Figure 11(b3)), where the transformer is loaded up to 95% of its rated capacity (Figure 11(a3)). The maximum level of cold water in a storage tank and the minimum temperature of hot water for case II and III represents the initial condition. However, for case I, the temperature can drop down to 54 • C for thermal unit connected at terminal T162 of feeder F6 (Figure 11(e1)). The minimum temperature is attained due to the mixing of hot and cold water in the storage tank due to heat exchange during operation of HP as discussed in Section 2. Table 7 indicates that the controller limits the terminal voltages to go below 0.92 pu for case I. However, for case II, there are many HPs storage unit with the level of cold water or temperature of hot water below the critical level as an initial condition. During such a condition the controllers are not effective because of their LV ride through property during the critical stage of the storage tank. However, in the long run, this problem is solved by the controllers which manage to flexibly operate the HPs and storage as seen in day 2 and day 3 analysis.

Day 2 and Day 3 Analysis
Similar observations as in Figure 11 are shown in Figure 12 for day 3. Day 2 and Day 3 analysis shows that the maximum transformer loading (60%) and line loading (69%) are well within the limit. In the long run, there are no more than 31% of HPs being operated at any time. Furthermore, the minimum temperature inside the tank is 54 • C at the time of operation of HP and the maximum level of cold water in the storage tank is 32% suggesting that the HP and storage are able to operate flexibly fulfilling the thermal demand simultaneously.
In the long run, despite any initial condition, the HPs coordinates to support peak shaping and line loading (Figure 12a,b), while maintaining the grid voltage (Figure12c,d) as well as thermal demand (Table 6-T hot ) together. Feeder F3 is the strongest feeder among all with voltage drop down to 0.96 pu, and the weakest one is F6. It is observed from Tables 6 and 7, and Figure 12 that the simulation results for all case studies have the similar outcome with respect to line loading, grid voltages, percentage operation of HP, and transformer loading. The grid parameters are well within their operational limit and satisfy the thermal demand simultaneously.

HPs with Lagging Power Factor
Some inverter-based HPs also operates at lagging 0.95 power factor. Thus, it is interesting to observe how the LV grid behaves for HPs with lagging pf. When the pf of all HPs is set to 0.95 lagging, with the initial condition as in the case I, the result is presented in Figure 13. Column I, II, and III of Figure 13 are for the day I, II, and III. Transformer loading and the total percentage of HP in operation is illustrated in Figure 13a. Compared to the case I (where pf is leading), maximum transformer loading and maximum percentage of HPs in operation are reduced (see Table 6), signifying that with lagging pf, the hosting capacity of LV distribution grid decreases. Figure 13b shows that line loading decreases due to reduction in hosting capacity as compared to case I. Terminal voltages of the terminal attaining minimum voltage in each feeder are shown in Figure 13c,d along with the temperature of hot water and level of cold water in storage tank associated with these terminals are shown in Figure 13e,f respectively. It is interesting to see that the terminal voltage of feeder F6, terminal T162 (Figure 13d), goes below 0.92 pu when the HPs associated with this terminal tries to operate. Thus, the HP in T162 is disconnected to support grid voltage. It again tries to connect after the time delay of 30 min and disconnects due to poor voltage, until the level of cold water in the storage tank attains critical limit (greater than 40%) as seen from Figure 13f. Figure 13e shows that the temperature of hot water in the storage tanks in different feeders at the terminal attaining minimum voltage. It shows that the supply temperature of hot water is well above 50 • C. In the long run, as seen from Figure 13 for day II and day III, HP coordinates to support peak shaping and maintain grid voltage simultaneously.

Discussion
The effectiveness of the proposed controller is well explored in this section for strong as well as in weak feeders. The effect of leading and lagging power factor of HPs in grid voltage fluctuation is well observed with the effectiveness of controller performance to set them within the operating limit of terminal voltage between 1-0.9 pu, together with the fulfillment of thermal demand. This enhances application for integrating renewables and synergy with P2H, supporting operation and planning of active grids. More voltage drop is observed with lagging power factor, compared to leading. By controlling its power factor, inverter-based HPs have potential to support grid voltage with many PV or wind turbines. The set points for the controller (control variables) can be specified for a particular network based on its capacity. If required, the set point of controllers can be changed individually (when there are some changes in the grid infrastructure) and is not needed frequently.
Flexibility is shared among the different HPs in the same feeder based on measurement and controlled variables. The controller fails to operate only when many storage tanks are below the critical level as seen in case II and case III initial conditions. In the long run, it manages to stabilize the system flexibility and satisfying demand simultaneously without any need for communication infrastructure, as a cost-effective solution for flexible load management of heating systems, minimizing grid reinforcements.
Along with the necessity to maintain thermal comfort of users at households, the reliable and stable operation of LV distribution grid are equally important when integrating sizeable loads such as HPs. With an increase in renewable energy generation, these HPs and storages are considered to be one of the major technologies that can provide flexibility to the power system providing efficient heating as well as storage of energy. Therefore, for the successful integration and use of HPs in distribution grids, it is important to understand the overall scenario of how the system is affected. Analysis of the distribution grid from the perspective of HP control and its effect on the grid voltage, line loading and transformer loading along with user comfort (as seen from level and temperature of hot water in the storage tank) presents the effectiveness of the proposed controller. The effect of selecting initial condition regarding the status of HP (on/off) and storage (level and temperature of hot water) in the LV grid are thoroughly analyzed based on its effect on grid voltage, and transformer and cable loadings. Same initial condition as in study case I and II introduces operating bottle necks in the LV grid on day one, whereas in long run despite any initial condition, the total percentage of HP in operation at a time is between 25-30% as seen for the analysis of day 3 analysis in Table 6. This reflects the distributed operation of HP over time that enables flexible peak shaping of the electricity demand as well as preventing under-voltage (day 3- Table 7) in the far end of the radial feeder. The peak loading of the transformer is just below 50% (day 3- Table 6) indicating that with proper demand control, more sizeable electrical loads, such as electric vehicles, can also be integrated in the future by adopting such simple and effective local controls. However, considering grid voltage of feeder F6 (Figures 11, 12 and 13c,d) there is a necessity of grid reinforcement in feeder F6 to ensure terminal voltage at the end nodes are within operating thresholds.
The results show that the proposed control architecture of HPs can reduce the need for costly grid reinforcement as well as a separate communication infrastructure for handling of big data, and complex control architecture. Hence, it is beneficial to both the networks/systems and the involved stakeholders. It not only reduces the network impact but also provides flexibility in production, conversion, storage, and end-use through the efficient coordination of the HP operation, helping to lower the electricity consumption at peak times. Synergies among the thermal and electrical systems with the interaction between local and distributed generation units such as wind solar have potential in the future. Thus, the results and findings from the work are in line with the contemporary integrated sustainable energy system concept. This work can be further carried out with economical aspects, scheduling and aggregation services based on forecasting and energy price.

Conclusions
HPs with storage is a good source of P2H in smart energy systems. Higher COP compared to EBs, and storage tank as energy buffer increases its potential as a flexible load and storage solution in future smart grid systems. This paper has analyzed the flexibility of individual HPs in every household, particular to the input voltage in the LV distribution network. The use of a simple local controller can satisfy the consumer demand while coordinating among the various HPs in the same electrical feeder, minimizing transformer loading, line loading, and voltage drop in the long run. Different scenarios based on leading and lagging power factor of HPs, and initial operating condition of HP and storages are demonstrated in a LV network with 164 households. The potential cost-effective solution towards P2H synergy in smart energy system has been presented with the use of the local controller for the individual HP.
Author Contributions: R.S. developed the methodology, corresponding models, and wrote the paper. B.B.-J. and J.R.P. supervised the work and reviewed the paper.
Funding: This work is supported by the Innovation fund Denmark through the DiCYPS project.