Vesselin Yamakov

Dynamics of nanoscale grain-boundary decohesion in aluminum by molecular-dynamics simulation

The dynamics and energetics of intergranular crack growth along a flat grain boundary in aluminum is studied by a molecular-dynamics simulation model for crack propagation under steady-state conditions. Using the ability of the molecular-dynamics simulation to identify atoms involved in different atomistic mechanisms, it was possible to identify the energy contribution of different processes taking place during crack growth. The energy contributions were divided as: elastic energy—defined as the potential energy of the atoms in fcc crystallographic state; and plastically stored energy—the energy of stacking faults and twin boundaries; grain-boundary and surface energy. In addition, monitoring the amount of heat exchange with the molecular-dynamics thermostat gives the energy dissipated as heat in the system. The energetic analysis indicates that the majority of energy in a fast growing crack is dissipated as heat. This dissipation increases linearly at low speed, and faster than linear at speeds approaching 1/3 the Rayleigh wave speed when the crack tip becomes dynamically unstable producing periodic dislocation bursts until the crack is blunted.

Interatomic potential for the ternary Ni–Al–Co system and application to atomistic modeling of the B2–L10 martensitic transformation

Ni–Al–Co is a promising system for ferromagnetic shape memory applications. This paper reports on the development of a ternary embedded-atom potential for this system by fitting to experimental and first-principles data. Reasonably good agreement is achieved for physical properties between values predicted by the potential and values known from experiment and/or first-principles calculations. The potential reproduces basic features of the martensitic phase transformation from the B2-ordered high-temperature phase to a tetragonal CuAu-ordered low-temperature phase. The compositional and temperature ranges of this transformation and the martensite microstructure predicted by the potential compare well with existing experimental data. These results indicate that the proposed potential can be used for simulations of the shape memory effect in the Ni–Al–Co system.

Multifunctional Electroactive Nanocomposites Based on Piezoelectric Boron Nitride Nanotubes

Space exploration missions require sensors and devices capable of stable operation in harsh environments such as those that include high thermal fluctuation, atomic oxygen, and high-energy ionizing radiation. However, conventional or state-of-the-art electroactive materials like lead zirconate titanate, poly(vinylidene fluoride), and carbon nanotube (CNT)-doped polyimides have limitations on use in those extreme applications. Theoretical studies have shown that boron nitride nanotubes (BNNTs) have strength-to-weight ratios comparable to those of CNTs, excellent high-temperature stability (to 800 °C in air), large electroactive characteristics, and excellent neutron radiation shielding capability. In this study, we demonstrated the experimental electroactive characteristics of BNNTs in novel multifunctional electroactive nanocomposites. Upon application of an external electric field, the 2 wt % BNNT/polyimide composite was found to exhibit electroactive strain composed of a superposition of linear piezoelectric and nonlinear electrostrictive components. When the BNNTs were aligned by stretching the 2 wt % BNNT/polyimide composite, electroactive characteristics increased by about 460% compared to the nonstretched sample. An all-nanotube actuator consisting of a BNNT buckypaper layer between two single-walled carbon nanotube buckypaper electrode layers was found to have much larger electroactive properties. The additional neutron radiation shielding properties and ultraviolet/visible/near-infrared optical properties of the BNNT composites make them excellent candidates for use in the extreme environments of space missions.

Boron nitride nanotube: synthesis and applications

Scientists have predicted that carbon’s immediate neighbors on the periodic chart, boron and nitrogen, may also form perfect nanotubes, since the advent of carbon nanotubes (CNTs) in 1991. First proposed then synthesized by researchers at UC Berkeley in the mid 1990’s, the boron nitride nanotube (BNNT) has proven very difficult to make until now. Herein we provide an update on a catalyst-free method for synthesizing highly crystalline, small diameter BNNTs with a high aspect ratio using a high power laser under a high pressure and high temperature environment first discovered jointly by NASA/NIA/JSA. Progress in purification methods, dispersion studies, BNNT mat and composite formation, and modeling and diagnostics will also be presented. The white BNNTs offer extraordinary properties including neutron radiation shielding, piezoelectricity, thermal oxidative stability (> 800°C in air), mechanical strength, and toughness. The characteristics of the novel BNNTs and BNNT polymer composites and their potential applications are discussed.

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Time scale for point-defect equilibration in nanostructures

Molecular dynamics simulations of high-temperature annealing are performed on nanostructured materials enabling direct observation of vacancy emission from planar defects (i.e., grain boundaries and free surfaces) to populate the initially vacancy-free grain interiors on a subnanosecond time scale. We demonstrate a universal time-length scale correlation that governs these re-equilibration processes, suggesting that nanostructures are particularly stable against perturbations in their point-defect concentrations, caused for example by particle irradiation or temperature fluctuations.

Multiscale modeling of sensory properties of Co–Ni–Al shape memory particles embedded in an Al metal matrix

The concept of utilizing ferromagnetic shape memory alloys as embedded sensory particles in aluminum alloys for damage detection is discussed. When embedded in a material, a shape memory particle can undergo an acoustically detectable solid-state phase transformation when the local strain reaches a critical value. The emitted acoustic signal can be used for real-time damage detection. To study the transition behavior of the sensory particle inside a metal matrix under load, a simulation approach based on a coupled atomistic-continuum model is used. The simulation results indicate a strong dependence of the particle’s pseudoelastic response on its crystallographic orientation with respect to the loading direction. These results serve as a basis for understanding the efficacy and variability in the sensory particle transformation to detect damage processes.

Multiscale Modeling of Grain-Boundary Fracture: Cohesive Zone Models Parameterized from Atomistic Simulations

A multiscale modeling strategy is developed to study grain boundary fracture in polycrystalline aluminum. Atomistic simulation is used to model fundamental nanoscale deformation and fracture mechanisms and to develop a constitutive relationship for separation along a grain boundary interface. The nanoscale constitutive relationship is then parameterized within a cohesive zone model to represent variations in grain boundary properties. These variations arise from the presence of vacancies, interstitials, and other defects in addition to deviations in grain boundary angle from the baseline configuration considered in the molecular-dynamics simulation. The parameterized cohesive zone models are then used to model grain boundaries within finite element analyses of aluminum polycrystals.

Effects of microstructural inhomogeneity on dynamic grain growth during large-strain grain boundary diffusion-assisted plastic deformation

Mesoscale simulations of grain-boundary (GB) diffusion creep in which GB migration-induced static grain growth is suppressed were carried out based on the variational principle of dissipated power. Assuming that the boundaries exhibit no sliding resistance in response to shear stress, the variation of the normal-stress distribution and the diffusive fluxes along the grain boundaries during Coble creep were analysed. The effects of microstructural inhomogeneity, including topological and physical inhomogeneity (i.e. distributions in the grain sizes and GB diffusivities) were investigated. We find that because of the lack of GB migration as an accommodation process to relax the stress concentrations in the microstructure, a topologically inhomogeneous microstructure becomes unphysical at high strains (of typically between 50–100%). Consistent with earlier simulations by Pan and Cocks (1993 Comput. Mater. Sci. 1 95), we find that even in the absence of static grain growth an inhomogeneous microstructure exhibits dynamic grain growth induced by grain-switching induced grain-disappearance events. Our simulations also reveal that in a topologically inhomogeneous microstructure, the diffusive fluxes along any given GB can be in the same direction at both delimiting triple points; i.e. qualitatively different from a homogeneous (i.e. regular hexagonal) microstructure in which, according to Spingarn and Nix (1978 Acta Metall. 26 1389), these fluxes oppose each other.

A Continuum-Atomistic Analysis of Transgranular Crack Propagation in Aluminum

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain is used to study plastic processes at a crack tip in a single crystal of aluminum. The case of mode I loading is studied. A transition from deformation twinning to full dislocation emission from the crack tip is found when the crack plane is rotated around the [ 111 ] crystallographic axis. When the crack plane normal coincides with the [112] twinning direction, the crack propagates through a twinning mechanism. When the crack plane normal coincides with the [011] slip direction, the crack propagates through the emission of full dislocations. In intermediate orientations, a transition from full dislocation emission to twinning is found to occur with an increase in the stress intensity at the crack tip. This finding confirms the suggestion that the very high strain rates, inherently present in MD simulations, which produce higher stress intensities at the crack tip, over-predict the tendency for deformation twinning compared to experiments. The present study, therefore, aims to develop a more realistic and accurate predictive modeling of fracture processes.