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.
Abstract We prove, under mild conditions, the convergence of a Riemannian gradient descent method for a hyperbolic neural network regression model, both in batch gradient descent and stochastic gradient descent. We also discuss a Riemannian version of the Adam algorithm. We show numerical simulations of these algorithms on various benchmarks.
Diffusion in NaFePO 4 can be enhanced through amorphization. Based on computations using DFT and machine learning potentials, we ascribe this phenomenon to the formation of less constrained Na-ion environments upon disordering. Most modern battery technologies depend on solid-state crystalline cathode materials. However, some of these materials are constrained by the low ionic conductivity of their most stable phases. An example of this is maricite (NaFePO 4 ). Interestingly, experiments have shown that maricite can improve its rate capability through disordering (amorphization). However, experimental characterization of amorphous cathode materials remains a major challenge, hindering a clear understanding of the structural origin of the disorder-induced improvement in sodium-ion mobility. To address this, we here employ molecular dynamics simulations by first training a machine learning potential for NaFePO 4 based on the atomic cluster expansion approach and a batch active learning potential parameterization scheme. This potential is then applied to explore the structural and dynamical properties of NaFePO 4 glasses as cathode materials. Specifically, we investigate the effect of glass structure on sodium-ion diffusion, revealing the relative influences of short-range and medium-range order features. We find significant heterogeneity in sodium-ion diffusivity in the glass, with fast-conducting ions residing in less constrained atomic environments with fewer P and Fe neighbors. These more mobile ions are also surrounded by larger ring-type structures. Overall, the results and developed approach present promising avenues for developing high-performance glassy cathodes for next-generation batteries.
Liquid-fuel Molten Salt Reactors (MSRs) allow for continuous online fuel treatment and processing using a variety of subsystems. Among these subsystems, the off-gas system (OGS) continuously removes fission products from the primary fuel salt through a helium sparging process. Accurately modeling fission product behavior is paramount for the development of MSRs; while there are a variety of depletion modeling and simulation tools that allow for removal of material from a defined system, most of them work in a semicontinuous, batch-wise manner. The new tools implemented into TRITON, in the SCALE suite for reactor analysis and design, aim to address this issue by implementing truly continuous material feeds and removals. In addition, material accountancy is enhanced in these tools, which enable tracking of removed material in user-specified mixtures. These continuous depletion tools provide a way to model material transport throughout systems more efficiently and accurately. This report validates the new SCALE MSR modeling capabilities using data from the Molten Salt Reactor Experiment (MSRE). After a framework was established for deriving realistic removal rates based on design specifications, these rates were used by TRITON to continuously remove specific elements while performing burnup calculations. Results were compared to those obtained in MSRE, with a focus on the 135 Xe poison fraction. After these capabilities were validated, the learned concepts and best practices were applied to model the Molten Salt Breeder Reactor (MSBR) and the Molten Salt Demonstration Reactor (MSDR) to study the effects of continuous fission product removal in fuel salt reactivity, as well as the lifetime effects.
TDA Research, Inc. has developed a low-cost sorbent (alkalized alumina) based technology for post-combustion CO 2 capture. The sorbent runs in an isothermal process for adsorption and regeneration. Neither temperature swing nor pressure swing is needed. TDA designed a system which operated 10 fixed beds to simulate a moving bed process. The flow pattern was optimized to minimize the steam usage for regeneration. An excellent sorbent (Batch 1) was also developed with collaboration with our partner. The sorbent degraded after long term operation. We developed a reprocessing procedure, which can be conducted in situ, to restore the sorbent capture capacity as the fresh sorbent after running for about 1250 hours. The reprocessing extended the sorbent longevity significantly. TDA designed a pilot unit (40’ x 32’ x 11.5’), processing flue gas equivalent to 0.5 MW of power generation. It had 10 beds and each bed could hold 1.75 m 3 of sorbent. We worked with Springs Fabrication, Inc. to complete the construction of the pilot unit. It was then shipped and installed at NCCC. The pilot unit was kept running in 24-hour mode when the flue gas and steam were available. It could run automatically without operator on site. Parametric and long-term (2 months) tests were carried out successfully. The results showed that TDA’s process can achieve 90% capture and 95% CO 2 purity for both coal and NG flue gases. For coal flue gas, the system reached performance target when processing up to 0.62 MW flue gas, 24% higher than the design capacity. Thus, capital cost could be saved on the reactors and sorbent. The strip air flow was designed to be 0.25 of that of the flue gas. The test data showed the strip/flue ratio can be reduced to as little as 0.18, which saves the power consumption for the strip air blower. The pressure drop was found to be lower than what the empirical equation calculated. The sorbent still had 91.7% of original capacity after 3-month test. With the data from the pilot test, we updated the TEA. For a 550-MW e supercritical coal fired power plant with CO 2 capture, the capture cost for TDA’s process is $34.9/tonne CO 2 captured, which meets DOE’s goal of $40/tonne and is 17.1% less than the DOE baseline Case 12. Therefore, TDA’s process has a good potential for commercialization.
Microbial production of succinic acid (SA) at an industrially relevant scale has been hindered by high downstream processing costs arising from neutral pH fermentation for over three decades. Here, we metabolically engineer the acid-tolerant yeast Issatchenkia orientalis for SA production, attaining the highest titers in sugar-based media at low pH (pH 3) in fed-batch fermentations, i.e. 109.5 g/L in minimal medium and 104.6 g/L in sugarcane juice medium. We further perform batch fermentation using sugarcane juice medium in a pilot-scale fermenter (300×) and achieve 63.1 g/L of SA, which can be directly crystallized with a yield of 64.0%. Finally, we simulate an end-to-end low-pH SA production pipeline, and techno-economic analysis and life cycle assessment indicate our process is financially viable and can reduce greenhouse gas emissions by 34–90% relative to fossil-based production processes. We expect I. orientalis can serve as a general industrial platform for production of organic acids.
Designing a superconducting qubit to realize specific Hamiltonian parameters typically requires iterating through a time and compute-intensive forward loop in which the designer chooses a layout geometry, simulates it, extracts circuit parameters such as capacitances, and refines the geometry. We study the inverse version of this task using a neural-network workflow that maps target Hamiltonian parameters directly to component-level layout parameters, which we subsequently demonstrate on a planar transmon layout. During training, we pair the inverse model with a frozen forward surrogate model and evaluate the loss in Hamiltonian space rather than in layout-parameter space. In validation against a conventional EM solver, 97% of generated designs produce usable geometries, and the inverse-plus-surrogate pipeline reaches mean percent errors of 0.73% for qubit frequency and 1.58% for anharmonicity, comparable to or below the fabrication and simulation-to-measurement uncertainty expected for academic-process transmon devices of this type. A single pipeline query takes ~60 ms on CPU, versus ~2 min for a conventional EM capacitance extraction on the same hardware, a speedup of approximately 2,000x. Batching minimizes the AI model inference overhead, reducing the runtime to 3.1 microseconds per sample on CPU and 2.6 microseconds per sample on GPU at a batch size of 2048, resulting in speedups of 3.9 x 10^7 and 4.6 x 10^7, respectively, relative to a single conventional CPU EM extraction. Our results indicate that component-level inverse design usefully extends and complements conventional EM simulation, including for small datasets on the order of 1,000 samples.
MCwrapper is a set of systems that manages the entire Monte Carlo production workflow for GlueX and provides standards for how that Monte Carlo is produced. MCwrapper was designed to be able to utilize a variety of batch systems in a way that is relatively transparent to the user, thus enabling users to quickly and easily produce valid simulated data at home institutions worldwide. Additionally, MCwrapper supports an autonomous system that takes user’s project submissions via a custom web application. The system then atomizes the project into individual jobs, matches these jobs to resources, and monitors the jobs status. The entire system is managed by a database which tracks almost all facets of the systems from user submissions to the individual jobs themselves. Users can interact with their submitted projects online via a dashboard or, in the case of testing failure, can modify their project requests from a link contained in an automated email. Beginning in 2018 the GlueX Collaboration began to utilize the Open Science Grid (OSG) to handle a bulk of simulation tasks; these tasks are currently being performed on the OSG automatically via MCwrapper. This talk will outline the entire system of MCwrapper, its use cases, and the unique challenges facing the system.
Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.
Double shell targets are an alternative ignition platform for inertial confinement fusion. One design consideration for double shell targets is the choice of inner shell material to help trap radiation emitted by the hot fuel to aid ignition. Materials such as molybdenum and tungsten are of interest for the inner shell layer of the targets. While molybdenum has a lower density that could inhibit instability growth and allow for radiography and code benchmarking, tungsten has a higher density that could provide more compression and confinement. These tradeoffs have been explored using optimized designs for each material. Our previous work [Vazirani et al., “Coupling 1D xRAGE simulations with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 28, 122709 (2021); Vazirani et al., “Coupling multi-fidelity xRAGE with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 30, 062704 (2023); and Vazirani et al., “Bayesian batch optimization for molybdenum versus tungsten inertial confinement fusion double shell target design,” Stat. Anal. Data Min. 17, e11698 (2024)] resulted in a multi-fidelity Bayesian optimization framework to find yield-optimized double shell target geometries. By leveraging simulations of varying fidelities (one-dimensional and two-dimensional) to inform one another, the multi-fidelity optimization was able to optimize a design in the highest fidelity with significantly fewer simulations than would be used in a systematic parameter scan. In this work, we apply the multi-fidelity Bayesian optimization to explore the optimized designs of double shell targets with molybdenum and tungsten inner shells as well as the physics producing the high performing implosions. A physics exploration of all the simulations used in this study shows trends in designs that contribute to high yields, ion temperatures, and fuel areal densities. Comparison of molybdenum and tungsten simulations shows that they can produce similar implosion conditions with different geometries, which would be important to study in experiments. Graded density layers produce varying performances with the two materials but continue to be of interest for future studies along with studies of doped inner shell materials and applied surface roughness.
The primary operational costs for existing nuclear reactors are plant operation costs, maintenance costs, and fuel costs, all of which are influenced by the materials used and the design of the reactor core. Optimizing core design parameters—including burnup limits and enrichment levels—can lengthen cycles, reduce outages, reduce reload batch fractions and spent fuel storage requirements, and lower maintenance and operating expenses, thereby enhancing economic viability. Furthermore, developing higher-fidelity tools to simulate these parameters enables better identification of the available margin, improves overall plant safety, and improves the understanding a given plant’s responses to accident scenarios. Here, in the US, much of the research and development focus has traditionally been on pressurized water reactors (PWRs), but boiling water reactors (BWRs) comprise approximately one-third of the US reactor fleet. Modeling and simulation advances for BWRs and PWRs—particularly those achieved through the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program—are crucial to the long-term viability of the light–water reactor industry. A key research area of the high burnup and increased enriched fuel initiative is focused on addressing issues related to postulated loss-of-coolant accident (LOCA) scenarios. NEAMS has dedicated significant effort to enhancing tools to better support BWRs. A current focus is showcasing the BWR framework for high-burnup LOCA analysis. This high-fidelity steady-state analysis is a first step toward demonstrating a best-estimate, pin-by-pin high-burnup BWR LOCA analysis to assess full-core cladding rupture behavior for a representative BWR. The objective of this effort is to provide a modeling capability that will help elucidate and provide a best-estimate evaluation for cladding rupture susceptibility in BWRs. This modeling capability could then be used to prevent and/or mitigate cladding ruptures in postulated accident scenarios without penalizing operational parameters. Additionally, the results of this work will help identify strategies for finding additional margins or potentially limiting cladding ruptures through core design optimizations to enable more efficient core designs.
Gamma Detector Response and Analysis Software (GADRAS) is used by the radiation detection and emergency response community to perform modeling and spectral analysis for gamma detector systems. Built into GADRAS is the ability to define a detector, geometry, background characteristics and source composition to generate synthetic spectra for drills and exercises (injects). Consequence Management is currently in development of a sample result data simulator tool in which a deposition model is probed for source conditions at moments in time and locations in space. These values are used to generate realistic sample results for use in drills and exercises. In addition to sample results, there is a need to simulate the actual spectra that would be observed in the field by downlooking HPGe instruments given a deposition activity. This way, the FRMAC Gamma Spectroscopist can practice their process of generating quantified results from spectra on realistic data as well. Recognizing the decades of work done in GADRAS to accurately generate synthetic spectra, this team decided to build a link between the new simulator and GADRAS to generate these spectra quickly and easily. The simulator tool will generate a file that specifies the name of the spectra, its location, date/time of measurement, duration of measurement, height off the ground, and the deposition activity and age for every radionuclide in the simulation. Then, a new tool within the Inject Tab of GADRAS was developed to read in this file given a detector selection and generate In-Situ spectra for each row in the file in any file format the user chooses. This way, simulation cell staff can take these files and then upload them to the appropriate data system (RAMS or RadResponder) for use during drills and exercises. An advanced feature of this tool allows for generating any spectra given an appropriate model and mapping of source to model layer in the batch inject tool. This way, spectra from field sample counts, mobile laboratories, or even fixed laboratories can be generated in bulk given an estimate of the radioactivity concentration or total radioactivity in an import file. This expands the capabilities of this tool a great deal and will make it a more useful tool for CM and others to help estimate detector response for nearly any situation. This user guide will explain the steps needed to perform a batch inject file generation.
Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.
Zeolites are one of the most widely used materials in the chemical industry due to their nanometer-sized pores that can adsorb and react upon molecules selectively. With hundreds of known framework topologies and hundreds of thousands of computationally predicted structures, the ability to rapidly predict zeolite performance allows researchers to prioritize their efforts on the most promising structures for a given application. Although the accuracy of forcefield-based atomistic simulations has advanced significantly in the past two decades, these simulations can be computationally expensive, especially for long-chain, complex molecules. Here, we present ZeoNet, a representation learning framework using convolutional neural networks (ConvNets) and 3D volumetric representations for predicting adsorption in zeolites. ZeoNet was trained on the task of predicting Henry's constants for adsorption, k H , of n-octadecane in more than 330 000 known and predicted zeolite materials. Employing a 3D grid based on the distances to solvent-accessible surfaces, a volumetric representation that can be generated efficiently, the best-performing ZeoNet achieved a correlation coefficient r 2 = 0.977 and a mean-squared error MSE = 3.8 in ln k H , which corresponds to an error of 9.3 kJ mol -1 in adsorption free energy. In comparison, a model based on hand-designed geometric features has values of r 2 = 0.783 and MSE = 35.7. ZeoNet is also relatively efficient and can process ≈8 structures per second on an Nvidia RTX 2080TI GPU, orders of magnitude faster than forcefield-based simulations. A systematic analysis was conducted to investigate how the choice of ConvNet architectures, the linear dimension (L) and spatial resolution (Δd) of the distance grids, batch size, optimizer, and learning rate impact the model performance. We found that ConvNets based on the ResNet architecture offer the best tradeoff between expressiveness and efficiency. The performance for all models reaches a plateau at L = 30–45 Å and depends less sensitively on grid resolution, with a small benefit around Δd = 0.30–0.45 Å. Finally, saliency maps were visualized to identify which regions of the materials contributed the most to model predictions. It was found, interestingly, that the predictions are driven primarily by the accessible pore volume rather than the region occupied by the framework atoms.
For solidification rates at or near equilibrium solidification conditions the effects of chrome (Cr) and nickel (Ni) on stainless steel solidification modes and microstructures are well detailed. However, fusion-based additive manufacturing (FBAM) processes that rely on faster, more rapid solidification rates call for a more in-depth understanding of the effects of rapid solidification on the solidification behavior and how these change with variations in alloy composition. Eleven custom stainless steel (SS) alloys with unique compositions based on 316L were made using targeted alloying element additions to generate feedstock with Cr/Ni eq ratios ranging from 0.9 to 2.1. Two-piston splat quenching (SQ) was used to produce rapid solidification conditions similar to those achieved in powder bed fusion processes in a fraction of the time and cost of traditional methods. Employing heat transfer simulations, SEM, TEM, and STEM characterization techniques, the SQ process was found to consistently produce solidification rates estimated to be between ~ 0.4 and 1.6 m/s. Five unique solidification microstructures were identified within the rapidly solidified SQ samples (primary austenite, massively transformed ferrite to austenite, primary austenite and massively transformed ferrite to austenite, primary ferrite, and primary ferrite and massively transformed ferrite to austenite). Both the primary ferrite solidification mode and primary ferrite phase were found to form at lower Cr/Ni eq values than previously predicted. SQ samples with compositions that were within the compositional specifications for 316L SS, an alloy which is commonly used in FBAM processes, demonstrated a wider range of potential solidification modes and microstructures that could form at rapid solidification rates than expected. Finally, by using SQ as a means to simulate rapid solidification conditions similar to those observed in FBAM processes, potential new alloy compositions can be screened faster and more cost effectively than purchasing and running test batches of the metallic powder.
Integrating genome-scale metabolic networks with reactive transport models (RTMs) provides a detailed description of the dynamic changes in microbial growth and metabolism. Despite promising demonstrations in the past, computational inefficiency has been pointed out as a critical issue to overcome because it requires repeated application of linear programming (LP) to obtain flux balance analysis (FBA) solutions in every time step and spatial grid. To address this challenge, we propose a new simulation method where we train and validate artificial neural networks (ANNs) using randomly sampled FBA solutions and incorporate the resulting surrogate FBA model (represented as algebraic equations) into RTMs as source/sink terms. We demonstrate the efficiency of our method via a case study of Shewanella oneidensis MR-1. During aerobic growth on lactate, S. oneidensis produces metabolic byproducts (such as pyruvate and acetate), which are subsequently consumed as alternative carbon sources when the preferred nutrients are depleted. To effectively simulate these complex dynamics, we used a cybernetic approach that models metabolic switches as the outcome of dynamic competition among multiple growth options. In both zero-dimensional batch and one-dimensional column configurations, the ANN-based surrogate models achieved substantial reduction of computational time by several orders of magnitude compared to the original LP-based FBA models. Moreover, the ANN models produced robust solutions without any special measures to prevent numerical instability. These developments significantly promote our ability to utilize genome-scale networks in complex, multi-physics, and multi-dimensional ecosystem modeling.
Not Available
Abstract Flowing electrode capacitive deionization systems (FE-CDI) have recently garnered attention because of their ability to prevent cross contamination and operate in uninterrupted cycles ad infinitum. Typically, FE-CDI electrodes suffer from low conductivity, reducing deionization performance. Utilization of higher mass loadings to combat this leads to poor rheological properties. Herein, Ti 3 C 2 T x MXene was introduced as 1 mg mL −1 slurry electrodes in an FE-CDI system for the removal and recovery of ammonia from simulated agricultural wastewater. The electrode performance was evaluated by operating the FE-CDI system with a feed solution of 500 mg L −1 NH 4 Cl running in batch mode at a constant voltage of 1.20 and −1.20 V in charging and discharging modes, respectively. Despite the low loading, Ti 3 C 2 T x flowing electrodes showed markedly improved performance, achieving 60% ion removal efficiency in a saturation time of 115 min with an adsorption capacity of 460 mg g −1 . To understand the high adsorption performance of the electrodes, physiochemical and structural analysis was done via a variety of characterization techniques such as SEM, TEM, XRD, DLS, and Raman spectroscopy. Cyclic voltammetry and galvanostatic charge/discharge profiles were obtained to evaluate the electrochemical properties of the electrodes. The system proved to be an energy-saving technology by exhibiting a charge efficiency of 58–70% while operating at an energy consumption of 0.45 kWh kg −1 . A 92% regeneration efficiency showed that the electrodes were stable and suitable for long term and scalable usage. The results demonstrate that MXenes have the potential to improve the FE-CDI process for energy-efficient removal and recovery of ammonia.
This report summarizes the data collected during the batching and melting of the Direct Feed High-Level Waste APPS Model Glass Matrix (DFHLW APPS) to serve as a quality-assured validation of the Aspen Process Performance Simulation (APPS) formulation method. Of 15 glasses tested, 12 satisfied all target property constraints. Two glasses, APPS-05 and -06, formed nepheline on canister centerline cooling heat-treatment and failed the Product Consistency Test response limits. Glass APPS-07-2 formed unacceptably high concentrations of crystals (primarily Na3Nd(PO4)2) when heat treated at 950 °C. All other glasses were found to be satisfactory. The measured property values were compared to predicted values from a set of current models. In many cases the current models were found to be inadequate for design of DFHLW glasses. These models are being adjusted to correct for mispredictions. Other models, e.g., density, toxicity characteristic leaching procedure, and sulfur solubility, are adequate for formulation of DFHLW glasses.