Search NASA⌕ Search

Engineering topics

Mikhail Mendelev

Publications and source records attributed to Mikhail Mendelev.

Molecular Dynamics Simulations of Microstructural Effects on Austenite-Martensite Interfaces in NiTi

Formation and migration of austenite-martensite interfaces plays the key role in reversible martensitic transformations of shape memory alloys (SMAs). How these interfaces interact with the SMA microstructure is a primary determining factor in important functional properties such as hysteresis and transformation span. As such, successful microstructural engineering of SMAs requires in-depth knowledge of interface behavior. The rapid nature of martensitic transformations makes experimental observations of moving austenite-martensite interfaces challenging. Molecular dynamics (MD) simulation is a unique tool which can probe the atomic-scale details of austenite-martensite interfaces as they migrate through different microstructures. However, in focusing on the entire transformation process, including the nucleation of new phases, MD studies are usually performed so far from equilibrium that their relevance to experiment is questionable. Here, we demonstrate new MD techniques to generate energetically preferred austenite-martensite interfaces in NiTi under near-equilibrium conditions. The interfaces are semi-coherent, exhibiting a series of structural disconnections, and they can migrate rapidly through single crystals under only small thermodynamic driving forces. In contrast, when interfaces migrate in polycrystals, their motion is impeded by thermoelastic effects as well as changes in orientation relationships at grain boundaries. Microstructures which accumulate large amounts of elastic energy tend to release some fraction through irreversible, hysteresis-inducing mechanisms. We demonstrate that engineering microstructures with less constraints is a viable strategy to produce SMAs with reduced hysteresis and transformation span. Similar thermoelastic and hysteresis-inducing mechanisms also arise when austenite-martensite interfaces encounter precipitates and can be controlled by tuning characteristics of the precipitate distribution.

Gabriel Plummer↗

Molecular Dynamics Simulations of Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

The unique properties of shape memory alloys (SMAs) arise from a reversible martensitic transformation. The nucleation and migration of austenite-martensite interfaces are the key to understanding the SMA properties. Molecular dynamics (MD) simulations can provide important atomic-scale information about these aspects of the transformation, but their time scales prevent the interface formation under near-equilibrium conditions relevant to experiment. We present a new MD methodology which allows for the natural formation of energetically preferred austenite-martensite interfaces under near-equilibrium conditions. Our simulation demonstrates that the interfaces in NiTi are semi-coherent, composed of a series of terrace planes and structural disconnections, and they migrate rapidly through single crystals with only a small thermodynamic driving force. In bi-crystals and polycrystals, the migration of these same interfaces is significantly impeded by grain boundaries and stored elastic energy. This behavior can result in SMA hysteresis via several mechanisms associated with nucleation and non-elastic strain accommodation.

Gabriel Plummer↗

Twin nucleation and growth mechanism in Ni-based superalloys

While micro-twinning is the dominant creep deformation mechanism in Ni-based superalloys at temperatures above 700 C, many aspects of twin nucleation and growth remain unexplored. Kolbe mechanism for micro-twinning, based on thermally activated reordering, is probably the only concept currently widely accepted by the scientific community. We propose a qualitatively different mechanism for nucleation and growth of twins. The mechanism can be briefly described as follows. Penetration of a   precipitate by two 1/2(110) edge dislocations travelling on adjacent {111} glide planes triggers nucleation (at the interface of the precipitate and the matrix) and emission of Shockley partial of screw character on the glide plane of the edge dislocation, which entered the precipitate first (generating trailing high-energy anti-phase boundary (APB)). Propagation of this Shockley partial into the precipitate converts the APB into super intrinsic stacking fault (SISF). The recurring arrival of additional edge dislocations on the glide planes adjacent to the configuration described above leads to formation of super extrinsic stacking fault SESF, subsequent micro-twin formation, and growth of the twinned region. We demonstrate the proposed mechanism via molecular dynamics simulations.

Ni-based superalloys↗

Effect of Alloying Additions on Twinning in Ni-Based Superalloys

Micro-twinning is the dominant creep deformation mechanism in Ni-based superalloys at temperatures above 700 C. We use atomistic simulations to study two mechanisms of twin nucleation and growth that are characterized by qualitatively different rate limiting processes. In case of the mechanism proposed by Kolbe, the rate limiting process is diffusion-mediated atomic reshuffling. In case of the other mechanism, we proposed recently, the rate limiting process is nucleation of Shockley partial dislocation. We demonstrate the effects of alloying additions on functionality of these mechanisms.

Ni-based superalloys↗

Influence of Microstructural Features on Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

Shape memory alloys (SMAs) exhibit several unique thermomechanical properties due to a reversible martensitic phase transformation. A key aspect of this transformation is the austenite-martensite interface. Here we present a molecular dynamics (MD) simulation methodology to study the atomic-scale features of austenite-martensite interface migration under near-equilibrium conditions. In single crystals, the interfaces migrate rapidly with only a small thermodynamic driving force. In polycrystals, however, interface migration is significantly impeded due to the change in orientation relationship at grain boundaries and the stored elastic energy resulting from microstructural constraints. This behavior can be linked to several mechanisms associated with transformation width and hysteresis in SMAs, properties of great importance for applications involving actuation. Additionally, we will present preliminary MD simulation results on the influence of precipitates on the formation and migration of austenite-martensite interfaces.

Gabriel Plummer↗

Solid-state Precipitation in Molecular Dynamics: KMC-MD Hybrid Simulations

A critical drawback of molecular dynamics (MD) simulations when modeling physical properties is the limitation to short timescales. Solid state precipitation is one behavior that is not significantly observable at timescales accessible to MD, and thus has been difficult to simulate at an MD scale. Taking inspiration from Kinetic Monte-Carlo (KMC) methods for simulating diffusion on a rigid lattice, we have developed a hybrid KMC-MD approach to enable simulation of solid-state precipitation during MD simulation implemented in LAMMPS. This method incorporates the atomic neighbor hopping behavior utilized by KMC simulations into a MD simulation with a variable local lattice configuration. By accelerating nearest-neighbor atomic hopping with this technique, solid-state precipitation can be observed at an MD timescale. Utilizing theory based on persistent-embryo methods for liquid-solid precipitation in MD we demonstrate the calculation of critical nucleus size for the gamma’ precipitation out of gamma phase in the Ni-Al system.

Molecular Dynamics↗

Effects of Cr and Nb Alloying Additions on High Temperature Creep Deformation of Ni-Based Superalloys

Micro-twinning is the major creep deformation mechanism in Ni-based superalloys at temperatures above 700 C°. Recent experiments suggest that superlattice stacking faults in phase may serve as the precursors to twin formation. Segregation of alloying elements to these precursors may have a significant effect on formation and extension of micro-twins. Using atomistic modeling we investigate and explain the effects of Cr and Nb alloying additions on these processes. We consider the experimentally obtained creep curves for several Ni-based superalloys with different elemental compositions and relate the observed trends to the simulation results presented in this study. Our results help to rationalize the experimentally observed puzzling effects of elemental composition of the alloy on creep resistance.

Ni-based superalloys↗

Effect of Cr Twinning in Ni-Based Superalloys

Micro-twinning is the major creep deformation mechanism in Ni-based superalloys at temperatures above 700 °C. Recent experiments indicate that formation and extension of micro-twins is reliant on the segregation of gamma phase formers (e.g., Co and Cr) along the twin interface in gamma prime phase. In this study we use molecular dynamics simulation to investigate the effect of site preference of Cr in Ni 3 Al gamma prime phase on deformation twinning. We demonstrate that, depending on the site preference (Cr on Ni sites, or Cr on Al sites), two drastically different deformation behaviors can be expected.

Ni-based superalloys↗

Special Grain Boundaries in NiTi Shape Memory Alloys as Sites for Preferential Martensite Nucleation

Shape memory alloys (SMAs) exhibit unique thermomechanical properties due to a reversible martensitic transformation which can be controlled via chemical composition and microstructural features. Among the latter, grain boundaries (GBs) are key in determining how the transformation nucleates and propagates. We performed molecular dynamics simulations to examine the roles of a few special GBs in the austenite phase of NiTi. The GBs can act as preferential host sites for martensite nuclei which can substantially lower the nucleation barrier for transformation and thereby result in reduced thermal hysteresis, an important SMA property for cyclic actuation applications. Free energy calculations show that a characteristic of these GBs is a negative entropy which drops sharply close to the transformation temperature, an anomalous behavior for GBs with fixed composition which is a direct result of the martensitic transformation. We discuss implications of these results with respect to SMA processing techniques for achieving improved properties.

Gabriel Plummer↗

Effects of Alloying Elements on Twinning in Ni-Based Superalloys

Micro-twinning is the major creep deformation mechanism in Ni-based superalloys at temperatures above 700 °C. Recent experiments suggest that superlattice stacking faults in y' phase may serve as the precursors to twin formation. Segregation of alloying elements to these precursors may have a significant effect on formation and extension of micro-twins. Using atomistic modeling we investigate and explain the effects of Nb and Cr alloying additions on these processes. The simulation shows that Nb increases the creep resistance which is mostly associated with impeding the reordering of the high energy double complex stacking fault. Cr, on the other hand, promotes twin growth, degrading the high temperature creep properties. These results can help to understand the effects of elemental composition of the alloy on creep resistance.

Ni-based superalloys↗

Scalable Variable Charge Molecular Dynamics Simulations of Metal-Oxide Systems

Interfaces between metals and oxides are an important feature in many technologically relevant materials, e.g., oxidation of metal surfaces, oxide-dispersion strengthened (ODS) alloys, dielectric components, and thermal barrier coatings among others. Experimental studies of such interfaces are challenging since the majority are buried within the bulk, making computational modeling an attractive alternative. Molecular dynamics (MD) simulations operate at the length scales relevant to many interface-mediated mechanisms, but the requisite interatomic potentials for metal-oxide systems require computationally expensive variable charge schemes to account for the disparate bonding types, thus often limiting their effectiveness. Here we introduce several improvements to the charge transfer interatomic potential (CTIP) model which enable greater computational efficiency for large scale MD simulations. Then, using a new CTIP parametrization for the Ni-O system, we demonstrate its capabilities to capture critical atomic scale mechanisms associated with metal-oxide interfaces. Long time scale simulations (>10 ns) are used to investigate high temperature oxidation and oxide precipitation from the melt, and large length scale simulations (> 1 million atoms) are used to study the interaction of dislocations with oxide particles. We have implemented the new CTIP model in the widely used, open-source MD code LAMMPS.

Gabriel Plummer↗