Search NASA⌕ Search

SEARCH · Search NASA

Results for “material response 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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Impact of Integration Scheme on Performance of Anisotropic Plasticity Models

Given the prevalent role of metals in a variety of industries, schemes to integrate corresponding constitutive models in finite element applications have long been studied. A number of formulations have been developed to accomplish this task; each with their own advantages and costs. Often the focus has been on ensuring the accuracy and numerical stability of these algorithms to enable robust integration. While important, emphasis on these performance metrics may often come at the cost of computational expense potentially neglecting the needs of individual problems. In the current work, the performance of two of the most common integration methods for anisotropic plasticity -- the convex cutting plane (CCP) and closest point projection (CPP) -- across a variety of metrics is assessed; including accuracy and cost. A variety of problems are considered ranging from single elements to large representative simulations including both implicit quasistatic and explicit transient dynamic type responses. The relative performance of each scheme in the different instances is presented with an eye towards guidance on when the different algorithms may be beneficial.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Predicting Failure Using Deep Learning SAND Report

Accurate prediction of ductile failure is critical to Sandia’s NW mission, but the models are computationally heavy. The costs of including high-fidelity physics and mechanics that are germane to the failure mechanisms are often too burdensome for analysts either because of the person-hours it requires to input them or because of the additional computational time, or both. In an effort to deliver analysts a tool for representing these phenomena with minimal impact to their existing workflow, our project sought to develop modern data-driven methods that would add microstructural information to business-as-usual calculations and expedite failure predictions. The goal is a tool that receives as input a structural model with stress and strain fields, as well as a machine-learned model, and output predictions of structural response in time, including failure. As such, our project spent substantial time performing high-fidelity, three-dimensional experiments to elucidate materials mechanisms of void nucleation and evolution. We developed crystal-plasticity finite-element models from the experimental observations to enrich the findings with fields not readily measured. We developed engineering length-scale simulations of replicated test specimens to understand how the engineering fields evolve in the presence of fine-scale defects. Finally, we developed deep learning convolutional neural networks, and graph-based neural networks to encode the findings of the experiments and simulations and make forward predictions in time for structural performance. This project demonstrated the power of data-driven methods for model development, which have the potential to vastly increase both the accuracy and speed of failure predictions. These benefits and the methods necessary to develop them are highlighted in this report. However, many challenges remain to implementing these in real applications, and these are discussed along with potential methods for overcoming them.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Refining a Novel Process Monitoring Method to Safeguard Continuously Cycling Designs Using Isotopic Ratios

In advanced reactor (AR) designs, a common feature is continuous chemical processing and circulation of the nuclear material. This work bridges a significant measurement gap in safeguarding reactors with circulating fuel or continuous refueling by leveraging and building on the isotope ratio method first developed by our team under an FY21 Advanced Reactors International Safeguards Engagement (ARISE) project (Uribe et al. 2021). In circulating fuel designs, the radioisotope inventory changes from traditional effects (e.g., radioactive decay, fission) but also includes material transport due to pressure and temperature gradients. Such designs may also require regular or continuous additions or removals during operation, which significantly increases the rate of inventory change compared to traditional pressurized water reactor (PWR)s. Thus, directly tracking the nuclear inventory is ineffective since the isotopes are continuously added and removed. The isotope ratio method instead focuses on detecting changes to the input and output flows of radioisotopes. Previous work showed that for well-chosen pairs of isotopes, the isotopic ratio provides a sensitive and lasting indicator of deviation from normal conditions (e.g., startup, shutdown, diversion). The isotope ratio method is a process monitoring method with potential for application in for forward-looking approaches to International Atomic Energy Agency (IAEA) safeguards. The original process monitoring method was developed for a specific case—the decay tank of a thorium-fueled molten salt breeder reactor. In this expanded work, we explored other types of reactors and processes with nonstationary (e.g., flowing) nuclear material, which are difficult to safeguard with traditional methods because of the transient nature of the systems. The goal of the present work is to generalize the isotope ratio method for use in processes with continuously flowing nuclear material. All continuous processes have an average time for isotopes to be replaced in the system. The isotope ratio method works by choosing radioisotopes with half-lives both above and below the processing time. The present work seeks to explore which isotopes are suitable for the method by simulating the nuclear inventory, radioisotope emissions, and detector responses for several classes of advanced reactors. While the method can in principle be applied to other processes (e.g., enrichment or reprocessing facilities), the present work limits scope to ARs with continuously flowing fuel. Section 2 details the mathematics supporting the isotope ratio method, and Section 3 introduces the representative ARs selected for this work. Section 4 discusses how each reactor was analyzed, and Section 5 showcases the results for each representative reactor. Finally, Section 6 provides concluding remarks and suggests pathways for further analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
NASA NTRS1985

Studies of erosion of solar max samples of Kapton and Teflon

Several samples of Kapton and Teflon which was exposed to solar radiation were examined. The samples represent material behavior in near Earth space. Clues to the identity of erosive processes and the responsible species were searched for. Interest centered around oxygen atoms which are ubiquitous at these altitudes and are known to erode some metal surfaces. Three diagnostic methods were employed: optical microscopy, scanning electron microscopy, and fourier transform infrared spectroscopy. Two types of simulation were used: a flow containing low energy oxygen atoms and bombardment with 3000 volt Ar ions. Results and conclusions are presented.

Fristrom, R. M.↗
NASA NTRS1998

Finite Element and Plate Theory Modeling of Acoustic Emission Waveforms

A comparison was made between two approaches to predict acoustic emission waveforms in thin plates. A normal mode solution method for Mindlin plate theory was used to predict the response of the flexural plate mode to a point source, step-function load, applied on the plate surface. The second approach used a dynamic finite element method to model the problem using equations of motion based on exact linear elasticity. Calculations were made using properties for both isotropic (aluminum) and anisotropic (unidirectional graphite/epoxy composite) materials. For simulations of anisotropic plates, propagation along multiple directions was evaluated. In general, agreement between the two theoretical approaches was good. Discrepancies in the waveforms at longer times were caused by differences in reflections from the lateral plate boundaries. These differences resulted from the fact that the two methods used different boundary conditions. At shorter times in the signals, before reflections, the slight discrepancies in the waveforms were attributed to limitations of Mindlin plate theory, which is an approximate plate theory. The advantages of the finite element method are that it used the exact linear elasticity solutions, and that it can be used to model real source conditions and complicated, finite specimen geometries as well as thick plates. These advantages come at a cost of increased computational difficulty, requiring lengthy calculations on workstations or supercomputers. The Mindlin plate theory solutions, meanwhile, can be quickly generated on personal computers. Specimens with finite geometry can also be modeled. However, only limited simple geometries such as circular or rectangular plates can easily be accommodated with the normal mode solution technique. Likewise, very limited source configurations can be modeled and plate theory is applicable only to thin plates.

Prosser, W. H.↗
NASA NTRSDate not supplied

ISS Payload Operations Training During the COVID-19 Pandemic: Impacts and Solutions

Introduction: The onset and the protracted duration of the COVID-19 pandemic brought a dramatic and rapid transformation to almost every aspect of humanity in 2020. The world’s space agencies and their missions were not immune to the unpredictable and wide-sweeping changes. One discipline principally affected was mission operations and the various groups supporting that critical function. Mission support teams, especially for complex and crewed missions like the International Space Station (ISS) were forced to rapidly rethink how and where control center staff performed their vital work. Operations training – an essential element to mission ops, had unique hurdles to overcome. Operations training is responsible for preparing astronaut crews for their scientific missions, providing recurring proficiency and currency material to flight controllers, as well as ensuring new team members are ready to join their colleagues on console. Every element of training was impacted by COVID restrictions. From orientation and introductory classes for new controllers, training for scientists and payload developers, to simulations and advanced lessons, critical On the Job Training (OJT) and final evaluations; all aspects faced new challenges. Operations trainers at NASA’s Marshall Space Flight Center (MSFC) in Huntsville, Alabama, were forced to rethink, retool, and become more efficient with trainees and training resources in order to continue supporting payload operations on the ISS. Impacts and Barriers to Operations Training: The pandemic started in the United States in early March 2020. Immediately, NASA mandated that the support for ISS real-time operations was critical. As a result, physical access to the ISS Payload Operations Center was restricted to only those people supporting console operations. Thus, use of on-site NASA facilities and any trainer to trainee face-to-face interaction was not authorized. Trainers and trainees had to quickly shift to 100% remote learning. In the short term, this was not a problem since NASA had already prepared Information Technology (IT) resources that facilitated remote access to mission resources like classes, documents and teleconferences. However, instructors quickly discovered the classes that they were accustomed to delivering in a large classroom environment did not translate well to remote teaching. The lack of non-verbal communication with the class became a major barrier to instruction. Another significant hurdle to operations training was the mandate that all simulations could only be held remotely. The logistics of planning and executing even small simulations proved to be challenging, almost impossible, due to IT technology and public internet limitations. With simulations essentially stopped, as well as the restrictions on most OJT training, trainees were essentially stopped in their advancement towards certification. An additional barrier that had to be overcome was the formal documentation of the curriculum objectives. Even though NASA had been deliberately transitioning to electronic documents for years, some organizations were still relying on “pen and ink” documents to formalize training completion. The logistics of implementing secure and controlled e-signature documents in a standardized, NASA approved format had to be addressed. Lesson Learned and Efficiencies: Once the various limitations were identified, ISS trainers at MSFC discussed, proposed, and prioritized new options and techniques. Transitioning to all electronic learning materials and certifications was a relevantly easy fix. Teaching courses to large and diverse groups took additional shifts in the NASA training paradigm. Likewise, methods and procedures for preparing astronauts for their missions were forced to be revised. Simulation supervisors and their teams were also forced to find more efficient techniques to provide realistic training experience. Communication and approvals from management was essential. In every case, the ISS payload operations instructors found novel and efficient solutions to all functions listed above. This paper will discuss the factors and solutions payloads operations trainers found to keep scientific research on the ISS flying forward to mission success.

ISS Operations↗
DOE OSTI2021

Reliability Barriers to Enhanced Solar PV Deployment: Selected Research Findings - State Outreach Technical Sessions

This presentation provides summaries of three reports selected from a body of work that examines potential reliability challenges related to planning future power systems with higher solar photovoltaic (PV) penetrations in the western United States. The three topics covered are 1) Power system flexibility requirements and supply; 2) Resource adequacy and the capacity credit of solar; 3) Simulating distributed energy resource responses to transmission system-level faults considering IEEE 1547 performance categories on three major WECC transmission paths. The full reports and related resources are available at the project website, https://www.westernenergyboard.org/western-interstate-energy-board/barrier-mitigation-to-enhanced-distributed-solar-photovoltaic/. This presentation is targeted toward a technical audience. A companion policy presentation provides higher-level summaries of the same material.

distributed energy resource↗
DOE OSTI2017

Molecular Dynamics Simulation of Fission Fragment Damage in Nuclear Fuel and Surrogate Material

We have performed classical molecular dynamics simulations of swift heavy ion damage, typical of fission fragments, in nuclear fuel (UO2) for energy deposition per unit length of 3.9 keV/nm. We did not observe amorphization. The damage mainly consisted of isolated point defects. Only about 1% of the displacements occur on the uranium sublattice. Oxygen Frenkel pairs are an order of magnitude more numerous than uranium Frenkel pairs in the primary damage state. In contrast, previous results show that the ratio of Frenkel pairs on the two sublattices is close to the stoichiometric ratio in ceria. These differences in the primary damage state may lead to differences in radiation response of UO2 and CeO2.

Molecular dynamics, radiation damage in materials,↗
DOE OSTI2021

Thermo-mechanical characterization of shale using nanoindentation

Abstract Shale can be a potential buffer for high-level radioactive nuclear wastes. To be an effective buffer while subject to waste heat, shale's mechanical response at elevated temperature must be known. Many researchers have experimentally characterized the mechanical behavior of various shales at different length scales in adiabatic conditions. However, its mechanical performance at elevated temperatures at the nano-scale remains unknown. To investigate the temperature dependency of nanomechanical properties of shale, we conducted both experimental and numerical studies. In this study, we measured mechanical and fracture properties of shale, such as hardness, elastic modulus, anisotropy, and fracture toughness from 25 °C up to 300 °C at different bedding planes. Statistical analysis of the results suggests that hardness and fracture toughness significantly increased at temperatures from 100 to 300 °C; while, temperature does not have a significant impact on elastic modulus. Data also shows that the bedding plane orientations have a substantial impact on both mechanical and fracture properties of shale at the nano-scale leading to distinct anisotropic behavior at elevated temperature below 100 °C. Additionally, we numerically investigated the mechanical performance of the shale samples at room temperature to gain an insight into its mechanical response through the thickness. Numerical results were validated against the experimental results, confirming the simulation can be used to predict shale deformation at the nano-scale or potentially be used in multi-scale simulations.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Simulating the Formation of Lunar Floor-Fracture Craters Using Elastoviscoplastic Relaxation

Lunar floor-fracture craters formed during the height of mare basalt emplacement. Due to a general temporal and spatial relation with the maria, these craters, numbering some 200, may be diagnostic of the thermal structure of the crust during this time. As the name suggests, these craters exhibit brittle failure, generally limited to the central floor region. That, and a shallower depth than fresh lunar craters, has led to two main theories as to their formation: laccolith emplacement under the crater and viscous relaxation. The implications of each model for the state of the Moon's crust during this time are quite different, so the viability of each model must be checked. Laccolith emplacement has been treated elsewhere. However, previous attempts to study the relaxation of the craters have assumed only a uniform, Newtonian viscous response of the near surface to the topographic driving forces, and simply postulated that the fractures resulted from tensile stresses associated with floor uplift. Here, we use a more sophisticated rheological model that includes not only non-Newtonian viscous behavior (i.e., the viscosity is stress-dependent), but also incorporates elastic behavior and a plastic component to the rheology to directly simulate the formation of the floor fractures. The results of our simulations show that while elastoviscoplastic relaxation is potentially viable for larger floor-fracture craters, it is not viable for craters with diameters < or = 60 km, the size of the majority of floor-fracture craters. We employ the finite element method, a numerical technique well suited for boundary-value problems, via the commercially available MARC software package. To test the viability of topographic relaxation, our goal is to prepare the simulations as to maximize the amount of relaxation. We take advantage of the natural axisymmetry of craters, simulating one radial plane. Initial shapes are based on data for fresh craters from Pike. To simplify implementation, a fourth order polynomial is used for the basin, while a third order inverse function is used for the rim. This form closely approximates the long-wavelength behavior of complex craters, while ignoring higher-frequency topography, save the rim. This approximation is appropriate because crater relaxation is strongly controlled by long-wavelength topography. Loading is accomplished assuming a uniform gravity field (1.62 m/s-square) and a uniform density of 2900 kg/cubic m. The initial stress state is set to be hydrostatic, with an additional pressure term to account for any overlying topography. This additional pressure term is tapered exponentially with depth. While the simulations quickly settle on a preferred stress state, and while the final solution is fairly insensitive to the choice of the e-folding depth of the taper, selecting an e-folding depth close to the diameter of the crater sets the initial stress state near the preferred state. We assume a diuranally averaged surface temperature of -20 C, and allow temperature to increase with depth at a rate of 50 K//km. Assuming a thermal conductivity of 2 W/in/K, this gradient translates to a heat flow of 100mW/square m, an extremely high value for the Moon. Temperature, of course, will not increase without bound. To maximize relaxation, we allow our temperature profile to increase linearly until it reaches the solidus (assumed to be 1200C) at a depth of 24.4 km, at which point it is kept constant. The presence of melt will drop the bulk viscosity; however, we have no rheological control for partial melts. Therefore, we make no attempt to simulate this situation. Elastoviscoplastic rheological model. In general, geologic materials can behave in three main ways: elastically, viscously (via solid-state creep), and brittly (plasticity is a continuum approach to simulate this phenomenon). We combine these three deformation mechanisms in an extended Maxwell solid, where the total strain can be broken down into a simple summation of the elastic, creep, and plastic strains. In relaxation phenomena in general, the system takes advantage of any means possible to eliminate deviatoric stresses by relaxing away the topography. Previous analyses have only modeled the viscous response. Comparatively, the elastic response in our model can augment the relaxation, to a point. This effect decreases as the elastic response becomes stiffer; indeed, in the limit of infinite elastic Young's modulus (and with no plasticity), the solution converges on the purely viscous solution. Igneous rocks common to the lunar near-surface have Young's modulii in the range of 10-100 GPa. To maximize relaxation, we use a Young's modulus of 10 GPa. (There is negligible sensitivity to the other elastic modulus, the Poisson's ratio; we use 0.25.) For the viscous response, we use a flow law for steady-state creep in thoroughly dried Columbia diabase, because the high plagioclase (about 70 vol%) and orthopyroxene (about 17 vol%) content is similar to the composition of the lunar highland crust as described by remote sensing and sample studies: noritic anorthosite. This flow law is highly non-Newtonian, i.e., the viscosity is highly stress dependent. That, and the variability with temperature, stands in strong contrast to previous examinations of lunar floor-fracture crater relaxation. To model discrete, brittle faulting, we assume "Byerlee's rule," a standard geodynamical technique. We implement this "rule" with an-angle of internal friction of about 40 deg, and a higher-than-normal cohesion of about 3.2 MPa (to approximate the breaking of unfractured rock). The actual behavior of geologic materials is more complex than in our rheological model, so the uncertainties in the plasticity do not represent the state-of-the-art error. Additional information is contained in the original.

Dombard, A. J.↗
DOE OSTI2017

Kinetic Monte Carlo simulations of GaN homoepitaxy on c- and m-plane surfaces

The surface orientation can have profound effects on the atomic-scale processes of crystal growth and is essential to such technologies as GaN-based light-emitting diodes and high-power electronics. We investigate the dependence of homoepitaxial growth mechanisms on the surface orientation of a hexagonal crystal using kinetic Monte Carlo simulations. To model GaN metal-organic vapor phase epitaxy, in which N species are supplied in excess, only Ga atoms on a hexagonal close-packed (HCP) lattice are considered. The results are thus potentially applicable to any HCP material. Growth behaviors on c-plane (0001) and m-plane (011¯0) surfaces are compared. We present a reciprocal space analysis of the surface morphology, which allows extraction of growth mode boundaries and direct comparison with surface X-ray diffraction experiments. For each orientation, we map the boundaries between 3-dimensional, layer-by-layer, and step flow growth modes as a function of temperature and growth rate. Two models for surface diffusion are used, which produce different effective Ehrlich-Schwoebel step-edge barriers and different adatom diffusion anisotropies on m-plane surfaces. Simulation results in agreement with observed GaN island morphologies and growth mode boundaries are obtained. These indicate that anisotropy of step edge energy, rather than adatom diffusion, is responsible for the elongated islands observed on m-plane surfaces. As a result, island nucleation spacing obeys a power-law dependence on growth rate, with exponents of –0.24 and –0.29 for the m- and c-plane, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Effect of Grain Boundary Facet Junctions on Segregation and Embrittlement

Junctions are discontinuities in flat grain boundaries that arise in all polycrystalline materials and are thought to play important roles in the response of a grain boundary network to thermal and mechanical loads. A key open question concerns the mechanisms by which solute segregation to junctions impacts properties of the grain boundary. Here, in this work, we investigate the influence of grain boundary facet junctions on solute embrittlement, and we present an analytical model that uses the hydrostatic stress field contributed by dislocations at multiple junctions to describe these effects. Specifically, we study junctions between {112} facets of various lengths in Au $\langle111\rangle$ Σ3 tilt grain boundaries. Copper and silver solutes are employed to determine if the effect of junctions on solute segregation and embrittlement is dependent on size relative to the host. Combined, atomistic simulation data and the analytical model show that Cu and Ag have opposite segregation responses to junctions due to the sign of the hydrostatic stress field induced by junctions. However, a positive shift in the embrittling potency is computed near junctions regardless of solute type or the stress state of the segregation site. Hence, for the conditions studied, junctions consistently shift the energetic landscape towards embrittlement.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Mechanical and electrochemical response of a LiCoO 2 cathode using reconstructed microstructures

As LiCoO 2 cathodes are charged, delithiation of the LiCoO 2 active material leads to an increase in the lattice spacing, causing swelling of the particles. When these particles are packed into a bicontinuous, percolated network, as is the case in a battery electrode, this swelling leads to the generation of significant mechanical stress. In this study we performed coupled electrochemical-mechanical simulations of the charging of a LiCoO 2 cathode in order to elucidate the mechanisms of stress generation and the effect of charge rate and microstructure on these stresses. Energy dispersive spectroscopy combined with scanning electron microscopy imaging was used to create 3D reconstructions of a LiCoO 2 cathode, and the Conformal Decomposition Finite Element Method is used to automatically generate computational meshes on this reconstructed microstructure. Replacement of the ideal solution Fickian diffusion model, typically used in battery simulations, with a more general non-ideal solution model shows substantially smaller gradients of lithium within particles than is typically observed in the literature. Using this more general model, lithium gradients only appear at states of charge where the open-circuit voltage is relatively constant. While lithium gradients do affect the mechanical stress state in the particles, the maximum stresses are always found in the fully-charged state and are strongly affected by the local details of the microstructure and particle-to-particle contacts. These coupled electrochemical-mechanical simulations begin to yield insight into the partitioning of volume change between reducing pore space and macroscopically swelling the electrode. Lastly, preliminary studies that include the presence of the polymeric binder suggest that it can greatly impact stress generation and that it is an important area for future research.

25 ENERGY STORAGE↗
NASA NTRS1987

A modular BLSS simulation model

The coordination of material flows in earth's biosphere is largely made possible by the buffering effect of huge material reservoirs. Without similarly-sized buffers, a bioregenerative life support system (BLSS) for extraterrestrial use will be faced with coordination problems more acute than those in any ecosystem found on earth. A related problem in BLSS design is providing an interface between the various life-support processors, one that will allow for their coordination while still allowing for system expansion. A modular model is presented of a BLSS that interfaces system processors only with the material storage reservoirs, allowing those reservoirs to act as the principal buffers in the system and thus minimizing difficulties with processor coordination. The modular nature of the model allows independent development of the detailed submodels that exist within the model framework. Using this model, BLSS dynamics were investigated under normal conditions and under various failure modes. Partial and complete failures of various components, such as the waste processor or the plants themselves, drive transient responses in the model system, allowing examination of the effectiveness of the system reservoirs as buffers. The results from simulations of this sort will help to determine control strategies and BLSS design requirements. An evolved version of this model could be used as an interactive control aid in a future BLSS.

Rummel, John D.↗
DOE OSTI2015

Tailored multivariate analysis for modulated enhanced diffraction

Modulated enhanced diffraction (MED) is a technique allowing the dynamic structural characterization of crystalline materials subjected to an external stimulus, which is particularly suited for in situ and operando structural investigations at synchrotron sources. Contributions from the (active) part of the crystal system that varies synchronously with the stimulus can be extracted by an offline analysis, which can only be applied in the case of periodic stimuli and linear system responses. In this paper a new decomposition approach based on multivariate analysis is proposed. The standard principal component analysis (PCA) is adapted to treat MED data: specific figures of merit based on their scores and loadings are found, and the directions of the principal components obtained by PCA are modified to maximize such figures of merit. As a result, a general method to decompose MED data, called optimum constrained components rotation (OCCR), is developed, which produces very precise results on simulated data, even in the case of nonperiodic stimuli and/or nonlinear responses. The multivariate analysis approach is able to supply in one shot both the diffraction pattern related to the active atoms (through the OCCR loadings) and the time dependence of the system response (through the OCCR scores). When applied to real data, OCCR was able to supply only the latter information, as the former was hindered by changes in abundances of different crystal phases, which occurred besides structural variations in the specific case considered. To develop a decomposition procedure able to cope with this combined effect represents the next challenge in MED analysis.

Caliandro, Rocco↗
DOE OSTI2022

Limiting velocities and transonic dislocations in Mg

To accurately predict the mechanical response of materials, especially at high strain rates, it is important to account for dislocation velocities in these regimes. Under these extreme conditions, it has been hypothesized that dislocations can move faster than the speed of sound. However, the presence of such dislocations remains elusive due to challenges associated with measuring these experimentally. In this work, molecular dynamics simulations were used to investigate the dislocation velocities for the basal edge, basal screw, prismatic edge, and prismatic screw dislocations in Mg in the sub-, trans-, and supersonic regimes. Also, our results show that only prismatic edge dislocations achieve supersonic velocities. Furthermore, this work demonstrates that the discrepancy between the theoretical limiting velocity and the MD results for Mg is due to its sensitivity to large hydrostatic stress around the dislocation core, which was not the case for fcc metals such as Cu.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Deformation and failure in extreme regimes by high-energy pulsed lasers: A review

The use of high-power pulsed lasers to probe the response of materials at pressures of hundreds of GPa up to several TPa, time durations of nanoseconds, and strain rates of 10 6 –10 1° s -1 is revealing novel mechanisms of plastic deformation, phase transformations, and even amorphization. This unique experimental tool, aided by advanced diagnostics, analysis, and characterization, allows us to explore these new regimes that simulate those encountered in the interiors of planets. Fundamental Materials Science questions such as dislocation velocity regimes, the transition between thermally-activated and phonon drag regimes, the slip-twinning transition, the ultimate tensile strength of metals, the dislocation mechanisms of void growth are being answered through this powerful tool. In parallel with experiments, molecular dynamics simulations provide modeling and visualization at comparable strain rates (10 8 –10 10 s -1 ) and time durations (hundreds of picoseconds). Finally, this powerful synergy is illustrated in our past and current work, using representative face-centered cubic (fcc) copper, body-centered cubic (bcc) tantalum and diamond cubic silicon as model structures.

36 MATERIALS SCIENCE↗