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.
CoSMo is a comprehensive shell model code suite designed for nuclear shell model and configuration interaction calculations. Its highly structured and templated nature allows for flexibility and ease in applications, including those to open quantum systems with non-Hermitian Hamiltonians, clustering, and time-dependent dynamics.
Abstract We present a data assimilation package for use with ocean circulation models in analysis, forecasting, and system evaluation applications. The basic functionality of the package is centered on a multivariate linear statistical estimation for a given predicted/background ocean state, observations, and error statistics. Novel features of the package include support for multiple covariance models, and the solution of the least squares normal equations either using the covariance matrix or its inverse—the information matrix. The main focus of this paper, however, is on the solution of the analysis equations using the information matrix, which offers several advantages for solving large problems efficiently. Details of the parameterization of the inverse covariance using Markov random fields are provided and its relationship to finite-difference discretizations of diffusion equations are pointed out. The package can assimilate a variety of observation types from both remote sensing and in situ platforms. The performance of the data assimilation methodology implemented in the package is demonstrated with a yearlong global ocean hindcast with a 1/4° ocean model. The code is implemented in modern Fortran, supports distributed memory, shared memory, multicore architectures, and uses climate and forecasts compliant Network Common Data Form for input/output. The package is freely available with an open source license from www.tendral.com/tsis/ .
Dynamic and adaptive mesh refinement is pivotal in high-resolution, multi-physics, multi-model simulations, necessitating precise physics resolution in localized areas across expansive domains. Today’s supercomputers’ extreme heterogeneity presents a significant challenge for dynamically adaptive codes, highlighting the importance of achieving performance portability at scale. Our research focuses on astrophysical simulations, particularly stellar mergers, to elucidate early universe dynamics. Here, we present Octo-Tiger, leveraging Kokkos, HPX, and SIMD for portable performance at scale in complex, massively parallel adaptive multi-physics simulations. Octo-Tiger supports diverse processors, accelerators, and network backends. Experiments demonstrate exceptional scalability across several heterogeneous supercomputers including Perlmutter, Frontier, and Fugaku, encompassing major GPU architectures and x86, ARM, and RISC-V CPUs. Parallel efficiency of 47.59% (110,080 cores and 6880 hybrid A100 GPUs) on a full-system run on Perlmutter (26% HPCG peak performance) and 51.37% (using 32,768 cores and 2048 MI250X) on Frontier are achieved.
Use of cellulosic biofuels from non-feedstocks are modeled using the BEPAM (Biofuel and Environmental Policy Analysis Model) model to quantifying the uncertainties about induced land use change effects, net greenhouse gas saving potential, and economic costs. The code is in GAMS, general algebraic modeling language. NOTE: Column 3 is titled "BAU" in "merged_BAU.gdx", "merged_RFS.gdx", and "merged_CEM.gdx", but contains "RFS" data in "merged_RFS.gdx" and "CEM" data in "merged_CEM.gdx".
The dataset consists of: (1) The replication codes and data for the BEPAM model are contained in the "BEPAM_Supplementary Environment Policy Analysis.zip" (2) Simulation results from the BEPAM model are contained in "ModelOutputs.zip" under the "BEPAM_Supplementary Environment Policy Analysis.zip"
This dataset provides raw millikelvin scanning tunneling microscopy/spectroscopy (STM/S) grid spectroscopy data and Python analysis scripts supporting the manuscript “Deciphering Majorana Zero Modes in Topological Superconductor FeTe0.55Se0.45 with Machine-Learning-Assisted Spectral Deconvolution.” The dataset includes a raw grid spectroscopy file acquired on FeTe0.55Se0.45 at 40 mK under magnetic field, together with Python/Jupytext analysis scripts used for STM/S data processing, visualization, spectral deconvolution, Lorentzian peak fitting, feature extraction, machine-learning-assisted clustering, and figure generation. These files support the analysis of vortex-core local density of states and the identification of zero-bias-peak-related spectral components from complex in-gap states. The dataset is intended to provide a citable archival record of the data and analysis code associated with the published manuscript and to support transparency and reproducibility of the reported STM/S and machine-learning workflow.
This paper discusses two cross section interpolation methods commonly found in popular nodal codes; the partial derivatives and multiple tables models. The motivation for choosing a model, and thus a case matrix structure, is a trade off between accuracy and computational cost. Due to decades of experience, there are default structures that are sufficient for current light water reactor analysis. However, this is not necessarily the case for advanced reactor designs. Therefore, it is advantageous to understand the sources of error in cross section interpolation models so that the quality of a case matrix may be improved. A mathematical framework for these models is presented in this work that provides a more rigorous connection between the nuclear engineering field's cross section interpolation methods and the broader mathematical field of function approximation. The two cross section models examined in this paper were found to utilize Lagrange interpolation and are a subset of Lagrange tensor products. Classical results of Lagrange polynomial error analysis were then applied to the partial derivative and multiple tables models to derive expressions for the total point-wise error. The analytical results classify the total error into two parts: the model form error and interpolation error. Finally, based on our observations, a better foundation for improving the quality of a case matrix is proposed. (authors)
Frequency-modulation spectroscopy (FMS) is generally suited to code-division multiplexing, and we demonstrate that capacity in a form of continuous-wave LIDAR, utilizing a sharp CO 2 absorption transition at 1.6 µm in simple ranging setups. The approach retains the advantages of FMS, including coherent detection and good rejection of broad absorption backgrounds. Extensions of this multiplexed approach to the continuous, simultaneous detection of several transitions would come by transmitting an encoded combination of frequency-modulated carriers, each tuned to detect a unique absorption transition. Signal analysis at the receiver involves a simple process of de-multiplexing that, in a general application, reveals targets at various distances and the absorption-related FMS signals in between.
High-resolution phase-contrast 3D imaging using nano-holotomography typically requires collecting multiple tomograms at varying sample-to-detector distances, usually 3 to 4. This multi-distance approach limits temporal resolution, making it impractical for operando studies. Moreover, shifting the sample complicates reconstruction, requiring precise alignment, registration, and interpolation to correct for shift-dependent magnification on the detector. In response, we propose and validate through simulations a novel, to the best of our knowledge, single-distance approach that leverages coded apertures to structure beam illumination while the sample rotates. Finally, this approach relies on a joint reconstruction scheme, which integrates phase retrieval with 3D tomography, ensuring data consistency and achieving artifact-free reconstructions from a single distance, unlocking dynamic experiments.
Data, code, and figures from our meta-analysis of the effect of precipitation change on soil respiration, with an emphasis on ecosystem variability and the effect of study duration. We combined the results of 81 separate studies to determine the effect of altered rainfall on the amount of carbon that leaves soil in the form of carbon dioxide, a major component of global carbon cycling. In addition we looked at how long the changes lasted, and how different soil properties and the intensity of precipitation changes affected the study results. We found that more precipitation resulted in greater amounts of carbon dioxide leaving the soil, and less precipitation resulted in less. However, the changes weakened over time in ecosystems that typically receive plenty of rainfall (e.g., forests), in contrast to ecosystems that typically receive little rainfall (e.g., deserts) where changes strengthened over time. Changes in the amount of carbon dioxide that left the soil were also affected by the amount of carbon in the soil, which impacts how much water soil can hold.
This dataset serves as supplementary information for the paper by Ciulla F. and Varadharajan C. A Network Approach for Multiscale Catchment Classification using Traits (see reference 1). It contains environmental and physical catchment traits, such as temperatures, precipitation, land use and human interference, from 9067 sites across the contiguous United States (CONUS). The purpose of this dataset is to provide information for a better trait-based categorization of river catchments in the CONUS using networks as an analytical tool. The traits variables match the ones present in the GAGES-II dataset and the preprocessing steps are described in the Methods section (processed_dataset.csv). Additionally we include the topologies (nodes, edges and clusters, also referred as classes) of the catchment network and traits network generated by said dataset (csv and json files). A series of tables support the information carried by the network providing more detailed descriptions of cluster components (SI1.pdf). A summary of all the plots of clusters of catchments with at least 50 nodes is provided (SI2.pdf). The characteristic traits for each cluster of catchments is presented as z-score (traits_categories_zscores_per_catchment_class.csv). The link to the hydrological behavior of clusters of catchments is displayed by boxplots, each describing a particular river discharge index (SI3.pdf). Both csv and json files can be read by common text editors but the data contained into them can be better handled using programming languages like python and database oriented libraries like pandas. Pdf files can be read by any pdf reader software.[02-23-2024] Update: The code and datasets necessary to reproduce the results of the study are available as a zipped repository (code_datasets_catchments_similarity.zip).
Atmospheric water vapor pressure is an essential meteorological control on land surface and hydrologic processes. It is not as frequently observed as other meteorologic conditions, but often inferred through the August–Roche–Magnus formula by simply assuming dew point and daily minimum temperatures are equivalent or by empirically correlating the two temperatures using an aridity correction. The performance of both methods varies considerably across different regions and during different time periods; obtaining consistently accurate estimates across space and time remains a great challenge. We applied an interpretable Long Short-Term Memory (iLSTM) network conditioned on static, location specific attributes to estimate daily vapor pressure for 83 FLUXNET sites in the United States and Canada. This data package includes all raw data of the 83 FLUXNET sites, input data for model training/validation/test, trained models and results, and python codes for the manuscript "Improving the Estimation of the Atmospheric Water Vapor Pressure Using an Interpretable Long Short-term Memory Network". Specifically, it consists of five parts. - First, "1_Daymet_data_83sites.zip" includes raw data downloaded from Daymet for the 83 sites used in the paper according to their longitude and latitude, in which vapor pressure is used. It also includes a pre-processed CSV data file combining all data from the 83 sites which is specifically used for the paper. - Second, "2_Fluxnet2015_data_83sites.zip" includes raw half hourly data of the 83 sites downloaded from FLUXNET2015 data portal, pre-processed daily data of the 83 sites, a CSV file including combined pre-processed daily data of the 83 sites, and a CSV file including the information (site ID, site name, latitude, longitude, data available period) of the 83 sites. - Third, "3_MODIS_LAI_data_83sites_raw.zip" includes raw leaf area index (LAI) data downloaded from the AppEEARs data portal. - Fourth, "4_Scripts.zip" includes all scripts related to model training and post-processing of a trained model, and a jupyter notebook showing an example for model post-processing. Two typo errors in files titled "run2get_args.py" and "postprocess.py" were corrected on March 27, 2024 to avoid confusions. - Finally, "Trained_models_and_results.zip" includes three folders and three files with suffix ".npy", and each folder corresponds to one file with suffix ".npy" with the same title. Each of the three folders include all trained models associated with one iLSTM model configuration (35 models for each configuration, details are described in the paper). Each file with suffix ".npy" includes the post-processed results of the corresponding 35 models under one iLSTM model configuration.
This dataset and R script accompany the published paper Bruns et al. (2024) in Geophysical Research Letters. The data are from the first six years of a field manipulation of whole-ecosystem warming and elevated CO2 experiment (Salt Marsh Accretion Response to Temperature eXperiment, or SMARTX) in the Smithsonian's Global Change Research Wetland (GCReW), a brackish, microtidal wetland site on a subestuary of the Chesapeake Bay. These data were generated to understand how warming and elevated CO2 interact to structure ecosystem-level responses to global change, particularly in terms of carbon sequestration. The dataset covers 2017-2022 and includes peak annual above ground biomass, annual belowground fine root productivity, and porewater NH4 for each experimental plot. The overall experiment is replicated in two locations on the marsh, a lower elevation zone dominated the C3 sedge S. Americanus and a higher elevation plot dominated by the C4 species. This paper and its data release is only for the C3 plot. Variable descriptions for data file is available in variable_descriptions.pdf. The R script Bruns_et_al_2024_GRL_make_figures.Rmd contains model code and other scripts used to generate all paper figures.
This data package provides inputs, codes, and outputs for a comprehensive analysis of projected changes in extreme precipitation across 10 regions of the continental United States, using 34 downscaled Earth System Models (ESMs) from the NASA Earth Exchange Global Daily Downscaled Projections, Coupled Model Intercomparison Project Phase 6 (NEX-GDDP-CMIP6) dataset. These models are part of the Coupled Model Intercomparison Project Phase 6 (CMIP6), a coordinated climate modeling framework widely used to assess climate change impacts. The analysis applies a stochastic storm transposition method to quantify changes in extreme rainfall under two Shared Socioeconomic Pathway (SSP) climate scenarios—SSP2-4.5 (moderate emissions) and SSP5-8.5 (high emissions)—compared to historical conditions (1995–2014 vs. 2081–2100). The dataset includes rainfall depth estimates for extreme events with return periods from 2 to 500 years across multiple storm durations (1, 3, and 5 days) for each of the 10 U.S. regions. Weighted ensemble statistics are derived from individual ESM performance against historical precipitation patterns, enabling robust uncertainty quantification through both sign-based and permutation-test-based model agreement assessments. Key analyses address: (1) relative changes in extreme precipitation for each climate scenario, (2) differences between SSP scenarios (SSP5-8.5 vs. SSP2-4.5), (3) contrasts between rare and frequent events, and (4) variations between multi-day and daily storm durations. The workflow produces ensemble statistics—median, 5th, 25th, 75th, and 95th percentiles—along with model agreement metrics that identify regions and event types with robust climate change signals. The dataset includes: processed rainfall depth outputs (netCDF format) from the RainyDay Python package, ESM weights from historical performance evaluation using DayMet observations, ensemble statistics across all storm dimensions, and figures summarizing key findings.
This data and code package supports the analysis presented in “Beyond Surface Cooling: Comprehensive Field Assessment of Reflective Pavement Thermal Performance in Phoenix, Arizona” and provides fully reproducible workflows for evaluating the thermal performance of cool pavement treatments in a hot urban environment. The dataset integrates multi-modal field measurements collected across residential and nonresidential settings, including mobile air temperature traverses, stationary air temperature monitoring, residential mean radiant temperature (MRT) measurements, subsurface temperature profiles, and controlled testbed observations. The data package contains raw and processed datasets in comma-separated value (CSV) format, accompanying metadata files describing site characteristics and measurement protocols, and R scripts (.R files) used for data cleaning, time synchronization, spatial and temporal matching, quality control filtering, statistical comparison, and figure generation. All analyses were conducted using R (version ≥ 4.2.0) with commonly available packages (e.g., tidyverse, lubridate, data.table, ggplot2). No proprietary software is required to reproduce results. Field campaigns were designed to quantify the effects of high-reflectance pavement coatings on surface temperature, near-surface air temperature, subsurface heat propagation, and radiative heat exposure. Temporal alignment procedures include standardized timestamp conversion and nearest-neighbor matching of high-frequency sensor measurements to stop-based metadata within defined tolerance windows to ensure comparability across instruments. The workflows generate summary statistics, treatment–control contrasts, depth-dependent thermal gradients, and time-series visualizations used in the associated publication. By integrating mobile, stationary, radiative, and subsurface measurements within a unified and transparent processing framework, this package enables comprehensive evaluation of cool pavement performance across multiple thermal exposure pathways and supports reuse in future urban heat mitigation and climate resilience studies.
The purpose of this data was to generate a scientific article that describes the responses of tropical roots to a warming experiment and to the effect of two consecutive hurricanes in Puerto Rico (Yaffar et al. in review). This data is from 10 months of minirhizotron images taken every 2 weeks at the experimental warming Tropical Responses to Altered Climate Experiment (TRACE) plots in Puerto Rico before and after Hurricanes Irma and Maria. This project has 3 warmed plots (plot 2,4,6) and 3 control plots (plot 1,3,5) with 2 minirhizotron tubes at each plot. As part of the data, there is also root data taken from cores and in-growth cores. Additionally, there is soil nutrient concentration data, soil microclimate, total leaf area, and canopy openness taken by Reed et al. 2020, and the TRACE census. Data files are in CSV format and the R code included in this package can be used with R 3.4.4 (R Core Team and contributors worldwide).
Jacketed Superconducting Radio Frequency (SRF) cavities structurally comprise of an inner niobium vessel surrounded by a liquid helium containment vessels. The pressure of the helium bath and/or its volume might be such that a jacketed SRF cavity shall be considered a system of pressure vessels. Thus, methods described in the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code (BPVC) should be used to analyze the structural soundness of jacketed SRF cavities. This paper will report the use of the set of rules developed at Fermilab for the design of SRF cavities, such as jacketed 1.3 GHz cavities for LCLS-II HE and jacketed Single Spoke Resonator type~2 (SSR2) for PIP-II, to ensure a similar level of safety as prescribed by the ASME BPVC.
We tested the performance of the HM-5 units and the NaIGEM code to determine the uranium enrichment. The uranium materials consist of two groups, the NBS series and the UISO series. The NBS items are U 3 O 8 in aluminum containers with an inner diameter of 7.0 cm and a frontal surface thickness of 2.0 mm. The UISO items are U 3 O 8 stored in doubly encapsulated stainless steel containers with a total thickness of 0.82 mm. The diameter of the inner UISO container is 10.0 cm. Four separate tests were performed.