Capacity-Constrained Contraflow Adaption for Lane Reconfiguration in Evacuation Planning

This paper presents a heuristic contraflow-based reconfiguration evacuation algorithm, which is named Capacity-Constrained Contraflow Adaption (CC-Adap). First, it effectively calculates optimal candidate routes for evacuation. Second, an evaluation method is proposed for estimating these candidate routes. Third, CC-Adap utilizes a contraflow-based method to reconfigure the evacuation routes to improve capacity constraints. Fourth, traffic conditions are updated in real time. Fifth, CC-Adap reuses historical evacuation routes to reduce the computational cost and accelerate the evacuation process. Experimental results show that CC-Adap generates high-performing evacuation strategies and can be used to tackle large-scale evacuation planning.


Introduction
In recent years, both natural and man-made disasters have posed serious threats to humans, such as Hurricane Andrew [1], Hurricane Katrina [2], the "9/11" terrorist attack event [3], and 2011 Tohoku Earthquake Tsunami [4].Effective evacuation plans are necessary for saving lives and minimizing casualties.The most important objective of an effective evacuation plan is to transfer evacuees to safe areas as quickly as possible.However, dynamic factors make this objective complicated to achieve.For example, the original intention in designing a transportation network is not to address the sharp increase of traffic flow when an emergency occurs [5].The number of evacuees might far exceed the road capacity, thereby making evacuation route planning a computational challenge.
A well-known routing algorithm, namely, Capacity-Constrained Route Planner (CCRP) [1], was proposed for tackling evacuation route planning.CCRP models the capacity constraints of roads and intersections and selects the shortest paths for evacuation.Later, two heuristics, namely, Intelligent Load Reduction (ILR) and Incremental Data Structure (IDS), were introduced to improve the algorithm's performance [6].Although routing algorithms can effectively solve mid-size traffic evacuation problems [6], they might not be effective for large-scale evacuation scenarios with massive number of evacuees, which might lead to heavy traffic congestion, owing to capacity constraints.Contraflow-based methods are considered effective techniques for alleviating traffic congestion [2].Contraflow-based methods reconfigure the edges in the ideal direction and reallocate edge capacity to reduce the evacuation time.For example, Kim et al. presented a heuristic approach that combines a Greedy method and a contraflow-based method for implementing transportation network reconfiguration and, to some extent, minimizing the total evacuation time [2].However, finding an optimal contraflow plan for reconfiguring network is computationally challenging because we must enumerate all possible combinations of edges in the exhaustive search space [2].Moreover, it takes a considerable amount of time to evaluate these contraflow candidates by calculating the evacuation time.
In this paper, we propose a Capacity-Constrained Contraflow Adaption (CC-Adap) algorithm for effective evacuation.CC-Adap consists of five steps.First, we utilize an evacuation algorithm to generate optimal candidate routes for evacuees, which finds the routes with maximal flow rate among all the available routes from source vertices to sink vertices.Second, we propose an evaluation method for heuristically evaluating the candidate routes' capabilities based on the traffic conditions and determining the appropriate routes for evacuation.Then, we implement a contraflow-based method for optimizing the routes' performances and reducing traffic congestion.Next, traffic conditions are updated in real time to make the route planning more practical.Finally, CC-Adap reuses the historical evacuation routes before calculating the available routes at each new time step.Experimental results indicate that CC-Adap boosts performance in evacuation route planning.
The remainder of this paper is organized as follows: Section 2 introduces some related works.Section 3 defines the transportation network and contraflow problem.Section 4 details the CC-Adap evacuation algorithm.Experimental results and analysis are presented in Section 5. Section 6 presents the conclusions of this paper.

Related Work
To the best of our knowledge, the existing methods can be divided into descriptive and prescriptive approaches [5].Descriptive approaches aim at mimicking real emergency evacuations based on the simulation of the traffic situations and drivers' behaviours, while prescriptive methods mainly focus on providing a high-performing evacuation plan with minimum evacuation time.

Descriptive Methods.
The aim of descriptive methods is to simulate traffic evacuation situations as vividly as possible.Many traffic evacuation simulation tools have been proven to be capable of solving complex traffic evacuation problems, including MITSIMLab [7], VISSIM [8], MATSim [9], DynusT and DYNASMART [10], and AIMSUN [11].In addition, some simulation methods utilize multimodel integration methods to simulate evacuees' behaviours during evacuation and combine various existing evacuation methods for more realistic traffic simulation.For example, Pan et al. implemented a simulation framework for classifying evacuees' behaviours into locomotion, steering, and social for evacuation analysis [12].Wu and Huang combined a control volume model and a flow merging hypothesis to simulate evacuees' behaviours in high-rise building evacuation situations [13].Noh et al. considered different behaviours of heterogeneous evacuees in building evacuation [14].They presented a partially dedicated evacuation strategy that divides evacuees into heterogeneous groups, assigns corresponding routes based on a flow model and a simulation optimization approach, and, to some extent, minimizes the average evacuation time [14].Beloglazov et al. studied the behaviours of people in wildfire evacuation and integrated traffic simulator, wildfire simulator, and behaviour simulator to define the risk metric and provided detailed plan of how evacuation unfold [15].Yuan et al. constructed a multilevel agent decision model for simulating driver's behaviours and determining actions for each agent [16].These descriptive methods can achieve satisfactory performance in traffic evacuation for local communities.However, most exhibit poor generalization performance owing to high computational complexity.

Prescriptive Methods.
Prescriptive methods typically provide suggested schedules for evacuation planners, aiming at reducing traffic congestion and minimizing total evacuation time.Zeng and Wang made a small modification to CCRP to improve the performance, which gives priority to longer evacuation routes for evacuating evacuees [17].Kang et al. calculated the Dijkstra shortest evacuation paths and chose the best one for evacuees based on current traffic conditions [18].Shahabi and Wilson proposed a Capacity-Aware Shortest Path Evacuation Routing (CASPER) framework for emergency traffic evacuation, which utilizes an advanced traffic model and an intelligent routing algorithm to calculate the optimal evacuation plan for evacuees [5].Khan et al. used the Intelligent Transportation System (ITS) to compute the maximum traffic flow and the routes with least traffic congestion towards safe places for evacuees based on realtime traffic conditions [19].Pourrahmani et al. considered the uncertain evacuee demand at the pick-up source points by utilizing fuzzy credibility theory to present a genetic algorithm (GA) for handling the stated evacuation problem [20].Chen et al. implemented a distribution Load-balancing Emergency Guiding System (LEGS) to assign the fastest paths to transfer evacuees to the exits according to the capacity constraints and concurrent movement of evacuees [21].Ikeda and Inoue proposed a Multi-Objective Genetic Algorithm (MOGA) to evaluate evacuation routes by considering route distance and evacuation time to select optimal evacuation routes [22].Heydar et al. divided the transportation network into pedestrian network and vehicular network, where pedestrians were assigned the shortest routes to the safe areas in pedestrian network and vehicles transferred evacuees to designated shelters in vehicular network [23].These prescriptive methods have demonstrated the feasibility of calculating evacuation strategies for traffic evacuation.However, the performances of these prescriptive methods are usually limited owing to the capacity constraints of road segments and intersections.
Previous works [24][25][26][27] have considered the benefits of contraflow-based methods on evacuation planning.In addition, in the literatures, the optimality of evacuation planning by integrating contraflow-based methods has been validated.For example, Xie et al. designed a bilevel framework for optimizing network evacuation performance, in which the lower level solves the traffic assignment problem and the upper level integrates a contraflow method with crossing elimination strategies to optimize the traffic assignment [28].Wang et al. implemented a multiobjective optimization model to determine the evacuation priorities and setup time for doing the contraflow operations [3].Even et al. implemented a conflict-based path-generation algorithm to simultaneously guide the evacuation and select the to-bereversed roads [29].Kim et al. presented a Bottleneck Relief (BR) approach that iteratively flips min-cut edges based on the theorem in [30] for solving massive traffic evacuation problems [2].
Above-mentioned contraflow-based methods are selected to generate contraflow strategy that can tackle capacity constraints effectively [2,3,[24][25][26][27][28][29].However, there are two difficulties in calculating the appropriate contraflow strategy for the complex traffic system.The first difficulty is that it is a time-consuming task to determine the appropriate to-be-reversed edges in the complex traffic network [2,24,25].The second difficulty is that it is difficult to determine the ideal direction of these edges after being reversed [2,3,28,29].Thus, it is vital for the planners to find an appropriate contraflow strategy that could overcome these two difficulties.
So, this paper presents a prescriptive algorithm to generate appropriate contraflow strategy and optimal evacuation plan.For the first difficulty, the presented algorithm utilizes Greedy method and iterative optimization technique to select the to-be-reversed edges on the evacuation routes.For the second difficulty, the ideal direction of these to-be-reversed edges is determined by the traffic flow on the evacuation routes.

Network Modeling.
Suppose we are given a multisource and multisink transportation network (, ), where  = {  ,   ,   } and , respectively, represent the sets of vertices and edges.Here,   is the source vertex set,   is the transition vertex set, and   is the sink vertex set.The definitions of the vertices and edges are given as follows: (i) Each vertex   ∈  has a maximal vertex capacity  , and each source vertex   ∈   has an extra initial occupancy  , , which represents the number of evacuees.For each sink vertex   ∈   , the value of  , can be set according to the actual demand.In this paper,  , of each sink vertex is set to infinity.
(ii) Each edge  ⟨  ,  ⟩ ∈  has a maximal edge capacity  ⟨  ,  ⟩  , a travel time  ⟨  ,  ⟩  , and an initial direction. ⟨  ,  ⟩  is the number of evacuees (e.g., residents or vehicles) that can pass per unit period.The edge direction is initialized according to the real road conditions. ⟨  ,  ⟩  is the time cost when evacuees pass edge  ⟨  ,  ⟩ .Various factors should be considered when formulating edge travel time, such as road length, traffic flow, and drivers' behaviours.Many works have studied the relationship between these factors and edge travel time [13][14][15]17].This paper focuses on the impact of road length and traffic flow on the edge travel time.Considering the emergency situations, the traffic flow on the edges should be limited by the maximal edge capacity.Thus, the evacuees can use maximal speed passing one given edge when the traffic flow does not exceed the maximal edge capacity.In other words, more evacuees can be evacuated to safe area using this strategy in the minimum time.According to this, the edge travel time in the transportation network model is determined by the length and the maximal edge capacity of one given road.
There are different traffic elements in transportation network including buildings, crossroads, playgrounds, shelters, and parks.Buildings, crossroads, and playgrounds are the most common traffic elements in transportation network.Figure 1 illustrates the modeling results of the buildings, the crossroad, and the playground.Firstly, the buildings are modeled as source vertex.Secondly, the crossroad is modeled as transition vertex.Thirdly, the playground is modeled as sink vertex.The edges of the vertices are modeled according to the traffic configuration.As for the other traffic elements, they are modeled similar to these three common traffic elements shown in Figure 1.
The objective of transportation evacuation is to find an optimal plan with minimum evacuation time.The objective function is defined in V is the total evacuation time,  is the initial time step when a group of evacuees leaves the source vertex,  is the to-be-solved evacuation route set, and  is an evacuation route in .
() is used to calculate the travel time of route  when evacuees pass through the edges along route  and is expressed in (2).For route where   is the th vertex of route ,  0 is the source vertex of route , and   is the sink vertex of route . denotes the available residual capacity of edge  ⟨ −1 ,  ⟩ on route  at its arrival time  −1 .
Formula (1) illustrates the capacity constraints along the time dimension during evacuation.Constraint A ensures that the outflow of source vertex   does not exceed the initial occupancy of source vertex   .Constraint B enforces that the inflow of the sink vertices must be equal to the initial occupancy of the source vertices.In addition, the model capacity of vertices and edges cannot be negative during evacuation due to Constraint C.

Contraflow Problem.
In a large transportation network, long distance traffic jams usually cover the roads when emergencies occur.Contraflow-based methods are considered effective techniques for remedying traffic congestion [2].However, an appropriate contraflow reconfiguration is difficult to calculate.Generally, two main difficulties are encountered in calculating optimal contraflow reconfiguration: (1) it is NP-hard to find all possible combinations of edge directions and calculate the evacuation time of these contraflow combinations and (2) it is time-consuming to evaluate the effects of dynamic traffic flow on these contraflow reconfiguration candidates.
Figure 2 illustrates an evacuation scenario and two contraflow-based evacuation planning results.The evacuation time is 8 time units based on the original network configuration in Figure 2(a).Figures 2(b) and 2(c) illustrate two contraflow reconfigurations, which combine all the twoway edge directions and merge the edge capacities to increase the upper bound on the evacuation capacity.The contraflow reconfiguration in Figure 2(b) reduces the evacuation time to 5 time units, while the contraflow reconfiguration in Figure 2(c) reduces the evacuation time to 7 time units.The only difference between the two contraflow reconfigurations is the edge direction between vertex  2 and vertex  3 .Thus, it is vital for an optimal contraflow reconfiguration to select the critical edges that affect evacuation performance and determine their ideal direction after reversal.

Proposed Algorithm
This paper proposes the CC-Adap algorithm, which heuristically evaluates the evacuation routes, and integrates it with a contraflow-based method for distributing the evacuation routes with the least traffic congestion and maximal flow rate to evacuees.In addition, CC-Adap reuses historical evacuation routes to reduce computational costs and accelerate the evacuation process.Algorithm 1 gives the pseudocode of the CC-Adap algorithm.
Step 0 (initialization).For a given network (, ), let  be the to-be-reversed edge set,  be the initial time step end if (18) end if (19) end if (20) end while Algorithm 1: CC-Adap algorithm.
when evacuees leave source vertices,  be the to-be-solved evacuation planning route set, and  be the historical evacuation route set.
Step 1 (GetEvaRoute).CC-Adap utilizes the evacuation route generation method in the Max-Flow Rate Priority (MFRP) algorithm [31] to select optimal candidate route  from among the available routes from source vertices to the nearest sink vertices.MFRP introduces parameter  to describe the route capacity.The formula of  is shown in (4).
The main strategy of the MFRP algorithm [31] is to calculate the available candidate routes by Dijkstra algorithm and balance the number of transferred evacuees with the travel time to select evacuation routes with maximal .The performance of MFRP has been analyzed in [31].CC-Adap guarantees the shortest possible evacuation time by employing MFRP [31] to select the optimal candidate routes.
Step 2 (RouteEvaluation).The use of low-quality routes might aggravate traffic congestion and increase the total evacuation time.CC-Adap introduces a parameter  for quantizing route , which is generated in Step 1, and historical evacuation route to avoid selecting lower-quality evacuation routes.The formula of  is presented in where   is the historical evacuation route in , which has the same source vertex as .Equation ( 5) represents the evacuation capability gap between routes  and   when they evacuate evacuees at same times under the same traffic conditions.
The RouteEvaluation function operates as follows: First, we introduce parameter  for reserving the evaluation result and set  = −1 initially.For each historical evacuation route   in , we evaluate  between  and   by employing (5), where  and   have the same source vertex, and reserve the evaluation result using  as follows, which is expressed in (6): If  is less than zero, then we set  = 0 and calculate the next historical evacuation route   in .Otherwise,  is suitable for evacuation and we set  = 1 and terminate the evaluation procedure.If there is no route with the same source vertex as  in , we assume that  can also be used for evacuation and set  = 1.
Step 3 (ContraflowReconfiguration).Greedy method can be used in evacuation route planning [2].The appropriate contraflow strategy can be determined using Greedy method which has the advantage of obtaining the global optimal solution by local optimal solution.Thus, CC-Adap integrates Greedy method and iterative optimization technique to reverse the two-way edges on the evacuation routes iteratively in this paper.However, the edges on the evacuation routes cannot be reversed repeatedly as long as the edges have been flipped.The modified edge direction of the to-be-reversed edges depends on the traffic flow along the evacuation routes.The upper bound on the capacity of the flipped edges is increased by merging the two-way edges.The travel time of the reversed edges is equal to that of the two-way edges, which have the same edge direction as the flipped edges.Figure 3 illustrates the steps for implementing the proposed contraflow reconfiguration method.Figure 3(a) is the original network configuration.In Figure 3(b), we select evacuation route ⟨ 1 −  2 −  4 −  6 ⟩.Its traffic flow is min{10, 3, 10, 4, 7, 2, ∞} = 2.As shown in Figure 3(c), the two-way edges on this route are merged and the maximal capacity of the flipped edges is increased.After contraflow reconfiguration, the traffic flow is increased to min{10, 7, 10, 7, 7, 5, ∞} = 5.Thus, the evacuation capability of route ⟨ 1 −  2 −  4 −  6 ⟩ is improved.
Step 4 (NetworkUpdate).After evacuees are transferred from source vertex to sink vertex, the available residual capacities of transition vertices and edges along route  are updated in (7) and (8), respectively, and the source vertex's residual evacuees of route  are updated in (9).
Step 5 (HRouteReuse).When there is no available evacuation route at the current time step, CC-Adap proceeds to the next time step (++) and returns to Step 1 to calculate routes.Before calculating routes, CC-Adap reuses historical evacuation routes in  to reduce the computational cost and accelerate the evacuation process.To avoid evacuating along the same evacuation routes continually, CC-Adap introduces a parameter  that marks historical evacuation route ℎ in  if ℎ is reserved at the last time step. is defined in The HRouteReuse function operates as follows: For each historical evacuation route ℎ in , if (ℎ) = 0 under the current traffic conditions, it will be excluded from .Otherwise, if the corresponding  is equal to 0, ℎ is reserved at the last time step and is not suitable for evacuation at the current time step, we set  = 1 for the next iteration.If  is equal to 1, then ℎ is used for evacuation at the current time step and traffic conditions are updated by Step 4.
In each iteration, the generated route  in Step 1 is reserved in  and  when  is used to transfer evacuees and  of  is set to 0. Meanwhile, the to-be-reversed edges on route  are reserved in .To complete the evacuation, Steps 1-5 can be repeated until there is no evacuee who needs to be evacuated.

Experimental Design.
There are many evacuation algorithms that can provide evacuation route planning, such as the classical CCRP algorithm [1], LRP algorithm [17], MR algorithm [31], MFRP algorithm [31], and Greedy algorithm [2].Reference [31] has validated that the performance of MFRP algorithm is better than LRP and MR.Besides, Greedy algorithm is a classical contraflow evacuation algorithm.CC-Adap algorithm integrates MFRP algorithm and Greedy algorithm in this paper.Thus, CCRP, MFRP, and Greedy algorithms are selected to validate the performance of CC-Adap algorithm.

Experiment Data.
In this paper, we use the following notations to describe evacuation situations:  is the number of vertices,  is the number of source and sink vertices,  is the number of edges, and  is the number of evacuees.The experimental data that are used in this paper is generated by the network simulator tool NETGEN [32].However, NETGEN does not generate the maximal vertex capacity for transition vertices.This paper uses a random number generation function, namely, rand(), to generate the maximal vertex capacity for transition vertices in (11).The seed of the rand() function is set to the vertex number to avoid generating the same networks. and  are two coefficients in (11), which imply that the lower bound on the edge capacity is  and the upper bound on the edge capacity is (()%(( * )/)) + . and  can be reset according to the practical requirements.In this paper, we set  = 3 and  = 20 [31].

Experimental Result
5.3.1.Feasibility Verification.The purpose of this section is to evaluate the feasibility of using CC-Adap in evacuation planning.CC-Adap is compared with CCRP and MFRP to determine whether CC-Adap can be used for evacuation planning and to optimize evacuation plans by applying a contraflow-based method.
In this section, two experiments are carried out.Each experiment has five test groups and the number of vertices, edges, source, and sink vertices is fixed.We vary the number of evacuees from 20000 to 100000.Figures 4 and 5 depict the results on the evacuation time and run time of the three algorithms in small-scale and large-scale transportation networks, respectively.According to Figure 4(a), CC-Adap has the smallest and nearly linearly growing evacuation time with the increase of the number of evacuees, while CCRP and MFRP have steeply increasing evacuation time and the evacuation time of CCRP is almost double that of CC-Adap.The run time of CCRP and MFRP grows rapidly as the number of evacuees increases, while CC-Adap has a stable run time and the run time of CCRP and MFRP is more than twice that of CC-Adap when the number of evacuees is 100000 in Figure 4(b).As shown in Figure 5, CC-Adap has the shortest evacuation time with the increase in the number of evacuees in Figure 5(a), followed by CCRP, and MFRP has the longest evacuation time.In terms of the run time, CC-Adap, CCRP, and MFRP all show nearly linear growth with the increase in the number of evacuees, and CC-Adap has a run time that is nearly half those of CCRP and MFRP.

Scalability
Verification.The purpose of this section is to evaluate the scalability of CC-Adap to complex evacuation situations and CC-Adap is compared with CCRP, MFRP, and Greedy algorithms.

Comparison on Networks with Different Number of Edges.
Figure 6 plots the experimental results in terms of evacuation time and run time of the four algorithms against the number of edges.This experiment is carried out with five subexperiments.The number of vertices is fixed at 256, and    time with the increase in the number of evacuees, comparing with CCRP, MFRP, and Greedy (which have growth ratios of nearly 50 percent).

Summary of Experiments.
There are many dynamic factors that are critical to the transportation network, such as number of evacuees and affected areas, which makes evacuation management a daunting task.These dynamic factors can be evaluated systematically by building appropriate evacuation models [33].Only by taking these factors into consideration we can optimize evacuation management.This section will discuss these dynamic factors and evaluate the optimality of the CC-Adap algorithm based on the experimental results.
(i) Number of edges: the more edges there are, the more difficult it is to find appropriate routes to evacuate evacuees.CC-Adap can intelligently evaluate the routes to choose the optimal one for evacuation, thereby making it more scalable with the number of edges.
(ii) Number of source and sink vertices: an increase in the number of source and sink vertices in a transportation network means that the evacuees are      more scattered in the network [1].It might lead to difficulty in calculating the evacuation routes.CC-Adap reuses the historical evacuation routes to reduce the computational costs, thereby making CC-Adap more scalable with the number of source and sink vertices.
(iii) Number of evacuees: capacity constraints are the major limiting factor in optimizing an evacuation strategy.Traffic congestion occurs easily when there are too many evacuees.CC-Adap applies a contraflow-based method to improve the performance of evacuation routes, which makes CC-Adap more scalable with the number of evacuees.
(iv) Traffic scale: large and complex urban evacuation problems cause most evacuation route planning approaches to fail.CC-Adap integrates an iterative optimization technique and a contraflow-based method to optimize an evacuation plan, thereby making CC-Adap more effective in addressing urban evacuation planning.

Conclusions
This paper presents a heuristic contraflow reconfiguration evacuation algorithm, which is called Capacity-Constrained Contraflow Adaption algorithm.CC-Adap can heuristically evaluate evacuation routes to assign high-quality routes to evacuees.Moreover, the proposed contraflow-based method is applied to reconfigure evacuation routes to reduce traffic congestion and improve evacuation performance.Meanwhile, CC-Adap reuses historical evacuation routes to reduce computational costs and accelerate the evacuation process.
Experimental results have shown that the CC-Adap algorithm can optimize evacuation route planning in terms of evacuation time and run time.
Journal of Advanced Transportation 13 However, more research is required for CC-Adap.First, a more effective evacuation route generation method should be constructed.Second, the evaluation method should be improved to evaluate the evacuation routes more effectively.

Figure 1 :Figure 2 :
Figure 1: Illustration of the vertices and edges.

Figure 3 :
Figure 3: Example of the proposed contraflow-based method.

Figure 4 :
Figure 4: Evacuation quality in a small transportation network.

Figure 5 :
Figure 5: Evacuation quality in a large transportation network.

Figure 6 :
Figure 6: Evacuation quality with respect to the number of edges.

Figure 8
plots the experimental results on the evacuation time and run time of the four algorithms against the number of source and sink vertices.This experiment is carried out with five subexperiments.The number of vertices is fixed at 256, the number of edges is fixed at 1528, and the number of evacuees is varied from 20000 to 100000.Each subexperiment has two test groups and the number of source and sink vertices is varied from 26 to 52.According to Figure8, CC-Adap produces the solutions with shortest evacuation time comparing with the CCRP, MFRP, and Greedy algorithms when the number of source and sink vertices is varied.In addition, the evacuation time of CC-Adap has the slowest increase with the increase in the number of source and sink vertices comparing with the other three algorithms.For different number of source and sink vertices, CC-Adap has the fastest run time comparing with the CCRP, MFRP, and Greedy algorithms.As the number of source and sink vertices varies from 26 to 52, the run time of the MFRP and Greedy algorithms rises steeply, CCRP shows a small growth in run time, and CC-Adap shows the smallest increase in run time.Comparison on Large Transportation Networks.Figure9depicts the experimental results of the four algorithms in terms of evacuation time and run time on a large-scale transportation network.The experiment is carried out with five test groups and has a fixed network configuration of 1024 vertices, 5258 edges, and 102 source and sink vertices.We vary the number of evacuees from 20000 to 100000 to test the algorithms' performances.As shown in Figure9(a), CC-Adap can generate high-quality evacuation solutions with the shortest evacuation time when the number of evacuees grows in a large-scale transportation network.Moreover, the run time of CC-Adap shown in Figure9(b) also has the smallest and most stable growth with different number of evacuees in large-scale transportation network evacuation scenarios.In addition, CC-Adap has the minimum growth ratio of run (a)  = 768;  = 256;  = 26 (b)  = 768;  = 256;  = 26 (c)  = 1528;  = 256;  = 26 (d)  = 1528;  = 256;  = 26 (e)  = 1528;  = 256;  = 52 (f)  = 1528;  = 256;  = 52

Figure 7 :
Figure 7: Evacuation quality with respect to the number of evacuees.

Figure 8 :
Figure 8: Evacuation quality with respect to the number of source and sink vertices.

Figure 9 :
Figure 9: Evacuation quality of large-scale transportation network evacuation.