Tobias Andersson Granberg, Ta0ana Polishchuk, Valen0n Polishchuk, Chris&ane Schmidt
Automatic Design of Aircraft Arrival Routes with Limited Turning - - PowerPoint PPT Presentation
Automatic Design of Aircraft Arrival Routes with Limited Turning - - PowerPoint PPT Presentation
Automatic Design of Aircraft Arrival Routes with Limited Turning Angle Tobias Andersson Granberg, Ta0ana Polishchuk, Valen0n Polishchuk, Chris&ane Schmidt Introduction: Air transportation, SIDs + STARs Grid-based IP formulation Experimental
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 2
Introduction: Air transportation, SIDs + STARs Grid-based IP formulation Experimental Study: Arlanda Airport Conclusion/Outlook
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 3
Air transportation:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 3
Air transportation:
- Significant growth over the last decades
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 3
Air transportation:
- Significant growth over the last decades
- International Air Transport Association (IATA) projected that the
number of passengers will double to reach 7 billion/year by 2034
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 3
Air transportation:
- Significant growth over the last decades
- International Air Transport Association (IATA) projected that the
number of passengers will double to reach 7 billion/year by 2034
- Terminal Maneuvering Area (TMA), i.e., the area surrounding one or
more neighboring aerodromes, is particularly affected by congestion
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 3
Air transportation:
- Significant growth over the last decades
- International Air Transport Association (IATA) projected that the
number of passengers will double to reach 7 billion/year by 2034
- Terminal Maneuvering Area (TMA), i.e., the area surrounding one or
more neighboring aerodromes, is particularly affected by congestion
- Design Arrival and Departure procedures ➜ higher throughput
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 3
Air transportation:
- Significant growth over the last decades
- International Air Transport Association (IATA) projected that the
number of passengers will double to reach 7 billion/year by 2034
- Terminal Maneuvering Area (TMA), i.e., the area surrounding one or
more neighboring aerodromes, is particularly affected by congestion
- Design Arrival and Departure procedures ➜ higher throughput
- In Air Traffic Management (ATM): humans-in-the-loop!
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 3
Air transportation:
- Significant growth over the last decades
- International Air Transport Association (IATA) projected that the
number of passengers will double to reach 7 billion/year by 2034
- Terminal Maneuvering Area (TMA), i.e., the area surrounding one or
more neighboring aerodromes, is particularly affected by congestion
- Design Arrival and Departure procedures ➜ higher throughput
- In Air Traffic Management (ATM): humans-in-the-loop!
๏ Planes constantly monitored/guided by air traffic controllers (ATCOs)
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 3
Air transportation:
- Significant growth over the last decades
- International Air Transport Association (IATA) projected that the
number of passengers will double to reach 7 billion/year by 2034
- Terminal Maneuvering Area (TMA), i.e., the area surrounding one or
more neighboring aerodromes, is particularly affected by congestion
- Design Arrival and Departure procedures ➜ higher throughput
- In Air Traffic Management (ATM): humans-in-the-loop!
๏ Planes constantly monitored/guided by air traffic controllers (ATCOs) ๏ Safe separation between aircraft
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 3
Air transportation:
- Significant growth over the last decades
- International Air Transport Association (IATA) projected that the
number of passengers will double to reach 7 billion/year by 2034
- Terminal Maneuvering Area (TMA), i.e., the area surrounding one or
more neighboring aerodromes, is particularly affected by congestion
- Design Arrival and Departure procedures ➜ higher throughput
- In Air Traffic Management (ATM): humans-in-the-loop!
๏ Planes constantly monitored/guided by air traffic controllers (ATCOs) ๏ Safe separation between aircraft ➡ Route design should:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 3
Air transportation:
- Significant growth over the last decades
- International Air Transport Association (IATA) projected that the
number of passengers will double to reach 7 billion/year by 2034
- Terminal Maneuvering Area (TMA), i.e., the area surrounding one or
more neighboring aerodromes, is particularly affected by congestion
- Design Arrival and Departure procedures ➜ higher throughput
- In Air Traffic Management (ATM): humans-in-the-loop!
๏ Planes constantly monitored/guided by air traffic controllers (ATCOs) ๏ Safe separation between aircraft ➡ Route design should: ๏ lead to traffic patterns with “low complexity”
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 3
Air transportation:
- Significant growth over the last decades
- International Air Transport Association (IATA) projected that the
number of passengers will double to reach 7 billion/year by 2034
- Terminal Maneuvering Area (TMA), i.e., the area surrounding one or
more neighboring aerodromes, is particularly affected by congestion
- Design Arrival and Departure procedures ➜ higher throughput
- In Air Traffic Management (ATM): humans-in-the-loop!
๏ Planes constantly monitored/guided by air traffic controllers (ATCOs) ๏ Safe separation between aircraft ➡ Route design should: ๏ lead to traffic patterns with “low complexity” ๏ avoid creating conflict points
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 4
At most airports predesigned standard routes for departure and arrival:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 4
At most airports predesigned standard routes for departure and arrival: Standard Instrument Departures (SIDs) and Standard Terminal Arrival Routes (STARs)
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 4
017° 30’ AMDT 80 29 SEP 2005 CHANGE: SID to ARS, DKR, NOSLI and TRS Swedish Civil Aviation Authority AD 2–ESSA–4–21 RWY 01L AIP-SVERIGE/SWEDEN STOCKHOLM/ARLANDA AERODROME FMS/RNAV SID 018° 00’ 60° 00’ 59° 30’ 018° 00’ 017° 30’ 59° 30’ 60° 00’ 018° 30’ 018° 30’ K O G A V 3 C ( K O G A 3 C ) ( 3 2 . 7 ° T ) < 2 6 . 2 > SA403 RESNA 3 C (RESN 3C) KOGAV SA850 BABAP DKR 116.80 DVOR DUNKER 591225.8N 0170043.5E NOSLI LEGEND See GEN 2.3 Tracks are in MAG. Tracks within brackets are in True. ELEV and ALT in ft MSL VAR 3.5°-4.5° E 2005 Fly-over wpt Fly-by wpt INITIAL CLIMB CLEARANCE Common to all SIDs published on this chart. Unless otherwise specified, climb to 5000 ft. 585616.5N 0173008.0E TROSA DVOR/DME TRS 114.30 90X elev 213 ft 593510.3N 0163901.4E AROS DVOR/DME ARS 112.80 75X elev 50 ft 594138.3N 0180335.6E ANE 113.30 80X elev 108 ft 594247.8N 0175109.2E DME ANW 112.05 57Y elev 163 ft 593515.7N 0174910.9E DME ASW 113.75 84Y elev 232 ft 594459.3N 0184600.6E NORTEL VOR/DME NTL 116.30 110X elev 68 ft DME 593154.1N 0181212.0E TEBBY DVOR/DME TEB 117.10 118X elev197 ft 593912.4N 0175451.9E ARLANDA DVOR/DME ARL 116.00 107X elev141 ft RESNA (003.4° T) < 31.2 > SA401 ( 7 9 . 9 ° T ) < 6 . 2 > ( 9 . 3 ° T ) < 1 8 . 7 > SA402 N O R T E L 3 C ( N T L 3 C ) (003.4° T) < 6.8 > (010.4° T) < 4.0 > SA421 (ARL DME 1.3) SA851 < 4 . > < 6.2 > SA723 < 4 . 5 > (224.1° T) AROS 4 C (ARS 4C) SA701 2 6 3 ° ( 2 6 6 . 9 ° T ) MENGA 1 C (MENA 1C) ( 9 9 . 3 ° T ) < 2 8 . 7 > TROSA 4 C (TRS 4C) 180° (183.1° T) < 35.2 > ( 1 9 7 . 3 ° T ) < 7 . > (118.4° T) < 6.9 > SA702 SA703 BABAP NEKLA < 30.6 > (145.7°) 593813.9N 0175726.5E DME ASE 114.45 91Y elev 141 ft 2 9 9 ° 359° 359° 7 6 ° 006° 2 4 9 ° (252.6° T) 221° 1 4 1 ° 9 5 ° 115° 194° 8 6 ° SA422 MENGA BABAP 3 C (BABA 3C) 153° (157.1°T) < 21.5 > < 3 . > < 27.5 > DUNKER 4C (DKR 4C) 218° (221.9° T) < 25.4 > <2.2> N O S L I 4 C ( N O S L 4 C ) 1 9 6 ° ( 1 9 9 . 3 ° T ) < 2 8 . 6 > SA724At most airports predesigned standard routes for departure and arrival: Standard Instrument Departures (SIDs) and Standard Terminal Arrival Routes (STARs)
AIRAC AMDT 5/2012 23 AUG 2012 CHANGE: New STAR NILUG 1M. Deleted TRS 3M, VAR, pagina LFV AD 2–ESSA–4–5 RWY 01L/01R AIP SWEDEN STOCKHOLM/ARLANDA AERODROME STAR Instrument LEGEND See GEN 2.3 BRG are MAG ELEV and ALT in ft MSL km 10 0 10 20 30 km NM 5 0 10 20 NM Description of ELTOK, NILUG and XILAN see AD 2–ESSA–4–6 V A R 4 °- 5
- M
- 2
STAR Stockholm, RWY 01L/01R SID Stockholm, RWY 01L
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 5
SIDs/STARs:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 5
SIDs/STARs:
- Designed manually
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 5
SIDs/STARs:
- Designed manually
- No optimal routes for any specific criteria
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 5
SIDs/STARs:
- Designed manually
- No optimal routes for any specific criteria
- here: mathematical programming framework for finding
- ptimal STAR merge trees
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA location and direction of the airport runway
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA location and direction of the airport runway Output:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA location and direction of the airport runway Output: arrival tree that merges traffic from the entries to the runway,
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA location and direction of the airport runway Output: arrival tree that merges traffic from the entries to the runway, i.e., a tree that has the entries as leaves and the runway as the root
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA location and direction of the airport runway Output: arrival tree that merges traffic from the entries to the runway, i.e., a tree that has the entries as leaves and the runway as the root (arborescence oriented differently than usual)
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA location and direction of the airport runway Output: arrival tree that merges traffic from the entries to the runway, i.e., a tree that has the entries as leaves and the runway as the root (arborescence oriented differently than usual)
- 1. No more than two routes merge at a point: in-degree ≤ 2
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA location and direction of the airport runway Output: arrival tree that merges traffic from the entries to the runway, i.e., a tree that has the entries as leaves and the runway as the root (arborescence oriented differently than usual)
- 1. No more than two routes merge at a point: in-degree ≤ 2
- 2. Merge point separation: distance threshold L
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA location and direction of the airport runway Output: arrival tree that merges traffic from the entries to the runway, i.e., a tree that has the entries as leaves and the runway as the root (arborescence oriented differently than usual)
- 1. No more than two routes merge at a point: in-degree ≤ 2
- 2. Merge point separation: distance threshold L
- 3. No sharp turns: angle threshold 𝛽, minimum edge length L
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA location and direction of the airport runway Output: arrival tree that merges traffic from the entries to the runway, i.e., a tree that has the entries as leaves and the runway as the root (arborescence oriented differently than usual)
- 1. No more than two routes merge at a point: in-degree ≤ 2
- 2. Merge point separation: distance threshold L
- 3. No sharp turns: angle threshold 𝛽, minimum edge length L
- 4. Obstacle avoidance
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA location and direction of the airport runway Output: arrival tree that merges traffic from the entries to the runway, i.e., a tree that has the entries as leaves and the runway as the root (arborescence oriented differently than usual)
- 1. No more than two routes merge at a point: in-degree ≤ 2
- 2. Merge point separation: distance threshold L
- 3. No sharp turns: angle threshold 𝛽, minimum edge length L
- 4. Obstacle avoidance
- 5. STAR–SID separation:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA location and direction of the airport runway Output: arrival tree that merges traffic from the entries to the runway, i.e., a tree that has the entries as leaves and the runway as the root (arborescence oriented differently than usual)
- 1. No more than two routes merge at a point: in-degree ≤ 2
- 2. Merge point separation: distance threshold L
- 3. No sharp turns: angle threshold 𝛽, minimum edge length L
- 4. Obstacle avoidance
- 5. STAR–SID separation:
STAR–SID crossings far from the runway,
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
6
Input: locations of the entry points to the TMA location and direction of the airport runway Output: arrival tree that merges traffic from the entries to the runway, i.e., a tree that has the entries as leaves and the runway as the root (arborescence oriented differently than usual)
- 1. No more than two routes merge at a point: in-degree ≤ 2
- 2. Merge point separation: distance threshold L
- 3. No sharp turns: angle threshold 𝛽, minimum edge length L
- 4. Obstacle avoidance
- 5. STAR–SID separation:
STAR–SID crossings far from the runway,
where arriving and departing planes sufficiently separated vertically (difference of descend and climb slopes)
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
7
Objective function:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
7
Objective function: ๏ Short flight routes for aircraft
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
7
Objective function: ๏ Short flight routes for aircraft ➡ Minimize total length of the routes
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
7
Objective function: ๏ Short flight routes for aircraft ➡ Minimize total length of the routes ๏ STAR tree should "occupy little space"
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
7
Objective function: ๏ Short flight routes for aircraft ➡ Minimize total length of the routes ๏ STAR tree should "occupy little space" ➡ Minimize total length of the edges
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
7
Objective function: ๏ Short flight routes for aircraft ➡ Minimize total length of the routes ๏ STAR tree should "occupy little space" ➡ Minimize total length of the edges paths length
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
7
Objective function: ๏ Short flight routes for aircraft ➡ Minimize total length of the routes ๏ STAR tree should "occupy little space" ➡ Minimize total length of the edges paths length tree weight
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Optimal STAR merge trees
7
Objective function: ๏ Short flight routes for aircraft ➡ Minimize total length of the routes ๏ STAR tree should "occupy little space" ➡ Minimize total length of the edges paths length tree weight Pareto frontier of multicriteria optimization problem: set of Pareto optimal solutions (cannot be improved with respect to one of the objectives without sacrificing on the other)
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 8
Grid-based IP formulation
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA ๏ Snap locations of the entry points and the runway onto the grid
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA ๏ Snap locations of the entry points and the runway onto the grid ๏ P: set of (snapped) entry points
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA ๏ Snap locations of the entry points and the runway onto the grid ๏ P: set of (snapped) entry points ๏ r: runway
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA ๏ Snap locations of the entry points and the runway onto the grid ๏ P: set of (snapped) entry points ๏ r: runway ๏ Side of the grid pixel: L (➜merge point separation)
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA ๏ Snap locations of the entry points and the runway onto the grid ๏ P: set of (snapped) entry points ๏ r: runway ๏ Side of the grid pixel: L (➜merge point separation) ๏ G = (V,E):
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA ๏ Snap locations of the entry points and the runway onto the grid ๏ P: set of (snapped) entry points ๏ r: runway ๏ Side of the grid pixel: L (➜merge point separation) ๏ G = (V,E): ๏ Every grid node connected to its 8 neighbors
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA ๏ Snap locations of the entry points and the runway onto the grid ๏ P: set of (snapped) entry points ๏ r: runway ๏ Side of the grid pixel: L (➜merge point separation) ๏ G = (V,E): ๏ Every grid node connected to its 8 neighbors ๏ G is bi-directed
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA ๏ Snap locations of the entry points and the runway onto the grid ๏ P: set of (snapped) entry points ๏ r: runway ๏ Side of the grid pixel: L (➜merge point separation) ๏ G = (V,E): ๏ Every grid node connected to its 8 neighbors ๏ G is bi-directed ๏ Only exceptions:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA ๏ Snap locations of the entry points and the runway onto the grid ๏ P: set of (snapped) entry points ๏ r: runway ๏ Side of the grid pixel: L (➜merge point separation) ๏ G = (V,E): ๏ Every grid node connected to its 8 neighbors ๏ G is bi-directed ๏ Only exceptions: ๏ Entry points (no incoming edges)
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA ๏ Snap locations of the entry points and the runway onto the grid ๏ P: set of (snapped) entry points ๏ r: runway ๏ Side of the grid pixel: L (➜merge point separation) ๏ G = (V,E): ๏ Every grid node connected to its 8 neighbors ๏ G is bi-directed ๏ Only exceptions: ๏ Entry points (no incoming edges) ๏ r (no outgoing edges)
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA ๏ Snap locations of the entry points and the runway onto the grid ๏ P: set of (snapped) entry points ๏ r: runway ๏ Side of the grid pixel: L (➜merge point separation) ๏ G = (V,E): ๏ Every grid node connected to its 8 neighbors ๏ G is bi-directed ๏ Only exceptions: ๏ Entry points (no incoming edges) ๏ r (no outgoing edges) ๏ length of an edge (i, j)
`i,j
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
9
๏ Square grid in the TMA ๏ Snap locations of the entry points and the runway onto the grid ๏ P: set of (snapped) entry points ๏ r: runway ๏ Side of the grid pixel: L (➜merge point separation) ๏ G = (V,E): ๏ Every grid node connected to its 8 neighbors ๏ G is bi-directed ๏ Only exceptions: ๏ Entry points (no incoming edges) ๏ r (no outgoing edges) ๏ length of an edge (i, j) ๏ IP formulation is based on flow IP formulation for Steiner trees (Min Cost Flow Steiner arborescence)
`i,j
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
10
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
10
decision variables: edge e participates in the STAR.
xe
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
10
decision variables: edge e participates in the STAR. flow variables: gives the flow on edge e = (i, j) (i.e., from i to j )
xe fe
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
10
decision variables: edge e participates in the STAR. flow variables: gives the flow on edge e = (i, j) (i.e., from i to j )
xe fe
X
k:(k,i)∈E
fki − X
j:(i,j)∈E
fij = 8 > < > : |P| i = r −1 i ∈ P i ∈ V \ {P ∪ r} xe ≥ fe N ∀e ∈ E fe ≥ 0 ∀e ∈ E xe ∈ {0, 1} ∀e ∈ E
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
11
Objective functions: min X
e∈E
`efe (1) min X
e∈E
`exe (2)
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
11
Objective functions: min X
e∈E
`efe (1) min X
e∈E
`exe (2)
paths length
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
11
Objective functions: min X
e∈E
`efe (1) min X
e∈E
`exe (2)
paths length tree weight
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
12
Degree constraints:
X
k:(k,i)∈E
xki ≤ 2 ∀i ∈ V \ {P ∪ r} X
j:(i,j)∈E
xij ≤ 1 ∀i ∈ V \ {P ∪ r} X
k:(k,r)∈E
xkr = 1 X
j:(r,j)∈E
xrj ≤ 0 X
k:(k,i)∈E
xki ≤ 0 ∀i ∈ P X
j:(i,j)∈E
xij = 1 ∀i ∈ P
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
13
Turn angle constraints:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
13
Turn angle constraints:
Ae
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
13
Turn angle constraints:
Ae
ae = |Ae|
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
13
Turn angle constraints: aexe + X
f∈Ae
xf ≤ ae ∀e ∈ E
Ae
ae = |Ae|
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Grid-based IP formulation
14
SID constraints:
We disallow STAR edges to intersect SID edges within distance d from the runway.
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 15
Experimental Study: Arlanda Airport
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
16
Stockholm TMA:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
16
Stockholm TMA: ๏ Arlanda’s runway 19L
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
16
Stockholm TMA: ๏ Arlanda’s runway 19L ๏ Four main entry points: NILUG, XILAN, HMR, and ARS
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
16
Stockholm TMA: ๏ Arlanda’s runway 19L ๏ Four main entry points: NILUG, XILAN, HMR, and ARS ๏ Square grids of size 14x20 and 25x30
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
16
Stockholm TMA: ๏ Arlanda’s runway 19L ๏ Four main entry points: NILUG, XILAN, HMR, and ARS ๏ Square grids of size 14x20 and 25x30 ๏ Solve IP with
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
16
Stockholm TMA: ๏ Arlanda’s runway 19L ๏ Four main entry points: NILUG, XILAN, HMR, and ARS ๏ Square grids of size 14x20 and 25x30 ๏ Solve IP with ❖ Both objective functions
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
16
Stockholm TMA: ๏ Arlanda’s runway 19L ๏ Four main entry points: NILUG, XILAN, HMR, and ARS ๏ Square grids of size 14x20 and 25x30 ๏ Solve IP with ❖ Both objective functions ❖ Degree constraints
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
16
Stockholm TMA: ๏ Arlanda’s runway 19L ๏ Four main entry points: NILUG, XILAN, HMR, and ARS ๏ Square grids of size 14x20 and 25x30 ๏ Solve IP with ❖ Both objective functions ❖ Degree constraints ❖ Turn Angle constraints
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
16
Stockholm TMA: ๏ Arlanda’s runway 19L ๏ Four main entry points: NILUG, XILAN, HMR, and ARS ๏ Square grids of size 14x20 and 25x30 ๏ Solve IP with ❖ Both objective functions ❖ Degree constraints ❖ Turn Angle constraints ๏ 8 grid directions
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
16
Stockholm TMA: ๏ Arlanda’s runway 19L ๏ Four main entry points: NILUG, XILAN, HMR, and ARS ๏ Square grids of size 14x20 and 25x30 ๏ Solve IP with ❖ Both objective functions ❖ Degree constraints ❖ Turn Angle constraints ๏ 8 grid directions ๏ Postprocessing for smoother paths:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
16
Stockholm TMA: ๏ Arlanda’s runway 19L ๏ Four main entry points: NILUG, XILAN, HMR, and ARS ๏ Square grids of size 14x20 and 25x30 ๏ Solve IP with ❖ Both objective functions ❖ Degree constraints ❖ Turn Angle constraints ๏ 8 grid directions ๏ Postprocessing for smoother paths: shortcuts by removing vertices as long as the turn angle constraint is not violated
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
16
Stockholm TMA: ๏ Arlanda’s runway 19L ๏ Four main entry points: NILUG, XILAN, HMR, and ARS ๏ Square grids of size 14x20 and 25x30 ๏ Solve IP with ❖ Both objective functions ❖ Degree constraints ❖ Turn Angle constraints ๏ 8 grid directions ๏ Postprocessing for smoother paths: shortcuts by removing vertices as long as the turn angle constraint is not violated
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
17
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
18
222 224 226 228 230 232 234 286 288 290 292 294 296 298 300 302 304
Tree Weight Paths Length
Pareto frontier:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
18
222 224 226 228 230 232 234 286 288 290 292 294 296 298 300 302 304
Tree Weight Paths Length
Pareto frontier: Pareto optimal solutions:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
19
Obstacle avoidance:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
19
Obstacle avoidance:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
19
Obstacle avoidance:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
19
Obstacle avoidance:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight serve the airlines’ request for short trajectories best
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight serve the airlines’ request for short trajectories best quite dense network
- f routes ➜
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight serve the airlines’ request for short trajectories best quite dense network
- f routes ➜
hard to control the traffic
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight most merge points are located close to TMA boundary serve the airlines’ request for short trajectories best quite dense network
- f routes ➜
hard to control the traffic
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight most merge points are located close to TMA boundary serve the airlines’ request for short trajectories best quite dense network
- f routes ➜
hard to control the traffic
➜helpful to use linear combination of these two functions
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight
➜helpful to use linear combination of these two functions
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight
➜helpful to use linear combination of these two functions
Solutions for large number of entry points could be used to suggest the number and location
- f entry points for a
design from scratch.
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight
➜helpful to use linear combination of these two functions
Solutions for large number of entry points could be used to suggest the number and location
- f entry points for a
design from scratch. 2 entry points
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
20
Increased Number of Entry Points:
paths length tree weight
➜helpful to use linear combination of these two functions
Solutions for large number of entry points could be used to suggest the number and location
- f entry points for a
design from scratch. 2 entry points 16 entry points
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
21
SID constraints:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25Each tree: within approximately 2 CPU hours (105 B&B nodes)
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
22
SID constraints:
200 210 220 230 240 250 260 270 280 290 1 2 3 4 5 6 7 8 9 10 11 12 13
Paths Length Radius
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
22
SID constraints:
200 210 220 230 240 250 260 270 280 290 1 2 3 4 5 6 7 8 9 10 11 12 13
Paths Length Radius
d 1
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle
Experimental Study: Arlanda Airport
22
SID constraints:
200 210 220 230 240 250 260 270 280 290 1 2 3 4 5 6 7 8 9 10 11 12 13
Paths Length Radius
d 1
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 23
Conclusion/Outlook
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 24
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 24
Proof of concept for our grid-based IP approach for finding aircraft arrival routes with limited turning angle
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 24
Proof of concept for our grid-based IP approach for finding aircraft arrival routes with limited turning angle Easily integrates constraints from the departure routes
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 24
Proof of concept for our grid-based IP approach for finding aircraft arrival routes with limited turning angle Easily integrates constraints from the departure routes Static obstacles, e.g., no-fly zones, can be added
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 24
Proof of concept for our grid-based IP approach for finding aircraft arrival routes with limited turning angle Easily integrates constraints from the departure routes Static obstacles, e.g., no-fly zones, can be added Might choose to minimize weighted version: minimize the sum of trajectory lengths flown by all arriving aircraft (easily integrated by changing right-hand side of first equation)
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 24
Proof of concept for our grid-based IP approach for finding aircraft arrival routes with limited turning angle Easily integrates constraints from the departure routes Static obstacles, e.g., no-fly zones, can be added Might choose to minimize weighted version: minimize the sum of trajectory lengths flown by all arriving aircraft (easily integrated by changing right-hand side of first equation) Outlook:
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 24
Proof of concept for our grid-based IP approach for finding aircraft arrival routes with limited turning angle Easily integrates constraints from the departure routes Static obstacles, e.g., no-fly zones, can be added Might choose to minimize weighted version: minimize the sum of trajectory lengths flown by all arriving aircraft (easily integrated by changing right-hand side of first equation) Outlook: Simultaneous design of both SIDs and STARs
ATMOS, 25.08.2016 Automatic Design of Aircraft Arrival Routes with Limited Turning Angle 24
Proof of concept for our grid-based IP approach for finding aircraft arrival routes with limited turning angle Easily integrates constraints from the departure routes Static obstacles, e.g., no-fly zones, can be added Might choose to minimize weighted version: minimize the sum of trajectory lengths flown by all arriving aircraft (easily integrated by changing right-hand side of first equation) Outlook: Simultaneous design of both SIDs and STARs 3D routes