Search NASA⌕ Search

SEARCH · Search NASA

Results for “randomized methods”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗
NASA NTRS1982

Advanced Multipurpose Rendezvous Tracking System Study

Rendezvous and docking (R&D) sensors needed to support Earth orbital operations of vehicles were investigated to determine the form they should take. An R&D sensor must enable an interceptor vehicle to determine both the relative position and the relative attitude of a target vehicle. Relative position determination is fairly straightforward and places few constraints on the sensor. Relative attitude determination, however, is more difficult. The attitude is calculated based on relative position measurements of several reflectors placed in a known arrangement on the target vehicle. The constraints imposed on the sensor by the attitude determination method are severe. Narrow beamwidth, wide field of view (fov), high range accuracy, and fast random scan capability are all required to determine attitude by this method. A consideration of these constraints as well as others imposed by expected operating conditions and the available technology led to the conclusion that the sensor should be a cw optical radar employing a semiconductor laser transmitter and an image dissector receiver.

Laurie, R. J.↗
NASA NTRS1985

Renormalization group formulation of large eddy simulation

Renormalization group (RNG) methods are applied to eliminate small scales and construct a subgrid scale (SSM) transport eddy model for transition phenomena. The RNG and SSM procedures are shown to provide a more accurate description of viscosity near the wall than does the Smagorinski approach and also generate farfield turbulence viscosity values which agree well with those of previous researchers. The elimination of small scales causes the simultaneous appearance of a random force and eddy viscosity. The RNG method permits taking these into account, along with other phenomena (such as rotation) for large-eddy simulations.

Yakhot, V.↗
DOE OSTI2018

Probability and Cumulative Density Function Methods for the Stochastic Advection-Reaction Equation

We present a cumulative density function (CDF) method for the probabilistic analysis of $d$-dimensional advection-dominated reactive transport in heterogeneous media. We employ a probabilistic approach in which epistemic uncertainty on the spatial heterogeneity of Darcy-scale transport coefficients is modeled in terms of random fields with given correlation structures. Our proposed CDF method employs a modified Large-Eddy-Diffusivity (LED) approach to close and localize the nonlocal equations governing the one-point PDF and CDF of the concentration field, resulting in a $(d + 1)$ dimensional PDE. Compared to the classsical LED localization, the proposed modified LED localization explicitly accounts for the mean-field advective dynamics over the phase space of the PDF and CDF. To illustrate the accuracy of the proposed closure, we apply our CDF method to one-dimensional single-species reactive transport with uncertain, heterogeneous advection velocities and reaction rates modeled as random fields.

advection-reaction equation↗
DOE OSTI2019

Radioablation by Image-Guided (HDR) Brachytherapy and Transarterial Chemoembolization in Hepatocellular Carcinoma: A Randomized Phase II Trial

Background and Aims: The aim of this single-center, open-label phase II study was to assess the efficacy of image-guided high-dose-rate (HDR) brachytherapy (iBT) compared with conventional transarterial embolization (cTACE) in unresectable hepatocellular carcinoma. Methods: Seventy-seven patients were treated after randomization to iBT or cTACE, as single or repeated interventions. Crossover was allowed if clinically indicated. The primary endpoint was time to untreatable progression (TTUP). Eligibility criteria included a Child–Pugh score of ≤ 8 points, absence of portal vein thrombosis (PVT) at the affected liver lobe, and ≤ 4 lesions. Survival was analyzed by using the Cox proportional hazard model with stratification for Barcelona Clinic Liver Cancer (BCLC) stages. Results: Twenty patients were classified as BCLC-A (iBT/cTACE 8/12), 35 as BCLC-B (16/19), and 22 as BCLC-C (13/9). The 1-, 2-, and 3-year TTUP probabilities for iBT compared with cTACE were 67.5% versus 55.2%, 56.0% versus 27.4%, and 29.5% versus 11.0%, respectively, with an adjusted hazard ratio (HR) of 0.49 (95% confidence interval 0.27–0.89; p = 0.019). The 1-, 2-, and 3-year TTPs for iBT versus cTACE were 56.0% versus 28.2%, 23.9% versus 6.3%, and 15.9% versus 6.3%, respectively, with an adjusted HR of 0.49 (0.29–0.85; p = 0.011). The 1-, 2-, and 3-year OS rates were 78.4% versus 67.7%, 62.0% versus 47.3%, and 36.7% versus 27.0%, respectively, with an adjusted HR of 0.62 (0.33–1.16; p = 0.136). Conclusions: This explorative phase II trial showed a superior outcome of iBT compared with cTACE in hepatocellular carcinoma and supports proceeding to a phase III trial.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2020

Data-driven predictive probability density function control of fiber length stochastic distribution shaping in refining process

Pulp is the most important raw material in paper in paper industries, whose Fiber Length Stochastic Distribution (FLSD) shaping directly determines the energy consumption and paper quality of the subsequent papermaking processes. However, the mean and variance are insufficient to describe the output FLSD shaping, which display non-Gaussian distribution properties. Therefore, the traditional control, method based on the mean and variance of fiber length is difficult to control the FLSD shaping effectively. In this paper, a novel data-driven predictive probability density function (PDF) control method is proposed for output FLSD shaping in refining process. Primarily, in order to improve the approximation accuracy of the PDF employing RBF neural network, the parameters (center value and width) of each RBF basis functions are tuned via utilizing iterative learning control (ILC) law, and the corresponding estimations of weights law can be obtained. Secondly, considering that the conventional linear model of weights vector has drawbacks of low accuracy and weak generalization ability, random vector functional link networks (RVFLNs) based data driven nonlinear modelling method is employed to characterize the prediction model between the input variables and weights vector. Finally, in order to reduce the randomness of the output FLSD, a minimum entropy control method under mean constraint is employed to design the predictive PDF controller for the output FLSD shaping with the help of the established nonlinear stochastic distribution model. Both simulations and experiments are given to illustrate the effectiveness and practicability of the proposed method.

Probability debsity function, Fiber length distrib↗
DOE OSTI2020

Comparing traditional and Bayesian approaches to ecological meta‐analysis

Abstract Despite the wide application of meta‐analysis in ecology, some of the traditional methods used for meta‐analysis may not perform well given the type of data characteristic of ecological meta‐analyses. We reviewed published meta‐analyses on the ecological impacts of global climate change, evaluating the number of replicates used in the primary studies ( n i ) and the number of studies or records ( k ) that were aggregated to calculate a mean effect size. We used the results of the review in a simulation experiment to assess the performance of conventional frequentist and Bayesian meta‐analysis methods for estimating a mean effect size and its uncertainty interval. Our literature review showed that n i and k were highly variable, distributions were right‐skewed and were generally small (median n i = 5, median k = 44). Our simulations show that the choice of method for calculating uncertainty intervals was critical for obtaining appropriate coverage (close to the nominal value of 0.95). When k was low (<40), 95% coverage was achieved by a confidence interval (CI) based on the t distribution that uses an adjusted standard error (the Hartung–Knapp–Sidik–Jonkman, HKSJ), or by a Bayesian credible interval, whereas bootstrap or z distribution CIs had lower coverage. Despite the importance of the method to calculate the uncertainty interval, 39% of the meta‐analyses reviewed did not report the method used, and of the 61% that did, 94% used a potentially problematic method, which may be a consequence of software defaults. In general, for a simple random‐effects meta‐analysis, the performance of the best frequentist and Bayesian methods was similar for the same combinations of factors ( k and mean replication), though the Bayesian approach had higher than nominal (>95%) coverage for the mean effect when k was very low ( k < 15). Our literature review suggests that many meta‐analyses that used z distribution or bootstrapping CIs may have overestimated the statistical significance of their results when the number of studies was low; more appropriate methods need to be adopted in ecological meta‐analyses.

Pappalardo, Paula↗
NASA NTRS1984

Advanced reliability method for fatigue analysis

When design factors are considered as random variables and the failure condition cannot be expressed by a closed form algebraic inequality, computations of risk (or probability of failure) may become extremely difficult or very inefficient. This study suggests using a simple and easily constructed second degree polynomial to approximate the complicated limit state in the neighborhood of the design point; a computer analysis relates the design variables at selected points. Then a fast probability integration technique (i.e., the Rackwitz-Fiessler algorithm) can be used to estimate risk. The capability of the proposed method is demonstrated in an example of a low cycle fatigue problem for which a computer analysis is required to perform local strain analysis to relate the design variables. A comparison of the performance of this method is made with a far more costly Monte Carlo solution. Agreement of the proposed method with Monte Carlo is considered to be good.

Wu, Y.-T.↗
NASA NTRS1978

A comparative study of Conroy and Monte Carlo methods applied to multiple quadratures and multiple scattering

An efficient numerical method of multiple quadratures, the Conroy method, is applied to the problem of computing multiple scattering contributions in the radiative transfer through realistic planetary atmospheres. A brief error analysis of the method is given and comparisons are drawn with the more familiar Monte Carlo method. Both methods are stochastic problem-solving models of a physical or mathematical process and utilize the sampling scheme for points distributed over a definite region. In the Monte Carlo scheme the sample points are distributed randomly over the integration region. In the Conroy method, the sample points are distributed systematically, such that the point distribution forms a unique, closed, symmetrical pattern which effectively fills the region of the multidimensional integration. The methods are illustrated by two simple examples: one, of multidimensional integration involving two independent variables, and the other, of computing the second order scattering contribution to the sky radiance.

Deepak, A.↗
NASA NTRS2003

Evaluation of Treadmill Exercise in a Lower Body Negative Pressure Chamber as a Countermeasure for Weightlessness-Induced Bone Loss: a Bed Rest Study with Identical Twins

Counteracting bone loss is required for future space exploration. We evaluated the ability of treadmill exercise in a LBNP chamber to counteract bone loss in a 30-day bed rest study. Eight pairs of identical twins were randomly assigned to sedentary control or exercise groups. Exercise within LBNP decreased the bone resorption caused by bed rest and may provide a countermeasure for spaceflight. INTRODUCTION: Bone loss is one of the greatest physiological challenges for extended-duration space missions. The ability of exercise to counteract weightlessness-induced bone loss has been studied extensively, but to date, it has proven ineffective. We evaluated the effectiveness of a combination of two countermeasures-treadmill exercise while inside a lower body negative pressure (LBNP) chamber-on bone loss during a 30-day bed rest study. MATERIALS AND METHODS: Eight pairs of identical twins were randomized into sedentary (SED) or exercise/LBNP (EX/LBNP) groups. Blood and urine samples were collected before, several times during, and after the 30-day bed rest period. These samples were analyzed for markers of bone and calcium metabolism. Repeated measures ANOVA was used to determine statistical significance. Because identical twins were used, both time and group were treated as repeated variables. RESULTS: Markers of bone resorption were increased during bed rest in samples from sedentary subjects, including the collagen cross-links and serum and urinary calcium concentrations. For N-telopeptide and deoxypyridinoline, there were significant (p < 0.05) interactions between group (SED versus EX/LBNP) and phase of the study (sample collection point). Pyridinium cross-links were increased above pre-bed rest levels in both groups, but the EX/LBNP group had a smaller increase than the SED group. Markers of bone formation were unchanged by bed rest in both groups. CONCLUSIONS: These data show that this weight-bearing exercise combined with LBNP ameliorates some of the negative effects of simulated weightlessness on bone metabolism. This protocol may pave the way to counteracting bone loss during spaceflight and may provide valuable information about normal and abnormal bone physiology here on Earth.

Twin Study↗
DOE OSTI2022

Parameter estimation for X-ray scattering analysis with Hamiltonian Markov Chain Monte Carlo

Bayesian-inference-based approaches, in particular the random-walk Markov Chain Monte Carlo (MCMC) method, have received much attention recently for X-ray scattering analysis. Hamiltonian MCMC, a state-of-the-art development in the field of MCMC, has become popular in recent years. It utilizes Hamiltonian dynamics for indirect but much more efficient drawings of the model parameters. We described the principle of the Hamiltonian MCMC for inversion problems in X-ray scattering analysis by estimating high-dimensional models for several motivating scenarios in small-angle X-ray scattering, reflectivity, and X-ray fluorescence holography. Hamiltonian MCMC with appropriate preconditioning can deliver superior performance over the random-walk MCMC, and thus can be used as an efficient tool for the statistical analysis of the parameter distributions, as well as model predictions and confidence analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Girsanov's transformation based variance reduced Monte Carlo simulation schemes for reliability estimation in nonlinear stochastic dynamics

The study considers the problem of simulation based time variant reliability analysis of nonlinear randomly excited dynamical systems. Attention is focused on importance sampling strategies based on the application of Girsanov's transformation method. Controls which minimize the distance function, as in the first order reliability method (FORM), are shown to minimize a bound on the sampling variance of the estimator for the probability of failure. Two schemes based on the application of calculus of variations for selecting control signals are proposed: the first obtains the control force as the solution of a two-point nonlinear boundary value problem, and, the second explores the application of the Volterra series in characterizing the controls. The relative merits of these schemes, vis-à-vis the method based on ideas from the FORM, are discussed. Illustrative examples, involving archetypal single degree of freedom (dof) nonlinear oscillators, and a multi-degree of freedom nonlinear dynamical system, are presented. The credentials of the proposed procedures are established by comparing the solutions with pertinent results from direct Monte Carlo simulations. - Highlights: • The distance minimizing control forces minimize a bound on the sampling variance. • Establishing Girsanov controls via solution of a two-point boundary value problem. • Girsanov controls via Volterra's series representation for the transfer functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Graphite Characterization: Baseline Variability Analysis Report

The Idaho National INL’s Baseline Material Property program provides many measurements of different material properties on a variety of nuclear graphite grades being investigated for irradiation-induced changes in material properties. The baseline material properties provide an unirradiated reference for comparison with measurements of those properties after irradiation. This study describes how uncertainty, or variance, in estimates of important material properties, like compressive strength, may depend on sample size, sample distribution and parameter estimation methods. Compressive strength data from two large collections of specimens from a billet of PCEA graphite (sample size = 230) and a billet of IG-110 graphite (sample size = 48) are used to illustrate these dependencies. The grades of graphite used represent approximate end members in terms of relative variability in strength measurements. Monte Carlo simulation and bootstrapping methods were used to examine variance in Weibull-parameter estimates as a function of sample size. Results were consistent with what previous statistical analyses have shown about the relationship between sample size and width of the confidence interval for maximum likelihood (ML) estimated parameters when samples are repeatedly drawn from a Weibull distribution. As a function of sample size, the precision of the modulus estimate increases faster than that of the characteristic strength. The modulus maximum likelihood estimates (MLEs) are biased for low sample sizes. The American Society for Testing and Materials document ASTM D7846-16 provides equations that describe the dependence of the variance of parameter MLEs on sample size. Using the same methodology, we show that the polynomials provided for variance dependence on sample size are inaccurate for sample sizes greater than ~120. We provide approximations to that relationship that illustrate relative constancy for sample size > ~10, in the logspace relationship of the modulus variance and the log transformed characteristic strength variance. Bootstrapping provided a means of comparing Weibull parameters obtained using ASME’s guidance for specimen collection to those obtained using different weighting in a case where strength is related to position in a billet and to grain orientation and a single distribution is used to characterize the billet. ASME specifies that samples be collected with equal representation in all locations and grain orientations. Parameter values obtained using ASME’s method were different from those obtained using random sampling of the complete data sets, and this was particularly true for the PCEA data. In addition, while the ASME approach yielded more conservative parameters for the PCEA data, the opposite was true for the IG-110 data. The difference observed for the PCEA data set likely reflects the larger number of against-grain specimens collected, which had generally higher strength than the with-grain specimens. This emphasizes the importance of the dependence of strength on location or grain orientation and suggests that it may be prudent to weight specimen collection toward the weaker subgroups if a conservative estimate as desired, as the equal weighting method recommended by ASME would will not always provide the most conservative Weibull parameters. Because spatial correlations commonly exist in certain grades of graphite, future work may involve developing a methodology for parameter inference that relaxes the assumption of independent specimens. ASTM and ASME codes recommend different parameter estimation methods and models of the Weibull distribution in different situations. Descriptions of the effects on parameter estimates for the example billets used in this study provide an example of how the method may affect the parameter estimate, and how—in this case—that compares to differences due to other factors, such as sample size and dependence on grain orientation or sample location.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Resolution of 100 photons and quantum generation of unbiased random numbers

Macroscopic quantum phenomena, such as observed in superfluids and superconductors, have led to promising technological advancements and some of the most important tests of fundamental physics. At present, quantum detection of light is mostly relegated to the microscale, where avalanche photodiodes are very sensitive to distinguishing single-photon events from vacuum but cannot differentiate between larger photon-number events. Beyond this, the ability to perform measurements to resolve photon numbers is highly desirable for a variety of quantum information applications including computation, sensing, and cryptography. True photon-number resolving detectors do exist, but they are currently limited to the ability to resolve on the order of 10 photons, which is too small for certain proposals. In this work, we extend photon measurement into the mesoscopic regime by implementing a detection scheme based on multiplexing highly quantum-efficient transition-edge sensors to accurately resolve photon numbers between zero and 100. Further, we then demonstrate the use of our system by implementing a quantum random number generator with no inherent bias. This method is based on sampling a coherent state in the photon-number basis and is robust against environmental noise, phase and amplitude fluctuations in the laser, loss and detector inefficiency as well as eavesdropping. Beyond true random number generation, our detection scheme serves as a means to implement quantum measurement and engineering techniques valuable for photonic quantum information processing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Multi-element flow-driven spectral chaos (ME-FSC) method for uncertainty quantification of dynamical systems

The flow-driven spectral chaos (FSC) is a recently developed method for tracking and quantifying uncertainties in the long-time response of stochastic dynamical systems using the spectral approach. The method uses a novel concept called enriched stochastic flow maps as a means to construct an evolving finite-dimensional random function space that is both accurate and computationally efficient in time. In this paper, we present a multi-element version of the FSC method (the ME-FSC method for short) to tackle (mainly) those dynamical systems that are inherently discontinuous over the probability space. In ME-FSC, the random domain is partitioned into several elements, and then the problem is solved separately on each random element using the FSC method. Subsequently, results are aggregated to compute the probability moments of interest using the law of total probability. To demonstrate the effectiveness of the ME-FSC method in dealing with discontinuities and long-time integration of stochastic dynamical systems, four representative numerical examples are presented in this paper, including the Van-der-Pol oscillator problem and the Kraichnan-Orszag three-mode problem. Results show that the ME-FSC method is capable of solving problems that have strong nonlinear dependencies over the probability space, both reliably and at low computational cost.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Scattering of electromagnetic waves from random media with strong permittivity fluctuations

By taking into account the singularity of the dyadic Green's function in the renormalization method, a theory is derived for vector electromagnetic wave propagation in a random medium with large permittivity fluctuations and with anisotropic correlation function. The strong fluctuation theory is then applied to a discrete scatterer problem in which the permittivity can assume only two values. The results are found to be consistent with those derived from discrete scatterer theory for all values of dielectric constants of the scatterers.

Tsang, L.↗
DOE OSTI2018

Metropolis–Hastings Monte Carlo Method for Neutron Emissivity Tomographic Inversion in Tokamak Plasma

The article presents a new approach to plasma neutron emissivity reconstruction based on the Metropolis–Hastings Monte Carlo algorithm. The algorithm is based on a biased random walk. A dedicated computer code generates pseudo-random samples within the domain that contains the solution. The properly chosen objective function ensures convergence to the desired solution. The method allows to overcome the problem of trapping in local minima by accepting, with some finite probability, changes in the direction opposite to that indicated by the objective function. The developed method was tested using a synthetic data set. The results show that the method provides reliable reconstructions and thus can be used as an alternative or complementary method to the commonly applied ones. The method capabilities were also tested against noise present in the line-integrated projections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Laser Confocal Microscopy Uncertainty Quantification Study

At Los Alamos National Laboratory (LANL), the Storage Safety and Engineering (SSE) team completes annual surveillance on a subset of in-use interim nuclear material storage containers in fulfilment of requirements outlined in DOE Manual M 441.1-1. The containers are selected through several methods, such as subject matter expert judgement, random selection, and trending items. Following these selections, the SSE team has the capacity to complete surveillance on 15-20 containers each fiscal year, composed of a combination of SAVY-4000 and Hagan storage containers. Through previous work, the stainless-steel components of the containers have been identified as life limiting components, with an emphasis on the thin-walled bodies. The team is focused on understanding the extent of general and pitting corrosion, due to observations of extensive corrosion from stored contents and bag-out-bag degradation. Quantifying corrosion effects on the thin-walled stainless steel container bodies, and understanding potential impacts to the respective design release rates and design qualification release rates is paramount to the team. To date, destructive examination (DE) has proven to be the most insightful method for developing an understanding on the extent of corrosion on used containers. To standardize this process, the SSE team developed a destructive examination guide for analyzing stainless steel components of the containers. Corroded containers of interest are identified during surveillance activities and set aside for sectioning and characterization. Following sectioning, a major step in the DE workflow is the utilization of laser confocal microscopy for scanning corroded samples of interest and extracting data on pits, such as count, depth, and equivalent diameter. Adhering to the techniques outlined in the DE guide, analysis has been completed on two Hagans and one SAVY-4000 container, with the maximum pit depth recorded as 139.1 ± 22.82 μm on a 17.5 year old Hagan. The findings from the completed destructive examinations will be utilized to support lifetime extension efforts of the SAVY-4000 as the team can better estimate corrosion rates and effects over time based on stored contents and age. Due to the implications of observing extreme pit depths that approach the nominal container body thickness of .0299 inches (0.759 mm) or minimum container thickness of 0.236” (0.6 mm), high confidence in the LCM measurements is desired. Through testing outlined in, it was concluded that the total error ascribed to the 20x objective when conducting large image mapping on the Keyence VK-X3050 laser confocal microscope (LCM) relative to a 50x objective (reference) is 16.4% (± 8.73%). For shallow features on the order of pristine SAVY surface defects (i.e. 5 μm), this uncertainty is appropriate. However, this conservative estimate of total error poses a fundamental concern for pit depths that approach the thickness of the measured samples. That is, with the measurement uncertainty currently employed on all measurements, the LCM would be unable to resolve if a pit with a depth of 515 μm is through wall. Standard step height samples were procured and used in the present study to assess the resolution and repeatability of height measurements. Understanding the resolution and repeatability of height measurements was the first focus of the team as it relates directly to pit depth, which is of primary concern. Calibration gratings were procured to evaluate the resolution and repeatability of measurements in the X and Y axes of the LCM stage. The results of the depth uncertainty study were conducted first and presented in the subsequent sections. The planar uncertainty study is appended to the depth study with conclusions from both summarized at the end of the report.

36 MATERIALS SCIENCE↗