A Hybrid Dynamic System Assessment Methodology for Multi-Modal Transportation-Electrification

In recent years, electrified transportation, be it in the form of buses, trains, or cars have become an emerging form of mobility. Electric vehicles (EVs), especially, are set to expand the amount of electric miles driven and energy consumed. Nevertheless, the question remains as to whether EVs will be technically feasible within infrastructure systems. Fundamentally, EVs interact with three interconnected systems: the (physical) transportation system, the electric power grid, and their supporting information systems. Coupling of the two physical systems essentially forms a nexus, the transportation-electricity nexus (TEN). This paper presents a hybrid dynamic system assessment methodology for multi-modal transportation-electrification. At its core, it utilizes a mathematical model which consists of a marked Petri-net model superimposed on the continuous time microscopic traffic dynamics and the electrical state evolution. The methodology consists of four steps: (1) establish the TEN structure; (2) establish the TEN behavior; (3) establish the TEN Intelligent Transportation-Energy System (ITES) decision-making; and (4) assess the TEN performance. In the presentation of the methodology, the Symmetrica test case is used throughout as an illustrative example. Consequently, values for several measures of performance are provided. This methodology is presented generically and may be used to assess the effects of transportation-electrification in any city or area; opening up possibilities for many future studies.


The Emergence of Electrified Transportation
Electrified transportation, be it in the form of buses, trains, or cars are an emerging form of mobility.While electric trains have been an important form of urban transport for decades, it is the adoption of electric vehicles (EVs) as cars and buses that is set to expand the importance of electrified transportation in terms of distance traveled and energy consumed [1].The shift to EVs as an enabling technology supports CO 2 emissions reduction targets [2,3].Relative to their internal combustion vehicle (ICV) counterparts, EVs are more energy efficient and consume less energy per unit distance [4].They also have the added benefit of not emitting any carbon dioxide in operation and rather shift their emissions to the existing power generation fleet.As the power generation portfolio gains a greater penetration of renewable energy sources, the "well-to-wheel" carbon footprint of the EV approaches neutrality [5].ICVs, in contrast, must emit CO 2 not just for gasoline refining but also in the combustion of the fuel during transportation [2].While the reduced well-to-wheel emissions of EVs drives renewable energy penetration, EVs have also been envisioned as an enabling technology to support renewable energy penetration [6].These "Vehicle-to-Grid" applications provide demand side management of battery charging/discharging for enhanced control of the power grid balance [7][8][9][10][11].

Infrastructure Considersations in Electrified Transportation
Despite these achievements, the true success of EVs depends on their successful integration with the infrastructure systems that support them.To that effect, EVs interact with three interconnected systems: the (physical) transportation system, the electric power grid [9], and their supporting information systems.In the transportation domain, such information systems are often called intelligent transportation systems [39][40][41][42][43][44][45][46], while in the electric domain they are often called energy management systems [47,48].In this paper, their combined functionality is referred to as Intelligent Transportation-Energy System (ITES) [49,50].Successful EV integration requires an assessment with regards to all three systems and their associated interdependencies.
With respect to the transportation system, EVs behave differently from ICVs in three regards.First, EVs typically have a travel range of approximately 150km [27].Second, while ICVs can refuel in a matter of minutes, a typical EV may require 6-8 hours in order to fully recharge [51].Finally, these two aspects of EVs can be further exacerbated by the geographic sparsity of charging stations causing drivers to go out of their way to charge [49].These three aspects can lead to significantly altered user driving patterns and potentially different aggregate traffic behavior.From a private user perspective, the lack of perfect availability may be inconvenient.From a commercial user perspective, such driving patterns may erode the business case for adoption [25,49,52].
With respect to the electric power system, multiple aspects have to be considered for successful integration [53][54][55][56][57][58][59][60][61].It is generally accepted that most EV adoption scenarios will not place excessive demands on the national power generation capacity [6,50].Nevertheless, it is very likely that EV penetration can place excessive power demands on electric distribution which may exceed transformer ratings or cause undesirable line congestion and voltage deviations.These demands can be further exacerbated if EV adoption is dense geographically at the neighborhood length scale.Furthermore, if users adopt similar charging patterns, driven perhaps by similar work and travel lifestyles, the power required for charging can be temporally concentrated.Behr [62], Deutch and Moniz [22], and Markel [23] all state that one central challenge in the upgrade of the two physical systems is the EV charging infrastructure.EV charging stations, the analog of gasoline stations, serve as origin-destination nodes in a transportation system while simultaneously acting as load nodes in an electrical grid [63].This coupling requires careful implementation [63,64].In contrast, "online EVs" have been advanced as a new technology in which vehicles are inductive charged wirelessly while still in motion.This novel solution couples transportation links to electrical grid nodes in a new electrified transportation infrastructure that has the potential to assure vehicle availability [50].
Finally, with respect to information systems, Galus et al. [59], Junaibi [52], and Farid [50,63] state that the demands that EVs place on the physical infrastructure networks imply the need for informatic coordination.A parking lot operator, for example, may wish to use its energy management system to apply coordinated charging for its private EVs [65][66][67][68][69][70].Furthermore, Gong et al. [71] states that a public transport operator might integrate the EV state of charge (SOC), its remaining available range, or even the real-time electricity spot price into its intelligent transportation system so as to co-optimize the dispatch of its transportation and charging services.Finally, electric utilities may wish to implement a vehicle-to-grid control scheme to optimize power grid balance, losses, congestion, and voltage deviations in the distribution management system [7,72,73].In all, this information integration can be classified into five distinct ITES decisions summarized in Table 1 [50].These functionalities, together, can serve to minimize operating costs, enhance power grid reliability, and enhance traffic congestion.
Traffic congestion and the effect it may have on the electric power grid requires dynamic modelling.A review of modeling approaches assessing the integration of electrified transportation has been found to, in general, use static analysis methods [74].Although a static analysis may provide some insights into impacts of electrification, it does not acknowledge the interdependence of vehicles within a transportation system and how this would translate to grid impacts.

•
Vehicle Dispatch: When a given EV should undertake a trip (from origin to destination) • Route Choice: Which set of roads and intersections it should take along the way

Original Contribution and Scope
This papers contributes a hybrid dynamic system assessment methodology for multi-modal transportation.It represents a natural evolution of the method originally developed in [75].At the heart of the assessment methodology is the recently developed hybrid dynamic system (mathematical) model for multi-modal transportation electrification [63].Such a system is called a transportation-electricity nexus (TEN).Definition 1. TEN [76,77]: A system-of-systems composed of a system with the artifacts necessary to describe at least one mode of electrified transport united with an interdependent system composed of the artifacts necessary to generate, transmit, distribute, and consume electricity.
The usage of the TEN model facilitates the identification of input data, system structure, system behavior, decision variables, and output data relevant to the TEN.The mathematical model also facilitates the presentation of the methodology in a general fashion which furthermore facilitates its application on newly planned as well as existing infrastructure.Given the emergence of EVs in transportation electrification, the methodology inclines towards supporting EV integration studies.However, given that the underlying mathematical model is multi-modal, the methodology is viewed as such as well.The research towards this methodology is conducted by first investigating the most recent published EV integration studies; Second, the most appropriate methods to conduct such a study were analyzed, and subsequently extended to incorporate dynamic decision-making and performance measures.To illustrate the concepts explained in the paper, the Symmetrica test case [76] is used as an example throughout.The methodology is validated through a simulation of the Symmetrica test case with 50% adoption of plug-in EVs.
This methodology, to our knowledge, is the first of EV integration studies that presents a dynamic model.It takes into account the dynamic behavior present in any transportation system.This is a significant development from the methodologies within published literature, as generally static analysis is used [74].The importance of such an assessment methodology can be understood by analogy from the precedent of renewable energy integration.At first, the integration of solar PV and wind energy was viewed from the perspective of small demonstration projects that had little or no effect on holistic power grid performance.As desired penetration rates grew into double-digit percentages, many renewable energy integration (case) studies [78][79][80][81][82] were developed; often with inconsistent methodological formulations.This work seeks to develop this methodological foundation for electrified transportation as it becomes an integral part of the coming sustainable energy transition.

Paper Outline
The remainder of the paper is devoted to the explanation of the methodology in Section 2. A summary of results achieved with this methodology using the Symmetrica test case [76] is presented in Section 3. The paper is brought to a conclusion in Section 4.

Transportation Electrification Assessment Methodology
This section forms the body of the paper and presents the transportation electrification assessment methodology.After a brief overview of the methodology in Section 2.1, Section 2.2 establishes the TEN structure, Section 2.3 presents the behavior of the TEN, Section 2.4 presents the ITES decision-making within the TEN dynamic system, and Section 2.5 presents performance measures to assess the performance of the TEN.As mentioned in Section 1.3, the methodology relies on a hybrid dynamic system model for multi-modal transportation electrification [63].

Overview
To begin the assessment methodology, Definition 1 of a TEN can be depicted graphically as in Figure 1.From a structural perspective, it consists of a transportation topology, a power system topology and a charging topology that couples them energetically and spatially.Each of these has their dynamic behaviors.For example, the transportation system may have mesoscopic or microscopic traffic behaviors [83,84].It responds to exogeneous use case driven traffic demand in the form of itineraries that vehicles take over the course of the day [85].The power system obeys Kirchoff's Laws through power flow analysis [85].It responds to exogeneous power demands (unrelated to transportation).The charging system draws electricity from the power system to charge vehicles in the transportation system.The structure and dynamics of such a TEN has been studied in detail as a hybrid dynamic model [63].These physical dynamics also respond to ITES decision-making.The types of decisions have been identified in Table 1.While these decisions are spatially and energetically coupled, they may be taken in an integrated or independent fashion depending on the circumstances of the system being studied.Such decisions send signals from the ITES functionality down to the TEN's physical system.Furthermore, such decisions can be conducted statically in an open-loop fashion or dynamically in a closed-loop fashion where the states and outputs of the physical dynamics are taken into account in the decision making.For example, dynamic routing [86] can be distinguished from static routing in that the former may take into account traffic congestion to optimize travel time.The presence of optimization objectives (e.g., travel time) suggests that the states and output of the TEN can be used to calculate performance measures of interest.
Consequently, the transportation electrification assessment methodology executes a numerical simulation of the TEN as a hybrid dynamic mathematical model [63] with its integrated ITES [49,50] within an operations time scale.This methodology requires four steps: (1) Establish the TEN structure.
(3) Establish the TEN ITES decision-making.(4) Assess the TEN performance by numerical simulation.
Each of these steps is discussed in the following associated subsections.To facilitate the general discussion, each step is also illustrated using the Symmetrica test case [76].It consists of a multi-modal electrified transportation system topology, an electric power topology, and activity-based use case data that spans transportation and charging [76].

Methodology
The first step in the assessment methodology is to establish the TEN's physical structure.In transportation and power systems, graph theoretic models [87][88][89] are often proposed as a modeling analytical framework.A graph (be it directed or undirected) D = {B, E} consists of nodes B and edges E. In transportation systems, the nodes often physically represent intersections and stations while edges represent roads, rails or transportation routes [85,86,[90][91][92].In power systems, the nodes often physically represent generators, substations, and loads while the edges represent the power lines [93,94].The relationships between these nodes and edges are then analyzed using incidence and adjacency matrices.Despite its strength in modeling individual engineering systems, graph theory has significant limitations when two network systems are connected together.A recent review on multi-layer networks shows that such models often place limiting assumptions on how the networks may be coupled together [95].From a practical electrified-transportation engineering perspective, such modeling limitations would constrain the charging topology in a TEN.To overcome this challenge, TEN's structure is modeled using Axiomatic Design for Large Flexible Engineering Systems [96][97][98][99] as a type of hetero-functional graph theory.Its advantage relative to more traditional graph theoretic works has been previously motivated [63,98,99].Definition 2. Large Flexible Engineering System (LFES) [98,100]: an engineering system with many system processes (P) that not only evolve over time, but also can be fulfilled by one or more system resources (R).
The establishment of the TEN's structure is captured in an LFES Knowledge Base.[96][97][98]: A binary matrix J S of size σ(P) × σ(R) whose element J S (w, v) ∈ {0, 1} is equal to one when action e wv (in the SysML sense [101]) exists as a system process p w ∈ P being executed by a resource r v ∈ R. σ() gives the size of a set.Essentially, the system knowledge base itself forms a bipartite graph which maps the set of system processes to their resources.Each filled element in the system knowledge base is called a structural degree of freedom and their number describes how many "capabilities" exist in the system [96][97][98].

Definition 3. LFES Knowledge Base
Axiomatic Design for LFES also classifies system processes [96][97][98] into three varieties: transformation (P µ ), transportation (P η ), and holding(P γ ).It also classifies resources into three varieties: transforming resources (M), independent buffers (B), and transporting resources (H).Table 2 identifies the various types of system processes and resources in a TEN [96][97][98].In order to construct the TEN knowledge base, the knowledge bases for the electrified transportation system and the electric power system must first be constructed.An overview of this process is presented here and the interested reader is referred to [63] for further mathematical details.

Establish TEN Structure
The system processes and resources in an electrified transportation system are similar to what one might expect for a conventional transportation system.As highlighted in Table 2, the main difference is the introduction of charging processes that occur while a vehicle is being transported from one location to another (not necessarily distinct) location.Definition 4. Charging Process [50,77]: A resource-independent process p γ ∈ P γ that positively or negatively affects an EV's SOC.These processes may draw or inject the required energy into the interdependent electricity grid.
This definition of the system processes and resources in the electrified transportation system allows the definition of the associated knowledge base [63]: where: where J METS , J HETS and J γETS are smaller knowledge bases that reflect transformation, transportation, and charging functionality [63].The ⊗ and • are the Kronecker and Hadamard products respectively and the notation 1 n reflects a column vector of ones of length n [96][97][98].In effect, the electrified transportation system knowledge base is able to distinguish a collectively exhaustive and mutually exclusive set of processes out of the feasible combinations of processes and resources.For example, these can include "staying in place at a charging station b y1 while charging by wire" or "moving from b y1 to b y2 while charging wirelessly in online EV" [50,63].The construction of an electric power system knowledge base follows similarly.Table 2 presents the system's processes and resources [63].Consequently, the electric power system knowledge base is defined as [63]: where: where J MEPS , J HEPS and J γEPS are smaller knowledge bases that reflect transformation, transportation, and holding functionality [63].
The construction of the (full) TEN knowledge base J TEN then follows straightforwardly [63].
It recognizes that any resource in the electrified transportation knowledge base may potentially inject or withdraw power from the grid depending whether it is capable of either charging process (i.e., wirelessly charging or charging by wire).Also note that while electrified roads and rail lines appear as transportation resources in the electrified transportation system, they appear as transforming resources in the larger nexus system.The full mathematical derivation of the knowledge base(s) is available in [63].

Symmetrica Example
The Symmetrica test case [98] is now used as an example to illustrate how the TEN structure can be captured in the TEN knowledge base.As shown in Figure 2, it consists of interlinked electrified transportation system and electric grid topologies.The transportation topology consists of a 12 × 12 km grid with intersections at every kilometer.Five charging stations are placed in the city center at coordinates (4,4), (4,8), (8,4), (8,8), and (6,6).Additionally, the peripheral nodes represent home charging.Finally, 13 groups of 2 km × 2 km electrified road crosses are uniformly distributed across the area.The power grid topology is a modified version of the 201-bus distribution system test case [107,108].This test case will be referred to at each step of the assessment methodology.
Even for a moderately-sized test case such as Symmetrica, the system knowledge base becomes quite tedious to write manually and automated approaches instead are required.Here, it is important to recognize that the system knowledge base is a quantitative representation of an allocated architecture [98].Several systems engineering languages (e.g., UML/SysML) represent an allocated architecture graphically [109][110][111].Automated approaches can be used to transform the graphical representation into a quantitative one.
Consider Figure 3.For simplicity of discussion, four resources have been extracted from the Symmetrica test case topology.On the left, they are represented as class blocks in UML/SysML.Node #1 is a charging station.It has two methods (or system processes), park and charging, at Node #1; with "Null Charging" (p c1 ) or with "Charging by Wire" (p c3 ).Node #2 is a parking lot.Like the charging station, it has a parking process.Link #1 is an electrified road.It has a system process of Transport Vehicles from Node #1 to #2 with battery discharging (p c2 ).It also has the process "charge wirelessly" (p c4 ) during the transportation between nodes.Similarly, Link #2 has two system processes in the opposite direction.Figure 3 on the right shows the associated portion of the TEN knowledge base.There are 16 rows to account for the four possible transportation processes combined with four possible charging processes.There are 4 columns; one for each of the classes mentioned above.The matrix is then filled with ones where there are feasible combinations of system processes and resources (i.e., structural degrees of freedom) [96][97][98].This small subsection has 7 degrees of freedom.Naturally, the system knowledge base is highly sparse and must be computationally implemented as such [50].

Methodology
The TEN knowledge base serves as a structural skeleton of the TEN behavior.As expected, this behavior consists of the behavior related to the electrified transportation system and that of the electric power system.Furthermore, as these systems include both continuous-time as well as discrete-event dynamics, it requires a hybrid dynamic model [63].The full definition of the model is provided in Definition 5 below.It assumes sufficient background in graph theory [87][88][89], timed Petri nets [112][113][114], microscopic traffic simulation [83,84,115] and power system dynamics [47] which is otherwise obtained from the provided references.Full mathematical details of the development of the hybrid dynamic model is presented in [63].Its essential features are highlighted here.

Definition 5. TEN hybrid dynamic model
is the set of arcs represented as the difference of two incidence matrices.It represents the logical relationship from the events to the places and from the places to the events.• W : M → {0, 1} is the weighting function on the arcs.w(p y , t ωv ) = 1 if and only if is the timed Petri net discrete state vector for all discrete event times k.
• U is a binary vehicle firing matrix for all discrete event times k.• X = [X T ; X E ] is a continuous-time state vector representing the kinematic and electric state of the TEN.
It describes the continuous-time evolution of state vector X (Equations ( 10) and ( 13)) domain is a set of invariant conditions [112] which associates a discrete state Q to an interval of X and U within which X and U must remain in order to also remain in the discrete state Q.
evolves by the Petri-net transition function [63]: where: The state vectors Q B and Q E represent the number of vehicles at a Petri net place or undergoing a timed Petri net transition respectively [63].The incidence matrices M + and M − are derived from the TEN knowledge base, and the firing vectors U − k and U + k are derived from the vehicle firing matrix U [63].
The continuous time functions f = f T ∪ f E are classified into differential equation of the electrified transportation system f T and algebraic equations of electric power system f E [63].The differential equations of the transportation system f T are implemented in state space form [63]: They represent microscopic traffic dynamics of each vehicle traveling within the electrified transportation system [63].In this regard, both free traffic flow and car-follower models have been previously studied [83,116,117].f E (X E ) = 0 are the well-known power flow analysis equations [47]: where P is active power injection, Q is reactive power injection, V is bus voltage, and Y is the bus admittance matrix.Furthermore, the discrete state of the timed Petri-net discrete events is coupled to the active power injection in the electric power grid [63]: where α(ψ) would represent the charging rate per vehicle of charging station or electrified road [63].This hybrid dynamic mathematical model captures the discrete-event and continuous-time dynamic, describes the electrified transportation and electric power systems, and couples them together via the discrete state.

Symmetrica Example
The hybrid dynamic model takes on an intuitive character when applied to a test case such as Symmetrica.As in the previous section, a small portion of the Symmetrica transportation and power system topologies are taken for explanation.Consider Figure 4.It shows a Petri net with two places B ETS (shown as circles); one for each intersection in the transportation.It has seven discrete events (shown as rectangles) E ETS .These correspond to the seven degrees of freedom discussed in the previous subsection.The arcs between the places M reflect the logical sequence of discrete states a vehicle may take and can be directly calculated from the knowledge base [63].The weighting function W is restricted to one where arcs exist.The discrete state vector Q represents the number of vehicles at each Petri net place and transition.In Figure 4, three vehicles are shown with distinct colors red, green, and blue.The Petri net transition function would then describe how these vehicle markings would evolve from discrete event k to the next.In the case of Figure 4, the red and blue vehicles were "fired" through transitions T5 and T6 respectively to mean both were transported from Node #1 to Node #2.The first with inductive charging and the other without.Naturally, the binary vehicle firing matrix indicates which vehicle is "fired" through a given transition at a given discrete event index k.Between discrete events, the continuous state X = [X T ; X E ] evolves in time.The continuous electrified transportation system state X T describes as a minimum the position, velocity, and SOC of the vehicles going through the roads.It evolves by the continuous differential equation f T .The continuous electric power system state X E describes as a minimum the active and reactive power injections and voltages at the power system buses.It evolves by the continuous algebraic equation f E .Finally, it is important to recognize the continuous state X T adheres to a domain of applicability.Physically speaking, as vehicles reach the end of their respective road segments, or finish charging, the continuous state X T ceases to evolve, "pops" the vehicle out of its associated Petri net transition, and sends it to the downstream Petri net place.In such a way, the hybrid dynamic model describes the TEN behavior.This relatively complex behavior requires numerical simulation software to implement the mathematical model.Nevertheless, its primary advantage is that it supports an arbitrary TEN topology while leveraging established models of transportation and power system behavior.

Methodology
Once the behavior of the TEN has been established, the focus shifts to assessing the ITES decision-making in the TEN.Here, it is important to distinguish between designing new ITES decision-making algorithms and simply assessing the decisions that already exist.While the former has been identified as an extensive area of future work [49,50,76], this work simply seeks to address the latter so as to describe TEN performance without change.
The TEN hybrid dynamic model presented in Definition 5 provides a structured approach for such an assessment.As expected, the decision-making can be classified into continuous-time and discrete-event decision-making.At a faster time scale, when the continuous time dynamics in Equation ( 10) are enhanced with real-time vehicle controller signals U T which affect the vehicles' position, velocity, and SOC, they may be replaced with a new closed-loop real-time dynamic of the form: Such an approach may be useful in the case when studying the integration of new EV concepts.For example, EVs with regenerative breaking increase the SOC as the vehicle decelerates.In this regard, the assessment resembles the assessment of different vehicle concepts under varying drive cycle conditions [118].
At a slower time scale, the timed Petri net in the hybrid dynamic model of the physical TEN system takes the binary vehicle firing matrix U at each discrete event k as an input.Therefore, it may be viewed as the output of the ITES (as shown in Figure 1).Furthermore, as mentioned in the introduction, the ITES is assumed to take exogeneous use-case driven traffic demand in the form of itineraries U I that vehicles take over the course of the day [76].Figure 5 shows two such home-work-home itineraries.With this understanding, the ITES decision-making can be viewed as a function Φ ITES that transforms U I into the much richer vehicle firing matrices c.
The ITES decision-making may also reflect dynamic decisions that take into the discrete-event state of the vehicles.In such a case: The transportation-electrification assessment methodology can therefore be viewed as writing explicit mathematical decision-making models for Equations ( 14) and (15).While it is possible to imagine a single such mathematical model, the reality of transportation and electric power infrastructure is that several, perhaps entirely uncoordinated, decisions are taken.In that regard, the ITES decisions highlighted in Table 1 are now discussed in the context of both private and commercial uses cases.
Vehicle Dispatch: is defined here as when a vehicle should undertake a trip from origin to destination.As mentioned previously, this work assumes that the demand for trips between origins and destination is taken as input data and therefore acts as a decision constraint.Nevertheless, the timing and choice of vehicle can be decided.In the private use case, individuals may decide to expedite or postpone their departure so as to avoid congestion [119].In the taxi use case, the choice of which taxi is dispatched when and where while optimizing battery capacity can enhance vehicle fleet utilization and therefore operating revenue [120].This decision takes on a new dimension when considering that the demand for certain origin-destination trips can be shared within a single vehicle (e.g., Park & Ride) [121].
Route Choice: is defined as the set of roads and intersections that a vehicle should take along the course of an itinerary.The earliest route choice algorithms determined the shortest path with respect to distance for a given vehicle [122] and have since developed to include static as well as dynamic traffic assignment models [123][124][125].Later the algorithms evolved to minimize time; taking into account road speed limits (and not just length).Dynamic route choice algorithms consider the congestion on the roads as a state.For EVs where the range may be limited, shortest path routing requires to also take into account locations of charging stations [126][127][128].Regulations on high occupancy vehicles (HOV) can shift route choice behavior [37,38].
Charging Station Queue Management: is defined as when and where an EV should charge as queues develop in real-time.Because of their relatively lengthy charging time, EVs drivers can be informed with real-time state of charging queues.Instead of idly waiting to charge, drivers may change their choice of charging station.This choice is particularly important in the commercial use case after a given origin-destination trip has been completed.Charging station queues can also be managed in infrastructure planning.A location-sizing model to optimally allocate charging spots without exceeding waiting times has been proposed [129].
Coordinated Charging: is defined as when the vehicles at a given charging station should charge to meet drive departure times and power grid constraints.Typically, a charging station might work on a first-in-first-out heuristic.However, some vehicles in charging station parking lots may need to leave first.This presents a challenge when a given charging station is limited in the number of vehicles it can charge at a time.Finally, charging stations may wish to participate in demand response programs where dynamic prices can provide incentives to reshape charging load curves [9,[53][54][55][56][57].
Vehicle-2-Grid Stabilization: is defined as when EVs are used as energy storage to enhance the stability of the power grid.Such a decision can be viewed as a generalization of coordinated charging in that it allows EVs to both inject electric power as well as withdraw it from the grid.This may be achieved at a slower time scale so as to be part of price-based demand response program.It may also occur as a real-time feedback control loop that enhances either the local voltage or frequency at the point of connection [9,[53][54][55][56][57].

Symmetrica Example
The Symmetrica test case is meant to simulate a naive business-as-usual EV integration scenario.In recent years, there have been many EV pilots that integrate EVs without consideration for holistic TEN system performance and the need for enhanced ITES decision-making.As such, the ITES decisions in the Symmetrica test case are entirely uncoupled and reflect typical industrial practice.How each decision was implemented in simulation is now described sequentially.The traffic demand data found in [76] reflects the vehicle firing matrices U in the hybrid dynamic model and is the output from the sequential execution of these five decisions.
Vehicle Dispatch: As expected from Equation ( 14), the generation of the Symmetrica vehicle firing matrices U began with use-case driven traffic demand in the form of itineraries U I .For all vehicles, these represent a home-commute-work-commute-home use case, as a simplification of an average work day [76].Returning to the transportation system topology in Figure 1, vehicles enter from the peripheral nodes, go to five work locations coinciding with charging stations and then after 8 hours return back to the periphery [76].This input data has a spatial distribution so that all five work locations are equally important [76].Meanwhile, the timing of commute departures was exponentially distributed around peak morning and afternoon timescite [76].This input data was taken as exogeneous and no further vehicle dispatch decisions were made to change travel patterns.
Route Choice: While itinerary data specifies the location of the home-work-home use case, it does not specify the exact route taken and its associated firing vectors.Routing was performed by calculating shortest routes.Because there are many such routes, the routes were constrained so as to pass through the centers of the electrified road segments.Furthermore, vehicles having exactly the same home-work-home triplet were evenly distributed amongst the multiple paths that met these two criteria.In this way, no single road segment received an undue portion of the traffic.
Charging Station Queue Management: In keeping with the business-as-usual EV integration scenario, no effort was made to manage charging station queues.The sizes of the queues is effectively an emergent property.
Coordinated Charging: Similarly, the Symmetrica test case employed uncoordinated charging.As online EVs drive on electrified roads, they charge accordingly.Meanwhile, plugin-in EVs charge for a pre-calculated time as soon as they arrive at a charging station.The time spent charging was calculated based upon the expected value of lost charge during a commute.Each charging station has a capacity limit of simultaneously charging vehicles.This is 30 vehicles for the peripheral charging stations and 25 for the central charging stations.Together, these basic rules augment the vehicle route-itineraries to include charging functionality.This information is sufficient to generate a complete set of vehicle firing matrices.
Vehicle-2-Grid Stabilization: Similarly, no vehicle-2-grid stabilization decisions were made.Taken together, these five decisions provide all the necessary information to numerically simulate the TEN hybrid dynamic model.
To complete the simulation of ITES decision-making, energy management decisions for the electric power grid need to be included.For an uncoordinated charging scenario, the EVs present an aggregate charging power demand P EV .This demand, taken in concert with the overall electric power demand P D , must be integrated into a electric power system energy management market [47].The simplest of these is an economic dispatch that minimizes the production cost function C(P G ): subject to: P min Gi ≤ P Gi ≤ P max Gi (18) where P Gi power plant generation levels, C i are their associated cost levels, and P min Gi and P max Gi are the minimum and maximum generation levels respectively [47].The economic dispatch is performed for each discrete event time k, dispatching generation for each demand level P D + P EV .Naturally, such an energy management market can be enhanced to include unit commitment, ramping limits, and transmission limits and losses [47].In this uncoordinated charging scenario, the power system energy market acts "downstream" of EV charging and must make sure to dispatch sufficient electricity supply.

Methodology
Once the TEN dynamic behavior and the ITES decision-making have been established, the TEN's performance can assessed by numerical simulation.In that regard, a number of quantitative performance measures are required.As infrastructure systems, both the electrified transportation system and the electric power system have many stakeholders each with their requirements and performance measures.While an exhaustive identification of performance measures is intractable here, several technical, economic and environmental concerns are highlighted as an overview on the subject.
The first group of performance measures address the electrified transportation infrastructure.These include congestion, queues, quality of service, asset utilization, and EV availability.Definition 6. Road Congestion G road : The number of vehicles on the conventional and electrified roads at a given discrete time k: where C road is an appropriately sized constant row vector with ones corresponding to the elements associated with roads in Q E .

Definition 7.
Parking Congestion G parking : The number of parked vehicles at a given discrete time k: where C parking is an appropriately sized constant row vector with ones corresponding to the elements associated with parking lots in Q E .
Definition 8. Charging Congestion G road : The number of parked charging vehicles at a given discrete time k: where C charging is an appropriately sized constant row vector with ones corresponding to the elements associated with charging stations in Q E .
Definitions 6-8 may be normalized by the total number registered vehicles N V to give relative measures of electrified transportation demand.Definition 9. System Queues G queue : The number of queued vehicles at a given discrete time k: where 1 is an appropriately size column vectors of ones.
Definition 10.Quality of Service S quality : The fraction of time that the vehicle fleet is receiving a useful transportation service (i.e., motion, parking, charging).During these times, the vehicle fleet is not queued: where K is the total number of discrete events over the course of the day.
Definition 11.Resource (Asset) Utilization Υ R : Given a resource R, the fraction of a given resource's capacity that is being utilized averaged over time: where C E is a column vector containing the capacities of each of the structural degrees of freedom and C R is an appropriately sized row vector with ones corresponding to the elements associated with resource R.
Consequently, the asset utilization of the full electrified transportation system is: Definition 12. Vehicle Fleet Utilization Υ VF : The fraction of time that the vehicle fleet is being driven: Definition 13.Vehicle Fleet Availability A EV : The fraction of time that the vehicle fleet is neither charging nor queued: A EV = 1 − Total Charging Congestion and System Queues Total Time*Fleet Size (28) In commercial and sharing use cases, the vehicle fleet utilization is relatively high limiting the time for charging.At some point, the vehicle fleet is effectively utilized to a maximum level of 100%.Definition 14. Effective Vehicle Fleet Utilization A EV : The fraction of the available time that the vehicle fleet is being driven: The second group of performance measures address the electric power system.Traditionally, these include balancing performance, line congestion, and voltage security limits [47].CO 2 emissions is added as an environmental measure.
Power System Balancing Performance: Several existing power system standards introduce balancing performance measures [48].More recently, advancements in renewable energy integration studies suggest the standard deviation of imbalances normalized by peak load [130,131].These values depend greatly on the design of power system energy markets and the quantities of different types of operating reserves [132][133][134].This work assumes that the balancing performance is held constant as needed with a potential increase in the required operating reserves.
Power System Line Congestion: Power line rating places a physical limit on the amount of transferred active power P. Definition 15.Line Safety Criterion SC L [75]: Given a set of N l lines, a given line i may have a line limit P * i , over all discrete events K, the line safety criterion SC L is defined as the average amount of excess active power in all the lines: where: Power System Voltage Security The IEEE Standard 519 [135] places bus voltage limits to be within 0.95 and 1.05 volts per unit.Definition 16.Bus Safety Criterion SC B [75]: Over all discrete events K, the bus safety criterion SC B is defined as the average amount of insufficient or excess voltage v in all the buses b: where: Power System CO 2 Emissions: The additional CO 2 emissions ∆E CO 2 arises from their quantities before (E CO 2 ) and after (E EVCO 2 ) electrification: The emissions associated with either scenario are calculated from the economic dispatch in the power system market.Each power plant i has a carbon emissions intensity curve I i which is a function of the amount of power generated P Gi .Consequently: Naturally, the difference ∆E CO 2 is caused by the difference between the load curve P D (κ) and the net load curve P D (κ) + P EV (κ).
The third group of performance measures address the investment and operating costs as a result of coupling the transportation and electric power systems together.
TEN Investment Costs: From the previous discussion, the investment cost from electrified transportation IC total comes from three sources: the purchase of EVs, the installation of charging infrastructure, the expansion of electric power generation, and distribution: The purchase of EV is simply their number N EV times their per unit cost c EV : The installation of charging infrastructure is the sum of charging stations N CS and electrified roads N ER weighted by their per unit costs (c CS and c ER ): Finally, the cost of additional power generation depends on the newly installed capacity ∆C gen and the per unit cost of new generation c gen : The quantity of required additional capacity is determined by simulation from the power system's load duration curve before and after electrification [50].Mathematically: TEN Operating Costs: The additional operating costs from electrified transportation ∆OC total comes from two sources: the additional production costs in the energy markets, and the additional costs of maintaining power generation operating reserves: The additional production costs in the energy markets is the optimal value of the objective function (Equation ( 16)) before and after electrification: The additional costs of maintaining power generation operating reserves arises from their quantities before (Q res ) and after (Q EVres ) electrification:

Results: Symmetrica Test Case
The TEN performance of the Symmetrica test case is now discussed, for the scenario with 50% adoption of plug-in EVs.Two previous works [50,63] have already investigated this performance for two different scenarios, and the interested reader is encouraged to find more detailed discussion there.Here, the discussion first focuses on the performance measures described in the methodology and subsequently on the quantitative results.CO 2 emissions, investment costs, and operating costs are not presented here as this would require additional data regarding electricity demand and power plant characteristics, which is outside the scope of this methodology paper.

Performance Measures
From the simulation data, the performance measures mentioned in Definitions 10-14 are calculated, and summarized in Table 3.The resource utilization and EV fleet utilization are particularly low as only two, relatively short, trips are undertaken between 6:00 a.m. and 8:00 p.m., requiring charging at only two moments during the day; the rest of the day all charging capacity is unused.Availability of the EV fleet is restricted by the time required to charge each EV, but remains particularly high indicating more effective use of each vehicle would be possible as is also indicated by the low effective utilization.Increased utilization can be achieved by car sharing strategies.

Quantitative Results
Figure 6 shows the results of the transportation electrification Symmetrica case with 50% plug-in EVs.It includes six subplots spanning the electrified transportation and electric power systems, each corresponding to a quantitative performance measure presented previously.
Road Congestion: In Figure 6a, the number of driving vehicles on a given road segment is shown with a distinct color, corresponding to Definition 6.Its shape reflects the decision-making for Symmetrica, as the traffic is dispatched in an exponentially distributed fashion in both the morning and the afternoon commute.
Parking Congestion: In Figure 6b, the number of parked vehicles at any given station is shown, corresponding to Definition 7. It rises sharply as the early-morning commute arrives at their destination, and then continues to rise over the course of a few hours.
Charging Congestion: This slow rise in the late afternoon is due to charging station congestion as shown in Figure 6c.Here, the quantity of stationary charging vehicles is shown, corresponding to Definition 8.The charging station capacity is quickly reached during the morning commute resulting in a built-up of queues.
System Queues: The size of these queues at each charging station are shown in Figure 6d.Here, the quantity of queued vehicles is shown, corresponding to Definition 9. Given this queue, the quality of service will not be 100%.The initial queue is due to the way Symmetrica data is presented, representing queuing before entering the cities from the periphery.
Feeder Lines Active Power: In Figure 6e, the active power in the three feeder lines are shown, corresponding to lines originating from bus 201.Depending on the limits on each line P * i set in defining the line safety criterion, as introduced in Definition 15, these values may exceed limits.Here, the structure of the electrical grid and the charging station topology must be carefully rationalized so as to avoid underutilized and overutilized lines.In this particular case, the morning charging load is provided by just two of the three feeders but does not peak because of the saturation of charging stations.In the afternoon though, the charging activity peaks around 4:00 pm as all peripheral charging stations provide sufficient capacity, but is supplied almost exclusively by two feeder lines.Leveraging the topological flexibility of the electric distribution system may reduce reduce this peaking by dividing the total required power more evenly over all available feeder lines.
Terminal Bus Voltages: In Figure 6f, the bus voltage magnitudes of a selection of buses is shown upon completion of a power flow analysis at each discrete event k.The buses shown correspond to the terminals associated with the longest branches of the power system topology as shown in Figure 2. The safety criterion SC B , as defined in Definition 16, is met in this case, as all bus voltages shown are within 0.95 and 1.05 volts per unit.Nonetheless, the effects of EV charging on the bus voltages coincide with the peak in charging demand in the afternoon, but show a different effect in response to the morning charging load.The morning voltages remain relatively stable throughout most of the charging.During the afternoon commute on the other hand, bus voltages are depressed much more substantially to almost 0.98 pu.While this is partially caused by the magnitude of the charging load, it is particularly amplified by the placement of several charging stations on the same feeder branch near its terminal.This raises questions as to the most appropriate nature of a power distribution system topology.In conclusion, the ITES decision-making has a significant impact on the TEN's holistic performance.Charging coordination, vehicle dispatch, charging station capacity and topology affect both the electrified transportation and electric power systems.New decision-making techniques are required to manage these simultaneous objectives.

Conclusions
This paper provides a new hybrid dynamic system assessment methodology for multi-modal transportation-electrification.Fundamentally, EVs interact with three interconnected systems: the (physical) transportation system, the electric power grid, and their supporting information systems.Coupling of the two physical systems essentially forms a TEN.The assessment methodology, presented here, provides clear and descriptive ways to establish the TEN structure, establish the TEN behavior, establish the TEN decision-making, and assess TEN performance.The transportation system model describes microscopic discrete-time traffic operations so as to study the kinematic and electric state of each EV at all times.The electric power system model studies balancing line and voltage performance.In this work, the Symmetrica test case has been used as an illustrative example.The presented hybrid dynamic system model lends itself to the design of performance measures to quantitatively assess technical, economic, and environmental characteristics.This methodology may be used to assess the effects of transportation-electrification in any city or area, as it is developed to be as general as possible, opening up possibilities for many future studies.

Figure 1 .
Figure 1.A systems overview of the TEN with connected vehicles and ITES.

Figure 2 .
Figure 2. The Symmetrica TEN topology: a test case for transportation electrification sesearch [98]: (a) on left, the transportation system topology; (b) on right, the electric power system topology.

Figure 3 .
Figure 3. Construction of part of the Symmetrica knowledge base from UML-notated classes for unique nodes and edges (charging station nodes and electric roads).

Figure 4 .
Figure 4.A small subsection from the Symmetrica transportation topology is taken, and illustrated into a Petri-net.The Petri-net tracks the movement of three unique cars, represented by three colored tokens, through the network.

Figure 6 .
Figure 6.Transportation Electrification Symmetrica test case results with 50% plug-in EV adoption.(a) Number of driving vehicles in time; (b) Number of parked vehicles in time; (c) Number of stationary charging vehicles in time; (d) Number of queued vehicles; (e) Active power through feeders (MW) in time; and (f) Bus voltages (p.u.) in time.
domain) where: • B ETS is the set of timed Petri net places.It represents transportation independent buffers (e.g., stations & intersections).• E ETS is the set of timed Petri net discrete events.It represents the structural degrees of freedom (as defined in the previous section).•

Table 3 .
Performance results for Symmetrica in 50% plug-in EVs scenario.