Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Detect and Avoid”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,207 records · Page 67

Operating displays - new concept: high performance displays - 118

New control rooms in Nuclear Power Plants using digital Distributed Control System (DCS) operated with software displays, are providing thousands of signals to the operating crew. When an event happens, the Flood of information provided to operators could obscure the initiator event increasing the operator response time, work load and complicating the selection of the actions to be performed. There are some initiatives that are trying to minimize the 'avalanche effect' as the alarm engineering filtering and prioritization or diagnostic displays. Likewise, new control rooms are lost the availability to know the plant status just with a quick glance to the panel indicators. Regarding Small Modular Reactors, where in some cases four operators should monitor and control twelve reactors, these problems are even more significant. The challenge is to find a solution, based on current digital systems, that supply information to operators in order to prevent events and take early actions to minimize or avoid any unexpected incident. In order to minimize human errors, Tecnatom has developed a specific set of displays, called 'High Performance Displays' (HPD) based on objectives, work environment limitation, and non-anticipated events. These displays are being design following the HFE standards: applying the NUREG-0711 elements to identify the information that they have to content and the NUREG-0700 to design the interface. The function of these displays is monitoring the plant situation in a way that the operators are able to recognize any potential event at a simple glance but controlling thousands of signals through pattern recognition, symmetries, regular polygons and colors. The objective of the HPD is to provide an integrated function minimizing the tunnel effect and reduce the workload during the diagnosis phase. Moreover, HPD enhance operators to take the right actions to anticipate the response to any event, even avoiding it, if possible This paper analyzes how HPD's have been design as well as presents the results of a preliminary validation of their capabilities. The preliminary validation is a reduce validation in a full scope simulator running two representative scenarios of abnormal and emergency operation. The analysis is focused on the first stage of the event, i.e. the detection phase, HPD's have been installed in a generic full scope simulator of a nuclear power plant, digital controlled. Both scenarios have been performed by two different operator crews. The first operator shift has worked with the conventional human machine interface, that is, displays similar to the ones currently installed in existing nuclear power plant with digital control. The second crew, reproduce the same scenarios but with HPD available. The studio collects preliminary data in order to compare results, draw lesson learned from the experience to improve HPD design, and find solutions in terms of reducing human errors, human workload, and increase the availability of the plant. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optimized, Budget-constrained Monitoring Well Placement Using DREAM

Defining the ideal suite of monitoring technologies to be deployed at a carbon capture and storage (CCS) site presents a challenge to project developers, financers, insurers, regulators and other stakeholders. The monitoring, verification, and accounting (MVA) toolkit offers a suite of technologies to monitor an extensive range of parameters across a wide span of spatial and temporal resolutions, each with their own degree of sensitivity to changes in the parameter being monitored. Understanding how best to optimize MVA budgets to minimize the time to leak detection could help to address issues around project risks, and in turn help support broad CCS deployment. This paper presents a case study demonstrating an application of the Designs for Risk Evaluation and Management (DREAM) tool using an ensemble of CO 2 leakage scenarios taken from a previous study on leakage impacts to groundwater. Impacts were assessed and monitored as a function of pH, total dissolved solids (TDS), and trace metal concentrations of arsenic (As), cadmium (Cd), chromium (Cr), and lead (Pb). Using output from the previous study, DREAM was used to optimize monitoring system designs based on variable sampling locations and parameters. The algorithm requires the user to define a finite budget to limit the number of monitoring wells and technologies deployed, and then iterates well placement and sensor type and location until it converges on the configuration with the lowest time to first detection of the leak averaged across all scenarios. To facilitate an understanding of the optimal number of sampling wells, DREAM was used to assess the marginal utility of additional sampling locations. Based on assumptions about monitoring costs and replacement costs of degraded water, the incremental cost of each additional sampling well can be compared against its marginal value in terms of avoided aquifer degradation. Applying this method, DREAM identified the most cost-effective ensemble with 14 monitoring locations. Here, while this preliminary study applied relatively simplistic cost and technology assumptions, it provides an exciting proof-of-concept for the application of DREAM to questions of cost-optimized MVA system design that are informed not only by site-specific costs and technology options, but also by reservoir simulation results developed during site characterization and operation.

47 OTHER INSTRUMENTATION↗

Accuracy evaluation of the optical surface monitoring system on EDGE linear accelerator in a phantom study

Frameless stereotactic radiosurgery (SRS) requires dedicated systems to monitor the patient position during the treatment to avoid target underdosage due to involuntary shift. The optical surface monitoring system (OSMS) is here evaluated in a phantom-based study. The new EDGE linear accelerator from Varian (Varian, Palo Alto, CA) integrates, for cranial lesions, the common cone beam computed tomography (CBCT) and kV-MV portal images to the optical surface monitoring system (OSMS), a device able to detect real-time patient's face movements in all 6 couch axes (vertical, longitudinal, lateral, rotation along the vertical axis, pitch, and roll). We have evaluated the OSMS imaging capability in checking the phantoms' position and monitoring its motion. With this aim, a home-made cranial phantom was developed to evaluate the OSMS accuracy in 4 different experiments: (1) comparison with CBCT in isocenter location, (2) capability to recognize predefined shifts up to 2° or 3 cm, (3) evaluation at different couch angles, (4) ability to properly reconstruct the surface when the linac gantry visually block one of the cameras. The OSMS system showed, with a phantom, to be accurate for positioning in respect to the CBCT imaging system with differences of 0.6 ± 0.3 mm for linear vector displacement, with a maximum rotational inaccuracy of 0.3°. OSMS presented an accuracy of 0.3 mm for displacement up to 1 cm and 1°, and 0.5 mm for larger displacements. Different couch angles (45° and 90°) induced a mean vector uncertainty < 0.4 mm. Coverage of 1 camera produced an uncertainty < 0.5 mm. Translations and rotations of a phantom can be accurately detect with the optical surface detector system.

61 RADIATION PROTECTION AND DOSIMETRY↗

Reconstruction of the Bi 2 Sr 2 CaCu 2 O 8 + δ Fermi surface

The effects of structural supermodulation with the period $λ ≈ 26$ Å along the $b$ axis of Bi 2 Sr 2 CaCu 2 O 8 + δ have been observed in photoemission studies from the early days as the presence of “diffraction replicas” of the intrinsic electronic structure. Although predicted to affect the electronic structure of the Cu-O plane, the influence of supermodulation potential on Cu-O electrons has never been observed in photoemission. In the present study, we clearly see the effects on the Bi 2 Sr 2 CaCu 2 O 8 + δ electronic structure; we observe a hybridization of the intrinsic bands with the supermodulation replica bands in the form of avoided crossings and a corresponding reconstruction of the Fermi surface. We estimate the hybridization gap, $2Δ_h ~ 25$ meV, in the slightly underdoped samples. The hybridization weakens with doping and the anticrossing can no longer be resolved in strongly overdoped samples. In contrast, the “shadow” replica, shifted by ($π , π$) , is found not to hybridize with the original bands within our detection limits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mt. Washington Icing Sensors Project (MWISP) Tests Technologies for Remote Sensing of Icing Conditions

The Mt. Washington Icing Sensors Project (MWISP) was an intensive multiagency field test of remote-sensing technologies and methods for in-flight icing detection. This field test was the first comprehensive examination of multiple remote-sensing technologies, and it included ground and airborne in situ measurements for comparison. The data gathered will allow the assessment of icing remote-sensing technologies leading toward the development of ground-based and airborne systems. These systems will provide flight crews with the information required for avoiding or exiting hazardous icing conditions.

Reehorst, Andrew L.↗

Magnetic Calorimeter Option for the Lynx X-Ray Microcalorimeter

One option for the detector technology to implement the Lynx x-ray microcalorimeter (LXM) focal plane arrays is the metallic magnetic calorimeter (MMC). Two-dimensional imaging arrays of MMCs measure the energy of x-ray photons by using a paramagnetic sensor to detect the temperature rise in a microfabricated x-ray absorber. While small arrays of MMCs have previously been demonstrated that have energy resolution better than the 3 eV requirement for LXM, we describe LXM prototype MMC arrays that have 55,800 x-ray pixels, thermally linked to 5688 sensors in “hydra” configurations, and that have sensor inductance increased to avoid signal loss from the stray inductance in the large-scale arrays when the detectors are read out with microwave superconducting quantum interference device multiplexers, and that use multilevel planarized superconducting wiring to provide low-inductance, low-crosstalk connections to each pixel. We describe the features of recently tested MMC prototype devices and simulations of expected performance in designs opti- mized for the three subarray types in LXM.

Stevenson, Thomas R.↗

Electrochemical Oxidation of Lithium Carbonate Generates Singlet Oxygen

Abstract Solid alkali metal carbonates are universal passivation layer components of intercalation battery materials and common side products in metal‐O 2 batteries, and are believed to form and decompose reversibly in metal‐O 2 /CO 2 cells. In these cathodes, Li 2 CO 3 decomposes to CO 2 when exposed to potentials above 3.8 V vs. Li/Li + . However, O 2 evolution, as would be expected according to the decomposition reaction 2 Li 2 CO 3 →4 Li + +4 e − +2 CO 2 +O 2 , is not detected. O atoms are thus unaccounted for, which was previously ascribed to unidentified parasitic reactions. Here, we show that highly reactive singlet oxygen ( 1 O 2 ) forms upon oxidizing Li 2 CO 3 in an aprotic electrolyte and therefore does not evolve as O 2 . These results have substantial implications for the long‐term cyclability of batteries: they underpin the importance of avoiding 1 O 2 in metal‐O 2 batteries, question the possibility of a reversible metal‐O 2 /CO 2 battery based on a carbonate discharge product, and help explain the interfacial reactivity of transition‐metal cathodes with residual Li 2 CO 3 .

25 ENERGY STORAGE↗

CD27{sup -}CD45{sup +} γδ T cells can be divided into two populations, CD27{sup -}CD45{sup int} and CD27{sup -}CD45{sup hi} with little proliferation potential

In addition to the majority of T cells which carry the αβ T cell receptor (TCR) for antigen, a distinct subset of about 1–5% of human peripheral blood T cells expressing the γδ TCR contributes to immune responses to infection, tissue damage and cancer. T cells with the Vδ2{sup +} TCR, usually paired with Vγ9, constitute the majority of these γδ T cells. Analogous to αβ T cells, they can be sorted into naive (CD27{sup +}CD45RA{sup +}), central memory (CD27{sup +}CD45RA{sup −}), effector memory (CD27{sup −}CD45RA{sup −}), and terminally-differentiated effector memory (CD27{sup −}CD45RA{sup +}) phenotypes. Here, we found that CD27{sup −}CD45RA{sup +} γδ T cells can be further divided into two populations based on the level of expression of CD45RA: CD27{sup −}CD45RA{sup int} and CD27{sup −}CD45RA{sup hi}. Those with the CD27{sup −}CD45RA{sup hi} phenotype lack extensive proliferative capacity, while those with the CD27{sup −}CD45RA{sup int} phenotype can be easily expanded by culture with zoledronate and IL-2. These CD27{sup -}CD45RA{sup hi} potentially exhausted γδ T cells were found predominantly in cancer patients but also in healthy subjects. We conclude that γδ T cells can be divided into at least 5 subsets enabling discrimination of γδ T cells with poor proliferative capacity. It was one of our goals to predict the feasibility of γδ T cell expansion to sufficient amounts for adoptive immunotherapy without the necessity for conducting small-scale culture tests. Fulfilling the ≥1.5% criterion for γδ T cells with phenotypes other than CD27{sup −}CD45RA{sup hi}, may help avoid small-scale culture testing and shorten the preparation period for adoptive γδ T cells by 10 days, which may be beneficial for patients with advanced cancer. - Highlights: • CD27{sup -}CD45{sup +} γδ T cells can be divided into two populations, CD27{sup -}CD45{sup hi} and CD27{sup -}CD45{sup int}. • CD27{sup -}CD45{sup hi} γδ T cells are detected in the peripheral blood of both cancer patients and healthy subjects. • CD27{sup -}CD45{sup hi} γδ T cells are distinct from so-called terminally differentiated CD27{sup -}CD45{sup int} γδ T cells. • CD27{sup -}CD45{sup hi} γδ T cells lack vigorous proliferative potential.

60 APPLIED LIFE SCIENCES↗

Induced photoemission from driven nonadiabatic dynamics in an avoided crossing system

When vibrational dynamics on an ionic state (large dipole moment) is coupled to that on a neutral state (small dipole moment) such as at an avoided crossing in the alkali halide system, the population transfer between the states cause oscillation of the molecular dipole, leading to dipole emission. Such dynamics may be driven by an external field. We study how the coupled wavepacket dynamics is affected by the parameters (intensity, frequency) of the driving field with the aim of making use of the photoemission as an alternative detection scheme of femtosecond and subfemtosecond vibrational and electronic dynamics or as a characteristic optical source.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Corrigendum to: Volatiles in Lunar Felsite Clasts: Impact-related Delivery of Hydrous Material to an Ancient Dry Lunar Crust

Hydrous components were measured in nominally anhydrous minerals, primarily lunar feldspar, and reported by Simon et al. (2020). In the feldspar mineral structure, these components are contained in the form of OH, H2O, and/or NH4+molecules (Johnson, 2006 and references therein). The advantages of the ion microprobe technique employed by Simon et al. (2020) include high spatial resolution, apparent insensitivity to crystal orientation, high precision, and low detection limits, but unlike complementary spectroscopy techniques (i.e., FTIR) they preclude assignment of molecular species(Mosenfelder et al. 2015). Convention in the lunar sample community is to report major elemental abundances and hydrogen as oxides, H2O is the oxide for H. This convention was not followed by Simon et al. (2020), but rather the measured hydrogen component was reported as H to avoid the appearance of assigning the measured hydrogen concentration to a specific molecular species. The speciation of the hydrous component in terrestrial feldspars varies predominantly between OH and H2O, with no apparent correlation to feldspar composition, and only slight correlation to volcanic-to-plutonic igneous rock type (Johnson and Rossman, 2004).

J. I. Simon↗

Semantic Complex Event Processing over End-to-End Data Flows

Emerging Complex Event Processing (CEP) applications in cyber physical systems like SmartPower Grids present novel challenges for end-to-end analysis over events, flowing from heterogeneous information sources to persistent knowledge repositories. CEP for these applications must support two distinctive features - easy specification patterns over diverse information streams, and integrated pattern detection over realtime and historical events. Existing work on CEP has been limited to relational query patterns, and engines that match events arriving after the query has been registered. We propose SCEPter, a semantic complex event processing framework which uniformly processes queries over continuous and archived events. SCEPteris built around an existing CEP engine with innovative support for semantic event pattern specification and allows their seamless detection over past, present and future events. Specifically, we describe a unified semantic query model that can operate over data flowing through event streams to event repositories. Compile-time and runtime semantic patterns are distinguished and addressed separately for efficiency. Query rewriting is examined and analyzed in the context of temporal boundaries that exist between event streams and their repository to avoid duplicate or missing results. The design and prototype implementation of SCEPterare analyzed using latency and throughput metrics for scenarios from the Smart Grid domain.

Zhou, Qunzhi↗

Hybridization-driven strong anharmonicity in Yb-filled skutterudites

A high-pressure study of the structural and thermal transport properties is carried out on one of the most efficient filled skutterudites, $\mathrm{Yb_{0.3}Co_4Sb_{12}}$, to understand the relatively low thermal conductivity behavior in this family. By combining x-ray diffraction and Raman scattering measurements, we detect a phase transition at around 12.4 GPa. The mode Grüneisen parameters of the observed phonon modes are obtained from the determined bulk modulus and the phonon frequency shifts with pressure. The strong anharmonicity in this material is demonstrated by the obtained large average Grüneisen parameter. We also find the depressed group velocity within the low-frequency range, the flat guest mode avoided crossing with the acoustic-phonon mode, and the significant contribution of optical phonons. The hybridization of the guest atom and host lattice and the related enhanced anharmonicity are suggested to account for the low lattice thermal conductivity in the studied system. So, these findings provide new insight into how phonon-phonon interactions lower lattice thermal conductivity in this important thermoelectric family.

36 MATERIALS SCIENCE↗

Turbocharging Quantum Tomography

Quantum tomography is used to characterize quantum operations implemented in quantum information processing (QIP) hardware. Traditionally, state tomography has been used to characterize the quantum state prepared in an initialization procedure, while quantum process tomography is used to characterize dynamical operations on a QIP system. As such, tomography is critical to the development of QIP hardware (since it is necessary both for debugging and validating as-built devices, and its results are used to influence the next generation of devices). But tomography suffers from several critical drawbacks. In this report, we present new research that resolves several of these flaws. We describe a new form of tomography called gate set tomography (GST), which unifies state and process tomography, avoids prior methods critical reliance on precalibrated operations that are not generally available, and can achieve unprecedented accuracies. We report on theory and experimental development of adaptive tomography protocols that achieve far higher fidelity in state reconstruction than non-adaptive methods. Finally, we present a new theoretical and experimental analysis of process tomography on multispin systems, and demonstrate how to more effectively detect and characterize quantum noise using carefully tailored ensembles of input states.

97 MATHEMATICS AND COMPUTING↗

Performance Modeling of a Pilot in a Free Flight Mode: A Free Flight Self-Separation Cancellations Due to the Requirement for Procedural Intervention - 1

Developments are being made that allow pilots to have more flexibility over the control of their aircraft. This new concept is called Free Flight. Free Flight strives to move the current air traffic system into an age where space technology is used to its fullest potential. Self-separation is one part of the Free Flight system. Self-separation provides pilots the opportunity to choose their own route to reach a specified destination provided that they maintain the 'minimum required separation distance between airplanes. In the event that pilots are unable to maintain separation, controllers will need to have the aircraft separation authority passed back to them. This situation is known as a procedural intervention point. This project attempted to examine and diagnose those particular situations in an effort to avoid reaching a procedural intervention point in the near future. Crews that reached procedural intervention points were compared with crews that made similar maneuver types in the same scenario, but did not reach procedural intervention points. Results showed that there were no significant differences between crews in a high-density acute angle flight conditions. However, significant differences in maneuver times, following the detection of an intruder aircraft and following the time the intruder aircraft came into view, were found in a low-density, acute angle scenario.

Ntuen, Celestine A.↗

Three-Dimensional High-Resolution Optical/X-Ray Stereoscopic Tracking Velocimetry

Measurement of three-dimensional (3-D) three-component velocity fields is of great importance in a variety of research and industrial applications for understanding materials processing, fluid physics, and strain/displacement measurements. The 3-D experiments in these fields most likely inhibit the use of conventional techniques, which are based only on planar and optically-transparent-field observation. Here, we briefly review the current status of 3-D diagnostics for motion/velocity detection, for both optical and x-ray systems. As an initial step for providing 3-D capabilities, we nave developed stereoscopic tracking velocimetry (STV) to measure 3-D flow/deformation through optical observation. The STV is advantageous in system simplicity, for continually observing 3- D phenomena in near real-time. In an effort to enhance the data processing through automation and to avoid the confusion in tracking numerous markers or particles, artificial neural networks are employed to incorporate human intelligence. Our initial optical investigations have proven the STV to be a very viable candidate for reliably measuring 3-D flow motions. With previous activities are focused on improving the processing efficiency, overall accuracy, and automation based on the optical system, the current efforts is directed to the concurrent expansion to the x-ray system for broader experimental applications.

Cha, Soyoung S.↗

Three-Dimensional High-Resolution Optical/X-Ray Stereoscopic Tracking Velocimetry

Measurement of three-dimensional (3-D) three-component velocity fields is of great importance in a variety of research and industrial applications for understanding materials processing, fluid physics, and strain/displacement measurements. The 3-D experiments in these fields most likely inhibit the use of conventional techniques, which are based only on planar and optically-transparent-field observation. Here, we briefly review the current status of 3-D diagnostics for motion/velocity detection, for both optical and x-ray systems. As an initial step for providing 3-D capabilities, we have developed stereoscopic tracking velocimetry (STV) to measure 3-D flow/deformation through optical observation. The STV is advantageous in system simplicity, for continually observing 3-D phenomena in near real-time. In an effort to enhance the data processing through automation and to avoid the confusion in tracking numerous markers or particles, artificial neural networks are employed to incorporate human intelligence. Our initial optical investigations have proven the STV to be a very viable candidate for reliably measuring 3-D flow motions. With previous activities focused on improving the processing efficiency, overall accuracy, and automation based on the optical system, the current efforts is directed to the concurrent expansion to the x-ray system for broader experimental applications.

Cha, Soyoung S.↗

High Speed Inlet Distortion Test for the X 59 Low Boom Flight Demonstrator in the NASA Glenn 8-by 6-Foot Supersonic Wind Tunnel

The Low Boom Flight Demonstrator (LBFD) High Speed Inlet Distortion Test (HSIDT) was completed in the NASA Glenn Research Center 8- by 6-Foot Supersonic Wind Tunnel (8x6 SWT) using a 9.5 percent-scale propulsion model of the X–59 to evaluate the compatibility and operability of the inlet and engine integration, focusing on the inlet dynamic effects at the engine face. The model was instrumented with steady-state and high-response total pressure probes in the conventional 40-point, equal-area array at the aerodynamic interface plane (AIP). The HSIDT successfully collected high response distortion data points for the X–59 from Mach = 0 to 1.55 for a wide range of angle-of-attack, angle-of-sideslip, and inlet flow. Analyses showed that inlet dynamic distortion and inlet planar wave (or inlet buzz) were within engine limits through most of the flight envelope. The test data showed that inlet planar wave was present at Mach 1.35 and above, at airflows between idle and maximum engine airflow. Engine throttle back will need to be limited when decelerating from supersonic flight Mach to avoid these regions of inlet planar wave. The HSIDT data was also used to compare the metric SUM40, computed from the 40 high-response AIP total pressure probes, with the metric SUM4, computed from four high-response inlet static pressure probes. The data showed that SUM4 can be used as a surrogate for SUM40 during the X–59 flight test to detect inlet planar wave.

X-59↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗