Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Graphite modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

SCALE modeling of foil irradiations at WSU’s TRIGA with sensitivity/uncertainty analysis

We report the Washington State University TRIGA reactor was modeled in the SCALE system’s KENO neutron transport code. The model includes 119 fuel rods, a central boral control rod, 4 control blades, 5 irradiation ports, and 20 graphite reflector blocks within a rectangular aluminum grid box. The model accurately simulated the irradiation of depleted uranium and highly enriched uranium foils and predicted reaction rate ratios for 238 U(n,2n) 237 U / 238 U(n,γ) 239 U and 99 Mo / fission, which compared favorably to empirical values previously measured at Los Alamos. Finally, a sensitivity/uncertainty analysis was performed using the continuous-energy TSUNAMI-3D module of SCALE.

237Np↗

FY22 Status Report on the ART-GCR CMVB and CNWG International Collaborations

This work presents the numerical model for the high temperature test reactor (HTTR) loss of forced cooling (LOFC) experiment with the INL codes Griffin, BISON, and RELAP-7, based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Promising results were obtained, with the overall behavior of the reactor successfully captured. Changes in the heat transfer coupling, as compared to the fiscal year (FY)-21 model, enabled drastic improvement of the steady-state solution, both in terms of computational time and global energy discrepancies. The former was reduced by a factor of roughly 60, whereas the latter decreased from 11% to less than 2%. Furthermore, the discrepancy in the steady-state multiplication factor was improved from +2,300 and +2,900 pcm (for the 30 and 9 MW cases) to -700 and +1,200 pcm, respectively, and now falls well within the large measurement uncertainties stemming from graphite impurities. Validation of the Monte Carlo model used to generate cross sections was also performed against available measurements. Though significant, the discrepancies remain acceptable overall in light of the large uncertainty. Numerous improvements are still needed to better compare with the experiment involving the 9 MW case and to instill greater confidence in the model’s ability to accurately predict the 30 MW behavior. Specifically, the power levels predicted by the 9 MW transient simulation following re-criticality remain low, pointing to an underestimation of the passive cooling of the core. A key aspect of future work will be to better understand the flow pattern during the LOFC event, particularly to determine if natural or forced convection is occurring inside the reactor pressure vessel (RPV). More generally, additional validation data would be immensely useful for further enhancing the numerical model and better matching the experiments. In addition, a more sophisticated thermal-hydraulics model that simulates all the channels as a single system model should be considered to take into account the rest of the primary loop. Finally, even if the results are in better agreement with the experiments, sensitivity analysis and uncertainty quantification will be necessary to evaluate the model uncertainty.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Laser ablation for structuring Li-ion electrodes for fast charging and its impact on material properties, rate capability, Li plating, and wetting

Laser ablation is a scalable technique for decreasing the effective tortuosity of electrodes by selectively removing material with high precision. Applied to ≈ 110 um thick electrode coatings, this work focuses on understanding the impact of laser ablation on electrode material properties at the beginning of life and synergistic impacts of ablated channels on cell performance throughout their cycle life. Post laser ablation, local changes in chemistry, crystallography, and morphology of the laser-impacted electrode regions are investigated. It is shown that femtosecond pulsed laser ablation can achieve high-rate material removal with minor material damage locally at the interface of the impacted zones. The capacity achieved during a 6C (10 min) constant-current constant-voltage charge to 4.2 V improved from 1 mAh cm -2 for the non-ablated electrodes to almost 2 mAh cm -2 for the ablated electrodes. This benefit is attributed to a synergistic effect of enhanced wetting and decreased electrode tortuosity. The benefit was maintained for over 120 cycles, and upon disassembly decreased Li-plating on the graphite anode was observed. Finally, multi-physics modeling in conjunction with wetting analyses showed that laser ablating either one of the electrodes led to substantial improvements in wetting and rate capability, indicating that substantial performance benefits can be achieved by ablating only the graphite anode as apposed to both electrodes.

25 ENERGY STORAGE↗

Using porous random fields to predict the elastic modulus of unoxidized and oxidized superfine graphite

Nuclear graphite is a candidate material for Generation IV nuclear power plants. Porous materials such as graphite can contain complex networks of pores that influence the material's mechanical and irradiation response. A methodology known as the random finite element method (RFEM) was adapted to create synthetic microstructures and predict the influence of porosity on the elastic properties of graphite during oxidation. RFEM combines random field theory and the finite element method in a Monte Carlo framework to estimate the mechanical response of a given grade of graphite. In this research, the random fields were verified through experimental characterization to predict the elastic response of three nuclear graphite grades, ETU-10, IG-110, and 2114. Finite element models (FEM) were generated using segmentations of x-ray computed tomography (XCT) data known as image-based models (IBMs) to validate and compare with the RFEM results and better understand the effects of uniform oxidation in these graphite grades. The RFEM predictions appear to correlate well with the experimental values of the measured Young’s modulus of the three graphite grades and display the same trends as IBMs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Examining Graphite Degradation in Molten Salt Environments: A Chemical, Physical, and Material Analysis

Molten-salt reactors (MSRs) are Generation IV nuclear reactors that use liquid salt as a coolant and/or fuel. In several MSR designs, graphite serves as a moderator and/or reflector. However, due to limited experimental data and operational experience, our understanding of graphite behavior in molten salt environments remains incomplete. This report aims to identify the degradation mechanisms of nuclear graphite in MSRs, detail the mechanisms of each factor, and provide an initial assessment of their impact on the structural integrity of graphite components. This assessment is based on an extensive literature review and insights from subject matter experts. Furthermore, given the limited data, a modeling strategy using existing Grizzly software is proposed for a more thorough analysis where appropriate. Additionally, it presents mitigation strategies where applicable. The report covers physical degradation mechanisms such as infiltration, erosion, and abrasion, as well as chemical degradation mechanisms including fluorination, intercalation, corrosion, and oxidation. Molten salt can infiltrate the porous structure of graphite, leading to several detrimental effects. Entrapment of fissile products within the graphite pores can cause radiation damage and could pose challenges in the handling and disposal of contaminated components. The differential thermal expansion between the infiltrated salt and graphite, along with internal stress from pressurized molten salt and volumetric heating, can compromise the structural integrity of graphite. To mitigate these effects, employing ultra-fine graphite grades and applying sealants and coatings are effective strategies. A computational model based on coupled solid mechanics and heat transfer phenomena could be used to predict the internal stresses using Grizzly software. In pebble-bed MSRs, graphite fuel pebbles can cause abrasion against reactor components due to friction and wear. The severity of wear is influenced by various factors such as temperature, environment, and the presence of lubricants. Tribological studies reveal that higher temperatures and molten salt environments, such as FLiBe, significantly reduce wear rates compared to dry conditions. Additionally, the chemical composition of the salt can further optimize graphite's tribological performance. Long-term wear effects can be modeled by incorporating surface defects into the geometry and predict stresses under thermal and radiation effects using Grizzly software. Chemical degradation of graphite in a molten salt environment can occur through fluorination and intercalation. Fluorination can occur via replacement of hydrogen or oxygen atoms, or at the active sites, but does not cause structural degradation. Intercalation, on the other hand, can lead to exfoliation, where layers of graphite separate and peel away, damaging the graphite. Protective coatings can enhance graphite's resistance to intercalation. Graphite generally exhibits good chemical stability in molten salt environments, though it can corrode under specific conditions, particularly in the presence of impurities or oxidants. Studies have shown that protective coatings, such as plasma-sprayed partially stabilized zirconia (PSZ), can effectively prevent such degradation. Corrosion behavior varies significantly with different graphite grades and coating applications, underscoring the need for detailed studies on uncoated and coated graphite to understand and mitigate corrosion mechanisms in MSRs. Research indicates that the presence of oxidants and impurities can accelerate graphite degradation in molten salts, making it essential to explore acceptable impurity limits. Oxidation is another critical degradation mechanism, leading to weight loss and structural damage due to the formation of CO and CO 2 from the reaction of carbon atoms with oxygen. This process creates new porosity and compromises graphite's integrity. While extensive research on graphite oxidation has been conducted for gas-cooled reactors, studies specific to MSRs are limited. Findings from the coal industry suggest that molten alkali metal salts can significantly accelerate graphite oxidation, a hypothesis worth exploring for fluoride salts in MSRs. Understanding oxidation behavior in MSRs is vital for developing protective measures. The analysis of post-irradiated graphite from the MSRE experiment demonstrated exceptional chemical compatibility with molten fluoride salt, suggesting that the extent of chemical attack on graphite largely depends on the salt's infiltration capability. Therefore, the use of ultra-fine grade graphite could help mitigate chemical degradation effects. Existing oxidation modeling capabilities in Grizzly, which use reaction-diffusion equations to model graphite-air interactions, could be adapted to simulate the chemical degradation effects of graphite in molten salt environments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Characterization of thermally heat-treated polyacrylonitrile carbon fibers

This study investigates the graphitization process of polyacrylonitrile (PAN) carbon fibers by subjecting commercial fibers to thermal heat treatment at temperatures ranging from 1400 to 2100 °C in 100 °C increments, using either argon or nitrogen gas atmospheres. Changes in crystallinity, surface morphology, and lattice parameters were analyzed for two commercial carbon fibers using X-ray diffraction, scanning electron microscopy, and Raman spectroscopy. Results indicated minimal changes in surface morphology with increasing heat-treatment temperature; however, crystallinity significantly increased. Crystallinity changes were more strongly dependent on temperature rather than gas atmosphere or fiber type. At intermediate heat-treatment temperatures (1600–1800 °C), fibers treated in argon showed a slight preference for graphitization. The highest level of graphitization was measured at 2100 °C. Crystallite size increased as the intensity ratio of the D1 to G Raman peaks increased, reaching a peak around ~1800 °C, after which the ratio started to decrease. This behavior aligns with Ferrari's three-stage model of carbon crystallization and is consistent with both the Marie-Mering degree of graphitization and Brubaker's Integrated Absolute Differential models, all of which describe the transformation from an amorphous to a more graphitic structure. At the higher heat-treatment temperatures, the changes between atmospheres and fiber types were measured to converge to similar levels of graphitization. In conclusion, this study evaluates the progressive change in commercial grade carbon fibers when heat-treated.

Characterization↗

Modelling dust transport in DIII-D with DTOKS-Upgrade

Comprehensive upgrades to the dust transport code Dust in TOKamaks (DTOKS) that extend the plasma-dust interaction model are presented and compared with recent measurements of dust transport in DIII-D. Simulations incorporating variation in physical properties of graphite dust with temperature and size in a stationary plasma background suggest a substantial decrease in lifetimes due principally to thermal expansion. The trajectories of 53 dust grains identified from analysis of visible camera data taken across two similar shots were used to measure the dust particle velocity distributions. Dust tracks terminated mostly at the outer divertor strike point having a mean observation time of 2.1 ± 0.4 ms. Stochastic modelling of 200 graphite dust particles in the DIII-D tokamak performed with DTOKS-U using plasma simulations generated by OEDGE found similar behaviour, with particles ablating rapidly after acquiring a positive charge in the region close to the outer strike point, creating an acute source of neutral carbon atoms. The simulated mean lifetime, 11 ± 2 ms, showed approximate agreement with experimental observation when corrected by accounting for dust visibility and ignoring the longest trajectories 1.5 ± 0.2 ms. Synthetic diagnostic data generated from coupling the results of DTOKS-U with the visualisation software Calcam offers a powerful new tool for validation of simulations and predictive calculations of dust dynamics.

Physics↗

Status memorandum on ORNL support of vendor irradiation capsule design and ASME irradiation code development and new ASTM test standard development activities

This report is in submission of completion of the Level 4 milestone number M4TG-24OR0501111 within the larger Advanced Reactor Technologies Gas Cooled Reactor program at Oak Ridge National Laboratory (ORNL). The focus of this year’s efforts is to continue support of industry needs related to the development of future industry-funded irradiation programs. At ORNL these efforts included the development of a generalized three-dimensional CAD model of an irradiation creep capsules for graphite, cost estimates of the expansion of the Materials Irradiation Facility (MIF), involvement with ASTM D02.F0 “Manufactured Carbon and Graphite Products” and the ASME Boiler and Pressure Vessel Code, and publication of papers supporting these efforts. This report will document the status of these activities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Secondary electron emission for reticulated carbon foam surfaces using direct measurements and spectroscopic analysis

This study investigates secondary electron emission (SEE) characteristics of reticulated foams using direct measurements and analytical modeling. Total SEE was quantified, revealing suppression of up to 44% in carbon foam structures compared to planar graphite surfaces. An optimal geometric configuration was identified and supported by analytical models. SEE angular dependence experiments showed diverse behaviors: fiber-like behavior and directional dependence for pore and ligaments on the mm scale, with fuzz-like characteristics when the foam features are between 10–100 µm. Electron energy analyzer measurements showed that carbon foams preferentially suppress inelastic backscattered electrons (BSEs) more so than true secondary electrons (SEs). The analysis indicated a larger fraction of low-energy SE generation in foams compared to flat surfaces due to increased emission from curved fiber ligaments and tertiary SEs from high-energy BSEs. These findings have implications for design and optimization of materials with tailored electron emission properties for applications like plasma-facing components, spacecraft materials, and accelerator surfaces.

Auger↗

Development and Validation of MALAMUTE model for Electric Field Assisted Sintering of Structural Materials

Fusion power plant designs feature extreme material performance requirements for structural material candidates. In addition to conventional alloys, more advanced composites and oxide dispersion strengthened (ODS) alloys are being explored, however, achieving the desired microstructures to maximize performance using traditional manufacturing methods can be challenging. The advanced manufacturing (AM) electric field-assisted sintering (EFAS) technique offers improved control over the final microstructure through higher heating and cooling rates and moderate pressures. Modeling and simulation tools show promise in elucidating the process-structure-property-performance (PSPP) correlation for AM-produced parts, including the EFAS process. An inherently multiscale process, the EFAS technique aligns well with the multiscale modeling capability of the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE)[cite]. We present here an electro-thermo-mechanical approach to modeling the EFAS process using the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) code. Prediction of the field and gradient distributions across the EFAS tooling is required to accurately describe the conditions for the lower-scale microstructural evolution models. In this work we present the MALAMUTE model developed to predict the electrical potential, temperature, and mechanical stress distribution across the EFAS graphite tooling and part at the larger engineering-scale. Validation of the MALAMUTE engineering-scale model is completed using data from experimental densification and pre-densified runs of iron powder via EFAS at 1000oC. These runs were conducted using a Thermal Technology DCS-5 EFAS system. Data collected during the experiment runs include the direct current (DC) supplied to the graphite tooling, the temperature of the graphite tooling as measured with a pyrometer, and the force applied to the top of the graphite tooling stack, and the data were recorded every 10 seconds. Our validation approach used the current and force data from the EFAS run as boundary condition inputs to the MAMALUTE simulation; the temperature data were used to evaluate the MALAMUTE EFAS model prediction. Results of the MALAMUTE simulations are employed to connect the external pyrometer temperature measurement to the temperature profile across the part undergoing consolidation. We investigate the impact of material property variation and mesh deformation on the temperature profile as predicted by MALAMUTE. We conclude by highlighting projects where the MALAMUTE EFAS modeling and simulation capabilities will be used to assist experimental design.

36 - MATERIALS SCIENCE↗

PSA 2025 Presentation: "Modeling and Sensitivity Analysis of a Generation IV Pebble Bed Reactor Using MELCOR 2.2"

Accompanying the advancement of reactor technologies is the need for computational modeling and simulation to predict their behavior under normal operating conditions and accident scenarios. New Generation IV reactor designs which employ non-conventional fuel have a particular need for modeling the behavior and release of radionuclides and other material from the fuel. In this work, MELCOR version 2.2, a system-level safety and accident scenario code developed by Sandia National Laboratories, was used to model a 200-MWth pebble bed modular reactor and calculate the inventories of circulating and deposited graphite, metal dust, and elemental components released from the fuel elements. A base case modeling the reactor under standard operating conditions was calculated using MELCOR and the inventories were extrapolated to 30 years of operation time using a logarithmic regression fit. A sensitivity analysis was also performed in which several key parameters for the base case model were modified to explore the effect of these changes on the inventories calculated by MELCOR. A set of transient scenario simulations for a depressurized loss of forced cooling (DLOFC) accident were also performed. The results of the sensitivity analysis and transient simulations are reported and discussed in relation to the modeling techniques used for this study.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Chemical Interplay of Silicon and Graphite in a Composite Electrode in SEI Formation

In this study, we investigated the effect of the Si/graphite weight ratio in half-cells on the solid electrolyte interphase (SEI) layer's chemistry. The nominal concentrations of active materials were (wt % Si/wt % Gr) 15/73, 30/58, 60/28, and 80/0. The electrolyte in the cells consisted of either 1.2 M LiPF 6 in ethylene carbonate/ethyl methyl carbonate (3:7 by wt) or 1.2 M LiPF 6 in ethylene carbonate:ethyl methyl carbonate (3:7 by wt) + 10 wt % fluoroethylene carbonate. These coin cells were cycled five times at the C/10 rate. As expected, the addition of silicon to the electrode significantly increased the measured capacity. Examination of the aged composite material showed that the electrolyte influenced the concentration of chemical environments on the surface. Depth profiling revealed that these concentrations of surface environments changed with sputtering time. A statistics-of-mixtures model was used to deconvolute how silicon and graphite interacted during the formation of these species and how the interaction changed with depth.

25 ENERGY STORAGE↗

A “Master Curve” Describing Reaction Inhomogeneity and Plating Onset during Fast-Charging of Graphite Electrodes

Lithium plating in porous graphite electrodes is a major limitation for fast charging. Theoretical evidence suggests plating during fast charging is largely due to inhomogeneous intercalation through the electrode thickness (caused by ionic/mass transfer limitations and inherent thermodynamic properties of the graphite material). Numerical and analytical predictions of plating onset during fast charging have been proposed but not compared directly with experiments. This work validates these model predictions against plating onset measurements via the “dOCV” method in graphite half-cells, for various electrode thicknesses and C-rates. Remarkably, it is shown that experimental and theoretical trends in plating onset (vs C-rate and electrode thickness) collapse to a single “master curve” using a nondimensional “reaction inhomogeneity” parameter. This observation supports the hypothesized reaction inhomogeneity mechanism for accelerating plating onset and provides practical guidance for electrode design. Furthermore, this work develops theory to quantify reaction inhomogeneity in situ directly from voltage V vs capacity Q data, using dQ / dV analysis. Here again, experiments and numerical predictions show good agreement, where peaks that correspond to LiC 6 correlate with plating onset. This work provides experimental validation of theoretical tools that can predict plating onset, aid electrode design, and give insight for plating onset mechanisms during fast charging.

25 ENERGY STORAGE↗

Influence of Polymorphs and Local Defect Structures on NMR Parameters of Graphite Fluorides

In this study, the role of local molecular structure on calculated 13 C and 19 F NMR chemical shifts for graphite fluoride materials was explored by using gauge-including projector augmented wave (GIPAW) computational methods for different periodic crystal polymorphs and density functional theory (DFT) gauge-including atomic orbital (GIAO) computational methods for individual graphite fluoride platelets, i.e., fluorinated graphene (FG). The impact of stacking sequences, d -spacing, and ring conformations on fully fluorinated graphite fluoride structures was investigated. A range of different defects including Stone–Wales, F and C vacancies, void formation, and F inversion were also evaluated using FG structures. These calculations show that distinct chemical shift signatures exist for many of these polymorphs and defects, therefore providing a basis for spectral assignment and development of models describing the mean local CF structure in disordered graphite fluoride materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational investigations of Dienes defect- and vacancy-induced changes in the electronic and vibrational properties of carbon fiber structural units

Carbon fiber (CF) is a promising lightweight alternative to steel and is of significant interest for energy applications. As CF continues to find new uses and is exposed to new external conditions, a noninvasive method of monitoring its structural integrity is critical. Raman spectroscopy is a commonly used method for this monitoring; however, it is highly inferential, and the interpretation of the data is not always straightforward. In this work, we perform density functional theory (DFT) calculations to investigate changes in the vibrational properties of CF structural units (i.e., graphene and graphite) caused by monovacancy and Dienes defects as a foundation for modeling more complex defects that move our model toward that of realistic CF. Using large computational supercells, we can understand how these defects change the electronic structure and vibrational properties of graphene and graphite for interdefect distances near those of the lower experimental limit. The monovacancy opens an electronic bandgap at the K point. Although no such electronic gap is opened by the Dienes defect, both defects introduce flat defect bands near the Fermi energy. The Dienes defect creates long-range deviations of the phonons, leading to substantial broadening of the highest frequency optical modes in the band structure compared to that of the pristine material. In contrast, the phonon changes caused by the monovacancy are short range, and only minor changes in the band structure or phonon density of states were observed. These findings can assist in the interpretation of experimental results by providing atomic-scale insight into key electronic and vibrational features.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the Representativity of Electrode Microstructure Parameters and Their Electrochemical Response for Lithium Ion Batteries

Lithium-ion battery electrochemical models require an accurate description of the electrodes microstructures to be predictive that can be achieved through nanoscale imaging. Such observations are however limited by their field of view (FOV), as they provide only a subset of the whole electrode volume that does not necessarily represent the whole electrode microstructure heterogeneity, and therefore can bias the microstructure analysis. A microstructure scale electrochemical model was used to investigate lithium plating onset, material non-uniform utilization, and in-plane heterogeneities for an NMC-graphite full cell. To evaluate the representativeness, and thus relevance, of these model predictions, a coupled representativity analysis has been performed on the microstructure parameters and, in a novel way, on the full cell electrochemical response. Electrode microstructure parameters representativeness has been first quantified using the representative volume element (RVE) methodology. The RVE major flaw is that ultimately it can only conclude if a FOV contains representative subvolumes of the FOV, but not if the FOV itself is representative of the electrode volume. Analysis can conclude negatively ('FOV is not representative'), but not positively ('FOV is representative'). One major contribution of this work was to quantify the convergence of the RVE size with the FOV, to actually investigate the FOV representativeness and thus partly remedy this intrinsic limitation. The analysis determined that performing a standard RVE calculation, without exploring its FOV convergence, is likely to strongly underestimate the actual RVE size. The new RVE methodology has been automated in the NREL open-source Microstructure Analysis Toolbox (MATBOX) and is available to the battery community. Representativeness of microstructure parameters is however only an intermediate step, as the end-results of an electrochemical model are performances predictions. Indeed, what is the practical consequence of a given deviation for a microstructure parameter? The microstructure parameter deviation propagations to the 3D microstructure scale electrochemical response have been then quantified for different charge rates. This defines a threshold for the microstructure parameters FOV for a desired maximum deviation of the electrochemical response. Such deviation propagation analysis is analogous to error propagation analysis and is necessary to determine the relevance of microstructure scale model predictions for macroscale predictions. Electrochemical model shows cell representative section areas are increasing with C-rate, due to higher in-plane heterogeneities, indicating larger FOVs are required specifically for fast charge modeling. Therefore, we introduced the novel concept of electrochemical RVE (eRVE) that is a function of the operating conditions (thus defined as a dynamic RVE), with an increasing dependence with the C-rate. Representativity analysis of the investigated cell determined a FOV of 144.4 x 54.4 m2 is large enough to establish a convergence on the representative section areas for low to intermediate C-rate (=2.5C), but not large enough to conclude for higher rates. This work aims to emphasize the importance of representativity analysis for LIB electrode microstructures, as it is required to estimate the error, and thus the relevance, of microstructure parameters intended to be used in macroscale models. The methodology and results can help researchers to select the relevant imaging and associated FOV required to provide accurate enough microstructure parameters.

ADVANCED PROPULSION SYSTEMS↗

AGR-5/6/7 Thermal Model with Non-uniform Gas Gaps

Fuel compact temperatures are a crucial factor in assessing the irradiation performance of tri-structural isotropic fuel particles. In the absence of direct measurement, fuel compact temperatures were calculated using a three-dimensional finite element thermal model, which is subject to simulation uncertainty. The most dominant factor in the uncertainty of calculated fuel temperatures is the gas gap uncertainty due to the nub-to-shell clearance caused by a design error of AGR-5/6/7 capsules. The thermal model was revised to examine the most probable graphite offset position for six different days during the irradiation for Capsules 1 and 2. The analysis varied the offset distance and azimuthal direction at both the top and bottom of the holder. The best-fit offset was estimated based on the minimum root mean square error of the residuals (measured minus calculated) for the operational thermocouples (TCs). From these results, the following conclusions were made: (1) The holder offsets led to slightly lower average temperatures but wider temperature variations (lower minimum and higher peak fuel temperatures) for both Capsule 1 and Capsule 2. (2) During earlier cycles (162A–164B), when numerous TCs were still operational, the best-fit offset distance varied over a specific range for both the top and bottom ([0.002–0.0035 in.] for Capsule 1 and [0.003-0.004 in] for Capsule 2). In contrast, the offset azimuthal direction varied widely, especially for the offset at the bottom of the Capsule 1 holder. This is because holder movement was somewhat constrained at the top of Capsule 1 by the TC leads running through the capsule head and into the holder and by the through tubes in Capsule 2, but the Capsule 1 bottom did not have this type of constraint. (3) During later cycles, when all TCs failed, applying the maximum possible offset of 0.006 in. to the northwest direction for both the top and bottom resulted in a calculated peak fuel temperature of 1557? in Capsule 1 (i.e., a 135? increase from 1422? with zero offset on September 20, 2019 (166A)); the maximum offset of 0.0068 in. to the south for both top and bottom resulted in a calculated peak fuel temperature of 1110°C in Capsule 2 on April 20, 2020 (i.e., a 116? increase from 994? with zero offset (168A)). High peak fuel temperatures in Capsule 1 during Cycle 166A could be the cause of massive particle failure near the end of this cycle. (4) Even though the highest temperature at the tip of Type-N TCs, such as TC-1-7, slightly exceeded 1000?, the temperature along the TC wire reached as high as 1335? assuming an offset of 0.006 in. in the northwest of Capsule 1 holder near the end of Cycle 166A. This temperature significantly exceeds the temperature threshold at which TC degradation is expected to occur, ultimately contributing to considerable particle failures in Capsule 1. For eight Type-N TCs in Capsule 2, the peak TC line temperature was much lower (i.e., 1029°C for TC-2-5), assuming maximum offset of 00068 in. to the south during cycle 168A.

advanced gas reactors↗