Assessment of Electric Vehicle Charging Costs in Presence of Distributed Photovoltaic Generation and Variable Electricity Tariffs

In this paper a general model for the estimation of the uncoordinated charging costs of Electric Vehicles (EVs) in the presence of distributed and intermittent generation, and variable electricity tariffs is presented. The proposed method aims at estimating the monthly average cost of uncoordinated charging of a single EV depending on the hour at which the EV is plugged into the EV Supply Equipment (EVSE). The feasibility and relevance of the proposed model is verified by applying the considered cost estimation method to a suitable use case. A single EV charging service offered at a public building equipped with a Photovoltaic (PV) system has been considered as reference case. The proposed model has been applied to the PV production and loads consumption data collected during one year, and the results of the study compared with the Time-Of-Use (TOU) electricity tariff. The application of the proposed model identified noticeable deviations among the computed EV charging costs and the reference TOU profile, with differences up to 40%, depending on the considered month and on the time of charging during the day. It can be concluded that such model could be used to properly detect opportunities of energy savings, and to define dedicated EV price signals that could help to promote the optimal use of distributed energy resources.


Introduction
The increasing pressure towards decarbonized and sustainable energy systems is heavily affecting the transportation sector, which is progressively called for a radical shift from traditional internal combustion engines (ICEs) towards greener solutions.The growing concerns on the anthropogenic greenhouse gas (GHG) emissions related to the use of road vehicles [1], and on their effect on the air pollution in urban centers [2], are in fact calling for a rapid decarbonization of the urban mobility.From this perspective, electric vehicles (EVs) are generally considered one of the most promising solution, thanks to their independence from the primary energy source, and to the total absence of direct GHG and pollutant emissions.Recent studies demonstrated that the large penetration of EVs could in fact help to significantly reduce both direct and indirect GHG emissions and air pollution in urban areas [3,4].EVs, in fact, avoid the direct combustion of fossil fuels, and might also help to mitigate the drawbacks caused by intermittent renewable energy sources (RESs), e.g., by potentially acting as storage solutions [5][6][7].
Energies 2019, 12, 499 2 of 20 Even though relevant technology developments and rapid cost declines, particularly in the EV battery market, are supporting the prospects of an affordable and effective diffusion of EVs in urban areas, several obstacles are still present, particularly concerning the impact of a large penetration of EVs on the management and operation of existing power networks [8][9][10][11].Several research streams have recently addressed this concern, by proposing advanced optimization strategies for the integration of EVs in smart grids [12][13][14][15][16][17].In particular, the application of the Internet of Things (IoT) paradigm and the use of industrial-grade IoT services [18,19] will make available a large amount of information on the use of road vehicles [20], which could be exploited for improving the integration of EVs into smart grids.At the same time, pervasive communication infrastructures, exploiting also the use of power line communication (PLC) systems [21], and of distributed measurement equipment [22], are expected to allow improved observability and enhanced control functions over EV charging networks, including the provisioning of demand response (DR) capabilities.
A particular interest has been paid to the integration and optimization of EV charging processes in the presence of distributed generation (e.g., from photovoltaic and wind power plants), and variable electricity tariffs, by forming the so-called smart charging concept [16].From this perspective, the smart charging concept, which is usually counterposed to the traditional uncoordinated charging mode, encompasses all of the communication technologies and control strategies which allow the coordinated charging of EVs.In the uncoordinated charging mode, in fact, the power supply provided by the grid to the EVSE depends only on the time of arrival of the EV at the charging station and on the requests of the EV battery charger.Conversely, in the smart charging mode the amount of power provided by the grid and the time of activation of the charging process are controlled by a supervisor, and coordinated with the constraints and operational plans set by distribution grid operators.In Reference [23], for instance, an optimal power flow based optimization framework has been proposed for the coordination of the charging cycles of plug-in hybrid EVs (PHEVs) by minimizing feeder losses and charging costs under time-of-use (TOU) electricity tariffs.In Reference [24], a forecast algorithm for the electrical load of households, based on k-means clustering and an artificial neural network, is used for the optimal scheduling of EV charging in residential applications.In Reference [25], a real-time control strategy for EV charging under dynamic price tariffs is proposed, with the objective of minimizing the charging costs.In Reference [26] a system dynamics (SD) model [27] is proposed to develop a real-time charge pricing mechanism of EVs.In Reference [28] an energy management system is presented to minimize EV charging costs while reducing the energy demand from the grid, and by increasing the self-consumption of Photovoltaic (PV) systems.In Reference [29] a generalized Nash equilibrium problem approach for the management of PHEV charging is presented, by allowing each PHEV to respond price tariff signals to minimize the charging costs, while also considering grid facility constraints.
As suggested by the brief literature review presented above, which represents only an excerpt of the large number of manuscripts published on the matter, the study of advanced EV smart charging strategies has been largely investigated.However, most of the existing works in the literature investigated the topic by mainly focusing on the supervisory control perspective, i.e., by proposing optimal EV charging schedules or variable electricity tariffs with the scope of minimizing the operating cost of the distribution grid.However, even though some of the examined studies took into account the estimation of EV charging costs in coordinated smart grid operations, none of them specifically investigated the modelling of uncoordinated EV charging costs from the prosumers' perspective in the presence of distributed RES generation and under variable electricity tariffs.
As a matter of fact, most of current EV supply equipment (EVSE), particularly if installed at end-users' premises (i.e., not for commercial services, such as public charging stations), operates in uncoordinated charging mode.In addition, following the increasing requests of EV charging services at work places or at public buildings, private companies and public administrations have to deal with the proper evaluation of the operating costs related to the requested services, particularly in the presence of distributed RESs, which typically causes both seasonal and intraday variations of prosumers' energy flows.The aim of the study is thus to provide a general model for the estimation of the uncoordinated charging costs of EVs in presence of distributed and intermittent generation, and variable electricity tariffs.The proposed method aims at estimating the monthly average cost of a single EV uncoordinated charging depending on the hour at which the EV is plugged into the EVSE.
The remainder of the paper is organized as follows: In Section 2, the mathematical definition of the cost estimation method is described in detail, by defining all of the system variables and parameters of the model.Afterwards, the feasibility and relevance of the proposed model is verified by applying the presented cost estimation method to a suitable use case, and by comparing the results of the study with the TOU tariff applied to the energy purchased from the grid.In this study, the evaluation of the operating costs of a single EV charging service offered at a public building equipped with a PV system has been considered as reference case, by applying the proposed model to the PV production and loads consumption data collected during one year.The detailed description of the use case is provided in Section 3. In particular, in Section 3.1 the description of the use case is provided by introducing the main characteristics of the system, and by defining the operation parameters and simulation constraints.In Section 3.2 the PV production and loads consumption data used for the simulation are presented and briefly discussed, while the simulation parameters are defined in Section 3.3.The results of the simulation are then discussed in Section 4, and the conclusion of the study are finally presented in Section 5.

Method
The proper evaluation of the charging cost of EVs must take into account all the energy costs related to the considered process, namely: the direct costs due to the increase of energy consumption from the grid, the loss of revenues due to the decrease of the amount of energy sold to the grid, and the costs related to the increase of the demand charge from the utility.
The average cost of an EV uncoordinated charging cycle starting at hour s of month m, C EV s,m , can be computed as follows: where s,m is the term due to the increase of the demand charge from the utility related to the specific EV charging cycle starting at hour s of month m, C ∆E FU s,m is the term related to the increase of energy consumption from the grid, C ∆E TU s,m is the loss of revenues due to the decrease of the amount of energy sold to the utility (which corresponds to the increase of self-consumed energy due to the specific EV charging cycle), and E EV s,m is the expected energy consumption of the EV charging cycle starting at hour s of month m.
The term C ∆P FU s,m is computed as follows: where: is the increase of the demand charge from the utility during month m, due to the EV charging cycle starting at hour s, expressed in kW; • T P FU m is the cost of the demand charge from the utility during month m, expressed as €/kW; • n s,m is the number of charging cycles starting at hour s during month m; In particular, the term ∆P FU,EV s,m is computed as follows: where P FU,EV i,d,m,s and P FU i,d,m are the average value of power absorbed from the utility for each i-th time interval of 15 min during the d-th day of month m, for an EV charging cycle starting at hour s, and for the reference case, respectively.
It must be noted that the term C ∆P FU s,m is computed taking into account the weight of the increase of power from the utility due to the EV charging cycle occurring at hour s, with respect to the sum of the increases of power due to each simulated charging cycle.In this way, the overall cost related to the increase of the demand charge of the system can be distributed to each charging cycle, depending on its specific contribution to that cost.

The term C
∆E FU s,m is computed as follows: where: • where E FU,EV h,d,m,s and E FU h,d,m are the values of energy absorbed from the grid during hour h of the d-th day of month m, for an EV charging cycle starting at hour s, and for the reference case, respectively.

The term C
∆E TU s,m is computed as follows: where: • T E TU h,m is the price of electricity sold to the grid during hour h of month m, expressed in €/kWh; • ∆E TU,EV h,d,m,s is the decrease of the energy injected to the grid (expressed in kWh) during hour h of the d-th day of month m, for an EV charging cycle starting at hour s.In particular, ∆E TU,EV h,d,m,s is computed as follows: where E TU,EV h,d,m,s and E TU h,d,m are the values of energy injected into the grid during hour h of the d-th day of month m, for an EV charging cycle starting at hour s, and for the reference case, respectively.

Use Case
In the present study, the evaluation of the operating costs of a single EV charging service offered at a public building located in the engineering campus of the University of Brescia, Italy, and equipped with a 64 kWp PV system, has been considered as reference use case.A brief description of the system and the definition of the operation parameters and simulation constraints of the simulated EV charging service are described in Section 3.1.The estimation of the EV charging costs for the considered case is carried out by applying the proposed model to the PV production and loads consumption data collected during 2017, as discussed in Section 3.2.In this study, it is assumed that a single point EV charging service is offered to the public during the opening hours of the building.The service calendars and EV charging profile used for the simulation are finally described in Section 3.3.

Description of the Use Case System
A public building located in the engineering campus of the University of Brescia, Italy, has been chosen as reference use case.The building includes both department facilities (i.e., offices and classrooms), and student services, comprising dormitories, study halls, a cafeteria, a gym, and a baseball field.The building is connected to the local low voltage (LV) distribution grid, and is equipped with a 64 kWp PV system.
In this study, it is assumed that a single point EVSE is installed at the LV point of delivery (POD) of the system, and that an EV charging service is offered publicly during the opening hours of the building.In this case, the results of the simulation are expected to provide useful information to the energy and facility manager of the University, by estimating the operating costs of the proposed EV charging service, by taking into consideration the variability and intermittency of the PV power generation, and the energy consumption of the building.Such information could be, for instance, compared to the current electricity tariffs applied by the local energy retailer, and then used to better estimate the price that must be charged to final users, depending on the time of EV charging.

PV Production and Loads Consumption Data
In the following, the PV production and loads consumption data used for the simulation are briefly introduced, by presenting the monthly average and peak power data of the PV production and of user's loads consumption recorded from January to June 2017, and from July to December 2017, as depicted in Figures 1 and 2, respectively.
The presented monthly average and peak data are computed from a set of data recorded every 5 min.Each of the recorded value is computed as the average of data sampled every 10 s.The yellow and red area represent the standard deviation of the measured PV power production and of the measured loads power consumption, respectively.Further information on the measurement system used to collect the herein presented data is provided in Reference [30].
As it can be observed by the graphs of Figures 1 and 2, the power consumption profile of loads shows a typical behavior during almost all the months of the year, by forming a sort of "reverse tilde" shape.The power demand of loads reaches its maximum value among 8 and 9 p.m. (with monthly averages ranging from 26.4 to 50.5 kW), and then rapidly decreases up to 1 am, by afterwards slowly reaching its minimum around 4 and 6 a.m. in the morning (with monthly averages ranging from 15.4 to 26.3 kW).After 6 a.m. in the morning, the power consumption rapidly increases up to 8 a.m., by then remaining almost constant up to 4 p.m. in the afternoon, and then by further increasing to reach the peak around 8 and 9 p.m.It must be also noted that this particular shape is generally replicated each month, with only few differences during July and August (apparently due to summer vacations), where the power consumption profile shows a sort of flattening, i.e., a reduction between peaks and valleys.Moreover, the typical shape of the daily power consumption profile also reveals a sort of "drift" effect during the year, by progressively showing a global down shift from January to August, and then a global up shift from September to December.Finally, it can also be observed that also the computed standard deviation of the mean shows noticeable variations depending on the time of the day and on the considered month.Generally, it can be noted that, for each observed month, the power demand profile seems quite stable, i.e., by showing very small deviations (ranging from 0.5 to 2.6 kW), while during peaks, or during months characterized by unstable weather conditions (such as April), the standard deviation of the average power consumption can be very high, by ranging from 5.1 to 8.9 kW.

PV Production and Loads Consumption Data
In the following, the PV production and loads consumption data used for the simulation are briefly introduced, by presenting the monthly average and peak power data of the PV production and of user's loads consumption recorded from January to June 2017, and from July to December 2017, as depicted in Figures 1 and 2, respectively.The presented monthly average and peak data are computed from a set of data recorded every 5 min.Each of the recorded value is computed as the average of data sampled every 10 s.The yellow and red area represent the standard deviation of the measured PV power production and of the measured loads power consumption, respectively.Further information on the measurement system used to collect the herein presented data is provided in Reference [30].
As it can be observed by the graphs of Figures 1 and 2, the power consumption profile of loads shows a typical behavior during almost all the months of the year, by forming a sort of "reverse tilde" shape.The power demand of loads reaches its maximum value among 8 and 9 p.m. (with monthly averages ranging from 26.4 and 50.5 kW), and then rapidly decreases up to 1 am, by afterwards slowly reaching its minimum around 4 and 6 a.m. in the morning (with monthly averages ranging from 15.4 and 26.3 kW).After 6 a.m. in the morning, the power consumption rapidly increases up to 8 a.m., by then remaining almost constant up to 4 p.m. in the afternoon, and then by further increasing to reach the peak around 8 and 9 p.m.It must be also noted that this particular shape is generally Data depicted in Figures 1 and 2 also show that, as expected, the power production of the PV system exhibits the typical daily and seasonal variations during the year, and that its standard deviation is noticeably greater than that of loads.The comparison between the PV and loads power profiles also reveals that during months characterized by high solar radiation (i.e., from April to September), the PV production may be greater than the loads consumption, thus producing reverse power flows into the grid.
The power balance at the user's POD has been then computed for each 5 min time interval, by taking into account the power consumption of the loads and the PV production, and the monthly average, standard deviation, and peak values have been plotted for each month, from January to June 2017, and from July to December 2017, as depicted in Figures 3 and 4, respectively.also the computed standard deviation of the mean shows noticeable variations depending on the time of the day and on the considered month.Generally, it can be noted that, for each observed month, the power demand profile seems quite stable, i.e., by showing very small deviations (ranging from 0.5 to 2.6 kW), while during peaks, or during months characterized by unstable weather conditions (such as April), the standard deviation of the average power consumption can be very high, by ranging from 5.1 to 8.9 kW.September), the PV production may be greater than the loads consumption, thus producing reverse power flows into the grid.
The power balance at the user's POD has been then computed for each 5 min time interval, by taking into account the power consumption of the loads and the PV production, and the monthly average, standard deviation, and peak values have been plotted for each month, from January to June 2017, and from July to December 2017, as depicted in Figures 3 and 4, respectively.The maximum power consumption from the grid is computed for each time interval of 15 min, and then used to determine the monthly charged power from the utility (i.e., equal to the highest value of all of the computed 15 min maximum values).
As it can be observed by the graphs depicted in Figures 3 and 4, the power profile at the user's POD shows a typical "tilted S" shape, which results from the combination of the typical profiles of the loads consumption and of the PV production.In particular, it can be noted that the computed standard deviations of the power flows at the user's POD are noticeably greater than that of user's loads, with Energies 2019, 12, 499 10 of 20 highest monthly values ranging from 10.4 to 17.4 kW, while the same values recorded by user's loads rage from 5.1 to 8.9 kW.This is apparently due to the variability and uncertain production of the PV system, which causes relevant disturbances to the power flows measured at the user's POD.
The monthly and yearly energy balance of the system is then computed based on all of the data recorded during 2017 every 5 min, as depicted in Figure 5.In sub-figure (a) the monthly energy flows of the system are reported for each month, by specifying the amount of energy absorbed from the grid, injected into the grid, and self-consumed by the system (i.e., produced by the photovoltaic system and consumed by the user's loads).In sub-figures (b), (c), and (d), the interquartile analyses of the daily energy absorbed from the grid, injected into the grid, and self-consumed by the system are reported for each month, respectively.and then used to determine the monthly charged power from the utility (i.e., equal to the highest value of all of the computed 15 min maximum values).
As it can be observed by the graphs depicted in Figures 3 and 4, the power profile at the user's POD shows a typical "tilted S" shape, which results from the combination of the typical profiles of the loads consumption and of the PV production.In particular, it can be noted that the computed standard deviations of the power flows at the user's POD are noticeably greater than that of user's loads, with highest monthly values ranging from 10.4 and 17.4 kW, while the same values recorded by user's loads rage from 5.1 to 8.9 kW.This is apparently due to the variability and uncertain production of the PV system, which causes relevant disturbances to the power flows measured at the user's POD.
The monthly and yearly energy balance of the system is then computed based on all of the data recorded during 2017 every 5 min, as depicted in Figure 5.In sub-figure (a) the monthly energy flows of the system are reported for each month, by specifying the amount of energy absorbed from the grid, injected into the grid, and self-consumed by the system (i.e., produced by the photovoltaic system and consumed by the user's loads).In sub-figures (b), (c), and (d), the interquartile analyses of the daily energy absorbed from the grid, injected into the grid, and self-consumed by the system are reported for each month, respectively.If we look at charts of Figure 5, it can be noted that the amount of energy absorbed from the utility decreases from January to August, and then increases during the following months, with values ranging from 8.6 MWh to 23.5 MWh.Conversely, the amount of energy injected into the utility, due to the overproduction of the PV system, and the amount of self-consumed energy increase from January to August, and then decrease during the following months, with values raging from almost zero to 3.6 MWh, and from 2.5 MWh to 7.7 MWh, respectively.Finally, the energy key performance indicators (KPIs) of the system have been computed, based on data measured during 2017, as shown in Figure 6.In subfigures (a) and (b) the self-sufficiency and the self-consumption of the system are reported for each month, respectively.In sub-figure (c), the Grid-Dependency (GD), Grid-Feed (GF), Self-Sufficiency (SS), and Self-Consumption (SC) indices are reported.

Simulation Parameters
In this study, a typical EV charging power profile is assumed as reference charging cycle.The considered EV charging profile, depicted in Figure 7, was measured during a real charging process of a commercial EV (Renault ZOE), with an initial state of charge (SOC) of 35%, up to a 100% SOC.In sub-figures (a) and (b) the self-sufficiency and the self-consumption of the system are reported for each month, respectively.In sub-figure (c), the grid-dependency (GD), grid-feed (GF), self-sufficiency (SS), and self-consumption (SC) indices are reported for the year 2017.
The self-sufficiency index of the system estimates the independency of the system from the power grid and is computed as the ratio between the energy provided by the PV to loads, and the overall loads consumption.In particular, the self-sufficiency index can be computed as follows: where: • SS m is self-sufficiency index of month m; The self-consumption index estimates the amount of energy produced by the PV plant, which is directly consumed by the system loads, and is computed as the ratio between the energy provided by the PV to loads, and the overall PV production.In particular, the self-consumption index can be computed as follows: where SC m is self-consumption index of month m.
The grid-dependency index estimates the dependency of the user's loads from the grid and is computed as the ratio between the grid consumption and the overall loads consumption.In particular, the grid-dependency index can be computed as follows: where: • GD m is grid-dependency index of month m; • P FU i,d,m is the average power absorbed from the grid by the user's loads during the i-th time interval of day d of month m.
Finally, the grid-feed index estimates the overproduction ratio of the PV system and is computed as the ratio between the energy injected into the grid by the PV system, and the overall PV production.In particular, the grid-feed index can be computed as follows: where: • GF m is grid-feed index of month m; • P TU i,d,m is the average power injected into the grid by the PV system during the i-th time interval of day d of month m.

Simulation Parameters
In this study, a typical EV charging power profile is assumed as reference charging cycle.The considered EV charging profile, depicted in Figure 7, was measured during a real charging process of a commercial EV (Renault ZOE), with an initial state of charge (SOC) of 35%, up to a 100% SOC.

Simulation Parameters
In this study, a typical EV charging power profile is assumed as reference charging cycle.The considered EV charging profile, depicted in Figure 7, was measured during a real charging process of a commercial EV (Renault ZOE), with an initial state of charge (SOC) of 35%, up to a 100% SOC.A three-phase power meter UPM209 produced by Algodue Elettronica s.r.l.(Maggiora, Italy), equipped with a Rogowski coil current transducer, was used to measure the power consumption of

Results and Discussion
The average cost of the EV uncoordinated charging starting at hour s for each month m, C EV s,m , was computed by using the measured loads consumption and PV production data described in Section 3.For each day of month m, and for each hour s, an EV charging process starting at hour s was simulated by adding the EV charging power profile to the measured loads power profile.The simulated POD power profile was then built by taking into account the measured PV power production, and the related monthly average charging cost was computed according to the model described in Section 2. The simulation has been implemented by means of a dedicated C++ code.
In the following, the results of the simulation are presented and discussed in detail.
In Figures 8 and 9, the computed monthly average cost of EV charging (reported as €/kWh) from January to June 2017, and from July to December 2017 are reported, respectively.The computed EV charging cost for each hour (i.e., depending on the time of arrival of the EV at the charging station), represented with the blue line, is compared to the time-of-use (TOE) electricity tariff applied to the amount of energy purchased from the grid (red line).
The results show that the EV charging costs computed by the proposed model may noticeably differ from the reference TOU profile, depending on the considered month and on the time of charging during the day.In particular, it can be noted that during seasons characterized by low PV production (i.e., from October to March), the computed EV charging cost is very close or equal to the TOU tariff applied by the utility to the amount of energy purchased by the grid.Only during early morning and late afternoon hours the computed EV charging cost is slightly higher than the TOU tariff, due to the sensible increase of rated power demand from the utility.
Conversely, during seasons with higher PV production (i.e., from April to September), the computed EV charging cost noticeably differ from the TOU tariff, particularly during late mornings and early afternoons, where the PV production reaches its higher values.It must be also noted that the highest reported difference (occurring at 1 p.m. in August) is on the order of 0.07 €/kWh, which corresponds to a difference of 1.24 € per each EV charging cycle, thus representing a possible estimated saving of about 40% with respect the considered TOU tariff.
To better estimate the effect of each tariff component on the overall computed EV charging cost, in Figures 10 and 11, the EV charging cost components (reported in €) from January to June 2017, and from July to December 2017 are reported, respectively.The cost components are presented for each charging cycle, depending on the time of arrival of the EV at the charging station.
As it can be noted by the analysis of the cost components, during seasons with low PV production (i.e., from October to March), the cost due to the amount of energy purchased by the grid is predominant over the other costs, thus leading to overall costs close to the applied TOU tariff.Conversely, during seasons with higher PV production (i.e., from April to September), the costs due to the loss of revenues related to the amount of energy self-consumed and not sold to the utility become more relevant, such as during late mornings and early afternoons of July and August.This clearly explains why during such periods the computed EV charging cost is noticeably lower than the considered TOU tariff.Finally, it must be noted that the effect of the demand charge cost component is almost negligible because of the low value currently applied for the rated power tariff.Section 3.For each day of month , and for each hour , an EV charging process starting at hour  was simulated by adding the EV charging power profile to the measured loads power profile.The simulated POD power profile was then built by taking into account the measured PV power production, and the related monthly average charging cost was computed according to the model described in Section 2. The simulation has been implemented by means of a dedicated C++ code.
In the following, the results of the simulation are presented and discussed in detail.In Figures 8 and 9, the computed monthly average cost of EV charging (reported as €/kWh) from January to June 2017, and from July to December 2017 are reported, respectively.The computed EV  The results show that the EV charging costs computed by the proposed model may noticeably differ from the reference TOU profile, depending on the considered month and on the time of charging during the day.In particular, it can be noted that during seasons characterized by low PV production (i.e., from October to March), the computed EV charging cost is very close or equal to the TOU tariff applied by the utility to the amount of energy purchased by the grid.Only during early morning and late afternoon hours the computed EV charging cost is slightly higher than the TOU tariff, due to the sensible increase of rated power demand from the utility.corresponds to a difference of 1.24 € per each EV charging cycle, thus representing a possible estimated saving of about 40% with respect the considered TOU tariff.
To better estimate the effect of each tariff component on the overall computed EV charging cost, in Figures 10 and 11, the EV charging cost components (reported in €) from January to June 2017, and from July to December 2017 are reported, respectively.The cost components are presented for each charging cycle, depending on the time of arrival of the EV at the charging station.Conversely, during seasons with higher PV production (i.e., from April to September), the costs due to the loss of revenues related to the amount of energy self-consumed and not sold to the utility become more relevant, such as during late mornings and early afternoons of July and August.This clearly explains why during such periods the computed EV charging cost is noticeably lower than the considered TOU tariff.Finally, it must be noted that the effect of the demand charge cost component is almost negligible because of the low value currently applied for the rated power tariff.

Conclusions
The results of the simulation show that, as expected, the contribution of distributed PV generation to the charging costs of EVs strongly depend on the period of the year and on the time of connection of the EV to the EVSE during the day.At the same time, the results also show that the application of the proposed model is able to better estimate the EV charging costs with respect to the TOU energy tariffs applied by the utility.In particular, during seasons with higher PV production (i.e., from April to September), the computed EV charging cost resulted noticeably lower from the considered TOU tariff, with a highest difference of about 0.07 €/kWh.It must be noted that this value corresponds to a difference of 1.24 € per each EV charging cycle, which leads to a possible saving estimation of about 40% with respect the TOU tariff.

Figure 1 .
Figure 1.Monthly average and peak power data of photovoltaic (PV) production and user's loads recorded from January to June 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.The yellow and red area represent the standard deviation of PV and loads data, respectively.UTC: Universal Time Coordinated.

Figure 1 .Figure 2 .
Figure 1.Monthly average and peak power data of photovoltaic (PV) production and user's loads recorded from January to June 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.The yellow and red area represent the standard deviation of PV and loads data, respectively.UTC: Universal Time Coordinated.

Figure 2 .
Figure 2. Monthly average and peak power data of photovoltaic (PV) production and user's loads recorded from July to December 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.The yellow and red area represent the standard deviation of PV and loads data, respectively.UTC: Universal Time Coordinated.

Figure 3 .
Figure 3. Monthly average and peak power data of the user's point of delivery (POD), recorded from January to June 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.Positive and negative values represent power flows from the grid and to the grid, respectively.Data are sampled every 10 s and the averages are recorded every 5 min.The red area represents the standard deviation of the measured POD power flows.UTC: Universal Time Coordinated.

Figure 3 .
Figure 3. Monthly average and peak power data of the user's point of delivery (POD), recorded from January to June 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.Positive and negative values represent power flows from the grid and to the grid, respectively.Data are sampled every 10 s and the averages are recorded every 5 min.The red area represents the standard deviation of the measured POD power flows.UTC: Universal Time Coordinated.

Figure 4 .
Figure 4. Monthly average and peak power data of the user's point of delivery (POD), recorded from July to December 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.Positive and negative values represent power flows from the grid and to the grid, respectively.Data are sampled every 10 s and the averages are recorded every 5 min.The red area represents the standard deviation of the measured POD power flows.UTC: Universal Time Coordinated.

Figure 4 .
Figure 4. Monthly average and peak power data of the user's point of delivery (POD), recorded from July to December 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.Positive and negative values represent power flows from the grid and to the grid, respectively.Data are sampled every 10 s and the averages are recorded every 5 min.The red area represents the standard deviation of the measured POD power flows.UTC: Universal Time Coordinated.

Figure 5 .
Figure 5. Energy balance of the system based on data measured during 2017.In sub-figure (a) the monthly energy flows of the system are reported for each month, by specifying the amount of energy absorbed from the grid, injected into the grid, and self-consumed by the system (i.e., produced by the photovoltaic system and consumed by the user's loads).In sub-figures (b), (c), and (d), the interquartile analyses of the daily energy absorbed from the grid, injected into the grid, and selfconsumed by the system are reported for each month, respectively.

Figure 5 .
Figure 5. Energy balance of the system based on data measured during 2017.In sub-figure (a) the monthly energy flows of the system are reported for each month, by specifying the amount of energy absorbed from the grid, injected into the grid, and self-consumed by the system (i.e., produced by the photovoltaic system and consumed by the user's loads).In sub-figures (b), (c), and (d), the interquartile analyses of the daily energy absorbed from the grid, injected into the grid, and self-consumed by the system are reported for each month, respectively.

Energies 2018 ,Figure 6 .
Figure 6.Energy key performance indicators (KPIs) of the system, based on data measured during 2017.In sub-figures (a) and (b) the self-sufficiency and the self-consumption of the system are reported for each month, respectively.In sub-figure (c), the grid-dependency (GD), grid-feed (GF), self-sufficiency (SS), and self-consumption (SC) indices are reported for the year 2017.

Figure 7 .
Figure 7. EV charging profile used as reference charging cycle in the simulations.EV: Electric Vehicle.

Figure 6 .
Figure 6.Energy key performance indicators (KPIs) of the system, based on data measured during 2017.In sub-figures (a) and (b) the self-sufficiency and the self-consumption of the system are reported for each month, respectively.In sub-figure (c), the grid-dependency (GD), grid-feed (GF), self-sufficiency (SS), and self-consumption (SC) indices are reported for the year 2017.
,m and P L i,d,m are the average power produced by the PV system and consumed by the user's loads during the i-th time interval of day d of month m, respectively; • ∆t is the duration of the considered time interval.

Figure 6 .
Figure 6.Energy key performance indicators (KPIs) of the system, based on data measured during 2017.In sub-figures (a) and (b) the self-sufficiency and the self-consumption of the system are reported for each month, respectively.In sub-figure (c), the grid-dependency (GD), grid-feed (GF), self-sufficiency (SS), and self-consumption (SC) indices are reported for the year 2017.

Figure 7 .
Figure 7. EV charging profile used as reference charging cycle in the simulations.EV: Electric Vehicle.

Figure 7 .
Figure 7. EV charging profile used as reference charging cycle in the simulations.EV: Electric Vehicle.

Figure 8 .
Figure 8. Monthly average cost of electric vehicle (EV) charging (€/kWh) from January to June 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.The computed EV charging cost for each hour (blue line), is compared to the time-of-use (TOE) electricity tariff applied to the amount of energy purchased from the grid (red line).UTC: Universal Time Coordinated.

Figure 8 .
Figure 8. Monthly average cost of electric vehicle (EV) charging (€/kWh) from January to June 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.The computed EV charging cost for each hour (blue line), is compared to the time-of-use (TOE) electricity tariff applied to the amount of energy purchased from the grid (red line).UTC: Universal Time Coordinated.

Figure 9 .
Figure 9. Monthly average cost of electric vehicle (EV) charging (€/kWh) from July to December 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.The computed EV charging cost for each hour (blue line), is compared to the time-of-use (TOE) electricity tariff applied to the amount of energy purchased from the grid (red line).UTC: Universal Time Coordinated.

Figure 9 .
Figure 9. Monthly average cost of electric vehicle (EV) charging (€/kWh) from July to December 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.The computed EV charging cost for each hour (blue line), is compared to the time-of-use (TOE) electricity tariff applied to the amount of energy purchased from the grid (red line).UTC: Universal Time Coordinated.

Figure 10 .
Figure 10.Electric vehicle (EV) charging cost components (€) from January to June 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.The cost components are presented for each charging cycle, depending on the time of arrival of the EV at the charging station.UTC: Universal Time Coordinated.

Figure 10 .
Figure 10.Electric vehicle (EV) charging cost components (€) from January to June 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.The cost components are presented for each charging cycle, depending on the time of arrival of the EV at the charging station.UTC: Universal Time Coordinated.

Figure 11 .
Figure 11.Electric vehicle (EV) charging cost components (€) from July to December 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.The cost components are presented for each charging cycle, depending on the time of arrival of the EV at the charging station.UTC: Universal Time Coordinated.

Figure 11 .
Figure 11.Electric vehicle (EV) charging cost components (€) from July to December 2017, as depicted from sub-figure (a) to sub-figure (f), respectively.The cost components are presented for each charging cycle, depending on the time of arrival of the EV at the charging station.UTC: Universal Time Coordinated.
m is the cost of electricity absorbed from the grid during hour h of month m (€/kWh);