Optimal Coordinated Management of a Plug-In Electric Vehicle Charging Station under a Flexible Penalty Contract for Voltage Security

Jip Kim 1, Seung Wan Kim 2,*, Young Gyu Jin 3, Jong-Keun Park 2 and Yong Tae Yoon 2 1 Power System Research Division, Korea Electrical Engineering & Science Research Institute, Bldg. 130, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Korea; jipkim@snu.ac.kr 2 Department of Electrical and Computer Engineering, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Korea; jkpark@snu.ac.kr (J.-K.P.); ytyoon@snu.ac.kr (Y.T.Y.) 3 School of Electrical Engineering, Korea University, 145 Anam-ro, Seongbuk-gu, Seoul 02841, Korea; ygjin@snu.ac.kr * Correspondence: pc9873@snu.ac.kr; Tel.: +82-2-880-9144


Introduction
Plug-in electric vehicles (PEVs) have been attracting interest as an alternative solution to manage the increasing amount of CO 2 emission.Authorities in a number of countries have set ambitious goals for integrating a large number of PEVs [1,2].However, the charging demands for PEVs are relatively larger than the household loads, and are inherently uncoordinated due to the unpredictable behavior of the consumers.Consequently, the high penetration of PEVs may threaten the stable and reliable operation of the distribution network.Specifically, there can be severe voltage drops, fluctuation problems, and degradation in the system efficiency [3][4][5].
To address these concerns, various PEV charging methods have been suggested from the perspective of a distribution network operator (DNO).In [6], the DNO manages the charging demands of PEVs to minimize the line losses of the distribution network.A similar charging method is suggested to minimize the line losses in [7], but it uses stochastic programming to consider the uncertainty of loads in the optimization problem.In [8], a real-time charging strategy for managing the uncoordinated and uncertain charging demands of PEVs is presented to minimize the total operating cost of the distribution network, including the line losses.In order to minimize the operating costs, the authors of [9] propose a method in which the inputs for optimization are continuously updated, and the optimization results are periodically renewed at each time step to reflect the difference between the forecasted value and the actual inputs.In the literature [6][7][8][9], the DNO directly manages the charging schedules of multiple PEVs considering the network constraints to pursue efficient and reliable operation of the distribution network.However, due to operational complexity and burden from the increasing numbers of PEVs, it seems to be difficult for the DNO to directly manage the charging of multiple PEVs in the distribution network [10,11].Moreover, various concepts of intermediate entities that are independent of the DNO have been emerging in practice [12,13].For these reasons, many studies assume that an intermediate entity is delegated to manage PEV charging demands.The intermediate entity does not have any responsibility for the stable and reliable operation of the distribution network, unlike the DNO, but the entity simply aims to maximize its profit with a proper method for managing the PEV charging demand [14][15][16][17].This conflict of interest can be mitigated by various types of coordination measures.The coordination measure should enable the DNO to inform the intermediate entity of the network status or to motivate the intermediate entity to shift its charging demand to another time.However, methods for coordination between the DNO and the intermediate entity have not been thoroughly investigated in the literature.
This study assumes that the PEV charging station is located at a specific bus in the distribution network and that hundreds of PEVs request to be provided with charging service at the same station.Furthermore, it is assumed that the intermediate entity, the so-called charging station operator (CSO), is delegated to manage the charging demands of PEVs.In this circumstance, the CSO only seeks to profit from providing charging service to the PEVs without considering the stability and reliability of the distribution network.Therefore, the hundreds of PEVs in the charging station can significantly drop the voltage of the following bus in the same feeder as well as the bus where the station is located.For this reason, this study suggests a flexible penalty contract (FPC) as a coordination measure, which is designed to induce the CSO not to compromise the voltage conditions of the distribution network.In this scheme, the DNO notices hourly the recommended charging power that does not cause voltages below the lower limit for the CSO, and imposes a penalty price for exceeding the maximum power as an addition to the current electricity tariff, such as the time-of-use (TOU) rate.Furthermore, a formulation for the optimal coordinated charging management of the CSO under the FPC is proposed in this paper.
The remainder of this paper is structured as follows: Section 2 presents the needs for coordination between the CSO and the DNO and details the coordination under the FPC; Section 3 describes the mathematical formulation for the coordinated management of a PEV charging station under the proposed FPC.; in Section 4, the proposed optimal coordinated management method under the FPC is demonstrated through case studies using simulations with the Pacific Gas and Electricity Company (PG&E) 69-bus distribution network; and concluding remarks are given in Section 5.

Need for Coordination between the Charging Station Operator and the Distribution Network Operator
In this study, it is assumed that the stakeholders in the PEV charging business are the DNO, the CSO, and a number of PEVs.The DNO generally takes responsibility for stable and reliable operation of the distribution network, and the CSO provides charging services to PEVs at a charging station in a commercial/office building.Unlike PEV charging services for households, the CSO provides a simultaneous charging service at a commercial/industrial building for a number of PEVs.Therefore, the charging service provided to PEVs mostly operates during office hours rather than at night.For this reason, the charging demands at the station may double the current electricity loads at the peak time or when a number of PEVs simultaneously try to charge at the same location.Therefore, this type of charging service can threaten the stability of the distribution network, mostly by causing severe Energies 2016, 9, 538 3 of 15 voltage drops in the distribution network.However, the CSO currently has no responsibility for stable network operation, but simply provides charging service via the network as a large-size customer.For this reason, there is little or no motivation for the CSO to shift its charging demands to other time slots without an additional penalty or incentive scheme.Among a number of possible solutions, this study suggests a penalty contract for voltage security.

Concept of Coordinated Operation with the Flexible Penalty Contract
This study considers an FPC between the CSO and the DNO to supplement the weakness of the TOU tariff.Although the TOU is surely one solution that can induce PEV charging demands to be shifted to off-peak times from peak times, it cannot flexibly cover the voltage problems in the distribution network according to the network status.The FPC can be useful for the DNO to induce the CSO to manage its charging demands by providing information about the network status.
Specifically, under the contract scheme, the DNO calculates the maximum available power P i MAX ptq at bus i at the time slot t that ensures all bus voltages remain within the limits using day-ahead forecasted demand data.For simplicity, it is assumed that the charging service is provided in units of hours, and each time slot t is an element of set T, where T " t0, 1, 2, ¨¨¨, 23u.Additionally, it is assumed that the DNO can notify the day-ahead recommended power consumption signal P i R ptq instead of P i MAX ptq with the security margin (SM) considering the uncertainty in the charging demands of PEVs and other electricity demands.The relationship between P i MAX ptq and P i R ptq is: where SM can be set according to the DNO's attitude regarding the risk from uncertainty, and P i MAX ptq can vary depending on the hourly network conditions.
Additionally, in the case when the power demand of the charging station exceeds P i R ptq, the DNO will impose an additional penalty price β to the CSO for the amount that exceeds P i R ptq by the contract.On the contrary, if the power demand does not exceed the noticed P i R ptq, the CSO is charged only at the current TOU rate.With the FPC having the described structure, the rigidity of the TOU can be supplemented and the DNO can partly charge the CSO for the cost of upgrading the network.

Specifying the Flexible Penalty Contract by the Distribution Network Operator
This study proposes the FPC as a coordination scheme for managing the conflicts of interest between the DNO and CSO.Within the FPC scheme, the DNO notifies the CSO of the day-ahead recommended power consumption signal and imposes a penalty price β, if the charging demands at each time slot exceed the recommended value at that time.The day-ahead recommended power consumption signal provides the threshold of the FPC, and it can be considered as a consumption guideline for the CSO to avoid a low-voltage problem in the network.This signal can be calculated by obtaining the decision variable vector P MAX , which can be securely accommodated in the network and it varies according to the day-ahead forecast load data P i L ptq, Q i L ptq in Equation (5).Since this study does not focus on the day-ahead forecast of demand, it is assumed that P i L ptq, Q i L ptq are given data.The optimal solution for P MAX can be determined by choosing the best solution among feasible vectors satisfying the constraint Equations ( 5) and (6).In this study, in terms of voltage security, the criteria for choosing the best solution is designed to maximize the voltage margin (VM) of the most vulnerable bus to the low-voltage problem.In addition, the decision variable vector TP, which means the position of tap on the main transformer, and its corresponding constraint Equations (7) and (8) are considered since the DNO can also adjust the tap position of the main transformer equipped with an on-load tap changer (OLTC) to control voltages of the network.OLTC controls the secondary side voltage magnitude of the substation transformers by adjusting the tap position.Voltage magnitudes are controlled discretely according to tap step, and OLTC can generally adjust the voltage by ˘10% in 16 steps, with 17 positions [18][19][20].Therefore, the optimization problem of the DNO on specifying the FPC can be described with the objective function described in Equation (2) subject to the Equations ( 5)- (8) as follows: max where the decision variable vector P MAX , TP can be represented as: TP " r TPp0q, ¨¨¨, TPp23q s 1 where TPptq P t0, ¨¨¨, 16u subject to: ‚ Power flow constraints: ‚ Tap operational constraints of the main transformer: P l F ptq and Q l F ptq are the active and reactive power flows on the distribution line l at t, and TPptq is tap position of the main transformer.Voltage magnitude of the secondary side of the transformer is determined by the tap position TPptq, and the relationship between tap position and the voltage is represented in Equation ( 7) where V UL , V LL are upper and lower voltage limits, and N TP is the number of tap levels.Equation ( 8) prevents the excessive operation of the tap changer by limiting the total number of tap movements to less than the maximum number of tap movements NT MAX .By including this integer variable and the operational constraints, the optimization problem for the objective function in Equation ( 2) subject to Equations ( 5)-( 8) can be considered as non-linear integer programming.
Over the 24 h, P i R ptq varies depending on hourly demands of P i L ptq and Q i L ptq, and can be expressed as: As mentioned in Section 2, the penalty rate by the proposed FPC is imposed only when P i CS ptq exceeds P i R ptq, as follows: Energies 2016, 9, 538 5 of 15 where β t should be set to be larger than the difference between the TOU rate at the peak time and at the part-peak time to be effective as a coordination measure.If not, the demand shift by the TOU overwhelms the effect from β t .

Optimal Charging Management of the Charging Station Operator under the Flexible Penalty Contract Scheme
The CSO pays the cost of using electricity to provide the charging service at the TOU tariff ρ t to the DNO.In addition to the current TOU rate ρ t , the CSO should pay the penalty price β t for the excess quantity when P i CS ptq exceeds P i R ptq under the FPC scheme.The signal P i R ptq plays a role as the connector between the DNO and the CSO and is an obtained value from Equation ( 9) in Section 3.1.On the other side, the CSO generates revenue by providing the charging service to PEVs at its charging station.The charging service fee is set to be a constant value µ t because it is not the main concern of this study.As a rational profit-seeking entity, the CSO aims to maximize its daily profit under the FPC scheme.Therefore, the objective function for the optimal charging management of the CSO can be expressed as Equation (11), considering the revenues and costs, subject to the constraints given in Equations ( 15)-( 17): where: N i R ptq " , txu is the largest integer less than or equal to x subject to: On the basis of the assumption that every PEV is provided electricity with the same charging cables, the charging demand of the charging station at t, P i CS ptq, can be calculated using Equation (12), where p PEV is the rated power of the charging cable, and N i C ptq is the number of PEVs that are being charged at t. N i C ptq is expressed in Equation ( 13) as the sum of charging statuses of the PEVs, ST k,i ptq, which is a binary variable.The recommended number of PEVs being charged N i R ptq, for the charging station on bus i at time slot t can be represented as Equation (14).Additionally, N i C ptq should not exceed the number of slots procured at the charging station described in Equation (15).The state of charge (SoC) of the kth PEV, SoC k,i ptq before its arrival should be maintained at the given initial SoC, SoC k,i INI and the relationship between the SoC of the adjacent time slot t and t ´1 can be formulated as Equation (16) where ∆t is a size of a time slot, which is 1 h, and cap PEV is the capacity of the PEV's battery.Equation (17) indicates that SoC k,i ptq should fulfill the required SoC, SoC k,i REQ before the PEV's departure and SoC k,i ptq must be smaller than the upper limit, SoC MAX .

Network Configuration and Data
The proposed optimal coordinated management method under the FPC was demonstrated through simulations with the PG and E 69-bus test distribution network in [21].The structure of the network is shown in Figure 1.The substation transformer located between buses 1 and 2 is equipped with OLTC, which controls the voltage magnitude of bus 2. The operating range of the OLTC was set to vary from ´5% to 5%, in 16 steps, and the interval of tap movement was set to 0.625% [18][19][20].As shown in Figure 1, three shunt capacitors are located at buses 14, 61, and 62, with capacities of 600 kVAR, 500 kVAR, and 450 kVAR, respectively.The base values of power and voltage were set to 10 MVA and 12.66 kV, respectively.The proposed optimal coordinated management method under the FPC was demonstrated through simulations with the PG and E 69-bus test distribution network in [21].The structure of the network is shown in Figure 1.The substation transformer located between buses 1 and 2 is equipped with OLTC, which controls the voltage magnitude of bus 2. The operating range of the OLTC was set to vary from −5% to 5%, in 16 steps, and the interval of tap movement was set to 0.625% [18][19][20].As shown in Figure 1, three shunt capacitors are located at buses 14, 61, and 62, with capacities of 600 kVAR, 500 kVAR, and 450 kVAR, respectively.The base values of power and voltage were set to 10 MVA and 12.66 kV, respectively.The peak demand of the network was set to 3791 kW; the forecasted hourly demands (%) in the network are shown in Figure 2 [22].The active/reactive power demand of each bus is assumed to be distributed according to the data in [21].It is assumed that the ratio of the forecasted hourly demand is the same for all buses.

Characteristic Data of the Plug-in Electric Vehicles
This simulation considered PEVs with a battery capacity of 21.6 kWh and a charging cable with a charging rate of 3.3 kW.Table A1 in the Appendix summarizes the randomly-generated information for 100 PEVs based on Table 1, including arrival and departure time, initial SoC, and required SoC.All data were generated assuming a normal distribution.The arrival time is the time when a PEV arrives at the charging station, and the departure time is the time when a PEV leaves the The peak demand of the network was set to 3791 kW; the forecasted hourly demands (%) in the network are shown in Figure 2 [22].The active/reactive power demand of each bus is assumed to be distributed according to the data in [21].It is assumed that the ratio of the forecasted hourly demand is the same for all buses.The proposed optimal coordinated management method under the FPC was demonstrated through simulations with the PG and E 69-bus test distribution network in [21].The structure of the network is shown in Figure 1.The substation transformer located between buses 1 and 2 is equipped with OLTC, which controls the voltage magnitude of bus 2. The operating range of the OLTC was set to vary from −5% to 5%, in 16 steps, and the interval of tap movement was set to 0.625% [18][19][20].As shown in Figure 1, three shunt capacitors are located at buses 14, 61, and 62, with capacities of 600 kVAR, 500 kVAR, and 450 kVAR, respectively.The base values of power and voltage were set to 10 MVA and 12.66 kV, respectively.The peak demand of the network was set to 3791 kW; the forecasted hourly demands (%) in the network are shown in Figure 2 [22].The active/reactive power demand of each bus is assumed to be distributed according to the data in [21].It is assumed that the ratio of the forecasted hourly demand is the same for all buses.

Others
The values of the other simulation parameters are listed in Table 3.If the penalty price is smaller than the gaps between the TOU rates at each time (off-peak, part-peak, and on-peak), the effect of the penalty price is overwhelmed by the effect of the TOU scheme.In this case, the FPC will not be an effective supplementary scheme for the TOU with respect to the voltage security.Therefore, the penalty price must be greater than the difference between the TOU rates.For the same reason, once the penalty price is greater than the rate gap between the TOU rates, the coordination results of the CSO will almost be the same, regardless of the penalty price value.However, the CSO's charging service cost will increase proportionally with the penalty price.Hence, the penalty price is determined to be larger than the rate gap ($0.216/kWh) between the on-peak rate ($0.361/kWh) and the part-peak rate ($0.145/kWh).

Simulation Cases
The simulation cases were configured as listed in Table 4.In order to verify the effectiveness of the proposed FPC, a comparison was conducted among three tariff conditions in each simulation case.Three simulation Cases, A, B, and C, were configured to investigate the effects of different locations of the charging station on the simulation results.In Case B, the charging station was located at bus 61, which is closer to the main substation at bus 1 compared to the location at bus 56 in Case A. In Case C, the charging station was located at bus 9, which is farther from the main substation at bus 1, compared to bus 56 in Case A. In each simulation case, Condition I is a flat-rate tariff, which means that there is no hour-to-hour rate variation.For Condition I, the DNO has no measure to induce the CSO to shift its charging demand from peak times to other times.Condition II is the PG and E EV-9A TOU tariff.This type of tariff indirectly helps the DNO maintain network stability but cannot flexibly reflect the network status.The authors consider Condition II as a "rigid" tariff condition that is not suitable for an environment with a high penetration of PEVs.Condition III is the proposed scheme, which applies the FPC described in Sections 2 and 3 on the basis of PG and E EV-9A TOU tariff to motivate the CSO to shift its charging demand at certain times when a voltage problem is expected.

Simulation Results
The DNO's optimization problem presented in Section 3.1 is solved using MATLAB's (The MathWorks, Natick, MA, USA) built-in genetic algorithm solver, and the CSO's optimization problem presented in Section 3.2 is solved using MATLAB's built-in fmincon solver.All of the simulations were conducted using a personal computer (PC) with 3.10-GHz Intel ® Core i5-4440 central processing unit (CPU), 8 GB of random-access memory (RAM), and 64-bit Windows ® 7 operating system.Both optimization routines are time-consuming mixed-integer nonlinear problems.The actual computation times of the simulation cases are listed in Table 5, and it is found that each case takes a few hours.Since the proposed FPC scheme is designed to be conducted a day ahead, the computation times in Table 5 seem to be implementable in practice.It should be noted that more time is consumed in case B than other cases, which means that the problem becomes more difficult as the charging station is located further from the main substation (bus 1).Under Condition I in Case A (Case A-I), where the flat rate is applied, the CSO is not motivated to adjust its charging schedule because there is no difference in the electricity tariff at each time.Therefore, the CSO only considers the required energy for the PEVs when it determines its charging schedule.The profile of N i C ptq for 24 h resulting from Case A-I is shown in Figure 3a.N i C ptq is observed to exceed N i R ptq in Figure 3a, and the minimum voltage of the distribution network is below the low voltage limit V LL during the same period as shown in Figure 4.
voltage in Case A-II, with the TOU rate, is less than the minimum voltage in Case A-I, with the flat rate, as seen in Figure 4.In particular, the minimum voltage in Case A-II occurs at 12 h, unlike in case A-I, which occurs at 11 h.The reason for this result is that the charging demand at 12 h increases by a strategic act of the CSO for avoiding the peak rate of the TOU rate.Furthermore, this result implies that the charging scheduling of the CSO considering the TOU rate can aggravate the voltage condition of the network.For Condition III in Case A (Case A-III), where the proposed FPC for voltage security is assumed, in addition to the PG and E EV-9A TOU tariff, the CSO evenly deconcentrates the charging demand to other times, as shown in Figure 3c, in order to avoid the potential penalty cost that might be imposed on the excess of C () Nt for all times do not exceed R () i Nt in Figure 3c, and the low voltage problems that occur in Case A-I and Case A-II are mitigated well in Case A-III, as presented in Figure 4.
In Case B, the location of the charging station (bus 61) is set to be farther from the main substation, compared to the location of the charging station (bus 56) in Case A. This means that Case B is configured to observe the changes in the simulation results when the voltage problem becomes more severe than in Case A. As a result of the calculation of R () i Nt with the method described in Section 3.1 for Case B, R () i Nt for 24 h is represented by the dashed lines in Figure 5.It is clearly observed that the values of R () i Nt for Case B are stricter than those for Case A by comparing Figures 3 and 5.In Energies 2016, 9, x 9 of 14 voltage in Case A-II, with the TOU rate, is less than the minimum voltage in Case A-I, with the flat rate, as seen in Figure 4.In particular, the minimum voltage in Case A-II occurs at 12 h, unlike in case occurs at 11 h.The reason for this result is that the charging demand at 12 h increases by a strategic act of the CSO for avoiding the peak rate of the TOU rate.Furthermore, this result implies that the charging scheduling of the CSO considering the TOU rate can aggravate the voltage condition of the network.For Condition III in Case A (Case A-III), where the proposed FPC for voltage security is assumed, in addition to the PG and E EV-9A TOU tariff, the CSO evenly deconcentrates the charging demand to other times, as shown in Figure 3c, in order to avoid the potential penalty cost that might be imposed on the excess of C () 3c, and the low voltage problems that occur in Case A-I and Case A-II are mitigated well in Case A-III, as presented in Figure 4.
In Case B, the location of the charging station (bus 61) is set to be farther from the main substation, compared to the location of the charging station (bus 56) in Case A. This means that Case B is configured to observe the changes in the simulation results when the voltage problem becomes more severe than in Case A. As a result of the calculation of R () i Nt with the method described in Section 3.1 for Case B, R () i Nt for 24 h is represented by the dashed lines in Figure 5.It is clearly observed that the values of R () i Nt for Case B are stricter than those for Case A by comparing Figures 3 and 5.In Applying Condition II in Case A (Case A-II), where the PG and E EV-9A TOU tariff is used, the CSO shifts the charging demand from the peak time to other times because the rate at the peak time is higher than the rates of the other times.For this reason, it can be observed that the charging demand at 13 h and 14 h in Case A-I is moved to other time slots before the peak time and after the peak time, as presented in Figure 3b.Additionally, it is observed in Figure 3b that N i C ptq violates N i R ptq between 9 h and 12 h, because there is no coordination measure for delivering information on the voltage condition of the network to the CSO from the DNO.A noteworthy point is that the minimum voltage in Case A-II, with the TOU rate, is less than the minimum voltage in Case A-I, with the flat rate, as seen in Figure 4.In particular, the minimum voltage in Case A-II occurs at 12 h, unlike in case A-I, which occurs at 11 h.The reason for this result is that the charging demand at 12 h increases by a strategic act of the CSO for avoiding the peak rate of the TOU rate.Furthermore, this result implies that the charging scheduling of the CSO considering the TOU rate can aggravate the voltage condition of the network.
For Condition III in Case A (Case A-III), where the proposed FPC for voltage security is assumed, in addition to the PG and E EV-9A TOU tariff, the CSO evenly deconcentrates the charging demand to other times, as shown in Figure 3c, in order to avoid the potential penalty cost that might be imposed on the excess of N i C ptq to N i R ptq.Additionally, it is shown that N i C ptq for all times do not exceed N i R ptq Energies 2016, 9, 538 10 of 15 in Figure 3c, and the low voltage problems that occur in Case A-I and Case A-II are mitigated well in Case A-III, as presented in Figure 4.
In Case B, the location of the charging station (bus 61) is set to be farther from the main substation, compared to the location of the charging station (bus 56) in Case A. This means that Case B is configured to observe the changes in the simulation results when the voltage problem becomes more severe than in Case A. As a result of the calculation of N i R ptq with the method described in Section 3.1 for Case B, N i R ptq for 24 h is represented by the dashed lines in Figure 5.It is clearly observed that the values of N i R ptq for Case B are stricter than those for Case A by comparing Figures 3 and 5.In addition, it is observed by comparing Figures 4 and 6 that the voltage problems are more effectively mitigated in Case B than in Case A. However, the minimum voltage values of a few time slots of Case B-III are just below the voltage lower limit V LL , which is a result of the CSO's excessive charging over N i R ptq in order to comply with contract with the PEVs.In other words, even if the proposed FPC improves the voltage security and mitigates the low-voltage problem, there is still a limit of total charging loads depending on the bus location of the charging station and the distribution network condition.Nt in order to comply with contract with the PEVs.In other words, even if the proposed FPC improves the voltage security and mitigates the low-voltage problem, there is still a limit of total charging loads depending on the bus location of the charging station and the distribution network condition.In Case C, the location of the charging station (bus 9) is set to be closer to the main station compared to the location of the charging station (bus 56) in Case A. This means that Case C is configured to analyze the changes in the simulation results if the voltage problem is less serious than that in Case A. Figure 7 shows that the values of   Nt in order to comply with contract with the PEVs.In other words, even if the proposed FPC improves the voltage security and mitigates the low-voltage problem, there is still a limit of total charging loads depending on the bus location of the charging station and the distribution network condition.In Case C, the location of the charging station (bus 9) is set to be closer to the main station compared to the the charging station (bus 56) in Case A. This means that Case C is configured to analyze the changes in the simulation results if the voltage problem is serious than that in Case A. Figure 7 shows that the values of  In Case C, the location of the charging station (bus 9) is set to be closer to the main station compared to the location of the charging station (bus 56) in Case A. This means that Case C is configured to analyze the changes in the simulation results if the voltage problem is less serious than that in Case A.
Figure 7 shows that the values of N i R ptq for 24 h have more generous values in Case C compared to those in Case A. As seen in Figure 8, the mitigation of the voltage problem is not as effective in Case C as in Case A or Case B. Consequently, considering the results of Case B and Case C, it can be said that the proposed FPC becomes more effective as stricter N i R ptq conditions are imposed.

Conclusions
This study proposes an FPC for voltage security as a coordination measure between the DNO and the CSO, which helps the DNO mitigate the low-voltage problem caused by the existence of a charging station.In this paper, a method for calculating R () i Nt , and a formulation for the optimal coordinate management of the CSO under the FPC scheme are proposed.Using the PG and E 69-bus test distribution network, the proposed scheme is verified to be effective as a coordination measure to handle the voltage security.In other words, the FPC can play a role as a measure to pass on the DNO's cost for maintaining stable operation of the network to the CSO.In addition, the effect of the FPC can be altered according to the location of the charging station.Specifically, it is more effective to use the FPC when a bus is more vulnerable to voltage drops.Although this study does not obtain the optimal penalty price to maximize the social welfare under the FPC scheme, this can be addressed in depth in future work.

Conclusions
This study proposes an FPC for voltage security a coordination between the DNO and the CSO, which helps the the low-voltage problem caused by the existence of a charging station.In this paper, a method for calculating R () i Nt , and a formulation for the optimal coordinate management of the CSO under the FPC scheme are proposed.Using the PG and E 69-bus test distribution network, the proposed scheme is verified to be effective as a coordination measure to handle the voltage security.In other words, the FPC can play a role as a measure to pass on the DNO's cost for maintaining stable operation of the network to the CSO.In addition, the effect of the FPC can be altered according to the location of the charging station.Specifically, it is more effective to use FPC when a bus is more vulnerable to voltage drops.Although study not obtain the optimal penalty price to maximize the social under the FPC scheme, this can addressed in depth in future work.

Conclusions
This study proposes an FPC for voltage security as a coordination measure between the DNO and the CSO, which helps the DNO mitigate the low-voltage problem caused by the existence of a charging station.In this paper, a method for calculating N i R ptq, and a formulation for the optimal coordinate management of the CSO under the FPC scheme are proposed.Using the PG and E 69-bus test distribution network, the proposed scheme is verified to be effective as a coordination measure to handle the voltage security.In other words, the FPC can play a role as a measure to pass on the DNO's cost for maintaining stable operation of the network to the CSO.In addition, the effect of the FPC can be altered according to the location of the charging station.Specifically, it is more effective to use the FPC when a bus is more vulnerable to voltage drops.Although this study does not obtain the optimal penalty price to maximize the social welfare under the FPC scheme, this can be addressed in depth in future work.

Figure 2 .
Figure 2. Forecasted hourly demand (%) over 24 h [22].4.1.2.Characteristic Data of the Plug-in Electric VehiclesThis simulation considered PEVs with a battery capacity of 21.6 kWh and a charging cable with a charging rate of 3.3 kW.TableA1in the Appendix summarizes the randomly-generated information for 100 PEVs based on Table1, including arrival and departure time, initial SoC, and required SoC.All data were generated assuming a normal distribution.The arrival time is the time when a PEV arrives at the charging station, and the departure time is the time when a PEV leaves the charging station.The initial SoC denotes the SoC of a PEV at arrival, and the required SoC indicates the level SoC that a PEV needs to reach until its departure time.

Figure 3 .
Figure 3. C () i Nt for (a) Case A-I; (b) Case A-II; and (c) Case A-III, with R () i Nt (dashed line) for each hours.

Figure 4 .
Figure 4. Minimum voltage of the test distribution network for Case A over 24 h.

Figure 3 .
Figure 3. N i C ptq for (a) Case A-I; (b) Case A-II; and (c) Case A-III, with N i R ptq (dashed line) for each hours.

Figure 3 .
Figure 3. C () i Nt for (a) Case A-I; (b) Case A-II; and (c) Case A-III, with R () i Nt (dashed line) for each hours.

Figure 4 .
Figure 4. Minimum voltage of the test distribution network for Case A over 24 h.

Figure 4 .
Figure 4. Minimum voltage of the test distribution network for Case A over 24 h.

Figure 5 .
Figure 5. C () i Nt for (a) Case B-I; (b) Case B-II; and (c) Case B-III, with R () i Nt (dashed line) for each hour.

Figure 6 .
Figure 6.Minimum voltage of the test distribution network for Case B over 24 h.

i
Nt for 24 h have more generous values in Case C compared to those in Case A. As seen in Figure8, the mitigation of the voltage problem is not as effective in Case C as in Case A or Case B. Consequently, considering the results of Case B and Case C, it can be said that the proposed FPC becomes more effective as stricter R () i Nt conditions are imposed.

Figure 5 .
Figure 5. N i C ptq for (a) Case B-I; (b) Case B-II; and (c) Case B-III, with N i R ptq (dashed line) for each hour.

Figure 5 .
Figure 5. C () i Nt for (a) Case B-I; (b) Case B-II; and (c) Case B-III, with R () i Nt (dashed line) for each hour.

Figure 6 .
Figure 6.Minimum voltage of the test distribution network for Case over 24 h.

i
Nt for 24 h have more generous values in Case C compared to those in Case A. As seen in Figure8, the mitigation of the voltage problem is not as effective in Case C as in Case A or Case B. Consequently, considering the results of Case B and Case C, it can be said that the proposed FPC becomes more effective as stricter R () i Nt conditions are imposed.

Figure 6 .
Figure 6.Minimum voltage of the test distribution network for Case B over 24 h.

Figure 7 .
Figure 7. C () i Nt for (a) Case C-I; (b) Case C-II; and (c) Case C-III, with R () i Nt (dashed line) for each hour.

Figure 8 .
Figure 8.The minimum voltage of the test distribution network for Case C over 24 h.

Figure 7 .
Figure 7. C () i Nt for (a) Case C-I; (b) Case C-II; and (c) Case C-III, with R ()

Figure 8 .
Figure 8.The minimum voltage of the test distribution for C over 24

Acknowledgments:
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP) Contributions: Jip Kim and Seung Wan Kim conceived and designed the study.Jip Kim developed the mathematical model, setup the simulation environment, and performed the simulations.Jip Kim drafted the

Figure 8 .
Figure 8.The minimum voltage of the test distribution network for Case C over 24 h.

Table 1 .
Basic data for generating the information of 100 plug-in electric vehicles (PEVs).SoC: state of charge.Time-of-use Rates for the Charging Station OperatorIn this study, it is assumed that the DNO charges the CSO according to the PG and E EV E-9A TOU pricing scheme, especially for PEVs in the summer season, as summarized in Table2[23].

Table 5 .
Computation times of the simulation cases (s).DNO: distribution network operator; and CSO: charging station operator.
is observed by comparing Figures4 and 6that the voltage problems are more effectively mitigated in Case B than in Case A. However, the minimum voltage values of a few time slots of Case B-III are just below the voltage lower limit VLL, which is a result of the CSO's excessive charging over R ()i is observed by comparing Figures4 and 6that the voltage problems are more effectively mitigated in Case B than in Case A. However, the minimum voltage values of a few time slots of Case B-III are just below the voltage lower limit VLL, which is a result of the CSO's excessive charging over R ()i Maximum available charging power on bus i the at time slot t P i R ptq Recommended charging power on bus i at the time slot t TPptq Tap position of the main transformer at the time slot t V i ptq Voltage on bus i at the time slot t P l F ptq, Q l F ptq Active/reactive power flow on distribution line l at the time slot t P i L ptq, Q i L ptq Day-ahead forecasted active/reactive power demand on bus i at the time slot t P i CS ptq Charging demand in the charging station on bus i at the time slot t δ i ptq Phase angle on bus i at the time slot t N i C ptq Number of PEVs being charged on bus i at the time slot t N i R ptq Recommended number of PEVs to be charged on bus i at the time slot t ST k,i ptq Charging status of the kth PEV on bus i at the time slot t SoC k,i ptq State-of-charge of the kth PEV on bus i at the time slot t Penalty price of the FPC at the time slot t ρ t TOU rate at the time slot t µ t Charging service fee at the time slot t V UL /V LL Upper/Lower voltage limit G ij , B ij Conductance/susceptance of the distribution line between the i th and jth buses REQInitial/required values of the state-of-charge for the kth PEV on bus i SoC MAX Maximum value of state-of-charge for PEV T k,i ARR , T k,i DEP Arrival/departure times of the kth PEV on bus i