Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dislocation dynamics”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A theoretical study of solid solution strengthening in the refractory medium entropy alloy Nb 45 Ta 25 Ti 15 Hf 15

The refractory medium-entropy alloy (RMEA) Nb 45 Ta 25 Ti 15 Hf 15 exhibits exceptional tensile ductility and fracture toughness at ambient temperature, but its engineering applications are limited by a lack of high temperature strength. Using a machine-learning interatomic potential (MLIP) with near-density functional theory (DFT) accuracy, we conducted molecular dynamics (MD) and statics simulations of the behavior of dislocations with both screw and edge characters. We also analyze experimentally measured yield strengths using the Rao-Suzuki model and the Maresca-Curtin model modified to include a temperature-dependent shear modulus and a bulk modulus-dependent misfit volume, thereby uncovering the mechanisms underlying the yielding of this RMEA. Compared with the published experimental yield strength, the models parameterized by the MLIP effectively reproduce the experimental results over a wide temperature range. The models and MD simulations indicate that yielding is governed by screw dislocations, with dipole dragging as the dominant mechanism. In MD simulations, we observed a potential softening mechanism not considered by the Rao-Suzuki screw model: slow migration of interstitial jogs along the dislocation core, which could lead to the annihilation of vacancy and interstitial jog pairs by their combination.

BCC complex concentrated alloys↗
DOE OSTI2025

Effect of reversible dislocation-based deformation on nanoparticle strain at failure

Molecular dynamics simulations are used to isolate and quantify the reversible and irreversible mechanisms that contribute to deformation in platinum nanoparticles under compression. Quantitative analysis reveals how the nucleation and entanglement of dislocations can lead to reversible dislocation-based deformation. Simulations run at different temperatures and loading conditions show that the formation of entangled dislocations is more likely at higher temperatures and is facilitated by loading orientations where dislocations nucleate on intersecting slip planes. The presence of entangled dislocations increases the strain at failure due to the ability of those dislocations to accommodate strain reversibly. The results are corroborated by the observation of similar entangled dislocation loops during in situ compression experiments on nanoparticles of the same material. Altogether, these findings provide insight into the role of dislocations in both reversible and irreversible deformation and their implications for nanoparticle stability and properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effect of helium bubbles on the mobility of edge dislocations in copper

Abstract Helium bubbles can form in materials upon exposure to irradiation. It is well known that the presence of helium bubbles can cause changes in the mechanical behavior of materials. To improve the lifetime of nuclear components, it is important to understand deformation mechanisms in helium-containing materials. In this work, we investigate the interactions between edge dislocations and helium bubbles in copper using molecular dynamics (MD) simulations. We focus on the effect of helium bubble pressure (equivalently, the helium-to-vacancy ratio) on the obstacle strength of helium bubbles and their interaction with dislocations. Our simulations predict significant differences in the interaction mechanisms as a function of helium bubble pressure. Specifically, bubbles with high internal pressure are found to exhibit weaker obstacle strength as compared to low-pressure bubbles of the same size due to the formation of super-jogs in the dislocation. Activation energies and rate constants extracted from the MD data confirm this transition in mechanism and enable upscaling of these phenomena to higher length-scale models.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dislocation nano-hydrides in nickel: Nucleation, evolution and effects on dislocation behaviors

Nano-hydrides have been predicted to precipitate at the core of edge dislocations in the Ni-H system, a mechanism that may promote hydrogen embrittlement. However, nano-hydride nucleation, growth, and effects on dislocation behavior have seldom been explored. This work combines molecular dynamics grand canonical Monte Carlo (MD-GCMC) simulations and continuum modeling to uncover a wide range of phenomena linked to dislocation nano-hydrides. Simulations reveal that nano-hydrides can be stabilized at dislocation cores with all character angles, including screw segments, due to the hydrostatic stresses around the cores of the Shockley partials. Nano-hydride nucleation takes place in these regions, and growth is dictated by the character angle θ of the perfect dislocation. The equilibrium stacking fault width d eq varies dynamically to increase the local hydrostatic stress field and facilitate the formation of the nano-hydride, forming a constriction-like feature and leading to three distinct behaviors: d eq decreases for θ > 30°, d eq remains unchanged for θ = 30°, and d eq increases for θ < 30°. Remote hydrostatic and Escaig stresses are also shown to influence the nucleation stage, implying stress concentrations such as those ahead of crack tips may facilitate nano-hydride precipitation. Moreover, we identify a new hydrogen-induced 60° dislocation reaction that emits a Shockley partial on a conjugate plane, with potential implications for twin nucleation. Testable predictions from this study are then used to reinterpret previous results from the literature. In conclusion, these findings provide a comprehensive framework to assess nano-hydride formation and evolution at dislocations in nickel and other face-centered cubic metals, with important implications to hydrogen embrittlement.

Dislocations↗
NASA NTRS1992

Prediction of dislocation generation during Bridgman growth of GaAs crystals

Dislocation densities are generated in GaAs single crystals due to the excessive thermal stresses induced by temperature variations during growth. A viscoplastic material model for GaAs, which takes into account the movement and multiplication of dislocations in the plastic deformation, is developed according to Haasen's theory. The dislocation density is expressed as an internal state variable in this dynamic viscoplastic model. The deformation process is a nonlinear function of stress, strain rate, dislocation density and temperature. The dislocation density in the GaAs crystal during vertical Bridgman growth is calculated using a nonlinear finite element model. The dislocation multiplication in GaAs crystals for several temperature fields obtained from thermal modeling of both the GTE GaAs experimental data and artificially designed data are investigated.

Tsai, C. T.↗
NASA NTRS1999

Study of Bulk and Elementary Screw Dislocation Assisted Reverse Breakdown in Low-Voltage (< 250 V) 4H-SiC p(sup +)n Junction Diodes--Part II: Dynamic Breakdown Properties: Dynamic Breakdown Properties - Part 2

This paper outlines the dynamic reverse-breakdown characteristics of low-voltage (<250 V) small-area <5 x 10(exp -4) sq cm 4H-SiC p(sup +)n diodes subjected to nonadiabatic breakdown-bias pulsewidths ranging from 0.1 to 20 microseconds. 4H-SiC diodes with and without elementary screw dislocations exhibited positive temperature coefficient of breakdown voltage and high junction failure power densities approximately five times larger than the average failure power density of reliable silicon pn rectifiers. This result indicates that highly reliable low-voltage SiC rectifiers may be attainable despite the presence of elementary screw dislocations. However, the impact of elementary screw dislocations on other more useful 4H-SiC power device structures, such as high-voltage (>1 kV) pn junction and Schottky rectifiers, and bipolar gain devices (thyristors, IGBT's, etc.) remains to be investigated.

Neudeck, Philip G.↗
NASA NTRS1990

The effect of frictional stress on the calculation of critical thickness in epitaxy

A model is developed for lattice-mismatched epitaxial layers, which makes it possible to calculate both the equilibrium critical thickness as well as the kinetic rate of misfit accommodations. The technique for calculating the equilibrium critical thickness is similar to that of Matthews (1966) in that it uses the force balance approach and the Volterra description of the dislocation line. For the kinetic component of misfit accommodation, the model uses a dynamic frictional force. A simple criterion is proposed for determining whether the generation of misfit dislocations will be controlled by the equilibrium or the kinetic factors.

Fox, Bradley A.↗
DOE OSTI2025

Dynamical Scaling Reveals Topological Defects and Anomalous Evolution of a Photoinduced Phase Transition

Nonequilibrium states of quantum materials can exhibit exotic properties and enable unprecedented functionality and applications. These transient states are inherently inhomogeneous, characterized by the formation of topologically protected structures, requiring nanometer spatial resolution on femtosecond timescales to resolve their evolution. Using ultrafast total x-ray scattering at a free electron laser and a sophisticated scaling analysis, we gain unique access to the dynamics on the relevant mesoscopic length scales. Our results provide direct evidence that ultrafast excitation of LaTe 3 leads to formation of topological vortex strings of the charge density wave. These dislocations of the charge density wave exhibit anomalous, subdiffusive dynamics, slowing the equilibration process, providing rare insight into the nonequilibrium mesoscopic response in a quantum material. Our findings establish a general framework to investigate properties of topological defects, which are expected to be ubiquitous in nonequilibrium phase transitions and may arrest equilibration and enhance competing orders.

Orenstein, Gal [SLAC National Accelerator Laborato↗
NASA NTRS2005

Growth from Solutions: Kink dynamics, Stoichiometry, Face Kinetics and stability in turbulent flow

1. Kink dynamics. The first segment of a polygomized dislocation spiral step measured by AFM demonstrates up to 60% scattering in the critical length l*- the length when the segment starts to propagate. On orthorhombic lysozyme, this length is shorter than that the observed interkink distance. Step energy from the critical segment length based on the Gibbs-Thomson law (GTL), l* = 20(omega)alpha/(Delta)mu is several times larger than the energy from 2D nucleation rate. Here o is tine building block specific voiume, a is the step riser specific free energy, Delta(mu) is the crystallization driving force. These new data support our earlier assumption that the classical Frenkel, Burton -Cabrera-Frank concept of the abundant kink supply by fluctuations is not applicable for strongly polygonized steps. Step rate measurements on brushite confirms that statement. This is the1D nucleation of kinks that control step propagation. The GTL is valid only if l* <Dk/vk, the diffusion path of a kink that has diffusivity Dk and average growth velocity vk. This is equivalent to supersaturations sigma less than approx. alpha/2l*, where alpha is the building block size. For lysozyme, sigma much less than (1%). Conventionally used interstep distance generated by screw dislocation, 19(omega)alpha/Delta(mu) should be replaced by the very different real one, approx.4l*. 2. Stoichiometry. Kink, and thus step and face rates of a non-Kossel complex molecular monocomponent or any binary, AB, lattice was found theoretically to be proportional to 1/(zeta(sup 1/2) + zeta(sup - 1/2)), where zeta = [B]/[A] is the stoichiometry ratio in solution. The velocities reach maxima at zeta = 1. AFM studies of step rates on CaOxalate monohydrate (kidney stones) from aqueous solution was found to obey the law mentioned above. Generalization for more complex lattice will be discussed. 3. Turbulence. In agreement with theory, high precision in-situ laser interferometry of the (101) KDP crystal face shows step bunching if solution flows parallel to the step flow. The bunch height increases with the distance the bunch travels, i.e. with the face size. However, when the flow rate, u, increases, at u greater than approx. 1 m / s , the average step bunch height decreases as 1/u. The pheonomenon is attributed to the turbulent rather than laminar viscous boundary layer where diffusivity Dt = 0.5u(sub tau),y, i.e. increases linearly with the distance y from the solid face. Friction velocity, u(sub tau) approx. u(sup 7/8). Dramatically larger rate of the mass/heat transport within the turbulent, as compared to the laminar, viscous layer will be discussed.

Chernov, A. A.↗
DOE OSTI2025

Interplay between hydrogen, temperature, and character angle on dissociated dislocation energies in Fe–Ni–Cr austenitic stainless steels

Dislocation energy has an important role in the mechanical performance of structural metals. While dislocation energies cannot be fully obtained from continuum theories due to the contribution of the dislocation core, they have been calculated via atomistic simulations in elemental metals. However, constraints on the local atomic environments have prevented the use of such approaches in systems that incorporate alloying or interstitial solutes. In this work, we develop robust molecular dynamics methods to resolve these issues through a geometric construction of dislocation dipoles and the calculation of time-averaged energies. Furthermore, we apply these methods to calculate dislocation energies (including core energies) in an Fe 70 Ni 11 Cr 19 austenitic steel at a variety of character angles, hydrogen concentrations, temperatures, and dipole spacings. The resulting highly converged energies show an excellent agreement with continuum expressions. Overall, hydrogen concentrations up to 1.0 % do not have a significant effect on the elastic parameters and dislocation energy. The methods and insights derived in this work have the potential to facilitate the calculation of dislocation energies in a wide range of systems, and to guide our understanding of hydrogen embrittlement.

Alloyed systems↗
DOE OSTI2025

Additively-manufactured Al-0.3Zr-0.2Ce-0.2Cu alloy with high creep resistance and electrical conductivity

Here, a new, solute-lean Al-0.3Zr-0.2Ce-0.2Cu (wt.%) alloy is developed for additive manufacturing that overcomes the classical tradeoff between conductivity and creep resistance. The rapid-cooling-enabled supersaturation of Zr, and its uniform distribution in α-Al matrix, along with formation of submicron (Ce,Cu)-rich intermetallic particles on solidification lead to unusually high creep resistance at 200 °C. Near-zero secondary creep rates are achieved up to the alloy yield stress (YS) of 65 MPa at 200 °C in as-fabricated state. The Zr-solute-induced dislocation-climb suppression mechanism underlying this improvement also restricts dynamic recovery above YS, as noted from appreciable primary creep and its transitioning to near-zero secondary creep rates. A combination of relatively coarse, epitaxially-grown α-Al grains, low Zr concentration in α-Al, and the impurity-scavenging effect of Ce to purify α-Al matrix produces high electrical conductivity of ∼48 %IACS. Aging precipitation of L1 2 -Al 3 Zr nanoprecipitates doubles the YS (to ∼150 MPa) at room temperature and increases alloy conductivity to ∼58 %IACS, but loss of solid-solution Zr out of α-Al matrix leads to activation of dislocation climb, degrading the creep properties as compared to the supersaturated Al-Zr solid solution in the as-fabricated state. Compared to L1 2 -Al 3 Zr nanoprecipitates, submicron (Ce,Cu)-rich particles formed on solidification are more effective at impeding dislocation climb, producing a threshold stress for dislocation creep of ∼ 50 MPa at 200 °C. The new alloy design concepts, especially solute-induced dislocation-climb suppression for creep resistance, explored here may pave way for the design of new metallic alloys for thermal/electrical conductors and other high-temperature applications.

Additive Manufacturing↗
DOE OSTI2025

Dislocation-Grain Boundary Interaction Dataset for FCC Cu

Interactions between dislocations and grain boundaries play a major role in controlling the strength and ductility of structural materials. Experimentally, assessing and probing geometric and stress-based criteria at the local level for dislocation transmission through grain boundaries remains challenging. Therefore, there have been many efforts to systematically generate datasets of dislocation-grain boundary interactions (DGI) via computational models such as molecular dynamics simulations. So far, most DGI datasets have focused only on the subset of nominal minimum-energy grain boundary structures, which limits their applicability, especially to materials processed far from equilibrium. We present a comprehensive database of dislocation-grain boundary interactions for edge, screw, and 60° mixed dislocation with 330 <110> and 257 <112> symmetric tilt grain boundaries (total of 587) in FCC Cu consisting of 73 minimum-energy grain boundary structures and 514 metastable structures. The dataset contains the outcomes for 5234 unique interactions for various dislocation types, grain boundary structures, and applied shear stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Anomalous entropy-driven kinetics of dislocation nucleation

The kinetics of dislocation reactions, such as dislocation multiplication, controls the plastic deformation in crystals beyond their elastic limit, therefore critical mechanisms in a number of applications in materials science. We present a series of large-scale molecular dynamics simulations that shows that one such type of reactions, the nucleation of dislocation at free surfaces, exhibit unconventional kinetics, including unexpectedly large nucleation rates under compression, very strong entropic stabilization under tension, as well as strong non-Arrhenius behavior. These unusual kinetics are quantitatively rationalized using a variational transition state theory approach coupled with an efficient numerical scheme for the estimation of vibrational entropy changes. These results highlight the need for a variational treatment of the kinetics to quantitatively capture dislocation reaction kinetics, especially at low-to-moderate strains where large deformations are required to activate reactions. These observations suggest possible explanations to previously observed unconventional deformation kinetics in both molecular dynamics simulations and experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Interpretation of Ion Irradiation and Neutron Irradiation Damage in Additively Manufactured 316 Stainless Steel using Multiscale Modeling

The accelerated adoption of nuclear energy necessitates advanced manufacturing technologies, such as additive manufacturing, to meet heightened supply chain requirements and support innovative reactor technologies. Due to the unique microstructural characteristics of additively manufactured materials under distinct solidification conditions, comprehensive evaluation of their performance in reactor environments is essential. The Advanced Materials and Manufacturing Technologies program under the Department of Energy's Office of Nuclear Energy focuses on understanding the irradiation performance and damage evolution of laser powder bed fusion 316 stainless steel, with an emphasis on integrating ion and neutron irradiation data to accelerate the development and qualification of materials for advanced nuclear reactor applications. While ion irradiation is a cost- and time-effective method, modeling and simulation are required to interpret the data for the broader range of irradiation conditions encountered in advanced reactors. In fiscal year 2025, integrated multiscale modeling and simulations were conducted to assess irradiation damage in additively manufactured 316 stainless steel. Key outcomes include predictions of chromium enrichment at grain boundaries, nickel enrichment at dislocation cell walls and void surfaces, and heterogeneous void evolution under ion and neutron irradiation conditions. Cluster dynamics simulations revealed the coarsening of voids at high irradiation temperatures and the suppression of void growth by high network dislocation density, while also demonstrating significant growth and coarsening of voids and self-interstitial atom loops at low dose rates. Machine learning-accelerated atomistic simulations highlighted the impact of the local environment and chromium concentration on vacancy diffusivity, providing key insights on the influence of composition on void swelling and radiation-induced segregation. Additionally, molecular dynamics simulations demonstrated the presence of defect production bias and a significant effect of carbon content on defect cluster behavior. These combined efforts aim to predict the performance of additively manufactured materials under various reactor conditions, supporting their qualification for nuclear reactor applications by interpreting ion irradiation data. This report underscores the potential of integrated multiscale modeling to analyze ion irradiation data in the effort to accelerate the qualification of additively manufactured materials for nuclear reactor components.

316 stainless steel↗
NASA NTRS2012

Upper-Ocean Heat Balance Processes and the Walker Circulation in CMIP5 Model Projections

Considerable uncertainty remains as to the importance of mechanisms governing decadal and longer variability of the Walker Circulation, its connection to the tropical climate system, and prospects for tropical climate change in the face of anthropogenic forcing. Most contemporary climate models suggest that in response to elevated CO2 and a warmer but more stratified atmosphere, the required upward mass flux in tropical convection will diminish along with the Walker component of the tropical mean circulation as well. Alternatively, there is also evidence to suggest that the shoaling and increased vertical stratification of the thermocline in the eastern Pacific will enable a muted SST increase there-- preserving or even enhancing some of the dynamical forcing for the Walker cell flow. Over the past decade there have been observational indications of an acceleration in near-surface easterlies, a strengthened Pacific zonal SST gradient, and globally-teleconnected dislocations in precipitation. But is this evidence in support of an ocean dynamical thermostat process posited to accompany anthropogenic forcing, or just residual decadal fluctuations associated with variations in warm and cold ENSO events and other stochastic forcing? From a modeling perspective we try to make headway on this question by examining zonal variations in surface energy fluxes and dynamics governing tropical upper ocean heat content evolution in the WCRP CMIP5 model projections. There is some diversity among model simulations; for example, the CCSM4 indicates net ocean warming over the IndoPacific region while the CSIRO model concentrates separate warming responses over the central Pacific and Indian Ocean regions. The models, as with observations, demonstrate strong local coupling between variations in column water vapor, downward surface longwave radiation and SST; but the spatial patterns of changes in the sign of this relationship differ among models and, for models as a whole, with observations. Our analysis focuses initially on probing the inter-model differences in energy fluxes / transports and Walker Circulation response to forcing. We then attempt to identify statistically the El Nino- / La Nina-related ocean heat content variability unique to each model and regress out the associated energy flux, ocean heat transport and Walker response on these shorter time scales for comparison to that of the anthropogenic signals.

Robertson, F. R.↗
DOE OSTI2025

Second harmonic Bessel-Gauss beam shaping with elliptic axicon aberrations

The second harmonic (SH) of an axicon generated Bessel-Gauss beam is created through the nonlinear interaction of photons with a crystal, resulting in the energy doubling of the output photons. In this work, we show experimentally that in addition to frequency doubling, the SH of Bessel-Gauss beams under asymmetric aberrations from an elliptic axicon exhibit intriguing beam formation. Particularly, the central region of the SH beam profile is composed of two central spots of various geometries surrounded by nested ellipses; one of which is the configuration of two central gamma dots with similar radius knotted by nested ellipses for a zeroth-order Bessel-Gauss pump. These SH beams consistently maintain their spatial profile throughout propagation, reminiscing the behavior of screw dislocations in wave patterns. Our numerical simulations produce beam dynamics consistent with that of experiments and further implicate the remarkable interweaving of bright spots with dark vortices. This is especially noticeable when the beams dynamically oscillate along the optical axis, resulting in the genesis of spatially polarized beams with a knotted framework. The insights gained from our study establish a novel paradigm for exploring interactions of Bessel-like beams with vortex dynamics. This, in turn, has the potential to spark innovations in optical applications, fostering new methodologies to harness and manipulate complex light structures. Our experimental findings could spur new methods of generating logical states of light and new opportunities for material processing control. Published by the American Physical Society 2025

Wang, Tianhong (ORCID:0000000169104201)↗
DOE OSTI2025

The interaction of edge dislocations with hydrogen-helium bubbles in tungsten

In fusion reactors, plasmas-facing components undergo degradation due to displacement damage and injection of gas impurities. Significant amounts of hydrogen (H) and helium (He) impurities can be introduced into materials through plasma exposure and nuclear transmutation, and their synergistic effects lead to the formation of mixed He/H bubbles under irradiation, changing their hardening mechanisms. In this study, molecular dynamics (MD) models were developed to study the interaction between ½<111> edge dislocations with mixed He/H-bubbles in W. The critical resolved shear stress (CRSS) that is required for dislocation breakaway increases with He/H-to-vacancy ratio at 600 K, indicating strong pinning effects of mixed He/H bubbles. However, as the temperature increases up to 1400 K, the mixed He/H bubbles become unstable as the H atoms are increasingly emitted from the mixed He/H bubbles and migrate into W matrix, which lowers bubble pressure and CRSS. Overall, these results highlight the synergistic effects of He and H on dislocation behavior in irradiated W during deformation.

36 - MATERIALS SCIENCE↗