Discrete and Dynamic Optimization Problems in Operations Management 2014
View this Special IssueResearch Article  Open Access
Macroscopic Model and Simulation Analysis of Air Traffic Flow in Airport Terminal Area
Abstract
We focus on the spatiotemporal characteristics and their evolvement law of the air traffic flow in airport terminal area to provide scientific basis for optimizing flight control processes and alleviating severe air traffic conditions. Methods in this work combine mathematical derivation and simulation analysis. Based on cell transmission model the macroscopic models of arrival and departure air traffic flow in terminal area are established. Meanwhile, the interrelationship and influential factors of the three characteristic parameters as traffic flux, density, and velocity are presented. Then according to such models, the macro emergence of traffic flow evolution is emulated with the NetLogo simulation platform, and the correlativity of basic traffic flow parameters is deduced and verified by means of sensitivity analysis. The results suggest that there are remarkable relations among the three characteristic parameters of the air traffic flow in terminal area. Moreover, such relationships evolve distinctly with the flight procedures, control separations, and ATC strategies.
1. Introduction
Air traffic management in terminal area is a knotty problem for controllers as this place is considered to be the air congestion, flight delay, and aviation accidentprone area. Researching on the fundamental operating features of terminal area traffic flow and deducing the macro emergence of traffic flow evolution may contribute to revealing the parameters with mutual relations, as well as the mechanism of spatiotemporal evolution, in terms of the traffic flow characteristic elements. By these means, we can move forward to exploring the objective law in air traffic in order to enrich the air traffic flow theory and to provide scientific basis for air traffic dispersion, which may have very important theoretical value and practical significance.
Traffic flow parameters are the physical variables that represent traffic flow characteristics including qualitative and quantitative features of operating states [1]. Basic theories of vehicle traffic flow have developed for decades and many results have been made by scholars. Lighthill and Whitham proposed the simulated dynamic model of traffic flow after researching on the evolution pattern of traffic flow under high traffic density circumstances [2]. Meanwhile, Richards proposed a first order continuum model of traffic flow, which has been integrated as the LWR theory [3]. Biham studied urban traffic flow based on a twodimensional cellular automaton [4], while Daganzo researched into dynamic traffic problems with a cellular transmission model [5–8]. Compared to the vehicle traffic, less research has been devoted to air traffic flow theory so far, not to mention that many studies focusing just on modeling. A simplified Eulerian network model of air traffic flow was proposed by Menon et al. [9, 10], and Bayen et al. studied the liner control problems derived from Eulerian network model [11–13]. Laudeman et al. noted a quantitative mathematical model on dynamic density of air traffic flow [14]. Complexity model based on traffic flow disturbance was advocated by Lee et al. [15]. Liu et al. proposed a onedimensional cellular transmission model specifically applicable to air route [16]. Wang et al. studied the microscopic planefollowing performance and built the air freeway flow model [17]. Primary discussion for the stability of air traffic flow operating system was advocated by Zhang and Wang, while some basic characteristics of air traffic flow were also involved in this paper [18]. Such research findings have made great foundations for further study on air traffic flow theory. However, no detailed studies have investigated the characteristic parameters with their objective evolution law of air traffic flow. In this paper, we will combine traditional mathematical formula derivation with modern simulation techniques. A macroscopic model of air traffic flow in airport terminal area will be proposed, with which we will simulate and analyze the interrelation and influential evolvement law of the characteristic parameters with the goal of providing theoretical basis for scientific air traffic management.
2. Macroscopic Model
2.1. Definitions
There are many kinds of definitions for air traffic flow parameters, because of different research intentions and different methods that can be used. Since we focus on airport terminal area, the definitions of velocity, density, and traffic flux on a segment of air route are as follows.
Traffic flux () is the number of aircraft passing a reference profile of the observation segment per unit of time. , where is the observation time and represents the number of aircraft passing in .
Density () is the number of aircraft per unit length of the observation segment. , where means length of the observation segment, while is the number of aircraft in .
Velocity can be divided into micro and macro definitions. Microscopic definitions including instantaneous velocity () and average velocity () focus on some point or profile of the observation segment; macroscopic ones including the space mean velocity () and time mean velocity () focus on some area extents or time ranges:
One flight () of all () uses time to pass the observation segment, where the actual distance is . And is the flying distance of in the time period .
Historically, the first macroscopic traffic flow model is a continuity equation, called the LighthillWhithamRichards (LWR) equation [2, 3]. Like the vehicle traffic flow, air traffic flux () and linear density () may also satisfy the corresponding relations for some given functions and with location and time , as follows:
On the right side of (3), denotes the number of aircraft that enter or exit the observation segment. With the basic equation , where denotes space mean velocity and assuming every flight has the same flying case, we can derive
There are two variables as density and velocity but only one equation in the LWR model. It cannot be solved as the differential equation is not closed. In response to this problem, a balancing velocitydensity functional relationship, as , was introduced into LWR theory by assuming traffic flow invariably in equilibrium state. To plug this into the equation, a hyperbolic equation of density can be derived, as follows:
2.2. Arrival Traffic Flow Model
LWR model described the propagation characteristic of the nonlinear density wave to find the evolution rule of traffic shock wave and rarefaction wave by characteristics method or numerical simulation [19]. However, air traffic flow differs from vehicle traffic, especially in airport terminal area (Figure 1). First, the arrival and departure traffic in terminal area will follow the designed STAR/SID (standard terminal arrival route/standard instrument departure), which means that aircraft in every route position must act in accordance with the operational flight program, that is, flying within a certain scope of designed flight level and speed. It is usually a small scope and can be reassigned by controllers. As a result of that, some balancing velocitydensity functional relationships in vehicle traffic may not exist in air traffic flow. Second, the density of air traffic flow is generally much lower than the vehicle traffic in real operations, which causes the interactions among aircraft to be weaker compared to vehicles. Therefore, it will be difficult to make statistical fit of the relationships according to the available radar data.
To solve the partial differential equation of the LWR theory, we use cell transmission model to discretize the continuity equation of macroscopic traffic flow. The method of discretization is applying a series of interconnected onedimensional cells to denote the air route and using the difference equation of time discretization to describe aircraft passing through every adjoining cell. It should be feasible to model any air traffic flow scenarios by using interconnected cells, as shown in Figure 2, putting all the aircraft flying along the arrival route while in several flight levels onto the same imaginary plane, hence the air traffic flow in one route can be approximately seen as onedimensional continuous flow. In addition, one air route will be divided into several segments by series of unit cells. To simplify the model we assume that there is a unified variable which we will introduce later to represent different kinds of control measures in one cell, such as speed control, maneuvering actions, or circling, that is, adjusting flows by changing speed or flight path of certain aircraft in the cell.
In the matter of arrival route, let be the number of arrival aircraft in cell at time ; then the changes in the number of arrival aircraft in one cell can be described by the following difference equation of time discretization:
In the equation above, is the number of arrival aircraft in cell at time , represents the flux of arrival aircraft entering cell from cell at time , while means the flux of arrival aircraft exiting cell at time , and is the time step.
Note. The number, flux et al. mentioned in this section, is specific to arrival aircraft, while departures will be included in following parts.
Air traffic flow in air route segment still satisfies the basic equation , where and stand, respectively, for segment liner density and space mean velocity of cell . The traffic flux of cell will be . Coefficient is the rate of outflow per unit time, which reflects the saturation level of cell . From Section 2.1., we know segment liner density , where is the length of cell . Since most of the time arrival aircraft in airport terminal area are in a deceleration process, we assume aircraft entering cell from some certain arrival route positions with an initial velocity that comes from the STAR and then uniformly decelerating along the cell (route) with a rate , so we can get the space mean velocity of cell , as follows:
To plug (7) into , we can derive
However, since the aforementioned situation is an ideal condition, it is necessary to take control measures for part of the arrival aircraft into consideration, such as speed control, maneuvering actions, or circling, because of the existence of traffic congestion, safety interval, and so forth in real conditions. For simplicity, we assume that there is a unified variable to stand for different kinds of control measures, since whichever measures controllers took, the actual effects caused by changing speed or flight path and so forth can all be seen as there to be aircraft flying along the cell (route) with a new constant velocity in cell . To be clear, the variable here represents the displacement velocity that is lower than the STAR designed velocity generally.
Considering these problems, we derived (9) after integrating control measures into (8), as follows:
After substituting in the difference equation (6) with (9), we got
Let ; ; ; we can have the simplified equation sets, as follows:
The aircraft in cell are taken as evenly distributed in the macroscopic traffic flow theory, so the number of aircraft should be equivalent number. Let be the mean nose interval of adjacent aircraft at the same direction in cell ; therefore the number of aircraft in the cell is . If becomes less than the control separation standard denoted by , that is to say, aircraft density exceeds the threshold level, the exceeded aircraft should be arranged to take control measures. In addition, arrival aircraft in airport terminal area also have to obey the safety separation according to STAR; that is, . The safety separation is evaluated by air traffic security department and usually gets more stringent compared to control separation as . The number of controlled aircraft at time in cell satisfies
We have discussed the onedimensional cell transmission model of air traffic flow before, while there may be multiple air routes in terminal area and some of them may cross with each other, in the way, in the same imaginary plane, as shown in Figure 3. So the converging and diverging situations are as shown in Figure 3.
According to the conservation of number, in converging situation the number of aircraft satisfies
On the contrary, in diverging situation we have
The coefficient stands for the proportion of traffic flow to different air route segments; meanwhile , .
2.3. Departure Traffic Flow Model
The model of departure flow in airport terminal area consists of two phases: taking off from runway and flying along air routes. In the takeoff phase, departure air traffic flow will make the most use of runway time slots based on the arrival priority, in order to guarantee a smooth arrival and landing process. Flight departure schedules are generated in stochastic cases, while the departure flights that disagree with the operation time interval of runway will be delayed to ground holding procedures until enough time slots come [20–22]. In the flying stage, departure aircraft diverge in various directions from the runway center, which differs from the arrival flow as there would be both diverging and converging situations in arrival processes. Under the condition of fully isolation between arrival and departure routes, flying stage will dispense with flow control if the runway interval problem is solved in the takingoff stage. And under the condition of semiisolation between arrival and departure routes, which means that some air routes may be overlapped in several segments (in the same imaginary plane), it is necessary to take both departure aircraft and arrival ones into consideration simultaneously to focus on the average nose interval of all aircraft in the overlapped segments (cells). If the average interval cannot meet safety needs, departure flow should be adjusted prior to the arrivals.
Let be the number of aircraft holding for departure on ground in time , and its change can be described as follows:
In (15), is the number of aircraft holding for departure on ground in time . The demand of departure per unit time that produced by flight schedules is represented by , while represents the actual number of aircraft taking off from the runway per unit time. The time step is still , and cannot be negative, so we can get
Consider maximizing the use of runway time slots; variable represents the maximum takeoff rate of the runway, which depends on the time interval of runway operation. If the time interval of arrival landing aircraft gets large, influence such as wake vortex on the runway will have no effect on takeoff aircraft, in which cases departure flows may take off by the standard time interval separately. Otherwise if arrival time interval gets intense, mutual interference between departures and arrivals will be strong, so that departure flows need to make use of the interspace of the time slots under the circumstances of an affected runway operation, as follows:
In (17), this is the actual time interval of arrival landing and is the standard takeoff time interval of departures, while and , respectively, stand for the irrelevant runway operation time interval and the modified operation time interval according to the mutual interference between landing and taking off. Among this, the irrelevant runway operation time will be a critical point where interference occurred.
Similar to the Arrival Model, we divide air routes into several segments by a series of unitinterconnected cells in the stage of flying along departure routes. Let be the number of departure aircraft in cell at time and its changing process can be described as (18). On the left side of the equation is the number of departures in cell at time ; represents the flux of departure aircraft entering cell from cell at time , while means that the flux of departure aircraft exiting cell at time , is the time step.
Since most of the time departure aircraft in airport terminal area are in an acceleration process, we assume aircraft entering cell from some certain departure route positions with an initial velocity that comes from the SID and then uniformly accelerating along the cell (route) with a rate . So according to the same modeling principle from the Arrival Model, we can get the following similar equations:
In the above equations, coefficient stands for the rate of departures outflow per unit time. The number of aircraft that need to be arranged to take departure flow control measures is denoted by the variable . As mentioned before, in circumstances of fully isolation between arrival and departure route segments there will be no extra control to the departures; that is, . On the other hand, in the routes overlapped condition part of the departures may take extra controls. We assume that these controlled departure flows would move with a new constant ATC velocity in cell and be the same with variable in the Arrival Model, in which both represent displacement velocity. Meanwhile the simplification form is
To determine the value of variable in overlapped segments, we will bring the average nose interval (not time interval) of arrivals noted by , the average nose interval of departures noted by , and the average interval between both arrivals and departures noted by in cell all into consideration. The fundamental aim of flow control in this stage is to make the entire average interval of arrivals and departures meet the ATC separation requirement by adjusting departure flows, as follows:
The average interval of arrivals and departures is .
3. Parameter Analysis
We focus on one single cell in the model as to deduce and analyze the interrelationship of air traffic flow characteristic parameters including flight flux, linear density, and traffic velocity. In the following sections, we discuss this problem from the two different aspects as arrival and departure, just like the models we had established before.
3.1. Arrival Routes
To plug (12) into (9), we can have
According to , we get
From (23), we can see that the basic parameters such as flux and density form a piecewise function and the segment point is the reciprocal value of ATC separation to the arrival aircraft. Assuming that coefficient , standard deceleration , initial speed , and length of one cell are all constant values, the slope and intercept of the relation curve will be determined by ATC velocity and the position of inflection point will be determined by . Let both sides of (23) be divided by :
If aircraft density comes lower than the critical value, the mean velocity of traffic flow in one cell will stay constant; otherwise, there will be an inverse relation between velocity and density . With an increase of density in cell , mean velocity will become lower gradually and the speed of reducing just gets more and more slow and eventually tends to the ATC velocity . Using the basic formula of fluid , we can get
From (25), we can see that the mean velocity of entire traffic flow in cell ranges from ATC velocity to standard velocity designed by STAR. If is equal to one of the critical values flux may be arbitrarysized data taking from zero to the maximum flux value. With an increase of mean velocity in cell , the flux value goes down gradually and the speed of reducing gets slow. When this mean velocity goes up to the peak value, it will cause a jump of traffic flux. Variables and are still the order parameters of the interrelationship between flux and velocity in cell.
3.2. Departure Routes
In the stage of flying along departure routes, aircraft keep accelerating and diverging along with the varying of orientation of air routes. There may be aircraft taken departure flow control measures in some of the overlapped segments. As an exceptional case of this situation, flying stage will dispense with flow control in the case of fully isolation between arrival and departure routes; that is, . We will focus on the overlapped segments and analyze the traffic flow characteristic parameters of departures in detail. To plug (21) into (18) we can have
Since , , substituting in (26), we get
From (27), we can see that flux and density form a piecewise function. Similar to the arrival segments (cells), we assume that coefficient , standard acceleration , initial speed , and length of one cell are all constant values. The density of arrival aircraft denoted by in cell changes over time, which has an effect on the intercept value of this linear equation. The ATC velocity of departures denoted by shows the slope while the reciprocal of ATC separation is the critical value of flow density, determining the position of inflection points. Then, letting both sides of (27) be divided by , we got
From (28), we can see that mean velocity and density also form a piecewise function. Two extreme cases are as follows: the entire density of arrivals and departures is lower than the ATC density that is denoted by ; density only considers arrivals that have already exceeded . These two cases will cause mean velocity of departures to be the standard velocity determined by , acceleration , initial speed , length , and departure ATC velocity , respectively. In a certain condition between such two extreme cases the slope of linear equation is greater than or equal to , while the intercept value is influenced by both and . According to the basic formula of fluid , we get
The form of (29) is much like (25) from the arrival parts. If mean velocity of departures lies between the standard velocity and the departure ATC velocity there will be an inverse proportional function between departure flux and mean velocity . Moreover, the coefficient of this inverse proportional function is determined by , acceleration , initial speed , and cell length all together and changes in value of the ATC velocity will make the function curve move horizontally. While if mean velocity of departures lies on either of the two critical positions departure flux may vary randomly.
4. Simulation Experiment
4.1. Simulation Sample
Based on the simulation platform NetLogo [23–25] we focused on each cell in the macroscopic model and designed traffic inflow/outflow behaviors to be each agent, so as to take control of the inflow/outflow volumes among all the cells. Applying the STAR/SID procedures to runway 02L of Guangzhou Baiyun International Airport (as shown in Figure 4) into the Netlogo system dynamic simulator (as shown in Figure 5) we can emulate the operation and evolution process of air traffic flow in the airport terminal area.
The network of arrival and departure routes in terminal can be separated into several single direction segments and converging or diverging segments that can be further decomposed into more unit cells. For the purpose of simplicity, we set one cell dimension to be equal to the length of the shortest segment in the network and the other air segments to be integer multiples of unit cell dimension. According to this, we made the following Tables 1 and 2 of the terminal network.


Each yellow rectangle in Figure 4 represents one cell, which means a “stock” of fluid in the simulator system. Each grey pipe represents flowing between adjacent cells, which means “flow” in system, and the direction of arrow is the same as the flow direction. Each black valve represents control measures, which means “strategy” in system, and the inflow and outflow of air traffic will be under control by giving the certain valves some corresponding rules. Take typical overlapped segments “TAN=AGVOS (Segment9&19)” as an example, the cell transmission model of this is like (30). In the model, respectively, represents the number of arrival aircraft and the simplified coefficients of cell in Segment . denotes the number of arrival aircraft taking control measures in cell of Segment . means the number of arrival aircraft that exit cell in Segment per unit time. respectively, represents the similar corresponding ones but specially for departures. Coefficient stands for the proportion of traffic flow into some segment when diverging occurs in air routes. is the designed time step for the simulator system:
4.2. Simulation Design
In this simulation sample there are 6 entry points for arrivals of the terminal network: MUBEL, OSIKA, BIPOP, P101, LONGTANG, and DANZHU. We assume that the arrival rate of each entry point obeys the negative exponential distribution [26]. To plug the average arrival rate that comes from flight historical statistics as expected value into distribution functions we can obtain the changes of traffic flux over time at each entry point as shown in Figure 6. It should be noted that we use equivalent traffic flow in this simulation and take each simulation time step as 1 minute.
It can be noticed that traffic flux changes significantly over time at each entry point. When aircraft keep going to convergent points, the peak of traffic wave may happen to meet another one, in which situation traffic density nearby will be too high to meet separation requirements and the probability of unsafe events will increase accordingly. Instead, the trough of traffic wave may also meet another trough that may lead to a low traffic density and a large aircraft interval at convergent points, thus reducing the time/space utilization of limited airspace. Based on the situation we designed a “valve”. Agent aims at balancing traffic flow in the simulator system to manage the inflow/outflow of adjacent cells. In addition, this is much like the principle of “Cutting peak and filling valley” in real air traffic flow management [27, 28].
The basic strategy of arrival “valve” control is as follows: first determine whether the mean arrival aircraft interval is lower than the arrival ATC separation in one cell at any time, which means whether the stock of aircraft in one cell exceeds. If yes then the excessive number of aircraft should operate with the assigned ATC velocity and the rest do not change their speed, while if no then all aircraft can still follow the STAR. Since is lower than normal velocity, the traffic flux will also be lower compared to no control situations. With outflow decreasing, the stock of aircraft in one cell will increase accordingly. At next time step inflow will be added to the original cell stock. If this total value still exceeds the standards, control measures should be taken like before. These cyclic steps keep going until sometime there comes a small inflow and the total value added with cell stock goes lower than the standards. In this circumstance all aircraft in the cell can operate with STAR and aircraft controlled before can exit cell normally. In the whole process, the total number of aircraft in one cell cannot exceed the safety value at any time; that is, mean interval should never be smaller than the safety separation.
Take typical arrival segment “GYAAGVOS (Segment5)” as an example and making a comparison of air traffic flux and stock in cells between before and after the “valve” control, it is easy to find that the inflow and outflow in this segment become smooth and steady when there is the “valve” control and the stock of aircraft in each cell always meets the safety requirement as shown in Figures 7(a) and 7(b). The unit is flights/minute.
(a) Before arrival control
(b) After arrival control
(c) Before departure control
(d) After departure control
The departure in terminal area consists of two major parts: taking off from runway and flying along departure routes. We designed two successive processes in this simulation accordingly called the airport surface part and the air routes part, which are complementary in departure process. Specifically, the airport surface part mainly consists of departure flights schedule generation, runway time occupation, takeoff slots allocation, and aircraft ground holding. We assume that the generation of departure flights schedule obeys a negative exponential distribution. The number of scheduled departure flights per unit time is as shown in Figure 8(a). The runway time occupation and takeoff slots allocation are actually how the landing aircraft and takingoff aircraft make the most use of runway slot resources within the limited and dynamic runway operation capacity. Changes in the number of landing and takingoff aircraft per unit time are as shown in Figures 8(b) and 8(c). Since landing aircraft have the priority, takingoff aircraft that are unable to go by flight schedules will be postponed to take groundholding processes. Changes in the number of ground holding aircraft over time are as shown in Figure 8(d).
(a) Scheduling
(b) Landing
(c) Taking off
(d) Gholding
In the process of flying along departure routes, there is no need to take control measures under the condition of full isolation between arrival and departure routes. The comparison of departure flux and stock over time before and after control measures is as shown in Figures 7(c) and 7(d). The basic strategy of departure “valve” control is as follows: first determine whether the mean arrival aircraft interval is lower than the ATC separation in one cell at any time. If yes then all the departures in this cell should take control measures, which means passing the segment with constant ATC velocity . If no then determine whether the mean interval of both arrivals and departures is lower than ATC separation in one cell, and if yes departures must be adjusted to increase the entire mean interval until it goes above .
4.3. Result Analysis
According to the export data derived from NetLogo simulation platform we got series of parameter scatter diagrams that could reflect the basic traffic flow characteristics through statistic and analysis, as shown in Figure 9. Based on this, we analyzed the mutual influence relationship among traffic flux , density , and velocity of the arrival and departure routes in airport terminal area. Limited by space this paper only made detailed discussion on relation of arrivals and departures; meanwhile the rest of and relations were listed in the form of statistical diagrams for reference.
(a)
(b)
Taking typical arrival segment, Segment5, as an example, we derived the basic tendency of traffic flux and density relationship for arrival routes in terminal area as shown in Figure 9(a1). From the diagram we can find that the relationship tendency consists of three main stages.
Stage I is the free flow state in which the number of aircraft in segment is very low and the mean interval exceeds arrival ATC separation , which means all the aircraft can follow STAR. The traffic flux is directly proportional to density in segment and the proportionality coefficient is equal to the mean velocity of segment. Stage II is the congestion flow state in which the number of aircraft in segment increases and then the mean interval is under arrival ATC separation , which means part of the aircraft should be assigned to take control measures including decelerating, maneuvering or holding, and so forth. The relationship between traffic flux and density occurs a inflection point. The traffic flux in segment is still directly proportional to density; however the new proportionality coefficient is equal to the space mean value of the controlled aircraft velocity and uncontrolled aircraft velocity; that is,
, see details in Section 2.1(4).
All of the control measures in the arrival process, including deceleration, maneuvering and holding, will lead to a decline in the displacement velocity. In the diagram it made the tendency of flux and density relationship to be leveling off. While the traffic flux will still increase with density, which differs from normal vehicle traffic since after inflection point the vehicle traffic flux decreases with density. The reason why this difference exists is that the carfollowing behavior between adjacent vehicles has significant influences when congestion occurs on the road while air traffic normally maintains a larger safe separation. Apart from this, some control measures as hold pattern make aircraft deviate from original air routes, which does not affect the other aircraft. Therefore, after the inflection point traffic flux will not decrease but increase with a low slope. Stage III is the block flow state in which the number of aircraft in segment exceeds safe value with the control adjustment in congestion state and large amount of traffic converging continuously. The mean interval is under the safe separation and insecurity factors surge that should be avoided as much as possible.
From the above theoretical derivation we can find that variables related to operational flight program including initial velocity , acceleration , and length of segment (cell) become constant values as the STAR/SID are established. This paper focuses on the rest of the variables especially flow control variables including ATC separation and ATC velocity . Such variables will become order parameters that influence the mutual relationships among three basic air traffic flow parameters.
As shown in Figure 9(a2), when arrival ATC separation rises to 25% the inflection point of traffic flux and density relationship moves forward and actually the horizontal axis value equals . The rise of ATC standard will lead to a decrease of free flow state in Stage I. More aircraft need to take control measures; meanwhile the frequency of high traffic density even block flow in Stage III also increases. But in whatever Stage I or Stage III, the slope of traffic flux and density relationship maintains constant; that is, before inflection point the increase tendency coincides and after inflection point the increase tendency still keeps parallel. As shown in Figure 9(a3), when arrival ATC separation remains the same and arrival ATC velocity turns down, the inflection point stays the same but aircraft pass the segment with a lower mean velocity after this point. Therefore, traffic flux has a decelerated growth with the density. Conversely, if it needs to let traffic flux reach the level before changing of ATC velocity , the density of air segment should be higher, which may lead to an early arrival of the block flow state in Stage III. The free flow state in Stage I has no change now since it is not affected by control measures.
Taking typical departure segment, Segment19, as an example, we derived the basic tendency of traffic flux and density relationship for departure routes in terminal area as shown in Figure 9(b1). The relationship tendency also consists of three main stages. Similar to arrival parts, Stage I is the free flow state in which the traffic flux of departures is directly proportional to density in segment and the proportionality coefficient is equal to the mean velocity of segment from SID. Stage II is the congestion flow state which is unlike arrival parts, since in departure routes we focus on the entire mean interval of both arrivals and departures, not just departures. With the density of departures rising up in segment, the entire mean interval can still meet the ATC separation if arrivals density is small enough. Then the tendency of traffic flux and density relationship lies on the extension line of Stage I. Stage III is the block flow state in which the number of arrivals and the number of departures increase simultaneously. Accordingly the entire mean interval becomes lower than the ATC separation . Considering the principle of arrivals priority part of the departures will accelerate to leave the heavytraffic segment in order to release the space resources for arrivals. The new slope of traffic flux and density relationship equals the departure ATC velocity . The number of adjusted departures is determined by arrivals, which is shown as a series of scatter values between the tendency line of Stage III and extension line of Stage II. The vertical distance from tendency line of Stage III to the scatter values varies inversely to the density of arrivals, that is, ( is constant coefficient); see details in (27).
As shown in Figure 9(b2), after raising up 25% of the ATC separation , there appear more departure aircraft that need to take control measures. Thus the Stage III also appears earlier while the slopes of tendency lines in congestion and block flow states have no change. However, there are more scatter values between these two lines. The Stage II congestion flow state in which entire mean interval (high departure, low arrival) still exceeds the ATC separation will reduce its frequencyofoccurrence. More scatter values lie in the Stage III block flow state. Conversely reducing the ATC separation departures may not need to be adjusted in most cases, which tend to be the Stage I free flow state. More scatter values of traffic flux and density lie more on the extension line of free flow state.
Keeping the ATC separation unchanged and changing the departure ATC velocity we can get new relationship tendencies as shown in Figure 9(b3). The Stage I free flow state keeps the same. But in Stage III block flow state the variation rate of traffic flux with density increases with the departure ATC velocity , which means departure flow has raised its outflow rate of heavytraffic segment. Meanwhile the impact of arrivals increases by ( is constant coefficient); see details in (27), which is shown as the scatter values spread a larger scope from center line of the free flow state extension direction.
5. Conclusions
Based on CTM we have proposed the macroscopic model of air traffic flow in airport terminal area and carried out a series of simulation experiments with the NetLogo platform. Through both of the theoretical and practical discussions, we could generally reveal the basic interrelationships and influential factors of air traffic flow characteristic parameters. The research findings are as follows.(1)The CTM could accurately reflect the macroevolution laws of air traffic flow in terminal air route network. Meanwhile, it may also be applied to different kinds of air traffic scenes including airways, sectors, or airspaces by modifying conditions correspondingly. Moreover, it has a remarkable operational efficiency in multiagent simulations.(2)There are obvious relationships among the three characteristic parameters as flux, density, and velocity of the air traffic flow in terminal area. And such relationships evolve distinctly with the flight procedures, control separations, and ATC strategies. The air traffic flow characteristics may take the specific changes through flight procedure optimization, control separation modification, or ATC strategy regulation, which could be part of the scientific basis for air traffic management in airport terminal area.(3)The default parameters we used in simulation experiments are from practical ATC rules. Automatic optimization of these parameters for desired traffic flow characteristics should be taken into consideration in further studies. In addition, discussions in this paper focus on macro perspectives, thus the research results seem rough anyway. To obtain more detailed and elaborate traffic flow characteristics it is necessary to combine such macro studies with micro perspectives that give full expressions to the following: overflying or turning and so forth of individual behaviors and interactive effects. It also should be an important direction for further study.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research is supported by “the National Natural Science Foundation of China (NSFC) no. 61104159” and “the Fundamental Research Funds for the Central Universities no. NJ20130019.”
References
 L. Li, R. Jiang, B. Jia, and X. M. Zhao, Modern Traffic Flow Theory and Application, vol. 1, Tsinghua University Press, Beijing, China, 2011.
 M. J. Lighthill and G. B. Whitham, “On kinematic waves. I. Floow movement in long rivers,” Proceedings of the Royal Society, London A, vol. 229, pp. 281–316, 1955. View at: Publisher Site  Google Scholar  MathSciNet
 M. J. Lighthill and G. B. Whitham, “On kinematic waves. II. A theory of traffic flow on long crowded roads,” Proceedings of the Royal Society: London A: Mathematical, Physical and Engineering Sciences, vol. 229, pp. 317–345, 1955. View at: Publisher Site  Google Scholar  MathSciNet
 O. Biham, A. A. Middleton, and D. Levine, “Selforganization and a dynamical transition in trafficflow models,” Physical Review A, vol. 46, no. 10, pp. R6124–R6127, 1992. View at: Publisher Site  Google Scholar
 C. F. Daganzo, “The cell transmission model: a dynamic representation of highway traffic consistent with the hydrodynamic theory,” Transportation Research B, vol. 28, no. 4, pp. 269–287, 1994. View at: Publisher Site  Google Scholar
 C. F. Daganzo, “The cell transmission model, part II: network traffic,” Transportation Research B, vol. 29, no. 2, pp. 79–93, 1995. View at: Publisher Site  Google Scholar
 C. F. Daganzo, “A variational formulation of kinematic waves: basic theory and complex boundary conditions,” Transportation Research Part B: Methodological, vol. 39, no. 2, pp. 187–196, 2005. View at: Publisher Site  Google Scholar
 C. F. Daganzo, “On the variational theory of traffic flow: wellposedness, duality and applications,” Networks and Heterogeneous Media, vol. 1, no. 4, pp. 601–619, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. K. Menon, G. D. Sweriduk, and K. D. Bilimoria, “Air traffic flow modeling, analysis and control,” in Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, pp. 11–14, 2003. View at: Google Scholar
 P. K. Menon, G. D. Sweriduk, and K. D. Bilimoria, “New approach for modeling, analysis, and control of air traffic flow,” Journal of Guidance, Control, and Dynamics, vol. 27, no. 5, pp. 737–744, 2004. View at: Publisher Site  Google Scholar
 A. M. Bayen, R. L. Raffard, and C. J. Tomlin, “Adjointbased control of a new Eulerian network model of air traffic flow,” IEEE Transactions on Control Systems Technology, vol. 14, no. 5, pp. 804–818, 2006. View at: Publisher Site  Google Scholar
 A. Bayen, R. Raffard, and C. Tomlin, Hybrid Control of PDE Driven Highway Networks, Lecture Notes in Computer Science, Springer, New York, NY, USA, 2004.
 R. Raffard, S. L. Waslander, A. Bayen, and C. Tomlin, “A cooperative, distributed approach to multiagent eulerian network control: application to air traffic management,” in Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, 2005, AIAA 20056050. View at: Google Scholar
 I. V. Laudeman, S. G. Shelden, and R. Branstrom, “Dynamic density: an air traffic management metric,” NASA/TM219982112226, 1998. View at: Google Scholar
 K. Lee, E. Feron, and A. Pritchett, “Air traffic complexity,” in Proceedings of the 44th Annual Allerton Conference, pp. 1358–1363, 2006. View at: Google Scholar
 Q. Liu, C. R. Bai, and Q. Lin, “Linearquadratic optimal control of air traffic flow,” Computer and Communications, vol. 26, no. 6, pp. 116–119, 2008. View at: Google Scholar
 L. Wang, X. Zhang, and Z. Zhang, “Following phenomenon and air freeway flow model,” Journal of Southwest Jiaotong University, vol. 47, no. 1, pp. 158–162, 2012. View at: Publisher Site  Google Scholar
 Z. N. Zhang and L. L. Wang, Air Traffic Flow Management Theory and Method, Science Press, Beijing, China, 2009.
 P. I. Richards, “Shock waves on the highway,” Operations Research, vol. 4, no. 1, pp. 42–51, 1956. View at: Publisher Site  Google Scholar  MathSciNet
 T. Vossen and M. Ball, “Optimization and mediated bartering models for ground delay programs,” Naval Research Logistics A: Journal Dedicated to Advances in Operations and Logistics Research, vol. 53, no. 1, pp. 75–90, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 T. W. M. Vossen and M. O. Bal, “Slot trading opportunities in collaborative ground delay programs,” Transportation Science, vol. 40, no. 1, pp. 29–43, 2006. View at: Publisher Site  Google Scholar
 T. Vossen, M. Ball, R. Hoffman, and M. Wambsganss, “A general approach to equity in traffic flow management and its application to mitigating exemption bias in ground delay programs,” in Proceedings of the 5th USA/Europe ATM R&D Seminar, Budapest, Hungary, 2003. View at: Google Scholar
 S. Tisue and U. Wilensky, “NetLogo: design and implementation of a multiAgent modeling environment,” in Proceedings of the Agent Conference on Social Dynamics: Interaction, Reflexivity and Emergence, Chicago, Ill, USA, 2004. View at: Google Scholar
 W. Rand and U. Wilensky, “Visualization tools for agentbased modeling in NetLogo,” in Proceedings of Agent Chicago, 2007. View at: Google Scholar
 G. Z. Jiao, C. Guo, X. J. Hu, Y. Qin, and X. Ou, “Research on selfgoverning traffic model based on micro behavior information decision,” Control and Decision, vol. 25, no. 1, pp. 64–68, 2010. View at: Google Scholar
 G. C. Ovuworie, J. Darzentas, and R. C. McDowell, “Free moves, followers and others: a reconsideration of headway distributions,” Traffic Engineering and Control, vol. 21, no. 89, pp. 425–428, 1980. View at: Google Scholar
 D. Bertsimas, G. Lulli, and A. Odoni, “The air traffic flow management problem: an integer optimization approach,” in The ATFM Problem: An Integer Optimization Approach, vol. 5035 of Lecture Notes in Computer Sciences, pp. 34–46, Springer, Berlin, Germany, 2008. View at: Google Scholar
 D. Bertsimas and S. S. Patterson, “The air traffic flow management problem with enroute capacities,” Operations Research, vol. 46, no. 3, pp. 406–422, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2014 Honghai Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.