SEARCH · Search NASA

Results for “fast time 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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Modeling and Simulation of Xe-100-type Pebble Bed Gas-Cooled Reactor with SCALE

The US Department of Energy (DOE) announced the Advanced Reactor Demonstration Program (ARDP) to accelerate the deployment of advanced reactor concepts. Awardees of ARDP funds are expected to demonstrate the operation of an advanced reactor within 7 years of receiving the award. X-Energy’s advanced reactor concept, the Xe-100, was selected as one of two advanced reactor concepts to receive funding to demonstrate the operation of its high-temperature gas-cooled pebble-bed reactor before the end of this decade. As a result of this push to bring advanced reactors to maturation and commercialization, transition and deployment scenario studies are being performed under the Systems Analysis and Integration (SA&I) campaign within the DOE Office of Nuclear Energy (DOE-NE) to evaluate the transition of the current US commercial fleet of light-water reactors (LWRs) to a future fleet of advanced reactors consisting of a mix of ARDP type reactor concepts and advanced LWRs. To accurately evaluate the front- and back-end resource requirements, it is important to perform reactor physics calculations to determine the discharge burnup and isotopic content, fuel residence time, as well as other parameters. For this purpose, a summer project funded by the SA&I campaign allowed for the setup of SCALE models for full-core Xe-100 type high-temperature gas-cooled pebble-bed reactor and a Xe-100 type slice using publicly available information. The core-averaged equilibrium compositions and zone-wise equilibrium compositions for the slice and 3D models, respectively, were obtained following an iterative depletion method developed by Bostelmann et al. using SCALE’s reactor physics sequence TRITON. The slice model was used with TRITON to generate burnup-dependent cross section libraries at different temperatures which can be used with SCALE’s ORIGAMI code to rapidly determine fuel inventory and therefore to perform quick sensitivity studies on parameters such as the pebble location in the core. The SCALE/TRITON transport and depletion calculation for the Xe-100 type slice model indicates that the isotopic concentrations are in good agreement at 1,300 effective full power days (EFPD) for 235 U. An analysis of 236 U results match 239 Pu results would seem to indicate a typographical error in Mulder and Boyes wherein the reported results of 236 U and 239 Pu are reversed. In addition to SCALE/TRITON calculations, a new capability within SCALE/ORIGAMI for the simulation of pebble-bed reactors was used to study the burnup sensitivity with respect to the pebble pathway through the core. The SCALE/ORIGAMI results show that pebbles that travel closer to the reflector for the entire depletion history have a higher burnup than pebbles that travel through the middle of the core because of the higher thermal to fast flux ratio near the reflector. Consequently, a pebble’s burnup is strongly affected by the pebble’s pathway for each pass. Additional phenomena such as temperature distributions in the core and different travel times of the pebbles in the individual radial zones further affect the burnup distribution. The sensitivity of the discharge vector to the pebble pathways taken during each pass can be evaluated in the future using SCALE/ORIGAMI now that the SCALE inputs have been established.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Assessing the vulnerability of solar inverters to EMPs: Port testing, PCI modeling, and protection strategies

Renewable energy sources are becoming an ever-larger contributor to the power grid. These renewable energy sources depend upon the power electronic devices, specifically inverters, being essential for connecting Photovoltaic (PV) generation to the grid. However, the Electromagnetic Pulses (EMPs) caused by the high-altitude nuclear explosions can generate fast broad-band pulses with nanosecond rise time, potentially causing damage or destruction to electronic components. To assess the vulnerability of PV inverters to high-altitude EMPs, the port testing and Pulsed Current Injection (PCI) modeling schemes are proposed based on the port impedance analysis. Wide-band frequency measurements are achieved by fusing impedance results from three vector network analyzers. Then, a PCI model is used to simulate the induced response to EMP, with two typical immunity levels of EC5 and EC8 tested. Here, the experiment successfully excites the induced voltage and current under EMP, where the voltage and current can reach 1500V/40A and 8000V/150Aunder EC5 andEC8, respectively. The port vulnerability analysis results demonstrate that only some ports can survive under EC5. To defend against the impact of EMP, three protection strategies are discussed.

42 ENGINEERING↗

DOE OSTI2025

Decoupled few-femtosecond phase transitions in vanadium dioxide

Abstract The nature of the insulator-to-metal phase transition in vanadium dioxide (VO 2 ) is one of the longest-standing problems in condensed-matter physics. Ultrafast spectroscopy has long promised to determine whether the transition is primarily driven by the electronic or structural degree of freedom, but measurements to date have been stymied by their sensitivity to only one of these components and/or their limited temporal resolution. Here we use ultra-broadband few-femtosecond pump-probe spectroscopy to resolve the electronic and structural phase transitions in VO 2 at their fundamental time scales. Our experiments show that the system transforms into a bad-metallic phase within 10 fs after photoexcitation, but requires another 100 fs to complete the transition, during which we observe electronic oscillations and a partial re-opening of the bandgap, signalling a transient semi-metallic state. Comparisons with tensor-network simulations and density-functional theory calculations show these features result from an unexpectedly fast structural transition, in which the vanadium dimers separate and untwist with two different timescales. Our results resolve the structural and electronic nature of the light-induced phase transition in VO 2 and establish ultra-broadband few-femtosecond spectroscopy as a powerful tool for studying quantum materials out of equilibrium.

Science & Technology - Other Topics↗

DOE OSTI2025

Status of Multiple Channel Fuel Performance Capabilities Within the SAS4A/SASSYS-1 Safety Analysis Software

SAS4A/SASSYS-1 (SAS) is a fast-running simulation tool used to perform deterministic analysis of anticipated events as well as design basis and beyond design basis accidents for advanced liquid-metal-cooled nuclear reactors. It is a critical element of safety analysis capabilities for the U.S. Department of Energy and is utilized within industry to perform the transient safety analyses required to support the licensing of Liquid Metal-cooled Fast Reactors (LMFRs). Although SAS is exceptionally fast for most transient scenarios, fuel performance calculations, along with the associated pre-transient characterization of the fuel pin, may be required for transient scenarios where fuel pin failure is hypothesized. Both the pre-transient characterization and the transient fuel performance calculation are necessary to properly quantify margins to potential fuel failure and assess the time spent potentially exceeding such margins during events. While safety analysis calculations with fuel performance models provide a more detailed characterization of the reactor during a transient, the pre-transient characterization can be time-consuming and computationally expensive. Often, large numbers of fuel pins have been exposed to similar pre-transient irradiation conditions. Similarly, the same pre-transient fuel characterization may be applicable to numerous transient conditions. This provides an opportunity to optimize the SAS computational framework such that pre-transient fuel characterization can be shared across multiple channels (fuel pins) and across multiple simulations, thus dramatically reducing overall computational costs. This report summarizes progress toward enhancing the SAS computational framework to support shared, multiple channel fuel performance characterizations intended to significantly reduce computational costs. Preliminary testing has shown that the computational time saved by using the pre-transient sharing capability is approximately equal to the time it takes to perform the pre-transient characterization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Integrating Advanced Modeling and Accelerated Testing for a Modernized Fuel Qualification Paradigm

With the increasing interest in sodium fast reactor technology, as seen by applications to the U.S. Nuclear Regulatory Commission for the OKLO Aurora plant, fuel testing for the TerraPower Traveling Wave Reactor, and the impending construction and startup of the versatile test reactor (VTR), a modernized, accelerated approach to fuel qualification is needed. To guide this effort, a Phenomena Identification Ranking Table–styled analysis was performed for a U-Pu-Zr sodium-free annular fuel system. This analysis evaluated a series of fuel design properties and parameters against their contributions to key fuel performance phenomena. The resulting priority parameters were then reviewed against existing modeling and experimental capabilities to support investigation of the highest-priority parameters. A pathway for qualification was then established using highthroughput, high-volume experiments from MiniFuel and FAST in parallel with advanced physics-based model development. This effort outlines how the first stages of qualification can be reduced from the typical 20+-year development cycle to 5 to 7 years by deploying accelerated irradiation testing platforms. As with any accelerated test, these methods are prototypic in some aspects and less so in others; however, by coupling with advanced fuel performance modeling and simulation capabilities, the larger space of irradiation parameters and material response provided offers advantages for the validation of physics-based models supporting the deployment of novel fuel designs. As a test case, this paper utilizes a proposed Mark II fuel system for the upcoming VTR. Thus, an accelerated qualification method can be tested for the development of MARK II driver fuel so that by the time of VTR startup, lead test assemblies for a Mark II fuel can be initiated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Avoiding Excess Computation in Asynchronous Evolutionary Algorithms

Asynchronous evolutionary algorithms are becoming increasingly popular as a means of making full use of many processors while solving computationally expensive search and optimization problems. These algorithms excel at keeping large clusters fully utilized, but may sometimes inefficiently sample an excess of fast-evaluating solutions at the expense of higher-quality, slow-evaluating ones. We introduce a steady-state parent selection strategy, SWEET (“Selection whilE EvaluaTing”), that sometimes selects individuals that are still being evaluated and allows them to reproduce early. This gives slow-evaluating individuals that have higher fitnesses an increased ability to multiply in the population. We find that SWEET appears effective in simulated take-over time analysis, but that its benefit is confined mostly to early in the run, and our preliminary study on an autonomous vehicle controller problem that involves tuning a spiking neural network proves inconclusive.

Scott, Eric↗

DOE OSTI2022

Redesign of the Coils for the 60T Controlled-Waveform Magnet at NHMFL

Driven by the1.4 GW generator, the 60T Controlled Waveform (60 TCW) magnet was the most powerful controlled waveform system in the world and had always been one of most important magnets to the National High Magnetic Field (NHMFL) and high-field research community because of its following unique features: (1) quasi-static field up to 60 T with 100 ms flat-top and total pulse-length of about 2000 ms, (2) variable magnetic field waveforms such as staircase and triangle with flat-top (3) relative large bore (32 mm) and (4) very fast cooling time (20 minutes) between pulses [Boebinger, 2001], [Crooker et al. 2001]. The magnet is composed of nine concentric coils, with each coil consisting of several conductor winding layers reinforced by a high-strength metallic shell. The magnet underwent a catastrophic failure in 2000 and all the coils had to be rebuilt. In late 2014 the second magnet version failed near the mid-plane of coil 7. The simulations afterward that incident indicated that the overall strength of the coil would be increased by replacing a section of the reinforcing shell with Zylon fiber-epoxy composite. This reduces the stress and thus significantly lowers the level of plastic deformation in the windings. The role of the metal and Zylon fiber reinforcing layers in bearing the axial and radial Lorentz forces has been studied to optimize the magnet design. Here, the results of the optimization will be discussed as well as challenges that have been presented in rebuilding the individual coils of the magnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Neural-Network-Enhanced COTSIM: Advancing Predictive Capabilities for Fast DIII-D Simulations

Sustaining fusion reactions in tokamaks requires heating plasma to thermonuclear temperatures while maintaining confinement and stability. Neutral beam injection (NBI) provides heating, current drive, torque, and fueling, while electron cyclotron (EC) waves are widely used for heating and current drive; together, these actuators shape the plasma current, temperature, and density profiles. The control-oriented tokamak simulator (COTSIM), a predictive, control-oriented code, has been enhanced with neural-network surrogates for transport and sources. Turbulent transport is predicted by MMMnet—a neural-network version of the updated multimode model (MMM 9.0.10)—with significantly reduced computation time relative to MMM; neoclassical transport follows the Chang–Hinton model. NUBEAMnet, a surrogate of the Monte Carlo NUBEAM module, predicts beam-driven heating, current, and torque. EC heating and current drive use a control-oriented, empirically scaled source model; plasma resistivity follows the Spitzer formulation; bootstrap current uses the Sauter model. Equilibrium is computed using both prescribed and fixed-boundary solvers (FBSs), and the pedestal structure is modeled with an empirical pedestal model. For a representative DIII-D discharge, COTSIM predicts electron and ion temperature and safety-factor profiles in close agreement with TRANSP predictive and interpretive simulations while extending predictions through the pedestal region to the plasma edge (versus 80% of the minor radius in TRANSP). Furthermore, the equivalent COTSIM simulation runs in under 3 min compared to about 2 h for TRANSP, enabling rapid scenario planning, optimization of tokamak operation, and between-pulse control design.

Control-oriented tokamak simulator (COTSIM)↗

DOE OSTI2025

Fast ion stabilization of tilt in large radius FRCs

The field reversed configuration (FRC) has been a curious case in plasma physics research in that early MHD analysis suggested FRCs should be grossly unstable, while experimental results contradicted that prediction. Later, this theory was able to resolve this contradiction by understanding that finite Larmor radius effects largely negated the MHD predictions. Similarly, previous theoretical studies of beam driven FRCs predicted that such system would be unstable to beam driven modes while, again, experimental results indicated the contradiction. In this paper, we reconcile the theoretical understanding of beam driven modes with experimental observations of stability in these systems. By self-consistently capturing fast ion generation from neutral beam injection and its impact on the plasma equilibrium, we show that low amplitude perturbations in the magnetic field, driven by betatron particles, modify the precession frequencies of the betatron particles such that the drive for compressional Alfvén waves in the thermal plasma is reduced. Finally, we are able to demonstrate, for the first time, stable beam driven FRC evolution at high S*/E in 3D kinetic simulations.

Betatrons↗

DOE OSTI2023

Synchrotron light source X-ray detection with Low-Gain Avalanche Diodes

Abstract Low Gain Avalanche Diodes (LGADs) represent the state-of-the-art in timing measurements and will instrument the future Timing Detectors of ATLAS and CMS for the High-Luminosity LHC. While initially conceived as a sensor for charged particles, the intrinsic gain of LGADs makes it possible to detect low-energy X-rays with good energy resolution and excellent time resolution (tens of picoseconds). Using the Stanford Synchrotron Radiation Lightsource (SSRL) at SLAC, several LGADs designs were characterized with energies from 5 to 70 keV. The SSRL provides 10 ps pulsed X-ray bunches separated by 2 ns intervals with an energy dispersion (ΔE/E) of 10 -4 . LGADs from Hamamatsu Photonics (HPK) and Brookhaven National Laboratory (BNL) with different thicknesses ranging from 20 μm to 50 μm and different gain layer designs were read out using fast amplification boards and digitized with a high bandwidth and high sampling rate oscilloscope. PIN devices from HPK and AC-LGADs from BNL were characterized as well. A systematic and detailed characterization of the devices' energy linearity, resolution, and time resolution as a function of X-ray energy was performed for different biasing voltages at room temperature and are reported in this work. The charge collection and multiplication mechanism were simulated using Geant4 and TCAD Sentaurus, providing an important handle for interpreting the data.

Instruments & Instrumentation↗

DOE OSTI2021

Characterization of Self-Magnetic Pinch (SMP) radiographic diode performance on RITS-6 at Sandia National Laboratories: 1) Diode Dynamics, DC Heating to extend Radiation Pulse

Radiographic diodes focus an intense electron beam to a small spot size to minimize the source area of energetic photons for radiographic interrogation. The self-magnetic pinch (SMP) diode has been developed as such a source and operated as a load for the RITS-6 Inductive Voltage Adder (IVA) driver. While experiments support the generally accepted conclusion that a 1:1 aspect diode (cathode diameter equals anode-cathode gap) delivers optimum SMP performance, such experiments also show that reducing the cathode diameter, while reducing spot size, also results in reduced radiation dose, by as much as 50%, and degraded shot reproducibility. Analyzation of the effective electron impingement angle on the anode converter with time made possible by a newly developed dose-rate array diagnostic indicates that fast-developing oscillations of the angle are correlated with early termination of the radiation pulse on many of the smaller-diameter SMP shots. This behavior as a function of relative cathode size persists through experiments with output voltages and currents up to 11.5 MV and 225 kA, respectively, and with spot sizes below ~ few mm. Since simulations to date have not predicted such oscillatory behavior, considerable discussion of the angle-behavior of SMP shots is made to lend credence to the inference. There is clear anecdotal evidence that DC heating of the SMP diode region leads to stabilization of this oscillatory behavior. This is the first of two papers on the performance of the SMP diode on the RITS-6 accelerator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Automatic Volume Balancing for Online Refueling in Molten Salt Reactor Simulations with SCALE

This work introduces recently implemented capabilities in the SCALE code system’s TRITON reactor physics sequence that improve molten salt reactor (MSR) modeling: (1) a volume balancing option, which automatically balances the volumes of the fed material with a corresponding material removal to maintain fixed mixture volumes in the neutron transport model and accurate densities, and (2) continuous feed from mixtures, which enables users to define material feed streams directly from salt mixture definitions rather than individual nuclides. The new capabilities were demonstrated via SCALE/TRITON simulations of two representative MSR concepts. Simulations of the 180 MWth molten chloride fast reactor—a system with a fixed fuel salt volume in the reactor core—applied continuous refueling with fresh fuel salt. The results confirm approximately constant reactivity when the new volume balancing capability is used. Additionally, excellent agreement with a manual feed-and-drain method for volume balancing further verified the new implementation. Simulations of the 400 MWth EIRENE reactor, an integral MSR in which the salt volume may grow over time within the reactor core, demonstrated the capability to represent growing salt volume within the TRITON depletion calculation. Compared with reference solutions, the results show consistent trends in reactivity and isotopic evolution, confirming that these new user-friendly capabilities provide accurate, physically consistent approaches for modeling MSRs in SCALE.

Elzohery, Rabab [ORNL] (ORCID:0000000160043633)↗

DOE OSTI2020

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Modeling the influence of the solid electrolyte interphase on the sand’s time and dendrite formation on lithium metal electrodes

Lithium metal is a sought after battery material for its high energy density due to the low electrochemical potential and density. However, lithium metal is also highly reactive, which results in a strong propensity for dendrite formation. The Sand’s time has previously been used to predict the time of dendrite initiation on metals that do not form a solid-electrolyte interphase (SEI), but it has been shown that the Sand’s time is not accurate for lithium electrodes when using transport parameters associated with the electrolyte. Thus, we built a numerical model to simulate lithium ion transport through a growing SEI to predict the Sand’s time. The numerical model is shown to be more accurate than previous analytical solutions, especially for low current densities. We then analyze the sensitivity of the Sand’s time to different SEI properties and the chemical potential gradients present in the SEI, driving lithium transport. The results showed that high lithium concentration has a greater impact at high current density, while fast diffusivity is more important at low current density. Lastly, we modeled the influence of surface roughness on the plating evolution and chemical potential gradients when an SEI is present in comparison to the electrolyte. As a result, we demonstrate that the SEI plays a critical role in lithium electrode stability, and that improved characterization techniques are needed to better understand transport through the SEI and increase lithium metal utilization in energy storage devices.

Chemistry↗

DOE OSTI2026

Accelerating laser ray tracing in high fidelity physics simulations of laser melting using squeeze U-net

Laser melting is a core component of the ongoing industrial revolution, dubbed Industry 4.0, as lasers facilitate fast and precise melting and fusion in advanced manufacturing. There is a strong need to optimize the laser process using simulations. However, this has proven challenging as high fidelity simulations are needed for predictive modeling and this is currently prohibitively expensive even when run on hundreds of processors on high performance computers. The challenge is capturing complex physics of laser material interaction, fluid dynamics, thermal physics and material phase transformations at various length and time scales. To close this technological gap, we modified a squeeze U-net to accelerate the laser ray tracing component of such high fidelity models by ~4x–40x while preserving the core physics principle of conservation of energy with 97% accuracy. This approach enables the accurate modeling of global laser energy absorption as a function of local surface temperatures and complex surface topologies, which govern the reflection directions and energy losses of laser rays upon interacting with the material surface.

Computer science↗

DOE OSTI2022

Characterization of self-magnetic pinch (SMP) radiographic diode performance on RITS-6 at Sandia National Laboratories. I. Diode dynamics, DC heating to extend radiation pulse

Radiographic diodes focus on an intense electron beam to a small spot size to minimize the source area of energetic photons for radiographic interrogation. The self-magnetic pinch (SMP) diode has been developed as such a source and operated as a load for the six-cavity radiographic integrated test stand (RITS-6) inductive voltage adder driver. While experiments support the generally accepted conclusion that a 1:1 aspect diode (cathode diameter equals anode–cathode gap) delivers optimum SMP performance, such experiments also show that reducing the cathode diameter, while reducing spot size, also results in reduced radiation dose, by as much as 50%, and degraded shot reproducibility. Analysis of the effective electron impingement angle on the anode converter with time made possible by a newly developed dose-rate array diagnostic indicates that fast-developing oscillations of the angle are correlated with early termination of the radiation pulse on many of the smaller-diameter SMP shots. This behavior as a function of relative cathode size persists through experiments with output voltages and currents up to 11.5 MV and 225 kA, respectively, and with spot sizes below approximately few millimeters. Since simulations to date have not predicted such oscillatory behavior, considerable discussion of the angle behavior of SMP shots is made to lend credence to the inference. We report there is clear anecdotal evidence that DC heating of the SMP diode region leads to stabilization of this oscillatory behavior. This is the first of two papers on the performance of the SMP diode on the RITS-6 accelerator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Efficient Ground-State Cooling of Large Trapped-Ion Chains with an Electromagnetically-Induced-Transparency Tripod Scheme

Here, we report the electromagnetically-induced-transparency (EIT) cooling of a large trapped 171 Yb + ion chain to the quantum ground state. Unlike conventional EIT cooling, we engage a four-level tripod structure and achieve fast sub-Doppler cooling over all motional modes. We observe simultaneous ground-state cooling across the complete transverse mode spectrum of up to 40 ions, occupying a bandwidth of over 3 MHz. The cooling time is observed to be less than 300 μ s , independent of the number of ions. Such efficient cooling across the entire spectrum is essential for high-fidelity quantum operations using trapped ion crystals for quantum simulators or quantum computers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Accelerated molecular dynamics simulations of dislocation climb in nickel

We report the mechanical behavior of materials operating under high temperatures is strongly influenced by creep mechanisms such as dislocation climb, which is controlled by the diffusion of vacancies. However, atomistic simulations of these mechanisms have traditionally been impractical due to the long time scales required. To overcome these time scale challenges, we use Parallel Trajectory Splicing (ParSplice), an accelerated molecular dynamics method, to simulate dislocation climb in nickel. We focus on modeling the activity of a vacancy near a jog on an edge dislocation in order to observe vacancy pipe diffusion and vacancy absorption at the jog. From rigorously constructed trajectories encompassing more than 2000 vacancy absorption events over a simulation time of more than 4μs at 900 K, a comprehensive sampling of available atomistic mechanisms is collated and analyzed further with molecular statics calculations. We estimate average rates for pipe diffusion and vacancy absorption into the jog using data from the dynamic and static calculations, finding very good agreement. Our results strongly suggest that the dominant mechanism for vacancy absorption by jogs is via biased diffusion to the dislocation core followed by fast pipe diffusion to the jog.

36 MATERIALS SCIENCE↗