Engineering topics
Casella, Andrew M.
Publications and source records attributed to Casella, Andrew M..
Effect of Hydrogen and Helium on Irradiation Damage Accumulation in Nickel
We report our results on the study of the effects of hydrogen and helium on irradiation damage accumulation in nickel, aiming to understand the TEM observations made by Edwards et al. regarding the Ni-coatings protecting the getter tube in a TPBAR. We hypothesized that differences in the fluxes of 3 H and/or He between the coatings might explain these microstructural variations. However, our simulation results show that during the initial stages of irradiation, the presence of H or He enhances void nucleation by stabilizing small vacancy clusters against thermal dissociation and does not enhance void growth. Nevertheless, we think that helium’s higher diffusivity, stronger binding with vacancy clusters, and ability to displace a nearby lattice atom—either to create a new vacancy or to enlarge an existing vacancy cluster—could result in a larger average vacancy cluster size. Therefore, it is likely that void growth will occur post-nucleation phase as more He accumulates at higher doses than those simulated in this study.
Thermodynamics of Tritium Trapping by Point Defects in Intermetallic Al 12 (TM) 2.35 Aluminide Coating Phase
Density functional theory simulations have been carried out to investigate the potential for tritium trapping by metal vacancies in intermetallic Al 12 (TM) 2.35 phase (TM = Fe, Cr, and Ni) as function of temperature and tritium partial pressure. It was found that tritium could be favorably trapped by Fe and Ni vacancies and not favorably trapped by Al and Cr vacancies. However, due to the presence of partially occupied Al sites in bulk Al 12 (TM) 2.35 , leading to the approximate number of ~255 Al atoms in the unit cell, 86 sites were found energetically favorable to the creation of an Al vacancy. While adding a tritium atom in an Al vacancy is not energetically favorable, the tritiated defect still has a negative Gibbs free energy because the energy gain for creating an Al vacancy overcome the energy cost of adding the tritium species. Based on the calculated Gibbs free energy, the first tritiation of a metal vacancy, at conditions relevant to in-reactor operations, should be more favorable for Al, followed Fe, Ni, and Cr vacancies. By comparing the behavior of tritium in Al 12 (TM) 2.35 with previously studied Fe-Al coating phases (i.e., FeNiAl 5 , Fe 4 Al 13 , and Fe 2 Al 5.6 ), we found that there is a correlation between interstitial tritium solubility and the potential for vacancy trapping. The current trend suggests that if the insertion of an interstitial tritium cost more than 0.3 eV, then trapping by metal vacancies should be preferred. By combining the simulations results obtained to date, we noticed different trapping mechanisms of tritium in the Al coating. Tritium is mostly trapped by Fe and Ni vacancies in the outer Fe-Al coating phase Al 12 (TM) 2.35 while tritium should be preferentially trapped by Al and Fe vacancies for the inner Fe-Al coating phases (FeNiAl 5 , Fe 4 Al 13 , Fe 2 Al 5.6 ). Altogether, these studies show that tritium interacts differently with the various Fe-Al aluminide phases, they also suggest that tritium trapping and retention could be more efficient if metal defects are present and if the solubility of interstitial tritium in the different phases is low.
Molecular dynamics study of interstitial He clusters in nickel
This study presents a molecular dynamics analysis focusing on the behavior of interstitial helium (He) clusters in nickel (Ni), examining their formation, stability, and migration energetics. Consistent with previous research, we found that the binding energies of interstitial helium within a helium cluster are positive and increase with the cluster size, indicating a preference for helium atoms to cluster together. However, our findings also reveal that while the formation energy increases monotonically with cluster size, the increase in binding energy is non-monotonic. Importantly, small He clusters were observed to be thermally unstable at reactor operational temperatures (approximately 600 K), with the He 2 cluster exhibiting instability even at room temperature. With a binding energy of 0.49 eV for a He 4 cluster, we hypothesize that for helium bubbles to form via homogeneous nucleation (i.e., through trap mutation) at reactor operating temperatures, the helium concentration must be high enough to facilitate the formation of helium clusters of at least size 4 or larger. As expected, interstitial helium and small helium clusters are highly mobile. This mobility was observed not only at room temperature but also at temperatures as low as approximately 200 K. Furthermore, the mean squared displacement method has been utilized to determine the migration barriers and the corresponding prefactors for clusters ranging from He 1 to He 6
Effects of irradiation temperature on the microstructure and deuterium retention in γ-LiAlO 2 pellets
This report presents the experimental findings obtained from November 2023 to September 2024. The study aims to investigate the effects of temperature on the microstructure, deuterium (D) retention, and lithium (Li) loss in γ-LiAlO 2 pellets subjected to sequential He + and D + ion irradiation to a high dose. Sequential ion irradiation has been a key method in our previous studies to simulate the behavior of γ-LiAlO 2 under neutron irradiation. In addition to dose and dose rate, irradiation temperature is a critical factor influencing microstructural and compositional changes. LiAl 5 O 8 precipitates have been observed in γ-LiAlO 2 pellets irradiated with reactor neutrons at ~300°C. These precipitates also form during ion irradiation at an elevated temperature of 500°C, but not at 300°C. In our ion irradiation experiments, the dose rate is typically three orders of magnitude higher than that of neutron irradiation, leading to a more rapid damage production. To better emulate the microstructural features in neutron-irradiated pellets using ion irradiation, a higher irradiation temperature is needed to accelerate the diffusion of point defects and enhance defect recovery rates, thereby compensating for the effects of the higher dose rate. The microstructural changes observed are the result of competing processes occurring during ion irradiation.
First-principles Study of the Tritium Formation in γ-LiAlO2 Pellets and Diffusion into Ni-Zircaloy-4 Getter
Presentation for American Physical Society (APS) March meeting.
Detection of tritium using ultrafast laser-induced breakdown spectroscopy
Detection of protium, deuterium and tritium using ultrafast LIBS with rapid, and standoff capability in addition to no sample preparation requirement which are crucial to nuclear nonproliferation, safeguards, and security
Cluster dynamics simulations of tritium and helium diffusion in lithium ceramics
Tritium (T) and He diffusion in LiAlO 2 and LiAl 5 O 8 phases influences the performance of tritium producing burnable absorber rods (TPBARs) by affecting the gas release, swelling and thermal conductivity of Li-bearing ceramic pellets. Frenkel pair defects and clusters created by irradiation can attract T and He interstitials and form clusters of the type He i x Li, He i x Al, He i x O, T i x Li, T i x Al, and T i x O, 1 ≤ x ≤ 4 in a Li, Al or O vacancy site (notation denotes x He or T atoms in a 1 Li, 1 Al or 1 O vacant site). The concentration and mobility of each of these clusters collectively contribute to the diffusion of the He and T gases in LiAlO 2 and LiAl 5 O 8 . Here, in this work, free energy cluster dynamics simulations implemented in the Centipede code, are used to obtain the concentration and diffusivities of these clusters which are then used to calculate the total diffusivity of T and He gases in LiAlO 2 and LiAl 5 O 8 . The results show that diffusivity of T is at least one order of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 whereas He diffusion is 2–13 orders of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 . There is a higher concentration of highly diffusive species (T interstitials and T i 03 Li for the case of tritium and He i 01 Li, He i 02 Li, and He i 03 Li for the case of He) in LiAlO 2 than in LiAl 5 O 8 which increase the total diffusion of T and He in LiAlO 2 .
Molecular dynamics simulations of displacement cascades in LiAlO2 and LiAl5O8 ceramics
Abstract Molecular dynamics was employed to investigate the radiation damage due to collision cascades in LiAlO 2 and LiAl 5 O 8 , the latter being a secondary phase formed in the former during irradiation. Atomic displacement cascades were simulated by initiating primary knock-on atoms (PKA) with energy values = 5, 10 and 15 keV and the damage was quantified by the number of Frenkel pairs formed for each species: Li, Al and O. The primary challenges of modeling an ionic system with and without a core–shell model for oxygen atoms were addressed and new findings on the radiation resistance of these ceramics are presented. The working of a variable timestep function and the kinetics in the background of the simulations have been elaborated to highlight the novelty of the simulation approach. More importantly, the key results indicated that LiAlO 2 experiences much more radiation damage than LiAl 5 O 8 , where the number of Li Frenkel pairs in LiAlO 2 was 3–5 times higher than in LiAl 5 O 8 while the number of Frenkel pairs for Al and O in LiAlO 2 are ~ 2 times higher than in LiAl 5 O 8 . The primary reason is high displacement threshold energies (E d ) in LiAl 5 O 8 for Li cations. The greater E d for Li imparts higher resistance to damage during the collision cascade and thus inhibits amorphization in LiAl 5 O 8 . The presented results suggest that LiAl 5 O 8 is likely to maintain structural integrity better than LiAlO 2 in the irradiation conditions studied in this work.
Microstructural and compositional evolutions in $γ$-LiAlO 2 pellets during ion irradiation at an elevated temperature
This study reports on the microstructural, phasic and compositional evolutions of γ-LiAlO 2 during ion irradiation. Polycrystalline γ-LiAlO 2 pellets were irradiated sequentially with He + and D 2 + ions to the same combined fluences of up to 3 × 10 17 (He + +D + )/cm 2 at 773 K. The irradiated pellets were characterized using scanning transmission electron microscopy and atom probe tomography. Surface amorphization likely due to radiolysis and planar defects as a possible precursor for formation of precipitates are created at 5 × 10 16 (He + +D + )/cm 2 , followed by the formation of nano-sized precipitates and fractures at higher doses. Spinel-like precipitates of non-stoichiometric LiAl 5 O 8 and gas-filled cavities are observed to grow with increasing dose. Faceted precipitates and rounded cavities appear at 2 × 10 17 (He + +D + )/cm 2 . With further increasing ion fluence to 3 × 10 17 (He + +D + )/cm 2 , amorphization of the precipitates takes place and micron-sized fractures appear. Surface exfoliation could occur at an extremely high ion fluence. There are compositional changes in the γ-LiAlO 2 pellets during the microstructural evolution. In the precipitate and amorphized regions, Li concentrations decrease to ~7 and 3.7 at.% from 25 at.% in γ-LiAlO 2 , respectively. This study reveals a full-cycle microstructural evolution with corresponding compositional changes in γ-LiAlO 2 pellets during ion irradiation at 773 K. In conclusion, the data could help model, assess, and predict the material performance during neutron irradiation in reactors.
Ab Initio Simulations of Tritium Diffusion in Al 2 O and Intermetallic Al 12 (TM) 2.34 Aluminide Coating Phases
Density functional theory simulations have been carried out to investigate the diffusion of interstitial tritium in Al 12 (TM) 2.34 and Al 2 O bulk phases. In Al 12 (TM) 2.34 the transition metal (TM) sites are occupied on average by 58.93 at.% Fe, 18.52 at.% Cr, and 22.54 at.% Ni. While the insertion of interstitial tritium in Al 2 O lead to a strong disordering of the structure, we only investigated interstitial tritium in Al 12 (TM) 2.34 and its iron end-member (i.e., Al 12 Fe 2.34 ). Nine diffusion pathways have been investigated along the a-, b-, and c -axis of Al 12 (TM) 2.34 . The direction of fastest diffusion for interstitial tritium is found to be along the b-axis, with a calculated diffusion coefficient of ≈10 -10 m 2 .s -1 at 600 K, which is at least one order of magnitude faster than those previously calculated for interstitial tritium in other Al-rich phases such as Fe2Alx, Fe4Al13, and FeNiAl5 for which D T =10 -11 m 2 .s -1 , ≈10 -12 m 2 .s -1 , and D T ≈10 -13 m 2 .s -1 respectively. The comparison of the energy landscape between Al 12 (TM) 2.34 and Al 12 Fe 2.34 for the fastest diffusion pathways in each axis direction, found that some pathways are not sensitive to transition metal mixing, while other are more affected, leading to higher energy barrier in Al 12 Fe 2.34 in part due to a more energetically favorable formation of Fe—T bond compared to Ni—T. However, we found that both materials have similar diffusion coefficients as the fastest diffusion pathways occurs along pathways that are not very sensitive to transition metal mixing.
Insights into radiation resistance of titanium alloys from displacement cascade simulations
Radiation damage in beam window materials limits the use of high power proton beams in high energy physics research. The alloy Ti-6Al-4V is presently used as a beam window material but a prospective alternative, Ti-15V-3Cr-3Sn-3Al has been proposed. Since both these alloys contain dual phases at room temperature, we compare the radiation damage in the α and ß-phases of these two materials via primary knock-on atom (PKA) cascade simulations in the 10-40keV energy range. At PKA energies 30 and 40keV, the number of Frenkel pairs in the ballistic stage is higher in the ß-phase of Ti-6Al-4V than that in the ß-phase of Ti-15V-3Cr-3Sn-3Al almost by a factor of 2. The α-phase, of both these alloys, by far outperforms the ß-phases of the two alloys, both in terms of damage during the ballistic stage and in terms of the surviving defects. The average displacement threshold energy (E d ) in the α-phase of both alloys was found to be 66eV while that in the ß-phase of Ti-15-3 was 55 eV and in the ß-phase of Ti-6-4 was 46. Further, while the number of surviving defects is almost equal in both alloys, the vacancy and interstitial clustering mechanisms differ notably, which can impact the degree of radiation hardening and loss of ductility. Our simulations show larger vacancy and interstitial clusters form in Ti-6Al-4V as compared to that in Ti-15-3-3-3 alloy. These results indicate that Ti-15-3-3-3 alloys may be a promising candidate for next generation beam window material with a higher radiation tolerance than the existing Ti-6-4 alloy.
Study of Tritium Diffusivity in Pure and Sn-Defective Zr: A First-Principles Density Functional Theory Approach
Zirconium alloys (e.g., zircaloy-4) are used as tritium ( 3 H) getter materials in tritium-producing burnable absorber rods (TPBARs) owing to their ability to capture 3 H and chemically convert them into metal hydrides. Understanding of 3 H diffusion mechanisms in zircaloy is crucial for the optimal design of material performance in nuclear technology. Here, for this work, we perform first-principles density functional theory calculations to study the 3 H diffusion mechanism in pure and impure Zr with a low concentration of tin (Sn) atoms to determine the impact of the presence of Sn on the movement of 3 H atoms through the material. First, we calculated the diffusion barriers for 3 H in pure Zr by taking different migration pathways. We then introduced a low concentration of Sn impurity and systematically explored the impurity effect on the diffusion barriers for 3 H. Using our calculated diffusion energy barriers, we further obtained the diffusion coefficients and analyzed the results by comparing them with the experimental and previously calculated values. A diffusion coefficient of the order of 10 –8 m 2 /s is predicted. We also found that the presence of a Sn impurity could reduce the diffusivity up to 4 orders of magnitude. In conclusion, our results could serve as guidelines for further experimental investigations.
First-principles study of the tritium diffusion and formation in γ-LiAlO2 pellets
TechConnect World Innovation Conference and Expo presentation
Ab Initio Simulations of Tritium Diffusion in Al-Rich Iron Aluminide Coating Phases
Surface coatings of steels used in extreme conditions and corrosive environments generally aim to provide protection and increased durability. In the case of tritium-producing burnable absorber rods (TPBARs) used in nuclear reactors, a 316 stainless steel has been coated with Al. Scanning transmission electron microscopy (STEM) characterization of the coating found three Al-rich (>60 atom % Al) iron aluminide alloys identified as hexagonal FeNiAl 5 , monoclinic Fe 4 Al 13 , and orthorhombic Fe 2 Al 5 . Density functional theory simulations using nudged elastic band have been performed to investigate the diffusion of interstitial tritium in each Al-rich iron aluminide phase. While FeNiAl 5 and Fe 4 Al 13 can be viewed as the stacking of two layers, the structural peculiarity of Fe 2 Al 5 is that channels of variable Al vacancy content are present along the c-axis. Therefore, three stoichiometries for Fe 2 Al x phase, namely, Fe 2 Al 4 , Fe 2 Al 5 , and Fe 2 Al 6 , have been considered to evaluate the impact of Al vacancy concentration on tritium diffusion behavior. Altogether, we found that at 600 K, tritium diffusion decreases from a faster rate in the channels of Fe 2 Al x phases (D T ≤ 10 –11 m 2 ·s –1 ) to Fe 4 Al 13 (D T ≈ 10–12 m2·s–1), and finally in FeNiAl5 (D T ≈ 10 –13 m 2 ·s –1 ). Here, we also find that interstitial tritium generally diffuses faster in Fe–Al coating phases than in the tritium breeding material γ-LiAlO 2 (D T ≈ 10 –14 m 2 ·s –1 ) but slightly slower than in 316 stainless steel (D T ≈ 10 –10 m 2 ·s –1 ).
Molecular dynamics study of primary damage in the near-surface region in nickel
In this work, we carried out a large-scale molecular dynamics (MD) simulation to elucidate the effect of the free surface on the defect production of displacement cascades in pure nickel. These MD simulations were performed both in the bulk and near-surface regions with primary knock-on atom (PKA) energies of E PKA = 1, 5, or 10 keV and at temperatures T = 300, 425, or 525 K. Additionally, for every (T, E PKA ), near-surface cascade simulations were performed as a function of depth. In the near-surface simulations, assuming isotropic neutron fluence, the PKAs were initiated in random directions, including toward and parallel to the free surface. For both the near-surface and bulk cascades, the effect of E PKA and T on the defect production is similar. In both cases, the defect production increases with E PKA , but the temperature has minimal effect. However, the production and clustering of vacancies are higher for near-surface cascades, and they decrease with increasing depth. In contrast, the production and clustering of self-interstitial atoms are lower and increase with depth. Eventually, the production and clustering of both defect types approach bulk cascade-like behavior, and the depth at which this occurs increases with E PKA , but is independent of temperature.
Influence of ambient gas on self-reversal in Li transitions relevant to isotopic analysis
Laser induced breakdown spectroscopy is a promising, rapid analysis method for the detection and quantification of Li and its isotopes needed in geochemical, nuclear, and energy storage applications. However, spectral broadening in laser produced plasmas, presence of fine and hyperfine structures, and self-reversal effects make Li isotopic analysis via laser induced breakdown spectroscopy challenging. The present study explores the influence of Ar, N 2 , and He ambient gases over the pressure range of 0.05 - 100 Torr on line broadening and self-reversal of the Li I transition with the greatest isotopic shift in the VIS spectral region (i.e., ≈670.8 nm, ≈15.8 pm isotopic shift). We perform spatially and temporally resolved optical emission spectroscopy of plasmas produced via laser ablation of LiAlO 2 substrates. Our results show that the self-reversal and linewidth is reduced at lower pressures for all gases, and using optimized plasma conditions with chemometric methods, the 6 Li/ 7 Li isotopic ratios can be predicted.
Molecular dynamics simulations of radiation response of $\mathrm{LiAlO_2}$ and $\mathrm{LiAl_5O_8}$
Here, we report the findings of a computational investigation of defect production and migration in two lithium aluminate ceramics: LiAlO 2 and LiAl 5 O 8 . Although the first ceramic, LiAlO 2 has been widely investigated computationally, especially via molecular dynamics (MD), there have been no MD studies done on LiAl 5 O 8 up to date, due to the lack of interatomic potentials for this system. This work implemented an existing set of Li-Al-O coulombic and Buckingham potential parameters for the LiAl 5 O 8 system from the literature, followed by validating it based on experimental values of density, crystallinity, bond lengths and melting point. The Li + diffusion in the two ceramics was examined and was found to be one order of magnitude slower in LiAl 5 O 8 (D 0 = 3.17 x 10 -11 m 2 /s) as compared to that in LiAlO 2 (D 0 = 4.02 x 10 -10 m 2 /s) at 600 K. The lithium vacancy migration barrier of 4.15 eV in LiAl 5 O 8 was more than three times that in LiAlO 2 (1.31 eV) possibly due to stronger ordering in the Li-O 6 octahedra of LiAl 5 O 8 than in Li-O 4 tetrahedra of LiAlO 2 . The Li displacement threshold energy (E d ) in LiAl 5 O 8 was found to have a median value of 68 eV which is much higher than that in LiAlO 2 (40 eV). The greater E d for Li signifies a lower likelihood for defect formation and clustering and thus lower tendency for amorphization in LiAl 5 O 8 . The presented results show the difference in defect dynamics in the two ceramics that may help us understand the susceptibility of the two ceramics to amorphization caused by irradiation.