SEARCH · Engineering Papers
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.
METTS Plot If you have set up latex form for matlibplot, you can run Plots.ipynb on your local machine. Otherwise, it will be good to run it on colab. METTS_Plot METTS Data METTS_data The plots and AVQMETTS data AVQMETTS data and plots
The ElectricityLCI v2 Python package (https://github.com/USEPA/ElectricityLCI/tree/v2.0) was used to generate the 2022 electricity baseline: a regionalized life cycle inventory model of U.S. electricity generation, consumption, and distribution using standardized facility and generation data. ElectricityLCI implements a local data store for downloading and accessing public data on an individual's computer. The data store follows the folder definition provided by USEPA's esupy Python package (https://github.com/USEPA/esupy), which utilized the appdirs Python dependency (https://pypi.org/project/appdirs/). This submission includes the background data used to generate the 2022 electricity baseline inventory. Each zip archive stores the source files as found in their data stores. Sub-folders in each of the data stores are archived separately. For example, stewi.zip contains the JSON files, while stewi.facility.zip is the 'facility' sub-folder of stewi data store that stores the parquet files. To reproduce the data store, extract each zip file and drag-and-drop sub-folders in to their appropriate root folders to recreate the data stores, then copy the root folders to your data store folder (as returned by running the following on the command line: `python -c "import appdirs; print(appdirs.user_data_dir())"`). The main five data stores include: 'electricitylci', 'facilitymatcher', 'fedelemflowlist', 'stewi', and 'stewicombo'. The log file generated by the 2022 model run is also included, which contains the statements at the DEBUG level and above.
The ElectricityLCI v2 Python package (https://github.com/USEPA/ElectricityLCI/tree/v2.0) was used to generate the 2021 electricity baseline: a regionalized life cycle inventory model of U.S. electricity generation, consumption, and distribution using standardized facility and generation data. ElectricityLCI implements a local data store for downloading and accessing public data on an individual's computer. The data store follows the folder definition provided by USEPA's esupy Python package (https://github.com/USEPA/esupy), which utilizes the appdirs Python dependency (https://pypi.org/project/appdirs/). An overview of the ElectricityLCI data stores may be found on the README (https://github.com/USEPA/ElectricityLCI/blob/v2.0/README.md#data-store). This submission includes the background data used to generate the 2021 electricity baseline inventory. Each zip archive stores the source files as found in their data stores. Sub-folders in each of the data stores are archived separately. For example, stewi.zip contains the JSON files, while stewi.facility.zip is the 'facility' sub-folder of stewi data store that stores the parquet files. To reproduce the data store, extract each zip file and drag-and-drop sub-folders in to their appropriate root folders to recreate the data stores, then copy the root folders to your data store folder (as returned by running the following on the command line: python -c "import appdirs; print(appdirs.user_data_dir())"). The main five data stores include: 'electricitylci', 'facilitymatcher', 'fedelemflowlist', 'stewi', and 'stewicombo'. The log file generated by the 2021 model run is also included, which contains the statements at the DEBUG level and above.
The ElectricityLCI v2 Python package (https://github.com/USEPA/ElectricityLCI/tree/v2.0) was used to generate the 2020 electricity baseline: a regionalized life cycle inventory model of U.S. electricity generation, consumption, and distribution using standardized facility and generation data. ElectricityLCI implements a local data store for downloading and accessing public data on an individual's computer. The data store follows the folder definition provided by USEPA's esupy Python package (https://github.com/USEPA/esupy), which utilizes the appdirs Python dependency (https://pypi.org/project/appdirs/). An overview of the ElectricityLCI data stores may be found on the README (https://github.com/USEPA/ElectricityLCI/blob/v2.0/README.md#data-store). This submission includes the background data used to generate the 2020 electricity baseline inventory. Each zip archive stores the source files as found in their data stores. Sub-folders in each of the data stores are archived separately. For example, stewi.zip contains the JSON files, while stewi.facility.zip is the 'facility' sub-folder of stewi data store that stores the parquet files. To reproduce the data store, extract each zip file and drag-and-drop sub-folders in to their appropriate root folders to recreate the data stores, then copy the root folders to your data store folder (as returned by running the following on the command line: python -c "import appdirs; print(appdirs.user_data_dir())"). The main five data stores include: 'electricitylci', 'facilitymatcher', 'fedelemflowlist', 'stewi', and 'stewicombo'. The log file generated by the 2020 model run is also included, which contains the statements at the DEBUG level and above.
Overview -------------------------- This repository contains datasets from the manuscript **"Enhanced Generalizability to Deep-Learning Quantification of 3D Microstructural Characteristics through Microstructurally Aware Augmentation of Scarce Data"** (*William F. Kent, Rochan Bajpai, Rachel C. Kurchin, William K. Epting, Harry W. Abernathy, Paul A. Salvador. Submitted 2026*). The methods are also described in the dissertation **Data Intensive Analysis of Solid Oxide Cell Microstructures** (*Doctoral dissertation, Carnegie Mellon University, 2025*). The datasets here are trained convolutional neural network (CNN) models for predicting key microstructural properties of solid oxide cell (SOC) electrodes from low-res, 2-channel 3D images, as well as some helpful code. The parameters for input images are provided in the paper. Sample data is provided in the file `Combined_anode_aug_dual_1k_examples` - that particular data was used to train `anode_all_aug.pth` and will work most accurately with that model. Please familiarize yourself with all caveats on accuracy and applicability, as detailed in the associated paper. Usage -------------------------- The basic usage is as follows, assuming `model_fn` is the path to the .pth file, and `X` is 2-channel input image(s) of the proper dimensions (either one image of shape `[2,12,24,24]`, or a batch of N input images of shape `[N,2,12,24,24]`): from CNN_inferencer import load_model_for_inference model = load_model_for_inference(model_fn) y_predicted = model(X) The model object automatically handles input scaling and output de-scaling based on the way the models were trained - in other words, pass in a 2-channel micro-CT image, and it will output microstructural property values in real units. ## Other model object attributes Note that model has useful attributes other than its forward pass model(X). * `model.output_descaler` - returns the output descaler object. Model does the de-scaling when generating inferences, but you may want to re-use this de-scaler on other values to e.g. compare predictions to ground truth from already-scaled training data. * `model.prop_names` - Gives the property names of the predicted y values, in order. Only exists if there's an output scaler as part of the model object, which there will be in the models provided here. ## Usage with sample data Here is a short script to use with the included sample data. from CNN_inferencer import display_predictions, load_model_for_inference, calculate_mape, parity_plot import h5py import numpy as np model_fn = 'anode_all_aug.pth' data_fn = 'Combined_anode_aug_dual_1k_examples.h5' N_samples = 200 figure_outdir = '.' model = load_model_for_inference(model_fn) with h5py.File(data_fn,'r') as f: XX = f['X'] #These are the 2-channel 3D images yy = f['y'] #These are the ground-truth microstructural properties, but they have been scaled for training - need to de-scale below N = XX.shape[0] #How many images total in the input data file #Run inferences on N_samples random samples from XX. #Run in a batch, much more efficient than one at a time. ii = np.random.choice(N,N_samples,replace=False) ii.sort() y_pred = model(XX[ii]) #Get the original/true (but normalized/scaled) values from the training dataset... #Because they were normalized, they are not in real units yet. So let's also de-scale them using model.output_scaler. y_true = model.output_scaler.transform(yy[ii]) #Let's display actual values for just 5 random ones for i in np.random.choice(N_samples,5,replace=False): display_predictions(y_true[i], y_pred[i], model.prop_names) #Make parity plots for each property (ground truth vs predicted values) #Also label each plot with the mean abs. percent error (MAPE) of the predicted values for i,key in enumerate(model.prop_names): mape = calculate_mape(y_true[:,i], y_pred[:,i]) parity_plot(y_true[:,i], y_pred[:,i], figure_outdir, key, extra_title=f' ({mape:.2f}% MAPE)')
After project completion on scope, on cost, and ahead of schedule, the Facility for Rare Isotope Beams began operations for scientific users in May of 2022. The ramp-up to a beam power of 400 kW is planned over a six-year period; 1 kW was delivered for initial user runs from in 2022, and 5 kW was delivered as of February 2023. Test runs with 10 kW 36Ar and 48Ca beams were conducted in July 2023. Upgrade plans include doubling the primary-beam energy to 400 MeV/nucleon for enhanced discovery potential (¿FRIB 400¿). This talk reports on the strategic plans towards high power operations emphasizing challenges and resolutions in beam-interception devices and targetry systems, radiation protection and controls, and legacy system renovation and integration.
A dual energy storage ring designed for beam cooling consists of two closed rings with significantly different energies: the cooling and damping rings. These two rings are connected by an energy recovering superconducting RF structure that provides the necessary energy difference. In our design, the RF acceleration has a main linac and harmonic cavities both running at crest that at first accelerates the beam from low energy E_{L} to high energy E_{H} and then decelerates the beam from E_{H} to E_{L} in the next pass. The purpose of the harmonic cavities is to extend the bunch length in a dual energy storage ring as such a longer bunch length may be very useful in a cooling application. Besides these cavities, a bunching cavity running on zero-crossing phase is used outside of the common beamline to provide the necessary longitudinal focusing for the system. In this paper, we present a preliminary lattice design along with the fundamental beam dynamics study in such a dual energy storage ring.
Nonlinear integrable optics (NIO) are a promising novel approach at improving the stability of high intensity beams. Implementations of NIO based on specialized magnetic elements are being tested at the Integrable Optics Test Accelerator (IOTA) at Fermilab. One method of verifying proper implementation of these solutions is by measuring the analytic invariants predicted by theory. The initial measurements of nonlinear invariants were performed during IOTA run in 2019/20, however the covid-19 pandemic prevented the full-scale experimental program from being completed. Several important improvements were implemented in IOTA for the 2022/23 run, including the operation at higher beam energy of 150 MeV, improved optics control, and chromaticity correction. This report presents on the improved calibrations of the NIO for nonlinear invariant measurements.
The nonlinear integrable optics studies at the integrable optics test accelerator (IOTA) demand fine control of the chromaticity using sextupole magnets. During the last experimental run undesirable misalignments and multipole composition in some sextupole magnets impacted operations. This report outlines the beam-based methods used to identify the nature of the misalignments and defects, and the subsequent magnetic measurements and remediation of the magnets for future runs.
The results of Run 24 experiments at Relativistic Heavy Ion Collider (RHIC) for improving luminosity using optics tuning are presented in this study. In the first experiment, MADx matching was used to output magnet strengths corresponding to specific s star movements around Interaction Region 8 (IR8). The corresponding Zero Degree Calorimeter (ZDC) signal was measured in place of luminosity, and Bayesian Optimization aids search of optimal movements. It was found that values retrieved from matching were inaccurate, resulting in negative feedback loops. The second experiment focused on calculating accurate s star movements. The matching method was replaced with a linear sensitivity matrix, directly relating optics to power supply, and its null space was used to fit constraints such as hysteresis effects. At the experiment, beam losses were observed at collimators around boundary of IR8, which were fixed for the third experiment. Dynamic mode decomposition was also introduced to improve quality of turn-by-turn (TBT) data as well as accuracy and consistency of optics measurements at IR8. These improvements will be tested in the experiment of next RHIC run for luminosity optimization.
The Integrable Optics Test Accelerator is a research electron and proton storage ring recently commissioned at the Fermilab Accelerator Science and Technology facility. Its research program is focused on testing novel techniques for improving beam stability and quality, notably the concept of non-linear integrable optics. In this paper, we report on run 1 results of experimental studies of a quasi-integrable transverse focusing system with one invariant of motion, a Henon-Heiles type system implemented with octupole magnets. Good agreement with simulations is demonstrated on key parameters of achievable tune spread, dynamic aperture, and invariant conservation. We also outline current simulation and hardware improvement efforts for run 2, planned for fall of 2019.
Understanding the interaction between humans and the Earth system is a computationally daunting task, with many possible approaches depending on resources available and questions of interest. For example, state-of-the-art impact models require decade-long time series of relatively high frequency, spatially resolved and often multiple variables representing climatic impact-drivers (Ruane et al., 2022). Most commonly these are derived from the outputs of detailed, computationally expensive Earth System Models (ESMs) run according to a standard, limited set of future scenarios, the latest being the SSP-RCPs run under CMIP6/ScenarioMIP (Eyring et al., 2016; O’Neill et al., 2016). At the time of writing, O’Neill et al. (2016) has been cited more than 1750 times and Eyring et al. (2016) more than 5000 times, highlighting the broad, general applications of this data. Often, however, impact modeling seeks to explore new scenarios that were not part of the ScenarioMIP protocol, and/or needs a larger set of initial condition ensemble members than are typically available to quantify the effects of ESM internal variability. In addition, the recognition that the human and Earth systems are fundamentally intertwined, and may feature potentially significant feedback loops, is making integrated, simultaneous modeling of the coupled human-Earth system increasingly necessary, if computationally challenging with most existing tools (Thornton et al., 2017). For impact modelers, climate model emulators can be the answer to meet both the needs of: 1) creating realizations for novel scenarios and 2) achieving a simplified, computationally tractable representation of ESM behavior in a coupled human-Earth system modeling framework. We proposed a new, comprehensive approach to such emulation of gridded, multivariate ESM outputs for novel scenarios without the computational cost of a full ESM, STITCHES (Tebaldi et al., 2022). The approach outlined in Tebaldi et al. (2022) should be extensible to future CMIP eras, although the STITCHES software at present is strictly focused on CMIP6/ScenarioMIP data hosted on Pangeo (https://gallery.pangeo.io/repos/pangeo-gallery/cmip6/). The corresponding STITCHES Python package uses existing archives of ESMs’ scenario experiments from CMIP6/ScenarioMIP to construct gridded, multivariate realizations of new scenarios provided by reduced complexity climate models (Hartin et al., 2015; Meinshausen et al., 2011; Smith et al., 2018), or to enrich existing initial condition ensembles. Its output provides the same characteristics as the emulated ESM output: multivariate (spanning potentially all variables that the ESM has saved), spatially resolved (down to the native grid of the ESM), and preserving the same high frequency as the original data. A new realization of multiple variables can be generated on the order of minutes with STITCHES, rather than the hours or sometimes days that ESMs require.
Elliptic partial differential equations (PDEs) can model many physical phenomena, such as electrostatics, acoustics, wave propagation, and diffusion. In scientific machine learning settings, a high-throughput PDE solver may be required to generate a training dataset, run in the inner loop of an iterative algorithm, or interface directly with a deep neural network. To provide value to machine learning users, such a PDE solver must be compatible with standard automatic differentiation frameworks, scale efficiently when run on graphics processing units (GPUs), and maintain high accuracy for a large range of input parameters. We have designed the jaxhps package with these use-cases in mind by implementing a highly efficient and accurate solver for elliptic problems with native hardware acceleration and automatic differentiation support.
The ATLAS experiment at CERN’s Large Hadron Collider (LHC) will be upgraded in two stages to prepare first for the Run~3 data-taking campaign, during which the integrated luminosity will roughly double, and then for the High-Luminosity LHC program with an ultimate integrated luminosity of up to 4 ab ^{{-1}} − 1 . These upgrades and high-statistics datasets will allow ATLAS to perform searches and precision measurements to constrain the Standard Model in yet-unexplored phase spaces and in the Higgs sector. This contribution summarizes the ATLAS detector upgrades and selected physics prospects for Run~3 and the HL-LHC.
This document describes the requirements for a software tool that will enable FRMAC to simulate large sets of sample result data that is based realistically on simulated radionuclide deposition grids from NARAC. The user of this tool would be scientists involved in exercise and drill planning or part of the simulation cell of an exercise controller team. A key requirement is that this tool must be able to be run with a reasonable amount of training and job aids by any person within the Assessment, Laboratory Analysis, or Monitoring and Sampling divisions of the FRMAC to support any level of exercise from the small IPX to the national level full scale exercise. This tool should be relatively lean and stand-alone so that the user can run it in the field with limited IT resources. This document will describe the desired architecture, design characteristics, order of operations, and algorithms that can be given to a software development team to assist them in project scoping, costing, and eventually, development.
The research program here combines Generic and Directed R\&D for Liquid Noble Gas detectors for neutrino physics and dark matter. In neutrino physics and dark matter research, we are addressing some of the most fundamental questions in particle physics today by studying these tiny particles in the electron family. To study them, detectors must be large, high precision, and ideally have the ability to collect detailed information from both charge and light from neutrino and dark matter interactions. Precision detection combined with large scales can be challenging. Liquid Argon detectors mitigate some of these challenges given the nature of its interaction medium and it’s relatively inexpensive cost per ton. Significant progress has been made in the last 15 years in developing these kinds of detectors to be built and operated at large scales. However there are still challenges and new ideas in moving these detectors from ton scale to kiloton scale. To be able to get the most out of these kinds of detectors directed questions with respect to specific components running and colleting data in the detectors must be understood. As well, new ideas on how to best combine charge and light measurements may lead to new ways to learn new things with these detectors. Experimental test stands to conduct this work, as is done in this proposal, to accomplish both of these goals are critical to address both of these questions. In the test stand enabled under this grant at Yale University’s Wright lab, new developments in efficient small scale setups were developed, instrumentation developed for running and future experiments (the SBND experiment at Fermilab at present, the DUNE experiment in the future) were and will be tested, and new ideas for charge and light determination for new measurements are under study. The long term impact of this work in neutrino physcis is both for the short and long baseline programs originating at Fermi National Accelerator Laboratory. At long baselines, the US flagship DUNE experiment will measure neutrino properties through neutrino oscillations using a beam originating at Fermilab and a massive LArTPC detector sited about a mile underground at the Sanford Underground Research Facility in Lead, South Dakota. Studies here on the high voltage for thes detectors, and ongoing work on charge and light production and collection in the TPC are relevant for the design and data taking of this and other LArTPC detectors to enable them to best measure the neutrinos created at Fermilab that then pass through the detector in South Dakota. What we learn in these detectors may ultimately change the course of how we understand neutrino’s impact in the early universe. What we learn about how to improve this detection technology and develop new ideas in Lar detectors in general may impact both future neutrino experiments and dark matter experiments.
The overall technical objective in phase I for the Surface Emitting Laser-Electro-Absorption Modulator (EAM) was to demonstrate feasibility of this approach by designing and building the laser-EAM chips, optimizing the fabrication process for the Laser/EAM in a series of 3 fabrication runs, irradiating Laser/EAM chips from each run and DC testing of the chips post radiation. 100 Gbps/channel testing would be done at the end of the phase II program before and after radiation. The overall technical objective for the driver circuit was to demonstrate the radiation tolerance of a circuit using the same technology as the driver, and to design the appropriate driver circuit.
This report provides a summary of notes for building and running the Sandia Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) code. SCEPTRE is a general purpose C++ code for solving the Boltzmann transport equation in serial or parallel using unstructured spatial finite elements, multigroup energy treatment, and a variety of angular treatments including discrete ordinates and spherical harmonics. Either the first-order form of the Boltzmann equation or one of the second-order forms may be solved. SCEPTRE requires a small number of open-source Third Party Libraries (TPL) to be available, and example scripts for building these TPL's are provided. The TPL's needed by SCEPTRE are Trilinos, boost, and netcdf. SCEPTRE uses an autoconf build system, and a sample configure script is provided. Running the SCEPTRE code requires that the user provide a spatial finite-elements mesh in Exodus format and a cross section library in a format that will be described. SCEPTRE uses an xml-based input, and several examples will be provided.