Search NASA⌕ Search

SEARCH · Search NASA

Results for “LOAD DISTRIBUTION - MONOCOQUE CONSTRUCTION”

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.

213 records · Page 12

A Multi-Model Framework for Assessing Long- and Short-Term Climate Influences on the Electric Grid

Climate change influences many aspects of the electric grid, but prior work and industry practices often ignore the potential effects of changing climate, or they only consider a single effect or individual effects in isolation. Challenges vary with each grid and include adapting to long-term trends such as changing temperature and precipitation or shorter-term events such as drought or storms that could increase in frequency or intensity. Here we present a multi-model framework designed to analyze the effects of long and short-term climate impacts in combination. This framework couples capacity expansion and production cost models with hydrologic models and future climate scenario data to analyze alternative climate and energy futures at high spatial, temporal, and process resolutions. Furthermore, we constructed and evaluated the results of a suite of simulated scenarios exploring climate impacts on capacity investment and stress-tested the resulting future infrastructures using hourly dispatch modeling under alternative drought and load conditions. We demonstrate the approach through a case study of the U.S. Western Interconnection, where climate impacts depend on interactions between temperature-induced load, water availability for hydropower, technology competitiveness, and demand flexibility. Changes in 2038 generating capacity range from -8.5-16.6 GW, and changes in 2038 transmission capacity range from -1-2 GW. Capacity increases are driven by higher load from higher temperatures, while capacity reductions can be achieved in scenarios with higher future hydropower availability and increased demand flexibility. Scenarios requiring additional capacity cost an additional $\$5$-$\$17$ billion (discounted) from 2018 to 2038; however, scenarios with capacity reductions cost $\$1$-$\$18$ billion less. Stress tests on four 2038 infrastructures demonstrated that the identified systems were able to serve at least 99.999% of load and 99.96% of reserves. However, drought and unexpected high-load conditions can result in reduced capacity to respond to contingency events we did not model. Although these results are system and scenario specific, they highlight the importance of considering multiple climate change impacts simultaneously in long-term planning efforts and demonstrate a multi-model, multiscale approach that can be flexibly applied to any system and set of climate change concerns.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Opportunities to Implement Solutions to Achieve Remedial Action Objectives in Consideration of Stakeholder Interests - 20365

A case study in the implementation of a soils excavation and removal project conducted at a Formerly Utilized Sites Remedial Action Program (FUSRAP) Maywood Superfund Site (FMSS) vicinity property; an active commercial business, in a densely populated area, with significant operational, technical and, logistical constraints that were expected to limit areas that could be remediated. The United States Army Corps of Engineers-led (USACE) team successfully navigated complex overlapping stakeholder interests to achieve a more efficient and complete removal of contaminated soils and debris while minimizing unnecessary excess costs to the Government. This paper centers on FUSRAP activities at the FMSS vicinity property located at 149-151 Maywood Avenue, Maywood, Bergen County, New Jersey. For most of its FUSRAP history, this ∼109,000 square meter (27-acre) property housed a now-demolished ∼26,000 square meter (6.5-acre) warehouse operated by Sears Logistics Services, a unit of the Sears Roebuck Company (Figure 2). Sears ended its property lease and vacated the warehouse in December 2016. Given that long history and for the purposes of this paper, the property will be referred to as the 'Sears property' or simply 'the property.' The Sears property is approximately 109,000 square meters (27-acres) and is currently zoned for commercial use. The property is bound to the north and northwest by 100 West Hunter Avenue (the Stepan Company), to the northeast by 205 Maywood Avenue, to the east by Maywood Avenue, to the south by 23 West Howcroft Road, and to the west by businesses on NJ Route 17 and the NJ Route 17 roadway. Until December 2016, the property housed a warehouse and distribution center operated by Sears Logistical Services. On-site structures were demolished by the property owner in 2017. The warehouse covered the north section of the property along the Stepan Company property line. A railroad spur from the adjacent property now known as the Maywood Interim Storage Site (MISS) ended at the northeast corner of the warehouse. Wetlands were located east of the warehouse; the rest of the property was covered with paved parking lots and grassed areas. During remediation at the property, opportunities frequently arose where the project team was able to coordinate effectively with stakeholders to sequence remediation activities to facilitate removal of otherwise inaccessible soils and identify opportunities for materials reuse when supported by residual radiological and chemical levels. Some specific opportunities include: - remediating individual truck loading dock bays to maintain tenant operations; - working along a busy highway and in a utility corridor containing a 30-inch high pressure gas main; - remediating Lodi Brook and associated wetlands requiring bypass pumping and other diversionary structures to maintain local community stormwater drainage; - protecting during remediation and supporting access to facilitate tenant/owner maintenance of critical fire protection and water supply system service lines buried in contaminated soils and necessary for safe warehouse operations; and - coordinating with stakeholders to ensure non-radiological contaminants of concern for the site (unaffiliated with FUSRAP) were addressed by the responsible party in a manner that maximized benefit to all parties while supporting an efficient overall site remediation program. Once the property was vacant and the warehouse structure and radiologically non-impacted above-grade structures were demolished, the project identified an opportunity to significantly reduce the volume of waste associated with the foundation of the former warehouse, a foundation suspected of being partially constructed in radiologically contaminated soils. The project developed and implemented supplemental radiological verification survey and sampling strategies based on radiological cross-contamination risk potential with process flow-charts and screening-level based decision points; a program that classified saw-cut sections of concrete based on observed residual surface radioactivity conditions using a combination of gamma sensitive sodium-iodide scans and beta sensitive Geiger-Muller direct measurements. The proposed approach was reviewed with stakeholders with feedback, including consideration of applicable State of New Jersey Site Remediation Program criteria and guidance before implementation. Once classified, additional sampling was performed at frequencies driven by screening. The sampling methods, stockpile sampling frequencies, criteria and, approach to data evaluation were developed in consideration of existing site cleanup standards, regional background ranges, State of New Jersey radiological remediation program guidance, and specific stakeholder input With the property remediation and survey efforts nearly completed, it is relevant to examine retrospectively the objectives, assumptions and limitations in the remedial design (i.e., what was planned) versus what was able to be accomplished through effective teaming. Benefits to other environmental remediation site programs include better understanding of how to work effectively with stakeholders to achieve win-win outcomes, and approaches to site remediation, waste minimization and materials reuse that were successful in their overall outcomes. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluation of an Interferometric Sensor for In-Space Detection of Gas Leaks

Space mission planning often involves long-term storage of volatile liquids or high-pressure gases. These may include cryogenic fuels and oxidizers, high-pressure gases, and life-support-critical consumables. The risk associated with the storage of fluids and gases in space systems has long been an issue and the ability to retain these fluids is often tied to mission success. A leak in the storage or distribution system can cause many different problems, including a simple, but mission endangering, loss of inventory or, in severe cases, unbalanced thrust loads on a flight vehicle. Cryogenic propellants are especially difficult to store, especially over a long duration. The propellant can boil off and be lost through the insulating walls of the tank or simple thermal cycling of the fittings, valves, and propellant feed lines may unseat seals allowing the fluid to escape. Current NASA missions call for long-duration in-space storage of propellants, oxidizers, and life support supplies. Leaks of a scale detectable through a pressure drop in the storage tank are often catastrophic and have long been the focus of ground-based mitigation efforts where redundant systems are often employed. However, there is presently no technology available for detecting and monitoring low-level, but still mission-endangering, gas leaks in space. Standard in-space gas detection methods either have a very limited pressure range over which they operate effectively or are limited to certain gases. Mass spectrometer systems are able to perform the detection tasks, but their size, mass and use of high voltage, which could potentially lead to an arc that ignites a combustible propellent, severely limit their usefulness in a space system. In this paper, we present results from testing of the light-based interferometric gas monitoring and leak detection sensor shown in Fig. 1. The output of the sensor is an interference fringe pattern that is a function of the gas density, and commensurate index of refraction, in the sample region. Changes in the density of gas cause the interference fringes to move across a photodiode detector, providing a temporal history of the leak. The sensor is fiber coupled and constructed from solid optics, allowing for placement almost anywhere on the spacecraft. It is also advantageous in that it consumes very little power and does not introduce an ignition source. Data are presented demonstrating the capability of the sensor to measure density variations in different gas species. In addition, the transient response of the sensor in vacuum is demonstrated. These data extend and improve upon the results previously presented by the authors in Ref. [1].

Polzin, Kurt A.↗

The Pursuit of K: Reflections on the Current State-of-the-Art in Stress Intensity Factor Solutions for Practical Aerospace Applications

The elastic stress intensity factor (SIF, commonly denoted as K) is the foundation of practical fracture mechanics (FM) analysis for aircraft structures. This single parameter describes the first-order effects of stress magnitude and distribution as well as the geometry of both structure/component and crack. Hence, the calculation of K is often the most significant step in fatigue analysis based on FM. This presentation will provide several reflections on the current state-of-the-art in SIF solution methods used for practical aerospace applications, including a brief historical perspective, descriptions of some recent and ongoing advances, and comments on some remaining challenges. Newman and Raju made significant early contributions to practical structural analysis by developing closed-form SIF equations for surface and corner cracks in simplified geometries, often based on empirical fits of finite element (FE) solutions. Those solutions (and others like them) were sometimes revised as new analyses were conducted or limitations discovered. The foundational solutions have exhibited striking longevity, despite the relatively "coarse" FE models employed many decades ago. However, in recent years, the accumulation of different generations of solutions for the same nominal geometry has led to some confusion (which solution is correct?), and steady increases in computational capabilities have facilitated the discovery of inaccuracies in some (not all!) of the legacy solutions. Some examples of problems and solutions are presented and discussed, including the challenge of maintaining consistency with legacy design applications. As computational power has increased, the prospect of calculating large numbers of SIF solutions for specific complex geometries with advanced numerical methods has grown more attractive. Fawaz and Andersson, for example, have been generating literally millions of new SIF solutions for different combinations of multiple cracks under simplified loading schemes using p-version FE methods. These data are invaluable, but questions remain about their practical use, because the tabular databases of key results needed to support practical life analysis can occupy gigabytes of storage for only a few classes of geometries. The prospect of using such advanced numerical methods to calculate in real time only those K solutions actually needed to support a specific crack growth analysis is also tempting, but the stark reality is that the computational cost is still so high that the approach is not practical except for specific, critical application problems. Some thoughts are offered about alternative paradigms. Compounding approaches are some of the earliest building blocks of SIF development for more complex geometries. These approaches are especially attractive because of their very low computational cost and their conceptual robustness; they are, in some ways, an intriguing contrast and complement to the brute-force numerical methods. In recent years, researchers at NRC-Canada have published remarkable results showing how compounding approaches can be used to generate accurate solutions for very difficult problems. Examples are provided of some successes--and some limitations--using this approach. These closed-form, tabulated numerical, and compounding approaches have typically been used for simple remote loading with simple load paths to the crack. However, many significant cracks occur in complex stress gradient fields. This is a job for weight function (WF) methods, where the arbitrary stress distribution on the crack plane in the corresponding uncracked body (typically determined using FE methods) is used to determine K. Several significant recent advances in WF methods and solutions are highlighted here. Fueled by advanced 3D numerical methods, many new solutions have been generated for classic geometries such as surface and corner cracks with wide ranges of geometrical validity. A new WF formulation has also be developed for part-through cracks considering the arbitrary stress gradients in all directions in the crack plane (so-called bivariant solutions). Basic WF methods have recently been combined with analytical expressions for crack plane stresses to develop a large family of accurate SIF solutions for corner, surface, and through cracks at internal or external notches with very wide ranges of shapes, sizes, acuities, and offsets. Finally, WF solutions are much faster than FE or boundary element solutions, but can still be much slower than simple closed-form solutions, especially for bivariant solutions that can require 2D numerical integration. Novel pre-integration and dynamic tabular methods have been developed that substantially increase the speed of these advanced WF solutions. The practical utility of advanced SIF methods, including both WF and direct numerical methods, is greatly enhanced if the FM life analysis can be directly and efficiently linked with digital models of the actual structure or component (e.g., FE models for stress analysis). Two recent advances of this type will be described. One approach directly interfaces the FM life analysis with the FE model of the uncracked component (including stress results). Through a powerful graphical user interface, simplified FM life models can be constructed (and visualized) directly on the component model, with the computer collecting the geometry and stress gradient information needed for the life calculation. An even more powerful paradigm uses expert logic to automatically build an optimum simple fracture model at any and every desired location in the component model, perform the life calculation, and even generate fatigue crack growth life contour maps, all with minimal user intervention. This paradigm has also been extended to the automatic calculation of fracture risk, considering uncertainty or variability in key input parameters such as initial crack size or location. Another new integrated approach links the engineering life analysis, the component model, and a 3D numerical fracture analysis built with the same component model to generate a table of SIF values at a specific location that can then be employed efficiently to perform the life calculation. Some attention must be given to verification and validation (V&V) issues and challenges: how good are these SIF solutions, how good is good enough, and does anyone believe the life answer? It is important to think critically about the different sources of error or uncertainty and to perform V&V in a hierarchal, building-block manner. Some accuracy issues for SIF solutions, for example, may actually involve independent material behavior issues, such as constraint loss effects for crack fronts near component surfaces, and can be a source of confusion. Recommendations are proposed for improved V&V approaches. This presentation will briefly but critically survey the range of issues and advances mentioned above, with a particular view towards assembling an integrated approach that combines different methods to create practical tools for real-world design and analysis problems. Examples will be selectively drawn from the recent literature, from recent enhancements in the NASGRO and DARWIN computer codes, and from previously unpublished research

CraigMcClung, R.↗

Software Testbed for Developing and Evaluating Integrated Autonomous Subsystems

To implement fault tolerant autonomy in future space systems, it will be necessary to integrate planning, adaptive control, and state estimation subsystems. However, integrating these subsystems is difficult, time-consuming, and error-prone. This paper describes Intelliface/ADAPT, a software testbed that helps researchers develop and test alternative strategies for integrating planning, execution, and diagnosis subsystems more quickly and easily. The testbed's architecture, graphical data displays, and implementations of the integrated subsystems support easy plug and play of alternate components to support research and development in fault-tolerant control of autonomous vehicles and operations support systems. Intelliface/ADAPT controls NASA's Advanced Diagnostics and Prognostics Testbed (ADAPT), which comprises batteries, electrical loads (fans, pumps, and lights), relays, circuit breakers, invertors, and sensors. During plan execution, an experimentor can inject faults into the ADAPT testbed by tripping circuit breakers, changing fan speed settings, and closing valves to restrict fluid flow. The diagnostic subsystem, based on NASA's Hybrid Diagnosis Engine (HyDE), detects and isolates these faults to determine the new state of the plant, ADAPT. Intelliface/ADAPT then updates its model of the ADAPT system's resources and determines whether the current plan can be executed using the reduced resources. If not, the planning subsystem generates a new plan that reschedules tasks, reconfigures ADAPT, and reassigns the use of ADAPT resources as needed to work around the fault. The resource model, planning domain model, and planning goals are expressed using NASA's Action Notation Modeling Language (ANML). Parts of the ANML model are generated automatically, and other parts are constructed by hand using the Planning Model Integrated Development Environment, a visual Eclipse-based IDE that accelerates ANML model development. Because native ANML planners are currently under development and not yet sufficiently capable, the ANML model is translated into the New Domain Definition Language (NDDL) and sent to NASA's EUROPA planning system for plan generation. The adaptive controller executes the new plan, using augmented, hierarchical finite state machines to select and sequence actions based on the state of the ADAPT system. Real-time sensor data, commands, and plans are displayed in information-dense arrays of timelines and graphs that zoom and scroll in unison. A dynamic schematic display uses color to show the real-time fault state and utilization of the system components and resources. An execution manager coordinates the activities of the other subsystems. The subsystems are integrated using the Internet Communications Engine (ICE). an object-oriented toolkit for building distributed applications.

Testbed↗

Evidence of Completion of Milestone 3: Experimental Testbed

Construction and deployment of the experimental testbed encountered scheduling delays, and completion was subsequently further delayed by curtailment of LANL operations due to the COVID-19 pandemic. Between March and October, LANL was in a state of “Limited Operations” such that only certain “mission-critical” work was performed with approximately 25% on-site staffing. In spite of this, much of the lab space preparation was executed and the five initial model house” units were constructed and installed during this period. Still, the pace of construction was impacted by difficulty in coordinating personnel and minimizing contact between workers. Currently, LANL is in the mode of “Normal Operations with Maximum Telework,” which continues to limit the availability of on-site personnel. Nevertheless, Milestone 3 is now complete. This document describes the experiment as it is deployed and provides status for each task for this milestone.

99 GENERAL AND MISCELLANEOUS↗

Hybrid Hydraulic-Electric Architecture (HHEA) for Mobile Machines (Final Report)

Traditionally, off-road mobile machines such as excavators and wheel loaders are primarily powered by hydraulics, and throttling valves are used to control their work circuits. In recent years, two general trends are towards more energy efficient systems and electrification. With electrification, both efficiency and control performance can be improved by the elimination of throttling losses and the use of high-bandwidth inverter control. Electrification is generally accomplished with Electro-hydraulic actuators (EHA) but they are limited to lower powered systems due to the high cost of electric machines capable of high power or high torque. This project proposes a new system architecture for off-road vehicles - Hybrid Hydraulic Electric Architecture (HHEA) to improve efficiency and control performance without requiring large electric machines. The widely applicable architecture combines hydraulic power and electric power in such a way that the majority of power is provided hydraulically while electric drives are used to modulate this power. In particular, HHEA utilizes multiple common pressure rails to transmit the majority of power and small electric machines to modulate the power. The energy-saving potential of the the HHEA has been validated for the work circuits of a variety of mobile machines, from small 5-ton excavators to medium sized 20-ton excavators and wheel loaders, and representative duty cycles to reduce energy input by 50-80% compared to the commercial state-of-art load-sensing systems. In addition, the corner power requirements of the electrical machines can be downsized by 85% compared to the EHA approach. Various tradeoff studies have also been conducted, including sensitivities to individual components performances, controllers, accumulator sizes, and variations of the system architecture etc. A control strategy has been developed to maintain or exceed the motion control precision compared to current systems. The motion control strategy consists of a nominal controller, based on a passivity-based backstepping design, and a transition controller, based on least-norn feedforward design. The nominal controller is used in between common pressure rail switches whereas the transition controller compensates for any disturbance that common pressure rail switchings inflict on the system. The control strategy has been experimentally validated on both a medium pressure (200bar) hardware-in-the-loop (HIL) testbed and a high pressure (300+bar) HIL testbed. An efficient and power-dense integrated electric-hydraulic machine consisting of an axial flux electric machine and a radial hydrostatic piston hydraulic machine has been designed, constructed and tested. The machine has an active material power density of 6.1kW/kg, a rated speed of 12500 RPM, and a design efficiency of 85%. This is among the highest power density electric machines using conventional materials. While the integrated machine was designed for modulating the hydraulic power within the HHEA, it can also be used in other applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Integrating AI Data Centers with the Power Grid

The rapid expansion of artificial intelligence (AI) has triggered an unprecedented surge in electricity demand, with US data center energy use projected to double or triple 2023 levels by 2028. This exponential growth places strain on grid infrastructure, which can hinder timely construction of desired computing capacity. To bridge this supply-demand gap, utilities and AI developers are increasingly turning to demand flexibility, a strategy that incentivizes shifting or reducing power use during peak periods of grid stress. Data centers are uniquely equipped for flexible operations due to their digital workloads, built-in redundancy, and onsite energy assets. This article outlines four primary mechanisms to enable data center flexibility: computational load flexibility (shifting tasks temporally or geographically), flexible use of core facility infrastructure adjustments, energy storage utilization, and onsite electricity generation. To encourage adoption, utilities are deploying new tariff designs, including voluntary interruptible service riders, mandated flexibility requirements, and streamlined interconnection processes for flexible loads. For the highly capitalized and rapidly growing AI industry, the primary motivators for embracing these strategies are expediting facility interconnection, satisfying emerging regulatory mandates, and mitigating community resistance. While demand flexibility cannot substitute the long-term need for new bulk power generation, it serves as an essential, immediate solution for enabling near-term deployment. By transforming data centers from grid stressors into stabilizing assets, flexible operations can ensure reliable grid integration, ease market pressures, and support a resilient power system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Results from Laboratory and Field Study of Thin Triple Pane Windows

Heat transfer through windows accounts for a significant percentage of a building’s energy use and adds substantially to the peak cooling load of a home. In recent years, improvements in glass manufacturing have enabled the use of a very thin central pane of glass similar to a cell phone screen to produce a thin triple-pane window, for finished insulated glass units (IGUs) with an overall thickness similar to standard double-pane windows. Because this highly insulating “thin triple” glass product can be incorporated into almost any existing window frame and can be fabricated at a modest added cost, the U.S. Department of Energy sponsored laboratory and field demonstration testing of thin triple-pane windows to validate thermal performance and installation requirements in real-life field settings. Thin triple pane windows were evaluated at the PNNL Lab Homes, a matched pair of manufactured houses located on PNNL’s campus in Richland, Washington and also at 16 different field study sites around the country. The experimental results include a comparison of heating, ventilation, and air-conditioning (HVAC) energy usage, condensation potential, occupant comfort, sound infiltration, and thermal performance. Field study data will be gathered through June of 2022; preliminary results are being shared in this paper. The lab houses are identical except that the reference house had standard double pane windows with assembly U-0.66 and the test house had thin triple pane windows with assembly U-0.19. Across the experimental test days, the daily HVAC savings ranged from 0.2 to 18.7 kWh (3%–18%) for the heating season and from 2.5 to 8.0 kWh (23%–41%) for the cooling season. The higher thermal performance of the thin triple-pane windows also reduced the condensation potential on the interior surface during winter months and provided more even distribution of temperatures throughout the home in comparison to the baseline. In addition to the added thermal performance, the thin triple-pane windows demonstrated significant acoustic benefits, reducing sound infiltration by 8 dB to 10 dB. For the field test portion of the project thin triple pane insulated glass units were produced by two different manufacturers, and then installed without modification into the ½” IGU pockets of the standard double pane frames of four other manufacturers. Field tests performed on existing homes in Washington, Montana, Colorado, and New York compared thin triple pane retrofits to original window conditions (before and after). Field tests at new construction sites in Minnesota, Michigan, and New York compared thin triple pane windows to commercially available solutions such as double pane or traditional triple pane (with a standard-thickness center pane). Field test work is ongoing, but preliminary results appear to follow the sound, surface temperature, and energy improvement results from the Lab Homes comparison. Additionally, reports from builders and installers indicate that thin triples require almost no added time or effort to install and look nearly identical to other windows, indicating the possibility of offering next-level performance with a product that requires very little modification to current production or installation practices, long considered a major barrier to technology uptake in the construction market.

energy efficiency, home retrofit, Windows, window ↗

Cloud Ozone Dust Imager (CODI): Investigation and Technical Plan - Volume 1

The Cloud Ozone Dust Imager (CODI) is proposed to investigate the current climatic balance of the Mars atmosphere, with particular emphasis on the important but poorly understood roles which dust and water ice aerosols play in this balance. The large atmospheric heating (20-50 K) resulting from global dust storms around Mars perihelion is well recognized. However, groundbased observations of Mars atmospheric temperatures, water vapor, and clouds since the Viking missions have identified a much colder, cloudier atmosphere around Mars aphelion that may prove as important as global dust storms in determining the interannual and long-term behavior of the Mars climate. The key climate issues CODI is designed to investigate are: 1) the degree to which non-linear interactions between atmospheric dust heating, water vapor saturation, and cloud nucleation influence the seasonal and interannual variability of the Mars atmosphere, and 2) whether the strong orbital forcing of atmospheric dust loading, temperatures and water vapor saturation determines the long-term balance of Mars water, as reflected in the north-south hemispheric asymmetries of atmospheric water vapor and polar water ice abundances. The CODI experiment will measure the daily, seasonal and (potentially) interannual variability of atmospheric dust and cloud opacities, and the key physical properties of these aerosols which determine their role in the climate cycles of Mars. CODI is a small (1.2 kg), fixed pointing camera, in which four wide-angle (+/- 70 deg) lenses illuminate fixed filters and CCD arrays. Simultaneous sky/surface imaging of Mars is obtained at an angular resolution of 0.28 deg/pixel for wavelengths of 255, 336, 502, and 673 nm (similar to Hubble Space Telescope filters). These wavelengths serve to measure atmospheric ozone (255 and 336 nm), discriminate ice and dust aerosols (336 and 673 nm), and construct color images (336, 502, and 673 nm). The CODI images are detected on four 512 x 512 pixel arrays, as partitioned on two 1024 x 1024 CCD's operated in frame transfer mode. The center of the CODI field-of-view is canted 40 deg from the zenith direction to obtain sky brightness measurements and a 20 deg surface field-of-view. Daily image observations will be conducted when the Sun is greater than or equal to 5 deg outside the edge of the CODI field-of-view, and twilight and nighttime imaging will obtained on a weekly basis. The 673 nm channel includes a polarizer wheel to obtain sky/surface polarimetry. A dust cover protects the entire lens assemblies of all four CODI channels. This opaque dust cover, which is normally opened for CODI imaging, includes a small fixed mirror and transparent window positioned above the 673 nm lens, to redirect the 673 nm field-of-view to the surface for descent imaging. Fixed pointing, internal data buffering, low operating power (2-4 W for less than or equal to 30 seconds), selective data transmission, and simple operational characteristics of the CODI experiment place minimum resource and operational demands on the Mars Surveyor 1998 lander. The CODI science goals are optimized for, but not restricted to, a low-latitude landing site (20 deg S-30 deg N). The primary CODI measurement objectives are the opacities, wave forms, particle properties (size, shape, and alignment), and heights of clouds; the opacities, particle properties, and vertical distribution of dust; and the opacity and vertical distribution of ozone. The variability of cloud, ozone, and dust opacities will be determined on diurnal, daily, and seasonal timescales. Wind velocities will be determined from cloud motions and wave characteristics; and the temporal variability of atmospheric water vapor, with limited altitude information, will be inferred from the CODI ozone observations. Secondary measurement objectives include limited descent imaging capability, surface uv-visible photometry and polarimetry, photochemistry, and meteorite infall rates.

Clancy, R. Todd↗

Grid-Integrated Production of Fischer-Tropsch Synfuels from Nuclear Power

Idaho National Laboratory (INL) investigates the relative economic profitability of an integrated energy system (IES) coupling an NPP with a synfuel production process at selected case study locations across the United States. In the synfuel IES, a high-temperature steam electrolysis (HTSE) plant is thermally and electrically coupled with an NPP to produce zero-carbon hydrogen. The synthetic fuel is produced from this H 2 combined with a CO 2 supply using the reverse water gas shift process followed by the Fischer-Tropsch (FT) reaction. This analysis considers a system in which the CO 2 is sourced from regional CO 2 emitters via the construction and operation of pipeline-based CO 2 supply networks. Locating the FT plant at the same site as the NPP and HTSE plants enables the NPP to provide zero-carbon heat and power to the HTSE plant and zero-carbon power to the FT plant as well as avoid the requirement for long-distance H 2 product transport from the HTSE plant to the FT plant. Hydrogen storage is used to enable the NPP to dispatch power to the electrical grid (instead of the HTSE plant) when grid demand increases, thus enabling the FT plant to continue to operate in a steady-state production mode. The ability to cease hydrogen production for several hours within each day enables the NPP to provide power to the grid to balance the electricity market during peak periods and maximize revenues for the nuclear synfuel IES. The FT process design considered has a 99% carbon conversion efficiency. The use of nuclear energy and nuclear energy-derived hydrogen enables synfuel production to achieve this high level of carbon utilization. Additionally, the life-cycle carbon emissions of the nuclear-based synfuel production process are very low, with WTW emissions of approximately 25 gCO 2 e/MJ, including steam credits (generated from FT process excess heat), and approximately 7 gCO 2 e/MJ, if steam credits are excluded. This compares favorably with the WTW emissions of 90.5 gCO 2 e/MJ for a compression-ignition, direct injection (CIDI) vehicle with a fuel economy of 31.6 miles per gallon gasoline equivalent (MPGGE), using low-sulfur diesel produced using conventional petroleum production and refining processes. Several NPPs in various regions of the U.S. are considered as case study analyses. Supply locations and transportation via pipeline of the CO 2 feedstock to the NPP site are analyzed through the National Energy Technology Laboratory (NETL) CO 2 Transport Cost model. The team finds that the amount of CO 2 generated by different sectors is sufficient for the synfuel production process at all locations considered. The CO 2 transportation costs are functions of the distance of the source to the NPP location, the CO 2 capture cost at the source, and the quantity of CO 2 transported. Historical electricity prices for the NPP case study locations are collected and analyzed. Monthly average prices, price range, and duration of negative-price periods vary among these locations. For each location, an auto-regressive moving average (ARMA) model is trained on historical electricity price data. ARMA validation is done to ensure the synthetic price distributions represent one of historical prices with high fidelity. Synthetic time series from these ARMA models are used in a coupled dispatch and system optimization in the Holistic Energy Resource Optimization Network (HERON) to compute the differential net present value (NPV) of the IES. The team finds that this econometric is positive, ranging from $14M–1.3bn (2020) depending on the location. The optimal synfuel IES configuration to obtain this increase in NPV often maximizes the size of the synfuel production process with regards to the size of the NPP. However, the team shows that the NPP still plays a stabilizing role for the grid: In periods of high prices and high loads, more electricity from the NPP is sent to the grid. A high variability of electricity prices and extreme maximum prices tend to drive up electricity production. While it requires significant investment, the synfuel IES could increase the economic profitability for the existing fleet of LWRs across the country while still maintaining the grid stabilizer role of NPPs. During its lifetime, the main costs for the nuclear synfuel IES are the carbon feedstock transportation costs, followed by the capital expenses (CAPEX) and operation and maintenance (O&M) costs while the revenue comes first from the IRA H 2 production tax credit (PTC) and then from the sales of synfuel products. The profitability of the synfuel IES is most sensitive to the value of the hydrogen PTC and the synfuel products as well as the cost of the carbon feedstock, highlighting the importance of governmental incentives regarding hydrogen, carbon emissions, and synfuel in driving the deployment of future nuclear synfuel IESs.

08 HYDROGEN↗

Direct Air Capture of CO 2 and Delivery to Photobioreactors for Algal Biofuel Production (Final Report)

A mobile DAC system was designed and constructed to pair with photobioreactors growing algae for biofuel production. The DAC system was designed as a versatile research system, rather than a compact production unit. The system was constructed and mounted on a mobile skid to facilitate transportation to the algae production site. Within the DAC system, CO 2 was captured using amine-loaded monoliths that allow for high CO 2 uptake with low pressure drop. The CO 2 is collected using a Global Thermostat patented temperature/vacuum swing adsorption (TVSA) process. Amine sorbents and process conditions were optimized to produce 10 to >20 g CO 2 .h-1. The stability of the amine sorbents was also studied, with sorbent modifications made to improve stability to degradation by oxidation. An Algenol-developed Spirulina strain (Arthrospira platensis AB2293) was selected as the production cyanobacterial strain. AB2293 cultured was inoculum for outdoor production following PBR installation by Algenol. The PBR system was composed of three independent PBRs, with each PBR composed of four hanging bags internally recirculated by a liquid turnover pump. The PBRs were operated outdoors in Atlanta, GA, and integrated with the DAC system. Algae were grown with similar productivity using DAC-CO 2 as algae grown using pure CO 2 obtained commercially (Airgas). Throughout the experimental duration, no discoloration was observed, and cellular morphology was consistent between the two experimental treatments. An LCA including lifecycle greenhouse gas emissions, full life cycle inventory of the Algenol system and the DAC system and integrated DAC+PBR system was developed. Lifecycle greenhouse gas emissions were calculated for capture of carbon dioxide using input from Global Thermostat and the National Renewable Energy Laboratory. Three scenarios for energy provision were evaluated: a natural gas combined heat and power system sized to meet the electricity requirement, a natural gas combined heat and power system sized to meet the process heat requirements, and a system without on-site power that procures the electricity from the grid. In all three cases, as expected, the major contributor to the emissions is the energy consumption associated with the desorption step of the DAC process. The LCA quantified the reduced potential energy and greenhouse gas emissions of heat and mass integration of DAC and Algenol compared to unintegrated DAC and Algenol systems. A life cycle assessment of the role of sorbent productivity and lifetime was also developed. The development of more robust, oxidation resistant DAC sorbents may enable small reductions in energy requirements and in lifecycle greenhouse gas emissions and other environmental impacts. NREL performed techno-economic analysis (TEA) to identify the integration scenario most likely to achieve a 15% cost reduction target versus the baseline. Heat and mass integration of DAC and the PBR is critical to minimizing the MFSP. The baseline case utilizes no heat and mass integration, and the DAC system provides 100% of the CO 2 required by the photobioreactors (20 tonnes/hr), operating for 12 hours/day capturing 40 tonnes CO 2 /operating hour. The minimum fuel selling price (MFSP) of ethanol calculated from the baseline case was $10.68/gal ethanol. This corresponds with a targeted MFSP of $9.07/gal ethanol (or 15% reduction). This target was achieved by integration Option 2a with the greatest cost reduction of 17.8% (or $8.78/gal) and integration Option 2b with a cost reduction of 16.4% (or $8.93/gal). Reductions in MFSP are attributed to two primary process considerations: (a) CO 2 storage at night reduces the capital expenses associated with DAC (i.e., increasing on-stream time); and (b) distributed DAC scenarios (DAC-PBR integration Options 2a and 2b) make use of boiler and DAC CHP flue gas CO 2 (free). Direct air capture on-stream time was one of the largest contributors to MFSP reduction.

09 BIOMASS FUELS↗

Ma'adim Vallis Estuarine Delta in Elysium Basin and Its Relevance as a Landing Site for Exobiology Exploration on Mars

The debouche of Ma'adim Vallis in the Elysium Basin generated a transitional transported sediment structure, which planimetric shape is controlled by the enclosing topography of a deep reentrant gulf of the Basin into the highland. We defined it as an estuarine delta. The location and the importance of this estuarine delta is supported by the theoretical model of graded profile constructed for Ma'adim Vallis, and by two approaches: (i) the reconstruction of Ma'adim Vallis downstream course from Gusev to Elysium Basin, and (ii) the survey of the sediment deposit in the alleged estuary. The longitudinal graded profile of Ma'adim Vallis finds its base-level in the Elysium Basin, at a about 1000 m elevation, which is in agreement with the observed Basin shoreline. This model is supported by observational evidence of flow between the northern rim of Gusev crater, and the Elysium Basin shoreline. This downstream course of Ma'adim Vallis can be divided into three hydrogeologic regions. into three hydrogeologic regions. (a) The first region is a flooded plain (Zephiria Mensae), consisting in chaotic terrain formed by highland rocks, and disintegrated lava of the western flank of Apollinaris. Morphologic indicators of the flood process are: (1) the sediment deposit over the Gusev crater northern rim that reflects the overspilling of the crater-lake water through a 40-km wide gap provided by an ancient impact crater, (2) the tear-drop shaped feature on the northeastern flank of Apollinaris Patera, and (3) the chaotic terrain that suggest the emergence of ground water generated by the seepage of the crater lake through high-permeable broken rampart material. This underground water circulation sustained by the hydrostatic pressure of the crater-lake has likely generated a hydrothermal system in the volcanic environment of Apollinaris Patera. The stratigraphy of the flooded area is identified as Hesperian age, with occurrences of Noachian hilly individual features, and as Amazonian flooded plain and chaotic material, (b) The second region is located on the western flank of Apollinaris Patera. It is surrounded by relics of deep valleys that suggest a former downstream course of Ma'adim Vallis. The geologic setting of this region (Lucus Planum) is interpreted to be an Amazonian formation composed by the middle and lower members of the Medusae Formation., c) The third region corresponds to the convergence of the west and east branches of Ma'adirn Vallis into a deep re-entrant wide gulf that penetrates about 100 km into the highland. This topographic depression is delineated by the 1000 in elevation contour. This gulf has formed an estuarine configuration centered at 3S/190W within the Elysium Basin. This configuration has favored the formation of a estuarine sedimentary delta, because of topographically controlled lateral migration. This estuarine structure is strongly dominated by the incoming supply of Ma'adim Vallis fluvial sediment extracted from Zephiria Mensae and Lucus Planum. The obtuse-angle geometry of the estuary increases the sedimentation rate, which is higher than in the course of the channel. The sediment deposition process is governed by the estuarine water circulation. The inflowing loaded fluvial water enters the estuary as a bottom current, and mixes with the relatively less-loaded water of the receiving basin. When they mixed. the inflowing fluvial material, and the landward basin circulating water generate an accumulation of highly-diversified estuarine deposit stratification. This accumulation of material is mostly centered in the transitional zone of the delta. The sediment trapping efficiency of the estuary is function of the energy balance between the inflowing fluvial water, and the ingoing basin current. The submergence of the delta by the rising of the water-level increases the estuary water-depth, and consequently the sediment entrapment is favored. The locus of sediment accumulation moves landward in the zone of inflowing fluvial water. This results in the rising of the channel base-level, thus in the increase of the length of the longitudinal graded-profile. The sediment deposit facies of the zone A shows a generally smooth surface. The longitudinal deposit is bordered by alluvial terraces that reflect the variations of the channel level. The waning of the Elysium Basin caused the erosion of the Basin estuarine zone by small channels, this episode being characterized by dissected tear-drop shaped mesa-like morphologies in the delta. Our estuarine delta model predicts a lithostratigraphic depositional sequence associated with the water submergence and the transgression of Elysium Basin. The thickness of the estuarine sediment corresponds to the Elysium Basin levels changes relatively to the bed floor of the estuary, The depositional sequence of Ma'adim Vallis are described: (1) a pro-current filled region (A), where fluvial are longitudinally accumulated by the inflowing water, (2) inverse current from Elysium Basin (B), where fluvial and lacustrine sediments are accumulated, and (3) zone of current equilibrium (C), where the sediments are distributed as a shoreline at the boundary of the estuarine delta. The estuary sedimentology dynamics collects and keeps the record of the geologic unit material crossed by Ma'adim Vallis, and those of the lakebed deposit of Elysium Basin. The predicted mixed stratigraphic sequence from fluvial and lacustrine sediment makes this site an exceptional environment to concentrate potential multi-origin biologic records. We envision four possible strategies to explore this sedimentologic record: (1) longitudinal surface and subsurface traverses in region A to investigate outcrop levees, (2) exploration of the mesa walls in region B, (3) deep drilling hole lodging of the sequential deposits in the zones A and B, and (4) surface and subsurface exploration of the shoreline delta. The expected results for each of these strategies are: (1) in the deepest layers of region A are predicted frequent and abundant coarse material, sandy lenses lamination grading downward from sand to cobbles. Volcanic debris from the Noachian crustal Plateau unit material, hydrothermal altered rocks, carbonates, Hesperian and possibly Amazonian volcanic material, from Apollinaris Patera, altered rocks and carbonates from Zephiria Mensae are expected. As a favorable environment for inception of life, possible biological records are expected in transported rock, (2) At the surface, and subsurface (<=100 m), large deposits sandy to silted material from Elysium paleolake basin mixed with fine-grained sediments from Ma'adim Vallis are expected mostly in the upstream part of region B, (3) on the shoreline of the estuarine delta, abundant fine material from Elysium paleolake basin (evaporites, carbonates), mostly Amazonian in age are expected. The Ma'adini estuary is a favorable landing site for all the above mentioned science aspects, and .for its location. The site lies near the equator, which is favorable for the rover solar power supply, and at 1000m elevation, which is a favorable configuration for the descent system braking. Another advantage is the extent of the area of high scientific interest (33,000 sq km), which provides a good ellipse, and potential long study traverses.

Grin, E. A.↗

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↗