Nickel and iron-based alloys are important structure and cladding materials for modern nuclear reactors due to their high mechanical properties and high corrosion resistance. To understand the radiative and corrosive environment influence on the mechanical strength, computer simulation works are conducted. In particular, this dissertation is focused on multiscale modeling of the effects of radiation-induced solute segregation and oxidation on grain boundary (GB) strength in nickel-based and iron-based alloys. Besides the atomistic scale density functional theory (DFT) based calculations of GB strength, continuum-scale cohesive zone model (CZM) is also used to simulate intergranular fracture at the microstructure scale.
First, the effects of solute or impurity segregation at GBs on the GB strength are studied. Thermal annealing or radiation induced segregation of solute and impurity elements to GBs in metallic alloys changes GB chemistry and thus can alter the GB cohesive strength. To understand the underlying mechanisms, first principles based DFT calculations are conducted to study how the segregation of substitutional solute and impurity elements (Al, C, Cr, Cu, P, Si, Ti, Fe, which are present in Ni-based X-750 alloys) influences the cohesive strength of Σ3(111),Σ3(112),Σ5(210) and Σ5(310) GBs in Ni. It is found that C and P show strong embrittlement potencies while Cr and Ti can strengthen GBs in most cases. Other solute elements, including Si, have mixed but insignificant effects on GB strength. In terms of GB character effect, these solute and impurity elements modify the GB strength of the Σ5(210) GB most and that of the Σ3(111) least. Detailed analyses of solute-GB chemical interactions are conducted using electron localization function, charge density map, partial density of states, and Bader charge analysis. The results suggest that the bond type and charge transfer between solutes and Ni atoms at GBs may play important roles on affecting the GB strength. For non-metallic solute elements (C, P, Si), their interstitial forms are also studied but the effects are weaker than their substitutional counterparts.
Nickel-base alloys are also susceptible to stress corrosion cracking (SCC), in which the fracture mainly propagates along oxidized grain boundaries (GBs). To understand how oxidation degrades GB strength, the next step is to use density functional theory (DFT) calculations to study three types of oxidized interfaces: partially oxidized GBs, fully oxidized GBs, and oxide/metal interface, using Ni as a model system. For partially oxidized GBs, both substitutional and interstitial oxygen atoms of different concentrations are inserted at three Ni GBs: Σ3(111) coherent twin, Σ3(112) incoherent twin, and Σ5(210). Simulation results show that the GB strength decreases almost linearly with the increasing oxygen coverage at all GBs. Typically, substitutional oxygen causes a stronger embrittlement effect than interstitial oxygen, except at the Σ3(111). In addition, the oxygen-induced mechanical distortion has a much smaller contribution to the embrittlement than its chemical effect, except for oxygen interstitials at the Σ3(111). For the fully oxidized GBs, three NiO GBs of the same types are studied. Although the strengths of Σ3(112) and Σ5(210) NiO GBs are much weaker than the Ni counterparts, the Σ3(111) NiO GB has a higher strength than that in Ni, indicating that Σ3(111) GB may be difficult to fracture during SCC. Finally, the strength of a Ni/NiO interface is found to be the weakest among all interfaces studied, suggesting the metal/oxide interface could be a favorable crack initiation site when the tensile stress is low.
Furthermore, the effects of co-segregation of oxygen and solute/impurity elements on GB strength are studied by DFT, using the 5(210) GB in an face-centered-cubic (FCC) Fe as a model system. Four elements (Cr, Ni, P, Si) that are commonly present in stainless steels are selected. Regarding the effects of single elements on GB strength, Ni and Cr are found to the increase the GB strength, while both P and Si have embrittlement effects. When each of them is combined with oxygen at the GB, the synergetic effect can be different from the linear sum of individual contributions. The synergetic effect also depends on the spatial arrangement of solute elements and oxygen. If they are aligned on the same plane at the GB, the synergetic effect is similar to the linear sum, although P and Si show stronger embrittlement potencies when they combine with both interstitial and substitutional oxygen. When they are arranged on a trans-plane structure, only nickel combined with oxygen show larger embrittlement potencies than the linear sum in all cases. Crystal Orbital Hamilton Populations analysis of bonding and anti-bonding states is conducted to interpret how the interaction between solutes and oxygen impacts GB strength.
Finally, the continuum-scale CZM method, which is based on the bilinear mixed mode traction separation law, is used to model SCC-induced intergranular fracture in polycrystalline Ni and Fe based alloys in the MOOSE framework. The previous DFT results are used to justify the input parameters for the oxidation-induced GB strength degradation. In this study, it is found that the crack path does not always propagate along the weak GBs. As expected, the fracture prefers to occur at the GB orientations perpendicular to the loading direction. In addition, triple junctions can arrest or deflect fracture propagation, which is consistent with experimental observations. Simulation results also indicate that percolated weak GBs will lead to a much lower fracture stress compared to the discontinuous ones. / Doctor of Philosophy / Iron and Nickel based alloys are important structural materials for nuclear reactors due to their good mechanical properties, corrosion resistance, and radiation resistance. Under radiation and corrosive conditions, those alloys are susceptible to radiation induced segregation (RIS) and stress corrosion cracking (SCC). This dissertation is mainly focused on understanding the influence of the two effects on grain boundary (GB) strength. Systematic atomistic scale density functional theory (DFT) simulations are applied for the nickel and iron systems. Based on the DFT results, cohesive zone model is utilized for the continuum scale fracture simulation in nickel and iron polycrystal.
First, DFT calculations are conducted for studying the RIS effect on the GB strength in nickel. Al, Cr, Cu, C, Si, P, Fe, and Ti are chosen as segregated element. Σ3(111), Σ3(112), Σ5(210), Σ5(310) four types of GBs are built for GB strength calculations. It is found that substitutional C and P always embrittle the GB, while substitutional Ti and Cr can strengthen the GB in most cases. Partial density of states (PDOS) analysis indicates the formation of C-Ni and P-Ni covalent bonds is the possible reason for their embrittlement effects. Bader charge analysis shows negatively charged elements likely reduce the GB strength. Interstitial element segregation is applied for non-metal elements (C, P, and Si). The results indicate interstitial elements have weaker effects than substitution ones.
On the next stage to study the SCC effect, DFT calculations are performed for nickel Σ3(111), Σ3(112), and Σ5(210) GBs with difference oxygen concentration and oxygen incorporation types. Besides partially oxidized GBs, fully oxidized GBs (NiO GBs) and metal-oxide interface are also constructed for comparison. Simulation results show that the GB strength decreases nearly monotonically as oxygen concentration goes up. Typically, substitution oxygen causes a larger embrittlement effect than interstitial oxygen except at Σ3(111). It is found that the large mechanical distortion in this coherent twin GB contributes significantly to the GB strength drop. NiO GBs can be weak (Σ3(112),Σ5(210)) or strong (Σ3(111)). NiO/Ni interface shows lowest strength compared with partially and fully oxidized GBs, indicating the importance of the metal-oxide interface in the SCC process.
Furthermore, the combined effects between segregated elements and oxygen are studied in face center cubic (FCC) iron system. In this part only Σ5(210) GB is selected with substitutional Cr, Ni, P, and Si as segregated elements. The results of single element effects show Cr can strength the GB while P has an opposite effect. Other two elements show little effect. For the co-segregation effects, the trans-plane structures have larger embrittlement potencies than in-plane ones for Ni, suggesting the GB strength can also be affected by the spatial arrangement of segregated elements.
Finally, cohesive zone model is applied for fracture simulations in polycrystalline nickel and iron under tensile loading condition. It is found that intergranular fracture depends on both GB strength and orientation. GBs perpendicular to the loading direction have higher chances to crack. It is also found the percolated weak GBs induce larger strength drop than the discontinuous ones.
Identifer | oai:union.ndltd.org:VTETD/oai:vtechworks.lib.vt.edu:10919/113173 |
Date | 13 January 2023 |
Creators | Xiao, Ziqi |
Contributors | Materials Science and Engineering, Bai, Xianming, Cai, Wenjun, Deshmukh, Sanket A., Corcoran, Sean G. |
Publisher | Virginia Tech |
Source Sets | Virginia Tech Theses and Dissertation |
Language | English |
Detected Language | English |
Type | Dissertation |
Format | ETD, application/pdf, application/pdf |
Rights | In Copyright, http://rightsstatements.org/vocab/InC/1.0/ |
Page generated in 0.0038 seconds