Search NASASearch

SEARCH · Search NASA

Results for “computer simulations”

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 phase diagram of quantum chromodynamics in one dimension on a quantum computer

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms.

Quantum information

Solvent-Mediated Control of Nanocellulose Dispersion: An Integrated Computational and Experimental Investigation

Fibrillated cellulose derived from forestry feedstocks represents a renewable and high-strength materials platform for circular bioeconomies. However, its practical implementation is hindered by the irreversible aggregation of nanocellulose architectures, including cellulose nanofibers (CNFs). Solvent-based dispersion offers a simple and practical route to prevent CNF aggregation. Here, in this work, we integrate classical and enhanced sampling molecular dynamics (MD) simulations with experimental suspension rheology and atomic force microscopy (AFM) to elucidate how solvent environments tune CNF–CNF interactions and dispersion stability. CNF–CNF contact free energies computed from MD simulations reveal reduced aggregation in acetone/water, γ-valerolactone (GVL)/water, and tetrahydrofuran (THF)/water and pure acetone compared with pure water, reflecting stronger CNF-solvent relative to inter-CNF interactions. Correspondingly, CNF-solvent suspensions in these solvent systems exhibit stronger inter-fibril network structures and enhanced recovery compared to water, indicating improved CNF-solvent affinity. Liquid cell AFM imaging in acetone–water mixtures and in pure acetone further confirm the presence of well-dispersed CNFs. By combining multiscale computation with targeted experiments, this study establishes a rational framework for solvent design to achieve stable nanocellulose dispersions for high-strength biobased materials and efficient bioenergy conversion.

cellulose

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE

Probabilistic Resource Adequacy Suite (PRAS) v0.8 Model Documentation

The Probabilistic Resource Adequacy Suite, or PRAS, is a software package for studying power system resource adequacy. It allows the user to simulate power system operations under a wide range of operating conditions, in order to study the system's risk of failing to meet demand due to a resource shortfall, and identify the time periods and regions in which that risk occurs. This reports documents version 0.8 of the tool.

29 ENERGY PLANNING, POLICY, AND ECONOMY

A Transient Hydrodynamic Model of Screen Channel Liquid Acquisition Devices for In-Space Cryogenic Propellant Management

Screen channel liquid acquisition devices (LADs) will play a crucial role in future deep space travel. It is essential that vapor-free delivery of propellants during tank-to-tank transfer is ensured to maximize yield from storage tanks and prevent potential combustion instabilities. The screen channel LAD utilizes a fine screen wire mesh that can separate phases in a low Bond number (i.e. microgravity) environment using surface tension forces. This study presents the development and verification of a new model for transient screen compliance, one of the influential factors for screen channel LAD design. Screen compliance is crucial during LAD channel outflow transients because the slight deflection of the screen can provide needed mass to satisfy rapid outflow demands and reduce the pressure difference across the screen. The model is successfully verified against CFD simulations. In addition, the characteristic speed for the governing screen compliance equations is derived which allows for numerical stability criteria to be established. As shown in this study, the transient maximum pressure difference across the screen can greatly exceed the steady state maximum pressure difference across the screen in many cases.

Hydrodynamics Simulations

Modeling Flow in a Microreactor Core: From High-Fidelity CFD to Subchannel Analysis

The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN

Computational Fluid Dynamics at the Lewis Research Center: An Overview

Lewis is a multidisciplinary Center with strong research and development programs in aeronautical and space propulsion, power, space communications, space experiments and materials. Computational fluid dynamics (CFD) is playing an important and growing role in most of these areas. Described here is how CFD is integrated into these programs and highlights elements of the CFD activities. Examples are presented of codes developed to predict flow fields in advanced propulsion systems and several of the code validation experiments are described. The CFD effort at Lewis ranges from basic research on new and improved algorithms through code development to the application of these codes to specific engineering problems. Because of the substantial improvement in CFD's predictive capability, its use at Lewis is on a steep growth path, spreading rapidly into new areas which had not traditionally taken advantage of the techniques of numerical simulation. Multidisciplinary codes and the future direction of CFD at Lewis are discussed.

Robert M Stubbs

Quantum Circuits for the Preparation of Spin Eigenfunctions on Quantum Computers

The application of quantum algorithms to the study of many-particle quantum systems requires the ability to prepare wave functions that are relevant in the behavior of the system under study. Hamiltonian symmetries are important instruments used to classify relevant many-particle wave functions and to improve the efficiency of numerical simulations. In this work, quantum circuits for the exact and approximate preparation of total spin eigenfunctions on quantum computers are presented. Two different strategies are discussed and compared: exact recursive construction of total spin eigenfunctions based on the addition theorem of angular momentum, and heuristic approximation of total spin eigenfunctions based on the variational optimization of a suitable cost function. The construction of these quantum circuits is illustrated in detail, and the preparation of total spin eigenfunctions is demonstrated on IBM quantum devices, focusing on three- and five-spin systems on graphs with triangle connectivity.

97 MATHEMATICS AND COMPUTING

Predicting Atomistic Transitions with Transformers

Accurate knowledge of the atomistic transition pathways in materials and material surfaces is crucial for many material science problems. However, conventional simulation techniques used to find these transitions are extremely computationally intensive. Even with large-scale, accelerated material simulations, the computational cost constrains the applicable domain in practice. Machine learning models, with the potential to learn the complex emergent behaviors governing atomistic transitions as a fast surrogate model, have great promise to predict transitions with a vastly reduced computational cost. Here, we demonstrate how transformers can be trained to predict atomistic transitions in nano-clusters. We show how we evaluate physical validity of the predictions and how a multitude of additional, different microstates can be generated by slightly varying the data provided to the model.

36 MATERIALS SCIENCE

Optimization of a 10–Step Reduced CH₄/O₂ Combustion Mechanism for RDRE Applications

A 10-step, 9-species reduced mechanism for methane/oxygen combustion is evaluated and optimized for use in rotating detonation rocket engine (RDRE) simulations at an operating pressure of 10 bar. The mechanism is assessed against GRI-Mech 3.0 for ignition delay times, adiabatic equilibrium temperature and species, Chapman-Jouguet (CJ) detonation parameters, and laminar flame speed. A sequential Nelder–Mead optimization procedure is applied to the rate constants of reactions R1, R4, and R5, targeting ignition delay agreement over the equivalence ratio range φ = 1.0–1.4 and temperature range T = 1400–2000 K. The optimized mechanism achieves ignition delay agreement within 28% of GRI-Mech 3.0 for the temperature range relevant to RDRE operation (τ ≥ 5 μs). Adiabatic equilibrium temperatures and CJ detonation parameters are reproduced within 0.15% of GRI-Mech 3.0 values, confirming thermodynamic consistency of the mechanism. Laminar flame speed predictions deviate significantly from GRI-Mech 3.0, which is attributed to the absence of HO₂ chemistry in the reduced mechanism; this limitation does not affect RDRE simulation fidelity since laminar flame propagation is of secondary importance in RDREs.

Computational Fluid Dynamics

A Laboratory Breadboard System for Dual-Arm Teleoperation

The computing architecture of a novel dual-arm teleoperation system is described. The novelty of this system is that: (1) the master arm is not a replica of the slave arm; it is unspecific to any manipulator and can be used for the control of various robot arms with software modifications; and (2) the force feedback to the general purpose master arm is derived from force-torque sensor data originating from the slave hand. The computing architecture of this breadboard system is a fully synchronized pipeline with unique methods for data handling, communication and mathematical transformations. The computing system is modular, thus inherently extendable. The local control loops at both sites operate at 100 Hz rate, and the end-to-end bilateral (force-reflecting) control loop operates at 200 Hz rate, each loop without interpolation. This provides high-fidelity control. This end-to-end system elevates teleoperation to a new level of capabilities via the use of sensors, microprocessors, novel electronics, and real-time graphics displays. A description is given of a graphic simulation system connected to the dual-arm teleoperation breadboard system. High-fidelity graphic simulation of a telerobot (called Phantom Robot) is used for preview and predictive displays for planning and for real-time control under several seconds communication time delay conditions. High fidelity graphic simulation is obtained by using appropriate calibration techniques.

A K Bejczy

Modeling Occupant Core Temperatures Across the Boston Building Stock to Advance Public Health

In recent history, extreme heat has been the cause of most deaths from a natural disaster. Exposure to extreme heat can aggravate preexisting conditions, increase hospitalization, and even cause death. Thus, modeling the thermal resilience of households across the United States will allow for a quantitative assessment of the health and safety risks posed by extreme heat. For the first time, we simulate occupant comfort in representative households across Boston by combining the granular results of the ResStock(TM) model with a two-node heat strain model. This model calculates occupants' core body temperature in each simulated household over a year. We compare the simulated core temperatures to two public health metrics: hyperthermia and heat stroke. Our results show for households in Boston that do not have or use aid conditioning, thousands potentially experience many dangerous heat events each summer and these events can last for more than a day at a time. These heat events peak during the late afternoon and evening, just as residents are coming home, cooking meals, and trying to go to sleep. We have found that multi-family buildings, renters, and low-income homes in more airtight and insulated homes are the most at risk for these events. These results demonstrate the scale and urgency of exposure to extreme heat.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)

Knowledge-guided learning with curated prior genetic biomarkers for robust model interpretation

Abstract Motivation Knowledge-guided learning offers effective and robust model training strategies in data-scarce settings by incorporating established domain knowledge, thereby enhancing generalization, robustness, and interpretability. By contrast, conventional deep learning approaches rely purely on data-driven learning, which can limit robust model interpretability, particularly in high-dimensional settings with limited size samples. In computational biology, knowledge-guided learning has primarily leveraged network- and structural-based knowledge, leading to biologically interpretable representations and enhanced predictive performance compared to conventional approaches. However, curated biomarkers, one of the most accessible forms of biological knowledge, remain largely unexplored within knowledge-guided paradigms. Results In this study, we propose a model-agnostic training paradigm, Biomarker-driven Explainable Prior-guided Learning (BioExPL), that can be applied to any neural networks that incorporates curated prior knowledge. BioExPL enforces neural networks to reflect curated biomarker priors in their latent representations through a novel knowledge-alignment loss. BioExPL consistently demonstrated significantly improved predictive performance and enhanced model interpretability with minimized computational overhead in simulation studies and intensive experiments on multiple cancer datasets. BioExPL not only integrates prior curated knowledge into the model but also accurately identifies unknown associated signals additionally. BioExPL is model-agnostic and domain-independent, enabling its integration into diverse neural network architectures. Availability and implementation The open-source is publicly available at: https://github.com/datax-lab/BioExPL.

Baek, Beomsu [Department of Computer Science, Univ

Micrometer: Micromechanics transformer for predicting full field mechanical responses of heterogeneous materials

Predicting mechanical responses of heterogeneous materials across scales remains a significant challenge. Traditional computational methods often struggle with complex and multiscale nature of these materials, limiting their effectiveness in real-world applications. Here, in this paper, we introduce Micrometer, a vision transformer based deep learning model designed to predict full field mechanical responses of heterogeneous materials, bridging the gap between computer vision and solid mechanics problems. We show that Micrometer, trained on a large-scale high-resolution dataset of 2D fiber-reinforced composites, can achieve state-of-the-art performance in predicting microscale strain fields across a wide range of material properties and loading conditions. Our model demonstrates accuracy and computational efficiency in applications such as computational homogenization and multiscale modeling, reducing computational time by up to two orders of magnitude compared to conventional numerical solvers while maintaining less than 1 % errors in predicting macroscale stress fields. Furthermore, we showcase Micrometer’s adaptability through transfer learning experiments on new materials with limited data, highlighting its potential to tackle diverse scenarios in computational solid mechanics. These results represent a significant step towards AI-driven innovation in materials science, addressing the limitations of traditional numerical methods and paving the way for more efficient simulations of heterogeneous materials across various industrial applications.

Composite materials

Entanglement properties of SU(2) gauge theory

Understanding how isolated quantum systems thermalize is central to both fundamental physics and the development of quantum technologies. In this Perspective, we present recent and new results on thermalization in nonabelian gauge theory from exact real-time simulations of two-dimensional lattice models. We discuss tests of the eigenstate thermalization hypothesis, entanglement entropy dynamics, the absence of persistent non-thermal states, spectral signatures, and measures of quantum complexity. Our new results include a two-step thermalization process in localized regions and the role of higher gauge field representations. These findings suggest that thermalization in gauge theories may offer a test case for quantum advantage.

information theory and computation

Nuclear Safety [Vol. 36, No. 2, July-December 1995]

Nuclear Safety is a journal that covers significant issues in the field of nuclear safety. Its primary scope is safety in the design, construction, operation, and decommissioning of nuclear power reactors worldwide and the research and analysis activities that promote this goal, but it also encompasses the safety aspects of the entire nuclear fuel cycle, including fuel fabrication, spent-fuel processing and handling, and nuclear waste disposal, the handling of fissionable materials and radioisotopes, and the environmental effects of all these activities. Table of Contents for this issue follows. THE CHORNOBYL ACCIDENT: 195 The Chornobyl Accident Revisited, Part III: Chernobyl Source Term Release Dynamics and Reconstruction of Events During the Active Phase, A. R. Sich; GENERAL SAFETY CONSIDERATIONS: 218 Second ANS Workshop on the Safety of Soviet-Designed Nuclear Power Plants, R. A. Bari; 234 Elements of a Nuclear Criticality Safety Program, C. M. Hopper; 243 Rickover, Excellence, and Criticality Safety Programs, R. E. Wilson; ACCIDENT ANALYSIS: 249 Transient Analysis of the PIUS Advanced Reactor Design with the TRAC-PF1/MOD2 Code, B. E. Boyack, J. L Steiner, S. C. Harmony, H. J. Stumpf, and J. F. Lime; 278 The Hierarchy-By-Interval Approach to Identifying Important Models that Need Improvement in Severe-Accident Simulation Codes, T. J. Heames, M. Khatib-Rahbar, J. E. Kelly, R. P. Jenks-Johnson, and Y.-S. Chen; 290 RELAP5/MOD3 Code Coupling Model, R. P. Martin; 299 Missiles Caused by Severe Pressurized-Water Reactor Accidents, R. Krieg; DESIGN FEATURES: 310 Validation of COMMIX with Westinghouse AP-600 PCCS Test Data, J. G. Sun, T. H. Chien, J. Ding, and W T. Sha; ENVIRONMENTAL EFFECTS: 321 Spent Nuclear Fuel Characterization for a Bounding Reference Assembly for the Receiving Basin for Off-Site Fuel, S. D. Kahook, R. L. Garrett, L. R. Canas, and M J. Beckum; OPERATING EXPERIENCES: 332 Reactor Shutdown Experience, Compiled by J. W. Cletcher; U.S. NUCLEAR REGULATORY COMMISSION INFORMATION AND ANALYSES: 335 Reactor Coolant System Blowdown at Wolf Creek on September 17, 1994, J. V. Kauffman and S. L. Israel; RECENT DEVELOPMENTS: 344 Reports, Standards, and Safety Guides, D. S. Oueener; 349 Proposed Rule Changes as of June 30, 1995; ANNOUNCEMENTS: 320 Symposium on Acceptability of Risk From Radiation—Application to Manned Space Flight; 320 24th DOE/NRC Nuclear Air Cleaning and Treatment Conference; 361 Radiation Biology and Radiation Protection— Modern Developments and Tendencies in Radiation Biology; 362 1997 IEEE Sixth Conference on Human Factors and Power Plants; 354 The Authors; 360 Reviewers of Nuclear Safety, Vol. 36.

05 NUCLEAR FUELS

Preparation of the Ice Certification of the Dornier 328 Regional Airliner By Numerical Simulation and By Ground Test

The Dornier 328, a new regional airliner, is to be qualified to FAR25/JAR25 requirements for operation into known icing conditions. All lifting surfaces are pneumatically deiced. Impingement limits were determined using a validated computational fluid dynamic (CFD) code. Wind tunnel tests with a model of the whole aircraft and the isolated empennage equipped with artificial ice shapes showed a degradation of handling characteristics due to ice accretion. Extensive two and three dimensional CFD calculations evaluated ice accretion at the deflection elevator horn. Icing tunnel tests were performed with a model of the horizontal tail with icing similitude fulfilled. The Dornier Do328 powerplant air induction system has successfully completed a program of icing tests in support of aircraft certification for operation into known icing conditions. Potential flight icing conditions were extensively analyzed using a CFD potential flow code and 3-D particle tracking routine to identify critical flight icing conditions. The tunnel test program verified the analytically predicted critical ice accretion surfaces of the air induction system, the adequacy of the ice protection provisions and demonstrated compliance with the applicable (JAR25) regulatory material.

D Welte