Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy usage”

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 955 records · Page 53

Vibration Testing of an Operating Stirling Convertor

The NASA John H. Glenn Research Center and the U.S. Department of Energy are currently developing a Stirling convertor for use as an advanced spacecraft power system for future NASA deep-space missions. As part of this development, a Stirling Technology Demonstrator Convertor (TDC) was recently tested to verify its survivability and capability of withstanding its expected launch random vibration environment. The TDC was fully operational (producing power) during the random vibration testing. The output power of the convertor was measured during the testing, and these results are discussed in this paper. Numerous accelerometers and force gauges were also present which provided information on the dynamic characteristics of the TDC and an indication of any possible damage due to vibration. These measurements will also be discussed in this paper. The vibration testing of the Stirling TDC was extremely successful. The TDC survived all its vibration testing with no structural damage or functional performance degradation. As a result of this testing, the Stirling convertor's capability to withstand vibration has been demonstrated, enabling its usage in future spacecraft power systems.

Hughes, William O.↗

US Department of Energy, Office of Science High Performance Computing Facility Operational Assessment 2021: Oak Ridge Leadership Computing Facility

Oak Ridge National Laboratory’s (ORNL’s) Leadership Computing Facility (OLCF) continues to surpass its operational target goals of supporting users; delivering fast, reliable computational ecosystems; creating innovative solutions for high-performance computing (HPC) needs; contributing to the community to build the next generation HPC workforce, and managing risks, safety, and security associated with operating some of the most powerful computers in the world. The results can be seen in the cutting-edge science conducted by users and the praise from the research community. Calendar year (CY) 2021 saw continued excellence in research supported by the OLCF’s leadership-class computing resources, including Summit (the nation’s most powerful supercomputer), the global scratch file system Alpine, the Scalable Protected Infrastructure (SPI), the Exploratory Visualization Environment for Research in Science and Technology (EVEREST), and the archival mass-storage resource High-Performance Storage System (HPSS). While maintaining access and exceptional user support for Summit, the OLCF continued to make progress on the installation and deployment of Frontier, which will be the nation’s first exascale system when it comes online at the start of CY 2023. Users have already begun running and optimizing scientific codes on Crusher, the OLCF test and development system equipped with Frontier’s architecture. Throughout the year, the OLCF maintained a strong culture of operational excellence, including risk management, workplace safety, and cybersecurity. The OLCF’s rigorous risk management strategy anticipated and mitigated risks, and at this time there are no high-priority operational risks. Similarly, ORNL and the OLCF were committed to operating under the US Department of Energy’s (DOE’s) safety regulations that ensure a safe workplace. Technical staff tracked and monitored existing threats and vulnerabilities within the OLCF while continually developing tools and practices to enhance operations without increasing the facility’s risk. CY 2021 was filled with outstanding results and accomplishments, including a very high rating from users on overall satisfaction for the eighth consecutive year; a tremendous number of node hours delivered to 1,671 researchers on Summit; and the successful delivery of the allocation split of roughly 60%, 20%, and 20% of core-hours offered for the Innovative and Novel Computational Impact on Theory and Experiment (INCITE), Advanced Scientific Computing Research Leadership Computing Challenge (ALCC), and Director’s Discretionary (DD) programs, respectively (Section 2). COVID-19 research remained a focus in 2021, and the ALCC and DD programs allocated over 1 million Summit hours to the COVID-19 High Performance Computing Consortium. These accomplishments, coupled with the high utilization rates (i.e., overall and capability usage), represent the fulfillment of the promise of leadership class machines: efficient facilitation of leadership-class computational applications.

97 MATHEMATICS AND COMPUTING↗

Fluid-Kinetic Coupling: Advanced Discretizations for Simulations on Emerging Heterogeneous Architectures (LDRD FY20-0643)

Plasma physics simulations are vital for a host of Sandia mission concerns, for fundamental science, and for clean energy in the form of fusion power. Sandia's most mature plasma physics simulation capabilities come in the form of particle-in-cell (PIC) models and magnetohydrodynamics (MHD) models. MHD models for a plasma work well in denser plasma regimes when there is enough material that the plasma approximates a fluid. PIC models, on the other hand, work well in lower-density regimes, in which there is not too much to simulate; error in PIC scales as the square root of the number of particles, making high-accuracy simulations expensive. Real-world applications, however, almost always involve a transition region between the high-density regimes where MHD is appropriate, and the low-density regimes for PIC. In such a transition region, a direct discretization of Vlasov is appropriate. Such discretizations come with their own computational costs, however; the phase-space mesh for Vlasov can involve up to six dimensions (seven if time is included), and to apply appropriate homogeneous boundary conditions in velocity space requires meshing a substantial padding region to ensure that the distribution remains sufficiently close to zero at the velocity boundaries. Moreover, for collisional plasmas, the right-hand side of the Vlasov equation is a collision operator, which is non-local in velocity space, and which may dominate the cost of the Vlasov solver. The present LDRD project endeavors to develop modern, foundational tools for the development of continuum-kinetic Vlasov solvers, using the discontinuous Petrov-Galerkin (DPG) methodology, for discretization of Vlasov, and machine-learning (ML) models to enable efficient evaluation of collision operators. DPG affords several key advantages. First, it has a built-in, robust error indicator, allowing us to adapt the mesh in a very natural way, enabling a coarse velocity-space mesh near the homogeneous boundaries, and a fine mesh where the solution has fine features. Second, it is an inherently high-order, high-intensity method, requiring extra local computations to determine so-called optimal test functions, which makes it particularly suited to modern hardware in which floating-point throughput is increasing at a faster rate than memory bandwidth. Finally, DPG is a residual-minimizing method, which enables high-accuracy computation: in typical cases, the method delivers something very close to the $L^2$ projection of the exact solution. Meanwhile, the ML-based collision model we adopt affords a cost structure that scales as the square root of a standard direct evaluation. Moreover, we design our model to conserve mass, momentum, and energy by construction, and our approach to training is highly flexible, in that it can incorporate not only synthetic data from direct-simulation Monte Carlo (DSMC) codes, but also experimental data. We have developed two DPG formulations for Vlasov-Poisson: a time-marching, backward-Euler discretization and a space-time discretization. We have conducted a number of numerical experiments to verify the approach in a 1D1V setting. In this report, we detail these formulations and experiments. We also summarize some new theoretical results developed as part of this project (published as papers previously): some new analysis of DPG for the convection-reaction problem (of which the Vlasov equation is an instance), a new exponential integrator for DPG, and some numerical exploration of various DPG-based time-marching approaches to the heat equation. As part of this work, we have contributed extensively to the Camellia open-source library; we also describe the new capabilities and their usage. We have also developed a well-documented methodology for single-species collision operators, which we applied to argon and demonstrated with numerical experiments. We summarize those results here, as well as describing at a high level a design extending the methodology to multi-species operators. We have released a new open-source library, MLC, under a BSD license; we include a summary of its capabilities as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Replacing photomultiplier tubes with silicon photomultipliers for nuclear safeguards applications

Photomultiplier tubes (PMTs) have been used for decades as the default light detection technology for scintillator-based radiation monitors. PMTs come with a handful of disadvantages, however, including large volume usage, fragility, high-voltage requirements, and susceptibility to magnetic fields. Arrays of silicon photomultipliers (SiPMs) are a possible alternative to PMTs providing similar performance while offering improvements in the areas listed. We focus on straightforward “drop-in replacement” evaluations by imposing a single channel output of the SiPM signal and a backend electronics data acquisition and analysis routine identical to that used for PMTs. In the realm of nuclear safeguards, the primary performance metric for gamma spectroscopy is detector resolution, and for neutron detection it is pulse shape discrimination (PSD) to separate neutron and gamma signals. In this work we present the results of replacing a PMT with a 2” x 2” SiPM array mounted, in turn, to the same 2” x 2” cylindrical sodium iodide (NaI) crystal. The ultimate comparison for this gamma spectrometry is confirmation of uranium enrichment standards, ranging from depleted uranium to 93% highly enriched uranium (HEU), by analyzing the resulting spectra with both a simple integral scaling in key energy regions of the spectrum, as well as with the NaIGEM software package. In addition to the gamma spectroscopy work, we will present the latest results on comparing the PSD capabilities of a 4” x 4” SiPM array to a 5” PMT. As with the gamma spectroscopy evaluations, the light detectors are mounted in turn to the same liter-scale organic scintillators. This neutron-focused work includes evaluation of prototype SiPM readout boards with a single output signal summed over 256 individual SiPM pixels, and evident tradeoffs between PSD capability and fast response.

42 ENGINEERING↗

Replacing photomultiplier tubes with silicon photomultipliers for nuclear safeguards applications

Photomultiplier tubes (PMTs) have been used for decades as the default light detection technology for scintillator-based radiation monitors. PMTs come with a handful of disadvantages, however, including large volume usage, fragility, high-voltage requirements, and susceptibility to magnetic fields. Arrays of silicon photomultipliers (SiPMs) are a possible alternative to PMTs providing similar performance while offering improvements in the areas listed. We focus on straightforward “drop-in replacement” evaluations by imposing a single channel output of the SiPM signal and a backend electronics data acquisition and analysis routine identical to that used for PMTs. In the realm of nuclear safeguards, the primary performance metric for gamma spectroscopy is detector resolution, and for neutron detection it is pulse shape discrimination (PSD) to separate neutron and gamma signals. In this work we present the results of replacing a PMT with a 2” x 2” SiPM array mounted, in turn, to the same 2” x 2” cylindrical sodium iodide (NaI) crystal. The ultimate comparison for this gamma spectrometry is confirmation of uranium enrichment standards, ranging from depleted uranium to 93% highly enriched uranium (HEU), by analyzing the resulting spectra with both a simple integral scaling in key energy regions of the spectrum, as well as with the NaIGEM software package. In addition to the gamma spectroscopy work, we will present the latest results on comparing the PSD capabilities of a 4” x 4” SiPM array to a 5” PMT. As with the gamma spectroscopy evaluations, the light detectors are mounted in turn to the same liter-scale organic scintillators. This neutron-focused work includes evaluation of prototype SiPM readout boards with a single output signal summed over 256 individual SiPM pixels, and evident tradeoffs between PSD capability and fast response.

Engineering - Instrumentation related to nuclear s↗

Particle hit clustering and identification using point set transformers in liquid argon time projection chambers

Liquid argon time projection chambers are often used in neutrino physics and dark-matter searches because of their high spatial resolution. The images generated by these detectors are extremely sparse, as the energy values detected by most of the detector are equal to 0, meaning that despite their high resolution, most of the detector is unused in a particular interaction. Instead of representing all of the empty detections, the interaction is usually stored as a sparse matrix, a list of detection locations paired with their energy values. Traditional machine learning methods that have been applied to particle reconstruction such as convolutional neural networks (CNNs), however, cannot operate over data stored in this way and therefore must have the matrix fully instantiated as a dense matrix. Operating on dense matrices requires a lot of memory and computation time, in contrast to directly operating on the sparse matrix. We propose a machine learning model using a point set neural network that operates over a sparse matrix, greatly improving both processing speed and accuracy over methods that instantiate the dense matrix, as well as over other methods that operate over sparse matrices. Compared to competing state-of-the-art methods, our method improves classification performance by 14%, segmentation performance by more than 22%, while taking 80% less time and using 66% less memory. Compared to state-of-the-art CNN methods, our method improves classification performance by more than 86%, segmentation performance by more than 71%, while reducing runtime by 91% and reducing memory usage by 61%.

calibration and fitting methods↗

BETO 2021 Peer Review - Analytical Development and Standardization for Biomass-Derived Thermochemical Liquids

This project began in FY14 to address the lack of standard chemical characterization analytical methods for bio-oils. Bio-oils are very complex and present numerous analytical challenges; yet reliable chemical information (quantification of both individual compounds and chemical functional groups) is needed to inform upgrading research and refinery co-processing. In this project, analysis needs are first determined from engaging the bioenergy community. Next, standard methods are developed to meet these needs, and then subsequently validated via inter-laboratory studies. Methods that are successfully validated (< 10% variability) are then shared as Laboratory Analytical Procedures (LAPs), which are free and publicly available. We have been tracking LAP use and have seen sustained usage as evidenced by an average of 500 pages views and 100 downloads per quarter, demonstrating the value of these methods to the bioenergy community. LAP methods that are particularly useful and reliable will be chosen for the next-level of standardization through ASTM. We have recently achieved approval by ASTM for our carbonyl titration method. This method (ASTM E3146) is the first example of an ASTM standard solely focused on the chemical characterization of bio-oils. Work in this project is meeting the analysis needs of the bioenergy community and will ultimately help enable the commoditization of bio-oils.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

GeoThermalCloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗

Transit Rider/Travel Behavior Inventory Survey - Minneapolis-St. Paul Metro - 2010

The Metropolitan Council administered a system-wide transit on-board survey between September and November 2010 as part of the 2010 Travel Behavior Inventory to provide detailed transit usage patterns and rider information to support modeling and planning efforts. The study’s focus was to capture only the most current and reliable data necessary to determine future public transportation needs in the Minneapolis region. The study was structured to collect detailed transit ridership data for different routes during different times of day to develop a disaggregate transit trip table that will support advanced travel demand modeling. The study was designed to leverage existing data sources, such as the 2005 on-board survey, to provide the best quality data to update the Travel Behavior Inventory in the Minneapolis-Saint Paul metropolitan area.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

2022 LANL Radionuclide Air Emissions Report (Rev. 2)

This report describes the emissions of airborne radionuclides from operations at Los Alamos National Laboratory (LANL) for calendar year 2022 and the resulting off-site dose from these emissions. This document fulfills the requirements established by the National Emissions Standards for Hazardous Air Pollutants in 40 CFR 61, Subpart H – Emissions of Radionuclides other than Radon from Department of Energy Facilities, commonly referred to as the Radionuclide NESHAP or Rad-NESHAP. Compliance with this regulation and preparation of this document is the responsibility of LANL’s Rad NESHAP compliance program, which is part of the Environmental Protection and Compliance (EPC) Division. The information in this report is required under the Clean Air Act and is being submitted to the U.S. Environmental Protection Agency (EPA) Headquarters and EPA Region 6. The highest effective dose equivalent (EDE) to an off-site member of the public was calculated using procedures specified by the EPA and described in this report. LANL’s EDE was 0.45 for 2022. The annual limit is 10 millirem per year, established by the EPA in 40 CFR 61 Subpart H. All measured air emissions are modeled to a single location, known as the Maximally Exposed Individual (MEI). During calendar year 2022, LANL continuously monitored radionuclide emissions at 27 “major” release points, or stacks. The Laboratory estimates emissions from an additional 34 “minor” release points using radionuclide usage source terms in lieu of stack monitoring. Also, LANL uses an EPA approved network of air samplers around the Laboratory perimeter to monitor ambient airborne levels of radionuclides. To provide data for dispersion modeling and dose assessment, LANL maintains and operates several meteorological monitoring towers. From these various systems, a comprehensive evaluation is conducted to calculate the MEI dose for the Laboratory. The MEI can be any member of the public at any off-site location where there is a residence, school, business, or office. In 2022, this MEI location was a business at 95 Entrada Drive, located in the eastern end of Los Alamos town site. The primary contributors to the off-site dose at this location are the ambient air data at that location combined with radioactive gas emissions from the LANSCE facility and the collected potential emissions from unmonitored (minor) sources. Overall, the MEI dose in 2022 is similar to that which has been observed in recent years, and it remains well below the EPA’s 10 millirem per year limit. Doses reported to the EPA for the past 10 years are shown in Table E1.

54 ENVIRONMENTAL SCIENCES↗

Notice of Submittal – 2023 Radionuclide Air Emissions Report for Los Alamos National Laboratory

This report describes the emissions of airborne radionuclides from operations at Los Alamos National Laboratory (LANL) for calendar year 2023 and the resulting off-site dose from these emissions. This document fulfills the requirements established by the National Emissions Standards for Hazardous Air Pollutants in 40 CFR 61, Subpart H – Emissions of Radionuclides other than Radon from Department of Energy Facilities, commonly referred to as the Radionuclide NESHAP or Rad-NESHAP. Compliance with this regulation and preparation of this document is the responsibility of LANL’s Rad NESHAP compliance program, which is part of the Environmental Protection and Compliance (EPC) Division. The information in this report is required under the Clean Air Act and is being submitted to the U.S. Environmental Protection Agency (EPA) Headquarters and EPA Region 6. The highest effective dose equivalent (EDE) to an off-site member of the public was calculated using procedures specified by the EPA and described in this report. LANL’s EDE was 0.43 for 2023. The annual limit is 10 millirem per year, established by the EPA in 40 CFR 61 Subpart H. All measured air emissions are modeled to a single location, known as the Maximally Exposed Individual (MEI). During calendar year 2023, LANL continuously monitored radionuclide emissions at 28 “major” release points, or stacks. The Laboratory estimates emissions from an additional 59 “minor” release points using radionuclide usage source terms in lieu of stack monitoring. Also, LANL uses an EPA approved network of air samplers around the Laboratory perimeter to monitor ambient airborne levels of radionuclides. To provide data for dispersion modeling and dose assessment, LANL maintains and operates several meteorological monitoring towers. From these various systems, a comprehensive evaluation is conducted to calculate the MEI dose for the Laboratory. The MEI can be any member of the public at any off-site location where there is a residence, school, business, or office. In 2023, this MEI location was a business at 129 New Mexico State Road 4 (NM-4), located in the northern end of White Rock. The primary contributors to the off-site dose at this location are the ambient air data at that location combined with the collected potential emissions from unmonitored (minor) sources. Overall, the MEI dose in 2023 is similar to that which has been observed in recent years, and it remains well below the EPA’s 10 millirem per year limit. Doses reported to the EPA for the past 10 years are shown in Table E1.

54 ENVIRONMENTAL SCIENCES↗

Real-Sim Interface: Enabling Multi-resolution Simulation and X-in-the-Loop Development for Connected and Automated Vehicles

We report connected and automated vehicles (CAVs) can bring safety, mobility, and energy benefits to transportation systems. Ideally, CAV applications would be fully evaluated and validated prior to real-world implementation. However, many technical challenges in both software and hardware hinder the process. To comprehensively evaluate all aspects of CAV applications, an integrated evaluation environment is needed with various simulation tools from different domains. In the current literature, there lacks a well-developed interface to enable multi-resolution simulation of vehicle, traffic, virtual environment, and hardware-in-the-loop (HIL) simulation. In this work, a modular and flexible interface is developed to enable multi-resolution vehicle and traffic co-simulation for CAV applications. This interface is built upon Oak Ridge National Laboratory’s (ORNL’s) Real-Sim approach, can support various simulation tools and real-time X-in-the-loop (XIL) simulations, and is based on network communication protocols. The network communication delays are analyzed for simultaneous connection of up to 20 simulators. Then, two example applications are studied to demonstrate the potential usage of the Real-Sim interface: (1) a centralized merging scenario where the dynamics of two ego vehicles are emulated with detailed vehicle and powertrain dynamics models and (2) a signal-controller-in-the-loop (SCIL) evaluation where one intersection of the simulated traffic scenario is controlled by a physical signal controller. This Real-Sim interface is a tool that allows researchers to bring real, tangible hardware and software into simulated environments to comprehensively evaluate and support a wide variety of CAV applications.

33 ADVANCED PROPULSION SYSTEMS↗

Vertiport Dynamic Density

Advanced Air Mobility (AAM) is envisioned to be another spoke in a region’s transportation system, supplementing ground-based travel with air-based travel. Eventually, air taxis will be pervasive, convenient, and affordable. As demand and operations tempo increase, congestion may arise. Several characteristics of AAM reduce the applicability of techniques used in today’s National Airspace System (NAS) to manage congestion. AAM will include both scheduled and non-scheduled, on-demand operations, challenging strategic deconfliction algorithms. Flights will be fairly short, traversing an urban area, not hundreds or thousands of miles, with corresponding low energy reserves, prohibiting excessive delays. Flight operators will need flexibility in operations, scheduling a flight only minutes before departure, or diverting to an alternate vertiport if it becomes advantageous, further adding to trajectory uncertainty that would challenge strategic deconfliction. Finally, operators will desire privacy to protect sensitive information or preserve competitive advantage, limiting early access to intent information. In this paper, we present an approach for managing congestion at vertiports by providing insight into the traffic situation to support operationally-advantageous and safe land or divert decisions. We propose a metric designed with usefulness and usability in AAM operations in mind. The metric uses the sociology concept of dynamic density (DD) that takes into consideration not only number of flights, but also the interaction of those flights with the vertiport’s limited resources, namely the landing pads and the parking spots. DD supports a Pilot in Command (PIC) with decisions about whether to proceed, expedite, delay, or divert; and supports air traffic control (ATC) and vertiport operators in airspace management and vertiport usage. We demonstrate the metric on a notional vertiport scenario. We also show that DD provides better insight into congestion and resulting flight delays than an aircraft count metric used in traditional air traffic management.

Lilly Spirkovska↗

Application of Monte Carlo techniques to optimization of high-energy beam transport in a stochastic environment

An algorithm employing a modified sequential random perturbation, or creeping random search, was applied to the problem of optimizing the parameters of a high-energy beam transport system. The stochastic solution of the mathematical model for first-order magnetic-field expansion allows the inclusion of state-variable constraints, and the inclusion of parameter constraints allowed by the method of algorithm application eliminates the possibility of infeasible solutions. The mathematical model and the algorithm were programmed for a real-time simulation facility; thus, two important features are provided to the beam designer: (1) a strong degree of man-machine communication (even to the extent of bypassing the algorithm and applying analog-matching techniques), and (2) extensive graphics for displaying information concerning both algorithm operation and transport-system behavior. Chromatic aberration was also included in the mathematical model and in the optimization process. Results presented show this method as yielding better solutions (in terms of resolutions) to the particular problem than those of a standard analog program as well as demonstrating flexibility, in terms of elements, constraints, and chromatic aberration, allowed by user interaction with both the algorithm and the stochastic model. Example of slit usage and a limited comparison of predicted results and actual results obtained with a 600 MeV cyclotron are given.

Parrish, R. V.↗

Tectonic-Climate Interactions in Action Orogenic Belts: Quantification of Dynamic Topography with SRTM data

This project was undertaken to examine the approach to steady state in collisional mountain belts. Although the primary thrust of this grant was to look at larger collisional mountain belts, such as the Himalaya, the Tien Shan, and Southern Alps, we began by looking at smaller structures represented by growing and propagating folds. Like ranges that are evolving toward a topographic steady state, these folds undergo a series of morphologic changes as they are progressively uplifted and eroded. We wanted to document the nature of these changes and to try to discern some of the underlying controls on them. We initially focused on the Wheeler Ridge anticline in southern California. Subsequently, we progressed to looking at the topographic development and the effects of differential uplift and glaciation on the Kyrgyz Range in the northern Tien Shan. This range is unusual inasmuch as it is transformed along its length from a simple uplift with a largely preserved Mesozoic erosion surface arching across it to a highly dissected and heavily glaciated uplift in the region where uplift has been sustained at higher rates over longer intervals. In efforts to understand the distribution of erosion rates at 10(exp 3) - 10(exp 5) year time scales, cosmogenic radionuclide (CRN) concentrations have been gaining increasingly widespread usage (Brown et al., 1995; Riebe et al., 2004; Riebe et al., 2001; Vance et al., 2003). Most studies to date, however, have been conducted in slowly eroding ranges. In rapidly eroding mountains where landslides deliver most of the sediments to the rivers, we hypothesized that CRN concentrations could be highly perturbed by the stochastic processes of landsliding. Therefore, we undertook the development of a numerical model that simulated the effects of both landsliding and grain-by-grain attrition within fluvial catchments. This modeling effort has shown the effects of catchment size and erosion rate on CRN concentrations and allows a prediction of where to sample to obtain the optimal erosion rate estimates using CRN techniques. Finally, we developed computational techniques to operate on DEMs to extract useful information that would enable quantification of climate-erosion interactions. In particular, we worked on rapid techniques to define catchments of any given range of sizes, to extract channel gradients, to combine precipitation information to calculate discharge, and to utilize various stream-power models to determine the erosional energy within any given catchment within a transect. We briefly describe results from Wheeler Ridge, the Kyrgyz Range, the Nepal Himalaya, and our numerical modeling.

Burbank, Douglas W.↗

Adoption of delivery services in light of the COVID pandemic: Who and how long?

A significant growth in demand for online shopping in light of the Coronavirus Disease (COVID) crisis has received attention from transportation practitioners, policy-makers, and researchers. However, an important question arises in this increase in online shopping and resulting deliveries: How long will this last? Very little is known whether this popularity would last a long time. Here, to address this question, the authors conducted a survey of 915 individuals residing in the U.S. and classified them into the four distinctive consumer types (i.e., the prior adopter, temporary adopter and permanent new adopter, and non-adopter) depending on their usage of delivery services before, during, and after (expected) the COVID crisis. This research aims to gain behavioral insight by exploring the differences between the four consumer types and investigating factors affecting the initial adoption and continuance intention of using delivery services. The descriptive analysis revealed that there are clear differences not only between the four types of consumers but also between the four product types (i.e., grocery, food, home goods, and other packages) considered in the survey. The models found that factors affecting the initial adoption and continuance intention are different from the previous studies conducted before the COVID pandemic. Implications for planning and policymaking are also discussed.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Mission and Navigation Design for the 2009 Mars Science Laboratory Mission

NASA s Mars Science Laboratory mission will launch the next mobile science laboratory to Mars in the fall of 2009 with arrival at Mars occurring in the summer of 2010. A heat shield, parachute, and rocket-powered descent stage, including a sky crane, will be used to land the rover safely on the surface of Mars. The direction of the atmospheric entry vehicle lift vector will be controlled by a hypersonic entry guidance algorithm to compensate for entry trajectory errors and counteract atmospheric and aerodynamic dispersions. The key challenges for mission design are (1) develop a launch/arrival strategy that provides communications coverage during the Entry, Descent, and Landing phase either from an X-band direct-to-Earth link or from a Ultra High Frequency link to the Mars Reconnaissance Orbiter for landing latitudes between 30 deg North and 30 deg South, while satisfying mission constraints on Earth departure energy and Mars atmospheric entry speed, and (2) generate Earth-departure targets for the Atlas V-541 launch vehicle for the specified launch/arrival strategy. The launch/arrival strategy employs a 30-day baseline launch period and a 27-day extended launch period with varying arrival dates at Mars. The key challenges for navigation design are (1) deliver the spacecraft to the atmospheric entry interface point (Mars radius of 3522.2 km) with an inertial entry flight path angle error of +/- 0.20 deg (3 sigma), (2) provide knowledge of the entry state vector accurate to +/- 2.8 km (3 sigma) in position and +/- 2.0 m/s (3 sigma) in velocity for initializing the entry guidance algorithm, and (3) ensure a 99% probability of successful delivery at Mars with respect to available cruise stage propellant. Orbit determination is accomplished via ground processing of multiple complimentary radiometric data types: Doppler, range, and Delta-Differential One-way Ranging (a Very Long Baseline Interferometry measurement). The navigation strategy makes use of up to five interplanetary trajectory correction maneuvers to achieve entry targeting requirements. The requirements for cruise propellant usage and atmospheric entry targeting and knowledge are met with ample margins.

Mission design↗

Quantum-Accurate Multiscale Modeling of Shock Hugoniots, Ramp Compression Paths, Structural and Magnetic Phase Transitions, and Transport Properties in Highly Compressed Metals

Fully characterizing high energy density (HED) phenomena using pulsed power facilities (Z machine) and coherent light sources is possible only with complementary numerical modeling for design, diagnostic development, and data interpretation. The exercise of creating numerical tests, that match experimental conditions, builds critical insight that is crucial for the development of a strong fundamental understanding of the physics behind HED phenomena and for the design of next generation pulsed power facilities. The persistence of electron correlation in HED materials arising from Coulomb interactions and the Pauli exclusion principle is one of the greatest challenges for accurate numerical modeling and has hitherto impeded our ability to model HED phenomena across multiple length and time scales at sufficient accuracy. An exemplar is a ferromagnetic material like iron, while familiar and widely used, we lack a simulation capability to characterize the interplay of structure and magnetic effects that govern material strength, kinetics of phase transitions and other transport properties. Herein we construct and demonstrate the Molecular-Spin Dynamics (MSD) simulation capability for iron from ambient to earth core conditions, all software advances are open source and presently available for broad usage. These methods are multi-scale in nature, direct comparisons between high fidelity density functional theory (DFT) and linear-scaling MSD simulations is done throughout this work, with advancements made to MSD allowing for electronic structure changes being reflected in classical dynamics. Main takeaways for the project include insight into the role of magnetic spins on mechanical properties and thermal conductivity, development of accurate interatomic potentials paired with spin Hamiltonians, and characterization of the high pressure melt boundary that is of critical importance to planetary modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗