SEARCH · Search NASA

Results for “Computational Physics and 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

3,226 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2022

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↗

DOE OSTI2025

Exploiting ν-dependence of projected energy correlators in HICs

We extend the recently derived factorization formula for energy-energy correlators to study the analytic structure of general ν-point projected energy correlators in heavy ion collisions. The ν-point projected energy correlators (or, ν-correlators) are an analytically continued family of the integer N-point projected energy correlators, which probe correlations between N final-state particles. By tracking the largest separation (χ) between the N particles, in vacuum, their structure is closely related to the DGLAP splitting functions and exhibits a classical scaling behavior ∼ 1/χ which is modified by resummation through the anomalous dimensions. We show that, in a thermal medium, the ν-correlators display non-trivial angular scaling already at the leading order in perturbation theory. We find that for non-integer values, particularly ν < 1, medium-induced jet function is enhanced compared to ν > 1. This is particularly manifested in the ratios of ν-correlators with respect to the two-point energy correlator, which encodes intrinsic angular information for ν < 1 when compared to large ν values. Moreover, for small-ν values, the ν-correlators appear to saturate at ν = 0.01 . We further confirm our leading-order numerical computations against simulated events from JEWEL for the parton-level production cross-section. Finally, we qualitatively discuss the effect of BFKL resummation for various values of ν.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Surface phase diagrams from nested sampling

From nested sampling, we compute the partition function and, from that, the phase diagram of gas adsorbates, including their anharmonic and configurational degrees of freedom, on flat and stepped surfaces of the Lennard-Jones solid.

Chemistry↗

DOE OSTI2026

Role of ion acoustic instability in magnetic reconnection

We report on a first-principles numerical study of magnetic reconnection in plasmas with different initial ion-to-electron temperature ratios. In cases where this ratio is significantly below unity, we observe intense wave activity in the diffusion region, driven by the ion-acoustic instability. Our analysis shows that the dominant macroscopic effect of this instability is to drive substantial ion heating. In contrast to earlier studies reporting significant anomalous resistivity, we find that anomalous contributions due to the ion-acoustic instability are minimal. These results shed light on the dynamical impact of this instability on reconnection processes, offering new insights into the fundamental physics governing collisionless reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2026

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↗

NASA NTRSDate not supplied

Natural Fiber Simulations and Analysis for Composites

As the U.S. government and the aviation industry strive to achieve net-zero carbon emissions by 2050, the Sustainable Manufacturing of Aircraft (SUMAC) project, supported by the Convergent Aeronautics Solutions (CAS) program of the National Aeronautics and Space Administration (NASA), aims to develop sustainably derived thermoplastic composites and manufacturing technologies for future aviation. A multidisciplinary team of experts drives this initiative through a collective effort. The team’s key focus areas are sustainably derived thermoplastic resins, natural fiber reinforcement, composite manufacturing, structural health monitoring, computational materials modeling, and systems analysis. This presentation will give an overview of the SUMAC project. Fiber-scale microscopy analysis characterized fiber cross sections finding that hemp and flax fibers are very deformable and non-circular (circularities of 0.2 were observed), and diameters range from submicron to 20 μm. Fiber shapes range from flattened, with aspect ratios Rlong/Rshort up to 5, to round, Rlong/Rshort =1. And Fiber core perimeters, which likely drives deformability, were calculated Pcore=7 μm. Explicit-fiber yarn simulations, parameterized from single-fiber measurements, resolved and identified slack as a cause for yarn strain-hardening tensile behavior. The bonded network model of PLA, PHA and PA11 resins captures tensile behavior, and can model an infused weave. Computational models of the weave unit cells matched the pattern (including plain, twill and triaxial braid), yarn ply (1to 2), and yarn twist from optical microscopy.

thermoplastic↗

NASA NTRSDate not supplied

Natural Fiber Simulations and Analysis for Composites

thermoplastic↗

NASA NTRS2026

NASA JSC’s Simulant Development Lab Capabilities and Artemis Testing

The Simulant Development Lab (SDL) is a multifunctional collaborative workspace that supports the development, curation, analysis, testing, and distribution of planetary regolith simulants – including lunar, Martian, asteroidal, and other granular materials. The lab provides a multidisciplinary setting for scientific characterization of simulant physical properties and for engineering evaluations conducted with simulant test beds. To enable this work, the SDL curates and maintains a stock of more than 35 metric tons of simulant material. To evaluate these materials and support testing goals, the lab is equipped with a comprehensive suite of processing tools and analytical instruments. These capabilities enable the SDL’s mission at NASA’s Johnson Space Center to distribute, develop, process, characterize, and test regolith simulants for mission relevant applications. Through controlled and repeatable testing environments that replicate the physical and compositional properties of lunar regolith, the SDL supports Artemis hardware maturation, providing safe, Earth‑based analogs for evaluating systems that must withstand regolith dust interactions, physical wear and abrasion, and operational loads. The facility’s extensive simulant inventory and integrated geological and engineering test infrastructure accelerate technology readiness for Artemis and future exploration campaigns (e.g., future crewed or robotic missions to Mars).

Simulant Development Lab↗

DOE OSTI2026

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗

DOE OSTI2025

Carbon in GaN as a nonradiative recombination center

Trap-assisted nonradiative recombination has been shown to limit the efficiency of optoelectronic devices. While substitutional carbon (⁠C N ) has been suggested to be a nonradiative recombination center in GaN devices, a complete recombination cycle including the two charge-state transition levels has not been previously described. In this work, we investigate the trap-assisted recombination process due to C N in GaN, including multiphonon emission, radiative recombination, trap-assisted Auger–Meitner (TAAM) recombination as well as thermal emission of holes. Our study shows the key role of TAAM processes at the high carrier densities relevant for devices. We also reveal the carrier-density regimes where thermal emission and radiative recombination are expected to play an observable role. Our results highlight that carbon concentrations exceeding ~10 17 cm −3 can have a noticeable impact on device efficiency, not just in GaN active layers but also in InGaN and AlGaN. Furthermore, our comprehensive formalism not only offers detailed results for carbon but also provides a general framework for assessing the multiple processes that participate in trap-assisted recombination in semiconductors.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗

NASA NTRS1990

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↗

NASA NTRS1995

A Case Study of View-Factor Rectification Procedures for Diffuse-Gray Radiation Enclosure Computations

The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.

Robert P Taylor↗

DOE OSTI2026

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↗

DOE OSTI2024

Copacabana: a probabilistic membership assignment method for galaxy clusters

Cosmological analyses using galaxy clusters in optical/near-infrared photometric surveys require robust characterization of their galaxy content. Precisely determining which galaxies belong to a cluster is crucial. In this paper, we present the COlor Probabilistic Assignment of Clusters And BAyesiaN Analysis (Copacabana) algorithm. Copacabana computes membership probabilities for all galaxies within an aperture centred on the cluster using photometric redshifts, colours, and projected radial probability density functions. We use simulations to validate Copacabana and we show that it achieves up to 89 per cent membership accuracy with a mild dependence on photometric redshift uncertainties and choice of aperture size. We find that the precision of the photometric redshifts has the largest impact on the determination of the membership probabilities followed by the choice of the cluster aperture size. We also quantify how much these uncertainties in the membership probabilities affect the stellar mass–cluster mass scaling relation, a relation that directly impacts cosmology. Using the sum of the stellar masses weighted by membership probabilities (⁠μ * ⁠) as the observable, we find that Copacabana can reach an accuracy of 0.06 dex in the measurement of the scaling relation at low redshift for a Legacy Survey of Space and Time type survey. These results indicate the potential of Copacabana and μ * to be used in cosmological analyses of optically selected clusters in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

NanoPSD: A software for automatic detection of Nano-Particle Shape Distribution in electron microscopy images

Accurate quantification of the size and morphology of nanoparticles from electron microscopy (EM) images is essential to understand growth mechanisms, surface reactivity, and functional behavior in nanoscale materials. Manual analysis remains slow, subjective, and difficult to reproduce in large datasets. We introduce NanoPSD (Nano-Particle Shape Distribution), an open-source and fully automated framework for quantitative particle detection and morphology analysis from EM images. NanoPSD integrates adaptive contrast enhancement, polarity-agnostic scale-bar detection, Optical Character Recognition (OCR)-based calibration, and classical segmentation via Otsu thresholding with morphological refinement. Particle contours are used to extract geometric descriptors, including equivalent circular diameter, aspect ratio, circularity, and solidity, enabling automated classification into spherical, rod-like, and aggregate morphologies. The framework supports both single-image and batch processing, generating publication-quality visualizations, LaTeX-ready tables, and structured comma-separated values (CSV) datasets. As a demonstration, we applied NanoPSD to plasma-synthesized nanoparticle samples diagnosed via transmission electron microscopy (TEM). The code produced statistically robust size and morphology distributions spanning a few to tens of nanometers with minimal user supervision. The pipeline demonstrates high reproducibility and scalability, processing large image collections with consistent calibration and output formatting. Its modular design enables seamless integration of future deep-learning-based segmentation models, providing a pathway toward intelligent, data-driven electron microscopy analysis.

36 MATERIALS SCIENCE↗