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.
Explore the source record for details and available documents.
Here, this work optimizes micro-prismatic high-temperature gas reactor (HTGR) designs to reduce the energy-normalized mass of spent nuclear fuel (SNF) and high-level waste (HLW) produced. The optimization was performed for the current graphite moderator and an inert matrix fuel (IMF) concept employing different composite moderators in a prismatic design architecture. The fuel matrix is magnesium oxide (MgO) with entrained tristructural-isotropic (TRISO) fuel. The moderator materials, including beryllium oxide (MgO-BeO) and beryllium (MgO-Be) at 40 vol % loading and yttrium hydride (MgO-YH x=1.9 ) and zirconium hydride (MgO-ZrH x=1.9 ) at 15 vol % loading, were entrained within the MgO host matrix. A generic graphite micro-prismatic HTGR is used as the baseline point design where the external dimensions are held constant. The composite moderator designs use 19.9% enriched uranium nitride TRISO fuel and hexagonal assemblies. For each IMF concept, an optimization study was performed to maximize the discharge burnup of the fuel by varying the TRISO packing fraction and the lattice pitch of the assemblies. The mass of SNF and HLW, other waste metrics, fuel cost, environmental impact metrics, and the activity of the SNF and HLW at 100 years and 100 000 years were calculated for the optimized IMF and graphite reference designs. The IMF results were subsequently compared to those of the graphite reference and the values for a light water reactor (LWR) and a small modular LWR. For the SNF and HLW, all the IMF concepts and the graphite reference produced less waste compared to the traditional LWR designs. However, the IMF concepts outperformed the graphite reference regarding the mass of SNF and HLW. For the other waste metrics, the IMF concepts showed reductions in fuel cost with improved environmental metrics relative to the graphite reference. Overall, the IMF concepts significantly reduced the SNF and HLW produced per unit of energy generated compared to traditional LWR designs.
The Minos Computing Library (MCL) is a task-based programming language and runtime for extremely heterogeneous systems. MCL facilitates writing program for heterogeneous devices and porting applications across different systems and devices. MCL supports asynchronous execution of computing tasks on all available heterogeneous devices, including GPUs, FGPAs, fixed-point accelerators, and AI accelerators. MCL provides a high-level programming interface and automatically and autonomously performs resource management, load balancing, and locality-aware scheduling.
matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.
The goal of the Scalable Thin Light-Emitting Diode (LED) Light Sheet Platform project is to derisk the manufacturing scalability for a high chip count heterogeneous smart lighting platform, with the potential to enable significant energy savings through dynamic directional light control and custom. Specifically, the project derisks the scalability of a disruptive new computer controlled microassembly fabrication process SRI International is developing, to address fundamental cost barriers to mass production of high chip count systems. This is key for enabling mass adoption and thus maximum societal energy savings impact. The LED light sheet technology is a smart illumination platform that can be applied to lighting, signage, and display. It has features and a form factor similar to bendable or conformal OLED light sheets but uses more efficient LEDs with a remote phosphor layer that is very close to the LED to deliver a luminous efficacy that exceeds 125 lumen/Watt and enable > 50% Lighting Application Efficiency (LAE) energy savings. During Budget Period 1 (BP1), all of the BP1 milestones were successfully completed: (M2.0.1) Design demonstrator details that will meet the final project goals, (M6.1.1) Show lighting output model can predict illuminance, spatial, spectral distribution for specific light-sheet designs and use case, (M3.2.1) Automated loading supports 200 chip arrays, and (M3.3.1) Finish first interconnect process run for 200 chip array. The BP1 Go/No-Go Decision Point (G/NG 1) 200 chip demonstrator sample was in the process of final assembly and test at the end of BP1 on December 31, 2023. The objective was successfully achieved on January 16, 2024, when a wired 200-chip sample underwent confirmation testing that demonstrated a 93.75% first pass electrical yield that successfully exceeded the 75% first-pass electrical yield BP1 G/NG requirement. In addition, all 93.75% of the sample devices lit up, further demonstrating the ability to assemble and transfer small LED chips without damaging their functionality.
It is of imperative interests for regional transmission organizations (RTOs) to effectively extract daily load profiles at transmission buses, which remains a gap in existing technology paradigm. This digest proposes an explicit yet efficient linear estimator, to disaggregate metered load profiles at buses with significant behind-the-meter (BTM) solar generations in a data driven manner. The proposed estimator is based on utility zonal load profiles and proxy solar irradiance profiles, which in reality is the aggregated waveform at each transmission bus and equivalent to the mix of summed load profiles minus actual BTM solar generation. To overcome technical challenges in the lack of “ground truth” and validate the performance of supervised learning algorithms, we propose semi-supervised mechanisms with parameter tuning, and leverage the unique characteristics of zero-crossing points in BTM solar peaking behaviors.
In Wendelstein 7-X (W7-X), stable plasma detachment in ECRH hydrogen plasmas is routinely achieved with the standard magnetic configuration, containing five magnetic islands in the scrape-off layer (SOL). Plasma detachment induced by intrinsic carbon impurities is characterized by a high radiation fraction (ƒ rad ) with significantly reduced divertor heat load and particle flux. In this work, we present the radiation dynamics during the detached plasma (DP) phase, such as the shift of radiation zones in radial and poloidal directions towards the X-points as ƒ rad increases (∼0.6–0.9), as well as the penetration of the radiation layer into the confinement region in the deep DP phase with ƒ rad approaching unity. In particular, a structure of multi-X-point radiation (multi-XPR) with an up/down asymmetry in the DP phase is highlighted, which is revealed by bolometer tomography and is further confirmed by video diagnostics. The multi-XPR structure forms helical 3D bands aligned with W7-X field periodicity. Field reversal experiments show that the brightest XPR displaces between the upper and lower SOL regions, suggesting the E× B drift effects. This paper presents the multi-XPR structure observed in the W7-X plasma in detail for the first time. A simplified model considering the influence of the poloidal E× B drift ( V d ) on the impurity flow in the SOL shows that the downstream drift toward the target or target shadow region ( V d < 0) decreases the impurity content, while the upstream drift toward the LCFS ( V d > 0) increases the impurity content. On this basis, the poloidal drift potentially leads to an up/down asymmetry of impurity density in the SOL despite the symmetry magnetic topology. The dynamics of the up/down asymmetry in the multi-XPR structure is also related to the magnitude V d /D (normalized to the impurity diffusivity), with an additional effect owing to the radial inward shift of the emission zone. These results provide new insights into impurity-induced detachment dynamics and provide a basis for improving 3D modeling of impurity transport, which typically does not consider drifts.
Here, this study introduces an airfoil-based refinement technique to enhance the Actuator Disk Model (ADM) for improved wind turbine aerodynamic load prediction and structural simulation in conjunction with Large Eddy Simulations of the wind flow. While ADM offers higher computational efficiency than the more detailed but resource-intensive Actuator Line Model (ALM), it traditionally lacks the resolution needed to capture the localized blade forces accurately. To address this limitation, we introduce a refinement technique that uses airfoil-specific data and employs interpolation-based grid point refinement, achieving ALM-comparable accuracy while preserving ADM's efficiency. Unlike conventional ADM that provides only rotor-disk averaged forces, our synthetic method tracks transient aerodynamic load variations over multiple blade revolutions, allowing us to calculate the distributions of maximum and minimum loads during typical cycles. Applied to the NREL 5 MW reference turbine, our enhanced ADM accurately predicts key aerodynamic parameters (angle of attack, axial velocity, lift, drag, axial and tangential forces along the blades) as well as structural responses (blade tip deflection, maximum stress, and stress concentration). Our results show that the tip deflection ranges from 2.33m (3.69 % of blade length) to 4.28m (6.79 %), with maximum stress concentration occurring near the blade root. This research demonstrates that a refined synthetic ADM approach can serve as a computationally efficient alternative for both aerodynamic analysis and structural simulation of wind turbine blades subjected to realistic wind fields.
Wind turbine wake models are essential tools for predicting power losses and structural loads in wind farms. Among these, the dynamic wake meandering (DWM) model, included as a recommended approach in the International Electrotechnical Commission design standard, is a widely used engineering-fidelity method that balances accuracy and computational cost. This study compares the performance of three DWM-based wake model implementations (from the Technical University of Denmark, the National Renewable Energy Laboratory, and the Institute for Energy Technology) under below-rated wind speed conditions. Model predictions of wake flow, power output, and structural loads for a four-turbine row are evaluated across different ambient turbulence levels and wind-direction misalignments and compared against high-fidelity large-eddy simulation results. All three models captured the overall wake evolution and mean turbine performance with reasonable accuracy; their predicted time-averaged thrust and power were typically within 5 %–10 % of the large-eddy simulation benchmark. However, notable differences emerged in wake structure and unsteady load predictions, with discrepancies increasing for turbines further downstream. These differences highlight the importance of modelling choices such as wake summation and turbulence treatment, which strongly influence power-deficit and fatigue-load predictions. Comparison with large-eddy simulations reveals each approach's strengths and weaknesses, indicating where improvements are needed. Overall, the findings point to specific refinements for DWM models to improve their fidelity, ultimately enabling more robust wake predictions for wind farm design and operation.
The technical and economic performance as well as the load-following capabilities of grid-connected geothermal hybrid systems were assessed in this work. The analyzed geothermal hybrid configuration is composed of a binary geothermal plant integrated with a concentrating solar-thermal system and underground thermal energy storage (UTES) through a primary heat exchanger. Physics-based models for the hybrid system for plant generation capacities of 1, 25, and 50 MW were developed from validated models for each subsystem. Also, an economic model was developed that accounts for different hybrid system capabilities, solar field sizes, and thermal storage duration. The advantage of the geothermal hybrid system was assessed by comparing the performance with the baseline benchmark geothermal plant with a similar configuration and generation capacity. It was found that hybridizing geothermal plants with concentrating solar and thermal energy storage not only improves the thermal efficiency by up to 8 percentage points when additional heat from the solar-UTES loop rises the evaporator temperatures from 70 to 125 °C, but also enhances the load-following capability for the geothermal plant, which can meet a typical residential load profile with a power rate of change 0.25 kW/s with an absolute error under 13 kW for a 1 MW plant. Other benefits of hybridization include resource preservation and a potential LCOE reduction of up to 56% for a 50 MW geothermal hybrid plant having a 50% solar share, a 1.4 solar multiple, and 24-h storage capacity. The results presented in this work demonstrate that hybridizing geothermal systems transforms them into a flexible and cost-effective solution for addressing the dynamic requirements of modern electric grids.
Toughening polymers has attracted significant interest. Traditionally, polymer toughness is enhanced by constructing polymer networks or introducing sacrificial bonds into the chains between crosslink points. These strategies, though, introduce pronounced energy dissipation and associated heat, both of which are undesirable under long-term cyclic loading, for example at the interface of implants in the human body. By incorporating single-chain nanoparticles (SCNPs) into linear polymer chains to generate all-polymer nanocomposites (APNCs), we have been able to achieve high strength, high toughness with low energy dissipation. Using a combination of simulation and experimental results, we are advancing a “SCNPs effect” where tightly cross-linked SCNPs produce a modulus contrast to achieve strengthening and toughening. Benefitting from the soft interface, the penetrable and deformable SCNPs cause the surrounding polymer chains to move in concert, significantly reducing the interfacial friction to achieve low energy dissipation. The intramolecular cross-linking of the SCNPs and adhesion between the SCNPs and polymer matrix are critical for realizing such high-performance systems. Based on a Gaussian regression model and back propagation (BP) neural network, the mechanical strength can be predicted and is supported by simulations. Finally, the APNC concept described can be applied to elastomers and gels, broadening its utilization in high-cycle and low-dissipation applications, like soft robots, flexible sensors and cartilage replacements, and artificial heart valves.
We present a novel distributed union-find algorithm that features asynchronous parallelism and k-d tree based load balancing for scalable visualization and analysis of scientific data. Applications of union-find include level set extraction and critical point tracking, but distributed union-find can suffer from high synchronization costs and imbalanced workloads across parallel processes. In this study, we prove that global synchronizations in existing distributed union-find can be eliminated without changing final results, allowing overlapped communications and computations for scalable processing. We also use a k-d tree decomposition to redistribute inputs, in order to improve workload balancing. We benchmark the scalability of our algorithm with up to 1,024 processes using both synthetic and application data. Here, we demonstrate the use of our algorithm in critical point tracking and super-level set extraction with high-speed imaging experiments and fusion plasma simulations, respectively.
In this paper, we propose a low-cost bidirectional four-switch inverter that transfers three-phase energy between a permanent magnet synchronous generator (PMSG) driven by a wave energy collector and the power grid or a local ac load. To achieve high control performance, the machine-side inverter uses a high-gain observer for rotor position estimation and a Lyapunov-based approach for torque control. For the gridside inverter control, we estimate the filter capacitor voltage with a model reference adaptive system and control the dc-link voltage with a Lyapunov-based energy function. We also track the maximum power point from the sea waves with an integral control law that actively matches the generator impedance with that of the wave collection device. The proposed system supports both grid-following/forming modes. We validate our control design with simulation results on a 3-kW hardware setup.
A long standing challenge in computational materials science is to establish a quantitative connection between the macroscopic properties of plastic deformation with the microscopic mechanisms of dislocations in crystalline materials. Although the discrete dislocation dynamics (DDD) simulation method has been developed for several decades with the goal of addressing this challenge, a one-to-one comparison between the DDD predictions on single crystal stress–strain curves and experimental measurements under identical conditions has not been possible to date. Such a comparison is an essential step towards establishing a dislocation-physics based theory of plasticity and a multiscale framework of the plastic behaviors of crystalline materials. Here we provide direct comparisons between the stress–strain curves of Cu single crystals under high strain rate loading in the [0 0 1] and [0 1 1] directions obtained from miniaturized desktop Kolsky bar experiments and those from DDD simulations under identical loading conditions. With an appropriate set of parameters, DDD simulations can produce stress–strain curves that are in reasonable agreement with the experimental results. However, the dislocation mobility values needed to achieve this agreement are an order of magnitude lower than expected based on previous measurements and atomistic simulations. We hypothesize that this discrepancy could be caused by drag forces from jogs and point defects produced during the plastic deformation. Cross-slip of screw dislocations is also found to be necessary to capture the experimental stress–strain behavior, especially for the [0 1 1] loading direction. Finally, this work provides an example of how direct comparisons between DDD simulations and experimental measurements can provide new insight into the fundamental mechanisms of plastic deformation.
This paper studies a power transmission system with both conventional generators (CGs) and distributed energy assets (DEAs) providing frequency control. We consider an operating condition with demand aggregating two dynamic components: one that switches between different values on a finite set, and one that varies smoothly over time. Such dynamic operating conditions may result from protection scheme activations, external cyber-attacks, or due to the integration of dynamic loads, such as data centers. Mathematically, the dynamics of the resulting system are captured by a system that switches between a finite number of vector fields - or modes -, with each mode having a distinct equilibrium point induced by the demand aggregation. To analyze the stability properties of the resulting switching system, we leverage tools from hybrid dynamic inclusions and the concept of ..omega.. -limit sets from sets. Specifically, we characterize a compact set that is semi-globally practically asymptotically stable under the assumption that the switching frequency and load variation rate are sufficiently slow. For arbitrarily fast variations of the load, we use a level-set argument with multiple Lyapunov functions to establish input-to-state stability of a larger set and with respect to the rate of change of the loads. The theoretical results are illustrated via numerical simulations on the IEEE 39-bus test system.
This work details the synthesis, materials characterization and absorption capacity of graphene aerogel (GA) loaded with UiO-66-NH 2 metal organic framework (MOF). Throughout the study, a series of density-tunable MOF/GA composites were synthesized by adjusting the MOF/graphene oxide (GO) mass ratio of the precursors before hydrothermal reduction and subsequent critical point drying to achieve MOF-x%/GA composites, with x, ranging from 0 to 100, denoting the weight percent ratio of MOF to GO in the starting solution. Scanning electron microscopy (SEM) images and Barrett-Joyner-Halenda (BJH) pore size calculations supported that MOF loading enlarged the macroporous (>50 nm) structure of the GA framework but had no influence on the mesoporous (2–50 nm) structure of GA. The degree of MOF loading first decreased and then increased the density of resulting composites, exhibiting a minimum for the 10–30% samples. Such behavior with respect to MOF loading has not been reported previously. MOF-x%/GA samples were further investigated for absorption capacity using various organic liquids. MOF-30%/GA proved to be the best absorbent sample for all solvents tested, achieving the highest capacity for chloroform at 147.0 ± 10.0 mg/mg. The behavior is attributed to the structural changes induced by the MOF incorporation as well as the interactions between the organic molecules and the MOF.
Machine learning methods have shown initial promise in constitutive modeling for single crystals or homogenized polycrystals, delivering notable computational efficiency. However, existing machine learning-based constitutive models often lack generalizability, limiting their application across diverse boundary value problems. This study introduces a thermodynamics-informed artificial neural network model to accelerate rate-tangent crystal plasticity fast Fourier transform simulations for cross-scale deformation behaviors of polycrystals under complex loading. Our model integrates microstructural variability and local interactions effectively. To address local effects in each grain, we employ K-means clustering to group Gauss points within the microstructure into clusters assumed to be in similar mechanical states. This approach, based on self-clustering analysis, extends model scope from macroscopic stress response to the granular level, capturing mechanical responses and orientation evolution across grains. This reduces the number of nonlinear problems to solve, with cluster responses propagated throughout each group. The thermodynamics-based artificial neural network-extracted features are further processed using local material state clusters to account for history-dependent deformation and evolving microstructures. Additionally, representative volume element simulations with rate-tangent crystal plasticity fast Fourier transform provide reliable datasets for model training. The proposed model demonstrates high efficiency, accuracy, self-consistency, and enhanced generalizability in predicting strain–stress responses and orientation evolution at both individual grain and aggregate scales under complex loading conditions, such as biaxial tension and arbitrary loading scenarios.
The accurate simulation of complex biochemical phenomena has historically been hampered by the computational requirements of high-fidelity molecular-modeling techniques. Quantum mechanical methods, such as ab initio wave-function (WF) theory, deliver the desired accuracy, but have impractical scaling for modeling biosystems with thousands of atoms. Combining molecular fragmentation with MP2 perturbation theory, this study presents an innovative approach that enables biomolecular-scale ab initio molecular dynamics (AIMD) simulations at WF theory level. Leveraging the resolution-of-the-identity approximation for Hartree-Fock and MP2 gradients, our approach eliminates computationally intensive four-center integrals and their gradients, while achieving near-peak performance on modern GPU architectures. The introduction of asynchronous time steps minimizes time step latency, overlapping computational phases and effectively mitigating load imbalances. Utilizing up to 9,400 nodes of Frontier and achieving 59% (1006.7 PFLOP/s) of its double-precision floating-point peak, our method enables us to break the million-electron and 1EFLOP/s barriers for AIMD simulations with quantum accuracy.