NetworkPartitioningAlgorithmsforSolvingtheTrafficAssignment .

3y ago
18 Views
2 Downloads
616.44 KB
18 Pages
Last View : 1m ago
Last Download : 3m ago
Upload by : Camryn Boren
Transcription

rk Partitioning Algorithms for Solving the Traffic Assignment ProblemUsing a Decomposition ApproachCesar N. YahiaGraduate Research AssistantDepartment of Civil, Architectural and Environmental EngineeringThe University of Texas at Austin301 E. Dean Keaton St. Stop C1761Austin, TX 78712-1172Ph: 779-804-4407, FAX: 512-475-8744Email: cesaryahia@utexas.edu(Corresponding author.)Venktesh PandeyGraduate Research AssistantDepartment of Civil, Architectural and Environmental EngineeringThe University of Texas at Austin301 E. Dean Keaton St. Stop C1761Austin, TX 78712-1172Ph: 737-222-8473, FAX: 512-475-8744Email: venktesh@utexas.eduStephen D. BoylesAssociate ProfessorDepartment of Civil, Architectural and Environmental EngineeringThe University of Texas at Austin301 E. Dean Keaton St. Stop C1761Austin, TX 78712-1172Ph: 512-471-3548, FAX: 512-475-8744Email: sboyles@mail.utexas.edu32Word count:4448 words text 188 words abstract 580 words references 5 tables/figures 250 words (each) 6466 words33Submission date: August 1, 201728293031TRB 2018 Annual MeetingPaper revised from original submittal.

1ABSTRACT2Recent methods in the literature to parallelize the traffic assignment problem consider partitioninga network into subnetworks to reduce the computation time for large-scale networks. In this article, we seek a partitioning method that generates subnetworks minimizing the computation timeof a decomposition approach for solving the traffic assignment (DSTAP). We aim to minimize thenumber of boundary nodes, the inter-flow between subnetworks, and the computation time whenthe traffic assignment problem is solved in parallel on the subnetworks. We test two different methods for partitioning. The first is an agglomerative clustering algorithm heuristic that decomposesa network with the objectives of minimizing subnetwork boundary nodes. The second is a flowweighted spectral clustering algorithm that uses the normalized graph Laplacian to partition thenetwork.345678910111213141516171819202122We assess the performance of both algorithms on different test networks. The results indicate that the agglomerative heuristic generates subnetworks with a low number of boundary nodes,which reduces the per iteration computation time of DSTAP. However, the partitions generated maybe heavily imbalanced leading to significantly higher computation time when the subnetworks aresolved in parallel separately at a certain DSTAP iteration. For the Austin network partitioned into4 subnetworks, the agglomerative heuristic requires 3.5 times more computational time to solvethe subnetworks in parallel. We also show that the spectral partitioning method is superior in termsof minimizing the inter-flow between subnetworks. This leads to a faster convergence rate of theDSTAP algorithm.Keywords— Network partitioning; Spectral partitioning; DSTAP; Regional modeling; Decentralized traffic assignment; Network contractionTRB 2018 Annual MeetingPaper revised from original submittal.

Yahia, Pandey, and Boyles1112The traffic assignment problem (TAP) is used to predict route choice and link flows for a giventravel demand. The static version of this problem can be formulated as a convex program (1) andsolved efficiently using modern specialized algorithms (2, 3). However, there are computationallydemanding problems that require solving TAP multiple times or on a large scale. Those problemsinclude bi-level mathematical programs with equilibrium constraints (4), statewide or nationalmodeling of the transportation network, and Monte Carlo simulations to account for forecastingerrors of TAP parameters 930INTRODUCTIONMethods for parallelizing the traffic assignment problem to decrease computation time havebeen studied recently e.g. (5, 6, 7). The decomposition approach to the static traffic assignmentproblem (DSTAP) was developed to decrease computation time by solving partitions of the transportation network in parallel (5). This approach creates subproblems for each partition and a masterproblem that equilibrates traffic across subnetworks. The master problem also includes regionaltraffic that has an origin or a destination outside a certain subnetwork or in two different subnetworks. To find equilibrium in this master-subproblem framework, the DSTAP algorithm exploitsthe equilibrium sensitivity analysis method developed in Boyles (8) to generate artificial links thatrepresent all paths between certain network nodes. DSTAP is shown to converge to equilibrium fora general network and its computation time is stated to depend on the subnetwork partitions (5).The objective of this study is to identify and test partitioning algorithms that can improvethe performance of a decomposition approach for solving the static traffic assignment problem.We seek to generate partitions that minimize the number of boundary nodes and the inter-flowbetween subnetworks. These requirements minimize the interactions between subnetworks whichinfluences the time needed to converge to a global equilibrium in a framework such as DSTAP. Inaddition, we seek partitions that minimize the computation time needed to solve the traffic assignment problem separately and in parallel for the subnetworks. This refers to the per iteration lowerlevel subproblems in DSTAP. With this motivation and objective in place, we test two algorithmsin our analysis. The first algorithm is the one proposed in Johnson et al. (12) for the objectivessimilar to those required in this paper. The second algorithm is based on flow weighted spectralclustering. We compare the performance of the algorithms on real-world networks against thestated objectives.35The remainder of this paper describes methods to parallelize TAP and evaluates the performance of partitioning algorithms when applied to decomposition methods for solving TAP. Section2 reviews current methods for solving TAP and partitioning networks. Section 3 presents the algorithms evaluated and their use in the DSTAP framework. Section 4 presents demonstrations fordifferent transportation networks. Section 5 concludes the paper.3623739This section summarizes existing literature in the following areas: the latest advancements in thetraffic assignment problem methods, a parallelization approach to the traffic assignment problem,and the need for efficient network partitioning algorithms.40Algorithms for solving TAP are generally classified as either link-based, path-based, or3132333438LITERATURE REVIEWTRB 2018 Annual MeetingPaper revised from original submittal.

Yahia, Pandey, and 7282930313233343536373839404142432bush-based. Link based methods work in the space of link flows and require less operationalmemory than path-based methods, but are much slower to converge (13, 14). Bush-based methodexploit the acyclic nature of paths that are used by origins at the equilibrium (2, 3, 15, 16). Recentwork also includes -optimal improved methods for solving TAP on large problems (17). Althoughthe recent advancements have improved the state-of-the-art for solving TAP fast and efficiently,there is still a need for faster methods. These are several instances which require TAP to be solvedon large scale networks or solved repeatedly. These include running large-scale statewide modelsfor evaluating multiple projects, solving TAP multiple times for network design problems withequilibrium constraints (4, 18), to name a few.To address the computational demands of large scale or iterative traffic assignment problems, methods that aim to parallelize TAP have been developed. Bar-Gera (6) describes a parallelization approach based on the paired-alternative segments. The algorithms proposed in Chenand Meyer (19) and Lotito (20) also parallelize TAP by decentralizing the computations for eachOD pair. A decomposition approach for static traffic assignment (DSTAP) developed by Jafari etal. parallelizes TAP by network geography instead of the traditional decomposition approach byOD pairs (5). This article aims to identify partitioning algorithms that minimize the computationtime per iteration of DSTAP and the total computation time required to reach convergence. Proofsof convergence and correctness of the algorithm are provided in Jafari et al. (5).The literature on network partitioning algorithms is extensive. These algorithms can bebroadly classified into agglomerative/divisive heuristics, integer programming based approaches,and spectral partitioning algorithms. Integer programming formulations for the partitioning problem are proven to be NP-hard (9) and approximation heuristics have been proposed (9, 21).Heuristics for generating partitions based on agglomerative and divisive clustering havebeen recently used in various transportation related applications. Saedmanesh and Geroliminis(11) used an agglomerative clustering heuristic for generating partitions based on “snake” similarities for applications of the macroscopic fundamental diagram. Etemadnia et al. (9) developedsimilar heuristics for distributed traffic management where a greedy agglomerative strategy wasadopted for combining adjacent nodes with high flow. Johnson et al. (12) developed anotherheuristic for decentralized traffic management. This heuristic aims to minimize boundary nodes insubnetworks and to create subnetworks of similar size. Their heuristic performed better than theMETIS algorithm proposed in (22).Spectral partitioning is an alternative approach for clustering and partitioning graph e.g.(23, 24, 25, 26, 27). Bell (10) applied a capacity-weighted form of the spectral partitioning methodsto investigate vulnerability. Other transportation applications include air traffic control and urbantraffic signal control systems (28, 29). The partitioning mechanism is based on the eigenvaluesassociated with the graph Laplacian. The partitions that result from the spectral partitioning havelow inter-cluster similarity (24). Additionally, using the normalized Laplacian generates graphsthat are balanced by the weight chosen. This is an important feature since ignoring the balancerequirement results in cuts that isolate a small number of peripheral nodes. For example, theminimum cut program that aims to minimize the weight between resulting partitions will oftenresult in separating one node from the rest of the network (26). However, incorporating balancerequirements causes the cut problems to become NP-hard. Spectral partitioning is an approximatemethod for obtaining a cut with minimal inter-flow while satisfying balance requirements (23, 26,TRB 2018 Annual MeetingPaper revised from original submittal.

Yahia, Pandey, and Boyles3127).23312We consider a directed network G defined by a set of nodes N and set of edges A. Let M bethe node-node adjacency matrix for the network. M is an N N matrix, with elements mijequal to 1 if there is a link connecting node i to j and zero otherwise. We also define the weightedadjacency matrix MGD with elements mG,D(i,j) equal to w(i,j) if (i, j) A and zero otherwise, wherew(i,j) is the weight assigned to link (i, j) A. In this article, we assume w(i,j) to be the flowon link (i, j) at the equilibrium. To construct a graph Laplacian, we use an undirected version ofMGD , denoted by MG , defined as the sum of MGD and its transpose. The elements of MG are mG(i,j) .The graph diagonal matrix DG is defined as a diagonal matrix with principal diagonal elements inrow i as the sum of elements in row i of MG : dii j mG(i,j) . The graph Laplacian is defined asLG D G M G .133.114We aim to partition a large scale network into subnetworks such that an algorithm based on thedecomposition approach to the static traffic assignment problem (DSTAP) is solved efficiently. Inorder to properly define the objectives of the partitioning algorithms evaluated, we review the mainelements of the DSTAP algorithm developed by Jafari et al. 33343536373839NETWORK PARTITIONING FOR DECENTRALIZED TRAFFIC ASSIGNMENTDecomposition Approach to the Static Traffic Assignment Problem (DSTAP)DSTAP is an iterative aggregation-disaggregation algorithm consisting of two levels, a master problem and a set of lower level subproblems corresponding to the respective subnetworks. Asubproblem corresponds to solving the traffic assignment problem for a specific subnetwork. Themaster problem is used to model interactions between the subproblems. In the master problem, thesubnetworks are aggregated using first order approximation methods based on equilibrium sensitivity analysis (8, 30). This results in artificial links representing the subnetworks in the masterlevel problem. The algorithm proceeds by solving the subproblems in parallel, aggregating thesubnetworks using artificial links, shifting flow towards equilibrium in the simplified master levelnetwork, obtaining subnetwork boundary flow from the master level iteration, and then proceeding to disaggregate the flow on subnetworks and solving the subproblems in parallel again. Thisprocedure is repeated until convergence to a global equilibrium as shown in Figure 1.The computational performance of DSTAP at each iteration depends on the number of artificial links. These links need to be updated at each iteration using equilibrium sensitivity analysisto incorporate the latest information on travel costs. To reduce the number of artificial links generated, the number of boundary nodes associated with the subnetworks needs to be minimized. Inaddition to the regional artificial links that approximate subnetworks at the master level, there aresubnetwork artificial links generated for each subproblem to represent flow that originates froma subnetwork then traverses other subnetworks before returning to the subnetwork. Therefore, toreduce subnetwork artificial links the flow traversing multiple subnetworks needs to be minimized.The computational performance at each iteration is also influenced by the time needed tosolve the traffic assignment problem in parallel for the subnetworks. This represents solving theK lower level subproblems in Figure 1. The computation time needed to solve the subproblems inTRB 2018 Annual MeetingPaper revised from original submittal.

Yahia, Pandey, and Boyles1234567891011121314151617181920214parallel is dominated by the subproblem that requires the greatest computational cost. Therefore,to reduce this computation time, the subproblems need to balanced in size. This can be achievedby balancing the flow distribution across subnetworks instead of having few subnetworks with themajority of flow.Consider the maximum excess cost termination criteria defined as the greatest differencebetween the longest used path and the shortest path for each OD pair, it was shown in Jafari etal. (5) that the maximum excess cost for the full network OD is bounded by the total numberof boundary points across subnetworks B̃ multiplied by the sum of the maximum excess cost forthe master level regional network rOD and the maximum excess cost for all subnetworks sOD asshown in Equation 1. Therefore, to reach convergence faster, we need to increase the rate at whichthe bound in Equation 1 tightens. The subproblem maximum excess cost sOD can be reduced bysolving the subproblems to a low gap level. After approximating the subnetworks with artificiallinks, the master level maximum excess cost rOD could be obtained. We note that if the interflow between subnetworks is minimized, then the artificial links representing the subnetworks willhave a similar cost structure across successive iterations since the influence of external flows onsubnetwork equilibrium is reduced. Therefore, the least cost path in the master level regionalnetwork would be relatively invariant across iterations. This implies that rOD could be reduced ata higher rate. In the extreme case where the master level least cost path is completely dominatedby constant costs on artificial links, the maximum excess cost could be reduced to zero by placingall the regional flow on the path with the least cost artificial link. Thus, faster convergence couldbe reached by minimizing the inter-flow between subnetworks. OD 2B̃( rOD sOD )22(1)233.22425We test the performance of two algorithms that aim to partition the network such that the computational time for a decomposition approach to solve traffic assignment is minimized.263.2.127The first heuristic algorithm used is the shortest domain decomposition algorithm (SDDA) proposed in Johnson et al. (12). This algorithm works in an agglomerative fashion and constructsa given number of partitions of a network such that the number of boundary nodes between thesubnetworks is minimized (primary objective) and the partitions are balanced in size (secondaryobjective). Achieving the algorithm objectives would minimize the computation time per iterationassociated with DSTAP. Minimizing the boundary nodes reduces the number of artificial links created, and generating balanced subnetworks would reduce the time needed to solve the subproblemsin parallel. In addition, this algorithm only depends on the topological properties of the graph. Thisfeature is desirable when limited information is available on link costs, flow between OD pairs, orother data that could form the basis of a partitioning algorithm.2829303132333435363738Partitioning AlgorithmsDomain decomposition algorithmThe sequential steps of the algorithm are shown in Algorithm 1. The algorithm constructsthe partitions by identifying source nodes which are “far” from each other given a distance mea-TRB 2018 Annual MeetingPaper revised from original submittal.

Yahia, Pandey, and Boyles5FIGURE 1 Algorithm for the decomposition approach to the static traffic assignment problem (5)TRB 2018 Annual MeetingPaper revised from original submittal.

Yahia, Pandey, and Boyles123456sure. We use the number of links on a breadth-first search tree between two nodes as the distancemeasure between the nodes instead of the P-LCA algorithm with unit costs proposed by Johnsonet al. (12). This distance measure indicates the extent of separation of two nodes and is used todetermine association of a node to the source nodes of the partitions. The reader is referred toJohnson et al. (12) for more information on this algorithm.Algorithm 1 Shortest domain decomposition algorithm (12)Step 1: InitializationLet ns be the number of subnetworks/partitions to be generatedSet Rsn : MAXStep 2: Determining first source nodeSet the rank of each node as the sum of the number of incoming and outgoing linksChoose the node with lowest rank s1 as the first source nodeStep 3: Updating the rank and determining other source nodesfor i in 2 : ns doPerform breadth-first search from every source node, sj 1 j iDetermine the rank of node n as a (i 1)-dimensional vector whose elements are thedistance of node n from source nodes sj where 1 j iChoose the node which has the highest total rank (sum of all elements in the rank vector)Resolve ties in favor of nodes which have minimum value of the sum of pair-wise difference between each element of the rank vectorAssign the chosen node as the i-th source node siStep 4: Populate subdomain associated with each source nodeFor each node, assign it to the source node to which it has the minimum distanceStep 5: Identify system boundary nodes and allocate the subnetworksfor (i, j) A doif i and j are assigned to different source nodes thenAdd i and j to the set of boundary nodes.Stop63.2.27Spectral graph theory is used to study graph properties using the Laplacian matrix of the graph.This theory can be used to identify cuts through the graph which have low cost by computing theeigenvalues and eigenvectors of the Laplacian matrix. The cost of a cut is defined as a ratio of theweights on cut links to the size of the smaller subnetwork separated by the cut (24, 26).89101112131415Spectral partitioningThe eigen

38 traffic assignment problem methods, a parallelization approach to the traffic assignment problem, 39 and the need for efficient network partitioning algorithms. 40 Algorithms for solving TAP are generally classified as either link-based, path-based, or TRB 2018 Annual Meeting Paper revised from original submittal.

Related Documents:

2.1 Anatomi Telinga 2.1.1 Telinga Luar Telinga luar terdiri dari daun telinga dan kanalis auditorius eksternus. Daun telinga tersusun dari kulit dan tulang rawan elastin. Kanalis auditorius externus berbentuk huruf s, dengan tulang rawan pada sepertiga bagian luar dan tulang pada dua pertiga bagian dalam. Pada sepertiga bagian luar kanalis auditorius terdapat folikel rambut, kelenjar sebasea .

Conducting a Comprehensive Skin Assessment Presented by Dr. Karen Zulkowski, D.N.S., RN. Montana State University

EA BA EWITA (Business Architecture Enterprise-wide IT Architecture) Recent work (e.g., Paul Harmon at Cutter, META Group, and a number of federal agency EA projects) emphasizes the importance of including BA in the EA definition.

Carlo DOMENICONI (b. 1947) Concerto Mediterraneo, Op. 67 Oyun Trilogy Toccata in blue Chaconne Concerto Mediterraneo, Op. 67 for two guitars and large orchestra The longing to live in a world in which nature is respected and its beauty recognised as vital to human life prompted Carlo Domeniconi to compose a concerto for two guitars and orchestra. His experience of witnessing the .

OCD in Children and Teens The information contained within this pack was correct at the time of sharing. We update this on a regular basis. If you notice any links are broken or information has changed please contact ShropshireFIS@shropshire.gov.uk and we will update the information. Further Family Information Services and Resource Packs are available through the Early Help website www .

Entry requirement for the CIPS Level 4 Diploma are 2 A Levels or 2 years business experience. Additional entry requirements for the Apprenticeship are Maths and English GCSE at A-C. (or their equivalents). Candidates without this may sit a literacy and/or numeracy test prior to beginning the course. 3.8 Taster sessions Weekly taster sessions from 6.30 – 7.30 at the Russell Square Study .

the risk of a new Cold War and are elaborated in Chapter 4. Recommendation 1: Reduce strategic tensions Avoiding a Cold War will require the US and China to strengthen, not reduce, the many areas of cooperation that once bound them, dampen down their hostile rhetoric and get serious about conducting a whole-

ways to fight Communism, the evil-s of Communism, and the false doctrines of Communism, it could become quite dangerous. Richard I. Miller, Associate Director of the National Education Association Project on Instruction, in re-porting on Communism to the Council of Chief State School Officers at Miami Beach said: