Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer experiments”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Expansion dynamics of femtosecond laser-induced plasmas: Influence of thermophysical plasma properties

This study investigates the expansion dynamics of femtosecond laser-induced plasmas, emphasizing the impact of plasma thermophysical properties and ambient gas composition. Through shadowgraphy experiments and multiphase computational fluid dynamics (CFD) simulations, the influence of parameters such as heat capacity, molecular weight, and thermal conductivity on plume morphology, shockwave evolution, and energy dissipation mechanisms is examined. A mixture multiphase model is implemented to capture the interaction between the plasma and the surrounding gas. Simulation results reveal that plasma expansion is strongly inertia-driven. Results show that differences in plasma properties and ambient conditions affect the shape and temperature distribution of the expanding plume. The early-stage dynamics are primarily dictated by pressure forces, whereas thermal and viscous effects play a growing role in the plume's behavior during later stages of expansion. The CFD findings show the necessity of accurate initial condition characterization, including crater geometry and plasma pressure and temperature, for reliable modeling of plasma evolution in laser ablation processes.

CFD modeling↗
DOE OSTI2022

Charting C–C coupling pathways in electrochemical CO 2 reduction on Cu(111) using embedded correlated wavefunction theory

The electrochemical CO 2 reduction reaction (CO 2 RR) powered by excess zero-carbon-emission electricity to produce especially multicarbon (C 2+ ) products could contribute to a carbon-neutral to carbon-negative economy. Foundational to the rational design of efficient, selective CO 2 RR electrocatalysts is mechanistic analysis of the best metal catalyst thus far identified, namely, copper (Cu), via quantum mechanical computations to complement experiments. Here, we apply embedded correlated wavefunction (ECW) theory, which regionally corrects the electron exchange-correlation error in density functional theory (DFT) approximations, to examine multiple C–C coupling steps involving adsorbed CO (*CO) and its hydrogenated derivatives on the most ubiquitous facet, Cu(111). We predict that two adsorbed hydrogenated CO species, either *COH or *CHO, are necessary precursors for C–C bond formation. The three kinetically feasible pathways involving these species yield all three possible products: *COH–CHO, *COH–*COH, and *OCH–*OCH. The most kinetically favorable path forms *COH–CHO. In contrast, standard DFT approximations arrive at qualitatively different conclusions, namely, that only *CO and *COH will prevail on the surface and their C–C coupling paths produce only *COH–*COH and *CO–*CO, with a preference for the first product. This work demonstrates the importance of applying qualitatively and quantitatively accurate quantum mechanical method to simulate electrochemistry in order ultimately to shed light on ways to enhance selectivity toward C 2+ product formation via CO 2 RR electrocatalysts.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Photoselective isotope fractionation dynamics of N 2 with cosmo and atmospheric chemistry perspectives

Stable isotope ratio measurements provide valuable insights into a broad range of natural processes, from planetary atmospheres and climate to interstellar chemistry. Nitrogen, which has two stable isotopes, exhibits varying isotope ratios across the solar system. To model these observations, the isotope fraction as a function of energy is essential. At the Advanced Light Source (ALS), we measured the photodissociation of molecular nitrogen (N 2 ) with vacuum UV photons where a single photon is sufficiently energetic to dissociate the strong bond. The nitrogen atoms produced are scavenged with H 2 to form ammonia, whose isotopic makeup is determined. Blending the experiments with dynamical computations that include the shielding of light, we examine the isotopic composition and electronic atomic states produced. The measured photodissociation of N 2 at a natural isotopic composition with a frequency broad light beam exceptionally strongly favors the formation of the heavier nitrogen isotope, 15 N. Computations concur and suggest that the maximum in the quantum yield reflects significant variations in the specific electronic quantum states of the product N atoms that have quite different reactivities. Our quantum computations show that at similar energies, photodissociation of 14 N 14 N and 15 N 14 N can lead to different product channels. The computed dynamics include extensive state-selective spin–orbit and nonadiabatic couplings affecting the light absorption and dissociation pathways that proceed via the triplet manifold of states. Our results are relevant for future exploration missions, both in situ and sample-return and for other molecules such as O 2 and CO.

Komarova, Ksenia [Hebrew Univ. of Jerusalem (Israe↗
DOE OSTI2025

Exploring the relation between transonic dislocation glide and stacking fault width in FCC metals

Theory predicts limiting gliding velocities that dislocations cannot overcome. Computational and recent experiments have shown that these limiting velocities are soft barriers and dislocations can reach transonic speeds in high rate plastic deformation scenarios. In this paper we systematically examine the mobility of edge and screw dislocations in several face centered cubic (FCC) metals (Al, Au, Pt, and Ni) in the extreme large-applied-stress regime using molecular dynamics simulations. Our results show that edge dislocations are more likely to move at transonic velocities due to their high mobility and lower limiting velocity than screw dislocations. Importantly, among the considered FCC metals, the dislocation core structure determines the dislocation’s ability to reach transonic velocities. This is likely due to the variation in stacking fault width due to relativistic effects near the limiting velocities.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Granger causal inference for climate change attribution

Abstract Climate change detection and attribution (D&A) is concerned with determining the extent to which anthropogenic activities have influenced specific aspects of the global climate system. D&A fits within the broader field of causal inference, the collection of statistical methods that identify cause and effect relationships. There are a wide variety of methods for making attribution statements, each of which require different types of input data and focus on different types of weather and climate events and each of which are conditional to varying extents. Some methods are based on Pearl causality (direct experimental interference) while others leverage Granger (predictive) causality, and the causal framing provides important context for how the resulting attribution conclusion should be interpreted. However, while Granger-causal attribution analyses have become more common, there is no clear statement of their strengths and weaknesses relative to Pearl-causal attribution and no clear consensus on where and when Granger-causal perspectives are appropriate. In this prospective paper, we provide a formal definition for Granger-based approaches to trend and event attribution and a clear comparison with more traditional methods for assessing the human influence on extreme weather and climate events. Broadly speaking, Granger-causal attribution statements can be constructed quickly from observations and do not require computationally-intesive dynamical experiments. These analyses also enable rapid attribution, which is useful in the aftermath of a severe weather event, and provide multiple lines of evidence for anthropogenic climate change when paired with Pearl-causal attribution. Confidence in attribution statements is increased when different methodologies arrive at similar conclusions. Moving forward, we encourage the D&A community to embrace hybrid approaches to climate change attribution that leverage the strengths of both Granger and Pearl causality.

Risser, Mark D. (ORCID:0000000319561783)↗
DOE OSTI2019

Electronic spin transition in FeO 2 : Evidence for Fe(II) with peroxide O 2 2 –

The discovery of FeO 2 containing more oxygen than hematite (Fe 2 O 3 ), which was previously believed to be the most oxygen rich iron compound, has important implications for the study of deep lower mantle compositions. Compared to other iron compounds, there are limited reports on FeO 2 , making studies of its physical properties of great interest in fundamental condensed matter physics and geoscience. Even the oxidation state of Fe in FeO 2 is the subject of debate in theoretical works and there have not been reports from experimental electronic and magnetic properties measurements. Here, we report the pressure-induced spin state transition from synchrotron experiments and our computational results explain the underlying mechanism. Using density functional theory and dynamical mean field theory, we calculated spin states of Fe with volume and Hubbard interaction U change, which clearly demonstrate that Fe in FeO 2 consists of Fe(II) and peroxide O 2– 2 . In conclusion, our paper suggests that the localized nature of both Fe 3d orbitals and O 2 molecular orbitals should be correctly treated for unveiling the structural and electronic properties of FeO 2 .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Increased localization of Majorana modes in antiferromagnetic chains on superconductors

Magnet-superconductor hybrid (MSH) systems are a key platform for custom-designed topological superconductors. Ideally, the ends of a one-dimensional MSH structure will host Majorana zero-modes (MZMs), the fundamental unit of topological quantum computing. However, some experiments with ferromagnetic (FM) chains show a more complicated picture. Due to tiny gap sizes and hence long coherence lengths, MZMs might hybridize and lose their topological protection. Recent experiments on a niobium surface have shown that both FM and antiferromagnetic (AFM) chains may be engineered, with the magnetic order depending on the crystallographic direction of the chain. While FM chains are well understood, AFM chains are less so. Here, we study two models inspired by the niobium surface: A minimal model to elucidate the general topological properties of AFM chains and an extended model to more closely simulate a real system by mimicking the proximity effect. Here, we find that, in general, for AFM chains, the topological gap is larger than for FM ones, and thus, coherence lengths are shorter for AFM chains, yielding more pronounced localization of MZMs in these chains. While for some parameters AFM chains may be topologically trivial, we find in these cases that adding an adjacent chain can result in a nontrivial system, with a single MZM at each chain end.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Brittle-ductile transitions in a metallic glass

Recent computational and laboratory experiments have shown that brittle-ductile transitions in the notch toughnesses of metallic glasses such as Vitreloy 1 are strongly sensitive to the initial effective disorder (or “fictive”) temperature. Glasses with lower effective temperatures are weak and brittle; those with higher effective temperatures are strong and ductile. The analysis of this phenomenon presented here examines the onset of fracture at the tip of a notch as predicted by the shear-transformation-zone theory of spatially varying plastic deformation. The central ingredient of this analysis is an approximation for the dynamics of the plastic zone formed by stress concentration at the notch tip. As a result, this zone first shields the tip but then, with increasing stress, expands suddenly, producing a discontinuous transition between brittle and ductile failure in satisfactory agreement with the numerical and experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Path integral Monte Carlo simulations of warm dense aluminum

Here, we perform first-principles path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of aluminum. Our equation of state (EOS) simulations cover a wide density-temperature range of $0.1-32.4$~g$\,$cm$^{-3}$ and $10^4-10^8$~K. Since PIMC and DFT-MD accurately treat effects of the atomic shell structure, we find two compression maxima along the principal Hugoniot curve attributed to K-shell and L-shell ionization. The results provide a benchmark for widely used EOS tables, such as SESAME, QEOS, and models based on Thomas-Fermi and average-atom techniques. A subsequent multi-shock analysis provides a quantitative assessment for how much heating occurs relative to an isentrope in multi-shock experiments. Finally, we compute heat capacity, pair-correlation functions, the electronic density of states, and Z to reveal the evolution of the plasma structure and ionization behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Shining light on cysteine modification: connecting protein conformational dynamics to catalysis and regulation

We report that cysteine is a rare but functionally important amino acid that is often subject to covalent modification. Cysteine oxidation plays an important role in many human disease processes, and basal levels of cysteine oxidation are required for proper cellular function. Because reactive cysteine residues are typically ionized to the thiolate anion (Cys-S - ), their formation of a covalent bond alters the electrostatic and steric environment of the active site. X-ray-induced photo-oxidation to sulfenic acids (Cys-SOH) can recapitulate some aspects of the changes that occur under physiological conditions. Here we propose how site-specific cysteine photo-oxidation can be used to interrogate ensuing changes in protein structure and dynamics at atomic resolution. Although this powerful approach can connect cysteine covalent modification to global protein conformational changes and function, careful biochemical validation must accompany all such studies to exclude misleading artifacts. Lastly, new types of X-ray crystallography experiments and powerful computational methods are creating new opportunities to connect conformational dynamics to catalysis for the large class of systems that use covalently modified cysteine residues for catalysis or regulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Towards a data-centric research and development roadmap for large-scale science user facilities

© 2017 IEEE. The U.S. Department of Energy (DOE) Office of Science (SC) operates approximately four dozen large-scale science user facilities (SUFs), each of which generates a tremendous amount of scientific data from experiments, observations and computations. To better understand the data needs and challenges, DOE has run many workshops in recent years to identify and articulate data-centric challenges and opportunities at varying resolution, from facility to community scale. Building on those workshop reports, as well as others from elsewhere in the community, this article goes beyond the findings-recommendations typical of workshop reports to consider how one might structure a broad, technology- and data-centric, coordinated research effort that would realize progress towards solutions that address the well documented challenges and opportunities. We focus on identifying practical issues of strategic relevance, along with offering a view about the focal points for a coordinated research and development effort that would target meeting data-centric needs of a broad set of science users and SUFs. These focal points would, by their nature, engage a spectrum of researchers from computer science, computational and experimental sciences, and data science in a coordinated fashion.

Bethel, EW↗
DOE OSTI2017

Processing-in-Memory Enabled Graphics Processors for 3D Rendering

The performance of 3D rendering of Graphics Processing Unit that convents 3D vector stream into 2D frame with 3D image effects significantly impact users’ gaming experience on modern computer systems. Due to the high texture throughput in 3D rendering, main memory bandwidth becomes a critical obstacle for improving the overall rendering performance. 3D stacked memory systems such as Hybrid Memory Cube (HMC) provide opportunities to significantly overcome the memory wall by directly connecting logic controllers to DRAM dies. Based on the observation that texel fetches significantly impact off-chip memory traffic, we propose two architectural designs to enable Processing-In-Memory based GPU for efficient 3D rendering.

Innovative hardware/software co-design↗
DOE OSTI2024

Future of plasma etching for microelectronics: Challenges and opportunities

Plasma etching is an essential semiconductor manufacturing technology required to enable the current microelectronics industry. Along with lithographic patterning, thin-film formation methods, and others, plasma etching has dynamically evolved to meet the exponentially growing demands of the microelectronics industry that enables modern society. At this time, plasma etching faces a period of unprecedented changes owing to numerous factors, including aggressive transition to three-dimensional (3D) device architectures, process precision approaching atomic-scale critical dimensions, introduction of new materials, fundamental silicon device limits, and parallel evolution of post-CMOS approaches. The vast growth of the microelectronics industry has emphasized its role in addressing major societal challenges, including questions on the sustainability of the associated energy use, semiconductor manufacturing related emissions of greenhouse gases, and others. The goal of this article is to help both define the challenges for plasma etching and point out effective plasma etching technology options that may play essential roles in defining microelectronics manufacturing in the future. The challenges are accompanied by significant new opportunities, including integrating experiments with various computational approaches such as machine learning/artificial intelligence and progress in computational approaches, including the realization of digital twins of physical etch chambers through hybrid/coupled models. These prospects can enable innovative solutions to problems that were not available during the past 50 years of plasma etch development in the microelectronics industry. To elaborate on these perspectives, the present article brings together the views of various experts on the different topics that will shape plasma etching for microelectronics manufacturing of the future.

Engineering↗
DOE OSTI2021

Extraction of the specific shear viscosity of quark-gluon plasma from two-particle transverse momentum correlations

The specific shear viscosity, η/s, of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions at RHIC and LHC is estimated based on the progressive longitudinal broadening of transverse momentum two-particle correlators, G 2 , reported as a function of collision centrality by the STAR and ALICE experiments. Estimates are computed as a function of collision centrality using the Gavin ansatz which relates the G 2 longitudinal broadening to the specific shear viscosity. Freeze out times required for the use of the ansatz are computed using a linear fit of freeze out times reported as a function of the cubic root of the charged particle pseudorapidity density (dN ch /dη) 1/3 . Estimates of η/s based on ALICE data exhibit little to no dependence on collision centrality at LHC energy, while estimates obtained from STAR data hint that η/s might be a function of collision centrality at top RHIC energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

CGSim: A Simulation Framework for Large Scale Distributed Computing Environment

Large-scale distributed computing infrastructures such as the Worldwide LHC Computing Grid (WLCG) require comprehensive simulation tools for evaluating performance, testing new algorithms, and optimizing resource allocation strategies. However, existing simulators suffer from limited scalability, hardwired algorithms, lack of real-time monitoring, and inability to generate datasets suitable for modern machine learning approaches. We present CGSim, a simulation framework for large-scale distributed computing environments that addresses these limitations. Built upon the validated SimGrid simulation framework, CGSim provides high-level abstractions for modeling heterogeneous grid environments while maintaining accuracy and scalability. Key features include a modular plugin mechanism for testing custom workflow scheduling and data movement policies, interactive real-time visualization dashboards, and automatic generation of event-level datasets suitable for AI-assisted performance modeling. We demonstrate CGSim’s capabilities through a comprehensive evaluation using production ATLAS PanDA workloads, showing significant calibration accuracy improvements across WLCG computing sites. Scalability experiments show near-linear scaling for multi-site simulations, with distributed workloads achieving 6 × better performance compared to single-site execution. The framework enables researchers to simulate WLCG-scale infrastructures with hundreds of sites and thousands of concurrent jobs within practical time budget constraints on commodity hardware.

Vatsavai, Sairam Sri [Brookhaven National Laborato↗
DOE OSTI2025

Relativistic Magnetic Reconnection in Astrophysical Plasmas: A Powerful Mechanism of Nonthermal Emission

Magnetic reconnection—a fundamental plasma physics process, where magnetic field lines of opposite polarity annihilate—is invoked in astrophysical plasmas as a powerful mechanism of nonthermal particle acceleration, able to explain fast-evolving, bright high-energy flares. Near black holes and neutron stars, reconnection occurs in the relativistic regime, in which the mean magnetic energy per particle exceeds the rest mass energy. This review reports recent advances in our understanding of the kinetic physics of relativistic reconnection (RR): ▪ Kinetic simulations have elucidated the physics of plasma heating and nonthermal particle acceleration in RR. ▪ The physics of radiative RR, with its self-consistent interplay between photons and reconnection-accelerated particles—a peculiarity of luminous, high-energy astrophysical sources—is the new frontier of research. ▪ RR plays a key role in global models of high-energy sources, in terms of both global-scale layers and reconnection sites generated as a by-product of local magnetohydrodynamic instabilities. We summarize themes of active investigation and future directions, emphasizing the role of upcoming observational capabilities, laboratory experiments, and new computational tools.

acceleration of particles↗
DOE OSTI2024

Cloud-Resolving ICON Simulations of Secondary Ice Production in Arctic Mixed-Phase Stratocumuli Observed during M-PACE

Field measurements and modeling studies suggest that secondary ice production (SIP) may close the gap between observed Arctic ice nucleating particle (INP) concentrations and ice crystal number concentrations ni. Here, we explore sensitivities with respect to the complexity of different INP parameterizations under the premise that ni is governed by SIP. Idealized, cloud-resolving simulations are performed for the marine cold air outbreak cloud deck sampled during the Mixed-Phase Arctic Cloud Experiment (M-PACE) with the Icosahedral Nonhydrostatic (ICON) model. The impact of the droplet shattering (DS) of raindrops and collisional breakup (BR) in addition to the existing Hallet–Mossop rime splintering mechanism were investigated. Overall, 12 different model experiments (12-h runs) were performed and analyzed. Despite the considerable amount of uncertainty remaining with regard to SIP mechanisms and their process representation in numerical models, we conclude from these experiments that (i) only simulations where DS dominates the SIP signal (potentially amplified by BR) capture observed ice-phase and liquid-phase cloud properties, and (ii) SIP events cluster around the convective outflow region and are structurally linked to mesoscale cloud organization. In addition, interactions with primary nucleation parameterizations of varied complexity were investigated. Here, our simulations show that (i) a stable long-lived mixed-phase cloud (MPC) can be maintained in the absence of primary nucleation once SIP is established, (ii) experiments using a computationally more efficient relaxation-based parameterization of primary nucleation are statistically invariant from simulations considering prognostic INP, and (iii) primary nucleation at cloud-top controls the areal extent of the mixed-phase cloud region, and reduces SIP efficacy via DS due to increased depletion of cloud liquid throughout the entire cloud column.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

A Laser-Based Optical Sensor for Broad-Band Measurements of Building Earthquake Drift

Accurate measurements of the time-dependent deformations of a building during earthquake excitation are essential for interpretation of the dynamic response of the as-built system and for quantifying the seismic demands. Traditional approaches for monitoring building systems are based on strong motion accelerometers mounted at selected elevations. However, accelerometer-based systems do not directly measure the deformations of the structure, and can have significant limitations that make it challenging to correctly measure deformations, particularly permanent deformations from inelastic response. In the study described herein, computational simulations and experiments were combined to evaluate the potential of a new optically based sensor to directly measure time-dependent deformations of a building, including inelastic deformations. The sensor methodology includes corrections for localized structural member rotations and can provide estimates of the absolute accelerations at each floor. A laser-based system utilizing a recently developed discrete diode position sensor (DDPS) is evaluated, and the ability of such a system to measure earthquake induced transient deformations characterized by building interstory drift is demonstrated.

58 GEOSCIENCES↗