Search NASA⌕ Search

SEARCH · Search NASA

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 163 records · Page 9

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.↗

Measuring, Cutting and Sorting Irradiated Control Rods to Limit the Volume of Intermediate Level Long-lived Waste - 17001

During operation, the PWR control rods are either completely extracted (scram or safety rods) or partially inserted in the operating zone of the core (fine control or operating rods). At the nominal rating conditions or close to them, the operating zone is limited to some steps of insertion in the order of a few tens of cm. In fact, only the lower part of the control rods are subjected to a significant neutron flux with an impact on the neutron activation and so on the waste classification (Low Level Short-lived Waste versus Intermediate Level Long-lived Waste). Our calculative as well as experimental experience confirms the significant presence of radionuclides in the control rods, but only in the lower parts. Measuring, cutting and finally sorting control rod scraps can reduce the volume of Intermediate Level Long-lived Waste planned for our storage deposit facility. This paper presents the EDF/DP2D methodologies followed to determine the sort criteria between the Low Level Short-lived Waste versus the Intermediate Level Long-lived Waste zones for the Silver-Indium-Cadmium control rods. The irradiated control rods are stored in cases without their head. Each case can contain several control rods, and the cases are stored under water in the fuel assembly pool. To calculate the radioactive inventories, we have developed a calculation scheme which is composed of a 3-dimensional mapping of the neutron flux, a calculation of the evolution of the isotopic compositions under neutron flux and finally a comparison between calculated and measured values. The main hypotheses of the control rods calculation take into account: - 3-dimensional geometry and 3-dimensional neutron sources distribution, - the isotopic compositions including impurities, - the history of irradiation resulting mainly from the daily power production, the positions of the control rods in the loading patterns and the positions in the axial operating zone. In order to validate the calculations, measurements were made 20 meters under water in the fuel assembly pool. The experimental results show the presence of many gamma emitter radionuclides. Only Ag-108m was retained for validating the calculations. This is due to the fact that silver is considered as the 'standard' chemical element for our control rods and silver produces Ag-108m which is a long-lived gamma emitter radionuclide. The analysis of the calculated results and the measured values allow us to build a simplified physics model for the simulation of the control rods activation. Finally three sort criteria can be proposed: - Without measurement devices, the sorting is based on an empirical sort criterion. Cuttings must then be made on the 2 back-ends of the control rod pins. This option is due to the difficulty in knowing the orientation of the pins in the storage cases. The cutting lines are placed at approximately 1 m from the 2 back-ends of the pins. This approach does not allow a great optimization of the sorting, less than 50% of the length of the control rods can be saved from Intermediate Level Long-lived Waste storage. - With a simple dose rate measurement device, the knowledge of the orientation of the pins in the storage cases allows us to limit cuttings to the lower part of the control rods. Depending on the 'top to tail' cases, between approximately 50% and 75% of the length of the control rods can be saved from Intermediate Level Long-lived Waste storage. - With an accurate measurement device, the sorting can be based on a specific activity of Ag-108m. Over 4400 Bq/g of Ag-108m (corresponding to the Intermediate Level Long-lived Waste limit), the control rod pieces can be considered as a long-lived waste. The measurements must be made over the 2 back-ends of control rod pins to examine the cases of 'top to tail' storage. This approach allows a great optimization of the sorting. Depending on the history of irradiation, more than 75% of the length of the control rods can be saved from Intermediate Level Long-lived Waste storage. (author)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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↗

Analyzing the Effects of Atmospheric Turbulent Fluctuations on the Wake Structure of Wind Turbines and Their Blade Vibrational Dynamics

In recent trends, a rising demand for renewable energy has driven wind turbines to larger proportions, where lighter blade designs are often adopted to reduce the costs associated with logistics and production. This causes modern utility-scale wind turbine blades to be inherently more flexible, and their amplified aeroelastic sensitivity results in complex multi-physics reactions to variant atmospheric conditions, including dynamic patterns of aerodynamic loading at the rotor and vortex structure evolutions within the wake. In this paper, we analyze the influence of inflow variance for wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex wake interactions for a range of flexible NRT blade variations, operating in differing conditions of variant atmospheric flow. CODEF solutions of turbine operation in Steady-In-The-Average (SITA) wind conditions are compared to SITA wind conditions featuring a controlled gust-like pulse overimposed, to isolate the effects of typical wind fluctuations. Finally, simulations of realistic time-varying wind conditions from SWiFT meteorological tower measurements are compared to the solutions of SITA wind conditions. These increasingly complex atmospheric inflow variations are tested to show the differing effects evoked by various patterns of spatiotemporal atmospheric flow fluctuations. An analysis is presented for solutions of wind turbine aeroelastic response and vortex wake evolution, to elucidate the consequences of variant inflow, which pertain to wind turbine dynamics at an individual and farm-collective scale. The comparisons of simulated farm flow for SITA and measured fluctuating wind conditions show that certain regions of the wake contain up to a 12% difference in normalized axial velocity, due to the introduction of wind fluctuations. The findings of this study prove valuable for practical applications in wind farm control and optimization strategies, with particular significance for modern utility-scale wind power plants operating in variant atmospheric conditions.

Farrell, Alayna (ORCID:000000023555720X)↗

Next Generation Hydrogen Storage Vessels Enabled by Carbon Fiber Infusion with a Low Viscosity, High Toughness System

Onboard hydrogen storage for fuel cell vehicles holds considerable commercial promise which could be accelerated by lower cost storage vessels. In past years, the use of lower cost reinforcing fiber had been studied as an intuitive avenue to reduce vessel costs. Although the matrix resin does not play a large role in bearing stress loads in many areas of a filament wound pressure vessel, the matrix is relied upon for stress transfer between fibers and appears to play an important role in certain regions of the vessel construction such as the “tangent region” where the cylinder section transitions to the dome. Also, the matrix resin can play a role in the damage tolerance of vessels after impact events. Furthermore, the effect of dry spots or voids in the laminate can affect stress transfer between fibers in a manner that is challenging to quantify. Therefore, excess fiber must be used to account for this variability in performance. Unfortunately, in the case hydrogen pressure vessels the cost of the vessel is strongly dependent on the amount of carbon-fiber used in the design and manufacturing. An extensive investigation on the use of non-traditional fabrication processes enabled by an ultra-low viscosity thermoset resin (<10 cP) with high fracture toughness was conducted to determine the opportunity to reduce carbon-fiber usage and the associated high cost. The hypothesis was based on enabling a lower variation in vessel performance from a reduction in void defects (often found with traditional wet-winding) and a more defect-tolerant vessel resulting from use of a matrix resin with high fracture toughness. These two features were anticipated to justify a more aggressive design approach to the vessel winding pattern. From the work in this project, the potential use of the vacuum infusion approach to produce high quality vessels was validated. However, the team recognized the need to re-optimize the design, the dry-process, and the vacuum infusion process at each scale-up stage. During the planning phase, the critical need to re-optimize at this frequency was under-estimated. Also, the team identified the lack of industry testing standards on subscale parts for providing predictive value of vessel performance. Therefore, significant effort is needed to produce statistically significant numbers of vessels. To produce such numbers of vessels the project needed in retrospect to partner closely with a series vessel manufactures to develop a proof of concept production line to make the vessels. However, without the appropriate data, a manufacturing partner was challenging to fully engage during this timeline. New technologies relying on changing raw materials and manufacturing processes carry inherent risks. This project did not reach its ultimate objectives of demonstrated a vessel with lower carbon-fiber. But some companies expressed interest in learning from the partial successes of this project for future challenges in pressure vessel applications.

08 HYDROGEN↗

Improving I/O Performance for Exascale Applications through Online Data Layout Reorganization

The applications being developed within the U.S. Exascale Computing Project (ECP) to run on imminent Exascale computers will generate scientific results with unprecedented fidelity and record turn-around time. Many of these codes are based on particle-mesh methods and use advanced algorithms, especially dynamic load-balancing and mesh-refinement, to achieve high performance on Exascale machines. Yet, as such algorithms improve parallel application efficiency, they raise new challenges for I/O logic due to their irregular and dynamic data distributions. Thus, while the enormous data rates of Exascale simulations already challenge existing file system write strategies, the need for efficient read and processing of generated data introduces additional constraints on the data layout strategies that can be used when writing data to secondary storage. We review these I/O challenges and introduce two online data layout reorganization approaches for achieving good tradeoffs between read and write performance. We demonstrate the benefits of using these two approaches for the ECP particle-in-cell simulation WarpX, which serves as a motif for a large class of important Exascale applications. Here, we show that by understanding application I/O patterns and carefully designing data layouts we can increase read performance by more than 80 percent.

97 MATHEMATICS AND COMPUTING↗

Coiled Tube Gas Heaters For Nuclear Gas-Brayton Power Conversion

This project developed an alternative design for heat exchangers for application to heating supercritical carbon dioxide (S-CO 2 ) or air for power conversion. We have identified an annular coiled tube bundle configuration–where hot sodium enters tubes from multiple vertical inlet manifold pipes, flows in a spiral pattern radially inward and downward, and then exits into an equal number of vertical outlet manifold pipes–as a potentially attractive option. The S-CO 2 gas or air flows radially outward through the tube bundle. Coiled tube gas heaters (CTGHs) are expected to have excellent thermal shock, long-term thermal creep, in-service inspection, and reparability characteristics, compared to alternative options. CTGHs have significant commonality with modern nuclear steam generators. Extensive experience exists with the design, manufacture, operation, in-service inspection and maintenance of nuclear steam generators. The U.S. Nuclear Regulatory Commission also has extensive experience with regulatory guidance documented in NUREG 0800. CTGHs leverage this experience and manufacturing capability. The most important difference between steam generators and gas-Brayton cycles such as the S-CO 2 cycle is that the heat exchangers must operate with counter flow with high effectiveness to minimize the pinch-point temperature difference between the hot liquid coolant and the heated gas. S-CO 2 -cycle gas heaters also operate at sufficiently elevated temperatures that time dependent creep is important and allowable stresses are relatively low. Designing heat exchangers to operate in this regime requires configurations that minimize stresses and stress concentrations. The cylindrical tubes and cylindrical manifold pipes used in CTGHs are particularly effective geometries. The first major goal of this research project was to develop and experimentally validate a detailed, 3-D multi-phase (gas-solid-liquid) heat transport model for CTGHs, using methods similar to earlier UCB multi-scale models for PCHEs, which will enable optimization of CTGH designs with respect to the number of manifold pipes, S-CO 2 and sodium circulating power, and other design parameters. CTGHs are relatively novel because the tubes loaded in compression as well as the manifold pipes which form the tube sheets. The second major goal of the experimental and modeling effort was to optimize the assembly and heat treatment of the tube-to-tubesheet joints, using a novel tapered joint configuration, and to develop diffusion-bonded joints with high shear strength and creep resistance, as has been done in other applications for tapered plugs. The project would also conduct creep testing of these tubesheet joints in pressurized furnaces under prototypical conditions of stress and temperature, however this was not feasible due to numerous hinderances. Detailed microstructural investigations was conducted on the diffusion bond interface, while the bond strength was evaluated by tensile testing. The mechanical testing quantified the overall strength of the bond and measured the pull out force as well as long term behavior, which complemented the microstructural analysis and gave a complete understanding of the diffusion bond.

30 DIRECT ENERGY CONVERSION↗

Final Report DE-FE0031785 Mohsen Ahmadian, Ph.D. Demonstration of Proof of Concept of a Multiphysics Approach for Real-Time Remote Monitoring of Dynamic Changes in Pressure and Salinity in Hydraulically Fractured Networks

Hydraulic fracturing has evolved into a multistep process with varying flow rates, carrier fluids (e.g., gel or slickwater), proppant loadings, and proppant grain sizes. As a result, primary recovery from a hydraulically fractured tight-oil reservoir is often a tiny fraction of the original oil in place, ranging between 5 and 10%. As stated in the FOA1990, “part of this problem is due to the inability of current well completion processes to effectively stimulate the entire reservoir volume in contact with the wellbore. Innovative technologies are needed that can help improve the effectiveness of reservoir completion methods, maximize stimulated reservoir volumes, and optimize recovery over the entire producing life span of a well”. We first need to enhance the current fracture diagnostic techniques to improve a well-completion design. However, detecting and delineating a subsurface hydraulic fracture is extremely difficult because the induced fracture network is only fractionally propped, and these propped fractures are generally very thin. Microseismic and tiltmeter monitoring techniques can provide information on the fracture extent but provide little or no information on the movement and final distribution of proppant or production fluids. On the other hand, electromagnetic (EM) imaging has shown the capability to monitor proppant distribution throughout the fracture area, especially in the presence of Electrically Active Proppants (EAPs). A previous EM survey of hydraulic fracturing at the Devine Fracture Pilot Site (DFPS) and subsequent EM code developments demonstrated this survey as a robust technique to remotely interrogate the extent of the EAP-filled hydraulic fracture during its propagation. The objectives of the project were threefold: (1) to capitalize on the material properties of an EAP to demonstrate remote monitoring of relative changes in pressure, pressure, and flow that are commonly encountered during production from a hydraulically fractured reservoir; (2) to evaluate EM imaging tools, to achieve Objective 1 in near real-time; and (3) to develop a multi-physics joint inversion approach to precisely predict flow patterns and physiochemical changes within an EAP-filled fracture network. This research project was built upon our previous work at the Devine Test Site managed by the Bureau of Economic Geology (BEG) at The University of Texas at Austin (UT-Austin). It also leveraged a significant investment from the Advanced Energy Consortium (AEC) to address the DOE's interest in subsurface flow, containment, and characterization by multiple signals. This three-year and three-month project succeeded in demonstrating the feasibility of a real-time dynamic fluid flow mapping technique at Technology Readiness Level 5 (TRL5) by utilizing a commercially available surface-based Controlled-Source Electromagnetic (CSEM) method (Objectives 1, 2). We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, the observed electric field in our study is influenced by bottomhole pressure, flow rate, and salinity, which is demonstrated by comparing these parameters with the electrical field potentials. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce the whole measured electric field magnitudes. Preliminary estimates showed that including Streaming Potential (SP) in our geophysical model is likely needed to reduce the simulation misfit.

02 PETROLEUM↗