Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Accelerator Simulation”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Combining Iteration-Free Polarization with Large Time Step Stochastic-Isokinetic Integration

In order to accelerate molecular dynamics simulations using polarizable force fields, we combine a new extended Lagrangian approach that eliminates the self-consistent field step (iEL/0-SCF) with a stochastic integration scheme that allows for a long time step using a multiple time scale algorithm (SIN(R)). We consider different algorithms for the combined scheme that places different components of the nonbonded forces into different time scales, as well as splitting individual nonbonded forces across time scales, to demonstrate that the combined method works well for bulk water as well as for a concentrated salt solution, aqueous peptide, and solvated protein. Depending on system and desired accuracy, the iEL/0-SCF and SIN(R) combination yields lower bound computational speed-ups of ~6–8 relative to a molecular dynamics Verlet integration using a standard SCF solver implemented in the reference program TINKER 8.1. Finally, the combined approach embodies a significant advance for equilibrium simulations in the canonical ensemble of many-body potential energy surfaces for condensed phase systems with speed-ups that exceed what is possible by either method alone.

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

Nitromethane Decomposition via Automated Reaction Discovery and an Ab Initio Corrected Kinetic Model

In the explore the systematic construction of kinetic models from in silico reaction data for the decomposition of nitromethane. Our models are constructed in a computationally affordable manner by using reactions discovered through accelerated molecular dynamics simulations using the ReaxFF reactive force field. The reaction paths are then optimized to determine reaction rate parameters. We introduce a reaction barrier correction scheme that combines accurate thermochemical data from density functional theory with ReaxFF minimal energy paths. We validate our models across different thermodynamic regimes, showing predictions of gas phase CO and NO concentrations and high-pressure induction times that are similar to experimental data. The kinetic models are analyzed to find fundamental decomposition reactions in different thermodynamic regimes.

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

Evidence for percolation diffusion of cations and reordering in disordered pyrochlore from accelerated molecular dynamics

Diffusion in complex oxides is critical to ionic transport, radiation damage evolution, sintering, and aging. In complex oxides such as pyrochlores, anionic diffusion is dramatically affected by cation disorder. However, little is known about how disorder influences cation transport. Here, we report results from classical and accelerated molecular dynamics simulations of vacancy-mediated cation diffusion in Gd 2 Ti 2 O 7 pyrochlore, on the microsecond timescale. We find that diffusion is slow at low levels of disorder, while higher disorder allows for fast diffusion, which is then accompanied by antisite annihilation and reordering, and thus a slowing of cation transport. Cation diffusivity is therefore not constant, but decreases as the material reorders. We also show that fast cation diffusion is triggered by the formation of a percolation network of antisites. This is in contrast with observations from other complex oxides and disordered media models, suggesting a fundamentally different relation between disorder and mass transport.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

The role of charge in microdroplet redox chemistry

In charged water microdroplets, which occur in nature or in the lab upon ultrasonication or in electrospray processes, the thermodynamics for reactive chemistry can be dramatically altered relative to the bulk phase. Here, we provide a theoretical basis for the observation of accelerated chemistry by simulating water droplets of increasing charge imbalance to create redox agents such as hydroxyl and hydrogen radicals and solvated electrons. We compute the hydration enthalpy of OH - and H + that controls the electron transfer process, and the corresponding changes in vertical ionization energy and vertical electron affinity of the ions, to create OH • and H • reactive species. We find that at ~ 20 - 50% of the Rayleigh limit of droplet charge the hydration enthalpy of both OH - and H + have decreased by >50 kcal/mol such that electron transfer becomes thermodynamically favorable, in correspondence with the more favorable vertical electron affinity of H + and the lowered vertical ionization energy of OH - . We provide scaling arguments that show that the nanoscale calculations and conclusions extend to the experimental microdroplet length scale. The relevance of the droplet charge for chemical reactivity is illustrated for the formation of H 2 O 2 , and has clear implications for other redox reactions observed to occur with enhanced rates in microdroplets.

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

Laboratory realization of relativistic pair-plasma beams

Relativistic electron-positron plasmas are ubiquitous in extreme astrophysical environments such as black-hole and neutron-star magnetospheres, where accretion-powered jets and pulsar winds are expected to be enriched with electron-positron pairs. Their role in the dynamics of such environments is in many cases believed to be fundamental, but their behavior differs significantly from typical electron-ion plasmas due to the matter-antimatter symmetry of the charged components. So far, our experimental inability to produce large yields of positrons in quasi-neutral beams has restricted the understanding of electron-positron pair plasmas to simple numerical and analytical studies, which are rather limited. We present the first experimental results confirming the generation of high-density, quasi-neutral, relativistic electron-positron pair beams using the 440 GeV/c beam at CERN’s Super Proton Synchrotron (SPS) accelerator. Monte Carlo simulations agree well with the experimental data and show that the characteristic scales necessary for collective plasma behavior, such as the Debye length and the collisionless skin depth, are exceeded by the measured size of the produced pair beams. Our work opens up the possibility of directly probing the microphysics of pair plasmas beyond quasi-linear evolution into regimes that are challenging to simulate or measure via astronomical observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Near-Field Corrosion Interactions between Glass and Corrosion Resistant Alloys

This study explores the corrosion interactions between model nuclear waste glass materials and corrosion resistant alloys, under conditions that simulate the near field of a nuclear waste repository. The interactions between the corrosion of stainless steel (SS) 316, alloy G30, or alloy 625, and international simple glass or soda-lime silica glass are systematically studied. The dissimilar materials were exposed in close proximity to each other in different electrolytes at 90 oC. After exposure, the glass surface exposed near metals showed different regimes of corrosion, with distinct surface morphologies and chemistries that were likely affected by the local environment created by the localized corrosion of metals. Surface and solution analyses showed that the corrosion rate of glass was enhanced by the presence of metals. Infrared spectroscopy data suggested the local build-up of stresses in the contact area of glass, which may lead to the mechanical instability of the glass alteration layer. On the other hand, the effect of glass on metal corrosion is strongly dependent on the leaching solution. In electrolytes containing abundant aggressive anions such as Cl-, glass seems to suppress the localized corrosion of SS by the precipitation of a Si-rich surface film that protects the SS substrate from solutions. However, in less aggressive electrolytes, the corrosion rate of SS was increased by the presence of glass corrosion products. Overall, our study showed that the hidden and localized damage on glass in contact with metals may enhance the release of incorporated radionuclides compared to typical uniform glass corrosion. Therefore, more attention should be paid to the corrosion interactions between metals and glass under repository conditions.

Guo, Xiaolei↗
DOE OSTI2017

The mobility of small vacancy/helium complexes in tungsten and its impact on retention in fusion-relevant conditions

Tungsten is a promising plasma facing material for fusion reactors. Despite many favorable properties, helium ions incoming from the plasma are known to dramatically affect the microstructure of tungsten, leading to bubble growth, blistering, and/or to the formation of fuzz. In order to develop mitigation strategies, it is essential to understand the atomistic processes that lead to bubble formation and subsequent microstructural changes. In this work, we use large-scale Accelerated Molecular Dynamics simulations to investigate small (N=1,2) V N He M vacancy/helium complexes, which serve as the nuclei for larger helium bubble growth, over timescales reaching into the milliseconds under conditions typical of the operation of fusion reactors. These complexes can interconvert between different I L V N+L He M variants via Frenkel pair nucleation (leading to the creation of a additional vacancy/interstitial pair) and annihilation events; sequences of these events can lead to net migration of these embryonic bubbles. The competition between nucleation and annihilation produces a very complex dependence of the diffusivity on the number of heliums. Finally, through cluster dynamics simulations, we show that diffusion of these complexes provides an efficient pathway for helium release at fluxes expected in fusion reactors, and hence that accounting for the mobility of these complexes is crucial.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Insights into dynamic processes of cations in pyrochlores and other complex oxides

Complex oxides are critical components of many key technologies, from solid oxide fuel cells and superionics to inert matrix fuels and nuclear waste forms. In many cases, understanding mass transport is important for predicting performance and, thus, extensive effort has been devoted to understanding mass transport in these materials. However, most work has focused on the behavior of oxygen while cation transport has received relatively little attention, even though cation diffusion is responsible for many phenomena, including sintering, radiation damage evolution, and deformation processes. Here, we use accelerated molecular dynamics simulations to examine the kinetics of cation defects in one class of complex oxides, A₂B₂O₇ pyrochlore. In some pyrochlore chemistries, B cation defects are kinetically unstable, transforming to A cation defects and antisites at rates faster than they can diffuse. When this occurs, transport of B cations occurs through defect processes on the A sublattice. Further, these A cation defects, either interstitials or vacancies, can interact with antisite disorder, reordering the material locally, though this process is much more efficient for interstitials than vacancies. Whether this behavior occurs in a given pyrochlore depends on the A and B chemistry. Pyrochlores with a smaller ratio of cation radii exhibit this complex behavior, while those with larger ratios exhibit direct migration of B interstitials. Similar behavior has been reported in other complex oxides such as spinels and perovskites, suggesting that this coupling of transport between the A and B cation sublattices, while not universal, occurs in many complex oxide.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Pair mobility functions for rigid spheres in concentrated colloidal dispersions: Force, torque, translation, and rotation

The formulation of detailed models for the dynamics of condensed soft matter including colloidal suspensions and other complex fluids requires accurate description of the physical forces between microstructural constituents. In dilute suspensions, pair-level interactions are sufficient to capture hydrodynamic, interparticle, and thermodynamic forces. In dense suspensions, many-body interactions must be considered. Prior analytical approaches to capturing such interactions such as mean-field approaches replace detailed interactions with averaged approximations. However, long-range coupling and effects of concentration on local structure, which may play an important role in, e.g., phase transitions, are smeared out in such approaches. An alternative to such approximations is the detailed modeling of hydrodynamic interactions utilizing precise couplings between moments of the hydrodynamic traction on a suspended particle and the motion of that or other suspended particles. For two isolated spheres, a set of these functions was calculated by Jeffrey and Onishi [J. Fluid Mech. 139, 261–290 (1984)] and Jeffrey [J. Phys. Fluids 4, 16–29 (1992)]. Along with pioneering work by Batchelor, these are the touchstone for low-Reynolds-number hydrodynamic interactions and have been applied directly in the solution of many important problems related to the dynamics of dilute colloidal dispersions [G. K. Batchelor and J. T. Green, J. Fluid Mech. 56, 375–400 (1972) and G. K. Batchelor, J. Fluid Mech. 74, 1–29 (1976)]. Toward extension of these functions to concentrated systems, here we present a new stochastic sampling technique to rapidly calculate an analogous set of mobility functions describing the hydrodynamic interactions between two hard spheres immersed in a suspension of arbitrary concentration, utilizing accelerated Stokesian dynamics simulations. These mobility functions provide precise, radially dependent couplings of hydrodynamic force and torque to particle translation and rotation, for arbitrary colloid volume fraction ϕ. The pair mobilities (describing entrainment of one particle by the disturbance flow created by another) decay slowly with separation distance: as 1/r, for volume fractions 0.05 ≤ ϕ ≤ 0.5. For the relative mobility, we find an initially rapid growth as a pair separates, followed by a slow, 1/r growth. Up to ϕ ≤ 0.4, the relative mobility does not reached the far-field value even beyond separations of many particle sizes. In the case of ϕ = 0.5, the far-field asymptote is reached but only at a separation of eight radii and after a slow 1/r growth. At these higher concentrations, the coefficients also reveal liquid-like structural effects on pair mobility at close separations. These results confirm that long-range many-body hydrodynamic interactions are an essential part of the dynamics of concentrated systems and that care must be taken when applying renormalization schemes.

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

The purported square ice in bilayer graphene is a nanoscale, monolayer object

The phase diagram of water is complex, and interfacial effects can stabilize unusual structures at the nanoscale. In this paper, we employ bond order accelerated molecular dynamics simulations to show that upon encapsulation within bilayer graphene, water can spontaneously adopt a two-dimensional (monomolecular) layer of “square ice” at ambient conditions, instead of an encapsulated water droplet. Free energy calculations show that this motif is thermodynamically stable up to diameters of approximately 15 nm due to enhanced hydrogen bonding and favorable binding to the graphene sheets. Entropic losses due to solidification and reduced graphene–graphene binding enthalpy are opposing thermodynamic forces that conspire to limit the maximum size, but modification of any of these thermodynamic factors should change the range of stability. Finally, simulated core-level spectroscopy reveals unambiguous orientation dependent signatures of square ice that should be discernable in experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Plateau–Rayleigh instability with a grain boundary twist

We demonstrate using theoretical, computational, and experimental studies a morphological instability, in which a polycrystalline nanorod breaks up at grain boundaries (GBs) into an array of isolated domains. Our theoretical model is used to establish a neutral stability surface demarcating stable and unstable perturbations. It is shown that GBs play a destabilizing role in which the critical wavelength for the instability decreases with the increase in the GB energy. We carry out phase field simulations, which reveal accelerated pinch-off kinetics with the increase in the GB energy and predict temporal evolution of interfacial profiles in quantitative agreement with experimental observations.

Hussein, Omar↗
DOE OSTI2024

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Mechanistic studies of small molecule ligands selective to RNA single G bulges

Abstract Small-molecule RNA binders have emerged as an important pharmacological modality. A profound understanding of the ligand selectivity, binding mode, and influential factors governing ligand engagement with RNA targets is the foundation for rational ligand design. Here, we report a novel class of coumarin derivatives exhibiting selective binding affinity towards single G RNA bulges. Harnessing the computational power of all-atom Gaussian accelerated molecular dynamics simulations, we unveiled a rare minor groove binding mode of the ligand with a key interaction between the coumarin moiety and the G bulge. This predicted binding mode is consistent with results obtained from structure-activity relationship studies and transverse relaxation measurements by nuclear magnetic resonance spectroscopy. We further generated 444 molecular descriptors from 69 coumarin derivatives and identified key contributors to the binding events, such as charge state and planarity, by lasso (least absolute shrinkage and selection operator) regression. Our work deepened the understanding of RNA-small molecule interactions and integrated a new framework for the rational design of selective small-molecule RNA binders.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Extracting symplectic maps for space-charge dominated beams

Symplecticity of transfer maps is important for reliable evaluation of space-charge dominated beams in accelerators. Unfortunately, most simulation codes that include collective effects, such as space charge, do not use canonical phase-space variables and therefore are not symplectic in the presence of electromagnetic fields. In this paper, we present a numerical method to extract local linear symplectic transfer maps using particle tracking simulation code for space-charge dominated beams. We demonstrate this method for the photoinjector (113 MHz SRF gun) section of the Coherent electron Cooling (CeC) Proof of Principle (POP) experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Ion acceleration from microstructured targets irradiated by high-intensity picosecond laser pulses

Structures on the front surface of thin foil targets for laser-driven ion acceleration have been proposed to increase the ion source maximum energy and conversion efficiency. While structures have been shown to significantly boost the proton acceleration from pulses of moderate-energy fluence, their performance on tightly focused and high-energy lasers remains unclear. Here, we report the results of laser-driven three-dimensional (3D)-printed microtube targets, focusing on their efficacy for ion acceleration. Using the high-contrast (~10 12) PHELIX laser (150 J, 10 21 W / cm 2 ), we studied the acceleration of ions from 1-μm-thick foils covered with micropillars or microtubes, which we compared with flat foils. The front-surface structures significantly increased the conversion efficiency from laser to light ions, with up to a factor of 5 higher proton number with respect to a flat target, albeit without an increase of the cutoff energy. An optimum diameter was found for the microtube targets. Our findings in this work are supported by a systematic particle-in-cell modeling investigation of ion acceleration using 2D simulations with various structure dimensions. Simulations reproduce the experimental data with good agreement, including the observation of the optimum tube diameter, and reveal that the laser is shuttered by the plasma filling the tubes, explaining why the ion cutoff energy was not increased in this regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Exascale workflow applications and middleware: An ExaWorks retrospective

Exascale computers offer transformative capabilities to combine data-driven and learning-based approaches with traditional simulation applications to accelerate scientific discovery and insight. However, these software combinations and integrations are difficult to achieve due to the challenges of coordinating and deploying heterogeneous software components on diverse and massive platforms. Here, we present the ExaWorks project, which addresses many of these challenges. We developed a workflow Software Development Toolkit (SDK), a curated collection of workflow technologies that can be composed and interoperated through a common interface, engineered following current best practices, and specifically designed to work on HPC platforms. ExaWorks also developed PSI/J, a job management abstraction API, to simplify the construction of portable software components and applications that can be used over various HPC schedulers. The PSI/J API is a minimal interface for submitting and monitoring jobs and their execution state across multiple and commonly used HPC schedulers. We also describe several leading and innovative workflow examples of ExaWorks tools used on DOE leadership platforms. Furthermore, we discuss how our project is working with the workflow community, large computing facilities, and HPC platform vendors to address the requirements of workflows sustainably at the exascale.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Virtual Observation System for Earth System Model: An Application to ACME Land Model Simulations

Investigating and evaluating physical-chemical-biological processes within an Earth system model (EMS) can be very challenging due to the complexity of both model design and software implementation. A virtual observation system (VOS) is presented to enable interactive observation of these processes during system simulation. Based on advance computing technologies, such as compiler-based software analysis, automatic code instrumentation, and high-performance data transport, the VOS provides run-time observation capability, in-situ data analytics for Earth system model simulation, model behavior adjustment opportunities through simulation steering. A VOS for a terrestrial land model simulation within the Accelerated Climate Modeling for Energy model is also presented to demonstrate the implementation details and system innovations.

58 GEOSCIENCES↗