Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “loading pattern optimization”

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

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

At least 145 records · Page 8

Datasets for Widespread Residential Space Heating Electrification in Texas

In this experiment, we explore long term patterns in electricity demand driven by the dual effects of full electrification of space heating in Texas (by adoption of electric heat pumps), and climate change. We use a predictive model of electricity demand, climate projections, and an open source nodal power system (DC Optimal Power Flow) model of the Electric Reliability Council of Texas (ERCOT) system. Heat pumps are a more energy efficient way of providing space heating and cooling in homes. We attempt to exhaustively investigate the impacts of full residential space heating electrification by adoption of heat pumps for the segment of Texas households that currently rely on fossil fuels (about 40%), while simultaneously incorporating climate change meteorological variables. We explore a range of scenarios of heat pump efficiency and climate uncertainty over a long period of future years (2020-2099). In total, the simulation experiment generates 1,280 simulation years of hourly data. We report and analyze results in form of impacts on residential load, total load, peak load, seasonality of peaking, and reliability measured by occurrence and frequency of loss of load events. While the experiment is for ERCOT, the insights and approach can be applied to other regions. The results from the analysis can inform system planners on a range of potential capacity requirements/ reliability implications and/or risks of full space heating electrification via the adoption of electric heat pumps, given the uncertainty in the scenarios/ climate futures. The dataset includes model output for residential, non residential and total load, and the results from the GO ERCOT model runs for 4 RCP Scenarios (RCP 4.5 Cooler, RCP 4.5 Hotter, RCP 8.5 Cooler, RCP 8.5 Hotter), 4 Heating electrification Scenarios (Base , Standard Efficiency HP, High Efficiency HP, Ultra-High Efficiency HP) over 80 years (2020-2099). The metrological variables at BA scale were weighted weighted using population projections consistent with the SSP3 scenario.

Climate Change↗

Chapter 2: Evaluating a Concentrating Solar Power Plantas an Extended-Duration Peaking Resource

We explore the ability of a concentrating solar power (CSP) plant with thermal energy storage (TES) to provide peaking capacity. We focus on future power systems, wherein net load patterns may be significantly different than they are today (e.g., due to higher renewable-energy penetrations). We examine 28 locations in the southwestern United States over an 18-year period. The hourly operation of the CSP plants are simulated to determine their potential to provide energy during an eight-hour peak-load window for each day up to 365 days per year. Our result shows that for the large majority of locations and years, CSP plants with certain configurations (i.e., in terms of solar field and TES sizes) can provide nearly 100% peak-load capacity. We examine also the amount of supplemental energy (e.g.,by using natural gas as a supplemental thermal-energy source) that would be required to ensure that a CSP plant could serve the eight highest-load hours of every day of the year. We find that in most cases, a CSP plant supplemented with natural gas would require less than 5% of the fuel that is used by a natural-gas fired power plant providing the same level of reliable capacity. A series of sensitivity analyses show that these results are robust to the number of peak-load hours and days that are considered and the configuration of the CSP plant.

capacity value↗

The power reliability event simulator tool (PRESTO): A novel approach to distribution system reliability analysis and applications

The growing interest in onsite solar photovoltaic and energy storage systems is partially motivated by customer concerns regarding grid reliability. However, accurately assessing the effectiveness of PVESS in mitigating these interruptions requires a comprehensive understanding of location-specific outage patterns and the ability to simulate realistic scenarios. To address the gap, we introduce the Power Reliability Event Simulation TOol (PRESTO), the first publicly available tool that simulates location-specific power interruptions at the county level. PRESTO allows for a more realistic assessment of system reliability by considering the unpredictability and location-specific patterns of power interruptions. We applied PRESTO in a case study of a single-family home across three U.S. counties, examining the performance of a solar photovoltaic system with 10kWh of battery storage during short-duration power interruptions. Our findings show that this system reliably met 93% of energy demand for essential non-heating and cooling loads, fully serving these loads in 84% of events, despite the constraints of daily time-of-use bill management which limits the battery's state-of-charge reserve. However, when heating and cooling loads were included, system performance decreased significantly, with only 70% of demand met and full service in 43% of events. These results highlight the challenges of using solar photovoltaic and energy storage systems for short-duration outages, emphasizing the need to consider factors like battery size and grid charging strategies to improve reliability. Our study demonstrates the practical applications of PRESTO, providing valuable insights into potential mitigation strategies including grid charging and optimizing battery size.

14 SOLAR ENERGY↗

Design and optimization of flexible decoupled high-temperature gas-cooled reactor plants with thermal energy storage

Advanced nuclear power plants are well-positioned for future zero-carbon grids, however, the need for flexible power generation will be required over the traditional emphasis on baseload generation for meeting historical demands. To achieve such flexibility, this work examines viable configurations for coupling nuclear energy production with thermal energy storage. Previous designs on nuclear-thermal energy storage configurations for advanced reactor designs, which utilized reactor steam as the heat source for charging the thermal energy storage, are restricted by the heat diversion ratio and efficiency losses, thus their impacts can be limited. In this context, this study proposes configurations for fully decoupling the nuclear reactor from the power cycle and positioning the storage as an intermediate loop, thereby achieving an unconstrained heat diversion ratio and improved efficiency. Compared with a standard high-temperature gas-cooled reactor’s power cycle, steady-state thermodynamic modeling and dispatch optimizations quantify the benefits of a steam reheat cycle within the fully decoupled thermal energy system to separate the plant cycle from the high-pressure primary side. These benefits are further detailed, compatible with required high-temperature and high-pressure conditions, through (1) open-source dynamic transient models that examine the impact of off-design operation on the systems, (2) the investigation of components design and costing and finally (3) sizing and dispatch optimization. The fully-decoupled design achieves a cycle efficiency of 43.1%, an enhancement over the vendor’s standard efficiency of 42.2% (Xe-100 design). Here, the proposed design offers strengthened physical barriers from the nuclear island as well as superior operational flexibility and power boosting. Dispatch optimization and market analysis reveal that thermal energy storage size is highly dependent on the peak patterns of electricity prices and the minimum generation level constraint imposed on the balance of plant. Evaluation of off-design operation demonstrates that the full decoupling design with the suggested fail-safe control mechanisms ensures a minimal impact on reactor parameters, even during rapid power ramping.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Probing the atmospheric boundary layer with integrated remote-sensing platforms during the American WAKE ExperimeNt (AWAKEN) campaign

The American WAKE ExperimeNt (AWAKEN) collaboration is an observational-based field campaign in northern Oklahoma intended to analyze the potential influence of onshore wind farms and their collective wakes on wind power production, turbine structural loads, and on the atmospheric boundary layer (ABL). Focusing on the ABL effects, the University of Oklahoma and the Lawrence Livermore National Laboratory collected continuous high-resolution kinematic and thermodynamic profile measurements during 2022 and Summer 2023. The deployment strategy for these campaigns is detailed first, followed by an initial comparison of data from two sites in the AWAKEN domain: a near-farm site to examine collective wake impacts on the ABL, and a far-field site remaining outside the wind farm-waked region. Here, we summarize the datasets available and demonstrate the benefits of these observations and multiple value-added products (VAPs) for investigation of ABL features observed during AWAKEN. We also highlight examples of preliminary analyses, including ABL height detection and nocturnal low-level jet examination, which are produced using novel VAPs based on optimal estimation to retrieve deeper Doppler lidar wind profiles than previously resolved, along with their uncertainty. By including the near-farm and far-field site in these analyses, we identified a pattern of stronger lower-atmospheric mixing at the near-farm site than the far-field site, motivating deeper investigation into the relationship between wind farms and general ABL characteristics. Future analysis will delve deeper into this relationship by examining other ABL characteristics, such as atmospheric stability and convection.

17 WIND ENERGY↗

Self-aligned concentrating immersion-lens arrays for patterning and efficiency recovery in scaffold-reinforced perovskite solar cells

In contrast to existing silicon, CIGS, and multi-junction cells that exhibit remarkable mechanical durability, perovskite solar cells have been shown to delaminate when subjected to the mechanical loads that occur during processing and field exposures, limiting their potential as a reliable solar technology. Mechanical reinforcement is thus essential for durable and reliable perovskite technologies. One current strategy to overcome the thermomechanical fragility of perovskite solar cells is to extrinsically shield them by introducing reinforcing scaffolds which partition the cell into many distinct microcells, but improvements in mechanical stability coincide with a reduced device efficiency due to parasitic absorption by the scaffold. We address this reduced efficiency by integrating concentrating immersion-lens arrays (CILAs) into scaffold-reinforced solar cells. These scaffolds are lithographically patterned by a maskless process, in which ultraviolet light is used to pattern the scaffolds through the CILAs, ensuring self-alignment and optical contact with the microcells. Perovskite devices are deposited into the microcells, and during operation, light is concentrated into the microcells and away from the insulating scaffolds, resulting in mechanically resilient solar cells with efficiencies comparable to planar devices. The CILAs also exhibit passive tracking of incident light—verified experimentally and with ray-tracing simulations—which is critical for optimal power output as the sun moves across the sky during the day. The simplicity of the fabrication process and the efficiency of the resulting scaffolded devices show that a lens-integrated scaffold-reinforced structure is a potential pathway to high-performance, robust, commercially viable perovskite solar cells.

14 SOLAR ENERGY↗

From Cell to System: Accelerated hpc Simulations of BESS Aging under Frequency Regulation and Arbitrage use cases

Lithium-ion battery energy storage systems (BESS) packs have emerged as a leading solution for grid-scale energy storage, enhancing resiliency and balancing load fluctuations. Yet, experimental characterization of large-format LIB packs-particularly to assess performance and degradation over hundreds of cycles - demands substantial hardware investment and multi-year testing campaigns. In this work, we couple a hierarchical, physics-based modeling framework agnostic to electrode chemistries with high-performance computing to accelerate systems level evaluation by upto two orders of magnitude. Building on the open-source liionpack platform, we implement cell, module, and pack-scale electrochemical models enriched with mechanistic aging mechanisms and deploy them on an HPC cluster to simulate 150−200kWh systems over 500 - 1,000 cycles with in days. We subject these virtual B ESS to both constant-current cycling and realistic grid service profiles spanning frequency regulation, ramp-rate support, and energy arbitrage-and quantify the resulting degradation patterns. Our results reveal that localized cell aging can induce substantial nonuniformity at module and pack levels, with service-specific cycling protocols driving distinct aging modes. This rapid, multiscale modeling approach provides a powerful design-space exploration tool for optimizing electrical architecture, control strategies, and operational schedules to prolong pack lifetime and lower total cost of ownership.

Ayalasomayajula, Surya [ORNL] (ORCID:0009000860788↗

Analysis of a Dynamic Multi-Track Airway Concept for Air Traffic Management

The Dynamic Multi-track Airways (DMA) Concept for Air Traffic Management (ATM) proposes a network of high-altitude airways constructed of multiple, closely spaced, parallel tracks designed to increase en-route capacity in high-demand airspace corridors. Segregated from non-airway operations, these multi-track airways establish high-priority traffic flow corridors along optimal routes between major terminal areas throughout the National Airspace System (NAS). Air traffic controllers transition aircraft equipped for DMA operations to DMA entry points, the aircraft use autonomous control of airspeed to fly the continuous-airspace airway and achieve an economic benefit, and controllers then transition the aircraft from the DMA exit to the terminal area. Aircraft authority within the DMA includes responsibility for spacing and/or separation from other DMA aircraft. The DMA controller is responsible for coordinating the entry and exit of traffic to and from the DMA and for traffic flow management (TFM), including adjusting DMA routing on a daily basis to account for predicted weather and wind patterns and re-routing DMAs in real time to accommodate unpredicted weather changes. However, the DMA controller is not responsible for monitoring the DMA for traffic separation. This report defines the mature state concept, explores its feasibility and performance, and identifies potential benefits. The report also discusses (a) an analysis of a single DMA, which was modeled within the NAS to assess capacity and determine the impact of a single DMA on regional sector loads and conflict potential; (b) a demand analysis, which was conducted to determine likely city-pair candidates for a nationwide DMA network and to determine the expected demand fraction; (c) two track configurations, which were modeled and analyzed for their operational characteristic; (d) software-prototype airborne capabilities developed for DMA operations research; (e) a feasibility analysis of key attributes in the concept design; (f) a near-term, transitional application of the DMA concept as a proving ground for new airborne technologies; and (g) conclusions. The analysis indicates that the operational feasibility of a national DMA network faces significant challenges, especially for interactions between DMAs and between DMA and non-DMA traffic. Provided these issues are resolved, sectors near DMAs could experience significant local capacity benefits.

Wing, David J.↗

Simulation of the Multi-Wake Evolution of Two Sandia National Labs/National Rotor Testbed Turbines Operating in a Tandem Layout

The future of wind power systems deployment is in the form of wind farms comprised of scores of such large turbines, most likely at offshore locations. Individual turbines have grown in span from a few tens of meters to today’s large turbines with rotor diameters that dwarf even the largest commercial aircraft. These massive dynamical systems present unique challenges at scales unparalleled in prior applications of wind science research. Fundamental to this effort is the understanding of the wind turbine wake and its evolution. Furthermore, the optimization of the entire wind farm depends on the evolution of the wakes of different turbines and their interactions within the wind farm. In this article, we use the capabilities of the Common ODE Framework (CODEF) model for the analysis of the effects of wake–rotor and wake-to-wake interactions between two turbines situated in a tandem layout fully and partially aligned with the incoming wind. These experiments were conducted in the context of a research project supported by the National Rotor Testbed (NRT) program of Sandia National Labs (SNL). Results are presented for a layout which emulates the turbine interspace and relative turbine emplacement found at SNL’s Scaled Wind Technologies Facility (SWiFT), located in Lubbock, Texas. The evolution of the twin-wake interaction generates a very rich series of secondary transitions in the vortex structure of the combined wake. These ultimately affect the wake’s axial velocity patterns, altering the position, number, intensity, and shape of localized velocity-deficit zones in the wake’s cross-section. This complex distribution of axial velocity patterns has the capacity to substantially affect the power output, peak loads, fatigue damage, and aeroelastic stability of turbines located in subsequent rows downstream on the farm.

Baruah, Apurva (ORCID:0000000252354068)↗

The Toughness of Interlocking Metasurfaces

Interlocking metasurfaces (ILMs) are arrays of autogenous latching unit cells patterned across a surface. These create structural joints similar to bioinspired suture joints but patterned over a 2D surface rather than a 1D seam. This enables ILMs to be an alternative to conventional joining technologies such as bolts, welds, and adhesives. However, compared to conventional joining methods, relatively little is known of the engineering considerations for designing structural ILMs. Herein, the interfacial toughness of an archetypal ILM is examined for the first time. Under the conditions studied here, the ILM is substantially tougher than the material from which it is made, in this case, exhibiting up to a 50% increase in interfacial crack initiation energy over the solid base material, a photocured 3D‐printed polymer. Through experimental tests using in‐situ digital image correlation along with complementary computational analyses, the mechanism of toughening in the ILM structure and the origins of toughness anisotropy are revealed. The increase in toughness is associated with cross‐cell interactions, that is, load‐sharing across unit cells, which give rise to a finite process zone length with different effective material properties. In this way, ILM toughening is analogous to crack blunting in ductile materials or fiber bridging in composites; yet here, the ILM is composed of a single‐phase base material and so the architected toughening is geometric in nature and hence amenable to future topological optimization.

fracture↗

Flexible Geothermal Power Generation utilizing Geologic Thermal Energy Storage (Final Seedling Project Report)

This report concludes that there is a cost-effective strategy for seasonal storage of heat that will provide firm, but dispatchable, electrical generating capacity in times when other renewable energy is not available to meet demand. Deployment of the technology appears to require no new technology, but instead combines solar, geothermal, and conventional oil and gas drilling technologies in a novel way. The study basis is the use of sedimentary geologic formations as a medium for thermal energy storage (TES), specifically for heat collected in concentrating solar collectors. The study identifies methodologies that could be used to transport this heat into and out of the subsurface in order to produce dispatchable electrical power, and reports on initial optimization results. The GeoTES system (heat input, storage, heat recovery, and heat to electric conversion) described in this analysis has the potential to provide a unique pathway for increasing the grid penetration of renewable energy in large blocks of power and across many states and regions. Further, the system can be used both to meet the nation’s flexible energy needs while also improving grid stability and reliability. The present study evaluated the use of a large number of dedicated wells to store and recover the heat, essentially creating a synthetic geothermal reservoir. The use of sedimentary geology allows the wells to be drilled at low cost. Dedicated hot and cold wells are used, arranged in a five-spot well pattern with each hot or cold well surrounded at an appropriate distance by the opposite type of well. In large numbers this becomes alternating rows of hot and cold wells. Each hot and cold well is operated using a push-pull strategy. This configuration provides the ability to immediately recover stored hot fluid from a GeoTES reservoir, or to store the heat over many months for recovery at low loss when needed. This is a practical approach for managing the system’s fluid inventory and reducing parasitic load. The production and injection power requirements are reduced because the rows of wells operating in “push” mode provide help to the wells operating in “pull” mode, and vice-versa. Initial charging of a GeoTES system increases the heat recovery temperature. Increasing the duration of the charging period decreases the magnitude of the temperature fluctuations that occur following prolonged system operation. Because of direct contact of the heated water with the reservoir formation, the production of both hot water and steam from the TES, and the temperature ranges of the recovered fluid (190 – 230°C or 375 – 445°F), conventional geothermal power cycles were used to convert the stored heat to electricity. A power cycle configuration for the GeoTES system was selected following a screening study of a number of flash, and flash/binary hybrid options. This analysis concluded that, of the configurations evaluated, a dual-stage flash steam cycle provides the lowest capital costs per unit net power generation with an acceptable hot brine inlet fluid flow rate. The evaluation included the power plant cost estimate, the cost and number of wells and the associated parasitic loads. Annual power generation performance was simulated to evaluate capacity factor and LCOE. The LCOE calculated for the inherently high capacity GeoTES system was $\$$0.13/kWhe. This value was calculated for the case where the solar thermal collector was sized in such a way that the solar collectors permitted an annual power plant capacity factor of up to 97%. The power cycle was able to provide power to the grid every night of the year, and flexible base-load power during the winter, if needed. This LCOE value compares favorably with reported values for solar photovoltaic plus battery energy storage (PV+BES) systems in the open literature, i.e. $\$$0.148/kWhe for a PV+BES system with 4 hours of electrochemical battery energy storage capacity. Addition of battery energy storage with more hours of storage would further increase PV+BES system LCOE and increase the separation between GeoTES and PV+BES. A GeoTES system would therefore provide superior economics for high capacity and long duration solar energy storage.

14 SOLAR ENERGY↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand while simultaneously maintaining air travel as one of the safest forms of transportation. One of the reasons for this success is the ability of the air traffic control system and the operators to adapt and accommodate to situations that routinely disrupt normal 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, complexity increases. This is because, under these conditions humans are required to make tactical decisions in response to external factors. This results in a departure from the original strategic plan where operations would be more efficiently managed. Human operators manage airspace complexity under rigid regulations but in a constantly changing environment. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. Some prior studies devised airspace complexity metrics in commercial aviation and related these metrics to controller workload. 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. 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 ours that identifies such contributing factors or precursor patterns.

Precursor↗

Evaluation of Graph Analytics Frameworks Using the GAP Benchmark Suite

The analysis of connected data is an increasingly important application in high-performance computing. Such analyses can reveal fraudulent patterns in financial transactions, optimize telecommunications networks, predict information flow in social networks, etc. However, the landscape of graph analytics is highly diverse. Graph algorithms stress processor architectures differently, and no one graph can represent all topologies. Consequently, no single approach or framework is expected to be optimal for all graph analytics problems. To help make sense of this diverse landscape, we evaluated four approaches to graph analytics: GraphBLAS, Galois, BGL17, GraphIt; and compare them against hand-tuned implementations that take advantage of hardware features on our test platform. Graph- BLAS formulates graph analytics as sparse linear algebra. Galois provides syntactic constructs for data parallelism over irregular data structures. BGL17 is a generic C++ template library for implementing graph algorithms. GraphIt provides a domain- specific language to describe and optimize graph algorithms. We use the GAP Benchmark Suite to establish baseline performance and guide the side-by-side evaluation of each framework. GAP consists of 30 tests: six graph analytics algorithms (breadth- first search, single-source shortest path, PageRank, betweenness centrality, connected components, and triangle counting) run on five graphs, each with different topological characteristics (e.g., high diameter, skewed degree distribution, high average degree). High-performance reference implementations are included for each benchmark algorithm. Because a graph can be loaded into memory a number of ways (e.g., flat file on disk, compressed sparse format, data frames, retrieved from SQL or NoSQL databases), our evaluation focused on computational performance rather than I/O. Our results show the relative strengths of each framework.

Graph algorithms, Benchmarking, shared-memory prog↗

Numerical Experiments of Counterflowiing Jet Effects on Supersonic Slender-Body Configurations

Previous studies have demonstrated that the use of counterflowing jets can greatly reduce the drag and heat loads on blunt-body geometries, especially when the long penetration mode jet condition can be established. Previously, the authors had done some preliminary numerical studies to determine the ability to establish long penetration mode jets on a typical Mach 1.6 slender configuration, and study its impact on the boom signature. The results indicated that a jet with a longer penetration length was required to achieve any impact on the boom signature of a typical Mach 1.6 slender configuration. This paper focuses on an in-depth parametric study, done using the space-time conservation element solution element Navier-Stokes flow solver, for investigating the effect of various counterflowing jet conditions/configurations on two supersonic slender-body models (cone-cylinder and quartic body of revolution). The study is aimed at gaining a better understanding of the relationship between the shock penetration length and reduction of drag and boom signature for these two supersonic slender-body configurations. Different jet flow rates, Mach numbers, nozzle jet exit diameters and jet-to-base diameter ratios were examined. The results show the characteristics of a short-to-long-to-short penetration-mode pattern with the increase of jet mass flow rates, observed across various counterflowing jet nozzle configurations. Though the optimal shock penetration length for potential boom-signature mitigation is tied to the long penetration mode, it often results in a very unsteady flow and leads to large oscillations of surface pressure and drag. Furthermore, depending on the geometry of the slender body, longer jet penetration did not always result in maximum drag reduction. For the quartic geometry, the maximum drag reduction corresponds well to the longest shock penetration length, while this was not the case for the cone-cylinder-as the geometry was already optimized for drag. Numerical results and assessments obtained from this parametric study along with the recommendation for future implementation of counterflowing jets as a means for drag and noise reduction are detailed in this paper.

Venkatachari, Balaji Shankar↗

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↗

Quantifying the Aerosol Semi-Direct Effect in the NASA GEOS-5 AGCM

Aerosols such as black carbon, dust, and some organic carbon species both scatter and absorb incoming solar radiation. This direct aerosol radiative forcing (DARF) redistributes solar energy both by cooling the surface and warming the atmosphere. As a result, these aerosols affect atmospheric stability and cloud cover (the semi-direct effect, or SDE). Furthermore, in regions with persistent high loadings of absorbing aerosols (e.g. Asia), regional circulation patterns may be altered, potentially resulting in changes in precipitation patterns. Here we investigate aerosol-climate coupling using the NASA Goddard Earth Observing System model version 5 (GEOS-5) atmospheric general circulation model (AGCM), in which we have implemented an online version of the Goddard Chemistry, Aerosol, Radiation and Transport (GOCART) model. GOCART includes representations of the sources, sinks, and chemical transformation of externally mixed dust, sea salt, sulfate, and carbonaceous aerosols. We examine a series of free-running ensemble climate simulations of the present-day period (2000-2009) forced by observed sea surface temperatures to determine the impact of aerosols on the model climate. The SDE and response of each simulation is determined by differencing with respect to the control simulation (no aerosol forcing). In a free-running model, any estimate of the SDE includes changes in clouds due both to atmospheric heating from aerosols and changes in circulation. To try and quantify the SDE without these circulation changes we then examine the DARF and SDE in GEOS-5 with prescribed meteorological analyses introduced by the MERRA analysis. By doing so, we are able to examine changes in model clouds that occur on shorter scales (six hours). In the GEOS-5 data assimilation system (DAS), the analysis is defined as the best estimate of the atmospheric state at any given time, and it is determined by optimally combining a first-guess short-term GCM forecast with all available observations. The Incremental Analysis Update (IAU) is added to the model forecast tendencies to align them with the analysis every six hours, thus preventing longer timescale feedbacks due to the aerosol forcing. We calculate the SDE by comparing model runs with and without aerosols, and the difference in the IAU between these runs is a useful metric with which to evaluate the impact of the SDE on the model atmosphere and clouds. Decreasing the IAU indicates that the aerosol direct and semi-direct effects act to reduce the bias between the model and observations and vice versa.

Randles, Cynthia A.↗

A BI-LEVEL APPROACH TO FUTURE POWER SYSTEM CO-OPTIMIZATION WITH HIGH PENETRATION OF RENEWABLE ENERGY AND RESPONSIVE DEMAND

Electricity generation accounts for nearly half of the total CO2 emissions in the United States. For this reason, the development and integration of renewable resources will play an essential role in achieving the societal objective of mitigating climate change through reduction of greenhouse gas emissions. In conjunction with the environmental benefits of renewable energy, the most common renewable sources, such as wind and solar, also increase the uncertainties surrounding generation in power systems, which adds significant challenges to the system operation and planning. The uncertainties and forecasting errors surrounding renewable generation are normally addressed through the use of reserves from traditional generators. At greater levels of renewable penetration, sufficient generator reserves may not be available or economically viable. In contrast, the promise of demand-side resources in this arena lies in the spatially widespread availability, rapid response potential, and lower cost of features that already exist in the system. However, the challenge of responsive demand arises from the need to understand, manage, and incentivize a very large number of resources to participate effectively in efficient operation of the complex power system. Since demand response comes from the distribution system, microgrids as the basic building blocks of future distribution systems will be a critical environment for the study of demand response. To support integration of microgrids with flexible loads in future power systems, the operational mode of power systems will need to evolve. Therefore, it is going to be critical to have new and efficient co-optimization methods for coordination of the various power market participants and the scheduling of resources in the power systems of the future. Motivated by the rapid increase in renewable penetration, the need for effective demand response programs, and a changing system structure, this dissertation seeks to define a new strategy that supports co-optimization of various participants in power systems with emphasis on high renewable penetration and demand response. This strategy has three components; 1) an exploration of the capabilities of different types of demand response programs in a microgrid, 2) development and implementation of a bi-level framework for co-optimization of the main grid with high renewable penetration and a microgrid with demand response capabilities, and 3) expansion of the bi-level framework from a microgrid to a general distribution system to explore the advantages of the bi-level co-optimization approach over the traditional optimization approach. Conclusions of this work illustrate that the stochastic rolling horizon approach could effectively manage the operation of a microgrid with various demand response programs. In addition, the bi-level approach is a promising co-optimization framework for the transmission and distribution levels that could increase system renewable penetration and reduce operation costs. Compared to the traditional framework, the bi-level framework yields more equitable cost sharing patterns among the market participants as well as better support for the power system evolution.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Efficient Parallelization of Irregular Applications on GPU Architectures

With the enlarging computation capacity of general Graphics Processing Units (GPUs), leveraging GPUs to accelerate parallel applications has become a critical topic in academia and industry. However, a wide range of irregular applications with the computation-/memory-intensive nature cannot easily achieve high GPU utilization. The challenges mainly involve the following aspects: first, data dependence leads to coarse-grained kernel and inefficient parallelism; second, heavy GPU memory usage may cause frequent memory evictions and extra overhead of I/O; third, specific computation patterns produce memory redundancies; last, workload balance and data reusability conjunctly benefit the overall performance, but there may exist a dynamic trade-off between them. Targeting these challenges, this dissertation proposes multiple optimizations to accelerate two real-world applications: many-body correlation functions to simulate nuclear physics in a large-scale scientific system; the other is the eALS-based matrix factorization recommendation system. To accelerate the calculations of many-body correlation functions, this dissertation presents three frameworks in GPU memory management and multi-GPU scheduling. Firstly, an optimized systematic GPU memory management framework, MemHC, utilizes a series of new memory reduction designs in GPU memory allocation, CPU/GPU communications, and GPU memory oversubscription. Secondly, an enhanced multi-GPU scheduling framework, MICCO, particularly by taking both data dimension (e.g., data reuse and data eviction) and computation dimension into account. MICCO designs a heuristic scheduling algorithm and a machine learning-based regression model to generate the optimal settings of a proposed new concept to manage the trade-off. Thirdly, a locality-aware multi-GPU scheduling framework. This scheduler leverages pipeline batch generation with a looking-ahead strategy by building local dependency graphs for memory transfer reduction and better data reuse, achieving up to 79.92% memory cost reduction and 1.67x speedup. To parallelize the eALS-based recommendation system, this dissertation proposes an efficient CPU/GPU heterogeneous recommendation system, HEALS. HEALS employs newly designed architecture-adaptive data formats to achieve load balance and good data locality on CPU and GPU. To mitigate the data dependence, HEALS presents a CPU/GPU collaboration model for both task parallelism and data parallelism with multiple kernel computation optimizations. In summary, this dissertation efficiently accelerates two typical irregular applications on GPUs by building four frameworks, including CPU/GPU collaboration, GPU memory management, and multi-GPU scheduling.

Wang, Qihan↗