Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer experiment”

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,027 records · Page 57

A Gaussian process based surrogate approach for the optimization of cylindrical targets

Simulating direct-drive inertial confinement experiments presents significant computational challenges, both due to the complexity of the codes required for such simulations and the substantial computational expense associated with target design studies. Machine learning models, and in particular, surrogate models, offer a solution by replacing simulation results with a simplified approximation. In this study, we apply surrogate modeling and optimization techniques that are well established in the existing literature to one-dimensional simulation data of a new cylindrical target design containing deuterium–tritium fuel. These models predict yields without the need for expensive simulations. We find that Bayesian optimization with Gaussian process surrogates enhances sampling efficiency in low-dimensional design spaces but becomes less efficient as dimensionality increases. Nonetheless, optimization routines within two-dimensional and five-dimensional design spaces can identify designs that maximize yield, while also aligning with established physical intuition. Optimization routines, which ignore constraints on hydrodynamic instability growth, are shown to lead to unstable designs in 2D, resulting in yield loss. However, routines that utilize 1D simulations and impose constraints on the in-flight aspect ratio converge on novel cylindrical target designs that are stable against hydrodynamic instability growth in 2D and achieve high yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

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↗

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↗

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↗

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)↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Fossil Energy Workshop on Quantum Information Science & Technology (Summary Report)

This report summarizes findings and discussion from the Fossil Energy Workshop on Quantum Information Science & Technology which took place on November 19-21, 2019. Quantum Information Science (QIS) is expected to profoundly change the practice of science and engineering in the coming decades. It is a rapidly progressing field, fueled by large investments from the private sector and governments. Its importance to the U.S. economy and national security is underscored by the National Quantum Initiative Act passed in December 2018. QIS includes quantum sensing, quantum communications, quantum simulation experiments, and quantum computation. QIS technology exploits quantum phenomena for performing tasks that are impossible to do today, such as finding prime factors of large numbers or elucidating reaction mechanisms in complex chemical systems. The opportunities for applying QIS to problems encountered in fossil energy technology development are not known today. This workshop brought together, for the first time, experts in these fields to exchange information and explore potential research opportunities for QIS to advance fossil energy. The goal of the workshop was to develop a set of priority research opportunities that can inform future research efforts in QIS and build a community of next-generation researchers at the intersection of QIS and fossil energy.

20 FOSSIL-FUELED POWER PLANTS↗

USAMP Low-Cost Mg Sheet Component Development and Demonstration Project

The overarching objective of this USAMP project was to develop and demonstrate door panels made from magnesium (Mg) sheet with a cost penalty over conventional steel stampings of no more than $\$5.50$/kg saved. The technical approach integrated experiments with advanced computational tools based on Integrated Computational Materials Engineering (ICME) methods to develop new alloy chemistries and their thermomechanical processing that promise improved formability and lower forming temperatures. A penultimate task before finally forming the stampings was to incorporate actual microstructure into models that would enable formability simulations. This approach would, for the first time, account for individual magnesium grains moving in an anisotropic fashion unlike that for aluminum or steel that have isotropic properties upon which the current simulation tools are based. In separate activities, new coatings and lubricants to facilitate forming and improved corrosion protection and joining strategies, were developed to ensure that the door could be produced with stated product requirements. A technical cost model, which included parts production, assembly, and paint for a door specifically designed for Mg sheet, showed the cost penalty to be between $\$4.26$ to $\$6.60$/kg saved, which enveloped the project’s cost targets. The cost of the coated Mg sheet was identified as the key driver for the cost penalty. The mass of the Mg-intensive door was 7.9 kg, which was 54% less than the baseline steel door.

36 MATERIALS SCIENCE↗

Low-Cost Magnesium Alloy Sheet Component Development and Demonstration Project

The overarching objective of this USAMP project was to develop and demonstrate door panels made from magnesium (Mg) sheet with a cost penalty over conventional steel stampings of no more than $5.50/kg saved. The technical approach integrated experiments with advanced computational tools based on Integrated Computational Materials Engineering (ICME) methods to develop new alloy chemistries and their thermomechanical processing that promise improved formability and lower forming temperatures. A penultimate task before finally forming the stampings was to incorporate actual microstructure into models that would enable formability simulations. This approach would, for the first time, account for individual magnesium grains moving in an anisotropic fashion unlike that for aluminum or steel that have isotropic properties upon which the current simulation tools are based. In separate activities, new coatings and lubricants to facilitate forming and improved corrosion protection and joining strategies, were developed to ensure that the door could be produced with stated product requirements. A technical cost model, which included parts production, assembly, and paint for a door specifically designed for Mg sheet, showed the cost penalty to be between 4.26 USD to 6.60 USD/kg saved, which enveloped the project’s cost targets. The cost of the coated Mg sheet was identified as the key driver for the cost penalty. The mass of the Mg-intensive door was 7.9 kg, which was 54% less than the baseline steel door.

36 MATERIALS SCIENCE↗

Nanomechanics and Electronic Structure of Organic Photovoltaics in Real Application Conditions by Advanced Scanning Probe Microscopy (Final Technical Report for Project DE-SC0018041)

This project report consists of the following three sections: A.) Investigation of Charge Transport Properties in Semiconducting Conjugated Polymer Blends Using Computational Chemistry and Experiments: This section provides unpublished but valuable research results that correlate the electronic structure of PCBM:PCDTBT organic photovoltaic systems with their current-voltage characteristics. The presentation is structured in the form of a scientific publication. B.) Scanning Probe Microscopy Methodology Development Efforts for Organic Photovoltaic Materials: This section summarizes various published works focusing on advancements in scanning probe microscopy methods. Each sub-section corresponds to a scientific publication for which the abstract and citation are provided. C.) Other Publications Supported by the Project: This section lists additional published works partially supported by the project, which are not included in the previous two sections.

36 MATERIALS SCIENCE↗

Spatially Resolved Near Field Spectroscopy of Vibrational Polaritons at the Small N Limit

Vibrational polaritons, which have been primarily studied in Fabry–Pérot cavities with a large number of molecules (N ~ 10 6 –10 10 ) coupled to the resonator mode, exhibit various experimentally observed effects on chemical reactions. However, the exact mechanism is elusively understood from the theoretical side, as the large number of molecules involved in an experimental strong coupling condition cannot be represented completely in simulations. This discrepancy between theory and experiment arises from computational descriptions of polariton systems typically being limited to only a few molecules, thus failing to represent the experimental conditions adequately. To address this mismatch, we used surface phonon polariton (SPhP) resonators as an alternative platform for vibrational strong coupling. SPhPs exhibit strong electromagnetic confinement on the surface and thus allow for coupling to a small number of molecules. As a result, this platform can enhance nonlinearity and slow down relaxation to the dark modes. In this study, we fabricated a pillar-shaped quartz resonator and then coated it with a thin layer of cobalt phthalocyanine (CoPc). By employing scattering-type scanning near-field optical microscopy (s-SNOM), we spatially investigated the dependency of vibrational strong coupling on the spatially varying electromagnetic field strength and demonstrated strong coupling with 38,000 molecules only–reaching to the small N limit. Through s-SNOM analysis, we found that strong coupling was observed primarily on the edge of the quartz pillar and the apex of the s-SNOM tip, where the maximum field enhancement occurs. In contrast, a weak resonance signal and lack of coupling were observed closer to the center of the pillar. This work demonstrates the importance of spatially resolved polariton systems in nanophotonic platforms and lays a foundation to explore polariton chemistry and chemical dynamics at the small N limit–one step closer to reconcile with high-level quantum calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗