Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational models”

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.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

DIRAC current, upcoming and planned capabilities and technologies

DIRAC is the interware for building and operating large scale distributed computing systems. It is adopted by multiple collaborations from various scientific domains for implementing their computing models. DIRAC provides a framework and a rich set of ready-to-use services for Workload, Data and Production Management tasks of small, medium and large scientific communities having different computing requirements. The base functionality can be easily extended by custom components supporting community specific workflows. DIRAC is at the same time an aging project, and a new DiracX project is taking shape for replacing DIRAC in the long term. This contribution will highlight DIRAC’s current, upcoming and planned capabilities and technologies, and how the transition to DiracX will take place. Examples include, but are not limited to, adoption of security tokens and interactions with Identity Provider services, integration of Clouds and High Performance Computers, interface with Rucio, improved monitoring and deployment procedures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Directional response of several geometries for reactor-neutrino detectors

Here we have modeled six abstracted detector designs, with the goal of determining their ability to resolve the direction to an antineutrino source, including two for which we have operational data for validating our computer modeling and analytical processes. We have found that the most promising options, regardless of scale and range, have angular resolutions on the order of a few degrees, which is better than any yet achieved in practice by a factor of at least two. We examine and compare several approaches to detector geometry for their ability not only to detect the inverse beta decay (IBD) reaction, but also to determine the source direction of incident antineutrinos. The information from these detectors provides insight into reactor power and burning profile, which is especially useful in constraining the clandestine production of weapons material. In a live deployment, a nonproliferation detector must be able to isolate the subject reactor, possibly from a field of much-larger power reactors; directional sensitivity can help greatly with this task. We also discuss implications for using such detectors in longer-distance observation of reactors, from a few kilometers to hundreds of kilometers.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Cooperative Education

Los Alamos National Laboratory (LANL) is a multidisciplinary national laboratory that conducts research and development in national security, engineering, materials science, computational modeling, and advanced manufacturing. The laboratory develops innovative technologies to address complex scientific and engineering challenges. This project focuses on the development and evaluation of high-performance absorbing structures through computational design, simulation, and engineering analysis. Absorbing structures are used in applications where damage mitigation, structural protection, and material efficiency are critical performance requirements. The increasing demand for lightweight, high-strength, and highly efficient structural systems has created a need for improved design methodologies capable of maximizing absorption while minimizing weight and material usage. The project utilizes advanced engineering software, including 3D CAD software and FEA, to generate and optimize structural concepts. Computational simulations are performed to evaluate structural behavior under loading conditions, while mathematical analyses are conducted using Python-based tools as well as established analytical equations from material and structural mechanics. The project benefits LANL by supporting the development of advanced design methodologies and improving the understanding of material and structural performance. During the internship term, a significant portion of the design development, simulation, and data analysis activities will be completed. Success of the project depends on collaboration among engineering mentors and technical staff members. Work will be conducted at Los Alamos National Laboratory using laboratory computing resources and engineering software.

42 ENGINEERING↗
DOE OSTI2025

Investigating Electron Conductivity Regimes in the Bacterial Cytochrome Wire OmcS

The anaerobic bacterium Geobacter sulfurreducens produces extracellular, electronically conductive cytochrome polymer wires that are conductive over micron length scales. Structure models from cryo-electron microscopy data show OmcS wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. However, the mechanism by which this heme arrangement supports long-range electronic conduction remains unknown. Structure models from cryo-electron microscopy data show these wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. Existing computational models using static heme redox potentials and coupling energies fail to explain experimental observations, predicting conductances 10,000 to 100,000 times lower than measured values. Here, we investigate how dynamic disorder affects site energies, interheme coupling, and long-range electronic conductivity within these cytochrome wires. We introduce an approach to extract charge carrier site information directly from Kohn–Sham density functional theory, without employing projector schemes, and show that site and coupling energies are highly sensitive to changes in interheme geometry and the surrounding electrostatic environment. Unlike models that incorporate dynamic disorder as a thermally averaged quantity, our quantum charge carrier model incorporates proxies for dynamic disorder through decoherence corrections, yielding predicted diffusion coefficient closer to what is expected from experiment and comparable with other organic-based electronic materials. Based on these simulations, we propose that the instantaneous fluctuations of the local electrostatic environment can transiently lift energy degeneracies and delocalize charge carriers. Furthermore, these studies reveal how incorporating dynamic fluctuations associated with the environment resolves the discrepancy between theory and experiment in microbial cytochrome wires and highlight design principles for bioinspired, heme-based conductive materials.

Bioinorganic chemistry↗
DOE OSTI2024

A Eulerian three-fluid flow framework for the study of fuel dispersal behavior under loss-of-coolant accident conditions

As the nuclear industry explores new fuel designs to accommodate increased burnup, studying fuel behavior during loss-of-coolant accidents is essential for ensuring the safe operation of light-water reactors. Here, this study presents a computational model to simulate the intricate three-phase flow dynamics of fuel dispersal following cladding breach, where solid fuel particles and carrier fission gases interact with the surrounding liquid or gas medium. The model utilizes a Eulerian three-fluid framework to depict the phase interaction within interpenetrating continua, treating fragmented fuel as granular material. Closure relationships for bulk-flow properties of the solid phase and detailed interfacial exchange terms in high-temperature, high-speed gas-solid-liquid flow scenarios are proposed. The model was validated using adiabatic experimental results in the literature, focusing on characterizing fuel dispersal behavior during loss-of-coolant accidents conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched probability density functions reported in experimental studies. Key highlights of this study include the theoretical description of a coupled three-phase flow with phase change, interpretation of the high-pressure boundary conditions and insights into the transient behavior of gas-solid-liquid phase dynamics during a simulated high-speed dispersal event.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Quantum utility in simulating the real-time dynamics of the Fermi–Hubbard model using superconducting quantum computers

The Fermi–Hubbard model is a fundamental model in condensed matter physics that describes strongly correlated electrons. On the other hand, quantum computers are emerging as powerful tools for exploring the complex dynamics of these quantum many-body systems. In this work, we demonstrate the quantum simulation of the one-dimensional Fermi–Hubbard model using IBM's superconducting quantum computers, employing over 100 qubits. We introduce a first-order Trotterization scheme and extend it to an optimized second-order Trotterization for the time evolution in the Fermi–Hubbard model, specifically tailored for the limited qubit connectivity of quantum architectures, such as IBM's platforms. Notably, both Trotterization approaches are scalable and maintain a constant circuit depth at each Trotter step, regardless of the qubit count, enabling us to precisely investigate the relaxation dynamics in the Fermi–Hubbard model by measuring the expectation value of the Néel observable (staggered magnetization) for time-evolved quantum states. Lastly, our successful measurement of expectation values in such large-scale quantum many-body systems, especially at longer time scales with larger entanglement, highlights the quantum utility of superconducting quantum platforms over conventional classical approximation methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

New Particle Formation and Growth in the Houston Atmosphere During TRACER (Final Report)

From 2020-2025, researchers from UC Irvine, UC Riverside, and Colorado State University collaborated on a Department of Energy-funded project to understand how airborne particles form and grow in urban atmospheres, conducting an intensive field campaign in Houston, Texas during summer 2022. Using advanced instruments to measure gas-phase chemicals, particle composition, and a specialized chamber to study particle growth, the team discovered that sulfur-containing compounds from industrial and power plant emissions are the dominant driver of new particle formation in Houston, with particles typically forming locally in the city and growing as air moves away in the urban plume. The research revealed an important methodological insight: measurements from fixed ground stations can be misleading when interpreting how particles actually evolve as air masses move, which has significant implications for how scientists worldwide interpret atmospheric observations. These findings improve understanding of urban air quality and help reduce uncertainties in climate models, since these particles play critical roles in cloud formation and Earth's radiation balance, while also providing detailed information about ultrafine particle composition relevant to public health. The project trained three doctoral students, developed enhanced computer models for urban particle formation, and made all data publicly available through the DOE Atmospheric Radiation Measurement data archive for use by the broader scientific community.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hierarchical Bayesian modeling for Inverse Uncertainty Quantification of system thermal-hydraulics code using critical flow experimental data

The best estimate plus uncertainty methodology in nuclear system thermal-hydraulic studies necessitates a comprehensive understanding of uncertainties in system code predictions. The forward uncertainty quantification (UQ) process involves the propagation of input uncertainties through the computational models to obtain uncertainties in the outputs. To this end, achieving an accurate estimation of input uncertainties is important, which is the focus of inverse UQ (IUQ). Traditionally, research in Bayesian IUQ within the nuclear engineering domain has largely relied on single-level Bayesian inference. While being effective for relatively small datasets, this approach encounters limitations for cases with large datasets. The use of a single-level model may prove inefficient, as the resultant posterior distributions can significantly differ when distinct subsets of data are employed. To address this issue, we employ an hierarchical Bayesian model for IUQ. Furthermore, this approach involves organizing observations into different groups based on the test conditions, thereby accommodating varying calibration parameters across these distinct groups. In this study, we developed and implemented a hierarchical Bayesian IUQ method to consider the grouping effect of critical flow measurement data from various geometries. Comparing the outcomes of IUQ under different selections of test data using hierarchical Bayesian IUQ against those obtained from single-level Bayesian IUQ, the forward propagation of hierarchical Bayesian IUQ results demonstrates a notably improved agreement with the experimental data.

42 ENGINEERING↗
DOE OSTI2026

Uncertainty-Driven Rapid Thermodynamic Assessment of Nb-Ta-Zr System and Effects of C impurities (L25GF9298S): Annual Progress Report

An integrated computational materials engineering (ICME) method is in development for rapid thermodynamic experimental investigation and high-fidelity computational modeling of refractory multi-principal element alloys (RMPEAs). These ultra-high-temperature (UHT) alloys are of interest for structural applications in extreme environments, but deficiency of reliable data, especially melting temperatures, impedes the prediction of alloys with favorable properties. The method leverages UHT capabilities and computational expertise of LLNL’s Materials Science Division and the McCormack Lab’s UHT conical nozzle levitation (CNL) system to iteratively map the Nb-Ta-Zr phase space, with focus on the liquidus surface, through targeted experiments selected by quantifying uncertainty in the thermodynamic model fitting parameters. This method will reduce the time to map uncharted RMPEA phase space and thereby accelerate discovery and development of advanced materials for applications in extreme environments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Accelerating kinetic plasma simulations with machine-learning-generated initial conditions

Computational models of plasma technologies often solve for the system operating conditions by time-stepping an initial value problem to a quasi-steady solution. However, the strongly nonlinear and multi-timescale nature of plasma dynamics often necessitate millions, or even hundreds of millions, of steps to reach convergence, reducing the effectiveness of these simulations for computer-aided engineering. We consider acceleration of kinetic plasma simulations via data-driven machine-learning-generated initial conditions, which initialize the simulations close to their final quasi-steady-state, thereby reducing the number of steps to reach convergence. Three machine-learning models are developed to predict the density and ion kinetic profiles of capacitively coupled plasma discharges relevant to the microelectronics industry. The models are trained on kinetic simulations over a range of device operating frequencies and pressures. Best performance was observed when simulations were initialized with ion kinetic profiles generated by a convolutional neural network, reducing the mean number of steps to reach convergence by 17.1× when compared to initialization with a zero-dimensional global model. We also outline a workflow for continuous data-driven model improvement and simulation speedup, with the aim of generating sufficient data for full device digital twins.

Artificial neural networks↗
DOE OSTI2026

Continuity of reaction kinetics across the pressure and materials gaps in CO oxidation on FeO–Pt interfaces

Translating atomic-scale insights from surface science studies of model catalysts to practical powder catalysts remains a persistent challenge in heterogeneous catalysis. Here, in this study, we demonstrate mechanistic continuity across the pressure and materials gaps during CO oxidation at the FeO-Pt interface using in situ microscopy, spectroscopy and computational modelling. Under reaction conditions, coordinatively unsaturated Fe (Fe cus ) sites at the interface enable selective O 2 activation on CO-saturated surfaces, circumventing the CO-poisoning limitation of platinum-group metals. We identify parallel reaction pathways involving the *O 2 -*CO intermediate. Remarkably, activation energies remain consistent at 12-15 kJ mol −1 (0.12-0.16 eV) from ultrahigh vacuum to atmospheric pressures and from FeO/Pt(111) model catalysts to FeO/Pt powder catalysts, validating mechanistic insights derived from surface science studies. Our findings show an example of bridging the long-standing divide between model and practical catalyst systems, establishing an effective approach to capture catalytic behaviours under operational conditions and advancing mechanism-driven catalyst design.

03 NATURAL GAS↗
DOE OSTI2024

Modeling-Based Design and Optimization of a Gradient Composite Transition Joint

An innovative additively manufactured gradient composite transition joint (AM-GCTJ) has been designed to join dissimilar metals, to address the pressing issue of premature failure observed in conventional dissimilar metal welds (DMWs) when subjected to increased cyclic operating conditions of fossil fuel power plants. The transition design, guided by computational modeling, developed a gradient composite material distribution, facilitating a smooth transition in material volume fraction and physical properties between different alloys. This innovative design seeks to alleviate structural challenges arising from distinct material properties, including high thermal stress and potential cracking issues resulting from the thermal expansion mismatch typically observed in conventional DMWs. In this study, we investigated the creep properties of transition joints comprising Grade 91 steel and 304 stainless steel through a combination of simulations and creep testing experiments. The implementation of a gradient composite design in the plate transition joint resulted in a significant enhancement of creep resistance when compared to the baseline conventional DMW. For instance, the creep rupture life of the transition joint was improved by > 400% in a wide range of temperature and stress testing conditions. Meanwhile, the failure location shifted to the base material of Grade 91 steel. Such enhancement can be primarily attributed to the strong mechanical constraint facilitated by the gradient composite design, which effectively reduced the stresses on the less creep-resistant alloy in the transition zone. Beyond examining plate joints, it is crucial to assess the deformation response of tubular transition joints under pressure loading and transient temperature conditions to substantiate and demonstrate the effectiveness of the design. The simulation results affirm that the tubular transition joint demonstrates superior resistance compared to its counterpart DMW when subjected to multiaxial stresses in tubular structures. In addition, optimization of the transition joint’s geometry dimensions has been conducted to diminish the accumulated deformation and enhance the service life. Lastly, the scalability and potential of the innovative transition joints for large-diameter pipe applications are addressed.

Zhang, Wei↗
DOE OSTI2024

Polynomial chaos expansions on principal geodesic Grassmannian submanifolds for surrogate modeling and uncertainty quantification

In this work we introduce a manifold learning-based surrogate modeling framework for uncertainty quantification in high-dimensional stochastic systems. Our first goal is to perform data mining on the available simulation data to identify a set of low-dimensional (latent) descriptors that efficiently parameterize the response of the high-dimensional computational model. To this end, we employ Principal Geodesic Analysis on the Grassmann manifold of the response to identify a set of disjoint principal geodesic submanifolds, of possibly different dimension, that captures the variation in the data. Since operations on the Grassmann require the data to be concentrated, we propose an adaptive algorithm based on Riemannian K-means and the minimization of the sample Fréchet variance on the Grassmann manifold to identify “local” principal geodesic submanifolds that represent different system behavior across the parameter space. Polynomial chaos expansion is then used to construct a mapping between the random input parameters and the projection of the response on these local principal geodesic submanifolds. Here, the method is demonstrated on four test cases, a toy-example that involves points on a hypersphere, a Lotka-Volterra dynamical system, a continuous-flow stirred-tank chemical reactor system, and a two-dimensional Rayleigh-Bénard convection problem.

42 ENGINEERING↗
DOE OSTI2025

A Generic FSC Wind Park EMT Model with IEEE Std 2800-Compliant Fault Ride-Through Capability

The modern power grid is seeing more and more electricity come from renewable sources like wind farms, which use sophisticated power electronics instead of traditional spinning generators. To keep everything running smoothly and meet industry standards such as IEEE 2800, these systems need smart control strategies. In our work, we built a flexible computer model of a full-scale wind farm converter that can handle grid disturbances without shutting down. When a fault or storm hits, the model’s built-in logic automatically adjusts the currents it sends to the grid and protects its internal energy storage, ensuring the wind farm stays connected and doesn’t damage its own equipment. Once the disturbance clears, the model restores normal operation seamlessly, so there’s no long interruption in power delivery. At the same time, it carefully regulates the voltage where the wind farm ties into the larger grid, helping to maintain safe voltage levels across the network. Our simulations show that this control setup not only meets all the requirements of IEEE Standard 2800 but also allows the wind farm to recover quickly and predictably, keeping the lights on no matter what happens on the grid.

17 WIND ENERGY↗
DOE OSTI2025

Computational Optimization of Room Temperature Usable Capacity for Hydrogen Storage in MFU-4-Type Metal–Organic Frameworks via Pairwise Metal Substitutions

The efficient storage of hydrogen is a critical challenge in the quest for sustainable energy solutions. Current adsorbent-based methods achieve satisfactory storage densities predominantly under cryogenic temperatures and/or high pressures, which imposes problems with cost-efficient and safe implementation of this technology. Materials that can bind hydrogen gas reversibly at ambient temperatures and more moderate pressures could play a pivotal role in enabling hydrogen-powered technologies. In this study, we use reliable computational modeling to investigate two synthetically feasible paths for tuning the enthalpy of H2 binding in MFU-4-type metal–organic frameworks (MOFs), aiming to maximize usable capacity. This study examines MIM4 IICl3(bta)6 (bta– = benzotriazolate) Kuratowski-type clusters as a model for strong binding sites in MFU-4l frameworks. We systematically evaluate the impact of separately tuning the central MII metal ion (which plays a structural role) and the peripheral MI metal ion (which binds the substrate) on the energetics of H2 binding. Our computational study reveals that H2 binding at an MI site mostly follows the trend AgI < CuI < NiI < CoI < AuI while a larger central MII site generally weakens the H2 binding at a MI site. Importantly, we have identified three new combinations of MI and MII to achieve high fractional usable capacities of the total H2 adsorbed under a pressure swing from 5 to 100 bar at room temperature. Additionally, we examine the nature of the binding interaction between the peripheral metal atom and the hydrogen molecule. While charge transfer predominantly induces this interaction, for several atom combinations, a change in the polarization (associated with variations in the ionic radius of the MI binding atom) is another important factor for adjusting the strength of the interaction. We suggest that the proposed compositions of Kuratowski-type clusters are highly desirable synthetic targets for future laboratory study.

Tkachenko, Nikolay V↗
DOE OSTI2025

Basic Physical Processes Involving Dust in Fusion Plasmas

This report presents the main outcomes of our research focused on understanding the behavior and effects of dust particles in fusion plasmas, particularly in the edge regions of tokamaks and stellarators. We conducted computer modeling studies of burst injections of carbon and tungsten dust particles in DIII-D like divertor plasmas. The studies investigated effects of transient influx of the low- and high-Z material plasma contaminants in form of dust on the edge plasma dynamics in a modern mid-size tokamak. We also performed theoretical and computational studies of the forces acting on non-spherical dust grains in magnetized and nonmagnetized plasmas. In addition, we cooperated with General Atomics slag management group on development of DIII-D dust injection experiments to measure trajectories of dust grains in the divertor plasma and compare them with DUSTT code predictions for validation of the dust modeling capabilities. We collaborated with JET experimentalists on evaluation of radiative losses induced by injection of mixed neon-deuterium ice pellets for disruption and run-away electron generation mitigation during thermal and current quench phases. We also cooperated with experimentalists at LHD fusion device (Japan) on modeling support of experiments on injection of boron granules in fusion plasma discharges to assess their dynamics and effectiveness for in situ dynamic wall conditioning.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Increasing aggregate size reduces single-cell organic carbon incorporation by hydrogel-embedded wetland microbes

Abstract Microbial degradation of organic carbon in sediments is impacted by the availability of oxygen and substrates for growth. To better understand how particle size and redox zonation impact microbial organic carbon incorporation, techniques that maintain spatial information are necessary to quantify elemental cycling at the microscale. In this study, we produced hydrogel microspheres of various diameters (100, 250, and 500 μm) and inoculated them with an aerobic heterotrophic bacterium isolated from a freshwater wetland (Flavobacterium sp.), and in a second experiment with a microbial community from an urban lacustrine wetland. The hydrogel-embedded microbial populations were incubated with 13C-labeled substrates to quantify organic carbon incorporation into biomass via nanoSIMS. Additionally, luminescent nanosensors enabled spatially explicit measurements of oxygen concentrations inside the microspheres. The experimental data were then incorporated into a reactive-transport model to project long-term steady-state conditions. Smaller (100 μm) particles exhibited the highest microbial cell-specific growth per volume, but also showed higher absolute activity near the surface compared to the larger particles (250 and 500 μm). The experimental results and computational models demonstrate that organic carbon availability was not high enough to allow steep oxygen gradients and as a result, all particle sizes remained well-oxygenated. Our study provides a foundational framework for future studies investigating spatially dependent microbial activity in aggregates using isotopically labeled substrates to quantify growth.

59 BASIC BIOLOGICAL SCIENCES↗