Electromigration (EM) has been a serious reliability concern in microelectronics packaging for close to half a century now. Whenever the challenges of EM are overcome newer complications arise such as the demand for better performance due to increased miniaturization of semiconductor devices or the problems faced due to undesirable properties of lead-free solders. The motivation for the work is that there exists no fully computational modeling study on EM damage in lead-free solders (and also in lead-based solders). Modeling techniques such as one developed here can give new insights on effects of different grain features and offer high flexibility in varying parameters and study the corresponding effects. In this work, a new computational approach has been developed to study void nucleation and initial void growth in solders due to metal atom diffusion. It involves the creation of a 3D stochastic mesoscale model of the microstructure of a polycrystalline Tin structure. The next step was to identify regions of current crowding or ‘hot-spots’. This was done through solving a finite difference scheme on top of the 3D structure. The nucleation of voids due to atomic diffusion from the regions of current crowding was modeled by diffusion from the identified hot-spot through a rejection free kinetic Monte-Carlo scheme. This resulted in the net movement of atoms from the cathode to the anode. The above steps of identifying the hotspot and diffusing the atoms at the hot-spot were repeated and this lead to the initial growth of the void. This procedure was studied varying different grain parameters. In the future, the goal is to explore the effect of more grain parameters and consider other mechanisms of failure such as the formation of intermetallic compounds due to interstitial diffusion and dissolution of underbump metallurgy.

Nanolaminate composite materials consist of alternating layers of materials at the nanoscale (≤100 nm). Due to the nanometer scale thickness of their layers, these materials display unique and tailorable properties. This enables us to alter both mechanical attributes such as strength and wear properties, as well as functional characteristics such as biocompatibility, optical, and electronic properties. This dissertation focuses on understanding the mechanical behavior of the Al-SiC system. From a practical perspective, these materials exhibit a combination of high toughness and strength which is attractive for many applications. Scientifically, these materials are interesting due to the large elastic modulus mismatch between the layers. This, paired with the small layer thickness, allows a unique opportunity for scientists to study the plastic deformation of metals under extreme amounts of constraint.

Previous studies are limited in scope and a more diverse range of mechanical characterization is required to understand both the advantages and limitations of these materials. One of the major challenges with testing these materials is that they are only able to be made in thicknesses on the order of micrometers so the testing methods are limited to small volume techniques. This work makes use of both microscale testing techniques from the literature as well as novel methodologies. Using these techniques we are able to gain insight into aspects of the material’s mechanical behavior such as the effects of layer orientation, flaw dependent fracture, tension-compression asymmetry, fracture toughness as a function of layer thickness, and shear behavior as a function of layer thickness.

The objective of this research is to develop robust, accurate, and adaptive algorithms in the framework of the extended finite element method (XFEM) for fracture analysis of highly heterogeneous materials with complex internal geometries. A key contribution of this work is the creation of novel methods designed to automate the incorporation of high-resolution data, e.g. from X-ray tomography, that can be used to better interpret the enormous volume of data generated in modern in-situ experimental testing. Thus new algorithms were developed for automating analysis of complex microstructures characterized by segmented tomographic images.

A centrality-based geometry segmentation algorithm was developed to accurately identify discrete inclusions and particles in composite materials where limitations in imaging resolution leads to spurious connections between particles in close contact.To allow for this algorithm to successfully segment geometry independently of particle size and shape, a relative centrality metric was defined to allow for a threshold centrality criterion for removal of voxels that spuriously connect distinct geometries.

To automate incorporation of microstructural information from high-resolution images, two methods were developed that initialize signed distance fields on adaptively-refined finite element meshes. The first method utilizes a level set evolution equation that is directly solved on the finite element mesh through Galerkins method. The evolution equation is formulated to produce a signed distance field that matches geometry defined by a set of voxels segmented from tomographic images. The method achieves optimal convergence for the order of elements used. In a second approach, the fast marching method is employed to initialize a distance field on a uniform grid which is then projected by least squares onto a finite element mesh. This latter approach is shown to be superior in speed and accuracy.

Lastly, extended finite element method simulations are performed for the analysis of particle fracture in metal matrix composites with realistic particle geometries initialized from X-ray tomographic data. In the simulations, particles fracture probabilistically through a Weibull strength distribution. The model is verified through comparisons with the experimentally-measured stress-strain response of the material as well as analysis of the fracture. Further, simulations are then performed to analyze the effect of mesh sensitivity, the effect of fracture of particles on their neighbors, and the role of a particles shape on its fracture probability.

Al 7075 alloys are used in a variety of structural applications, such as aircraft wings, automotive components, fuselage, spacecraft, missiles, etc. The mechanical and corrosion behavior of these alloys are dependent on their microstructure and the environment. Therefore, a comprehensive study on microstructural characterization and stress-environment interaction is necessary. Traditionally, 2D techniques have been used to characterize microstructure, which are inaccurate and inadequate since the research has shown that the results obtained in the bulk are different from those obtained on the surface. There now exist several techniques in 3D, which can be used to characterize the microstructure. Al 7075 alloys contain second phase particles which can be classified as Fe-bearing inclusions, Si-bearing inclusions and precipitates. The variation in mechanical and corrosion properties of aluminum alloys has been attributed to the size, shape, distribution, corrosion properties and mechanical behavior of these precipitates and constituent particles. Therefore, in order to understand the performance of Al 7075 alloys, it is critical to investigate the size and distribution of inclusions and precipitates in the alloys along with their mechanical properties, such as Young's modulus, hardness and stress-strain behavior. X-ray tomography and FIB tomography were used to visualize and quantify the microstructure of constituent particles (inclusions) and precipitates, respectively. Microscale mechanical characterization techniques, such as nanoindentation and micropillar compression, were used to obtain mechanical properties of inclusions. Over the years, studies have used surface measurements to understand corrosion behavior of materials. More recently, in situ mechanical testing has become more attractive and advantageous, as it enables visualization and quantification of microstructural changes as a function of time (4D). In this study, in situ X-ray synchrotron tomography was used to study the SCC behavior of Al 7075 alloys in moisture and deionized water. Furthermore, experiments were performed in EXCO solution to study the effect of applied stress on exfoliation behavior in 3D. Contrary to 2D measurements made at the surface which suggest non-uniform crack growth rates, three dimensional measurements of the crack length led to a much more accurate measurement of crack growth rates.

Hydrogen embrittlement (HE) is a phenomenon that affects both the physical and chemical properties of several intrinsically ductile metals. Consequently, understanding the mechanisms behind HE has been of particular interest in both experimental and modeling research. Discrepancies between experimental observations and modeling results have led to various proposals for HE mechanisms. Therefore, to gain insights into HE mechanisms in iron, this dissertation aims to investigate several key issues involving HE such as: a) the incipient crack tip events; b) the cohesive strength of grain boundaries (GBs); c) the dislocation-GB interactions and d) the dislocation mobility.

The crack tip, which presents a preferential trap site for hydrogen segregation, was examined using atomistic methods and the continuum based Rice-Thompson criterion as sufficient concentration of hydrogen can alter the crack tip deformation mechanism. Results suggest that there is a plausible co-existence of the adsorption induced dislocation emission and hydrogen enhanced decohesion mechanisms. In the case of GB-hydrogen interaction, we observed that the segregation of hydrogen along the interface leads to a reduction in cohesive strength resulting in intergranular failure. A methodology was further developed to quantify the role of the GB structure on this behavior.

GBs play a fundamental role in determining the strengthening mechanisms acting as an impediment to the dislocation motion; however, the presence of an unsurmountable barrier for a dislocation can generate slip localization that could further lead to intergranular crack initiation. It was found that the presence of hydrogen increases the strain energy stored within the GB which could lead to a transition in failure mode. Finally, in the case of body centered cubic metals, understanding the complex screw dislocation motion is critical to the development of an accurate continuum description of the plastic behavior. Further, the presence of hydrogen has been shown to drastically alter the plastic deformation, but the precise role of hydrogen is still unclear. Thus, the role of hydrogen on the dislocation mobility was examined using density functional theory and atomistic simulations. Overall, this dissertation provides a novel atomic-scale understanding of the HE mechanism and development of multiscale tools for future endeavors.

Increased priority on the minimization of environmental impacts of conventional construction materials in recent years has motivated increased use of waste materials or bi-products such as fly ash, blast furnace slag with a view to reduce or eliminate the manufacturing/consumption of ordinary portland cement (OPC) which accounts for approximately 5-7% of global carbon dioxide emission. The current study explores, for the first time, the possibility of carbonating waste metallic iron powder to develop carbon-negative sustainable binder systems for concrete. The fundamental premise of this work is that metallic iron will react with aqueous CO2 under controlled conditions to form complex iron carbonates which have binding capabilities. The compressive and flexural strengths of the chosen iron-based binder systems increase with carbonation duration and the specimens carbonated for 4 days exhibit mechanical properties that are comparable to those of companion ordinary portland cement systems. The optimal mixture proportion and carbonation regime for this non-conventional sustainable binder is established based on the study of carbonation efficiency of a series of mixtures using thermogravimetric analysis. The pore- and micro-structural features of this novel binding material are also evaluated. The fracture response of this novel binder is evaluated using strain energy release rate and measurement of fracture process zone using digital image correlation (DIC). The iron-based binder system exhibits significantly higher strain energy release rates when compared to those of the OPC systems in both the unreinforced and glass fiber reinforced states. The iron-based binder also exhibits higher amount of area of fracture process zone due to its ability to undergo inelastic deformation facilitated by unreacted metallic iron particle inclusions in the microstructure that helps crack bridging /deflection. The intrinsic nano-mechanical properties of carbonate reaction product are explored using statistical nanoindentation technique coupled with a stochastic deconvolution algorithm. Effect of exposure to high temperature (up to 800°C) is also studied. Iron-based binder shows significantly higher residual flexural strength after exposure to high temperatures. Results of this comprehensive study establish the viability of this binder type for concrete as an environment-friendly and economical alternative to OPC.

For decades, microelectronics manufacturing has been concerned with failures related to electromigration phenomena in conductors experiencing high current densities. The influence of interconnect microstructure on device failures related to electromigration in BGA and flip chip solder interconnects has become a significant interest with reduced individual solder interconnect volumes. A survey indicates that x-ray computed micro-tomography (µXCT) is an emerging, novel means for characterizing the microstructures' role in governing electromigration failures. This work details the design and construction of a lab-scale µXCT system to characterize electromigration in the Sn-0.7Cu lead-free solder system by leveraging in situ imaging.

In order to enhance the attenuation contrast observed in multi-phase material systems, a modeling approach has been developed to predict settings for the controllable imaging parameters which yield relatively high detection rates over the range of x-ray energies for which maximum attenuation contrast is expected in the polychromatic x-ray imaging system. In order to develop this predictive tool, a model has been constructed for the Bremsstrahlung spectrum of an x-ray tube, and calculations for the detector's efficiency over the relevant range of x-ray energies have been made, and the product of emitted and detected spectra has been used to calculate the effective x-ray imaging spectrum. An approach has also been established for filtering `zinger' noise in x-ray radiographs, which has proven problematic at high x-ray energies used for solder imaging. The performance of this filter has been compared with a known existing method and the results indicate a significant increase in the accuracy of zinger filtered radiographs.

The obtained results indicate the conception of a powerful means for the study of failure causing processes in solder systems used as interconnects in microelectronic packaging devices. These results include the volumetric quantification of parameters which are indicative of both electromigration tolerance of solders and the dominant mechanisms for atomic migration in response to current stressing. This work is aimed to further the community's understanding of failure-causing electromigration processes in industrially relevant material systems for microelectronic interconnect applications and to advance the capability of available characterization techniques for their interrogation.

Concrete is the most widely used infrastructure material worldwide. Production of portland cement, the main binding component in concrete, has been shown to require significant energy and account for approximately 5-7% of global carbon dioxide production. The expected continued increased use of concrete over the coming decades indicates this is an ideal time to implement sustainable binder technologies. The current work aims to explore enhanced sustainability concretes, primarily in the context of limestone and flow. Aspects such as hydration kinetics, hydration product formation and pore structure add to the understanding of the strength development and potential durability characteristics of these binder systems. Two main strategies for enhancing this sustainability are explored in this work: (i) the use of high volume limestone in combination with other alternative cementitious materials to decrease the portland cement quantity in concrete and (ii) the use of geopolymers as the binder phase in concrete. The first phase of the work investigates the use of fine limestone as cement replacement from the perspective of hydration, strength development, and pore structure. The nature of the potential synergistic benefit of limestone and alumina will be explored. The second phase will focus on the rheological characterization of these materials in the fresh state, as well as a more general investigation of the rheological characterization of suspensions. The results of this work indicate several key ideas. (i) There is a potential synergistic benefit for strength, hydration, and pore structure by using alumina and in portland limestone cements, (ii) the limestone in these systems is shown to react to some extent, and fine limestone is shown to accelerate hydration, (iii) rheological characteristics of cementitious suspensions are complex, and strongly dependent on several key parameters including: the solid loading, interparticle forces, surface area of the particles present, particle size distribution of the particles, and rheological nature of the media in which the particles are suspended, and (iv) stress plateau method is proposed for the determination of rheological properties of concentrated suspensions, as it more accurately predicts apparent yield stress and is shown to correlate well with other viscoelastic properties of the suspensions.