A Stochastic Methodology for EV Fast-Charging Load Curve Estimation Considering the Highway Traffic and User Behavior

: The theoretical impact of the electric vehicle (EV) market share growth has been widely discussed with regards to technical and socioeconomic aspects in recent years. However, the prospection of EV scenarios is a challenge, and the difficulty increases with the granularity of the study and the set of variables affected by user behavior and regional aspects. Moreover, the lack of a robust database to estimate fast-charging stations’ load curves, for example, affects the quality of planning, allocation, or grid impact studies. When this problem is evaluated on highways, the challenge increases due to the reduced number of trips related to the reduced number of charger units installed and the limited EVs range, which influence user anxiety. This paper presents a methodology to estimate the highway fast-charging station operation condition, considering regional and EV user aspects. The process is based in a block of traffic simulation, considering the traffic information and highway patterns composing the matrix solution model. Also, the output block estimates charging stations’ operational conditions, considering infrastructure scenarios and simulated traffic. A Monte Carlo simulation is presented to model entrance rates and charging times, considering the PDF of stochastic inputs. The results are shown for the aspects of load curve and queue length for one case study, and a sensibility study was conducted to evaluate the impact of model inputs.


Introduction
In recent years, public and private policies have been developed to improve global energy conditions, including transport market actions.Based on the reduction of Greenhouse Gas Emissions and the incentive for countries' energy autonomy, the electric vehicles (EV) market is developing [1].Large light-duty and heavy-duty EV occurrences are already observed in China, the USA, and European countries, reflecting a global EV stock reaching 20M during 2022, between battery EVs (BEV) and plug-in hybrid EVs (PHEV) [2].
Although the new "re-fueling" paradigm was introduced by EVs, where major charging events occur at home (50-80%), the development of public charging infrastructures (PCIs) has had impact on the EV market growth [3].Reference [4] compares the EV development and PCIs as a "chicken-or-the-egg fashion" problem, where the increase of EV stock enhances PCI implementation, which in turn impacts the EV market.
This relation is straightforward, observing the three main barriers to improving technology adherence.The first barrier, the EV battery range, is related to manufacturing technologies and traffic patterns.Second is the user anxiety range, which measures the fear level of EV users when the EV State-of-Charge (SoC) is near the minimum [5].Finally, the charging duration time impacts the total displacement time and EV user routine.According to [6], these conditions could be handled by increasing and distributing the PCIs and adopting fast-charging stations (FCSs) with power capacity greater than 50 kW.Furthermore, high-power charger units reduce the charging duration time for the same charging event energy condition, aligning it with the time to re-fuel internal combustion vehicles [7].When the PCI models are extrapolated to highways and roads, these points impact the mobility transition directly, where long-distance trips above 200 km could be concluded with energy and time safety [8,9].
However, due to the disruptive aspect of EV technologies and the development focusing on urban environments, FCSs for road or highway applications lack more specific studies about their operating conditions and business models to guarantee sustainability in this environment [10].Behind this approach's lack is the complexity to determine traffic flow patterns and, consequently, to project the FCS load curves impacted by spatialtemporal conditions [11], mainly on first-time implementation projects.
The FCS load curves estimation improves the results of an extensive planning field.For the development of FCSs, it results in a better definition of infrastructure placing and sizing [12].Moreover, due to the FCS power, new installations modify local power system load forecasting, and a disclosed forecasting model supports the FCS uncertainty representation [13].Also, grid impacts, like voltage drops or feeder overload, are better estimated and anticipated for future system actions [14].
Finally, from an economic perspective, the FCS load curves evaluate the optimized relation between FCS and the local grid on energy and power capacity aspects, influencing the operation cost and, consequently, charging prices and business models [15].The nominal power and number of charging units affect the initial costs of an FCS implementation plan.In turn, when it is evaluated from an operation view, it is related to power demand contracts according to the power available.These two parcels affect the FCS availability costs.Moreover, charging events in different day periods influence the average energy price and operational costs, mainly for time-of-use or real-time tariffs.Then, alternative load curve patterns could generate different charging prices, alternating amounts of energy consumed, or just the curve shape, influencing the relation between the user and the station and the user's own behavior.Reference [16] also highlights the impact of different tariffs' on-demand peaks, for example.When smart tariffs are adopted, the charging price must be adapted to new operational conditions.
Thus, this paper aims to develop a model to estimate the FCS hourly load curves, considering road or highway infrastructures and the spatial-temporal aspect evolving these applications.The methodology is composed of three steps.The first was modeled for user patterns recognition and determine charging duration times and FCS entry probabilities.These variables are inputs to the second step, a traffic simulation model, to evaluate the position and time that charging events will occur on the planning horizon.Finally, the third step is to sort each FCS operational routine, characterized by the charger unit uses and queue formation.
The results are FCS load and queue curves occurring in the time domain, able to identify and project regional and temporal seasonalities related to origin-destination patterns, EV market conditions, and road or highway occupancy according to hour and day type.This article presents advances in the methodology and discussion improvements from the first step of this research, considering Reference [17].

Literature Review
The infrastructure localization reflects different operational behaviors, comparing FCS infrastructures on highways or urban spaces.Generally, the uncertainty related to these load curves, mainly for low EV market share scenarios, is expressive.Some stochastic factors are related to the traffic dynamic due to random origin-destination (OD) pairs, highway seasonalities, user behaviors and their driving patterns, EV type, and the relation between user and EV SoC.For these reasons, FCS load curves for road or highway infrastructures have a spatial-temporal characterization [18].
Energies 2024, 17, 1764 3 of 27 Some theoretical studies define that the load curve projection must start from a traffic simulation result, mainly in the time domain, to include hourly and daily traffic variations [19].Moreover, it is necessary to identify FCS occupancy aspects since the number of chargers is finite, resulting in a queue composition [20].Due to the lack of user planning or anxiety range, the unbalanced use of subsequent FCSs results in queue variations [21].
The literature approaches related to PCIs and FCSs on highway load curve modeling follow two model groups.A bottom-up model is applied for scenario generation according to historical data [22].On the other hand, in top-down models, some empirical aspects need to be modeled, characterized by probability distributions functions to estimate operational elements [23].
About bottom-up methodologies, ref. [24] estimates probability distribution functions based on local data, defining the charging start and stop time and respective SoCs.Following a Big Data application, ref. [14] proposes load curves based on battery electric vehicles (BEV) data, from charging events or V2G mode.A cluster model aggregates EV data on day types, creating a decision tree solution.Reference [25] also uses Big Data, using secure information of 105,000 trips on urban ways, defining the direct relation between traffic and charging station occupation.On the other hand, following conventional load forecasting methods, ref. [13] uses historical temporal series and fARIMA methods to decompose these series to linear, seasonal, and random components, projecting the charging station load just for 120 min ahead.Finally, on consumer behavior standardization, ref. [26] proposes a model to estimate the daily load curve considering field research about EV use patterns.
A common aspect of these references is that they are applied to urban environments.The data availability enhances these studies.Moreover, traffic variations in urban zones are smaller than in road or highway zones.Due to higher EV units circulating, seasonalities, and exogenous factors, like the traffic condition, non-urban zones must be better explored to generate confident FCS load curves.
In turn, top-down models are applied to model the FCS operation condition based on probability distribution functions (PDFs) or traffic estimations indicating entry rates on FCSs.Reference [27] estimates the load curve based on EV traffic flow PDFs, applying a queue theory M/M/s model, also considering charging duration time PDFs based on historical data.On the other hand, ref. [19] uses the same queue theory M/M/s model; however, it estimates the traffic flow on a Fluid Dynamic Traffic model based on differential equations to correlate distance and average traffic speed.Reference [12] solves a very short-term FCS load forecasting based on Monte Carlo simulations (MCS) by applying a queue theory model.Also, ref. [28] estimates regional urban FCS demand-based MCS.
Evaluating these references and following a top-down approach, this paper's methodology estimates FCS on highway load curves, observing an equilibrium between the data lack, generation, and use of traffic flow information.It reduces computational efforts, similar to [19], and increases the FCS load curve granularity to observe the spatial-temporal effects on regional optics and the variations between FCSs installed along an electric route.
Reference [17] presents the background concepts of this methodology, discussing general EV user behavior and traffic flow impacts on FCS operational conditions.Departing from these preliminary discussions and the proposal of the three-step model for FCS operational conditions, this article explores all advances considering the maturity of electric mobility conditions.It includes new models for the initial SoC and anxiety range, estimation of origin-destination patterns, and improvements in Monte Carlo processing related to the problem stochasticity.Also, new parameters relating to EV users and the charging process are proposed, increasing model capacities to represent actual behaviors and increasing this proposal's result reliability.A real case study is analyzed, and the general results and sensibilities on empirical parameters are discussed.

Paper Contributions
Based on these aspects, the paper's contributions are: • An adaptative MCS model to estimate user behavior variables of FCS entry probability and charging duration times, considering regional OD patterns, regional influence zones, EV market, user anxiety range, and initial SoC condition.

•
Improvements on the traffic flow methodology, considering [17], estimating EVs traveling on the road or highway on a time domain.Since the inputs patterns are stochastic, a discrete matrix model is proposed, with adaptative sizes according to distances and average velocities.
• An FCS operation model, based on charger unit operation matrices, where it considers the charging duration frequency distribution.Like the queue theory model, but increases this approach with recognition of hourly seasonalities, using the traffic results lastly obtained and the relation between EV user and charger technologies.• To discuss variables affecting FCS load curves on highways, uncertainty points, and analysis about how to apply this load information in planning studies.Also, sensibilities on empirical parameter incentive discussions about secondary influences on FCS operation efficiency.

Methodology
The FCS on the road or highway load curves methodology is not a conventional process, mainly when load curve projections must be set on new EV environments, applying these results to expansion planning models.It occurs due to problems defining user behaviors related to the equipment (EV) and trips.In this form, besides considering the spatial-temporal characteristics and user behavior, the methodology must be highly adaptative in terms of parameters, being possible to apply to different regions and EV scenarios.
Thus, the proposed methodology is based on ten steps, divided into three blocks, corresponding to EV user behavior modeling, traffic simulation, and charging simulation.Figure 1 shows the interrelations between these steps and inputs/outputs block connections.
Firstly, the frequency distribution of charging duration time and FCS entry probabilities are processed on the EV user behavior block.The treatment of these variables is highly stochastic due to the extensive set of EV scenarios in an operational region.A Monte Carlo simulation is proposed to handle user uncertainty inputs.After, for the traffic simulation block, interest points of the road or highway are identified, looking for matrices modeling based on predictions of entries and exits at these points, distances, and average speed results in traffic flow.It is proposed as a discrete solution to the simulation horizon (T).However, the smaller the simulation step (∆t), the closer the continuity solution will be to the result.Finally, the charging simulation departs from FCS EV user entrances along the time steps, and considering local charger units, it is possible to estimate charger occupancies on the simulation horizon.Therefore, load curves and queue length are calculated assuming the temporal charger occupancy and equipment load capacity.Some boundary conditions were applied in this article's methodology.Regarding user behavior, EV users want to recharge the energy amount for the trip conclusion, minimizing the waiting time for a charging event.The small charging duration (1 to 5 min) is valid for this article.Moreover, it was proposed that all EVs connected have a power capacity equal to the charger's nominal power.Regarding traffic simulation, highway and road senses are simulated separately, considering a regular average speed and the highway speed limit.
The following subsections will better explore these three blocks, including internal steps.Moreover, Table 1 presents the article variable list, consisting of acronyms and the meaning of each variable used in this methodology.The following subsections will better explore these three blocks, including internal steps.Moreover, Table 1 presents the article variable list, consisting of acronyms and the meaning of each variable used in this methodology.

EV User Behavior
For a better comprehension of the problem treated in this paper, it is essential to define how the traffic pattern influences FCS operations.The primary key to solving this question is the EV user behavior.The user defines trip start and end points and EV battery management decisions.Therefore, modeling these aspects is fundamental for the FCS simulation.
EV user displacement patterns along a road or highway connecting n cities have some probabilities associated with departure and arrival points.This condition defines origindestination (OD) matrices composed of pair-of-points interconnection probabilities.Some places around the world have developed specific studies about OD patterns.However, due to the methodology's application range and low data availability premise, an OD calculation is proposed considering historical data of entrances/exits on roads or highways.
Directly, these entrance/exit data do not result in probabilities between points.Then, a concept of local influence (LI) is introduced.Figure 2 demonstrates these zones.Considering the city on P 2 , with a respective total population (Pop 2 ), it influences the population of P 1 and P 3 since people from towns look for services, workplaces, and leisure in big cities.This effect creates local influence zones, and displacement probabilities between P 2 and P 1 or P 3 are more significant than between P 1 and P 3 .This analysis is similar for the city on point P n−1 , according with the neighborhood points.

EV User Behavior
For a better comprehension of the problem treated in this paper, it is essential to define how the traffic pattern influences FCS operations.The primary key to solving this question is the EV user behavior.The user defines trip start and end points and EV battery management decisions.Therefore, modeling these aspects is fundamental for the FCS simulation.
EV user displacement patterns along a road or highway connecting  cities have some probabilities associated with departure and arrival points.This condition defines origin-destination (OD) matrices composed of pair-of-points interconnection probabilities.Some places around the world have developed specific studies about OD patterns.However, due to the methodology's application range and low data availability premise, an OD calculation is proposed considering historical data of entrances/exits on roads or highways.
Directly, these entrance/exit data do not result in probabilities between points.Then, a concept of local influence (LI) is introduced.Figure 2 demonstrates these zones.Considering the city on  , with a respective total population ( ), it influences the population of  and  since people from towns look for services, workplaces, and leisure in big cities.This effect creates local influence zones, and displacement probabilities between  and  or  are more significant than between  and  .This analysis is similar for the city on point  , according with the neighborhood points.Moreover, each pair of points has a connection distance (), and when the distance increases, the OD probability decreases.Inspired by Newton's theory of gravity, the population of each point and the distance between a pair , , generate an attraction OD force ( ), as follows: Based on Equation (1), the connection force is directly proportional to city's population and inversely proportional to the distance square.For example, in Figure 2, the city Moreover, each pair of points has a connection distance (d), and when the distance increases, the OD probability decreases.Inspired by Newton's theory of gravity, the population of each point and the distance between a pair x, y, generate an attraction OD force (F OD ), as follows: Based on Equation (1), the connection force is directly proportional to city's population and inversely proportional to the distance square.For example, in Figure 2, the city on P n−1 has the biggest population.However, by distance, the connection probability between P 2 and P 3 may be higher, due to the influence of city P 2 .
Considering this calculated force, the OD is divided into two matrices.The OD in matrix, whose elements indicate the entry probabilities on point x, observing next y point, and the OD out matrix, to estimate exit probabilities on the following y points, related with a previous point x.
Equations ( 2) and (3) show matrices element calculation.For the same point (x = y), both equations determine OD in or OD out elements proportional to own point historical average entrance or exit data.For (2), when x < y, the probability already estimated to y − 1, a point localized between the treated pair, is subtracted from the historical average entrance on x and the normalized interaction force (F OD ) value.For (3), the exit probability on point y reduces, when the trip starts from point x, with the distance increment between x and y, subtracting the matrix element in y + 1 by the historical average exit on y and the normalized F OD value.In addition to highway points, x and y are matrix column and row indexes, as seen in Equation ( 4).Moreover, an empirical parameter (γ x,y ) could be applied to adjust the far pair-of-points index, mainly above 600 km.If γ x,y = 1, just the OD force is considered.Then, Equations ( 5) and ( 6), represents the OD in and OD out global matrices.
x (column), In the OD in matrix (Equation ( 5)), each line is a sub-vector for each highway point.Considering the highway senses processed separately when in Point 1 (highway extreme), the probability of entrance in this point is equal to 100%.In turn, analyzing Point 2, there are entrance probabilities on Points 1 and 2. Thus, for each FCS point processing, this subvector is extracted from the matrix, and the entrance highway point is randomly defined for one MCS sample according to the cumulative distribution function of this sub-vector.The same is true for the OD out matrix (Equation ( 6)) since each line is also a sub-vector.The first line comprises all matrix terms because one EV may exit from the highway along all exits until the end of the path, considering Point 1.In this case, each point has different probabilities according to the historical data and city/town relevance.
Estimating the travel behavior characteristic, the second part of EV user behavior is characterized.Two conditions must compose the EV battery management.The first is the initial SoC, reflecting the EV range considering the first charge event before the trip.According to the distance, this value is near to 100% for highway use.Moreover, the second parameter is the anxiety range (AR).Each user has their anxiety feel, making this variable estimation difficult.In this paper, the AR is modeled as the EV SoC to stop on an FCS or at the endpoint.Based on literature reviews, the final SoC varies between 5% and 30%.
Because it is not a deterministic parameter, initial SoC and AR variables are modeled according to an equivalent probability distribution function (PDF).Since the initial SoC is near 90-100% and AR near 5-30%, a beta PDF is applied, modeling a peak near these extremes.Equation (7) presents this function based on two form parameters, α and β.
The EV type is the last variable applied to the Monte Carlo simulation.Different EV models are available, and the number will increase with the development of the EV market.Fortunately, each has its efficiency and battery size aspects, which are generally associated with price.For example, market conditions indicate that EVs with larger ranges should be more expensive due to the battery prices.Otherwise, the EV comfort and optional items also influence the cost.Thus, this variation must influence the charging necessities and duration for a highway environment.
It is proposed to identify the best-selling EVs in a study region to handle the EV types, determining the range, battery capacity, and efficiency.This ranking composes a frequency distribution and the respective cumulative distribution function, making it possible to generate one EV type for each sample.
Thus, based on these four variables, it is possible to estimate the entry probabilities of each highway FCS (EP FCS ) and the charging duration time ( CD FD ).Considering stochastic aspects around input variables, a Monte Carlo simulation (MCS) model generates these outputs, according to i scenarios.Based on random input draws, considering PDFs and FDs characteristics, an origin, destination, initial SoC, AR, and EV type is associated to each scenario.
About the charging duration time, Equation ( 8) represents the equivalent charging distance (Char dist i ).It is the sum of the distance traveled ( D tvd ) and the distance to go ( D togo , subtracted by the capable distance with initial SoC ( D start ).Then Char dist , from i, is related to the EV consume (kW/km), and charger nominal power (P char ) (kW), to estimate the charging duration time (CD), as seen in Equation (9).In (9), the sum of 1 represents that the minimum connection time is 1 min.It is subtracted in Equation ( 11) to estimate the event occurrence.
The result of processing all MCS scenarios composes the charging duration time FD.Considering the entry probability, Equations ( 10)-( 13) are applied to estimate entry probabilities (EP FCS ).
Equation ( 10) defines the entry probability as the relation between the number of charging events (N char ) subtracted from possible charging events on the following FCS in the path (N next ), divided by the total MCS scenarios (N sc ), subtracted from the not available scenarios (N na ).Equation (11) indicates that all i scenarios with equivalent charging duration time above 0 are associated with one charging event, discounting 1 in this balance, from one-minute minimum charging duration time as defined.Equation ( 12) compares the distance traveled until the FCS with the initial SoC distance to determine if the EV arrives at the FCS.Finally, Equation ( 13) estimates how many charging events of the current FCS could occur in the following FCS installed in the path.A Limit Factor (F lim ) defines how many users with capacity will leave the actual FCS to recharge at the next station.If F lim = 1, all EVs with enough range will go to the following FCS.
Algorithm 1 shows the pseudocode representation for the user behavior model, focusing on the MCS interactions aspect.Firstly, the SoC and AR PDFs, EV type ( EV FD ), and OD matrices are modeled and introduced.The parameter Sc max defines the maximum number of scenarios.Early stop is the number of consecutive interactions without output updates, related with an ε MCS error.
Energies 2024, 17, 1764 9 of 27 Algorithm 1: Pseudocode representation for the user behavior model, returning charging duration frequency distribution and entrance probability Input EV FD , SoC( f Beta (x, α, β)) , AR( f Beta (x, α, β)), OD in , OD out Input Sc max , Early stop , ε MCS While all FCS are not processed do Initialize CD FD j−1 = 0 and EP FCS j−1 = 0, j = 0, cont = 0, k = 0 While stopping criteria not satisfied do For i ← 1 to 100 do Estimate EV i , SoC i , AR i , OD in i , OD out i , based in random numbers f cud (x, 0, 1) Calculate CD i , N char i , N next i , N na i using (8,9,10,11,12) end k ← k + 100 Estimate CD FD j with CD FD j−1 and new one hundred samples CD i Estimate EP FCS j with EP FCS j−1 and new scenarios using ( 13) The number of samples (MCS scenarios) is adaptative according to the algorithm stop criteria.Each interaction comprises one hundred scenarios, and the algorithm stop occurs by two conditions: maximum scenario reached or early stop.After each interaction, the charging duration FD average, standard deviation, and FCS entrance probabilities result changes are compared with the last global result (one interaction behind).If these comparisons are behind the error (ε MCS ), it indicates low results variability, and stopping just with the Early stop criteria reaching.Thus, the higher the model's variable stochasticity, the higher the number of samples and interactions.
All FCSs on road or highway is independently treated, and the output is the CD FD matrix for the FCSs, and a vector with binomial functions with the entry probability for all FCSs.On roads or highways with two senses, it is also processed independently, aggregating the charging duration times at the end of the EV user behavior process.

Traffic Simulation Model
Since spatial-temporal aspects impact highway or road FCS load curves, the number of EVs crossing these FCS points modifies the occupation rates and the probability of a charging event.If traffic flow projections are hardly stochastic, applications of exact models, e.g., based on differential equations, can generate confident and similar results compared with reliable approximation model.However, the second application is mathematically less complex and requires less computational effort.
Thus, a discrete solution is proposed for the traffic flow on the operation planning horizon.This model solves the traffic flow problem by minute on a regular application, updating the EVs' position at each simulation step.However, it is possible to change this granularity, where smaller steps leave the model closer to a continuous solution.
Figure 3 shows a highway/road model, considering interest points.In this case, the interest points (P n ) are towns, cities, or FCSs installed on highways.They are localized on a relative position according to the highway origin (rp n ), separated by distances (d ), and traveled at an average velocity (v).
updating the EVs' position at each simulation step.However, it is granularity, where smaller steps leave the model closer to a contin Figure 3 shows a highway/road model, considering interest interest points ( ) are towns, cities, or FCSs installed on highway a relative position according to the highway origin ( ), separate traveled at an average velocity ().In Figure 3, the interval between  and  was divided Still, based on this approach, it is possible to simulate the vehicle' be demonstrated.Focusing on point  , the number of EVs depa proportional to the sum of EV arriving from the last subpoint ( the highway ( ) subtracted from EVs exiting the highway ( sents this relation, applied to the time step .In Figure 3, the interval between P n−1 and P n was divided on m subpoints (SP).Still, based on this approach, it is possible to simulate the vehicle's displacements, as will be demonstrated.Focusing on point P n , the number of EVs departing from P n (EV P n ) is proportional to the sum of EV arriving from the last subpoint ( EV SP m ), with EV entering the highway (EV in ) subtracted from EVs exiting the highway (EV out ). Equation ( 14) presents this relation, applied to the time step t.
Regarding Equation ( 14), considering a road or highway with N points, and T simulation steps (horizon), it is possible to organize this data on respective matrices.Equation (15) shows the EV P matrix, composed by N•T elements, and each column is one highway point P. The matrix has T rows, based on the simulation horizon.For example, considering a daily simulation, the matrix has 1440 rows on a minute step.
Equations ( 16) and ( 17), taking the exact dimension of (15), comprise highway EV entry and exit forecasting to all points.Like highway parameters of distance and speed, this information is a model input from external data acquisition technologies like those using "digital borders" data to estimate the relation between cities/towns and the highway.Equation (18) shows matrices sizes and Equation (19) the row dimension, estimated from the simulation horizon in hours, and ∆t, varying between 1 and 60.
. . . . . . . . . ) N Columns × T Rows (18) The last Equation ( 14) variable is the EV SP matrix, used to estimate EVs arriving at a point P. To simulate the highway flow, matrix Equation ( 20) is developed.This matrix has the same T rows, however columns (M n−1,n ) are adaptative for different pair of points.It varies as Equation (22), whose columns are proportional to the integer number from the relation between distance (d n−1,n ), pair-of-points average speed (v n−1,n ), and simulation time step (∆t).For example, considering d = 60 km, v = 60 km/h, and 1 min time step, M = 60.However, for v = 90 km/h, and the same time step and distance, M = 40.
Returning to the analysis of Equation (20), due to the relation between column and rows on each simulation step, the state of a determined matrix element m will be the same state of matrix element m − 1, in the last simulation step (t − ∆t).Thereby, for M = 60 and t = 0, the EVs on SP 1 will just arrive on interest point of P n on simulation step t = 60.
Finally, considering all traffic flow processed, EV P and the binomial PDF of FCS entry probability (EP FCS ) are g T function inputs to estimate the FCS EV entrance matrix (λ P ), as seen in Equation ( 23), representing all EV users in the simulation horizon.λ P and EP FCS have the same dimension, but only FCS points are non-zero values.All EVs on the EV P matrix are processed separately, generating a single random number from the binomial PDF.
λ P = g T (EV P , f Binomial (x, EP FCS , 1 − EP FCS )) Algorithm 2 shows the pseudocode of the traffic simulation model.Initially, two looping structures sweep all N points of the highway/road until the simulation horizon T.During the steps, EV SP is updated, resulting in the traffic flow aspect.Moreover, the EV P matrix is also updated, considering the last submatrix element from all pair of points.In the second part of the Algorithm 2 pseudocode, the FCS EV entrance matrix is calculated according to function g t .Each FCS and time step t composes one interaction of a looping structure.The algorithm estimates a true or false signal from the binomial PDF for all EVs on FCS points and increments an EV counter.After all are processed, the FCS EV user matrix is available for FCS simulation.

Fast-Charging Station Model
The FCS operation is processed on the last methodology block from Figure 1, the fastcharging station model, defining load curves and queue patterns.Firstly, the infrastructure composition is determined, mainly related to the number of chargers and their power.Since the charging duration time FD was defined considering a power parameter, this aspect must be the same for FCS processing.
Moreover, two already-processed variables are inputs of this step.A random charging duration time for each event comes from the MCS block.Otherwise, the FCS EV user matrix comprises the number of charging events for all highway FCSs modelled.Among the boundary conditions, a first-come first-serve pattern is used.Also, in this article, all EV charging events remain on chargers' power capacity and recharge the equivalent energy to end the trip, as previously mentioned in the boundary conditions.
The FCS model is composed by a matrix solution.Initially, Equation ( 24) defines the sum of FCS EV user matrices, since the road or highway is composed by two senses, direct and reverse.Then, to each FCS P n , it results in a global λ FCS P n matrix.Related to FCS infrastructure, a charger's binary matrix computes time-domain charging event, as seen in Equations ( 25)- (27).The Equation ( 25) matrix has T rows, equivalent to the simulation horizon with the same time step of previous simulation steps.In addition to this, the matrix has C columns, each one corresponding an FCS charger unit.Matrix elements assume 0, for a charger waiting status on t, or 1, for an operation condition on t, as indicated by Equation (26).
A cumulative frequency distribution is derived from CD FD P n and a random number generator returns the charging duration time from one event, in minutes, as shown in Equation (28).Moreover, a parameter introduces extra charger occupation times (δ FCS ), representing minimum intervals between two consecutive operations (EV remotion and new EV charging start).It varies with how complex the FCS operation or user technological affinity is.
After all time steps in simulation horizon processing, the FCS load curves are estimated as Equation (29).For each t time step, the sum of the product between chargers' status and chargers' power results in the time step power.Once all time steps are processed, Load FCS Pn represents FCS load curve data.In addition, during FCS processing, for a determined time step, if all chargers had been in use and one or more EVs are waiting, they are added on the next time step.Thus, Equation (30) shows that the FCS queue matrix comes from the FCS EV user matrix treated after processing ends.
This development is shown in Algorithm 3, and the road or highway FCS points are independently processed.Firstly, if an FCS is installed on P n , the chargers' matrix, FCS EV user's matrix, and charging duration PD are initialized.Then, an external looping process of all time steps on the simulation horizon and an internal looping control the EV charging events, defining charging duration and allocating them to the chargers' matrix or adding them to the next time step.Finally, based on Equations (29) and (30), the FCS load curve and the queue length matrix are obtained and available for operation or expansion planning studies.

Results and Discussion
This section develops and discusses the proposed methodology for implementing a case study, starting from base case results through sensibility analysis of model parameters.It focuses on parameters whose definitions involve empirical aspects and defines their impacts on FCS load curve results.All algorithms were implemented in Python v.3.9.12, using mainly NumPy and Pandas Libraries, and processed on an Intel Core I5 11th generation with 8 GB RAM.

Case Study
A case study from an actual highway localized in the south region of Brazil is proposed.The distance between the first and the end point of this highway is 629 km.The FCS placing is not the focus of this work, but studies like [27] consider approximately one FCS each 100 km for highway scenarios, according to EV range effects.Following these results, the case study highway starts with seven FCS.However, due to aspects such as local influences, and since this highway has two senses, more than one FCS is necessary around big cities.Therefore, 10 FCS were strategically positioned.
Moreover, thirteen heterogeneous cities exist along the highway, with different sizes, inhabitants, and economic relevance.The sum of cities and FCS points results in 23 interest points, indicated on the map of Figure 4.In this case, P5 and P15 are the main big cities, creating influence zones to the other points.For convention, the southern point of Figure 4 is the first, and the northern, P23, is the last.Also, the direct highway sense is from south to north (P1 to P23), and the opposite for the reverse sense.Regarding the city's entrance and exit database, information from the Transport Infrastructure National Department (TIND) was used [29], considering all 2019 data hourly.
Based on Figure 4, Table 2 shows the population, distances, and average speed for the proposed case study.About influence zones, P5 and P15, with 555,100 and 4,319,860 inhabitants, respectively, are the main points.Also, the traffic flow between P5 and P15 in both senses is expected to be higher than at highway extremes.For P19, P21, and P23, the number of inhabitants increases during the summer due to their being vacation destinations.This modifies patterns of traffic flow, influencing FCS occupation rates and, consequently, load curves and queue observations.Regarding the distance between FCS points, due to the allocation near existing services, mainly gas stations, for P2 and P16, they are located 5 km and 6 km above 100 km.However, due to EV range improvements, this is a low-impact addition.FCSs at P4 and P6, and P14 and P16, are near between them due to the influence of P5 and P15 cities, respectively.Finally, the average speed in case study simulations is 80 km/h between P1 and P5 and 100 km/h between P5 and P23, according to the highway limits.
In terms of case study traffic measures, Table 3 presents some 2019 highway entrance data information.Although it is impossible to estimate EV traffic patterns just by data, this methodology allows the composition of different market share scenarios and impacts on charging infrastructures.The base case, for example, departs from an EV market share of 1%.As seen in Table 3, the relevance of P5 and P15 is again clear.In either direct or reverse highway senses, main displacements occur between P5 and P15 and P15 and P23.For P18 and above, local maximums are observed in the reverse sense, with effect seen in the summer season.Conversely, P1 and P3 are towns with local relevance, so the traffic numbers decrease.FCSs behind P5 should have lower occupancy rates than those above P5.
Other used inputs for the MCS user behavior estimation are EV type PD, initial SoC, and anxiety range PDF.Table 4 and Figure 5 show these inputs considered in the case study.According to Brazilian regional sales, eleven EV models compose the average EV type model whose battery range indicates one to three charging events to complete the path between P1 and P23.In turn, the red curve of Figure 5 represents the anxiety range beta PDF, with α = 3 and β = 33, and a maximum near to 6%.The initial SoC curve (black) was parameterized with α = 28 and β = 2, and with a maximum of 96%.Both cases vary approximately 30% from the SoC extreme.These values were used for the base case, but parameter updates will be evaluated on sensibility studies.  1 Brazilian National Traffic Secretary Data. 2 ev-DataBase. 3Manufacturer's Information.

Base Case Results
Based on the proposed case study, a one-year simulation horizon in an hourly time step was conducted.FCS entry probabilities and charging duration time PD have independent results for each climate season and weekdays or weekends to estimate annual seasonalities better.Also, the base case results simulated consider a 1% EV market share (MS).Moreover, for results comprehension, FCSs were numbered according to the direct sense position, i.e., P2-FCS1, P4-FCS2, through P22-FCS10.
Figure 6 shows the FCS entry probabilities bar chart for highway senses and day types.High index variations due to traffic dynamics occur between highway senses.FCS1

Base Case Results
Based on the proposed case study, a one-year simulation horizon in an hourly time step was conducted.FCS entry probabilities and charging duration time PD have independent results for each climate season and weekdays or weekends to estimate annual seasonalities better.Also, the base case results simulated consider a 1% EV market share (MS).Moreover, for results comprehension, FCSs were numbered according to the direct sense position, i.e., P2-FCS1, P4-FCS2, through P22-FCS10.
Figure 6 shows the FCS entry probabilities bar chart for highway senses and day types.High index variations due to traffic dynamics occur between highway senses.FCS1 direct sense results are six times greater than the same point on reverse sense.This also appears in the intermediate stations, FCS4, FCS5, FCS8, and FSC9.Since the main trip start points are P5 and P15, the FCS probabilities around these points are expected to be smaller than in the central stations, except by FCS2, the last station until P1, in the reverse sense, separated 105 km from FCS1, near to P1.This simulation generally indicates a concentration of charging events on central FCSs, according to main displacement ways.About the charging duration times, Figure 7 presents FCS4 (P9) and FCS5 (P12) curves, representing the highest entry probabilities, and FCS1 and FCS10, located on the highway extremes.Except for FCS1, the other points show maximums in the first 5 min since, in some MCS scenarios, the battery range is near the displacement distance.Also, the user-modelled interest is to charge energy to conclude the trip.For example, 5 min results in 5 kWh and approximately 25 km extra for case study EVs.FCS4, 5, and 10 also present local peaks near to 10 min.On FCS4, it is related to the reverse sense since the concentration is near 25 min in the direct sense, mainly for after P15 destinations.FCS10 presents a charging duration time local peak of 50 min due to the reverse sense conditions.Although the most negligible probabilities are for FCS1, the charging duration times occur between 20 and 35 min, mainly due to small initial SoC (70%) MCS scenarios.All FCSs were processed.About the charging duration times, Figure 7 presents FCS4 (P9) and FCS5 (P12) curves, representing the highest entry probabilities, and FCS1 and FCS10, located on the highway extremes.Except for FCS1, the other points show maximums in the first 5 min since, in some MCS scenarios, the battery range is near the displacement distance.Also, the usermodelled interest is to charge energy to conclude the trip.For example, 5 min results in 5 kWh and approximately 25 km extra for case study EVs.FCS4, 5, and 10 also present local peaks near to 10 min.On FCS4, it is related to the reverse sense since the concentration is near 25 min in the direct sense, mainly for after P15 destinations.FCS10 presents a charging duration time local peak of 50 min due to the reverse sense conditions.Although the most negligible probabilities are for FCS1, the charging duration times occur between 20 and 35 min, mainly due to small initial SoC (70%) MCS scenarios.All FCSs were processed.
According to the second methodology block, Figure 8 shows the average and standard deviation results for the EV P matrix, considering the highest population points, P5 and P15, in summer (a) and winter (b) for the direct sense.Initially, based on the curve average, a low EVs/hour index is expected on 1% EV MS, considering the highway traffic data.On the presented points, the average EV curves vary between 0.05 and 1.69, and a standard deviation (shadow zones) between 3.4 and 3.1 for P15 in summer and winter, respectively.The highest EV traffic concentration occurs between 6:00 and 21:00, mainly due to sunlight displacements.
present local peaks near to 10 min.On FCS4, it is related to the reverse sense since the concentration is near 25 min in the direct sense, mainly for after P15 destinations.FCS10 presents a charging duration time local peak of 50 min due to the reverse sense conditions.Although the most negligible probabilities are for FCS1, the charging duration times occur between 20 and 35 min, mainly due to small initial SoC (70%) MCS scenarios.All FCSs were processed.According to the second methodology block, Figure 8 shows the average and standard deviation results for the  matrix, considering the highest population points, P5 and P15, in summer (a) and winter (b) for the direct sense.Initially, based on the curve  average, a low EVs/hour index is expected on 1% EV MS, considering the highway traffic data.On the presented points, the average EV curves vary between 0.05 and 1.69, and a standard deviation (shadow zones) between 3.4 and 3.1 for P15 in summer and winter, respectively.The highest EV traffic concentration occurs between 6:00 and 21:00, mainly due to sunlight displacements.Both graphs in Figure 8 show the earlier traffic increases for P5, compared with P15, since despite local influences, P15 also has regional influence on P5.For P15, concentrations between 15:00 and 17:00 occur during summer due to the vacation condition previously explained, another seasonal effect which affects standard deviation.It has similar behaviors between graphs.However, the SD is bigger for P15, increasing the uncertainty related to the traffic data.Traffic aspects for other highway points are generally similar, except for the highway region of P1 through P5, whose smaller traffic numbers, combined with a small market share, increase the day-to-day traffic curve variation.In terms of FCS load curves, it introduces higher uncertainty components.
Finally, the FCS operation model was processed considering the annual traffic simulation.Initially, since the FCS sizing is not an objective of this paper, one charger by FCS point was evaluated.An empirical study was conducted to expand this conclusion and found a maximum of 10 EVs waiting for charging criteria on quantile 0.99.That is, during 99% of the year, the queue must be smaller than 10.Thus, for a 1% EV MS, one 60 kW charger was sufficient, except FCS4, which has the highest entry probability, requiring two charger units.
Table 5 shows the total, maximum, and minimum consumption months for the FCS annual energy consumption simulated.Since the occupancy rates increase, there are direct relations between the FCS entry probabilities and the total FCS energy consumed.FCS4 and FCS5 presented the highest values, 274.58 MWh and 185.82MWh, respectively.Both graphs in Figure 8 show the earlier traffic increases for P5, compared with P15, since despite local influences, P15 also has regional influence on P5.For P15, concentrations between 15:00 and 17:00 occur during summer due to the vacation condition previously explained, another seasonal effect which affects standard deviation.It has similar behaviors between graphs.However, the SD is bigger for P15, increasing the uncertainty related to the traffic data.
Traffic aspects for other highway points are generally similar, except for the highway region of P1 through P5, whose smaller traffic numbers, combined with a small market share, increase the day-to-day traffic curve variation.In terms of FCS load curves, it introduces higher uncertainty components.
Finally, the FCS operation model was processed considering the annual traffic simulation.Initially, since the FCS sizing is not an objective of this paper, one charger by FCS point was evaluated.An empirical study was conducted to expand this conclusion and found a maximum of 10 EVs waiting for charging criteria on quantile 0.99.That is, during are at FCS4, reaching approximately 60 kW.The lowest result is for FCS1, whose average maximum varies between 3 kW and 4 kW, denoting low occupancy rates.Moreover, for all FCSs between P1 and P15 (FCS1-FCS6), the average hourly loads increase during the day until 21:00 due to the reverse sense flow after 15:00.The charging times on this highway path in a reverse sense are longer than on direct one.This explains the winter FCS4 peak near 19:00.In a global analysis, the FCS load curves follow regional traffic patterns.However, with representative queue occurrences, the load behavior stays near the maximum power for specific days despite a momentaneous traffic flow reduction.Finally, observing Figure 9 curves, occupancy rates and EV market share affect load curve uncertainty.The FCS1 load curve is deeply stochastic since it is exposed to some periods without any charging events.Then, thinking about load curve applications in planning studies, the condition introduces complex uncertainty components for the results.

Sensibilities Analysis
Despite the limited available historical data, reducing empirical influences was one aspect of this methodology proposal.However, some strategic decision-making is essential and could support applications such as load curve projections for future scenario variations.This section presents some sensibilities analysis related to these empirical parameters.
First, about the EV market share conditions.As estimated on the base case, the EV MS and the FCS entry probabilities affect occupancy rates.Thereby, as EV MS increases, current FCS infrastructure will saturate in terms of queues, making it an insoluble problem.Three MS variations were proposed: 0.1%, 1%, and 5%.Table 6 shows charger units defined for each FCS, and Figure 10 comprises a bar chart from queue quantiles 0.9 and 0.99.As seen in Table 6, in the 0.1% MS scenario, all FCSs have one charger.This also occurred for 1% MS, except FCS4, which considers two units.However, in the 5% MS scenario, only FCS1 and FCS2 remain with one charger, respecting the Q0.99 queue rule.For FCS4, it is necessary to have six chargers, and for FCS5 and FCS8, five units of 60 kW.Another conclusion is that FCS8 has smaller energy consumption than FCS9, on 1% MS, but needs more chargers for the 5% scenario, indicating a more significant number of short-time charging events.Observing Figure 10, for 1% MS, FCS9 will receive the next investment in a charger due to a limit queue on Q0.99, and the same for the 5% MS condition, together with FCS7.FCS1 did not process any queue on Q0.99 for 5% MS, becoming the last expansion point.As seen in Table 6, in the 0.1% MS scenario, all FCSs have one charger.This also occurred for 1% MS, except FCS4, which considers two units.However, in the 5% MS scenario, only FCS1 and FCS2 remain with one charger, respecting the Q0.99 queue rule.For FCS4, it is necessary to have six chargers, and for FCS5 and FCS8, five units of 60 kW.Another conclusion is that FCS8 has smaller energy consumption than FCS9, on 1% MS, but needs more chargers for the 5% scenario, indicating a more significant number of short-time charging events.Observing Figure 10, for 1% MS, FCS9 will receive the next investment in a charger due to a limit queue on Q0.99, and the same for the 5% MS Figure 11 shows the FCS4 average load curves, varying the EV MS scenarios.Comparing load curves maximum, for 0.1%, it reaches a maximum average of 10 kW, increasing to 60 kW and 220 kW approximately, for 5%.Moreover, the relation between maximum and minimum increases with MS variation, following traffic curve patterns.Also, considering the FCS sizing rule, average demands are far from the installed power for all MS scenarios.Highways localized on low variation traffic flow have better relations between sizing and queue length, avoiding early saturation.This improves the FCS load factor.
Energies 2024, 17, x FOR PEER REVIEW 23 of 28 condition, together with FCS7.FCS1 did not process any queue on Q0.99 for 5% MS, becoming the last expansion point.
Figure 11 shows the FCS4 average load curves, varying the EV MS scenarios.Comparing load curves maximum, for 0.1%, it reaches a maximum average of 10 kW, increasing to 60 kW and 220 kW approximately, for 5%.Moreover, the relation between maximum and minimum increases with MS variation, following traffic curve patterns.Also, considering the FCS sizing rule, average demands are far from the installed power for all MS scenarios.Highways localized on low variation traffic flow have better relations between sizing and queue length, avoiding early saturation.This improves the FCS load factor.The next sensibility relates to extra-time charger occupancy ( ) and how it influences the queue.Longer times between two consecutive charging events introduce idle periods.The sum of these times generates a cascade of charging event delays, affecting queue conditions.
Varying this parameter on the base case, for FCS4, the highest FCS entry probability, Table 7 returns queue quantiles for 0.85, 0.9, 0.95, and 0.99 at 1, 3, and 5 min of extra-time charger occupancy.Although these three scenarios respect the queue rule and have similar results for 1 and 3 min, the 5 min scenario shows one EV maximum on Q0.9, reaching a 5 EV maximum on Q0.99.Then, the idle time and queue length are directly proportional.The next sensibility relates to extra-time charger occupancy (δ FCS ) and how it influences the queue.Longer times between two consecutive charging events introduce idle periods.The sum of these times generates a cascade of charging event delays, affecting queue conditions.
Varying this parameter on the base case, for FCS4, the highest FCS entry probability, Table 7 returns queue quantiles for 0.85, 0.9, 0.95, and 0.99 at 1, 3, and 5 min of extra-time charger occupancy.Although these three scenarios respect the queue rule and have similar results for 1 and 3 min, the 5 min scenario shows one EV maximum on Q0.9, reaching a 5 EV maximum on Q0.99.Then, the idle time and queue length are directly proportional.Figure 12 shows the average load curves for these three scenarios.The load curve behavior has a similar pattern considering the same traffic flow simulation.Thus, occupancy rates increase on high-demand hours between 18:00 and 21:00.This period is highly susceptible to queue formation and consecutive charging events, introducing idle times.Longer idle times reduce the average hour load, consequently increasing the queue.It is possible to observe where the load maximum reduces from 60 kW to 55 kW.Technological facilities for autonomous charging coordination and user information campaigns about how to handle the chargers improve these idle times.Table 8 compares the results for both cases, considering FCS4.The new beta function maximum probability for the initial SoC was 0.85 and 0.15 for the anxiety range, updating  and  parameters.For all-day type and highway senses, results to FCS entry probability were duplicated, approximately.Thus, it affects the occupancy rate, resulting in 81% more energy consumed and increasing the queue quantiles, resulting in three EVs on Q0.85 and fifteen EVs on Q0.99, violating the FCS charger rule.Other empirical factors are the initial SoC and anxiety range PDFs.Changing α and β influences the maximum and tail forms according to the beta PDF parameters.A variation from the base case was simulated.It created a pessimistic scenario in which the user is carelessly related to the initial SoC and can assume a more apprehensive behavior in front of the minimum SoC.
Table 8 compares the results for both cases, considering FCS4.The new beta function maximum probability for the initial SoC was 0.85 and 0.15 for the anxiety range, updating α and β parameters.For all-day type and highway senses, results to FCS entry probability were duplicated, approximately.Thus, it affects the occupancy rate, resulting in 81% more energy consumed and increasing the queue quantiles, resulting in three EVs on Q0.85 and fifteen EVs on Q0.99, violating the FCS charger rule.
These results demonstrate the sensibility of FCS infrastructures to EV user behavior, mainly in the current market situation.The methodology developed can generate different operational results, being robust mainly for expansion planning studies.
The last methodology analysis is about results stochasticity.Applying the base case parameters, Figure 13 shows the FCS4 stochastic behavior load curve.It simulated ten load curve scenarios (narrow black curves), and it was possible to determine the average load curve (red line) for a random workday in summer.The average curve follows a similar pattern from previous results, with load hours mainly between 7:00 and 21:00 and increasing during the day.Due to the stochastic methodology factors, some simulations also estimated charging events during the early morning (0:00 until 6:00).Then, based on the stochastic study nature, for applications using this methodology results, it is essential to conduct an uncertainty treatment or an error impact evaluation, e.g., in cases for operation or expansion planning.The last methodology analysis is about results stochasticity.Applying the base case parameters, Figure 13 shows the FCS4 stochastic behavior load curve.It simulated ten load curve scenarios (narrow black curves), and it was possible to determine the average load curve (red line) for a random workday in summer.The average curve follows a similar pattern from previous results, with load hours mainly between 7:00 and 21:00 and increasing during the day.Due to the stochastic methodology factors, some simulations also estimated charging events during the early morning (0:00 until 6:00).Then, based on the stochastic study nature, for applications using this methodology results, it is essential to conduct an uncertainty treatment or an error impact evaluation, e.g., in cases for operation or expansion planning.

Conclusions
Worldwide mobility technological transitions have significant challenges, mainly related to current user behaviors, anxiety range, EV battery range, and charging times.Public charging infrastructures partially treat these aspects.Focusing on PCI applied to interurban displacements, trips longer than the EV range and without significant travel time impacts are possible by fast-charging station infrastructures.
The literature has been solving FCS highway placing and sizing problems.However, particularities from FCS operation are essential for operation and expansion planning in an infrastructure or power system context.The FCS load curves and queue pattern estimations are the main variables in this planning study, and a lack of methodologies to project these conditions, mainly in a low-available-data environment, was identified.
This article presented a methodology for the FCS operational projection, handling high stochastic parameters.The FCS operational condition is based on two main questions: the relationship between EV users and their vehicles and the relationship between EV users and the highway.The answer to these questions enables us to indicate when the charging event occurs, on a time scale, and how much energy will be used.

Conclusions
Worldwide mobility technological transitions have significant challenges, mainly related to current user behaviors, anxiety range, EV battery range, and charging times.Public charging infrastructures partially treat these aspects.Focusing on PCI applied to interurban displacements, trips longer than the EV range and without significant travel time impacts are possible by fast-charging station infrastructures.
The literature has been solving FCS highway placing and sizing problems.However, particularities from FCS operation are essential for operation and expansion planning in an infrastructure or power system context.The FCS load curves and queue pattern estimations are the main variables in this planning study, and a lack of methodologies to project these conditions, mainly in a low-available-data environment, was identified.This article presented a methodology for the FCS operational projection, handling high stochastic parameters.The FCS operational condition is based on two main questions: the relationship between EV users and their vehicles and the relationship between EV users and the highway.The answer to these questions enables us to indicate when the charging event occurs, on a time scale, and how much energy will be used.
The main discussions from the base case and method validation comprise EV infrastructure, EV market, and policy implications.

Impact on EV Infrastructure
The main supply or service in the EV infrastructure chain is the energy sold to the users.Then, the load curve projection input is elementary for the FCS implementation and operation planning.Obtained from this article's sensibility cases, the FCS occupation rate varies with EV user behaviors, increasing charging probabilities in an SoC-afraid region or without trip planning.Moreover, the FCS localization impacts entry probabilities, charging duration, or both.
The case study showed an FCS localized between two critical points had high probability for a charge event.Otherwise, stations on extreme points or far from cities in influence zones may have low entry probabilities but longer charging events.During the sizing plan step, regions with low entry probabilities needing fewer charger units are only partially valid.This must be regulated by the relation between charging events and charging duration.
Also, FCS load curves on highways follow similar patterns with traffic flow.This must be explored in implementation planning, considering the expected queue on a limited charger resource.In turn, on operation planning, this fact can increase the energy and power demand costs based on time-of-use or real-time energy tariffs models.The recognition of these load patterns allows for better business modeling.
Finally, FCS is deeply impacted by daily and monthly seasonalities.For example, an energy consumption variation between minimum and maximum months of 32% was identified for FCS4 in the base case.In this case, supply contracts must fill this seasonality to a global cost optimization besides the influence on FCS occupancy rates and queue length.

Impact on EV User Experience and Market
Optimizing the charging process can improve the competitiveness of EVs compared with internal combustion vehicles, mainly in this first EV wave.Guaranteeing good experiences for users can increase market expansion.Based on scenarios obtained from this article, this methodology can improve placing and sizing models, reducing charging durations and increasing infrastructure investments in high-demand regions, influencing queue formation.
Moreover, private and public campaigns could influence current or future EV users for better technology use.The results showed that a smaller initial SoC increases FCS entry probabilities, charging duration, and, consequently, queue length.Also, increasing FCS infrastructures can reduce the user anxiety range, thus better exploring the EV total range.Otherwise, facilities must be implemented for equipment and applications for better process efficiency, reducing, for example, charger idle times.

Impact on Policies
In the policy environment, FCS load curve estimation methodologies improve the reliability analysis of future EV scenarios.For example, according to sensibility results, infrastructure saturation occurs quickly, according to the EV market share growth.Then, the proposition of FCS policies for new infrastructures can be synchronized with future EV market necessities.
Regarding energy efficiency and Greenhouse Gas Emissions, the FCS load curves indicate future energy deficits and the expected growth for generation, mainly related to alternative sources or the development of new technologies for energy storage systems, due to the high power demand of FCS with various charging units.
Moreover, the results presented a high relation between FCS load curve shape and traffic patterns.Incentives on energy tariffs, passed on through the charging price by the FCS owner, can introduce new trip strategies for EV users, increasing the low demand factor of FCS and reducing the coincidence load factor with other power system loads once the power limitation for charging events on highways is not a good alternative, considering the user experience.
For power systems, the application of this methodology to project FCS placing and sizing, due to the high-power demand, can modify system expansion policies, including new transmission lines and substations to supply this new load type.
The last results demonstrated operational day stochasticity.Necessarily, all FCS load curves or post-processing models must handle uncertainties due to many EV-related noncontrollable parameters.Thus, it is a reliable tool for generating FCS operation conditions according to time step and horizon definitions in data restrictive environments.
The next step is the improvement of this methodology, monitoring and updating stochastic variables' treatment with previous experiences on FCSs already implemented around the world and expanding the problem boundary conditions.The main advances mapped are related to the charging duration, considering minimum-duration connections and the inclusion of EV-type charging load curves, as well as differences between SoC and power demand.Also, an application on FCS energy stochastic optimization related to the FCS load curves is being developed.

Energies 2024 , 28 Figure 1 .
Figure 1.Methodology block diagram proposed.Three main steps are represented: user behavior, traffic simulation, and charging simulation.

Figure 1 .
Figure 1.Methodology block diagram proposed.Three main steps are represented: user behavior, traffic simulation, and charging simulation.

Figure 2 .
Figure 2. Interaction between local and regional influence zones, related with the population along the highway/road points.

Figure 2 .
Figure 2. Interaction between local and regional influence zones, related with the population along the highway/road points.

Figure 3 .
Figure 3. Highway/road model demonstrating points, subpoints, and input/output relations for each point.

Figure 4 .
Figure 4. Case study map representation, positioning the twenty-three points modeled (thirteen cities and ten fast-charging stations).

Figure 4 .
Figure 4. Case study map representation, positioning the twenty-three points modeled (thirteen cities and ten fast-charging stations).

Figure 5 .
Figure 5. State-of-Charge and anxiety range curves modeled as a beta probability distribution.

Figure 5 .
Figure 5. State-of-Charge and anxiety range curves modeled as a beta probability distribution.

Energies 2024 , 28 Figure 6 .
Figure 6.Bar chart for the FCS entry probability index, considering both highway senses and weekday and weekend results.

Figure 6 .
Figure 6.Bar chart for the FCS entry probability index, considering both highway senses and weekday and weekend results.

Figure 8 .
Figure 8. EV_P results for P5 and P15 (highest population), considering the direct highway sense for summer (a) and winter (b) seasons.Shaded zones represent standard deviation of the curves.

Figure 8 .
Figure 8. EV_P results for P5 and P15 (highest population), considering the direct highway sense for summer (a) and winter (b) seasons.Shaded zones represent standard deviation of the curves.

Figure 9 .
Figure 9. FCS average load curve for all simulated FCS points in summer (a) and winter (b) seasons.

Figure 9 .
Figure 9. FCS average load curve for all simulated FCS points in summer (a) and winter (b) seasons.

Energies 2024 ,
17, x FOR PEER REVIEW 24 of 28Other empirical factors are the initial SoC and anxiety range PDFs.Changing  and  influences the maximum and tail forms according to the beta PDF parameters.A variation from the base case was simulated.It created a pessimistic scenario in which the user is carelessly related to the initial SoC and can assume a more apprehensive behavior in front of the minimum SoC.

Figure 13 .
Figure 13.FCS4 stochastic behavior of load curves, considering 1% of EV market share.The black lines are ten load curve scenarios, and the red line represents the scenarios' average load curve.

Figure 13 .
Figure 13.FCS4 stochastic behavior of load curves, considering 1% of EV market share.The black lines are ten load curve scenarios, and the red line represents the scenarios' average load curve.

Table 1 .
List of acronyms and variable meanings applied in this article methodology.

Table 1 .
List of acronyms and variable meanings applied in this article methodology.

Algorithm 2 :
Pseudocode of the traffic simulation model, composed by the traffic flow and the EV entrance matrix calculation.Input EV in P , EV out P , Simul h , ∆t, d, v Input f Binomial (x, EP FCS , 1 − EP FCS ) Calculate T according to (19) Generate EV P according to (15) Generate all EV SP P n−1 ,Pn according to (20) considering results from (22) For t ← 1 to T do n ← 1 While n < N do if n = 1 then EV SP P n−1 ,Pn (EV SP M (t)) ← 0 else if t = 1 then EV SP P n−1 ,Pn (EV SP M (t)) ← 0 else Update EV SP P n−1 ,Pn according to (20): EV SP P n−1 ,Pn (EV SP1 (t)) ← EV P n−1 (t − ∆t) For m ← 2 to M do EV SP P n−1 ,Pn (EV SP m (t)) ← EV SP P n−1 ,Pn (EV SP m−1 (t − ∆t)) end end Calculate EV Pn (t) according to (14), considering EV SP P n−1 ,Pn (EV SP m (t)) Update EV P according to Equation (15) considering EV Pn (t) For aux ← 1 to EV P (EV Pn (t)) do if f Binomial x, EP FCS Pn , 1 − EP FCS Pn = True then cont ← cont + 1 n ← n + 1 end end Estimate λ P according to (23): For t ← 1 to T do n ← 1 While n < N do if EV P (EV Pn (t)) > 0 then cont ← 0

Algorithm 3 :
Pseudocode of the fast-charging station model for the load curve and queue length estimation.

Table 2 .
Main parameters from the proposed case study.

Table 2 .
Main parameters from the proposed case study.

Table 3 .
Numbers from 2019 highway vehicle entrance used in the case study.

Table 4 .
Characteristics of the actual Brazilian EV market scenario.

Table 4 .
Characteristics of the actual Brazilian EV market scenario.

Table 5 .
Energy consumed for all FCS simulated in the base case.

Table 6 .
Charger units defined for each FCS point and market share scenarios.

Table 7 .
Maximum value of highest quantiles for queue length curves for FCS4.

Table 7 .
Maximum value of highest quantiles for queue length curves for FCS4.

Table 8 .
Results comparison between base and pessimist case related to user behavior parameters for FCS4.

Table 8 .
Results comparison between base and pessimist case related to user behavior parameters for FCS4.