## Abstract

Homogenization heat treatment is performed to attain uniformity in microstructure which is helpful to achieve the desired workability and microstructure in final products and, eventually, to gain predictive and consistent performance. Fabrication of low-enriched uranium alloys with 10 wt% molybdenum (U-10Mo) fuel plates involves multiple thermomechanical processing steps. It is well known that the molybdenum homogeneity in the final formed product affects the performance in the nuclear reactor. To ensure uniform homogenization, a statistical method is proposed to quantify and characterize the molybdenum concentration variation in U-10Mo fuel plates by analyzing the molybdenum concentration measurement data from scanning electron microscopy energy dispersive spectroscopy line-scan. Statistical tolerance intervals (TI) are employed to determine the qualification of the U-10Mo fuel plate. We formulate an argument for the minimum number of independent samples to define fuel plate qualification if no molybdenum measurement data are available in advance and demonstrate that the given TI requirements can be equivalently reduced to a sample variance criterion in this application. The outcome of the statistical analysis can be used to optimize casting design and eventually increase productivity and reduce fabrication costs. The statistical strategy developed in this paper can be implemented for other applications especially in the field of material manufacturing to assess qualification requirements and monitor and improve the process design.

## 1 Introduction

The US National Nuclear Security Administration drives the need to develop and design low-enriched uranium (LEU) fuels to progressively reduce and eventually replace high-enriched uranium (HEU) fuels for US high-power civil and research reactors [1]. Uranium alloys with 10 wt% molybdenum (U-10Mo) are identified as the most promising candidate for LEU fuel due to its high intrinsic density and good irradiation behavior [2–4]. The manufacturing of the U-Mo alloy requires multiple complex thermomechanical processes [5,6] including casting, homogenization, hot/cold rolling, annealing, and hot isostatic pressing. Every process is highly challenging and requires a tremendous amount of planning and preparation. Fuel plate quality after casting and homogenization is fundamental to ensure high quality of the final U-Mo product because it can greatly affect the subsequent microstructure evolution and eventually influence the fuel performance in a reactor. To determine the qualification of a given fuel plate, a statistical analysis of the chemical composition measured using an energy dispersive spectroscopy (EDS) technique is a reliable, efficient, and cost-effective. Statistical methods have been widely used to characterize the quality of processing products in many applications [7–10]. For this study, tolerance interval (TI) calculations will be used to determine the quality of as-cast plates, identify potential quality problems, and eventually enhance the quality and efficiency of the entire fabrication process, thereby reducing fabrication costs.

Previous works [6,11] reveal that the distribution of molybdenum in fuel plates has a significant impact on the quality and performance of U-Mo fuel. In principle, the overall average molybdenum concentration should be 10% throughout the entire plate. However, the variation in molybdenum concentrations after casting is significant due to molybdenum segregation and a homogenization treatment is necessary to produce a more uniform molybdenum distribution. Therefore, statistical analysis is needed to quantify the variation in molybdenum concentrations that determines the qualification of the fuel plate.

Two common measures of product quality compare the confidence interval (CI) for a parameter and/or a TI to specified limits [12,13]. A CI specifies the statistical interval that might contain the true value of population parameter, e.g., mean or variance, with a given probability. The width of the confidence interval will approach zero with the increasing sample size. A TI bounds a specified proportion of a sampled population at a given CI. The width of TI is also dependent on the number of samples and variance.

Extensive work [14] has been performed for quality control in a broad range of applications, where the TI is used as a quantitative metric for assessing the product quality. However, no such attempt has been published on characterizing the molybdenum concentration variation in U-Mo fuel plates for fuel qualification. A recent statement of qualification specification for the U-Mo fuel plate has been proposed [15] and quoted herein as follows: “A sufficient number of samples shall be randomly taken from a sufficient number of fuel plates randomly selected from a fuel plate lot to allow statistical determination of 95% confidence interval that 95% of the U-Mo has 10.0 ± 1.0 wt% molybdenum microscopically throughout the U-Mo.” This statement prompts one to use TI to quantify the molybdenum concentration variation and ensure product quality. This specification can be equivalently interpreted as follows: 95% of the total U-Mo population within a plate must fall between 9% and 11% molybdenum wt% with a confidence interval level of 95%. However, there is no guideline on the “sufficient number of samples” in the statement. As in other sampling scenarios, there is a tradeoff between the number of high-quality measurements that can be obtained and the total time cost for data collection. Specifically, a single molybdenum wt% measurement from EDS line-scans takes about 30 s to 1 min, and a large number of measurements are required to meet statistical specifications, given inherent molybdenum variability within a plate. In addition to determining the fuel plate qualification, this study also provides insights into an inverse problem, i.e., how to determine the optimal number of measurements required to qualify the given plate.

Fuel plate qualification can be specified in terms of the sample variance rather than the calculation of tolerance limits, which is a simplified metric that may be applied given reasonable assumptions. One of these assumptions is that the collected samples must be independent of one another, which ensures that the selected data properly represent the characteristics of the entire fuel plate. The product is considered qualified only if the simplified metrics are satisfied and the assumptions can be verified. The final results of the statistical analysis can be used to optimize the casting design, refine casting specifications, and eventually enhance the casting productivity. The statistical strategy developed in this paper can be implemented in other applications, especially in the field of material manufacturing, to define reasonable qualification requirements and provide insights to improve the process design.

The rest of the paper is organized as follows: Section 2 describes the statistical methodology including a comprehensive introduction of TI approaches and key mathematical formulas. Section 3 presents the TI implementation on the specified U-Mo qualification specification, and two examples are analyzed as demonstration. Section 4 introduces our recommendation on the minimum number of molybdenum wt% measurements and the sample variance criterion to determine the fuel plate qualification. The conclusions are summarized in Sec. 5.

## 2 Methodology

**= (**

*Y**Y*

_{1},

*Y*

_{2},…,

*Y*

_{N}) be

*N*-independent random samples drawn from a continuous cumulative distribution function,

*F*. In addition, let an interval [

*L*(

**),**

*Y**U*(

**)] based on the sample vector**

*Y***be a two-sided TI for**

*Y**F*such that for any given 0 <

*γ*< 1 and 0 <

*P*< 1, the following probability expression is satisfied:

*γ*and

*P*define a general two-sided TI that has a confidence level of $\gamma \xd7100%$ containing at least $P\xd7100%$ of the population. If

*L*(

**) = −∞ and**

*Y**U*(

**) < +∞, then the previous two-sided interval becomes [−∞,**

*Y**U*(

**)] and can be defined more properly as a one-sided TI (**

*Y**γ*,

*P*) with an upper tolerance bound. Similarly, if

*L*(

**) > −∞ and**

*Y**U*(

**) = +∞, the interval [**

*Y**L*(

**), +∞] is a one-sided TI (**

*Y**γ*,

*P*) with a lower tolerance bound only [16].

**= (**

*Y**Y*

_{1},

*Y*

_{2}, …,

*Y*

_{N}) is independently drawn from a standard normal distribution

*N*(

*μ*,

*σ*

^{2}) with the population mean

*μ*and the population variance

*σ*

^{2}. We can use $Y\xaf$ and

*s*to denote sample mean and sample variance, respectively. The mathematical expressions for $Y\xaf$ and

*s*can be written as

*k*

_{1}is the factor to ensure that at least a proportion $P\xd7100%$ of the entire population with a confidence interval level of $\gamma \xd7100%$ is satisfied based on the sample of size

*N*. The value of

*k*

_{1}can be calculated based on the inverse cumulative distribution function for a non-central

*t*distribution [17]:

*t*

_{N−1,γ,δ}is the critical value (

*γ*× 100 percentile), i.e.,

*F*(

*t*

_{N−1,γ,δ}) =

*γ*, for a non-central

*t*distribution with

*N*− 1 degrees-of-freedom and the non-centrality parameter $\delta =zPN$. Here,

*z*

_{P}is the critical value of the standard normal distribution associated with a cumulative probability of $P\xd7100%$. For example, Fig. 1

^{,}

^{2}shows that the accumulated probability for −∞ ≤

*x*≤ 1 is 84%, which is equivalent to $Pr[x\u2264zP=1]=84%$. The value of

*z*

_{P}can be found in Ref. [18] for any given

*P*. Some commonly used

*P*and the associated

*z*

_{P}values are listed in Table 1.

*k*

_{2}factor can be obtained by solving the following implicit nonlinear integral equation [21,22]:

*N*− 1 degrees-of-freedom; $\chi 1,P,\delta 2$ is the $P\xd7100%$ percentile of the chi-squared distribution with one degree-of-freedom and non-centrality parameter

*δ*=

*z*

^{2}, where

*z*is a dummy variable used for integration. Equation (9) is quite complicated to solve directly. Moreover, the following approximation can be used to simplify the calculation of the two-sided

*k*

_{2}factor [23]:

*N*versus the relative error for $P=\gamma =0.95$. The approximations first slightly underpredict and then overpredict the true value of $k2$ with the increasing

*N*, but the relative error is always within the range of ±0.36% and gradually approaches zero with increasing

*N*. These results indicate that the $k2$ values obtained from Eq. (10) can be considered as good approximations to represent the true solutions for

*N*approaching 100.

## 3 Statistical Analysis for Qualification of U-10Mo

Based on the specification statement for U-10Mo, we can set $P=\gamma =0.95$. In addition, this requirement also defines the lower and upper bounds for 95/95 tolerance interval as 9% and 11%, respectively. Note that an average concentration $Y\xaf=10%$ is assumed for all samples and is assumed known without error so Eqs. (4), (5), and (8) become a function of $k2$ and *s*. The number of samples, $N$, namely the number of measurements of molybdenum concentration required when performing EDS scan experiments for each fuel plate is not explicitly specified and will be determined in this study. Molybdenum concentration will drop significantly where the scan intersects a carbide or secondary phase spot. To avoid the bias due to carbides, only molybdenum weight percentage between 7% and 13% is considered as effective measurements and used in the plate qualification analysis. The effect of the distribution of carbides and/or secondary phase spots on the fuel plate quality is currently not considered in this study.

For any independently drawn molybdenum concentration vector of length $N$ which is close to the neighborhood of 100, Eq. (10) can be applied to calculate the $k2$ factor. With $P=0.95$, we can find $z(1+P)/2=z0.975=1.96$. The values of $\chi 1\u2212\gamma ,N\u22121,02=\chi 0.05,N\u22121,02$ and $k2(N,\gamma ,P)=k2(N,0.95,0.95)$ become a function of *N* only and are calculated once $N$ is provided. Equation (8) can be used to determine lower and upper bounds (*L* and *U*) for the observed vector of molybdenum concentrations. The fuel plate is determined to be qualified in terms of the molybdenum distribution if both the upper and lower bounds satisfy the specification requirement, i.e., $L\u22659%$ and $U\u226411%$. Otherwise, the plate will be rejected.

*r*for lag

*j*is defined as

A unit lag distance (*j* = 1) in our study is 1 *µ*m which represents the space between the adjacent molybdenum concentration measurements. By observing the line-scan data listed in Fig. 4(b), it is evident that the measurements vary in a fashion that is not purely random, and it is a strong indication of data correlation. This trend becomes even more obvious if carbides are excluded. The findings from Fig. 6, which plots the ACF results for Plate 1, are consistent with the observations from the measured data. In particular, the critical lag distance indicating no data intercorrelation increases from ∼15 *µ*m to ∼100 *µ*m after neglecting carbides. Critical lag distance is an important parameter for determining the repeating data pattern. If multiple data are collected within the critical lag distance for data analysis, then the behavior and the characteristic of the entire data set cannot be properly represented. In other words, the measurement must be taken at the interval of critical lag distance to ensure data independence and to avoid data similarity. The critical lag distance is approximated at where the ACF values start to vary between the dashed lines which specifies a 95% CI of data independence.

To simulate the distribution of *k*_{2} as a function of the number of samples, we randomly sampled *N _{s}* samples from the observed data vector, repeating the experiment 1000 times for each value of

*N*in Table 3. The last column in this table shows the results for

_{s}*N*=

_{s}*N*= 946. For any given

*N*, the

_{s}*k*

_{2}for 95/95 TI can be calculated from Eq. (10) and appear in the third row of Table 3. The mean value of sample standard deviation,

*s*, listed in the fourth row slightly increases with

*N*and then remains constant for

_{s}*N*greater than 20.

_{s}Equation (12) shows that sample variance will gradually increase and approach the population variance (constant) with increasing *N*—our findings are consistent. The next row of Table 3 shows the mean value of *k*_{2}*s* (multiplication of *k*_{2} and *s*), which decreases with increasing *N _{s}*. The last two rows list the corresponding lower and upper bounds for each sample size

*N*.

_{s}*L*< 9% and

*U*> 11% are observed for all cases, which indicates that this plate will not meet specified qualification and the molybdenum spatial distribution is not sufficiently homogeneous.

Next, Fig. 7 shows a second example plate (Plate 2) with its SEM and three EDS lines (250 measurements on each solid line). The statistics for varying *N _{s}* on plate qualification for Plate 2 is listed in Table 4. The minimum value of

*N*to make $L>9%andL<11%$ and reach the qualification requirement is 10 for Plate 2. The corresponding column is highlighted in bold. Therefore, the molybdenum distribution in Plate 2 is determined to be sufficiently homogeneous and hence qualified for the subsequent fabrication processes.

_{s}A spatial independence study using ACF is also conducted for Plate 2. The critical lag distance for both with and without particles scenarios is similar and close to 5 *µ*m (Fig. 8). The ACF results are in line with the observation from Fig. 7, which shows no clear pattern of data distribution except white noise.

## 4 Recommendation on the Minimum Number of Molybdenum wt% Measurements to Determine Fuel Plate Qualification

It was shown in the second example that ten measurements are already sufficient for Plate 2 to be qualified. However, this number can be different for other areas within a plate (or just a different plate) with different degrees of homogenization. Next, we provide a general recommendation for the optimal number of measurements to determine fuel plate qualification if no available measurement data exist in advance. As mentioned earlier, Eq. (10) indicates that *k*_{2}(*N*, *γ*, *P*) is only a function of sample size *N* with fixed *P* = *γ* = 0.95. A plot of *k*_{2} versus *N* in Fig. 9 shows that *k*_{2} gradually decreases (lower curve) but its slope slowly increases (top curve) with the increasing *N*. The value of *k*_{2} eventually approaches an asymptote of ∼2. A gradient of *k*_{2} greater than −0.1% can be considered as a mathematical indication of stability. The corresponding value *N* at $dk2/dN=\u22120.1%$ is 133, and the associated measurement time (∼2 h) is experimentally affordable. Therefore, it is reasonable to select *N* = 133 as the suggested number of measurements to determine the fuel plate qualification, resulting in *k*_{2} = 2.19.

Based on these findings, a suggested procedure of data collection is described as follows to make sure that the selected data are not autocorrelated and are independent:

Collect the initial 133 EDS line-scan data with a normal 1

*µ*m spacing between each measurement.Conduct the ACF analysis and then determine the critical lag distance.

Determine the number of valid collections among the original 133 data and then make additional collections based on the critical lag distance obtained from step 2.

Discard the measurements if the molybdenum weight percentage is outside the range of 7–13% to avoid the influence of carbides.

For any given fuel plate that meets the above-mentioned criterion, it can be considered as a qualified piece and be supplied to the subsequent processes only if all 133 collected data satisfy *s* ≤ 0.46, which is equivalent to *k*_{2}*s* ≤ 1. Otherwise, the plate is disqualified.

## 5 Conclusion

Characterizing the quality of U-Mo in casting is important and also challenging. A statistical approach based on the concept of TI is implemented to determine fuel plate qualification. This study also recommends a minimum number of required molybdenum wt% measurements to represent the distribution behavior of molybdenum for the entire fuel plate when other data are not available to determine the number of samples. Then, the given qualification specification can be simplified using only sample variance under reasonable assumptions. In particular, for the fuel plates discussed here, at least 133 data points across the line-scans are sufficient. Data points wherein the line-scan intersects a carbide or a secondary phase particle discarded. The collection of data must be independent without any correlation to assure that the collected measurements are a valid representation of the entire fuel plate. For sample standard deviation values equal to or less than 0.46 and means equal to 10 wt% molybdenum the given plate will ideally be qualified, whereas for standard deviation values higher than 0.46 the material is not qualified. By monitoring the fuel plate quality using the proposed method, the efficiency of the plate qualification is greatly enhanced, and fabrication costs are reduced. The developed methodology can be extended to other fields and serve as a metric for quality control.

## Footnote

## Acknowledgment

The authors would like to thank Shelley Carlson, Mark Rhodes, and Alan Schemer-Kohrn for helping with the metallography samples and capturing the SEM images. Pacific Northwest National Laboratory is operated by Battelle Memorial Institute for the United States Department of Energy under Contract DE-AC06-76RLO1830.

## Funding Data

U.S. Department of Energy (Funder ID: 10.13039/100000015).

National Nuclear Security Administration’s Office of Material Management and Minimization (Funder ID: 10.13039/100006168).