Algorithm For Scheduling Runway Operations Under Constrained Position Shifting
Algorithm For Scheduling Runway Operations Under Constrained Position Shifting
Algorithm For Scheduling Runway Operations Under Constrained Position Shifting
1. Introduction
The air transportation system is a tightly constrained system that is operating at (or close to) capacity at most
major airports. In 2005, terminal-area congestion accounted for only 13% of all delays at the 35 busiest
airports; that number had risen to 17% in 2008 and is currently at 21% over the first nine months of 2009.
The increasing delays, coupled with the expected increase in the demand for air transportation in the future,
have motivated several initiatives for the enhancement of terminal-area capacities (Arkind 2004, Boehme
1994). The runway system has been identified as the primary bottleneck in airport capacity, due to various
operational constraints on runway operations. Consequently, even small enhancements to runway throughput
can have a significant impact on system wide delays. The terminal area is a dynamic and uncertain
environment, with constant updates to aircraft states being obtained from surveillance systems and airline
reports (Atkins and Brinton 2002). The dynamic nature of the terminal area necessitates the development of
scheduling algorithms that are computationally efficient and therefore amenable to replanning when new
events occur, such as when a new aircraft enters the centre boundary or when data updates flexibility in
reordering aircraft, however large deviations from a nominal schedule may be unacceptable to airlines from a
fairness standpoint. This observation led to the constrained position shifting (CPS) approach for scheduling
aircraft, which stipulates that an aircraft may be moved up to a specified maximum number of positions from
its FCFS order. For example, if the maximum position shift (MPS) allowed was 2, an aircraft that is in the 8th
position in the FCFS sequence can be placed at the 6th, 7th, 8th, 9th or 10th position in the new sequence.
Several researchers have used CPS to model fairness and have worked toward developing fast solution
techniques for scheduling within the CPS framework
This report presents an algorithm for efficient runway scheduling on a single runway with CPS constraints,
while accounting for various operational considerations (including time-window restrictions and precedence
constraints. The core problem we consider is that of maximizing runway throughput (the landing time of the
last of a given set of aircraft) for arrivals-only or departures-only operations. Consequently, most practical
implementations resort to heuristic or approximate approaches that produce good solutions in a short-time.
The difficulty in solving these scheduling models arises primarily because the solution space allows for the
optimal sequence to deviate arbitrarily from the first-come-first-served (FCFS) sequence.
Keywords
FCFS, Constrained Position Shifting, Runway scheduling, Operations Research
2. Problem Definition
The runway scheduling problem is to find a sequence and corresponding arrival/departure times that optimize
Some objective of the schedule (for example, minimize the spending time in air or minimize a weighted
average of aircraft delay), subject to the following constraints.
1. Position Shift Constraint
As the airlines are major stakeholders in the air transportation system, it is important that an increase in
efficiency is not achieved at the expense of an equitable allocation of resources. This could happen if an
aircraft that would have had an early arrival or departure in the FCFS sequence is rescheduled to operate last,
thereby incurring a disproportionate amount of delay. CPS ensures some degree of fairness because it does
not allow the final sequence to deviate significantly from the FCFS order. The maximum number of position
shifts allowed is denoted by k, and the resultant scenario is referred to as a k-CPS scenario.
2. Minimum spacing requirements
An aircraft operating on a runway faces the risk of instability if it interacts with the wake-vortex of an aircraft
landing or taking off before it. To prevent this, the Federal Aviation Administration (FAA) mandates
minimum spacing requirements under instrument approach conditions (IAC) between aircraft operations on a
runway, which depend on the on the maximum take-off weight capacity of the aircraft. Although these
spacing requirements are specified in terms of distance, They can be converted to time requirements.
Representative values for these separations for three weight classes - small, large, and heavy
3. Time-window constraints
Limits on the levels of delay that can be incurred by an aircraft due to downstream traffic flow management
initiatives or constraints on possible manoeuvres that can be performed by the aircraft restrict the times at
which an aircraft can reach a runway. These constraints could possibly result in a set of disjoint time intervals
in which an aircraft can arrive/depart. For simplicity of notation, the case of a continuous-time interval
defined by an earliest and latest time, but the approach is applicable to disjoint intervals as well.
4. Precedence constraints
Precedence constraints are pairwise requirements on aircraft that stipulate whether one aircraft must land
before another. Sources of such constraints are the airlines themselves, which have precedence constraints
due to banking operations or priority flights. In addition, arrivals on the same jet route are constrained to not
overtake each other. Precedence constraints can also represent the restricted freedom available to taxiing
departures that are not allowed to overtake each other.
4. Complexity
The complexity of the algorithm for finding the minimum make span for n aircraft and maximum position
shift of k is O (n (2k + 1) ^ (2k+2)). The nodes in each stage of the network are generated by all combinations
of length 2k + 1, where each position in the sequence has at most 2k + 1 possible aircraft. The number of
nodes in each stage is therefore O ((2k + 1) ^ (2k+1)), because there are n stages, the total number of nodes in
the network is O (n (2k + 1) ^ (2k+1)). Each node can have at most 2k + 1 predecessors because the sequence
of a node differs from the sequence of its predecessor only in the first and last position. So the number of arcs
is O (n (2k + 1) ^ (2k+2)). Pruning the network requires looking at each arc at most twiceonce during the
forward pass and once during the backward pass. The dynamic programming recursion examines each arc at
most once, so total complexity is equal to the number of arcs in the network, which is O (n (2k + 1) ^
(2k+2)).In the presence of precedence constraints, each node would require O(k) work because we would
need to check if precedence is violated between the final aircraft in the node and each of the other aircraft in
the subsequence of the node. The complexity for pre-processing the entire network would thus be the number
of nodes times, which is the same as the running time of the pruning and dynamic programming recursion.
Although the complexity is exponential in k, it is of little consequence because k is typically small.
Therefore, the complexity of the algorithm is essentially linear in n.
6. Conclusion
The model represents a unified framework for runway scheduling under constrained position shifting and has
demonstrated that the problems of enhancing throughput, decreasing delay, and ensuring fairness can be
effectively modelled and solved in polynomial time (linearly in the number of aircraft) while accounting for
most operational constraints. An extension of the framework to include more general cost functions (by using
discrete-time models) and to mixed arrival and departure operations, including the merging of multiple
departure queues can also be developed. The algorithms can be easily implemented, and a prototype
implementation for arrival scheduling demonstrates that the run times are sufficiently small to enable
real-time deployment.
7. References
7.1 Anagnostakis, I., J.-P. Clarke, D. Bhme, U. Vlckers. 2001. Runway operations planning and control:
Sequencing and scheduling. J. Aircraft 38(6) 988996.
7.2 Lee, H. 2008. Trade-off evaluation of scheduling algorithms for terminal-area air traffic control. Masters
thesis, Massachusetts Institute of Technology, Cambridge, MA.
8. Attachments
-input.txt
-filerunway.py
-output.txt