Search NASA⌕ Search

SEARCH · Search NASA

Results for “Probability estimates”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Risk-Informed MultiPphysics BestEestimate Plus Uncertainties (BEPU) Application Development of RELAP5-3D Perturbation Model

The United States nuclear industry is facing a strong challenge to ensure maximum safety while enhancing economic benefit. Safety is a key parameter to all aspects related to light water reactor (LWR) nuclear power plants (NPPs), especially cost savings. Since the goal is to extend the lifetimes of these NPPs, the traditional deterministic safety concept may not guarantee a current economic asset. The Light Water Reactor Sustainability (LWRS) Program has been promoting a wide range of research and development (R&D) in this field to maximize the safety, economics, and performance of these NPPs through improved scientific understanding. One of the best practices to achieve this goal is to identify and optimize safety margins, which can lead to cost reduction. To do this, under the LWRS framework, the Risk-Informed Systems Analysis (RISA) Pathway will focus on the optimization of safety margin and minimization of uncertainties to ensure both safety and economics at the highest level. The RISA Pathway will provide enhanced capabilities for analyzing and characterizing LWR systems performance by developing and demonstrating methods, tools, and data to enable risk-informed margins management (RIMM). The goals of the RISA Pathway are twofold: (1) deploy the risk-informed tools and methods that enable better representation of safety margins and factors that contribute to cost and safety; and (2) conduct advanced risk assessment applications with industry to support margin management strategies that enable more cost-effective plant operation. The tools and methods provided by the RISA Pathway will support effective margin management for both active and passive safety systems, structures, and components (SSC) of an NPP. The tools and methods used in the RISA Pathway should have high confidence and highest technical maturity for and implementation to industry at its current setting. They should also have a capability to support risk-informed decision making for both probabilistic and deterministic elements of safety. The RISA Pathway will, therefore, perform a comprehensive assessment of verification and validation (V&V) status of RISA Toolkit to enhance credibility RISA Toolkit which be used by industry. This report summarizes RELAP5-3D development activities for best estimated and uncertainty (BEPU) capability to support ongoing RISA Pathway pilot projects for risk informed uncertainty quantification applications. Work scope includes (1) development of initial perturbation model for selected closure laws and (2) testing of suitable probability distribution functions (PDF) for the selected correlations. Developed BEPU modules are tested with selected separate effect cases and compared with conventional thermal-hydraulics model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Use of Tracer Elements for Estimating Community Exposure to Marcellus Shale Development Operations

Since 2009, unconventional natural gas development (UNGD) has significantly increased in Appalachia’s Marcellus Shale formation. Elevations of fine particulate matter <2.5 µm (PM2.5), have been documented in areas surrounding drilling operations during well stimulation. Furthermore, many communities are experiencing increased industrial activities and probable UNGD air pollutant exposures. Recent studies have associated UNGD emissions with health effects based on distances from well pads. In this study, PM2.5 filter samples were collected on an active gas well pad in Morgantown, West Virginia, and three locations downwind during hydraulic stimulation. Fine particulate samples were analyzed for major and trace elements. An experimental source identification model was developed to determine which elements appeared to be traceable downwind of the UNGD site and whether these elements corresponded to PM2.5 measurements. Results suggest that 1) magnesium may be useful for detecting the reach of UNGD point source emissions, 2) complex surface topographic and meteorological conditions in the Marcellus Shale region could be modeled and confounding sources discounted, and 3) well pad emissions may be measurable at distances of at least 7 km. If shown to be more widely applicable, future tracer studies could enhance epidemiological studies showing health effects of UNGD-associated emissions at ≥15 km.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

A Southern Photometric Quasar Catalog from the Dark Energy Survey Data Release 2

We present a catalog of 1.4 million photometrically selected quasar candidates in the southern hemisphere over the ~5000 deg 2 Dark Energy Survey (DES) wide survey area. We combine optical photometry from the DES second data release (DR2) with available near-infrared (NIR) and the all-sky unWISE mid-infrared photometry in the selection. We build models of quasars, galaxies, and stars with multivariate skew-t distributions in the multidimensional space of relative fluxes as functions of redshift (or color for stars) and magnitude. Our selection algorithm assigns probabilities for quasars, galaxies, and stars and simultaneously calculates photometric redshifts (photo-z) for quasar and galaxy candidates. Benchmarking on spectroscopically confirmed objects, we successfully classify (with photometry) 94.7% of quasars, 99.3% of galaxies, and 96.3% of stars when all IR bands (NIR YJHK and WISE W1W2) are available. The classification and photo-z regression success rates decrease when fewer bands are available. Our quasar (galaxy) photo-z quality, defined as the fraction of objects with the difference between the photo-z z p and the spectroscopic redshift z s , |Δz| ≡ |z s - z p |/(1 + z s ) ≤ 0.1, is 92.2% (98.1%) when all IR bands are available, decreasing to 72.2% (90.0%) using optical DES data only. Our photometric quasar catalog achieves an estimated completeness of 89% and purity of 79% at r < 21.5 (0.68 million quasar candidates), with reduced completeness and purity at 21.5 < r ≲ 24. Among the 1.4 million quasar candidates, 87,857 have existing spectra, and 84,978 (96.7%) of them are spectroscopically confirmed quasars. Finally, we provide quasar, galaxy, and star probabilities for all (0.69 billion) photometric sources in the DES DR2 coadded photometric catalog.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Projecting U.S. coastal storm surge risks and impacts with deep learning

Storm surge is one of the deadliest hazards posed by tropical cyclones (TCs), yet assessing its current and future risk is difficult due to the phenomenon’s rarity and physical complexity. Recent advances in artificial intelligence applications to natural hazard modeling suggest a new avenue for addressing this problem. We develop a deep learning storm surge model to efficiently estimate coastal surge risk in the United States from 900 000 synthetic TC events, accounting for projected changes in TC behavior and sea levels. The derived historical 100 year surge (the event with a 1% yearly exceedance probability) agrees well with historical observations and other modeling techniques. When coupled with an inundation model, we find that heightened TC intensities and sea levels by the end of the century result in a 50% increase in population at risk. Key findings include markedly heightened risk in Florida, and critical thresholds identified in Georgia and South Carolina.

RAFT↗
DOE OSTI2023

Stochastic Price Generation for Evaluating Wholesale Electricity Market Bidding Strategies

This work presents a novel method for generating electricity price scenarios from statistical properties of past electricity prices using a hybrid statistical and reduced-form stochastic model. Previous work in applying stochastic differential equations (SDE) to model electricity prices has focused on daily average prices. To extend stochastic price generation methods to hourly or sub-hourly pricing, we address several weaknesses in the state-of-the-art: (1) we replace the mean-reversion component of the SDE with an ARIMA process that is better able to characterize the daily and weekly trends; (2) we extend the price-spike, or jump process to account for conditional probabilities of price spikes occurring in consecutive time steps by replacing the traditional Poisson process for modeling jumps with a generalized point process model inspired by brain neuron models; and (3) we replace the traditional method of estimating spike intensity with empirical variance with a Markov process based on observed price spike intensity transitions. The method is demonstrated with electricity prices from the US ERCOT market and a use-case example is provided for bidding an energy storage unit into the day-ahead and real-time energy markets of ERCOT using stochastic optimization methods. Results show that the the synthetic price model out performs a (naive) persistence forecast model by resulting in 24% to 47% more in profits over 168 simulated days.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Radiological Impact of 2023 Operations at the Savannah River Site

This report presents the environmental dose assessment methods and the estimated potential doses to the public from 2023 Savannah River Site (SRS) air and liquid radioactive releases. Also documented are potential doses from special-case exposure scenarios, such as the consumption of wildlife or goat milk. Dose to the Offsite Representative Person The 2023 dose to the offsite representative person from SRS liquid releases was 0.14 mrem and from SRS air releases it was 0.016 mrem. To show compliance with the U. S. Department of Energy (DOE) all pathway dose standard of 100 mrem/yr, SRS conservatively adds these two doses for a total representative person dose of 0.16 mrem which is 0.16% of the DOE standard. Sportsman Doses Onsite Hunter: SRS conducts annual hunts to control onsite deer and feral hog populations. The estimated dose from consuming harvested deer or hog meat is determined for every onsite hunter. During 2023, the maximum potential dose an onsite hunter received was 9.42 mrem, or 9.42% of DOE’s 100 mrem/yr all pathway dose standard. Creek Mouth Fisherman: SRS estimated the maximum potential dose from fish consumption at 0.17 mrem from bass collected at the mouth of Lower Three Runs. This dose is 0.17% of the DOE standard. SRS bases this hypothetical dose on the low probability scenario that, during 2023, a fisherman consumed 24 kg (53 lbs) of bass caught exclusively from the mouth of Lower Three Runs. Release of Material Containing Residual Radioactivity SRS did not release any real property (land or buildings) in 2023. SRS unconditionally released a total of 13,324 items of personal property (such as tools) from radiological areas in 2023. Most of these items did not leave the Site. However, all of these items required no additional radiological controls post-survey as they met DOE Order 458.1 release criteria. Radiation Dose to Aquatic and Terrestrial Biota SRS conducts screening evaluations of plant and animal doses for aquatic and terrestrial ecosystems. For 2023, all SRS aquatic system locations passed the initial (Level 1) screenings and no further assessments were required at those locations. For the land-based systems evaluation, SRS performed initial screenings using concentration data from the five onsite radiological soil sampling locations. Typically, SRS collects and analyzes only one soil sample per year from each location. For 2023, all land-based locations passed their initial (Level 1) pathway screenings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluating the Effectiveness of a Detection and Deterrent System in Reducing Golden Eagle Fatalities at Operational Wind Facilities

The Renewable Energy Wildlife Institute (REWI) was appointed as the prime awardee of DOE award number DE-EE0007883 to lead a team of scientists, wind developers, and technology manufacturers toward the overarching goal of evaluating the effectiveness of the current DTBird system in minimizing the risk of golden eagles (Aquila chrysaetos) and other large soaring raptors from approaching the rotor-swept zone (RSZ) of operating wind turbines. As part of this goal, the team set out to 1) quantify the expected reduction in collision risk for golden eagles from operation of the detection and deterrence modules in a manner that supports the approach used by the U.S. Fish and Wildlife Service (USFWS) to assess and credit facility operators for their efforts to minimize predicted collision fatalities and 2) provide information to help improve the technology to maximize its effectiveness. DTBird is an automated detection and audio deterrent system created by the Spanish company Liquen, designed to discourage birds from entering the RSZ of spinning wind turbines. The system uses cameras to automatically detect airborne targets of interest, records each such event in an online database, and triggers a warning signal (loud sound) if the tracked object has moved close to the turbine. If the object moves even closer to the RSZ, a more aggressive dissuasion signal is broadcast. To meet our objectives, the team conducted a two-year experiment at the Goodnoe Hills wind facility in Washington state, in which 14 turbines were outfitted with DTBird units. Daily, each DTBird-equipped turbine was randomly assigned to a control or treatment group. Treatment turbines operated with DTBird running as intended—broadcasting warning or deterrent signals when DTBird detected a target within range. On control turbines, no sound signals were broadcast if a moving target triggered the DTBird system. The team also flew unmanned aerial vehicles (UAVs) designed to coarsely mimic the general size, weight, and coloration of golden eagles in programmed flight transects across DTBird detection ranges to quantify DTBird’s ability to detect intended targets and to evaluate factors that influence the probability of detection and DTBird’s response distances. Additionally, the team evaluated the behavioral responses of in situ eagles exposed to spinning turbines alone (visual and sound influences) versus spinning turbines plus broadcasted DTBird audio deterrents, to estimate the effectiveness of deterrence by the DTBird system. The data and results from these investigations were combined with those from a pilot study conducted at the Manzana Wind Power Project in California to better evaluate DTBird’s effectiveness across different landscapes.

17 WIND ENERGY↗
DOE OSTI2023

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

Fuel compact temperatures are a crucial factor in assessing irradiation performance of the Tristructural 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 uncertainty of calculated fuel temperature is the gas gap uncertainty due to the nub-to-shell clearance because of fabrication error. A revision to the thermal model was made to examine the most probable offset position for four different dates of interest. The offset position options varied in magnitude and azimuthal direction of offset of the holder top and bottom position. The best-fit offset of the holder is estimated based on the minimum root mean square error of residuals for all operational thermocouples (TCs). From these results, the following conclusions were made: 1) During earlier cycles (162A – 164A) when numerous TCs were still operational, the best-fit offset distance varied in range [0.002 – 0.0035 in.] for both the top and bottom holder while offset azimuthal direction varied widely, especially for the holder bottom. 2) Holder offset led to slightly lower Capsule 1 average temperature, but wider temperature variation (lower minimum and higher peak fuel temperatures). 3) The maximum offset of 0.006 in. to the northwest direction for both the top and bottom holder yielded a much higher temperature than the original peak fuel temperature, which increased from 1422 to 1557 ? (or 135 ? increase). This could be one of the reasons for massive particle failure near the end of Cycle 165A. 4) Even though, the highest temperature at the TC-1-7 tip reached slightly higher than 1000?, the temperature along the TC wire reached as high as 1335? assuming a holder offset of 0.006 in. to the northwest. Consequently, at this temperature nickel from this TC wire could cause particle failures in Capsule 1.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Energy Budget Constraints on the Time History of Aerosol Forcing and Climate Sensitivity

Abstract An observationally constrained time series of historical aerosol effective radiative forcing (ERF) from 1750 to 2019 is developed in this study. We find that the time history of aerosol ERFs diagnosed in CMIP6 models exhibits considerable variation and explore how the time history of aerosol forcing influences the probability distributions of present‐day aerosol forcing and emergent metrics such as climate sensitivity. Using a simple energy balance model, trained on CMIP6 climate models and constrained by observed near‐surface warming and ocean heat uptake, we derive estimates for the historical aerosol forcing. We find 2005–2014 mean aerosol ERF to be −1.1 (−1.8 to −0.5) W m −2 relative to 1750. Assuming recently published historical emissions from fossil fuel and industrial sectors and biomass burning emissions from SSP2‐4.5, aerosol ERF in 2019 is −0.9 (−1.5 to −0.4) W m −2 . There is a modest recovery in aerosol forcing (+0.025 W m −2 decade −1 ) between 1980 and 2014. This analysis also gives a 5%–95% range of equilibrium climate sensitivity of 1.8°C –5.1°C (best estimate 3.1°C) with a transient climate response of 1.2°C –2.6°C (best estimate 1.8°C).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Sparsifying priors for Bayesian uncertainty quantification in model discovery

We propose a probabilistic model discovery method for identifying ordinary differential equations governing the dynamics of observed multivariate data. Our method is based on the sparse identification of nonlinear dynamics (SINDy) framework, where models are expressed as sparse linear combinations of pre-specified candidate functions. Promoting parsimony through sparsity leads to interpretable models that generalize to unknown data. Instead of targeting point estimates of the SINDy coefficients, we estimate these coefficients via sparse Bayesian inference. The resulting method, uncertainty quantification SINDy (UQ-SINDy), quantifies not only the uncertainty in the values of the SINDy coefficients due to observation errors and limited data, but also the probability of inclusion of each candidate function in the linear combination. UQ-SINDy promotes robustness against observation noise and limited data, interpretability (in terms of model selection and inclusion probabilities) and generalization capacity for out-of-sample forecast. Sparse inference for UQ-SINDy employs Markov chain Monte Carlo, and we explore two sparsifying priors: the spike and slab prior, and the regularized horseshoe prior. UQ-SINDy is shown to discover accurate models in the presence of noise and with orders-of-magnitude less data than current model discovery methods, thus providing a transformative method for real-world applications which have limited data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Characterizing a supernova’s standing accretion shock instability with neutrinos and gravitational waves

Here, we perform a novel multi-messenger analysis for the identification and parameter estimation of the Standing Accretion Shock Instability (SASI) in a core collapse supernova with neutrino and gravitational wave (GW) signals. In the neutrino channel, this method performs a likelihood ratio test for the presence of SASI in the frequency domain. For gravitational wave signals we process an event with a modified constrained likelihood method. Using simulated supernova signals, the properties of the Hyper-Kamiokande neutrino detector, and O3 LIGO Interferometric data, we produce the two-dimensional probability density distribution (PDF) of the SASI activity indicator and calculate the probability of detection P D as well as the false identification probability P FI . We discuss the probability to establish the presence of the SASI as a function of the source distance in each observational channel, as well as jointly. Compared to a single-messenger approach, the joint analysis results in P D (at P FI = 0.1) of SASI activities that is larger by up to ≈ 40% for a distance to the supernova of 5 kpc. We also discuss how accurately the frequency and duration of the SASI activity can be estimated in each channel separately. Our methodology is suitable for implementation in a realistic data analysis and a multi-messenger setting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A data science approach for analysis and reconstruction of spinodal-like composition fields in irradiated FeCrAl alloys

A statistical method for the analysis of continuously distributed data representative of composition fluctuations in irradiated FeCrAl alloys acquired using Energy Dispersive X-ray Spectroscopy (EDS) method is presented. Using probability distribution functions, direct and cross-covariances between the elemental compositions, the effects of alloy composition and irradiation dose were investigated on the spatial distribution and length scale of composition fluctuations at the nanoscale. We have observed that, for neutron-irradiated FeCrAl alloys, the distribution of Fe and Cr followed a left-skewed and right-skewed distribution, respectively for all (average) alloy compositions and irradiation doses. The analysis also revealed enhanced spatial gradients in the elemental compositions at higher irradiation dose. Direct and cross-covariance estimates of the experimental data were also utilized for reconstruction of composition data through fitting it to a parametric form of the covariance functions. Linear Model of Coregionalization was used to determine the parameters of the covariance functions. Subsequently, a spectral method was utilized for simulating a realization of the alloy compositions. Close correspondence was observed between the experimental and the reconstructed data which was analyzed using probability distribution functions and covariance functions. Composition space of the experimental and reconstructed data and dislocation velocities as a function of applied stress and line directions over the entire composition maps were also examined.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Renewable Energy Potential Model: Hawaii Geothermal Supply Curves

This dataset extends the development of the Renewable Energy Potential (reV) model to include geothermal energy, with a specific focus on Hawaii. Provided here are the results of two scenarios that were modeled for geothermal energy in Hawaii: binary enhanced geothermal systems (EGS) at a depth of 2.5 km and hydrothermal binary systems at a depth of 1.5 km. The resource data for both scenarios were derived from Lautze and Haskins (2024) using an exponential method. The PFA probability of heat map was used as a look up table for which temperature gradient to use (Lautze and Haskins, 2024). The dataset provides geospatial and techno-economic details for evaluating geothermal energy potential. It includes spatial coordinates, estimated capacity factors, developable area, resource potential, and annual energy production metrics. Economic details such as levelized cost of electricity (LCOE), site development costs, transmission costs, and fixed-charge rates are also included. The reV model, originally developed for wind and solar energy, incorporates these variables to evaluate deployment constraints related to land use, environmental and cultural factors, and grid integration.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Defender Policy Evaluation and Resource Allocation against MITRE ATT&CK Data and Evaluations

Protecting against multi-step attacks of uncertain duration and timing forces defenders into an indefinite, always ongoing, resource-intensive response. To effectively allocate resources, a defender must be able to analyze multi-step attacks under assumption of constantly allocating resources against an uncertain stream of potentially undetected attacks. To achieve this goal, we present a novel methodology that applies a game-theoretic approach to the attack, attacker, and defender data derived from MITRE´s ATT&CK ® Framework. Time to complete attack steps is drawn from a probability distribution determined by attacker and defender strategies and capabilities. This constraints attack success parameters and enables comparing different defender resource allocation strategies. By approximating attacker-defender games as Markov processes, we represent the attacker-defender interaction, estimate the attack success parameters, determine the effects of attacker and defender strategies, and maximize opportunities for defender strategy improvements against an uncertain stream of attacks. This novel representation and analysis of multi-step attacks enables defender policy optimization and resource allocation, which we illustrate using the data from MITRE´ s APT3 ATT&CK ® Framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Natural Hazard Forecast Alert Grid Risk System

Weather events cause most power outages. Often, we even get notifications on our phones to take cover or be prepared for an imminent event. If electric grid utilities had a similar warning that also included probable scenarios and the equipment involved, they could prepare and minimize the effects. Idaho National Laboratory had a project with the U.S. Department of Energy’s Cybersecurity, Energy Security, and Emergency Response program to develop a grid alert application that receives messages from the existing emergency alert system, filters and determines components possibly affected by the emergency event, calculates probable scenarios using MASTERRI (Modeling And Simulation for Targeted Reliability and Resilience Improvement). For high-risk events, the application can then send alert links to subscribed electric distribution utility operations staff to allow them to see and evaluate the scenarios and the impact in a web based interactive map tool. This proof of concept application used data from utilities and organizations, such as the international regulatory body North American Electric Reliability Corporation, which have complied historical failure data of elements that comprise the U.S. electric grid. Nominal failure rates are obtained from this data. To make this tool possible, estimated failure rates were calculated for different component types given the alert type, severity, and location. Historic weather-related grid element failures were correlated with historic weather events from the Integrated Public Alert & Warning System. These correlated events and failures are used along with Bayesian updates from the historical norms to provide a modified failure rate for grid elements in the alert areas and calculate probable scenarios. Working with an industry collaborator, actual grid models and data were used for demonstration cases. This report outlines the work performed for this project.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Emergent quantum mechanics at the boundary of a local classical lattice model

We formulate a model in which quantum mechanics emerges from classical mechanics. Given a local Hamiltonian H acting on n qubits, we define a local classical model with an additional spatial dimension whose boundary dynamics is approximately—but to arbitrary precision—described by Schrödinger's equation and H. The bulk consists of a lattice of classical bits that propagate towards the boundary through a circuit of stochastic matrices. The bits reaching the boundary are governed by a probability distribution whose deviation from the uniform distribution can be interpreted as the quantum-mechanical wave function. Bell nonlocality is achieved because information can move through the bulk much faster than the boundary speed of light. Finally, we analytically estimate how much the model deviates from quantum mechanics, and we validate these estimates using computer simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A Bootstrap Approach to Quantifying Reliabilities and Uncertainties of Complex Systems

While much work has been done in analyzing margins and uncertainties at the component level, a gap exists in NNSA methodology relating component level reliabilities and uncertainties to system level reliability and uncertainty. This paper shows how component level reliability data can be combined via a bootstrap analysis to estimate system level reliabilities and uncertainties. The performance of the bootstrap for this problem is validated through simulation studies. This paper extends the original work of Crowder by including sensitivity analyses related to changes in samples sizes and component failure probabilities. The use of the bootstrap as a decision-making tool is thus developed by quantifying the effect of such changes on system reliability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

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↗