Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Despite the technological importance of body-centered cubic (BCC) metals, models of their plastic deformation are less common than those of face-centered cubic (FCC) metals, due in part to the complexity of slip in BCC crystals caused by the thermal activation of screw dislocation motion. This paper presents a physically based crystal plasticity model that incorporates atomistic models and experimental measurements of the thermally activated nature of screw dislocation motion. This model, therefore, reproduces the temperature, stress, and strain rate dependence of flow in BCC metals in a simple formulation that will allow for large, grain-scale simulations. Furthermore, the results illustrate the importance of correctly representing mechanistic transitions in materials with high lattice friction. © 2012 Elsevier Ltd. All rights reserved.

During large-strain plastic deformation, subgrain structures typically develop within the grains. At large enough equivalent strains above, say 0.5, recrystallization occurs via abnormal coarsening of the subgrain structure or abnormal (sub-) grain growth (AsGG). The fraction of subgrains that develop into new, recrystallized grains has been quantified as a function of texture spread (Grain Reference Orientation Deviation) using Monte Carlo simulation. When this fraction is combined with the known monotonic increase in mean misorientation with strain, the recrystallized grain size can be predicted as a function of von Mises strain. The prediction is in good agreement with experimental results drawn from the literature. © (2012) Trans Tech Publications, Switzerland.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In this project, we performed a preliminary set of sintering experiments to examine nanocrystal-enabled diffusion bonding (NEDB) in Ag-on-Ag and Cu-on-Cu using Ag nanoparticles. The experimental test matrix included the effects of material system, temperature, pressure, and particle size. The nanoparticle compacts were bonded between plates using a customized hot press, tested in shear, and examined post mortem using microscopy techniques. NEDB was found to be a feasible mechanism for low-temperature, low-pressure, solid-state bonding of like materials, creating bonded interfaces that were able to support substantial loads. The maximum supported shear strength varied substantially within sample cohorts due to variation in bonded area; however, systematic variation with fabrication conditions was also observed. Mesoscale sintering simulations were performed in order to understand whether sintering models can aid in understanding the NEDB process. A pressure-assisted sintering model was incorporated into the SPPARKS kinetic Monte Carlo sintering code. Results reproduce most of the qualitative behavior observed in experiments, indicating that simulation can augment experiments during the development of the NEDB process. Because NEDB offers a promising route to low-temperature, low-pressure, solid-state bonding, we recommend further research and development with a goal of devising new NEDB bonding processes to support Sandia's customers.

Legislated requirements and industry standards are replacing eutectic lead-tin (Pb-Sn) solders with lead-free (Pb-free) solders in future component designs and in replacements and retrofits. Since Pb-free solders have not yet seen service for long periods, their long-term behavior is poorly characterized. Because understanding the reliability of Pb-free solders is critical to supporting the next generation of circuit board designs, it is imperative that we develop, validate and exercise a solder lifetime model that can capture the thermomechanical response of Pb-free solder joints in stockpile components. To this end, an ASC Level 2 milestone was identified for fiscal year 2010: Milestone 3605: Utilize experimentally validated constitutive model for lead-free solder to simulate aging and reliability of solder joints in stockpile components. This report documents the completion of this milestone, including evidence that the milestone completion criteria were met and a summary of the milestone Program Review.

One of the tenets of nanotechnology is that the electrical/optical/chemical/biological properties of a material may be changed profoundly when the material is reduced to sufficiently small dimensions - and we can exploit these new properties to achieve novel or greatly improved material's performance. However, there may be mechanical or thermodynamic driving forces that hinder the synthesis of the structure, impair the stability of the structure, or reduce the intended performance of the structure. Examples of these phenomena include de-wetting of films due to high surface tension, thermally-driven instability of nano-grain structure, and defect-related internal dissipation. If we have fundamental knowledge of the mechanical processes at small length scales, we can exploit these new properties to achieve robust nanodevices. To state it simply, the goal of this program is the fundamental understanding of the mechanical properties of materials at small length scales. The research embodied by this program lies at the heart of modern materials science with a guiding focus on structure-property relationships. We have divided this program into three Tasks, which are summarized: (1) Mechanics of Nanostructured Materials (PI Blythe Clark). This task aims to develop a fundamental understanding of the mechanical properties and thermal stability of nanostructured metals, and of the relationship between nano/microstructure and bulk mechanical behavior through a combination of special materials synthesis methods, nanoindentation coupled with finite-element modeling, detailed electron microscopic characterization, and in-situ transmission electron microscopy experiments. (2) Theory of Microstructures and Ensemble Controlled Deformation (PI Elizabeth A. Holm). The goal of this Task is to combine experiment, modeling, and simulation to construct, analyze, and utilize three-dimensional (3D) polycrystalline nanostructures. These full 3D models are critical for elucidating the complete structural geometry, topology, and arrangements that control experimentally-observed phenomena, such as abnormal grain growth, grain rotation, and internal dissipation measured in nanocrystalline metal. (3) Mechanics and Dynamics of Nanostructured and Nanoscale Materials (PI John P. Sullivan). The objective of this Task is to develop atomic-scale understanding of dynamic processes including internal dissipation in nanoscale and nanostructured metals, and phonon transport and boundary scattering in nanoscale structures via internal friction measurements.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The kinetic Monte Carlo method and its variants are powerful tools for modeling materials at the mesoscale, meaning at length and time scales in between the atomic and continuum. We have completed a 3 year LDRD project with the goal of developing a parallel kinetic Monte Carlo capability and applying it to materials modeling problems of interest to Sandia. In this report we give an overview of the methods and algorithms developed, and describe our new open-source code called SPPARKS, for Stochastic Parallel PARticle Kinetic Simulator. We also highlight the development of several Monte Carlo models in SPPARKS for specific materials modeling applications, including grain growth, bubble formation, diffusion in nanoporous materials, defect formation in erbium hydrides, and surface growth and evolution.

Interfaces are a critical determinant of the full range of materials properties, especially at the nanoscale. Computational and experimental methods developed a comprehensive understanding of nanograin evolution based on a fundamental understanding of internal interfaces in nanocrystalline nickel. It has recently been shown that nanocrystals with a bi-modal grain-size distribution possess a unique combination of high-strength, ductility and wear-resistance. We performed a combined experimental and theoretical investigation of the structure and motion of internal interfaces in nanograined metal and the resulting grain evolution. The properties of grain boundaries are computed for an unprecedented range of boundaries. The presence of roughening transitions in grain boundaries is explored and related to dramatic changes in boundary mobility. Experimental observations show that abnormal grain growth in nanograined materials is unlike conventional scale material in both the level of defects and the formation of unfavored phases. Molecular dynamics simulations address the origins of some of these phenomena.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Many conventional continuum approaches to solid mechanics do not address the size sensitivity of deformation to microstructural features like grain boundaries, and are therefore unable to capture much of the experimentally observed behavior of polycrystal deformation. We propose a non-local crystal plasticity model, in which the geometrically necessary dislocation (GND) density is calculated using a non-local integral approach. The model is based on augmented FeFp kinematics, which account for the initial microstructure (primarily grain boundaries) present in the polycrystal. With the augmented kinematics, the initial GND and the evolving GND state are determined in a consistent manner. The expanded kinematics and the non-local crystal plasticity model are used to simulate the tensile behavior in copper polycrystals with different grain sizes ranging from 14 μm to 244 μm. The simulation results show a grain size dependence on the polycrystal's yield strength, which are in good agreement with the experimental data. © 2007 Elsevier Ltd. All rights reserved.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Like other interfaces, equilibrium grain boundaries are smooth at low temperature and rough at high temperature; however, little attention has been paid to roughening except for faceting boundaries. Using molecular dynamics simulations of face-centered cubic Ni, we studied two closely related grain boundaries with different boundary planes. In spite of their similarity, their boundary roughening temperatures differ by several hundred degrees, and boundary mobility is much larger above the roughening temperature. This has important implications for microstructural development during metallurgical processes.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We describe two geometric structures, the shortest path and the minimum cut, and show that these structures emerge at special threshold points in the highly non-linear electrical response of complex networks. Algorithms which find the shortest path and the minimum cut directly as well as methods for finding the full non-linear response of complex networks are outlined. Scaling laws for the behavior of the shortest path and minimum cut in random networks are then surveyed. Finally, applications of the shortest path and minimum cut to grain boundary controlled polycrystalline materials are elucidated. © 2005 Elsevier Ltd. All rights reserved.

Abstract not provided.

Abstract not provided.

Abstract not provided.

As current experimental and simulation methods cannot determine the mobility of flat boundaries across the large misorientation phase space, we have developed a computational method for imposing an artificial driving force on boundaries. In a molecular dynamics simulation, this allows us to go beyond the inherent timescale restrictions of the technique and induce non-negligible motion in flat boundaries of arbitrary misorientation. For different series of symmetric boundaries, we find both expected and unexpected results. In general, mobility increases as the grain boundary plane deviates from (111), but high-coincidence and low-angle boundaries represent special cases. These results agree with and enrich experimental observations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Percolation theory is now standard in the analysis of polycrystalline materials where the grain boundaries can be divided into two distinct classes, namely 'good' boundaries that have favorable properties and 'bad' boundaries that seriously degrade the material performance. Grain-boundary engineering (GBE) strives to improve material behavior by engineering the volume fraction c and arrangement of good grain boundaries. Two key percolative processes in GBE materials are the onset of percolation of a strongly connected aggregate of grains, and the onset of a connected path of weak grain boundaries. Using realistic polycrystalline microstructures, we find that in two dimensions the threshold for strong aggregate percolation c{sub SAP} and the threshold for weak boundary percolation c{sub WBP} are equivalent and have the value c{sub SAP} = c{sub WBP} = 0.38(1), which is slightly higher than the threshold found for regular hexagonal grain structures, c{sub RH} = 2 sin({pi}/18) = 0.347. In three dimensions strong aggregate percolation and weak boundary percolation occur at different locations and we find c{sub SAP} = 0.12(3) and c{sub WBP} = 0.77(3). The critical current in high T{sub c} materials and the cohesive energy in structural systems are related to the critical manifold problem in statistical physics. We develop a theory of critical manifolds in GBE materials, which has three distinct regimes: (1) low concentrations, where random manifold theory applies, (2) critical concentrations where percolative scaling theory applies, and (3) high concentrations, c > c{sub SAP}, where the theory of periodic elastic media applies. Regime (3) is perhaps most important practically and is characterized by a critical length L{sub c}, which is the size of cleavage regions on the critical manifold. In the limit of high contrast {open_square} {yields} 0, we find that in two dimensions L{sub c} {proportional_to} gc/(1-c), while in three dimensions L{sub c} {proportional_to} g exp[b{sub 0}c/(1-c)]/[c(1-c)]{sup 1/2}, where g is the average grain size, {open_square} is the ratio of the bonding energy of the weak boundaries to that of the strong boundaries, and b{sub 0} is a constant which is of order 1. Many of the properties of GBE materials can be related to L{sub c}, which diverges algebraically on approach to c=1 in two dimensions, but diverges exponentially in that limit in three dimensions. We emphasize that GBE percolation processes and critical manifold behavior are very different in two dimensions as compared to three dimensions. For this reason, the use of two dimensional models to understand the behavior of bulk GBE materials can be misleading.

In the present work the authors describe the adaptation of a standard SEM into a flexible microjoining tool. The system incorporates exceptional control of energy input and its location, environmental cleanliness, part manipulation and especially, part imaging. Beam energetics, modeling of thermal flow in a simple geometry, significant effects of surface energy on molten pools and beam size characterization are treated. Examples of small to micro fusion welds and molten zones produced in a variety of materials (Ni, tool steel, Tophet C, Si) and sizes are given. Future directions are also suggested.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Biological systems create proteins that perform tasks more efficiently and precisely than conventional chemicals. For example, many plants and animals produce proteins to control the freezing of water. Biological antifreeze proteins (AFPs) inhibit the solidification process, even below the freezing point. These molecules bond to specific sites at the ice/water interface and are theorized to suppress solidification chemically or geometrically. In this project, we investigated the theoretical and experimental data on AFPs and performed analyses to understand the unique physics of AFPs. The experimental literature was analyzed to determine chemical mechanisms and effects of protein bonding at ice surfaces, specifically thermodynamic freezing point depression, suppression of ice nucleation, decrease in dendrite growth kinetics, solute drag on the moving solid/liquid interface, and stearic pinning of the ice interface. Stearic pinning was found to be the most likely candidate to explain experimental results, including freezing point depression, growth morphologies, and thermal hysteresis. A new stearic pinning model was developed and applied to AFPs, with excellent quantitative results. Understanding biological antifreeze mechanisms could enable important medical and engineering applications, but considerable future work will be necessary.

Abstract not provided.

Abnormal subgrain growth has been proposed as the nucleation mechanism for recrystallization. To test this hypothesis, Monte Carlo Potts model simulations of subgrain growth were performed on single-phase, strain-free subgrain structures with experimentally validated microstructure, texture, boundary character, and boundary properties. Results indicate that abnormal growth events emerge spontaneously during evolution in such systems, and abnormal subgrains behave as predicted by mean field theory. An analysis predicts the frequency of abnormal growth events as a function of local neighborhood and the boundary misorientation distribution. A recrystallization model is derived based on the abnormal subgrain growth analysis. Using data for aluminum subgrain structures, the model predicts reasonable recrystallized grain sizes as a function of von Mises strain. The extension of these results to abnormal grain growth is discussed. © 2003 Published by Elsevier Ltd on behalf of Acta Materialia Inc.

The purpose of microstructural control is to optimize materials properties. To that end, they have developed sophisticated and successful computational models of both microstructural evolution and mechanical response. However, coupling these models to quantitatively predict the properties of a given microstructure poses a challenge. This problem arises because most continuum response models, such as finite element, finite volume, or material point methods, do not incorporate a real length scale. Thus, two self-similar polycrystals have identical mechanical properties regardless of grain size, in conflict with theory and observations. In this project, they took a tiered risk approach to incorporate microstructure and its resultant length scales in mechanical response simulations. Techniques considered include low-risk, low-benefit methods, as well as higher-payoff, higher-risk methods. Methods studied include a constitutive response model with a local length-scale parameter, a power-law hardening rate gradient near grain boundaries, a local Voce hardening law, and strain-gradient polycrystal plasticity. These techniques were validated on a variety of systems for which theoretical analyses and/or experimental data exist. The results may be used to generate improved constitutive models that explicitly depend upon microstructure and to provide insight into microstructural deformation and failure processes. Furthermore, because mechanical state drives microstructural evolution, a strain-enhanced grain growth model was coupled with the mechanical response simulations. The coupled model predicts both properties as a function of microstructure and microstructural development as a function of processing conditions.

The Microsystems Subgrid Physics project is intended to address gaps between developing high-performance modeling and simulation capabilities and microdomain specific physics. The initial effort has focused on incorporating electrostatic excitations, adhesive surface interactions, and scale dependent material and thermal properties into existing modeling capabilities. Developments related to each of these efforts are summarized, and sample applications are presented. While detailed models of the relevant physics are still being developed, a general modeling framework is emerging that can be extended to incorporate evolving material and surface interaction modules.

A kinetic, three-dimensional Monte Carlo model for simulating grain growth in the presence of mobile pores is presented. The model was used to study grain growth and pore migration by surface diffusion in an idealized geometry that ensures constant driving force for grain growth. The driving forces, pore size, and pore mobilities were varied to study their effects on grain-boundary mobility and grain growth. The simulations captured a variety of complex behaviors, including reduced grain-boundary velocity due to pore drag that has been predicted by analytical theories. The model is capable of treating far more complex geometries, including polycrystals. We present the capabilities of this model and discuss its limitations.

Abstract not provided.

Computational materials simulations have traditionally focused on individual phenomena: grain growth, crack propagation, plastic flow, etc. However, real materials behavior results from a complex interplay between phenomena. In this project, the authors explored methods for coupling mesoscale simulations of microstructural evolution and micromechanical response. In one case, massively parallel (MP) simulations for grain evolution and microcracking in alumina stronglink materials were dynamically coupled. In the other, codes for domain coarsening and plastic deformation in CuSi braze alloys were iteratively linked. this program provided the first comparison of two promising ways to integrate mesoscale computer codes. Coupled microstructural/micromechanical codes were applied to experimentally observed microstructures for the first time. In addition to the coupled codes, this project developed a suite of new computational capabilities (PARGRAIN, GLAD, OOF, MPM, polycrystal plasticity, front tracking). The problem of plasticity length scale in continuum calculations was recognized and a solution strategy was developed. The simulations were experimentally validated on stockpile materials.