Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulator study”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Activation of atomic transport via vibrational coupling-induced force fluctuations

Vibrational coupling, although well-explored in many fields, has seldom been investigated in the context of mass transport. In this Letter, we examine the impact of vibrational coupling on atomic transport using simple molecular dynamics simulations. Our study shows that the atomic transport can be activated when the natural frequency of the atomic slit is close to the natural frequency of the atom being transported. We uncover the presence of fluctuating forces induced by vibrational coupling with higher amplitudes observed when the coupling is strong. We show that the transport activation mechanism is due to the high force fluctuations that arise during strong vibrational coupling allowing the atom to temporarily surpass the transport barrier of the slit. Our findings will serve as a foundation for the continued examination of vibrational coupling in the realm of mass transport.

Physics↗
DOE OSTI2023

Critical behavior in a chiral molecular model

Understanding the condensed-phase behavior of chiral molecules is important in biology as well as in a range of technological applications, such as the manufacture of pharmaceuticals. Here, we use molecular dynamics simulations to study a chiral four-site molecular model that exhibits a second-order symmetry-breaking phase transition from a supercritical racemic liquid into subcritical D-rich and L-rich liquids. We determine the infinite-size critical temperature using the fourth-order Binder cumulant, and we show that the finite-size scaling behavior of the order parameter is compatible with the 3D Ising universality class. We also study the spontaneous D-rich to L-rich transition at a slightly subcritical temperature of T = 0.985Tc, and our findings indicate that the free energy barrier for this transformation increases with system size as N2/3, where N is the number of molecules, consistent with a surface-dominated phenomenon. The critical behavior observed herein suggests a mechanism for chirality selection in which a liquid of chiral molecules spontaneously forms a phase enriched in one of the two enantiomers as the temperature is lowered below the critical point. Furthermore, the increasing free energy barrier with system size indicates that fluctuations between the L-rich and D-rich phases are suppressed as the size of the system increases, trapping it in one of the two enantiomerically enriched phases. Such a process could provide the basis for an alternative explanation for the origin of biological homochirality. We also conjecture the possibility of observing nucleation at subcritical temperatures under the action of a suitable chiral external field.

Chemistry↗
DOE OSTI2023

Heat transport with a twist

Despite the desirability of polymers for use in many products due to their flexibility, light weight, and durability, their status as thermal insulators has precluded their use in applications where thermal conductors are required. However, recent results suggest that the thermal conductance of polymers can be enhanced and that their heat transport behaviors may be highly sensitive to nanoscale control. Here we use non-equilibrium molecular dynamics simulations to study the effect of mechanical twist on the steady-state thermal conductance across multi-stranded polyethylene wires. We find that a highly twisted double-helical polyethylene wire can display a thermal conductance up to three times that of its untwisted form, an effect which can be attributed to a structural transition in the strands of the double helix. We also find that in thicker wires composed of many parallel strands, adding just one twist can increase its thermal conductance by over 30%. However, we find that unlike stretching a polymer wire, which causes a monotonic increase in thermal conductance, the effect of twist is highly non-monotonic, and certain amounts of twist can actually decrease the thermal conductance. Finally, we apply the Continuous Chirality Measure (CCM) in an attempt to explore the correlation between heat conductance and chirality. The CCM is found to correlate with twist as expected, but we attribute the observed heat transport behaviors to structural factors other than chirality.

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

Kohn–Sham accuracy from orbital-free density functional theory via Δ-machine learning

Here, we present a Δ-machine learning model for obtaining Kohn–Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine-learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn–Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations and study its accuracy, performance, and sensitivity to parameters for representative systems. We find that the formalism not only improves the accuracy of Thomas–Fermi–von Weizsäcker orbital-free energies and forces by more than two orders of magnitude but is also more accurate than MLFFs based solely on Kohn–Sham DFT while being more efficient and less sensitive to model parameters. We apply the framework to study the structure of molten Al 0.88 Si 0.12 , the results suggesting no aggregation of Si atoms, in agreement with a previous Kohn–Sham study performed at an order of magnitude smaller length and time scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Mining experimental magnetized liner inertial fusion data: Trends in stagnation morphology

In magnetized liner inertial fusion (MagLIF), a cylindrical liner filled with fusion fuel is imploded with the goal of producing a one-dimensional plasma column at thermonuclear conditions. However, structures attributed to three-dimensional effects are observed in self-emission x-ray images. Despite this, the impact of many experimental inputs on the column morphology has not been characterized. We demonstrate the use of a linear regression analysis to explore correlations between morphology and a wide variety of experimental inputs across 57 MagLIF experiments. Results indicate the possibility of several unexplored effects. For example, we demonstrate that increasing the initial magnetic field correlates with improved stability. Although intuitively expected, this has never been quantitatively assessed in integrated MagLIF experiments. We also demonstrate that azimuthal drive asymmetries resulting from the geometry of the “current return can” appear to measurably impact the morphology. In conjunction with several counterintuitive null results, we expect the observed correlations will encourage further experimental, theoretical, and simulation-based studies. Finally, we note that the method used in this work is general and may be applied to explore not only correlations between input conditions and morphology but also with other experimentally measured quantities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Invariant discovery of features across multiple length scales: Applications in microscopy and autonomous materials characterization

Physical imaging is a foundational characterization method in areas from condensed matter physics and chemistry to astronomy and spans length scales from atomic to universe. Images encapsulate crucial data regarding atomic bonding, materials microstructures, and dynamic phenomena such as microstructural evolution and turbulence, among other phenomena. The challenge lies in effectively extracting and interpreting this information. Variational Autoencoders (VAEs) have emerged as powerful tools for identifying the underlying factors of variation in image data, providing a systematic approach to distilling meaningful patterns from complex data sets. However, a significant hurdle in their application is the definition and selection of appropriate descriptors reflecting local structures. Here, we introduce the scale-invariant VAE approach (SI-VAE) based on the progressive training of the VAE with the descriptors sampled at different length scales. The SI-VAE allows the discovery of the length scale-dependent factors of variation in the system. Here, we illustrate this approach using the ferroelectric domain images and generalize it to the movies of the electron-beam induced phenomena in graphene and topography evolution across combinatorial libraries. This approach can further be used to initialize the decision making in automated experiments including structure–property discovery and can be applied across a broad range of imaging methods. This approach is universal and can be applied to any spatially resolved data including both experimental imaging studies and simulations, and can be particularly useful for exploration of phenomena such as turbulence and scale-invariant transformation fronts.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Irradiation-induced gas production in REBCO-based magnet materials used for future compact fusion reactors

Nuclear fusion is an enticing alternative to current sources of energy, with multilayered Rare-Earth Barium Copper Oxide (REBCO) coated conductors deemed pivotal in the race toward fully realized, commercially viable, and magnetic confinement fusion reactors. In this study, we simulated the ion spectrum expected to evolve from REBCO's nickel-based Hastelloy C-276 substrate and copper stabilizer in an affordable robust compact-like reactor. We then emulated this gas production through helium implantation to investigate changes in materials and superconducting properties. Our results revealed that the substrate and stabilizer are capable of producing protons energetic enough to recoil throughout the tape thickness in appreciable doses, and alphas energetic enough to deposit 7.54 × 1014 ions/cm2 or 50.1 helium appm in the superconducting layer over a 30-year reactor lifetime. The superconducting layer of SuperPower® tapes exhibited at least double the swelling rate of the other major layers, and both SuperPower and Fujikura Ltd. tapes displayed microstructural changes in the REBCO layer not observed in isotropic metals. For the estimated lifetime fluence, the Fujikura tapes showed a ∼1 K reduction in critical temperature and a 32% degradation in critical current for compact reactor-relevant conditions (16 T, 20 K). Nuclear transmutation, low-temperature solder implantations, gas-ion evolution, the influence of gas production on vortex dynamics, and other related considerations are also discussed.

Reis, Chris↗
DOE OSTI2025

Single-shot in-line x-ray phase-contrast imaging of void-shockwave interactions in fusion energy materials

Recent breakthroughs in nuclear fusion, specifically the report of reactions exceeding scientific breakeven at the National Ignition Facility (NIF), highlight the potential of inertial fusion energy (IFE) as a sustainable and virtually limitless energy source. However, further progress in IFE requires characterization of defects in ablator materials and how they affect fuel capsule compression. Voids within the ablator can degrade energy yield, but their impact on the density distribution has primarily been studied through simulations, with limited high-resolution experimental validation. To address this, we used the x-ray free-electron laser (XFEL) at the matter in extreme conditions (MECs) instrument at the Linac coherent light source (LCLS) to capture 2D x-ray phase-contrast (XPC) images of a void-bearing sample with a composition similar to inertial confinement fusion (ICF) ablators. By driving a compressive shockwave through the sample using MEC's long-pulse laser system, we analyzed how voids influence shockwave propagation and density distribution during compression. To quantify this impact, we extracted phase information using two phase retrieval algorithms. First, we applied the contrast transfer function (CTF) method, paired with Tikhonov regularization and a fast optimization approach to generate an initial phase estimate. We then refined the result using a projected gradient descent (PGD) method that works directly with the sample's refractive index. Comparing these results with radiation adaptive grid Eulerian (xRAGE) radiation hydrodynamic simulations enables identification of model validation needs or improvements. By calculating phase maps in situ, it becomes possible to reconstruct areal density maps, improving understanding of laser-capsule interactions and advancing IFE research.

Hodge, D. S. [Colorado State Univ., Fort Collins, ↗
DOE OSTI2025

THz cyclotron resonance of a 2D hole gas in a GaN/AlN heterostructure

The recent discovery of highly conducting two-dimensional hole gases (2DHGs) in GaN/AlN heterojunctions has opened the door to efficient complementary GaN electronics, a long-standing challenge in wide-bandgap semiconductor device physics. Electrical transport studies and simulations indicate that both heavy- and light-hole valence bands are occupied in these 2DHGs, but direct experimental characterization of the fundamental parameters of the mobile holes remains at an early stage. Here, we use time-domain terahertz spectroscopy and pulsed magnetic fields up to 31 T to directly measure cyclotron resonance of the mobile 2D holes in these GaN-based 2DHGs at low temperature (8 K), revealing key material properties including effective masses, densities, scattering times, and mobilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Population inversion X-ray laser oscillator

Oscillators are at the heart of optical lasers, providing stable, transform-limited pulses. Until now, laser oscillators have been available only in the infrared to visible and near-ultraviolet (UV) spectral region. In this paper, we present a study of an oscillator operating in the 5- to 12-keV photon-energy range. We show that, using the Kα 1 line of transition metal compounds as the gain medium, an X-ray free-electron laser as a periodic pump, and a Bragg crystal optical cavity, it is possible to build X-ray oscillators producing intense, fully coherent, transform-limited pulses. As an example, we consider the case of a copper nitrate gain medium generating ~ 5 × 10 10 photons per pulse with 37-fs pulse length and 48-meV spectral resolution at 8-keV photon energy. Our theoretical study and simulation of this system show that, because of the very large gain per pass, the oscillator saturates and reaches full coherence in four to six optical-cavity transits. We discuss the feasibility and design of the X-ray optical cavity and other parts of the oscillator needed for its realization, opening the way to extend X-ray–based research beyond current capabilities.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Early-stage dynamics of chloride ion–pumping rhodopsin revealed by a femtosecond X-ray laser

Chloride ion–pumping rhodopsin (ClR) in some marine bacteria utilizes light energy to actively transport Cl – into cells. How the ClR initiates the transport is elusive. Here, we show the dynamics of ion transport observed with time-resolved serial femtosecond (fs) crystallography using the Linac Coherent Light Source. X-ray pulses captured structural changes in ClR upon flash illumination with a 550 nm fs-pumping laser. High-resolution structures for five time points (dark to 100 ps after flashing) reveal complex and coordinated dynamics comprising retinal isomerization, water molecule rearrangement, and conformational changes of various residues. Combining data from time-resolved spectroscopy experiments and molecular dynamics simulations, this study reveals that the chloride ion close to the Schiff base undergoes a dissociation–diffusion process upon light-triggered retinal isomerization.

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

Elastic forces drive nonequilibrium pattern formation in a model of nanocrystal ion exchange

The widely used process of nanocrystal ion exchange operates out of thermodynamic equilibrium and can require mixing components of varying sizes. Here we use theory and computer simulation to study a simple model which captures these two basic features of ion exchange reactions. We show that a strong driving force for exchange among different-sized species creates nonequilibrium patterns within model nanocrystals. We further demonstrate that such patterns can be thermodynamically stable within core/shell nanocrystals. These results help us understand the heterostructures formed in ion-exchanged nanocrystals and suggest strategies for leveraging elasticity to design patterns in nanoscale materials. Chemical transformations, such as ion exchange, are commonly employed to modify nanocrystal compositions. Yet the mechanisms of these transformations, which often operate far from equilibrium and entail mixing diverse chemical species, remain poorly understood. Here we explore an idealized model for ion exchange in which a chemical potential drives compositional defects to accumulate at a crystal’s surface. These impurities subsequently diffuse inward. We find that the nature of interactions between sites in a compositionally impure crystal strongly impacts exchange trajectories. In particular, elastic deformations which accompany lattice-mismatched species promote spatially modulated patterns in the composition. These same patterns can be produced at equilibrium in core/shell nanocrystals, whose structure mimics transient motifs observed in nonequilibrium trajectories. Moreover, the core of such nanocrystals undergoes a phase transition—from modulated to unstructured—as the thickness or stiffness of the shell is decreased. Our results help explain the varied patterns observed in heterostructured nanocrystals produced by ion exchange and suggest principles for the rational design of compositionally patterned nanomaterials.

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

Influence of nuclear quantum effects on the electronic properties of amorphous carbon

We carry out quantum simulations to study the physical properties of diamond-like amorphous carbon by coupling first-principles molecular dynamics with a quantum thermostat, and we analyze multiple samples representative of different defective sites present in the disordered network. We show that quantum vibronic coupling is critical in determining the electronic properties of the system, in particular its electronic and mobility gaps, while it has a moderate influence on the structural properties. We find that despite localized electronic states near the Fermi level, the quantum nature of the nuclear motion leads to a renormalization of the electronic gap surprisingly similar to that found in crystalline diamond. We also discuss the notable influence of nuclear quantum effects on band-like and variable-hopping mechanisms contributing to electrical conduction. Our calculations indicate that methods often used to evaluate electron–phonon coupling in ordered solids are inaccurate to study the electronic and transport properties of amorphous semiconductors composed of light atoms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Arrested coarsening and large density fluctuations in driven particle mixtures in two dimensions

Abstract Using molecular dynamics simulations, we study a driven, nonadditive binary mixture of spherical particles confined to move in two dimensions and immersed in an explicit solvent consisting of point particles with purely repulsive interactions. We show that, without a drive, the mixture of spherical particles phase separates and coarsens with kinetics consistent with an Ising-like conserved dynamics. Conversely, when the drive is applied, the coarsening is arrested and the system develops large density fluctuations. We show that the drive creates domains of a characteristic size which decreases with an increasing force. Furthermore, we find that these domains are anisotropic and can be oriented either parallel or perpendicular to the drive direction. Finally, we connect our findings to existing theories of strongly-driven systems, pointing out the importance of introducing the explicit solvent particles to break the Galilean invariance of the system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Progress in DIII-D towards validating divertor power exhaust predictions

In this study, UEDGE simulations highlight the role of cross-field drifts on the onset of detached conditions, and new calibrated divertor vacuum ultra violet (VUV) spectroscopy is used to challenge the predictions of radiative constituents in these simulations. UEDGE simulations for DIII-D H-mode plasmas with the open divertor with the ion ∇B-drift towards the X-point show a bifurcated onset of the low field side (LFS) divertor detachment, consistent with experimentally observed step-like detachment onset (Jaervinen A.E. et al 2018 Phys. Rev. Lett. 121 075001). The divertor plasma in the simulations exhibits hysteresis in upstream separatrix density between attached and detached solution branches. Reducing the drift magnitude by a factor of 3 eliminates the step-like detachment onset in the simulations, confirming the strong role of drifts in the bifurcated detachment onset. When measured local plasma densities and temperatures are within proximity of predicted values in the simulations, there is no shortfall of the local emission of the dominant resonant radiating lines. However, the simulations systematically predict a factor of two lower total integrated radiated power than measured by the bolometer with the difference lost through radial heat flow out of the computational domain. Even though there is no shortfall in the emission of the dominant lines, a shortfall of total radiated power can be caused by underpredicted spatial extent of the radiation front, indicating a potential upstream or divertor transport physics origin for the radiation shortfall, or shortfall of radiated power in the spectrum between the dominant lines. In addition to the underpredicted spatial extent, in detached conditions, the simulations overpredict the peak radiation and dominant carbon lines near the X-point, which can be alleviated by manually increasing divertor diffusivity in the simulations, highlighting the ad hoc cross-field transport as one of the key limitations of the predictive capability of these divertor fluid codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Experimental evaluation of the vapor box divertor concept with an open vapor box module in Magnum-PSI

A promising approach to handle the intense plasma heat flux in the divertor region of a tokamak is the Vapor Box Divertor (VBD). Here plasma-lithium interaction creates a dense lithium vapor cloud which interacts with the incoming plasma, effectively shielding the tungsten surface beneath, therefore preventing overheating and sputtering of the tungsten and increasing the component’s lifetime. Two key steps must be addressed to validate this concept: investigating plasma-Li interactions and the transport mechanism of the latter in the presence of plasma. To explore this, a Vapor Box Module (VBM) has been designed for use with the linear plasma generator Magnum-PSI. The VBM consists of a series of three cylindrical boxes, with Li being evaporated at a controlled temperature in the central box (CB). Divertor-like plasma enters the VBM from an upstream aperture, interacts with the Li vapor cloud and exits through a downstream aperture ultimately impacting a target equipped with a calorimetry system. Optical diagnostics Thomson scattering, Filtered fast camera Imaging and optical emission spectroscopy provided information on plasma parameters in terms of electron density n e , temperature T e and plasma composition before and after this interaction. A significant reduction in plasma power was observed upon the establishment of lithium vaporization determined via cooling water calorimetry and embedded thermocouples at the downstream target. This resulted in a drop of the target temperature from ∼800 °C to ∼350 °C (57%) at the highest applied power (11.1 MW m −2 ) and CB temperature (∼700 °C). Lithium condensation at the side boxes of the VBM in combination with strong plasma momentum transfer effectively prevented upstream migration of lithium, while enhancing Li transport toward the target. The achievement of the two main goals of reducing plasma power and confining the Li in the VBM, are consistent with earlier published preliminary SOLPS-ITER simulations. This study shows that the presence of lithium in a VBD-like configuration can result in a strong reduction of power to the target surface, while at the same time the lithium vapor is effectively confined by the VBM, aided by the incoming plasma, preventing its escape from the VBM geometry. This represents a valuable step toward validating the feasibility of the VBD configuration in future fusion reactors.

Magnum-PSI↗
DOE OSTI2025

The track-length extension fitting algorithm for energy measurement of interacting particles in liquid argon TPCs and its performance with ProtoDUNE-SP data

This paper introduces a novel track-length extension fitting algorithm for measuring the kinetic energies of inelastically interacting particles in liquid argon time projection chambers (LArTPCs). The algorithm finds the most probable offset in track length for a track-like object by comparing the measured ionization density as a function of position with a theoretical prediction of the energy loss as a function of the energy, including models of electron recombination and detector response. The algorithm can be used to measure the energies of particles that interact before they stop, such as charged pions that are absorbed by argon nuclei. The algorithm's energy measurement resolutions and fractional biases are presented as functions of particle kinetic energy and number of track hits using samples of stopping secondary charged pions in data collected by the ProtoDUNE-SP detector, and also in a detailed simulation. Additional studies describe the impact of the dE/dx model on energy measurement performance. The method described in this paper to characterize the energy measurement performance can be repeated in any LArTPC experiment using stopping secondary charged pions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗