SEARCH · Search NASA
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.
The objective of RAPIDS2 SciDAC Institute for Computer Science and Data is to assist Office of Science (SC) application teams in overcoming computer science and data challenges in the use of DOE supercomputing resources to achieve science breakthroughs.
The quantum computing market is entering a rapid growth phase. Growth was originally driven by private investment interest in high value applications such as drug discovery, materials science, and finance. Currently the main market drivers include technical advances, continued private investment, and rising public funding. The quantum computing market is poised for explosive growth — with its market size estimated at $1.42 billion USD in 2024, and an expected compound annual growth rate of 20.5% from 2025 to 2030 [4]. Superconducting qubits are the largest qubit type by revenue. However, alternative approaches such as trapped ion, neutral atom, spin based, and photonic qubits are also attracting investment.
The renaissance of nuclear energy has arrived, heralding an age of abundant inexpensive clean energy, and renewed space exploration. In nuclear-powered spacecraft and microreactors, safety and size are of utmost importance. Yttrium Hydride (YHx) is being researched for its utility as a neutron moderator in nuclear reactors; the hydrogen in YHx slows down neutrons, enabling a continuous nuclear reaction in the reactor. This has the benefit of allowing reactors to be more safe, compact, and efficient. The goal of this effort is to computationally predict the coefficient of temperature-dependent hydrogen diffusion within YHx, the Soret coefficient. This parameter is essential for determining the safe operating modes of YHx moderators. Zirconium Hydride (ZrHx) is used in the Training, Research, Isotopes, General Atomics (TRIGA) reactor, is the reference material for these calculations. In this work, nanoscale atomic modeling in the Vienna Ab initio Simulation Package (VASP) is combined with the mesoscale finite element phase-field module in the Multiphysics Object-Oriented Simulation Environment (MOOSE); this culminates in a new multiscale computational method to simulate Soret diffusion of hydrogen in YHx. This data is useful for predicting experimental outcomes. This workflow involves convergence testing followed by static, Nudged Elastic Band (NEB), Quasi-Harmonic Approximation (QHA), and Molecular Dynamics (MD) calculations - linked with phase field simulation. NEB simulates hydrogen migration, while QHA and MD predict temperature-dependent properties. The static calculations align with literature, and preliminary NEB and QHA calculations yield accurate results. Once the atomic calculations are complete, we will incorporate Electron Backscatter Diffraction (EBSD) images and VASP-generated parameters into the phase field module to simulate intra- and intergranular transport of hydrogen in ZrHx and YHx. Future research will extend our approach to fuel-moderator materials systems such as Uranium-Yttrium Hydride (U-YHx). This work contributes to the development of advanced nuclear energy solutions for space travel.
Omega version 0, the first step toward a new ocean model, has been developed to take full advantage of next-generation supercomputers, including those powered by GPUs. Omega is rebuilt from the ground up, redesigning computational aspects of ocean models to run faster and more efficiently across diverse high-performance computing systems.
While many researchers can design knockdown and knockout methodologies to remove a gene product, this is mainly untrue for new chemical inhibitor designs that empower multifunctional DNA Damage Response (DDR) networks. Here, we present a robust Goldilocks (GL) computational discovery protocol to efficiently innovate inhibitor tools and preclinical drug candidates for cellular and structural biologists without requiring extensive virtual screen (VS) and chemical synthesis expertise. By computationally targeting DDR replication and repair proteins, we exemplify the identification of DDR target sites and compounds to probe cancer biology. Our GL pipeline integrates experimental and predicted structures to efficiently discover leads, allowing early-structure and early-testing (ESET) experiments by many laboratories. By employing an efficient VS protocol to examine protein-protein interfaces (PPIs) and allosteric interactions, we identify ligand binding sites beyond active sites, leveraging in silico advances for molecular docking and modeling to screen PPIs and multiple targets. A diverse 3,174 compound ESET library combines Diamond Light Source DSI-poised, Protein Data Bank fragments, and FDA-approved drugs to span relevant chemotypes and facilitate downstream hit evaluation efficiency for academic laboratories. Two VS per library and multiple ranked ligand binding poses enable target testing for several DDR targets. This GL library and protocol can thus strategically probe multiple DDR network targets and identify readily available compounds for early structural and activity testing to overcome bottlenecks that can limit timely breakthrough drug discoveries. By testing accessible compounds to dissect multi-functional DDRs and suggesting inhibitor mechanisms from initial docking, the GL approach may enable more groups to help accelerate discovery, suggest new sites and compounds for challenging targets including emerging biothreats and advance cancer biology for future precision medicine clinical trials.
Packed columns are commonly used in post-combustion processes to capture CO 2 emissions by providing enhanced contact area between a CO 2 -laden gas and CO 2 -absorbing solvent. To study and optimize solvent-based post-combustion carbon capture systems (CCSs), computational fluid dynamics (CFD) can be used to model the liquid–gas countercurrent flow hydrodynamics in these columns and derive key determinants of CO 2 -capture efficiency. However, the large design space of these systems hinders the application of CFD for design optimization due to its high computational cost. In contrast, data-driven modeling approaches can produce fast surrogates to study large-scale physics problems. We build our surrogates using MeshGraphNets (MGN), a graph neural network framework that efficiently learns and produces mesh-based simulations. We apply MGN to a random packed column modeled with over 160K graph nodes and a design space consisting of three key input parameters: solvent surface tension, inlet velocity, and contact angle. Our models can adapt to a wide range of these parameters and accurately predict the complex interactions within the system at rates over 1700 times faster than CFD, affirming its practicality in downstream design optimization tasks. This underscores the robustness and versatility of MGN in modeling complex fluid dynamics for large-scale CCS analyses.
Since 2021, NREL has been developing a system to measure heliostats by measuring the deflection of printed computer vision targets, called the Reflected Target Non-intrusive Assessment (ReTNA) [6], [7]. While this system will have lower resolution than a fringe deflectometry system, it has several important advantages that make it a complimentary technology: 2D surface slope measurement can be generated from a single image, it can operate in ambient lighting, target points can be directly located in 3D space with photogrammetry allowing for a non-flat target, and it's well suited to using a smaller target, and multiple images to measure larger optical surfaces. ReTNA has undergone several significant changes and improvements, described below. This talk will summarize new system layouts designed for commercial use, new computer vision algorithms used to automate the analysis process and validation campaigns for the ReTNA software.
Knowledge of the transport properties of biomass particles such as porosity, tortuosity, and permeability is paramount for high-fidelity modeling of biomass pyrolysis due to the heat and mass transfer limitations imposed by particle microstructure. X-ray computed tomography (XCT) is a non-destructive imaging method that enables full 3D reconstructions of the biomass particle microstructure with high resolution, permitting direct calculation of porosity, tortuosity, and permeability from real particle geometries. In this study, XCT imaging revealed the 3D microstructures of particles and chars from pyrolytic conversion of cylindrically cut or milled/pelletized loblolly pine samples. The porosity, tortuosity, and permeability were calculated directly from the XCT geometries via open-source microstructural analysis tool MATBOX+TauFactor (https://github.com/NREL/MATBOX_Microstructure_analysis_toolbox) and computational fluid dynamics (CFD) simulations using our solver, Mesoflow (https://github.com/NREL/mesoflow). These properties were used in a reactor scale model developed in COMSOL of the single particle reactor at NREL to investigate the impact of feedstock pre-processing on biomass conversion during pyrolysis with rigorous experimental validation.
It is desirable for fourth-generation Small Modular Reactors to be passively cooled in standard and accident operations. Passive Containment Cooling Systems can reject heat from the containment structure, without using pumps or blowers. The targeted design containment structure is a large, domed, stainless steel, cylindrical vessel. In a postulated Design Basis Accident, steam will flash inside containment. Steam condensation occurs on the inner containment wall and transfers heat through the steel containment into a large body of water known as the annular reservoir (AR) surrounding the vessel serving as the ultimate heat sink. Natural circulation drives the flow in the AR and heat will be released to the environment by evaporation of water. Unique containment geometry requires a separate effects test (SET) facility for the verification and validation of the computer code and evaluation model development and assessment for reactor licensing efforts. In this study, STAR-CCM+, a computational fluid dynamics (CFD) code was used to inform the decision-making process on the design of the SET. The CFD simulation modeled, a two-phase turbulent flow with fluid film development and heat transfer for different containment geometries. The Reactor Excursion and Leak Analysis Program will also be used in a code-to-code verification against the CFD results.
Quantum computers are expected to outperform classical computers for specific problems in quantum chemistry. Such calculations remain expensive, but costs can be lowered through the partition of the molecular system. In the present study, partition was achieved with range-separated density functional theory (RS-DFT). The use of RS-DFT reduces both the basis set size and the active space size dependence of the ground state energy in comparison with the use of wave function theory (WFT) alone. The utilization of pair natural orbitals (PNOs) in place of canonical molecular orbitals (MOs) results in more compact qubit Hamiltonians. To test this strategy, a basis-set independent framework, known as multiresolution analysis (MRA), was employed to generate PNOs. Tests were conducted with the variational quantum eigensolver for a number of molecules. The results show that the proposed approach reduces the number of qubits needed to reach a target energy accuracy.
Most observables at particle colliders involve physics at a wide variety of distance scales. Due to asymptotic freedom of the strong interaction, the physics at short distances can be calculated reliably using perturbative techniques, while long distance physics is non-perturbative in nature. Factorization theorems separate the contributions from different scales, allowing to identify the pieces that can be determined perturbatively from those that require non-perturbative information, and if the non-perturbative pieces can be reliably determined, one can use experimental measurements to extract the short distance effects, sensitive to possible new physics. Without the ability to compute the non-perturbative ingredients from first principles one typically identifies observables for which the non-perturbative information is universal in the sense that it can be extracted from some experimental observables and then used to predict other observables. In this paper we argue that the future ability to use quantum computers to calculate non-perturbative matrix elements from first principles will allow to make predictions for observables with non-universal non-perturbative long-distance physics.
Abstract The Idaho National Laboratory (INL) has implemented laboratory-based micro-X-ray computed tomography in a laboratory equipped for the examination of highly radioactive samples. This capability provides nondestructive three-dimensional volumetric information on samples to inform subsequent traditional destructive examinations as well as real-world inputs for high-fidelity scientific modeling. Samples can be imaged with spatial resolutions ranging from several hundred nm/voxel up to ~ 100 µm/voxel. The best usable spatial resolution achieved to date is 384 nm/voxel with this instrument, while the highest radiological dose rate of a sample imaged is ~ 60 R/h β/γ on contact. Advanced data analysis, including custom tomographic reconstruction and segmentation methods, have also been developed to support this capability. In addition to traditional digital X-ray radiography and tomography, this instrument is also able to visualize in situ tensile and compression testing as well as perform diffraction contrast tomography. This work describes the X-ray computed tomography post-irradiation examination capabilities at INL, as well as detailing a variety of applications this instrument has examined.
Multicomponent plastics cannot be processed using mechanical recycling technologies, hindering efforts to deal with plastic waste. Multicomponent plastics include multilayer plastic films, which are widely used for food and healthcare packaging. Multilayer films combine several layers (potentially dozens) of different polymers to protect products from external factors (e.g., oxygen, water, temperature, shock, and light). Solvent-based separation processes have emerged as a promising alternative to recycle these complex materials. For instance, the Solvent-Targeted Recovery and Precipitation (STRAP TM ) process uses sequential solvent washes to selectively dissolve and separate constituent polymers from multicomponent plastic waste, including films. STRAP TM process design (separation sequence, type of solvents, and operating conditions) changes significantly depending on the design of the multilayer plastic film (e.g., number, types, and proportions of polymers). The ability to quickly quantify the economic and environmental benefits of diverse STRAP TM process designs is essential to accelerate the development of sustainable recycling processes and more recyclable multilayer film products. In this work, we present a fast computational framework that integrates molecular-scale models, process modeling, and techno-economic and life cycle analysis to quickly evaluate STRAP TM designs. The computational framework is general and can be used to study the processing of complex multilayer plastic waste streams that contain many layers. Furthermore, we highlight the different uses of the framework via targeted case studies.
Although machine learning (ML) has emerged as a powerful tool for rapidly assessing grid contingencies, prior studies have largely considered a static grid topology in their analyses. This limits their application, since they need to be re-trained for every new topology. Here, this paper explores the development of generalizable graph convolutional network (GCN) models by pre-training them across a range of grid topologies and contingency types. We found that a GCN model with auto-regressive moving average (ARMA) layers with a line graph representation of the grid offered the best predictive performance in predicting voltage magnitudes (VM) and voltage angles (VA). We introduced the concept of phantom nodes to consider disparate grid topologies with a varying number of nodes and lines. For pre-training the GCN ARMA model across a variety of topologies, distributed graphics processing unit (GPU) computing afforded us significant training scalability. The predictive performance of this model on grid topologies that were part of the training data is substantially better than the direct current (DC) approximation. Although direct application of the pre-trained model to topologies that are not part of the grid is not particularly satisfactory, fine-tuning with small amounts of data from a specific topology of interest significantly improves predictive performance. In general, this paper highlights the feasibility of training large-scale GNN models to assess the reliability of power grids by considering a wide variety of grid topologies and contingency types. With the advent of foundational models in ML and the exponential increase in GPU computing clusters, generalizable ML models will significantly enhance how utilities manage power systems and make decisions in real-time or near-real-time.
Here we propose a nodal stochastic generation and transmission expansion planning model that incorporates the output from high-resolution global climate models through load and generation availability scenarios. We implement our model in Pyomo and perform computational studies on a realistically-sized test case of the California electric grid in a high performance computing environment. We propose model reformulations and algorithm tuning to efficiently solve this large problem using a variant of the Progressive Hedging Algorithm. We utilize the parallelization capabilities and overall versatility of mpi-sppy, exploiting its hub-and-spoke architecture to concurrently obtain inner and outer bounds on an optimal expansion plan. Initial results show that instances with 360 representative days on a system with over 8,000 buses can be solved to within 5% of optimality in under 4 h of wall clock time, a first step towards solving a large-scale power system expansion planning problem across a wide range of climate-informed operational scenarios.
The vertical structure of vegetation canopies creates micro-climates. However, the land components of most Earth System Models, including the Energy Exascale Earth System Model (E3SM), typically neglect vertical canopy structure by using a single layer big-leaf representation to simulate water, CO 2 , and energy exchanges between the land and the atmosphere. In this study, we developed a Multi-Layer Canopy Model for the E3SM Land Model to resolve the micro-climate created by vegetation canopies. The model developed in this study re-implements the CLM-ml_v1 to support heterogeneous computing architectures consisting of CPUs and GPUs and includes three additional optimization-based stomatal conductance models. The use of Portable, Extensible Toolkit for Scientific Computation provides a speedup of 25–50 times on a GPU relative to a CPU. The numerical implementation of the model was verified against CLM-ml_v1 for a month-long simulation using data from the Ameriflux US-University of Michigan Biological Station site. Model structural uncertainty was explored by performing control simulations for five stomatal conductance models that exclude and include the control of plant hydrodynamics (PHD) on photosynthesis. The bias in simulated sensible and latent heat fluxes was lower when PHD was accounted for in the model. Additionally, six idealized simulations were performed to study the impact of three environmental variables (i.e. air temperature, atmospheric CO 2 , and soil moisture) on canopy processes (i.e. net CO 2 assimilation, leaf temperature, and leaf water potential). Increasing air temperature reduced net CO 2 assimilation and increased air temperature. Net CO 2 assimilation increased at higher atmospheric CO 2 , while decreasing soil moisture resulted in lower leaf water potential.
In this paper, we introduce a nonlinear stochastic model to describe the propagation of information inside a computer processor. In this model, a computational task is divided into stages, and information can flow from one stage to another. The model is formulated as a spatially-extended, continuous-time Markov chain where space represents different stages. This model is equivalent to a spatially-extended version of the M/M/s queue. The main modeling feature is the throttling function which describes the processor slowdown when the amount of information falls below a certain threshold. We derive the stationary distribution for this stochastic model and develop a closure for a deterministic ODE system that approximates the evolution of the mean and variance of the stochastic model. In conclusion, we demonstrate the validity of the closure with numerical simulations.
Questions of ‘how best to acquire data’ are essential to modelling and prediction in the natural and social sciences, engineering applications, and beyond. Optimal experimental design (OED) formalizes these questions and creates computational methods to answer them. This article presents a systematic survey of modern OED, from its foundations in classical design theory to current research involving OED for complex models. We begin by reviewing criteria used to formulate an OED problem and thus to encode the goal of performing an experiment. We emphasize the flexibility of the Bayesian and decision-theoretic approach, which encompasses information-based criteria that are well-suited to nonlinear and non-Gaussian statistical models. We then discuss methods for estimating or bounding the values of these design criteria; this endeavour can be quite challenging due to strong nonlinearities, high parameter dimension, large per-sample costs, or settings where the model is implicit. A complementary set of computational issues involves optimization methods used to find a design; we discuss such methods in the discrete (combinatorial) setting of observation selection and in settings where an exact design can be continuously parametrized. Finally we present emerging methods for sequential OED that build non-myopic design policies, rather than explicit designs; these methods naturally adapt to the outcomes of past experiments in proposing new experiments, while seeking coordination among all experiments to be performed. Throughout, we highlight important open questions and challenges.