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 [24]. 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 [710]. 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

The general definition of a TI can be described as follows. Let Y = (Y1, Y2,…, YN) be N-independent random samples drawn from a continuous cumulative distribution function, F. In addition, let an interval [L(Y), U(Y)] based on the sample vector Y be a two-sided TI for F such that for any given 0 < γ < 1 and 0 < P < 1, the following probability expression is satisfied:
(1)
where γ and P define a general two-sided TI that has a confidence level of γ×100% containing at least P×100% of the population. If L(Y) = −∞ and U(Y) < +∞, then the previous two-sided interval becomes [−∞, U(Y)] and can be defined more properly as a one-sided TI (γ, P) with an upper tolerance bound. Similarly, if L(Y) > −∞ and U(Y) = +∞, the interval [L(Y), +∞] is a one-sided TI (γ,P) with a lower tolerance bound only [16].
Now, suppose that the sample vector Y = (Y1, Y2, …, YN) is independently drawn from a standard normal distribution N(μ, σ2) with the population mean μ and the population variance σ2. We can use Y¯ and s to denote sample mean and sample variance, respectively. The mathematical expressions for Y¯ and s can be written as
(2)
(3)
The formulas to describe the one-sided TI can be expressed as
(4)
or
(5)
where k1 is the factor to ensure that at least a proportion P×100% of the entire population with a confidence interval level of γ×100% is satisfied based on the sample of size N. The value of k1 can be calculated based on the inverse cumulative distribution function for a non-central t distribution [17]:
(6)
where tN−1,γ,δ is the critical value (γ × 100 percentile), i.e., F(tN−1,γ,δ) = γ, for a non-central t distribution with N − 1 degrees-of-freedom and the non-centrality parameter δ=zPN. Here, zP is the critical value of the standard normal distribution associated with a cumulative probability of P×100%. For example, Fig. 1,2 shows that the accumulated probability for −∞ ≤ x ≤ 1 is 84%, which is equivalent to Pr[xzP=1]=84%. The value of zP can be found in Ref. [18] for any given P. Some commonly used P and the associated zP values are listed in Table 1.
If N is greater than 50 [19,20], k1 can be approximated using the following equation:
(7a)
(7b)
(7c)
Next, the upper and lower bounds for a two-sided TI are defined as
(8a)
and
(8b)
The exact solution for the two-sided k2 factor can be obtained by solving the following implicit nonlinear integral equation [21,22]:
(9)
where χN12 is the chi-squared probability distribution function with N − 1 degrees-of-freedom; χ1,P,δ2 is the P×100% percentile of the chi-squared distribution with one degree-of-freedom and non-centrality parameter δ = z2, 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 k2 factor [23]:
(10)
where χ1γ,N1,02 is the critical value (1γ percentile) of the central chi-squared distribution with N1 degrees-of-freedom and non-centrality parameter δ=0, i.e., F(χ1γ,N1,02)=1γ. Figure 2 shows an example chi-squared distribution with N1=3. The probability of having a value greater than χ10.57,3,02=2 is 57%, which is equivalent to Pr[xχ1γ,N1,02=2]=57%. The values of χ1γ,N1,02 for any given γ and N can be found in Ref. [18]. Some commonly used critical values of chi-square distribution are listed in Table 2. To evaluate the quality of the k2 approximation from Eq. (10), a comparison is conducted by solving the true solutions calculated from Eq. (9) using the R package developed in Ref. [22]. The relative error between approximation and true solution is defined as 1(k2(approximation))/(k2(true solution)). Figure 3 plots the number of sample N versus the relative error for P=γ=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=γ=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¯=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 χ1γ,N1,02=χ0.05,N1,02 and k2(N,γ,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., L9% and U11%. Otherwise, the plate will be rejected.

Two examples with EDS line-scan measurements of molybdenum concentration were used to demonstrate the procedure of qualification determination. Figure 4 shows the secondary electron image from the scanning electron microscope (SEM), together with its EDS line-scan in red for the first sample fuel plate (Plate 1). The domain size for Fig. 4 is 726 × 567 µm. Figure 5 shows the observed probability density function from a total of 999 molybdenum concentration measurements taken along the line. A bimodal distribution with two peaks is evident, the smaller of which represents the carbide regions. To exclude the anomalies from the carbides, only 946 observations are considered as effective measurements with molybdenum concentrations between 7 wt% and 13 wt%. Given the exclusion of carbides, molybdenum wt% is expected to be normally distributed and consistent with the assumptions described in Sec. 2. Next, the dependence of the molybdenum concentration measurements for Plate 1 is explored using the autocorrelation function (ACF), which is a commonly used method for detecting non-randomness in data. Given equally spaced measurements, the autocorrelation function r for lag j is defined as
(11)

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 k2 as a function of the number of samples, we randomly sampled Ns samples from the observed data vector, repeating the experiment 1000 times for each value of Ns in Table 3. The last column in this table shows the results for Ns = N = 946. For any given Ns, the k2 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 Ns and then remains constant for Ns greater than 20.

The relationship between the sample and population variance is expressed as [24]
(12)

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 k2s (multiplication of k2 and s), which decreases with increasing Ns. The last two rows list the corresponding lower and upper bounds for each sample size Ns. 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 Ns on plate qualification for Plate 2 is listed in Table 4. The minimum value of Ns 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.

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 k2(N, γ, P) is only a function of sample size N with fixed P = γ = 0.95. A plot of k2 versus N in Fig. 9 shows that k2 gradually decreases (lower curve) but its slope slowly increases (top curve) with the increasing N. The value of k2 eventually approaches an asymptote of ∼2. A gradient of k2 greater than −0.1% can be considered as a mathematical indication of stability. The corresponding value N at dk2/dN=0.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 k2 = 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:

  1. Collect the initial 133 EDS line-scan data with a normal 1 µm spacing between each measurement.

  2. Conduct the ACF analysis and then determine the critical lag distance.

  3. 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.

  4. 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 k2s ≤ 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).

References

1.
Wachs
,
D. M.
,
Clark
,
C. R.
, and
Dunavant
,
R. J.
,
2008
,
Conceptual Process Description for the Manufacture of Low-Enriched Uranium-Molybdenum Fuel
,
Idaho National Laboratory (INL)
,
Idaho Falls
.
2.
Van Den Berghe
,
S.
, and
Lemoine
,
P.
,
2014
, “
Review of 15 Years of High-Density Low-Enriched UMo Dispersion Fuel Development for Research Reactors in Europe
,”
Nucl. Eng. Technol.
,
46
(
2
), pp.
125
146
. 10.5516/NET.07.2014.703
3.
Meyer
,
M. K.
,
Hofman
,
G. L.
,
Hayes
,
S. L.
,
Clark
,
C. R.
,
Wiencek
,
T. C.
,
Snelgrove
,
J. L.
,
Strain
,
R. V.
, and
Kim
,
K. H.
,
2002
, “
Low-Temperature Irradiation Behavior of Uranium–Molybdenum Alloy Dispersion Fuel
,”
J. Nucl. Mater.
,
304
(
2–3
), pp.
221
236
. 10.1016/S0022-3115(02)00850-4
4.
Snelgrove
,
J. L.
,
Hofman
,
G. L.
,
Meyer
,
M. K.
,
Trybus
,
C. L.
, and
Wiencek
,
T. C.
,
1997
, “
Development of Very-High-Density Low-Enriched-Uranium Fuels 1 Work Supported by the US Department of Energy, Office of Nonproliferation and National Security, Under Contract No. W-31-109-ENG-38.1
,”
Nucl. Eng. Des.
,
178
(
1
), pp.
119
126
. 10.1016/S0029-5493(97)00217-3
5.
Hu
,
X.
,
Wang
,
X.
,
Joshi
,
V. V.
, and
Lavender
,
C. A.
,
2018
, “
The Effect of Thermomechanical Processing on Second Phase Particle Redistribution in U-10 wt% Mo
,”
J. Nucl. Mater.
,
500
, pp.
270
279
. 10.1016/j.jnucmat.2017.12.042
6.
Xu
,
Z.
,
Joshi
,
V.
,
Hu
,
S.
,
Paxton
,
D.
,
Lavender
,
C.
, and
Burkes
,
D.
,
2016
, “
Modeling the Homogenization Kinetics of as-Cast U-10 wt% Mo Alloys
,”
J. Nucl. Mater.
,
471
, pp.
154
164
. 10.1016/j.jnucmat.2015.11.026
7.
Garrett
,
C. E.
, and
Prasad
,
K.
,
2004
, “
The Art of Meeting Palladium Specifications in Active Pharmaceutical Ingredients Produced by Pd-Catalyzed Reactions
,”
Adv. Synth. Catal.
,
346
(
8
), pp.
889
900
. 10.1002/adsc.200404071
8.
Ren
,
H.
,
Hou
,
Z.
,
Huang
,
M.
,
Bao
,
J.
,
Sun
,
Y.
,
Tesfa
,
T.
, and
Ruby Leung
,
L.
,
2016
, “
Classification of Hydrological Parameter Sensitivity and Evaluation of Parameter Transferability Across 431 US MOPEX Basins
,”
J. Hydrol.
,
536
, pp.
92
108
. 10.1016/j.jhydrol.2016.02.042
9.
Ulewicz
,
R.
,
2003
, “
Quality Control System in Production of the Castings From Spheroid Cast Iron
,”
Metalurgija
,
42
(
1
), pp.
61
63
.
10.
Reid
,
R. D.
, and
Sanders
,
N. R.
,
2007
,
Operations Management: an Integrated Approach
,
John Wiley
,
New York
.
11.
Wang
,
X.
,
Xu
,
Z.
,
Soulami
,
A.
,
Hu
,
X.
,
Lavender
,
C.
, and
Joshi
,
V.
,
2017
, “
Modeling Early-Stage Processes of U-10 Wt.%Mo Alloy Using Integrated Computational Materials Engineering Concepts
,”
JOM
,
69
(
12
), pp.
2532
2537
. 10.1007/s11837-017-2608-z
12.
Fedorov
,
A.
,
Wells
,
W. M.
,
Kikinis
,
R.
,
Tempany
,
C. M.
, and
Vangel
,
M. G.
,
2014
, “
Application of Tolerance Limits to the Characterization of Image Registration Performance
,”
IEEE Trans. Med. Imaging
,
33
(
7
), pp.
1541
1550
. 10.1109/TMI.2014.2317796
13.
Rebafka
,
T.
,
Clémençon
,
S.
, and
Feinberg
,
M.
,
2007
, “
Bootstrap-Based Tolerance Intervals for Application to Method Validation
,”
Chemom. Intell. Lab. Syst.
,
89
(
2
), pp.
69
81
. 10.1016/j.chemolab.2007.06.001
14.
Sharma
,
G.
, and
Mathew
,
T.
,
2012
, “
One-Sided and Two-Sided Tolerance Intervals in General Mixed and Random Effects Models Using Small-Sample Asymptotics
,”
J. Am. Stat. Assoc.
,
107
(
497
), pp.
258
267
. 10.1080/01621459.2011.640592
15.
INL
,
2017
, “
Fuel Specification for MP-1, MP-2, and FSP-1
,” SPC-1691, Idaho National Laboratory, Idaho Falls.
16.
Guenther
,
W. C.
,
1972
, “
Tolerance Intervals for Univariate Distributions
,”
Nav. Res. Logist.
,
19
(
2
), pp.
309
333
. 10.1002/nav.3800190208
17.
2012
, NIST/SEMATECH e-Handbook of Statistical Methods.
18.
Loucks
,
D. P.
, and
van Beek
,
E.
,
2017
, “An Introduction to Probability, Statistics, and Uncertainty,”
Water Resource Systems Planning and Management
,
Springer
,
Cham
. pp.
213
300
.
19.
Lieberman
,
G. J.
,
1957
,
Tables for One-Sided Statistical Tolerance Limits
,
Applied Mathematics and Statistics Labs Stanford Univ
,
CA
.
20.
Natrella
,
M. G.
,
2013
,
Experimental Statistics
,
Courier Corporation
,
Mineola, NY
.
21.
Janiga
,
I.
,
Garaj
,
I.
, and
Witkovský
,
V.
,
2009
, “
On Exact Two-Sided Statistical Tolerance Intervals for Normal Distributions With Unknown Means and Unknown Common Variability
,”
Quality and Productivity Research Conference
, Yorktown Heights, New York, IBM Thomas J. Watson Research Center.
22.
Young
,
D. S.
,
2016
, “
Normal Tolerance Interval Procedures in the Tolerance Package
,”
R J.
,
8
(
2
), pp.
200
212
. 10.32614/RJ-2016-041
23.
Howe
,
W. G.
,
1969
, “
Two-Sided Tolerance Limits for Normal Populations, Some Improvements
,”
J. Am. Stat. Assoc.
,
64
(
326
), pp.
610
620
. 10.1080/01621459.1969.10500999
24.
Kenney
,
J. F.
, and
Keeping
,
E. S.
,
1951
,
Mathematics of Statistics, Part One
, 2nd ed.,
Van Nostrand
,
Princeton, NJ
.