Optimization of Electric Vehicle Charging Scheduling in Urban Village Networks Considering Energy Arbitrage and Distribution Cost

: Electric vehicles (EV) replacing the internal combustion engine vehicle may be the solution for the particulate matter (PM) 2.5 pollution issue. However, the uncontrolled charging of EVs would challenge the power system operation. Therefore, it is necessary to implement some level of control over the EV charging procedure, especially in the residential network. In this paper, an optimization of EVs charging scheduling considering energy arbitrage and the distribution network cost of an urban village environment is presented. The optimized strategy focuses on decreasing the loss of EV owners’ energy arbitrage beneﬁt, introduced as the penalty cost. Also, peak demand, power loss, and transformer aging are included in the estimation of the cost function for the distribution network. The optimization problem is solved using the genetic algorithm. As a case study, data from the urban village in Udon Thani, Thailand, are utilized to demonstrate the applicability of the proposed method. Simulation results show a reduction in the loss of energy arbitrage beneﬁt, transformer peak load, power loss and the transformer loss of life. Therefore, the application of the optimized EV charging can prolong transformer lifetime beneﬁting both the EV owner and the distribution system operator.


Introduction
Nowadays, climate change and global warming on account of greenhouse gas (GHG) emissions such as CO 2 , are negatively impacting the world environment. Moreover, CO 2 emissions in the transportation sector are also the cause of particulate matter (PM) 2.5 pollution which is generated by internal combustion vehicles. Many countries are looking into changing their policies to increase electric vehicle (EV) usage and there has been a consequent increase of EVs to over 5 million cars worldwide in 2018 [1]. However, a large increase in EVs can significantly impact the power network, especially the low voltage (LV) network, such as the urban village power system with a higher potential of energy consumption [2]. EV charging can affect the distribution network in many ways, for example, power quality issue and reliability problems [3], power fluctuation due to the uncertainty of charging [4], voltage drop due to charging demand exceeding prior design specifications [5] and network overloading because of dump charging [6,7]. Moreover, its impacts on grid asset deterioration have been reported in references [8][9][10]. The research in reference [8] presented that the impact of EVs charging demand in the city with a high annual average temperature environment can significantly reduce the loss of life of the distribution transformer. In reference [9], the simulation result illustrated

•
The loss of the EV owner's benefit from battery energy arbitrage is formulated and taken into account in the objective function, which has not been studied, in previous works.

•
To optimize EVs charging allocation in the urban village environment, the network operating limitation is considered in the optimization constraints and evaluated based on the power flow analysis.

•
The network peak demand, power loss and transformer loss of life are formulated to be the economic term which represents the network operating costs for including in the objective function. • Finally, the transformer aging is calculated using the winding thermal characteristic model where the real baseload and the actual local ambient temperature, of an urban village, are utilized to demonstrate the effectiveness of the proposed optimization.
The remainder of the paper is organized as follows. Section 2, describes the mathematical problem formulation. Section 3 explains the proposed method that includes optimization functions and the framework. Simulation results are discussed in Section 4, followed by the conclusions in Section 5.

Problem Formulation
On one hand, the EV owner likes to charge their EV with the lowest electricity costs and benefits by using the energy arbitrage; however, this may negatively impact the network as previously mentioned. On the other hand, the DSO prefers to sell electricity using distribution assets for a longer time and with minimal operational cost. In addition, factors affecting the operation cost of distributing electricity in the village network include peak demand power, power loss in the distribution system, and grid asset degradation such transformer aging. All these factors are included in our formulation as described in the next section.

Energy Arbitrage
The electricity price in the village network varies during each period of a day, therefore the price will be on the TOU. In addition, the basic need of EV owners prefers to minimize their electricity charging cost which can be achieved by scheduling their EV charging following the TOU rate. This charging method can benefit EV owners by them selling the residual energy of an EV's battery back to the distribution network by the vehicle to grid (V2G) mode at a high electricity price period during peak time. Then, the EV owner buys energy from the network by the grid to vehicle (G2V) mode with the lower electricity price in the off-peak period. This process is called energy arbitrage [21]. The profit from the energy arbitrage can be expressed as in Equation (1).
where C AB,i is the daily energy arbitrage benefit of the i-th EV in $. P d i,t and P c i,t are discharged and charged powers of the i-th EV in kW at time t, respectively. Also, ∆t is the time interval at the time t and c t is the electricity price at time t in $/kWh. T is the total number of time slots of a day.

Peak Demand
The power flow balance equations for the village network are expressed as follows: where in Equations (2) and (3), P t and Q t are total active power in kW and total reactive power in kVar at the time t, respectively. P EV i,t is charging power in kW of the i-th EV at time t. P base j,t and Q base j,t are the village baseload of bus j in kW and kVar at the time t, respectively. In addition, P loss k,t and Q loss k,t are the active power loss in kW and the reactive power loss in kVar of the distribution branch k at the time t, respectively. In Equation (4), S t is the total apparent power in kVA at the time t. P to k,t in Equation (5) and Q to k,t in Equation (6) are the power flow to the receiving end of branch k at time t in kW and kVar, respectively. P f r k,t and Q f r k,t are the power flow from the sending end of branch k at time t in kW and kVar, respectively. Also, N ev is the total number of electric vehicles while N bus and N bra are the total number of network buses and branches, respectively.
Since only one distribution transformer supplies the village network, the transformer loading at any time t can be evaluated using Equation (4). Thus, the cost of village daily peak demand can be formulated by multiplying the daily peak load with the demand charge fee, as follows: P peak = max{P 1 , P 2 , . . . , P T }, ∀t = 1, 2, . . . , T Equation (8) is provided to find the peak load power of 24 h of a day that is used to calculate the peak demand charge. Where C peak is the peak demand cost in $, P peak is the daily peak load power in kW and c d is a daily demand charge fee in $/kW which is obtained by dividing the monthly fee (C m ) by 30.

Power Loss
The second operation cost is the power loss in the distribution network that includes power loss in conductor lines and a distribution transformer. This cost can be formulated in an economic term as shown below: where, C loss is the daily power loss cost of the network in $.

Transformer Aging Model
According to the IEEE standard C57.91-2011 [22], the transformer aging is mostly estimated from the transformer insulation life. In fact, aging or deterioration of insulation is considered to be due to aging effects produced by the hottest-spot temperature. The winding hottest-spot temperature is the result of three heat components as shown in Equation (10). θ H temperature rise over the top-oil temperature at the end of the same time period. All temperatures are in degrees Celsius ( • C).
The top-oil temperature rise over the ambient temperature at the time t after a load change is given by the exponential expression function containing an oil time constant [19] as follows: where ∆θ To,U t is the ultimate top-oil temperature rise over the ambient temperature at time t and ∆θ To t−1 is the top-oil temperature rise over the ambient temperature at the previous period in • C. ∆t is the interval of evaluation in hour and τ o is the oil time constant of transformer for any load and for any specific temperature differential between the ultimate top-oil rise and the initial top-oil rise [22]. The calculation of the ultimate top-oil temperature rise over the ambient temperature (∆θ To,U t ) is based on the top-oil temperature rise over the ambient temperature at the transformer rated loading (∆θ To,rated ) as shown in Equation (12).
where k u is the ratio of transformer ultimate loading and the rated loading. R is the ratio of transformer power loss at rated load and no-load and n is an empirically derived exponent of the cooling parameter, whose value depends on the type of transformer cooling method and approximately accounted for considering the effects of change in transformer resistance with a change in load. Also, the winding hottest-spot temperature rise over the top-oil temperature at time t after load changes is given by the exponential expression function containing the winding time constant at hot spot location as in Equation (13) [19].
where ∆θ H,U t is the ultimate winding hottest-spot temperature rise over the top-oil temperature at the end of time t in • C. ∆θ H,U t−1 is the winding hottest-spot temperature rise over the top-oil temperature previously in • C and τ w is the transformer winding temperature time constant. The calculation of the ultimate winding hottest-spot rise temperature over top-oil temperature (∆θ H,U t ) is based on the winding hottest-spot rise over the top-oil temperature at rated loading (∆θ H,rated ) as shown in Equation (14).
where m is an empirically derived exponent of the cooling parameter of the transformer [22]. The relation between insulation deterioration with time and the winding hottest-spot temperature can be evaluated using an aging acceleration factor (F AA t ) [22]. The F AA t of a given hottest-spot temperature at time t, usually is based on the reference temperature of 110 • C and can be evaluated using Equation (15).
Thus, the F AA t is equal to 1, when the winding hottest-spot temperature is equal to the reference temperature. The F AA t is more than 1 when the winding hottest-spot temperature is higher than the reference temperature. This condition indicates that the transformer is experiencing accelerated aging and its lifetime is reduced. The equivalent of the aging acceleration factor and daily transformer percent loss of life over a time period of a day can be written, as in Equations (16) and (17), respectively.
Energies 2020, 13, 349 6 of 20 where F eqa is the equivalent of the aging acceleration factor and LoL d is the daily transformer percent loss of life in hours while T d is the studied time duration of a day which normally is considered to be 24 h. Also, β is the normal insulation life of the transformer, according to IEEE standard [22], the residual tensile strength (RTS) of the end point in sealed tube aging was considered at a level of 20% residual. The end of life criteria for the oil-immersed transformer is an RTS of 20% at the reference temperature of 110 • C and a normal life of 150,000 h (17.12 years). Finally, the daily transformer aging cost can be calculated as in Equation (18).
where, C TxLoL is the transformer aging cost in $. S Tx and Tx cap are the transformer power rating in kVA and the transformer installation cost per unit in $/kVA, respectively.

The Proposed Optimization Method
The goal of the proposed optimization is the minimization of the loss of EV owner's energy arbitrage benefit through optimal charging and reduction of the distribution network operation cost. The objective and constraint functions are formulated as shown below:

Objective Function
To combine the EV owner and the DSO perspectives mentioned above, control of EV charging in the distribution network is one way to benefit both the owner and the DSO. Consequently, the objective function for minimizing the related cost can be written by the following: where, C ABpen is the penalty cost of the EV owners' arbitrage benefit loss in $ while C Ab,i is the energy arbitrage of the i-th EV in $, which can be calculated using Equation (1). P ev i,t is the optimized charging power in kW for the i-th EV at time t when parking at home and c t is the electricity price rate in $/kWh at time t. The optimal charging power value can be either positive value when charging (G2V) or negative if discharging (V2G).

Electric Vehicle (EV) Charging Constraints
The optimization of EVs charging in urban village networks is performed after EV is parked at home. The ratio of the EV stored energy level and its battery capacity is defined as the state of charge (SoC), which the initial battery level of the i-th EV (SoC i,int ) at the beginning of charging which can be calculated using Equation (21). The high accuracy model of EV battery's SoC level estimation can be found in [23]. However, this model must be evaluated by using many parameters of the battery. Therefore, to reduce the complexity of this constraint, the SoC of EV battery energy formula as shown in Equation (22) is used in this work. This equation can be compared with the equation in [23] and is successfully used in reference [15].
Energies 2020, 13, 349 7 of 20 where, SoC i,depart is the battery level of the i-th EV when it departs from home and ε i is the electrical consumption rate of the i-th EV in kWh/km. d i is the daily traveling distance in km and c i,bat is the EV battery capacity in kWh. SoC i,t and SoC i,t−1 are the battery level of the i-th EV at time t and the previous time (t-1), respectively. Besides, η EV is equal to η c when EV battery operates at charging mode and η EV is equal to −1/η d if it operates at discharging mode, where η c and η d are charging and discharging efficiencies, respectively. The charging power rating of EVs in a household is limited by the nominal rated power of its on-board charger as described below: where P EV i,rated is the on-board charger power rating in kW for the i-th EV. To prolong the EV battery's lifetime, the state of charge (SoC) of each EV battery at time t, should be limited between the minimum (SoC min ) and the maximum (SoC max ) state of charge limits. Also, the final state of charge of each EV battery (SoC i, f inal ) before departing must not be below its minimum desired level (SoC i,des ) and not above its maximum level (SoC max ). Both constraints are shown in Equations (24) and (25), respectively. The minimum required charging time of the i-th EV (T ch i,min ) can be calculated using Equation (26).
In this study, both η c and η d are set to be 95%. The state of charge of each EV battery at the time t should be between SoC min and SoC max . To prolong the EV battery's lifetime, the SoC min is set to be 20% while the SoC max is set to be 95% in order to protect the overcharging (no more than 100%) of the battery. In addition, the SoC i,des of each EV is set to be SoC max when the parking period is longer than the T ch i,min . Otherwise, it must be continuously charged with a rated power and the final SoC depends on its available charging time at home.

Transformer Constraints
The transformer loss of life is related to the winding hottest-spot temperature that correlates with the transformer loading and the ambient temperature. Therefore, to protect the transformer's aging, one can control the transformer winding hottest-spot temperature to make sure that the temperature does not go over the limit. This condition can be achieved by optimizing the EV's charging load. Thus, the constraints we have used are as follows.
where, θ H max is the maximum limit of the winding hottest-spot temperature to avoid gassing in the solid insulation and the oil. According to the recommendation in [22], θ H max is set to be 140 • C for this study. In addition, S max is the maximum limit of the transformer loading to safeguard its connector and protecting switch from damage due to higher temperatures. To defer the upgrade cost of the distribution transformer, transformer overloading can be allowed up to some appropriated limit. Following the guide in IEEE standard C57.91-2011, the S max in this study is set to be 1.4 times the transformer nameplate rating in kVA.

Network Constraints
The final constraint term is the network inequality limitation as listed below: where, V j,t is the voltage of the bus j at time t while V min and V max are the lower and upper limits of the voltage level, respectively, which are set to be 0.9-1.1 p.u. of the nominal voltage. In addition, I k,t is the current of the branch k at time t and I k,max is the maximum current limit per specification of the conductor type at branch k.

Optimization Framework
The aim of the proposed charging algorithm is to find the optimal day-ahead planning of EVs charging schedule in the village network for increasing EV owner's energy arbitrage benefit and decreasing the power distribution operating cost. This concept can be applied in the real world by using the EV aggregation system [24] which is provided by the DSO/ISO of a village. The EV aggregator can communicate with an EV owner for obtaining daily traveling data such as the expected arrival time at home and departure time from home while the driven distance can be forecast by using the historical data or user provides. Also, the EV charging scheduling can be controlled through the home energy management system (HEMS) [25].
According to the objective function in Equation (19) which may have multiple locally optimal points. Also, the constraints in Equations (27)-(30) are nonlinear function which its value are evaluated based on power flow analysis (nonlinear problem). That means this is a non-convex optimization problem and cannot be solved by using the linear optimization method. While the metaheuristic optimization methods are widely known for their ability to solve the complexity problem. A genetic algorithm (GA) is a search heuristic algorithm that is inspired by the natural revolution. It is one of the popular intelligent algorithms because of its exploitation and robustness. The examples of success using the GA to solve the optimization of EV charging were presented in references [9,11,26,27]. Therefore, the GA method is selected to solve the optimization problem in this work. However, because the GA method does not guarantee the global optimal solution. The process of the proposed optimization method is operated 10 times in order to obtain the best result. The parameters of GA, such as the population size is set to be 400 with the number of maximum generations is 100 to increase the convergence rate. While, crossover and mutation rates have been determined by a trial and error method to balance the quality of the solution and computation time, which are found to be 0.8 and 0.1, respectively. A diagram of the proposed optimization process based on the GA method is presented in Figure 1. The mechanism to determine the optimal EV charging scheduling is given as follows: (1) The process begins by creating a random initial population of each individual which feasible in the region considering constraints in Equations (22) (1), peak demand (7), power loss (9) and the transformer loss of life (18) are calculated. Then these related costs are used to evaluate the fitness function in Equation (19). (4) The GA optimization process continually generates a new population using GA terminology operation as crossover, mutation and selection until one of the stopping criteria is met. In this Energies 2020, 13, 349 9 of 20 study, two stopping criteria are used. First, when the convergence tolerance is within 10 −6 , the process will be stopped. Second, the generation process reaches the maximum number. (5) Finally, the best solution of the GA shows the optimal charging power of each EV charging at any time interval. Then the EV charging and discharging schedule considering the energy arbitrage and distribution network cost for the village is used for the day ahead planning optimization.
Energies 2020, 13, x FOR PEER REVIEW 9 of 20 (3) In this step, the cost functions of the EV owners' arbitrage benefit (1), peak demand (7), power loss (9) and the transformer loss of life (18) are calculated. Then these related costs are used to evaluate the fitness function in Equation (19). (4) The GA optimization process continually generates a new population using GA terminology operation as crossover, mutation and selection until one of the stopping criteria is met. In this study, two stopping criteria are used. First, when the convergence tolerance is within 10 −6 , the process will be stopped. Second, the generation process reaches the maximum number. (5) Finally, the best solution of the GA shows the optimal charging power of each EV charging at any time interval. Then the EV charging and discharging schedule considering the energy arbitrage and distribution network cost for the village is used for the day ahead planning optimization.

A Studied Urban Village Case
To test the proposed optimization approach, we have selected the distribution radial network in Udon Thani, Thailand, which represents an urban village network. A typical single-phase distribution transformer with a rating of 50 kVA connects the main grid with the 22 kV voltage system. This transformer serves 12 households with a low voltage level of 230 V and delivers power using the 50 THW conductor line, which spans 40 m for the main feeder and 10 m for the service drop line. The network topology for this study is shown in Figure 2. In order to evaluate the transformer aging, the village baseload and ambient temperature of a day in summer (hottest) and winter (coldest) in 2018 were used for simulations. The village baseload profile with the time step of 15 min was obtained from the automatic meter recording (AMR) of the Provincial Electricity Authority (PEA) [28] with a constant power factor of 0.9 (lagging). In addition, the hourly ambient temperature data were obtained from the Thai Meteorological Department (TMD) [29], Thailand. Both village demand power and the environment temperature profile used for this study are presented in Figure  3.

A Studied Urban Village Case
To test the proposed optimization approach, we have selected the distribution radial network in Udon Thani, Thailand, which represents an urban village network. A typical single-phase distribution transformer with a rating of 50 kVA connects the main grid with the 22 kV voltage system. This transformer serves 12 households with a low voltage level of 230 V and delivers power using the 50 THW conductor line, which spans 40 m for the main feeder and 10 m for the service drop line. The network topology for this study is shown in Figure 2. In order to evaluate the transformer aging, the village baseload and ambient temperature of a day in summer (hottest) and winter (coldest) in 2018 were used for simulations. The village baseload profile with the time step of 15 min was obtained from the automatic meter recording (AMR) of the Provincial Electricity Authority (PEA) [28] with a constant power factor of 0.9 (lagging). In addition, the hourly ambient temperature data were obtained from the Thai Meteorological Department (TMD) [29], Thailand. Both village demand power and the environment temperature profile used for this study are presented in Figure 3.

EV Model
A report of the top 10 EV sales by model in the U.S. [30], is used to select 3 EV models for this study, namely the Chevy Volt [31], Nissan Leaf [32] and Tesla Model 3 [33]. Because the area of this study is focused on the urban village, therefore, we have assumed one EV per house in this urban village study. The specifications of the three selected EV types are listed in Table 1. In this work, we expect all participating EVs in a village must have completed their charging when they park at home. Since the interval time (∆t) is set to be 15 min, the total number of time slots in a day is 96. The EV charger equipment at home is AC level type 2, which operates with a unity power factor. The charging and discharging rates follow the electricity tariffs based on the TOU rate. The electricity tariff and the monthly demand charge fee for a residential network in this study are shown in Table  2. The energy rates are converted from THB to USD using the exchange rate of 1 USD equals 32 THB.  [33] 55 6.6 0.14 Table 2. Electricity tariff data from [34].

EV Model
A report of the top 10 EV sales by model in the U.S. [30], is used to select 3 EV models for this study, namely the Chevy Volt [31], Nissan Leaf [32] and Tesla Model 3 [33]. Because the area of this study is focused on the urban village, therefore, we have assumed one EV per house in this urban village study. The specifications of the three selected EV types are listed in Table 1. In this work, we expect all participating EVs in a village must have completed their charging when they park at home. Since the interval time (∆t) is set to be 15 min, the total number of time slots in a day is 96. The EV charger equipment at home is AC level type 2, which operates with a unity power factor. The charging and discharging rates follow the electricity tariffs based on the TOU rate. The electricity tariff and the monthly demand charge fee for a residential network in this study are shown in Table  2. The energy rates are converted from THB to USD using the exchange rate of 1 USD equals 32 THB.

EV Model
A report of the top 10 EV sales by model in the U.S. [30], is used to select 3 EV models for this study, namely the Chevy Volt [31], Nissan Leaf [32] and Tesla Model 3 [33]. Because the area of this study is focused on the urban village, therefore, we have assumed one EV per house in this urban village study. The specifications of the three selected EV types are listed in Table 1. In this work, we expect all participating EVs in a village must have completed their charging when they park at home. Since the interval time (∆t) is set to be 15 min, the total number of time slots in a day is 96. The EV charger equipment at home is AC level type 2, which operates with a unity power factor. The charging and discharging rates follow the electricity tariffs based on the TOU rate. The electricity tariff and the monthly demand charge fee for a residential network in this study are shown in Table 2. The energy rates are converted from THB to USD using the exchange rate of 1 USD equals 32 THB.  Table 2. Electricity tariff data from [34]. The uncertainties associated with EVs travel such as departure time, arrival time, and daily traveling distance are applied in this study through random generation of data based on the statistical probability of the traveling pattern from the National Household Travel Survey (NHTS) data [35]. The normal distribution is used to create the random set of the departure time and the arrival time. The departure times from home are randomly generated with mean µ = 7.0 and standard deviation σ = 1.5 while the arrival times at home are randomly generated with the mean µ = 18.0 and the standard deviation σ = 3.0. Based on reference [36], the random values for the daily driven distance is generated using a lognormal distribution with mean µ = 3.2 and standard deviation σ = 0.88. The probability distribution function (PDF) of the departure time, the arrival time and the daily driven distance are generated using the MATLAB probability density function shown in Figures 4 and 5, respectively.

Period
Energies 2020, 13, x FOR PEER REVIEW 11 of 20 The uncertainties associated with EVs travel such as departure time, arrival time, and daily traveling distance are applied in this study through random generation of data based on the statistical probability of the traveling pattern from the National Household Travel Survey (NHTS) data [35]. The normal distribution is used to create the random set of the departure time and the arrival time. The departure times from home are randomly generated with mean μ = 7.0 and standard deviation σ = 1.5 while the arrival times at home are randomly generated with the mean μ = 18.0 and the standard deviation σ = 3.0. Based on reference [36], the random values for the daily driven distance is generated using a lognormal distribution with mean μ = 3.2 and standard deviation σ = 0.88. The probability distribution function (PDF) of the departure time, the arrival time and the daily driven distance are generated using the MATLAB probability density function shown in Figures 4 and 5, respectively.  The village baseload at each house in summer and winter cases, EV model and its initial SoC of each EV battery are shown in Table 3. The baseload power of each house is obtained from the monthly electricity bill while the daily load profile of an individual house is generated using the similar load pattern of the transformer. The peak load of each house is multiplied with a normalized curve of the transformer daily load profile to represent the daily load curve of each house.  The uncertainties associated with EVs travel such as departure time, arrival time, and daily traveling distance are applied in this study through random generation of data based on the statistical probability of the traveling pattern from the National Household Travel Survey (NHTS) data [35]. The normal distribution is used to create the random set of the departure time and the arrival time. The departure times from home are randomly generated with mean μ = 7.0 and standard deviation σ = 1.5 while the arrival times at home are randomly generated with the mean μ = 18.0 and the standard deviation σ = 3.0. Based on reference [36], the random values for the daily driven distance is generated using a lognormal distribution with mean μ = 3.2 and standard deviation σ = 0.88. The probability distribution function (PDF) of the departure time, the arrival time and the daily driven distance are generated using the MATLAB probability density function shown in Figures 4 and 5, respectively.  The village baseload at each house in summer and winter cases, EV model and its initial SoC of each EV battery are shown in Table 3. The baseload power of each house is obtained from the monthly electricity bill while the daily load profile of an individual house is generated using the similar load pattern of the transformer. The peak load of each house is multiplied with a normalized curve of the transformer daily load profile to represent the daily load curve of each house. The village baseload at each house in summer and winter cases, EV model and its initial SoC of each EV battery are shown in Table 3. The baseload power of each house is obtained from the monthly electricity bill while the daily load profile of an individual house is generated using the similar load pattern of the transformer. The peak load of each house is multiplied with a normalized curve of the transformer daily load profile to represent the daily load curve of each house.

Simulation and Scenarios
Based on the village network topology, EV models and the studied data as mentioned before, the power flow analysis was computed using the Newton Raphson method. Next, MATLAB is used for simulation to demonstrate the ability of the proposed optimization. To evaluate the transformer loss of life, the transformer thermal parameters used in this study are listed in Table 4. For simulations, existing EV datasets for both summer and winter seasons are used. Finally, the numerical simulations were done for 3 scenarios as explained below: • Case I: Dump charging (uncontrolled charging), each EV starts charging immediately when it arrives at home with a rated charging power without V2G operation.

•
Case II: TOU charging (EV owner's perspective), every EV will be delayed for charging for low electricity TOU rates. In addition, the discharging mode will be used at the high TOU rate period.

•
Case III: Optimal charging, all EVs are aggregated to charge and discharge with the proposed method. Table 4. Transformer parameters data from [37].

Impact on Transformer Aging
Bar graphs in Figure 6 illustrate the EVs charging scheduling in 24 h for each scenario. The transformer power was calculated using Equation (4) which is shown its load profiles in Figure 7. In addition, the transformer aging acceleration factor (F AA t ) profile versus the EV's charging load is presented in Figure 8, and the correlation of the transformer loss of life (LoL) and the winding hottest-spot temperature (θ H t ) is shown in Figure 9. In EV charging with the dump charging method, results show that when the EVs are connected to the village network with uncontrolled charging, the transformer peak load is increased to 95.56 kVA (1.91 p.u.) due to dumping of the EVs charging load during a day in summer and 90.87 kVA (1.82 p.u.) in winter. Thus, the winding hottest-spot temperature exceeds 110 • C, as a result, F AA t is more than 1 for both the summer and winter days. Under this condition, the transformer experiences accelerated aging and its lifetime is reduced. In Case 2, when EVs are charged following TOU charging, this reflects the basic preference of EV owners who want to make profit by the energy arbitrage. Without any optimization, the EV that has the capability to discharge its stored energy is operated in V2G mode during the high rate period (e.g., 9 a.m. to 10 p.m.). Consequently, the large amount of EV discharging (V2G) can cause reverse power flow between 7 p.m.-10 p.m., which in turn, may create overvoltage related issues in the village network. Moreover, all EVs status will change to dump charge when the tariff rate is low, which occurs at 10 p.m. This leads to the transformer peak load at 112.78 kVA (2.26 p.u.) during a summer day and 104.62 kVA (2.09 p.u.) during a winter day. This condition rapidly increases the winding hottest-spot temperature over 250 • C during a summer day and over 190 • C during a winter day. Consequently, it leads to the extreme raising of the accelerated aging factor, which ultimately reduces the transformer lifetime. As a result, the insulation is gassed, which can lead to an explosion.
In the case of optimal charging, results show that the proposed EV charging method using optimal scheduling of both charging and discharging can minimize the transformer peak load to 68.81 kVA (1.38 p.u.) in summer and 65.22 kVA (1.30 p.u.) in winter. In addition, the transformer winding hottest-spot temperature is lower than 140 • C during the day of summer and lower than 110 • C during the day of winter. Accordingly, F AA t declines indicating that the proposed optimization can prolong the transformer aging. Results of transformer loading and its winding hottest spot temperature were evaluated and the summary of the transformer aging including the equivalent of the aging acceleration factor and the percent of daily transformer loss of life for each scenario in summer and winter season are presented in Table 5. Consequently, it leads to the extreme raising of the accelerated aging factor, which ultimately reduces the transformer lifetime. As a result, the insulation is gassed, which can lead to an explosion.
In the case of optimal charging, results show that the proposed EV charging method using optimal scheduling of both charging and discharging can minimize the transformer peak load to 68.81 kVA (1.38 p.u.) in summer and 65.22 kVA (1.30 p.u.) in winter. In addition, the transformer winding hottest-spot temperature is lower than 140 °C during the day of summer and lower than 110 °C during the day of winter. Accordingly,  Table 5.   The TOU charging method results in the highest of the equivalent of the aging acceleration factor and the percent transformer loss of life, followed by the dump charging method. Compared to TOU and dump charging, the proposed optimal charging method can reduce the aging acceleration factor with the lowest percent transformer loss of life.

Cost Evaluation
This subsection presents the simulation results estimating the penalty cost of the EV owners' arbitrage benefit loss, peak demand cost, power loss cost and the transformer loss of life cost under all three scenarios. All these costs are evaluated based on per day and the daily energy arbitrage ( ) values that are obtained from the TOU charging method case. In addition, the transformer capital cost ( cap Tx ) is set to be 166.1 $/kVA [19]. A summary of the economic cost evaluation is presented in Table 6. Also, the efficiency improvement of the proposed optimization method based on the cost evaluation compared with the dump charging method and TOU charging method is demonstrated in Figure 10.  The TOU charging method results in the highest of the equivalent of the aging acceleration factor and the percent transformer loss of life, followed by the dump charging method. Compared to TOU and dump charging, the proposed optimal charging method can reduce the aging acceleration factor with the lowest percent transformer loss of life.

Cost Evaluation
This subsection presents the simulation results estimating the penalty cost of the EV owners' arbitrage benefit loss, peak demand cost, power loss cost and the transformer loss of life cost under all three scenarios. All these costs are evaluated based on per day and the daily energy arbitrage (C AB,i ) values that are obtained from the TOU charging method case. In addition, the transformer capital cost (Tx cap ) is set to be 166.1 $/kVA [19]. A summary of the economic cost evaluation is presented in Table 6. Also, the efficiency improvement of the proposed optimization method based on the cost evaluation compared with the dump charging method and TOU charging method is demonstrated in Figure 10.  In the case of dump charging, the penalty cost and the power loss cost are the highest when compared to the other charging methods because the energy arbitrage opportunity by V2G operation is not considered. The TOU charging case shows none of the penalty costs, which means EV owners can profit from the energy arbitrage using the V2G mode, but the peak demand cost, the transformer loss of life cost and the total cost are the highest. For the EV optimal charging case, results show that it induces a penalty cost; however, this value is lower than the dump charging case. This is because, in the proposed optimization method, we have included objective and constraint functions to minimize peak load, power loss and slow down the transformer aging. Hence, the transformer loss of life by the optimized case is the lowest among the EV charging cases. It can be observed that in both summer and winter the costs of peak demand, power loss and the transformer loss of life of the proposed optimized EV charging are the lowest compared to other charging methods. Also, the overall cost of the optimal charging method in both season cases is the lowest. However, due to the addition of penalty cost being presented, the EV owners' arbitrage benefit loss. Thus, the DSO or ISO could manage this cost with the suitable electricity tariff program to motivate the EV owner to participate in this optimal charging strategy.  In the case of dump charging, the penalty cost and the power loss cost are the highest when compared to the other charging methods because the energy arbitrage opportunity by V2G operation is not considered. The TOU charging case shows none of the penalty costs, which means EV owners can profit from the energy arbitrage using the V2G mode, but the peak demand cost, the transformer loss of life cost and the total cost are the highest. For the EV optimal charging case, results show that it induces a penalty cost; however, this value is lower than the dump charging case. This is because, in the proposed optimization method, we have included objective and constraint functions to minimize peak load, power loss and slow down the transformer aging. Hence, the transformer loss of life by the optimized case is the lowest among the EV charging cases. It can be observed that in both summer and winter the costs of peak demand, power loss and the transformer loss of life of the proposed optimized EV charging are the lowest compared to other charging methods. Also, the overall cost of the optimal charging method in both season cases is the lowest. However, due to the addition of penalty cost being presented, the EV owners' arbitrage benefit loss. Thus, the DSO or ISO could manage this cost with the suitable electricity tariff program to motivate the EV owner to participate in this optimal charging strategy.

Conclusions
This paper presents the optimization of the EV charging method in an urban village network considering energy arbitrage and network operating cost. The focus of this optimization work is to minimize the loss of the EV owner's benefit, peak demand cost, power loss cost, and the transformer aging cost. The loss of the EV owner's benefit from battery energy arbitrage is introduced as the penalty cost function. The proposed optimization process includes distribution transformer aging affected by EV charging load. Also, the transformer loss of life model is formulated based on the winding hottest-spot temperature rising mechanism. Besides, the network operation limitations such as bus voltage, branch current and the transformer loading are set as constraints and are obtained from the power flow analysis result. To validate the outcomes of the proposed optimal EV charging method, a typical distribution radial network in Udon Thani, Thailand with real village baseload and the local ambient temperature was selected for simulations.
The numerical simulations were undertaken using MATLAB. Results show that the negative impacts due to the EVs dump charging and the TOU charging method on distribution network that consists of the network overloading, increased power loss, and the transformer aging, can be mitigated by the proposed optimization charging strategy. Also, the proposed optimization of the EVs charging method presents the lowest of the peak demand cost, power loss cost and the transformer loss of life cost when compared with other charging methods. In addition, the result shows that the loss of the EV owner's energy arbitrage benefit can be reduced and the overall network operation cost is minimized. The simulation results prove that the proposed optimization method based on the genetic algorithm works efficiently for a variety of EV models, different initial EV battery energies, arrival times and departure times. Moreover, the robustness of the algorithm is presented when it performs with the different seasons of the study village case. Where the proposed optimization EVs charging shows high performance in both season cases. In summary, the application of the proposed optimal EV charging strategy can benefit both the EV owner and the distribution system operator (DSO) or independent system operator (ISO).