Multi-Objective Evaluation of Airborne Self-Separation Procedure in Flow Corridors Based on TOPSIS and Entropy

This paper proposes a simulation-based framework for assessing airborne self-separation procedures in flow corridors with consideration of different performance metrics, including air traffic operations, corridor capacity, safety, and environmental impacts. Firstly, the airborne self-separation concept in flow corridors is introduced, followed by an agent-based flow corridor simulation model. Then, data were collected to initialize a parallel-lane flow corridor model connecting A461 upper air route from Beijing to Guangzhou in China which can also simulate aircraft self-separating in the flow corridor. The total control delay, flow corridor throughput, breakout rate, and the CO2 emissions of traffic flow were considered as the impact measurements, and the TOPSIS and entropy method was used to rank the performances of different self-separation procedures. We found that combining multiple objectives into one, the optimum scheme can be obtained to guide the design of self-separation procedures for flow corridors. The research results can be used by airspace managers to dynamically develop appropriate operational procedures and rules for flow corridors given different operational conditions and constraints. Also, the framework proposed in the research may be used to evaluate the design of airspace structure with consideration of multiple objectives.


Introduction
Flow corridors concept is a long tube-shaped new airspace structure introduced within Next Generation Air Transportation System (NextGen). It is designed for reducing current air traffic complexity and restructuring the airspace to provide more capacity. The essence of flow corridors refers to bundles of near-parallel four-dimensional trajectory (4DT) assignments, which consequently achieve a very high traffic throughput, while allowing traffic to shift as necessary to enable more effective weather avoidance and reduce congestion [1,2]. Flow corridors integrate many similar prototype concepts such as high-volume tube-shaped sectors (HTS), highways-in-the-sky, tubes, and dynamic multi-track airways [3][4][5][6], three common and prominent characteristics of these concepts that would distinguish them from today's airways are the multiple-parallel lanes design, self-separation operations, and dynamic activation rules [7]. With the anticipated flight demand increase in the future, flow corridors concept has become an important component of future air transportation systems.
Currently, the flow-corridors-related literature is mainly focused on the development of concept of operations, such as the geographical allocation, required equipage and capabilities, and

Description of Airborne Self-Separation in Flow Corridors
The airborne self-separation in flow corridors refers to the transfer of the responsibility of maintaining separation with other aircraft from air traffic controllers to the pilots of each aircraft in flow corridors. Different from the usual airborne self-separation concept for resolving complex encounter scenarios in crossing waypoints [24], the airborne self-separation in flow corridors mainly refers to aircrafts flying within multiple near-parallel lanes with appropriate separations and shifting over to other lanes for overtaking slower-moving aircrafts or resolving potential conflicts.
In [6], NASA proposed three airborne self-separation procedures for track configurations with spacing and passing capabilities that are aimed to enable aircrafts to fly more consistently at their optimal cruise speed. One of the promising procedures that permits aircrafts to switch and stay at the available lane is selected for multi-objective evaluation. In this procedure, aircrafts flying in flow corridors are much like the ground traffic moving in the highways. If an aircraft enters the flow corridor, it is responsible for self-separated with its leading aircraft by itself all the time. For safety and efficiency consideration, it may adjust its acceleration, velocity, and switch lane for overtaking and/or avoiding loss of the separation. Generally speaking, it is a procedure with more flexibility and efficiency, but it also comes with some potential risks if we cannot deal well with the self-separation procedures. Thus, we take this procedure as the starting point for our research.

Simulation Model of Self-Separation in Flow Corridors
To perform multi-objective evaluation for airborne self-separation procedures in flow corridors, simulation modeling is the best method which can capture various aircraft behaviors. An agent-based approach was proposed to develop a flow corridor simulation model with C++ language [16,17], in which each aircraft in the model is represented by an agent with both operational performance states and self-separation states. In this model, the flow corridor is designed with two closely parallel-lane structures located at the altitude of 35,000 ft. Aircrafts are assumed to fly in the same direction within different lanes in the flow corridor and will self-separate with others in the same lane by adjusting velocity and/or switching lanes. Some high-level description of the simulation model is illustrated as follows.
The aircraft dynamic model is adapted from Glover and Lygeros' research [25], which includes three control variables and six key operational state variables. Aircraft thrust (T), bank angle (ϕ), and angle of attack (α) are set as control variables in along-track, across-track, and vertical dimensions, separately. The along-track position (x), true airspeed (v), across-track position (y), heading (ψ), altitude (z), and flight-path angle (γ) are six key state variables used to simulate aircraft flying in the flow corridor.
For capturing different self-separation behaviors, five self-separation states were also defined in the model, which are the target velocity flying state (TVF), velocity adjusting state (VA), lane-changing state (LCS), breakout state (BS), and locked state (LS). The VA and TVF state refer to states in which an aircraft is flying with varied speed or constant speed for self-separation from lead aircraft. LCS and BS refer to states that an aircraft switching lane or breakout from the flow corridor for efficient or safety consideration. The LS is an auxiliary state works with VA and TVF states to prevent simultaneous lane changes or breakouts. Some aircraft self-separation behavior examples within flow corridors are shown in Figure 1. and efficiency consideration, it may adjust its acceleration, velocity, and switch lane for overtaking and/or avoiding loss of the separation. Generally speaking, it is a procedure with more flexibility and efficiency, but it also comes with some potential risks if we cannot deal well with the self-separation procedures. Thus, we take this procedure as the starting point for our research.

Simulation Model of Self-Separation in Flow Corridors
To perform multi-objective evaluation for airborne self-separation procedures in flow corridors, simulation modeling is the best method which can capture various aircraft behaviors. An agent-based approach was proposed to develop a flow corridor simulation model with C++ language [16,17], in which each aircraft in the model is represented by an agent with both operational performance states and self-separation states. In this model, the flow corridor is designed with two closely parallel-lane structures located at the altitude of 35,000 ft. Aircrafts are assumed to fly in the same direction within different lanes in the flow corridor and will self-separate with others in the same lane by adjusting velocity and/or switching lanes. Some high-level description of the simulation model is illustrated as follows.
The aircraft dynamic model is adapted from Glover and Lygeros' research [25], which includes three control variables and six key operational state variables. Aircraft thrust (T), bank angle (φ), and angle of attack (α) are set as control variables in along-track, across-track, and vertical dimensions, separately. The along-track position (x), true airspeed (v), across-track position (y), heading (ψ), altitude (z), and flight-path angle (γ) are six key state variables used to simulate aircraft flying in the flow corridor.
For capturing different self-separation behaviors, five self-separation states were also defined in the model, which are the target velocity flying state (TVF), velocity adjusting state (VA), lanechanging state (LCS), breakout state (BS), and locked state (LS). The VA and TVF state refer to states in which an aircraft is flying with varied speed or constant speed for self-separation from lead aircraft. LCS and BS refer to states that an aircraft switching lane or breakout from the flow corridor for efficient or safety consideration. The LS is an auxiliary state works with VA and TVF states to prevent simultaneous lane changes or breakouts. Some aircraft self-separation behavior examples within flow corridors are shown in Figure 1.  The key self-separation variables may influence the performance of airborne self-separation procedures include the initial traffic density, minimum separation, target velocity, target separation, lane-switch buffer, and velocity difference threshold. The initial traffic density refers to the interval arrival of aircraft in the flow corridor; the minimum separation is a distance-based interval standard between adjacent aircraft that no separation is allowed to be smaller than. Target velocity and target separation are two variables used for controlling the operational speed and density of traffic flows in the flow corridor. Lane-switch buffer and velocity difference threshold are important variables used for trigging flights to switch states. The main simulation procedures are as follows: (1) Build up the simulation framework and initialize parameters, including construction of the flow corridor model, aircraft model, atmosphere model, and initialization of simulation parameters.
Parameters for flow corridor model include the number of lanes, corridor length, width, and altitude. Important parameters used for aircraft modeling are aircraft type, reference mass, cruising Mach number, flight envelope, engine thrust and aerodynamics coefficients, etc. Atmosphere model will provide wind speeds, the standard pressure, temperature, density, and the speed of sound on the altitude of flow corridor. Simulation parameters include the number of simulated flights, replication times, time-step, etc. The proposed simulation model and procedures were checked and visually verified by Google Earth. However, the model does not consider the unknown failures of the equipment. For more details of the agent-based flow corridor simulation model, refer to [17].

Data and Methods
This section presents data and methods used for multi-objective evaluation of airborne self-separation procedures in the paper. Specifically, field data were collected first, followed by the introduction of impact measurements. Then, a multi-objective evaluation method, comprising TOPSIS and entropy, was used to compare the performances of different self-separation procedures.

Data Collection and Preprocess
To construct a practical flow corridor model for evaluation, field data were collected based on Chinese historical flight operational data in 2017. According the primary statistic results, the A461 upper route in China was selected as the basis for flow corridor modeling in this research. A461 is one of the busiest routes in China; it serves more than 200 flights per day and 486 million passengers in total. Firstly, we created the route structure model based on the Aeronautical Information Publication (AIP) data. Currently, the A461 segment connecting Beijing to Guangzhou is about 970 nmi, including 21 important navigation/fix points, as shown in Figure 2a. We designed a parallel-lane flow corridor for this high-altitude route using the great circle trajectory between RENOB and ATAGA points, which are set as the entry and exit points separately. The configuration of the model is supposed to be using RNAV Q-routes with two closely parallel lanes at the FL350 (Flight Level 35,000 ft). The flow corridor is 940 nmi long in total with two 8 nmi width near-parallel lanes, as shown in Figure 2b. Firstly, we created the route structure model based on the Aeronautical Information Publication (AIP) data. Currently, the A461 segment connecting Beijing to Guangzhou is about 970 nmi, including 21 important navigation/fix points, as shown in Figure 2a. We designed a parallel-lane flow corridor for this high-altitude route using the great circle trajectory between RENOB and ATAGA points, which are set as the entry and exit points separately. The configuration of the model is supposed to be using RNAV Q-routes with two closely parallel lanes at the FL350 (Flight Level 35,000 ft). The flow corridor is 940 nmi long in total with two 8 nmi width near-parallel lanes, as shown in Figure 2b. Then, the flights flying from ZBAA (Beijing Capital International Airport), ZBNY (Beijing Nanyuan Airport), and ZBTJ (Tianjin Binhai International Airport) to ZGGG (Guangzhou Baiyun International Airport), ZGSZ (Shenzhen Bao'an International Airport), and ZGSD (Zhuhai Jinwan Airport) between 26 March and 28 October 28 2017 were selected for traffic-flow-features extraction. After preprocessing, 20,434 flights records were available for data analysis; important data fields include the departure airport, arrival airport, aircraft type, planned departure/arrival time, actual departure/arrival time, route segments, passing time, etc.
After preliminary analysis, we can see that A320, A332, A333, A388, B737, B772, B77W, and B787 are eight primary aircraft types, as shown in Table 1. The maximum amount of flights is B737, which comprises about 43% of the flights. A320 and A333 are two other aircraft types that are frequently used, and they account for 19% and 16%, separately. The mean and standard deviation of flying time for all flights in A461 were calculated and collected, and they were later used for model verified and validation. The average flying time ranges from about 122 to 131 min. Then, the flights flying from ZBAA (Beijing Capital International Airport), ZBNY (Beijing Nanyuan Airport), and ZBTJ (Tianjin Binhai International Airport) to ZGGG (Guangzhou Baiyun International Airport), ZGSZ (Shenzhen Bao'an International Airport), and ZGSD (Zhuhai Jinwan Airport) between 26 March and 28 October 28 2017 were selected for traffic-flow-features extraction. After preprocessing, 20,434 flights records were available for data analysis; important data fields include the departure airport, arrival airport, aircraft type, planned departure/arrival time, actual departure/arrival time, route segments, passing time, etc.
After preliminary analysis, we can see that A320, A332, A333, A388, B737, B772, B77W, and B787 are eight primary aircraft types, as shown in Table 1. The maximum amount of flights is B737, which comprises about 43% of the flights. A320 and A333 are two other aircraft types that are frequently used, and they account for 19% and 16%, separately. The mean and standard deviation of flying time for all flights in A461 were calculated and collected, and they were later used for model verified and validation. The average flying time ranges from about 122 to 131 min. The statistics of aircraft types and their ratios are used as inputs for the proposed simulation model, the performance parameters of aircraft, including aircraft mass, Mach number, reference wing surface area, weight, thrust specific fuel flow, etc. are obtained from BADA (User Manual for the Base of Aircraft Data) database [26]. Another important input parameter is the interval arrival separation for the aircrafts in the flow corridor. We did a survey of the number of arrival flights per hour for A461 route, and the exponential distribution with mean of 81.76 nmi is fitted as the simulated function for random arrivals in our experiments. The initial values of other important parameters and variables are shown in Table 2.

Traffic Operational Metrics Used for Evaluation
The air traffic operational impact, corridor capacity, safety, and environmental impacts of different target separation standards were estimated. For a particular self-separation procedure, the total control delay was used to measure the operational impact of aircraft in the corridor. The corridor capacity was measured by using the average aircraft throughput per hour. The safety was measured by using the corridor breakout rate, while the aircraft CO 2 emissions were considered as the environmental impact measure.

Impacts on Traffic Operation and Corridor Capacity
The total control delay was used to measure the operational impact of aircraft in the parallel-lane corridor. The control delay refers to the time difference between the nominal travel time and the practical travel time that results from different self-separation procedures. Note that different type of aircraft has different nominal travel time, which equals to the length of the corridor divided by the aircraft standard cruise velocity provided by BADA at the corridor altitude. Given the simulation information from the flow corridor model, the operational impacts associated with various self-separation procedures for aircraft were estimated by standard Monte Carlo simulation. The average control delay was measured as follows: where W d [min] is the average control delay for the aircraft that passed the corridor; K is the simulation replications for the particular scenario. N is the number of the simulated aircraft for each replication; B j is the set of breakout aircraft and |B j | represents the breakout events happened in the j-th replication; L represents the flow corridor length; V CR,i [kt] represents the standard cruise velocity for aircraft i at the corridor altitude; t sim,i,j [min] represents the actual corridor travel time of aircraft i in the j-th replication. Similarly, the corridor capacity was measured using total passed aircraft divided by the total simulation time. The average capacity for the corridor was estimated as: where C [aircraft/h] represents the average capacity for the corridor; Γ j [h] represents the total simulation time in the j-th replication.

Impacts on Traffic Safety
The safety is measured by the breakout rate which represents the percentage that an aircraft has to break out from the flow corridor to avoid flying with the separation less than the minimum standard. In the simulation, the separations between adjacent aircrafts are calculated at each time step. With the updated performance data, each flight agent checks the potential loss of separation and resolves the conflict by changing its acceleration.
However, if the conflict cannot be resolved in time, and the current separation with the lead flight is already smaller than the minimum separation, it has to exit the corridor. Therefore, the safety metric may also refer to the lane-switch rate or potential conflict rate (sum of breakout rate and lane-switch rate). However, in the collision risk assessment, along with the impending near midair collision (NMAC), the breakout rate could further be used to estimate the overall rate of a collision rate [24], as shown in Equation (3). We decide to use the breakout rate as the metric.
Collision rate = breakout rate × Pr (NMAC | breakout rate) ×Pr (collision | NMAC, breakout rate) In addition, this paper only focuses on the potential non-nominal self-separation events that occur in the corridor. We will not further estimate the potential collision rate. The average breakout rate P r for a particular target separation standard was estimated as follows:

Impacts on Environment
The environmental performance was measured by using aircraft emissions, which are typically associated with the aircraft engines. The greenhouse gases are the primarily considered emissions at the higher altitudes for corridor (FL350), and the aircraft engines emissions are made up about 70% of CO 2 , a little less than 30% H 2 O, and less than 1% other compounds [27]. CO 2 is considered as the environmental impact measure in the paper. Aircraft engine emissions are measured by emission indexes which have units of grams per kilogram of fuel consumed. The emission index of CO 2 EI CO2 used is 3155 grams/kilogram fuel [26]. In order to estimate aircraft CO 2 emissions in the corridor, the fuel consumption was estimated first. For the jet engines, according to the Operations Performance File (OPF) of BADA, the thrust specific fuel consumption η [kg/(min·kN)] is specified as a function of where C f1 [kg/(min·kN)] is the 1st thrust specific fuel consumption coefficient; C f2 [knots] is the 2nd thrust specific fuel consumption coefficient. The corridor simulation model could calculate the true airspeed and thrust for all the aircraft at any time interval in the corridor. The average CO 2 emissions, E CO2 [kg], were estimated as a function of the CO 2 emission index EI CO2 , cruise fuel consumption and cruise time, as shown in Equation (7). Here, t represents the discrete time period; ∆t is the time-step and T is the last time period; f cr,i,t,j is the fuel flow of aircraft i in cruise phase at the simulation time t in the j-th replication.

Multi-Objective Evaluation Compared with the Ideal Solution
To evaluate the performance of different treatments of self-separation procedures, an ideal solution and a negative ideal solution are created as basis for comparison. The ideal solution has the best performance for each index among all alternatives, whereas the negative ideal solution has the worst performance for each index [28][29][30]. The relative closeness between each alternative and the ideal solution is calculated. The alternative with the shortest geometric distance comparing with the ideal solution, or with the largest geometric distance comparing with the negative ideal solution is regarded to have the best performance. Therefore, different alternatives can then be ranked. A key process of the multi-objective evaluation process is to determine the relative importance of different indexes, which is known as the weight. In this paper, the entropy method is used to determine the weight. It is regarded that the performance measure with greater dispersion may contain more information, thus have larger effects on the multi-objective evaluation results [31,32]. The larger the variation, the greater the weight is.

•
Step 1: Generate the evaluation matrix. Suppose that there are m scenarios to be evaluated (m equals ten for each case). Let a ij indicate the performance evaluation value for the jth performance measure for the ith scenario (i = 1, 2, . . . , m, j = 1, 2, . . . , n); then, the evaluation matrix is illustrated as A = [a ij ] m × n .

•
Step 2: Normalize the performance measures. As the performance measures are in different units, i.e., the throughput is measured in aircraft/h, the potential conflict rate is measured in probability, the average delay is measured in minutes, and the average fuel consumption is measured in ton, these values are required to be normalized by transforming them into a dimensionless value. The standardization equation is written as: where r ij represents the jth performance value for the ith scenario after normalization, 0 ≤ r ij ≤ 1. The evaluation matrix after normalization is written as R = [r ij ] m × n .

•
Step 3: Calculate the entropy of different indexes. Let E j represent the entropy for the jth performance measure. The following equation can be used to calculate E j : where 0 ≤ E j ≤ 1; k is an adjustment factor, which can be estimated as 1/ln m.

•
Step 4: Determine the entropy weight ω j for the jth performance measure. • Step 5: Calculate the geometric distance between each alternative and the ideal alternative. Let R + = {r 1 + , r 2 + , . . . , r n + } represent the decision matrix for the ideal solution and . . , r n − } represent the decision matrix for the negative ideal solution, where r j + = max{r ij |I = 1, 2, . . . , m }, r j − = min{r ij |I = 1, 2, . . . , m }, j = 1, 2, . . . , n. We can calculate the distance from each performance measure to D i + and D i − as follows: • Step 6: Estimate the relative closeness G i for the ith treatment to the ideal solution.
The ith treatment can be closer to the ideal solution with the increase of the value of G i .

Numerical Test
This section presents a case study for the established Beijing-Guangzhou flow corridor, which is arguably the busiest city-pair in China. According to the simulation model and the proposed TOPSIS and entropy methods, the comprehensive performance of airborne self-separation procedures were evaluated.

Results of Self-Separation Simulation in Flow Corridors
The simulation model was initialized by using Chinese domestic flights information of the A461 segment between Beijing and Guangzhou, and it was visually validated by the Google Earth software. To capture the performances of different traffic flow within flow corridors, firstly we simulated the airborne self-separation procedures with different traffic volume. For example, we increased the initial traffic density from 0.01 to 0.2 aircraft/nmi/lane based on the interval arrival of A461, and assuming that the minimum separation is set as 5 nmi, the target separation is set as 7 nmi, the lane-switch buffer is 1 nmi, the target velocity is 0.78 Mach, the distance threshold is 10 nmi, and the velocity difference threshold is 40 nmi. Figure 3 shows the simulation results. In general, with the increase of traffic density, the breakout rate curve shows a slight decrease and sharp increase trend, while the throughput, average delay, and CO 2 emissions curves show some increase and decrease trends. Similarly, Figure 4 illustrates the metrics-changing trends with the increase of target separation. It is worth noting that the curves have different turning points and performance trade-off, and it is difficult to decide the best procedure with consideration of multi-objective.
Sustainability 2020, 12, x FOR PEER REVIEW 10 of 16 trends. Similarly, Figure 4 illustrates the metrics-changing trends with the increase of target separation. It is worth noting that the curves have different turning points and performance tradeoff, and it is difficult to decide the best procedure with consideration of multi-objective.

Multi-Objective Evaluation Based on TOPSIS and Entropy
Taking the TOPSIS and entropy methods discussed above, the self-separation procedures can be evaluated with different treatments. For example, assume that both the target velocity and target separation have two options, which implies four alternative procedure could be used for the flow corridor. The relative closeness for the four treatments of self-separation is shown in Table 3. The procedure with standard cruise Mach number as the target velocity and 6 nmi as the target separation would be considered the best choice if the decision had only been made based on the average delay performance. However, when other objectives were incorporated, the procedure with 0.78 Mach as the target velocity and 6 nmi as the target separation was superior to the other three treatments of self-separation procedure. Table 3. Rank and relative closeness for four alternative self-separation procedures.

Multi-Objective Evaluation Based on TOPSIS and Entropy
Taking the TOPSIS and entropy methods discussed above, the self-separation procedures can be evaluated with different treatments. For example, assume that both the target velocity and target separation have two options, which implies four alternative procedure could be used for the flow corridor. The relative closeness for the four treatments of self-separation is shown in Table 3. The procedure with standard cruise Mach number as the target velocity and 6 nmi as the target separation would be considered the best choice if the decision had only been made based on the average delay performance. However, when other objectives were incorporated, the procedure with 0.78 Mach as the target velocity and 6 nmi as the target separation was superior to the other three treatments of self-separation procedure. Table 3. Rank and relative closeness for four alternative self-separation procedures. Moreover, a sensitivity analysis can be conducted to evaluate the impacts of alternative self-separation procedures in flow corridors with considering other self-separation parameters. Taking the lane-switch buffer and velocity difference threshold combination as example, we can assume that the target velocity ranged from 0.78 to 0.82 Mach, along with standard cruise Mach number, the target separation ranged from 6 to 12 nmi, and the traffic volume from 1X (statistical traffic demand, 81.76 aircraft/nmi/lane) to 4X (4 times traffic demand). Considering the lane-switch buffer (1 or 2 nmi) and the velocity difference threshold (20 or 40 knots), a total of 128 scenarios were considered, as shown in Table 4. For the 1X traffic density based on the interval arrival of A461, the ranks of different treatments of self-separation procedures with different lane-switch buffer (LB) and velocity difference threshold (VDT) are provided in Figure 5. Different combinations of LB and VDT are represented with different types of bars. For the self-separation procedure with the target velocity of 0.78 Mach as an example, if the target separation is 8 nmi, the 1 nmi LB with 40 knots VDT (treatment 2) has the best comprehensive performance than others. If the target separation is 6, 10, or 12 nmi, the 1 nmi LB with 40 knots VDT (treatment 4) is superior to other treatments. In general, with different target velocities, the rank of treatments changes with the target separation, but the 2 nmi LB with 40 knots VDT (treatment 4) obtains a greater number of the best comprehensive performances than others. Similar experiments were conducted for the 4X traffic density, as shown in Figure 6. For the selfseparation procedure with the target velocity of 0.78 Mach, if the Target Separation is 6, the 1 nmi LB with 40 knots VDT (treatment 4) is still superior to other treatments. However, if the target separation is 8 or 10, the 2 nmi LB with 20 knots VDT (treatment 3) is now superior to other treatments. In addition, if the target separation is 12, the 1 nmi LB with 20 knots VDT (treatment 1) becomes the best comprehensive performance. In general, with different target velocity, the 2 nmi LB with 20 knots VDT (treatment 4) cannot obtain as many the best comprehensive performances as 1X traffic density. Similar experiments were conducted for the 4X traffic density, as shown in Figure 6. For the self-separation procedure with the target velocity of 0.78 Mach, if the Target Separation is 6, the 1 nmi LB with 40 knots VDT (treatment 4) is still superior to other treatments. However, if the target separation is 8 or 10, the 2 nmi LB with 20 knots VDT (treatment 3) is now superior to other treatments. In addition, if the target separation is 12, the 1 nmi LB with 20 knots VDT (treatment 1) becomes the best comprehensive performance. In general, with different target velocity, the 2 nmi LB with 20 knots VDT (treatment 4) cannot obtain as many the best comprehensive performances as 1X traffic density.

Self-Separation Variables and
with 40 knots VDT (treatment 4) is still superior to other treatments. However, if the target separation is 8 or 10, the 2 nmi LB with 20 knots VDT (treatment 3) is now superior to other treatments. In addition, if the target separation is 12, the 1 nmi LB with 20 knots VDT (treatment 1) becomes the best comprehensive performance. In general, with different target velocity, the 2 nmi LB with 20 knots VDT (treatment 4) cannot obtain as many the best comprehensive performances as 1X traffic density.

Conclusions
This paper proposed a framework for multi-objective evaluation of airborne self-separation procedures in flow corridors, with consideration of the total control delay, flow corridor throughput, breakout rate, and the CO 2 emissions of traffic flow. Through combining multiple objectives into one with the TOPSIS and entropy method, we may obtain the optimum scheme to guide the design of self-separation procedures for flow corridors. The research results can be used by airspace managers to dynamically develop appropriate operational procedures and rules for flow corridors given different operational conditions and constraints. Also, the framework proposed in the research may be used to evaluate the design of airspace structure with consideration of multiple objectives.
One limitation of the study is that all variables evaluated for self-separation procedures are discrete from a practical standpoint. To achieve more efficient operational procedures for flow corridors, we can extend the problem to a multi-objective continuous-state optimization problem. In addition, since the safety should be the top consideration, we would like to filter the procedures with unacceptable safety performance in the future. In addition, our analysis was based on the data obtained from proposed microscopic simulation model. Although the self-separation behaviors have been visually validated, it is assumed that all information can be obtained timely and accurately without any time-lag, pilot blunders, and/or wind prediction error. The simulation model can be further improved by incorporating more realistic and stochastic conditions.