Uranium Dioxide (UO2) is a significant nuclear fission fuel, which is widely used

in nuclear reactors. Understanding the influence of microstructure on thermo-mechanical behavior of UO2 is extremely important to predict its performance. In particular, evaluating mechanical properties, such as elasticity, plasticity and creep at sub-grain length scales is key to developing this understanding as well as building multi-scale models of fuel behavior with predicting capabilities. In this work, modeling techniques were developed to study effects of microstructure on Young’s modulus, which was selected as a key representative property that affects overall mechanical behavior, using experimental data obtained from micro-cantilever bending testing as benchmarks. Beam theory was firstly introduced to calculate Young's modulus of UO2 from the experimental data and then three-dimensional finite element models of the micro-cantilever beams were constructed to simulate bending tests in UO2 at room temperature. The influence of the pore distribution was studied to explain the discrepancy between predicted values and experimental results. Results indicate that results of tests are significantly affected by porosity given that both pore size and spacing in the samples are of the order of the micro-beam dimensions. Microstructure reconstruction was conducted with images collected from three-dimensional serial sectioning using focused ion beam (FIB) and electron backscattering diffraction (EBSD) and pore clusters were placed at different locations along the length of the beam. Results indicate that the presence of pore clusters close to the substrate, i.e., the clamp of the micro-cantilever beam, has the strongest effect on load-deflection behavior, leading to a reduction of stiffness that is the largest for any location of the pore cluster. Furthermore, it was also found from both numerical and i

analytical models that pore clusters located towards the middle of the span and close to the end of the beam only have a very small effect on the load-deflection behavior, and it is concluded that better estimates of Young's modulus can be obtained from micro- cantilever experiments by using microstructurally explicit models that account for porosity in about one half of the beam length close to the clamp. This, in turn, provides an avenue to simplify micro-scale experiments and their analysis.

Shock loading is a complex phenomenon that can lead to failure mechanisms such as strain localization, void nucleation and growth, and eventually spall fracture. The length scale of damage with respect to that of the surrounding microstructure has proven to be a key aspect in determining sites of failure initiation. Studying incipient stages of spall damage is of paramount importance to accurately determine initiation sites in the material microstructure where damage will nucleate and grow and to formulate continuum models that account for the variability of the damage process due to microstructural heterogeneity, which is the focus of this research. Shock loading experiments were conducted via flyer-plate impact tests for pressures of 2-6 GPa and strain rates of 105/s on copper polycrystals of varying thermomechanical processing conditions. Serial cross sectioning of recovered target disks was performed along with electron microscopy, electron backscattering diffraction (EBSD), focused ion beam (FIB) milling, and 3-D X-ray tomogrpahy (XRT) to gain 2-D and 3-D information on the spall plane and surrounding microstructure. Statistics on grain boundaries (GB) containing damage were obtained from 2-D data and GBs of misorientations 25° and 50° were found to have the highest probability to contain damage in as-received (AR), heat treated (HT), and fully recrystallized (FR) microstructures, while {111} Σ3 GBs were globally strong. The AR microstructure’s probability peak was the most pronounced indicating GB strength is the dominant factor for damage nucleation. 3-D XRT data was used to digitally render the spall planes of the AR, HT, and FR microstructures. From shape fitting the voids to ellipsoids, it was found that the AR microstructure contained greater than 55% intergranular damage, whereas the HT and FR microstructures contained predominantly transgranular and coalesced damage modes, respectively. 3-D reconstructions of large volume damage sites in shocked Cu multicrystals showed preference for damage nucleation at GBs between adjacent grains of a high Taylor factor mismatches as well as an angle between the shock direction and the GB physical normal of ~30°-45°. 3-D FIB sectioning of individual voids led to the discovery of uniform plastic zones ~25-50% the size of the void diameter and plastic deformation directions were characterized via local average misorientation maps. Incipient transgranular voids revealed from the sectioning process were present in grains of high Taylor factors along the shock direction, which is expected as materials with a low Taylor factor along the shock direction are susceptible to growth due their accomodation of plastic deformation. Fabrication of square waves using photolithography and chemical etching was developed to study the nature of plasticity at GBs away from the spall plane. Grains oriented close to <0 1 1> had half the residual amplitudes than grains oriented close to <0 0 1>.

The United States Department of Energy (DOE) has always held the safety and reliability of the nation's nuclear reactor fleet as a top priority. Continual improvements and advancements in nuclear fuels have been instrumental in maximizing energy generation from nuclear power plants and minimizing waste. One aspect of the DOE Fuel Cycle Research and Development Advanced Fuels Campaign is to improve the mechanical properties of uranium dioxide (UO2) for nuclear fuel applications.

In an effort to improve the performance of UO2, by increasing the fracture toughness and ductility, small quantities of oxide materials have been added to samples to act as dopants. The different dopants used in this study are: titanium dioxide, yttrium oxide, aluminum oxide, silicon dioxide, and chromium oxide. The effects of the individual dopants and some dopant combinations on the microstructure and mechanical properties are determined using indentation fracture experiments in tandem with scanning electron microscopy. Indentation fracture experiments are carried out at room temperature and at temperatures between 450 °C and 1160 °C.

The results of this work find that doping with aluminosilicate produces the largest favorable change in the mechanical properties of UO2. This sample exhibits an increase in fracture toughness at room temperature without showing a change in yield strength at elevated temperatures. The results also show that doping with Al2O3 and TiO2 produce stronger samples and it is hypothesized that this is a result of the sample containing dopant-rich secondary phase particles.

Fission products in nuclear fuel pellets can affect fuel performance as they change the fuel chemistry and structure. The behavior of the fission products and their release mechanisms are important to the operation of a power reactor. Research has shown that fission product release can occur through grain boundary (GB) at low burnups. Early fission gas release models, which assumed spherical grains with no effect of GB diffusion, did not capture the early stage of the release behavior well. In order to understand the phenomenon at low burnup and how it leads to the later release mechanism, a microstructurally explicit model is needed. This dissertation conducted finite element simulations of the transport behavior using 3-D microstructurally explicit models. It looks into the effects of GB character, with emphases on conditions that can lead to enhanced effective diffusion. Moreover, the relationship between temperature and fission product transport is coupled to reflect the high temperature environment.

The modeling work began with 3-D microstructure reconstruction for three uranium oxide samples with different oxygen stoichiometry: UO2.00 UO2.06 and UO2.14. The 3-D models were created based on the real microstructure of depleted UO2 samples characterized by Electron Backscattering Diffraction (EBSD) combined with serial sectioning. Mathematical equations on fission gas diffusion and heat conduction were studied and derived to simulate the fission gas transport under GB effect. Verification models showed that 2-D elements can be used to model GBs to reduce the number of elements. The effect of each variable, including fuel stoichiometry, temperature, GB diffusion, triple junction diffusion and GB thermal resistance, is verified, and they are coupled in multi-physics simulations to study the transport of fission gas at different radial location of a fuel pellet. It was demonstrated that the microstructural model can be used to incorporate the effect of different physics to study fission gas transport. The results suggested that the GB effect is the most significant at the edge of fuel pellet where the temperature is the lowest. In the high temperature region, the increase in bulk diffusivity due to excess oxygen diminished the effect of GB diffusion.

Commercially pure (CP) and extra low interstitial (ELI) grade Ti-alloys present excellent corrosion resistance, lightweight, and formability making them attractive materials for expanded use in transportation and medical applications. However, the strength and toughness of CP titanium are affected by relatively small variations in their impurity/solute content (IC), e.g., O, Al, and V. This increase in strength is due to the fact that the solute either increases the critical stress required for the prismatic slip systems ({10-10}<1-210>) or activates another slip system ((0001)<11-20>, {10-11}<11-20>). In particular, solute additions such as O can effectively strengthen the alloy but with an attendant loss in ductility by changing the behavior from wavy (cross slip) to planar nature. In order to understand the underlying behavior of strengthening by solutes, it is important to understand the atomic scale mechanism. This dissertation aims to address this knowledge gap through a synergistic combination of density functional theory (DFT) and molecular dynamics. Further, due to the long-range strain fields of the dislocations and the periodicity of the DFT simulation cells, it is difficult to apply ab initio simulations to study the dislocation core structure. To alleviate this issue we developed a multiscale quantum mechanics/molecular mechanics approach (QM/MM) to study the dislocation core. We use the developed QM/MM method to study the pipe diffusion along a prismatic edge dislocation core. Complementary to the atomistic simulations, the Semi-discrete Variational Peierls-Nabarro model (SVPN) was also used to analyze the dislocation core structure and mobility. The chemical interaction between the solute/impurity and the dislocation core is captured by the so-called generalized stacking fault energy (GSFE) surface which was determined from DFT-VASP calculations. By taking the chemical interaction into consideration the SVPN model can predict the dislocation core structure and mobility in the presence and absence of the solute/impurity and thus reveal the effect of impurity/solute on the softening/hardening behavior in alpha-Ti. Finally, to study the interaction of the dislocation core with other planar defects such as grain boundaries (GB), we develop an automated method to theoretically generate GBs in HCP type materials.

This research examines several critical aspects of the so-called "film induced cleavage" model of stress corrosion cracking using silver-gold alloys as the parent-phase material. The model hypothesizes that the corrosion generates a brittle nanoporous film, which subsequently fractures forming a high-speed crack that is injected into the uncorroded parent-phase alloy. This high speed crack owing to its kinetic energy can penetrate beyond the corroded layer into the parent phase and thus effectively reducing strength of the parent phase. Silver-gold alloys provide an ideal system to study this effect, as hydrogen effect can be ruled out on thermodynamic basis. During corrosion of the silver-gold alloy, the less noble metal i.e. silver is removed from the system leaving behind a nanoporous gold (NPG) layer. In the case of polycrystalline material, this corrosion process proceeds deeper along the grain boundary than the matrix grain. All of the cracks with apparent penetration beyond the corroded (dealloyed) layer are intergranular. Our aim was to study the crack penetration depth along the grain boundary to ascertain whether the penetration occurs past the grain-boundary dealloyed depth. EDS and imaging in high-resolution aberration corrected scanning transmission electron microscope (STEM) and atom probe tomography (APT) have been used to evaluate the grain boundary corrosion depth.

The mechanical properties of monolithic NPG are also studied. The motivation behind this is two-fold. The crack injection depth depends on the speed of the crack formed in the nanoporous layer, which in turn depends on the mechanical properties of the NPG. Also NPG has potential applications in actuation, sensing and catalysis. The measured value of the Young's modulus of NPG with 40 nm ligament size and 28% density was ~ 2.5 GPa and the Poisson's ratio was ~ 0.20. The fracture stress was observed to be ~ 11-13 MPa. There was no significant change observed between these mechanical properties on oxidation of NPG at 1.4 V. The fracture toughness value for the NPG was ~ 10 J/m2. Also dynamic fracture tests showed that the NPG is capable of supporting crack velocities ~ 100 - 180 m/s.

The initial microstructure of oxide fuel pellets can play a key role in their performance. At low burnups, the transport of fission products has a strong dependence on oxygen content, grain size distribution, porosity and grain boundary (GB) characteristics (crystallography, geometry and topology), all of which, in turn depend on processing conditions. These microstructural features can also affect the fuel densification, thermal conductivity and microstructure evolution inside the reactor. Understanding these effects can provide insight into microstructure evolution of fuels in-pile. In this work, mechanical and ion beam serial sectioning techniques were developed to obtain Electron Backscatter Diffraction (EBSD) data, both in 2-D and 3-D, for depleted UO2+X pellets manufactured under different conditions. The EBSD maps were used to relate processing conditions to microstructural features, with emphasis on special GBs according to the Coincident Site Lattice (CSL) model, as well as correlations between pore size and location in the microstructure. Furthermore, larger grains (at least 2.5 times the average grain size) were observed in all the samples and studied. Results indicate that larger grains, in samples manufactured under different conditions, dominate the overall crystallographic texture and have a fairly strong GB texture. Moreover, it seems that the preferential misorientation axis for these GBs, regardless of the O/M, is {001}. These results might be related to GB energy and structure and, suggest that the mechanism that controls grain growth seems to be independent of both processing conditions and stoichiometry. Additionally, a sample was heat treated to relate grain growth and crystallography. The results indicate that at least two mechanisms were involved. Lengthening of GBs was observed for larger grains. Another mechanism of grain growth was observed, in this case, grains rotate to match a neighboring grain forming a larger grain. In the new grain, the misorientation between the two neighboring grains decreases to less than 5 degrees, forming a new larger grain. The results presented in this work indicate that detailed studies of the initial microstructure of the fuel, with emphasis on the crystallography of grains and GBs could help to give insights on the in-pile microstructural evolution of the fuel.

Shock loading is a complex phenomenon that can lead to failure mechanisms such as strain localization, void nucleation and growth, and eventually spall fracture. Studying incipient stages of spall damage is of paramount importance to accurately determine initiation sites in the material microstructure where damage will nucleate and grow and to formulate continuum models that account for the variability of the damage process due to microstructural heterogeneity. The length scale of damage with respect to that of the surrounding microstructure has proven to be a key aspect in determining sites of failure initiation. Correlations have been found between the damage sites and the surrounding microstructure to determine the preferred sites of spall damage, since it tends to localize at and around the regions of intrinsic defects such as grain boundaries and triple points. However, considerable amount of work still has to be done in this regard to determine the physics driving the damage at these intrinsic weak sites in the microstructure. The main focus of this research work is to understand the physical mechanisms behind the damage localization at these preferred sites. A crystal plasticity constitutive model is implemented with different damage criteria to study the effects of stress concentration and strain localization at the grain boundaries. A cohesive zone modeling technique is used to include the intrinsic strength of the grain boundaries in the simulations. The constitutive model is verified using single elements tests, calibrated using single crystal impact experiments and validated using bicrystal and multicrystal impact experiments. The results indicate that strain localization is the predominant driving force for damage initiation and evolution. The microstructural effects on theses damage sites are studied to attribute the extent of damage to microstructural features such as grain orientation, misorientation, Taylor factor and the grain boundary planes. The finite element simulations show good correlation with the experimental results and can be used as the preliminary step in developing accurate probabilistic models for damage nucleation.

The focus of this investigation includes three aspects. First, the development of nonlinear reduced order modeling techniques for the prediction of the response of complex structures exhibiting "large" deformations, i.e. a geometrically nonlinear behavior, and modeled within a commercial finite element code. The present investigation builds on a general methodology, successfully validated in recent years on simpler panel structures, by developing a novel identification strategy of the reduced order model parameters, that enables the consideration of the large number of modes needed for complex structures, and by extending an automatic strategy for the selection of the basis functions used to represent accurately the displacement field. These novel developments are successfully validated on the nonlinear static and dynamic responses of a 9-bay panel structure modeled within Nastran. In addition, a multi-scale approach based on Component Mode Synthesis methods is explored. Second, an assessment of the predictive capabilities of nonlinear reduced order models for the prediction of the large displacement and stress fields of panels that have a geometric discontinuity; a flat panel with a notch was used for this assessment. It is demonstrated that the reduced order models of both virgin and notched panels provide a close match of the displacement field obtained from full finite element analyses of the notched panel for moderately large static and dynamic responses. In regards to stresses, it is found that the notched panel reduced order model leads to a close prediction of the stress distribution obtained on the notched panel as computed by the finite element model. Two enrichment techniques, based on superposition of the notch effects on the virgin panel stress field, are proposed to permit a close prediction of the stress distribution of the notched panel from the reduced order model of the virgin one. A very good prediction of the full finite element results is achieved with both enrichments for static and dynamic responses. Finally, computational challenges associated with the solution of the reduced order model equations are discussed. Two alternatives to reduce the computational time for the solution of these problems are explored.

Dealloying induced stress corrosion cracking is particularly relevant in energy conversion systems (both nuclear and fossil fuel) as many failures in alloys such as austenitic stainless steel and nickel-based systems result directly from dealloying. This study provides evidence of the role of unstable dynamic fracture processes in dealloying induced stress-corrosion cracking of face-centered cubic alloys. Corrosion of such alloys often results in the formation of a brittle nanoporous layer which we hypothesize serves to nucleate a crack that owing to dynamic effects penetrates into the un-dealloyed parent phase alloy. Thus, since there is essentially a purely mechanical component of cracking, stress corrosion crack propagation rates can be significantly larger than that predicted from electrochemical parameters. The main objective of this work is to examine and test this hypothesis under conditions relevant to stress corrosion cracking. Silver-gold alloys serve as a model system for this study since hydrogen effects can be neglected on a thermodynamic basis, which allows us to focus on a single cracking mechanism. In order to study various aspects of this problem, the dynamic fracture properties of monolithic nanoporous gold (NPG) were examined in air and under electrochemical conditions relevant to stress corrosion cracking. The detailed processes associated with the crack injection phenomenon were also examined by forming dealloyed nanoporous layers of prescribed properties on un-dealloyed parent phase structures and measuring crack penetration distances. Dynamic fracture in monolithic NPG and in crack injection experiments was examined using high-speed (106 frames s-1) digital photography. The tunable set of experimental parameters included the NPG length scale (20-40 nm), thickness of the dealloyed layer (10-3000 nm) and the electrochemical potential (0.5-1.5 V). The results of crack injection experiments were characterized using the dual-beam focused ion beam/scanning electron microscopy. Together these tools allow us to very accurately examine the detailed structure and composition of dealloyed grain boundaries and compare crack injection distances to the depth of dealloying. The results of this work should provide a basis for new mathematical modeling of dealloying induced stress corrosion cracking while providing a sound physical basis for the design of new alloys that may not be susceptible to this form of cracking. Additionally, the obtained results should be of broad interest to researchers interested in the fracture properties of nano-structured materials. The findings will open up new avenues of research apart from any implications the study may have for stress corrosion cracking.