Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase Retrieval”

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.

872 records · Page 49

Elucidating proximity magnetism through polarized neutron reflectometry and machine learning

Polarized neutron reflectometry is a powerful technique to interrogate the structures of multilayered magnetic materials with depth sensitivity and nanometer resolution. However, reflectometry profiles often inhabit a complicated objective function landscape using traditional fitting methods, posing a significant challenge for parameter retrieval. In this work, we develop a data-driven framework to recover the sample parameters from polarized neutron reflectometry data with minimal user intervention. We train a variational autoencoder to map reflectometry profiles with moderate experimental noise to an interpretable, low-dimensional space from which sample parameters can be extracted with high resolution. We apply our method to recover the scattering length density profiles of the topological insulator–ferromagnetic insulator heterostructure Bi2Se3/EuS exhibiting proximity magnetism in good agreement with the results of conventional fitting. We further analyze a more challenging reflectometry profile of the topological insulator–antiferromagnet heterostructure (Bi,Sb)2Te3/Cr2O3 and identify possible interfacial proximity magnetism in this material. We anticipate that the framework developed here can be applied to resolve hidden interfacial phenomena in a broad range of layered systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reliable Determination of Pulses and Pulse-Shape Instability in Ultrashort Laser Pulse Trains Using Polarization-Gating and Transient-Grating Frequency-Resolved Optical Gating Using the RANA Approach

Devices that measure the presence of instability in the pulse shapes in trains of ultrashort laser pulses do not exist, so this task necessarily falls to pulse-measurement devices, like Frequency-Resolved Optical Gating (FROG) and its variations, which have proven to be a highly reliable class of techniques for measuring stable trains of ultrashort laser pulses. Fortunately, multi-shot versions of FROG have also been shown to sensitively distinguish trains of stable from those of unstable pulse shapes by displaying readily visible systematic discrepancies between the measured and retrieved traces in the presence of unstable pulse trains. However, the effects of pulse-shape instability and algorithm stagnation can be indistinguishable, so a never-stagnating algorithm—even when instability is present—is required and is generally important. In previous work, we demonstrated that our recently introduced Retrieved-Amplitude N-grid Algorithmic (RANA) approach produces highly reliable (100%) pulse-retrieval in the second-harmonic-generation (SHG) version of FROG for thousands of sample trains of pulses with stable pulse shapes. Further, it does so even for trains of unstable pulse shapes and thus both reliably distinguishes between the two cases and provides a rough measure of the degree of instability as well as a reasonable estimate of most typical pulse parameters. Here, we perform the analogous study for the polarization-gating (PG) and transient-grating (TG) versions of FROG, which are often used for higher-energy pulse trains. We conclude that PG and TG FROG, coupled with the RANA approach, also provide reliable indicators of pulse-shape instability. In addition, for PG and TG FROG, the RANA approach provides an even better estimate of a typical pulse in an unstable pulse train than SHG FROG does, even in cases of significant pulse-shape instability.

47 OTHER INSTRUMENTATION↗

Aerosol Observability and Predictability: From Research to Operations for Chemical Weather Forecasting. Lagrangian Displacement Ensembles for Aerosol Data Assimilation

A challenge common to many constituent data assimilation applications is the fact that one observes a much smaller fraction of the phase space that one wishes to estimate. For example, remotely sensed estimates of the column average concentrations are available, while one is faced with the problem of estimating 3D concentrations for initializing a prognostic model. This problem is exacerbated in the case of aerosols because the observable Aerosol Optical Depth (AOD) is not only a column integrated quantity, but it also sums over a large number of species (dust, sea-salt, carbonaceous and sulfate aerosols. An aerosol transport model when driven by high-resolution, state-of-the-art analysis of meteorological fields and realistic emissions can produce skillful forecasts even when no aerosol data is assimilated. The main task of aerosol data assimilation is to address the bias arising from inaccurate emissions, and Lagrangian misplacement of plumes induced by errors in the driving meteorological fields. As long as one decouples the meteorological and aerosol assimilation as we do here, the classic baroclinic growth of error is no longer the main order of business. We will describe an aerosol data assimilation scheme in which the analysis update step is conducted in observation space, using an adaptive maximum-likelihood scheme for estimating background errors in AOD space. This scheme includes e explicit sequential bias estimation as in Dee and da Silva. Unlikely existing aerosol data assimilation schemes we do not obtain analysis increments of the 3D concentrations by scaling the background profiles. Instead we explore the Lagrangian characteristics of the problem for generating local displacement ensembles. These high-resolution state-dependent ensembles are then used to parameterize the background errors and generate 3D aerosol increments. The algorithm has computational complexity running at a resolution of 1/4 degree, globally. We will present the result of assimilating AOD retrievals from MODIS (on both Aqua and TERRA satellites) from AERONET for validation. The impact on the GEOS-5 Aerosol Forecasting will be fully documented.

da Silva, Arlindo↗

BWR Spent Nuclear Fuel Acquisition and Testing to Support DOE-NE High Burnup Spent Fuel Data Project

The Office of Spent Fuel and Waste Disposition (SFWD) within the US Department of Energy (DOE) Office of Nuclear Energy (NE) established the Spent Fuel and Waste Science and Technology (SFWST) campaign to conduct research and development (R&D) activities related to the storage, transportation, and disposal of spent nuclear fuel (SNF) and high-level radioactive waste. The SFWST program was created within SFWD to address issues of extended or long-term SNF storage and transportation. Some near-term objectives of SFWST are to use a science-based, engineering-driven approach to: Support the enhancement of the technical bases to support the continued safe and secure dry storage of SNF for extended periods; Support the enhancement of the technical bases for retrieving SNF after extended dry storage; Support the enhancement of the technical bases for transporting high burnup (HBU) fuel and transporting low burnup fuel and HBU fuel after dry storage DOE-NE, in partnership with the Electric Power Research Institute, developed the High Burnup Spent Fuel Data Project to perform a large-scale demonstration and laboratory-scale testing of HBU pressurized water reactor (PWR) fuels (exceeding 45 gigawatt-days per metric ton of uranium [GWd/MTU]). Under this project, 25 sister rods—which are rods that have the same design, power histories, and other characteristics—were removed from assemblies at the North Anna Nuclear Power Station and sent to Oak Ridge National Laboratory (ORNL) in January 2016. ORNL performed detailed nondestructive examination (NDE) on all 25 rods. The NDE consisted of visual examinations, gamma and neutron scanning, profilometry and rod length measurements, and eddy current examinations. After completing the NDE, 10 of the sister rods were delivered to Pacific Northwest National Laboratory (PNNL) in a NAC International, Inc. legal-weight truck cask in September 2018 for destructive examination (DE). To date, SFWD work has focused on the PWR fuel that is part of the Sister Rod Test program. No boiling water reactor (BWR) fuel has been tested in the program, and the data needs that were identified for the PWR fuel have not been collected for BWR fuel. The goal to obtain six to nine BWR rods and test them at ORNL will support closing this important data gap. BWR fuel comprises approximately 56% of the total fuel assemblies currently in storage at nuclear power plants in the United States. BWR nuclear fuel and cladding designs and manufacturing are significantly different from PWRs. Differences include the following: BWR fuel pellets are larger than PWR pellets; Variations of Zircaloy-2 (including liners) are used instead of the Zircaloy-4 cladding materials used in PWRs; Clad manufacturing and stress-relief processes are different between PWRs and BWRs; The fuel rod dimensions are different because larger rod diameters and thicker cladding are used in BWRs; BWR fuel typically has lower internal rod pressures and sees vastly different operating conditions than PWR fuel (i.e., two-phase flow); BWR assemblies are “canned,” meaning each assembly is surrounded by a metal fuel channel; BWR cladding is often composed of an inner pure Zr liner that has widely different mechanical properties than the Zircaloy-2 alloy and exhibits a stronger affinity for hydrogen; The BWR SNF generally has more total hydrogen in the cladding/liner than typical PWR fuel; The construction of the PWR and BWR assemblies is vastly different; BWR rods are solidly attached to the assembly nozzles and experience a much different vibration and shock load than PWR rods, which are “floating” within a grid system attached to guide tubes, and the rods sit loosely on the bottom end plates. These numerous differences will affect the way the BWR SNF responds under dry storage preparation processes (e.g., vacuum drying) and during transportation. The results collected in the PWR experimental program must be compared with a subset of similar data collected on BWR SNF to establish a technical basis for whether the larger PWR database is sufficient to bound the BWR SNF end-of-life conditions as is currently assumed for several fuel/clad properties. Changes that occur in both fuel types at HBU could exacerbate any mechanical property differences. As the fuel burnup increases, several changes occur that might affect the performance of the fuel, cladding, and assembly hardware in storage and transportation. These changes include increased cladding corrosion layer thickness, increased cladding hydrogen content, increased cladding creep strains, increased fission gas release, and the formation of the HBU structure at the surface of the fuel pellets. The Nuclear Regulatory Commission (NRC) limits the current maximum rod-averaged burnup to 62 GWd/MTU due to these changes and the lack of data at higher burnups.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Microstructural Characterization of the Second High Fluence Baffle-Former Bolt Retrieved from a Westinghouse Two-loop Downflow Type PWR

As one of the pressurized water reactor (PWR) internal components, baffle-former bolts (BFBs) are subjected to significant mechanical stress and neutron irradiation from the reactor core during the plant operation. Over the long operation period, these conditions lead to potential degradation and reduced load-carrying capacity of the bolts. In support of evaluating long-term operational performance of materials used in core internal components, the Oak Ridge National Laboratory (ORNL), through the Department of Energy (DOE), Light Water Reactor Sustainability (LWRS) Program, Materials Research Pathway (MRP) has harvested two high fluence BFBs from a commercial Westinghouse two-loop downflow type PWR. The two bolts of interest, i.e. bolts # 4412 and 4416, were withdrawn from service in 2011 as part of a preventative replacement plan. No identification of cracking or potential damage was found for these bolts during their removal in 2011. However, the bolts required a lower torque for removal from the baffle structure than the original torque specified during installation. Irradiation displacement damage levels in the bolts range from 15 to 41 displacements per atom. The goal of this project is to perform detailed microstructural and mechanical property characterization of BFBs following in-service exposures. The information from these bolts will be integral to the LWRS program initiatives in evaluating end of life microstructure and properties. Furthermore, valuable data will be obtained that can be incorporated into model predictions of long-term irradiation behavior and compared to results obtained in high flux experimental reactor conditions. In this report, we present our latest study in FY22 on microstructural characterizations of the second high fluence baffle-former bolt, i.e., bolt # 4412. Analytical electron microscopy and atom probe tomography characterization were performed. The radiation-induced defects in the material add to the large wealth of knowledge for neutron-induced defects in 304/316 grades of stainless steels, specifically for radiation-induced precipitation after high fluence commercial PWR irradiation. The main findings are summarized as follows: 1) The cavity size was considerably larger in the bolt thread section than in the bolt head, with the bolt thread section having a bimodal distribution of cavities greater than ~6 nm in diameter and less than ~3 nm in diameter. The bolt head only had the small-sized cavities. In addition, there was a denuded zone of large cavities near grain boundaries in the thread section of the bolt. 2) Radiation-induced precipitation in the BFB #4412 was highly complex, with the volume fraction, size, and number density of Ni/Si and Cu-rich precipitates depending strongly on the radiation temperature/dose. In many cases, co-precipitates of adjoined clusters were found with Ni/Si-rich precipitates sandwiched between Cu-rich clusters and Mo/Cr/P-rich clusters. 3) Solute segregation out of solution was highest for most solutes in the thread section of the bolt #4412 with the exception of Cu, which experienced more separation out of solution into Curich clusters in the bolt head section. This highlights the difference in the mechanisms for precipitation of Ni/Si clusters, which have the Ni 3 Si phase composition, and precipitation of Cu-rich clusters. 4) There appear to be multiple simultaneous influences that affect the microstructural variation along the length of the bolt that overcomes the ~2X difference in irradiation dose between the bolt head and the bolt thread. The irradiation temperature, thermal/fast neutron ratio variation, potential strain gradient, and exposure to PWR coolant water that each section of the bolt sees may have more influence on the microstructural evolution than the total irradiation dose. The bolt thread and shank, with higher temperature, higher relative fast neutron flux, higher strain, and exposure to coolant but lower dose, underwent more enhanced cavity formation, precipitation, and solute segregation than the bolt head section.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Electrolytically Assisted Surface Decontamination (EASD{sup TM}) for POCO Operations - 20282

NNL in collaboration with C-Tech Innovation Ltd and Sellafield Ltd has been exploring innovative technologies to enable a significant reduction of radiological hazards within facilities during the Post Operational Clean Out (POCO) phase of a nuclear plant's life cycle. Reducing the hazard by effectively decontaminating plants in-situ during POCO delivers huge cost reductions for future decommissioning operations. These cost reduction benefits are realized by reducing the number and complexity of remote operations as well as lowering the long-term waste disposal costs. Whilst chemical decontamination can achieve the desired level of decontamination, applying aggressive chemical reagents is hazardous, potentially difficult to control and requires there to be complimentary effluent treatment and waste routes. This work has aimed to develop flexible and controllable decontamination processes which could be operated without the additional complexity and issues associated with chemical decontamination. The processes needed to be relatively fast and effective to minimize the time operators would spend in an active area. In addition, there was a driver to produce decontamination methodologies which generated a secondary waste compatible with current routes and which takes advantage of current waste capacity. The strategy therefore was to come up with a solution that works with, rather than against, the science and engineering behind process plants at Sellafield. Electrolytically Assisted Surface Decontamination (EASD{sup TM}) is an innovative electrochemical decontamination process (developed by NNL and C-Tech Innovation Ltd) that can remove activity from contaminated metal in very short time periods. Application of the patented electric waveform to the surface when treating contaminated metal with nitric acid has been shown to significantly enhance the decontamination performance. When compared to proposed baseline washout procedures at Sellafield, the only change to the process is the applied electrical waveform which is temporary and controllable. The current causes dissolution of the metal surface being decontaminated which leads to activity transferring from the plant item into the nitric acid effluent stream. This innovation has the potential to transform the POCO process and allow a nitric acid based washout approach to deliver POCO quickly and cost effectively within the existing infrastructure. Several devices have now been designed that incorporate this EASD{sup TM} technology and enable the decontamination process to be applied to a range items commonly identified as being contaminated during nuclear decommissioning programmes such as pipework, tanks and hotspots of walls/floors. Inactive and active laboratory-scale testing has been completed using both simulated and 'real-life' contaminated (low level waste) metallic items retrieved from different nuclear sites. Active trials demonstrated that contaminated items could be treated to free release levels within minutes. Engineering scale tests are currently being performed, with guidance from Sellafield's system engineers and plant managers, to provide the necessary re-assurance the designed devices could be successfully deployed in an on-plant scenario. The final stage of this development work is to carry out an active demonstration of the decontamination technology on Thorp at Sellafield. This paper aims to highlight progress made to date with the EASD{sup TM} technology and, more specifically, the development of an in-situ decontamination device for deployment within radioactive pipework. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Custom surface reflectance, shade mask, and equivalent water thickness maps for the Colorado Headwaters Ecological Spectroscopy Study (2025)

This dataset contains land surface reflectance estimates and additional derived products generated from NEON Imaging Spectrometer (NIS) data collected in the Upper Gunnison river basin during June and July of 2025. Data was collected over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). These products were derived from radiance and LiDAR data collected by the NEON Airborne Observation Platform (AOP) campaign funded by the Colorado Headwaters Ecological Spectroscopy Study (CHESS) (doi:10.15485/3017965). Products include per-pixel surface reflectance (rfl) and reflectance uncertainty (rfl_unc), observational data (obs), canopy equivalent water thickness (ewt), and shade masks. Atmospheric correction was performed per flightline using the ISOFIT (Imaging Spectrometer Optimal FITting) optimal estimation framework to estimate surface reflectance and the associated per-band reflectance uncertainty. Reflectance retrievals achieved a mean absolute error of 1.5% across diverse validation surfaces (see validation report.pdf). Equivalent water thickness was calculated from surface reflectance using the Beer–Lambert absorption of liquid water. Shade masks were generated based on the geometry between the sun angle, ground surface, and sensor at the time of flight. Data products are provided per-flightline and as mosaics for each domain. Flightline data products are provided as ENVI-formatted binary files (rfl, rfl_unc, ewt) and GeoTIFFs (shade). Reflectance and uncertainty mosaics are provided as tiled NetCDFs, while all other mosaicked products are provided as cloud-optimized GeoTIFFs. These formats are supported by common geospatial software (e.g., QGIS, ArcGIS, ENVI) and programmatic libraries in Python (e.g., rasterio, xarray, spectral, netCDF4) and R (e.g., terra, ncdf4). Processing workflows were designed to be equivalent to those used to generate the 2018 CHESS campaign airborne imaging spectroscopy data products (doi:10.15485/3013527). All outputs were co-registered to a common spatial grid to support time series analyses. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgment: Data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). Computational research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (80NM0018D0004) and was funded by EMIT Extended Mission Phase E Science.

2018 NEON and 2025 CHESS Campaigns↗

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↗