Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods:statistical”

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 973 records · Page 54

A physical-based gas–surface interaction model for rarefied gas flow simulation

Empirical gas–surface interaction models, such as the Maxwell model and the Cercignani–Lampis model, are widely used as the boundary condition in rarefied gas flow simulations. The accuracy of these models in the prediction of macroscopic behavior of rarefied gas flows is less satisfactory in some cases especially the highly non-equilibrium ones. Molecular dynamics simulation can accurately resolve the gas–surface interaction process at atomic scale, and hence can predict accurate macroscopic behavior. They are however too computationally expensive to be applied in real problems. In this work, a statistical physical-based gas–surface interaction model, which complies with the basic relations of boundary condition, is developed based on the framework of the washboard model. In virtue of its physical basis, this new model is capable of capturing some important relations/trends for which the classic empirical models fail to model correctly. As such, the new model is much more accurate than the classic models, and in the meantime is more efficient than MD simulations. Therefore, it can serve as a more accurate and efficient boundary condition for rarefied gas flow simulations. - Highlights: • The new physical-based gas–surface interaction model satisfies the basic relations of boundary condition so that it can be used in rarefied gas flow simulations. • The new physical model possesses several advantages comparing to the empirical gas–surface interaction models, hence it's a more accurate boundary condition for rarefied gas flow study. • The new physical model is a statistical model which can be efficiently implemented by Monte Carlo method. It's much more efficient comparing to the MD simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Recent advances, and unresolved issues, in the application of computational modelling to the prediction of the biological effects of nanomaterials

Nanomaterials research is one of the fastest growing contemporary research areas. The unprecedented properties of these materials have meant that they are being incorporated into products very quickly. Regulatory agencies are concerned they cannot assess the potential hazards of these materials adequately, as data on the biological properties of nanomaterials are still relatively limited and expensive to acquire. Computational modelling methods have much to offer in helping understand the mechanisms by which toxicity may occur, and in predicting the likelihood of adverse biological impacts of materials not yet tested experimentally. This paper reviews the progress these methods, particularly those QSAR-based, have made in understanding and predicting potentially adverse biological effects of nanomaterials, and also the limitations and pitfalls of these methods. - Highlights: • Nanomaterials regulators need good information to make good decisions. • Nanomaterials and their interactions with biology are very complex. • Computational methods use existing data to predict properties of new nanomaterials. • Statistical, data driven modelling methods have been successfully applied to this task. • Much more must be learnt before robust toolkits will be widely usable by regulators.

60 APPLIED LIFE SCIENCES↗

A method of surface marker location optimization for tumor motion estimation in lung stereotactic body radiation therapy

Purpose: Accurately localizing lung tumor localization is essential for high-precision radiation therapy techniques such as stereotactic body radiation therapy (SBRT). Since direct monitoring of tumor motion is not always achievable due to the limitation of imaging modalities for treatment guidance, placement of fiducial markers on the patient’s body surface to act as a surrogate for tumor position prediction is a practical alternative for tracking lung tumor motion during SBRT treatments. In this work, the authors propose an innovative and robust model to solve the multimarker position optimization problem. The model is able to overcome the major drawbacks of the sparse optimization approach (SOA) model. Methods: The principle-component-analysis (PCA) method was employed as the framework to build the authors’ statistical prediction model. The method can be divided into two stages. The first stage is to build the surrogate tumor matrix and calculate its eigenvalues and associated eigenvectors. The second stage is to determine the “best represented” columns of the eigenvector matrix obtained from stage one and subsequently acquire the optimal marker positions as well as numbers. Using 4-dimensional CT (4DCT) and breath hold CT imaging data, the PCA method was compared to the SOA method with respect to calculation time, average prediction accuracy, prediction stability, noise resistance, marker position consistency, and marker distribution. Results: The PCA and SOA methods which were both tested were on all 11 patients for a total of 130 cases including 4DCT and breath-hold CT scenarios. The maximum calculation time for the PCA method was less than 1 s with 64 752 surface points, whereas the average calculation time for the SOA method was over 12 min with 400 surface points. Overall, the tumor center position prediction errors were comparable between the two methods, and all were less than 1.5 mm. However, for the extreme scenarios (breath hold), the prediction errors for the PCA method were not only smaller, but were also more stable than for the SOA method. Results obtained by imposing a series of random noises to the surrogates indicated that the PCA method was much more noise resistant than the SOA method. The marker position consistency tests using various combinations of 4DCT phases to construct the surrogates suggested that the marker position predictions of the PCA method were more consistent than those of the SOA method, in spite of surrogate construction. Marker distribution tests indicated that greater than 80% of the calculated marker positions fell into the high cross correlation and high motion magnitude regions for both of the algorithms. Conclusions: The PCA model is an accurate, efficient, robust, and practical model for solving the multimarker position optimization problem to predict lung tumor motion during SBRT treatments. Due to its generality, PCA model can also be applied to other imaging guidance system whichever using surface motion as the surrogates.

60 APPLIED LIFE SCIENCES↗

Scaled Vecchia Approximation for Fast Computer-Model Emulation

Many scientific phenomena are studied using computer experiments consisting of multiple runs of a computer model while varying the input settings. Gaussian processes (GPs) are a popular tool for the analysis of computer experiments, enabling interpolation between input settings, but direct GP inference is computationally infeasible for large datasets. We adapt and extend a powerful class of GP methods from spatial statistics to enable the scalable analysis and emulation of large computer experiments. Specifically, we apply Vecchia’s ordered conditional approximation in a transformed input space, with each input scaled according to how strongly it relates to the computer-model response. The scaling is learned from the data by estimating parameters in the GP covariance function using Fisher scoring. Our methods are highly scalable, enabling estimation, joint prediction, and simulation in near-linear time in the number of model runs. In several numerical examples, our approach substantially outperformed existing methods.

97 MATHEMATICS AND COMPUTING↗

Expanded method for the determination of burnup in nuclear fuels using multiple neodymium isotopes

Established methods for the determination of burnup in nuclear fuels commonly rely on the measurement of 148Nd in the spent fuel in concert with known fission product yields to determine the atom percent of fissions in the fuel. This isotope of neodymium is used for various reasons, including chemical and radioactive stability, ease of measurement, and low rates of formation and destruction due to neutron flux apart from fission. However, careful calculation of effective cumulative fission yields and correction factors for (n,γ) capture reactions allows for additional stable and long-lived isotopes of neodymium to be used to provide additional independent measurements of burnup, reducing statistical uncertainty. This method was developed and successfully applied to measure the burnup of compacts from the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program. The mean burnup measured using 143Nd, 145Nd, 146Nd, 148Nd, and 150Nd was statistically observed to be the same as that measured using 148Nd alone, but the statistical uncertainty in the measurement was reduced by a factor of 2, providing a tighter confidence interval in the final results.

Helmreich, Grant [ORNL] (ORCID:0000000330464394)↗

Bayesian inference on the sphere beyond statistical isotropy

We present a general method for Bayesian inference of the underlying covariance structure of random fields on a sphere. We employ the Bipolar Spherical Harmonic (BipoSH) representation of general covariance structure on the sphere. We illustrate the efficacy of the method as a principled approach to assess violation of statistical isotropy (SI) in the sky maps of Cosmic Microwave Background (CMB) fluctuations. SI violation in observed CMB maps arise due to known physical effects such as Doppler boost and weak lensing; yet unknown theoretical possibilities like cosmic topology and subtle violations of the cosmological principle, as well as, expected observational artefacts of scanning the sky with a non-circular beam, masking, foreground residuals, anisotropic noise, etc. We explicitly demonstrate the recovery of the input SI violation signals with their full statistics in simulated CMB maps. Our formalism easily adapts to exploring parametric physical models with non-SI covariance, as we illustrate for the inference of the parameters of a Doppler boosted sky map. Our approach promises to provide a robust quantitative evaluation of the evidence for SI violation related anomalies in the CMB sky by estimating the BipoSH spectra along with their complete posterior.

79 ASTRONOMY AND ASTROPHYSICS↗

End-Use Savings Shapes Measure Documentation: Dispatch Schedule Generation for Demand Flexibility Measures

This supplemental document describes the methodology used for determining the dispatch timing of various EUSS demand flexibility measures. Demand flexibility measures are designed to reduce/dispatch electricity demand in buildings during especially beneficial/critical times. The method used in this work utilizes predictions of building loads to generate a schedule that reflects the periods when the building's daily peak load occurs to support decision making in demand flexibility measures. The dispatch schedule generation method described in this document creates an hourly schedule that includes a load dispatch (peak) window for each day for a whole year based on load prediction, with options using different prediction methods: perfect prediction, bin-sampling method, fixed schedule, and outdoor air temperature (OAT)-based prediction method. The perfect prediction method performs a simulation to obtain the annual load profile as predicted load, representing the scenario of perfect load prediction. The bin-sampling method (1) categorizes days into representative bins by temperature characteristics, (2) performs simulations on sample days from each of those bins to create representative (or predicted) load, and (3) assigns representative loads for all days in a year based on the bin categorization. The fixed schedule method defines uniform start and end time of peak window with assumed fixed daily peak time, for all days in a season or a year. The OAT-based prediction method uses the statistics of OAT (minimum and maximum) as the indicators of peak load, with specified delay response time from building loads to temperature. Given the load prediction, daily peak periods are determined as a time window with specified length in each day that include the predicted daily peak load and with a secondary rule such as maximizing energy saving potential. The dispatch schedule generation method is not a standalone measure and is intended to be combined with other demand flexibility measures that could leverage the peak schedule and apply demand controls on specific systems or devices for demand response, such as measures described in "Measure Documentation - Thermostat Control for Load Shedding" and "Measure Documentation - Thermostat Control for Load Shifting".

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Structural fluctuation of proteins induced by thermodynamic perturbation

A theory to describe structural fluctuations of protein induced by thermodynamic perturbations, pressure, temperature, and denaturant, is proposed. The theory is formulated based on the three methods in the statistical mechanics: the generalized Langevin theory, the linear response theory, and the three dimensional interaction site model (3D-RISM) theory. The theory clarifies how the change in thermodynamic conditions, or a macroscopic perturbation, induces the conformational fluctuation, which is a microscopic property. The theoretical results are applied, on the conceptual basis, to explain the experimental finding by Akasaka et al., concerning the NMR experiment which states that the conformational change induced by pressure corresponds to structural fluctuations occurring in the ambient condition. A method to evaluate the structural fluctuation induced by pressure is also suggested by means of the 3D-RISM and the site-site Kirkwood-Buff theories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design, Detection, and Countermeasure of Frequency Spectrum Attack and Its Impact on Long Short-Term Memory Load Forecasting and Microgrid Energy Management

This paper introduces a frequency-domain false data injection attack called Frequency Spectrum Attack (FSA) and explores its effects on load forecasting and the energy management system (EMS) in a microgrid. The FSA analyzes time-series signals in the frequency domain to identify patterns in their frequency spectrum. It learns the distribution of dominant frequencies in a dataset of healthy signals. Subsequently, it manipulates the amplitudes of dominant frequencies within this healthy distribution, ensuring a stealthy attack against statistical analysis of the signal spectrum. We evaluated the performance of FSA on LSTM, a state-of-the-art network for load forecasting. The results show that FSA can triple the Mean Absolute Error (MAE) of predictions compared to the normal case and increase it by 70% compared to noise injection attacks. Furthermore, FSA indirectly enhances battery utilization in the EMS by 45%. We then proposed a detection method that combines statistical analysis and machine-learning-based classification techniques with features. The model effectively distinguishes FSA from healthy and noisy signals, achieving an accuracy of 98.7% and an F1-score of 98.1% on a load dataset, covering healthy, FSA, and noisy load data. Finally, a countermeasure was introduced based on the statistical analysis of the frequency spectrum of healthy signals to mitigate the impact of FSA. This countermeasure successfully reduces the MAE of the attacked model from 0.135 to 0.053, validating its effectiveness in mitigating FSA.

Nazeri, Amirhossein↗

Measurement of the muon anomalous precession frequency in runs 2 \& 3 of the Muon $g-2$ Experiment at Fermilab

This dissertation presents a measurement of the muon anomalous precession frequency for the Runs 2 \& 3 data of the E989 Muon $g-2$ Experiment at Fermilab. The muon anomalous precession frequency is one of two key inputs, the other being the magnetic field, used to precisely determine the muon magnetic anomaly, $a_{\mu}$. In April 2021, the E989 collaboration reported its first measurement of $a_{\mu}$ to an unprecedented precision of 460 parts-per-billion (ppb). This result is in agreement with the previous measurement performed by the E821 collaboration at Brookhaven National Laboratory, and the combined experimental value is in tension with the Standard Model prediction at 4.2$\sigma$, a possible hint of new physics. The first result from E989 was based on the Run-1 data, which was collected in 2018 and comprises 6\% of the experiment’s target statistics; the Run-1 result was statistics-limited. The Runs 2 \& 3 data were collected in 2019-2020 and constitute a four-fold increase i n statistics compared to Run-1 and consequently, a factor of two reduction in the measurement’s statistical uncertainty. This reduced statistical uncertainty, as well as continued understanding of systematic effects, will result in an even more precise measurement of $a_{\mu}$ and will help clarify the observed tension between theory and experiment. The precession frequency analysis presented in this dissertation includes a number of improvements compared to Boston University's Run-1 analysis: the implementation of a more robust pileup-subtraction procedure, the implementation of a kernel ratio method, the adoption of the statistically optimal asymmetry-weighted method, and investigations that led to an improved understanding of an outstanding residual slow effect. This dissertation motivates a measurement of the muon magnetic anomaly, describes the experimental principle, gives an overview of the E989 experiment, and presents a precession frequency analysis with full systematic unc ertainty evaluation for the Runs 2 \& 3 data. The combined precession frequency measurement, using the ratio asymmetry-weighted method, has a statistical uncertainty of 201 ppb and a systematic uncertainty of 25 ppb, constituting the most precise determination of the muon anomalous precision frequency to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Testing cosmic geometry without dynamic distortions using voids

We propose a novel technique to probe the expansion history of the Universe based on the clustering statistics of cosmic voids. In particular, we compute their two-point statistics in redshift space on the basis of realistic mock galaxy catalogs and apply the Alcock-Paczynski test. In contrast to galaxies, we find void auto-correlations to be marginally affected by peculiar motions, providing a model-independent measure of cosmological parameters without systematics from redshift-space distortions. Because only galaxy-galaxy and void-galaxy correlations have been considered in these types of studies before, the presented method improves both statistical and systematic uncertainties on the product of angular diameter distance and Hubble rate, furnishing the potentially cleanest probe of cosmic geometry available to date.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimating region-specific fuel economy in the United States from real-world driving cycles

Here we describe a method for estimating region-specific real-world light-duty vehicle fuel economy in the United States that is unique in both the size and representativeness of real-world driving that was considered, and for its ability to model regional variations in driving patterns. Over one million miles of national driving data were used to select real-world cycles representative of observed trip categories. The six cycles were compared to U.S. legislative cycles, revealing some key differences. Finally, a set of cycle weighting factors for 533 separate U.S. regions was derived from annual traffic statistics. Applying this method, it was found that regional fuel economy varies due to differences in driving patterns alone and that rural driving patterns lead to improved fuel economy (for conventional vehicles). The driving cycles and regional weighting factors described here are useful for testing and simulation studies, specifically those sensitive to regional variations in driving patterns.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Study of the filter method for neutron pulse-height distributions measured with organic scintillators

Identification of neutron sources is of interest in several fields, such as homeland security or basic nu-clear physics. Organic scintillators offer the possibility to unfold fast-neutron energy spectra. Since organic scintillators are sensitive to neutrons and gamma-rays, a method to discriminate these two types of particles is required. This paper is focused on investigation of the discrimination properties of the filter method by using a new procedure. We compare the results of an experimental filter method to a digital pulse-shape discrimination (PSD) method based on charge integration. In addition, these methods are compared to the simulation results obtained with the MCNPX-PoliMi code. The experimental and nu-merical investigations were performed with a 252Cf spontaneous-fission neutron source. The integration of the number of counts gave the relative differences between the experimental filter and digital PSD neutron pulse-height distributions (PHDs) and simulated PHDs less than approximately 5% in the range between 60 keVee and 1.715 MeVee of light output. Above 1.715 MeVee, the PSD method has advantages over the filter method, due to the filter method having significantly worse counting statistics. The results show that the filter method has potential for robust neutron measurements when the PSD method cannot be applied, such as for ‘old’ organic plastic scintillators without PSD capability.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗

Opportunities in AI/ML for the Rubin LSST Dark Energy Science Collaboration

The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will produce unprecedented volumes of heterogeneous astronomical data (images, catalogs, and alerts) that challenge traditional analysis pipelines. The LSST Dark Energy Science Collaboration (DESC) aims to derive robust constraints on dark energy and dark matter from these data, requiring methods that are statistically powerful, scalable, and operationally reliable. Artificial intelligence and machine learning (AI/ML) are already embedded across DESC science workflows, from photometric redshifts and transient classification to weak lensing inference and cosmological simulations. Yet their utility for precision cosmology hinges on trustworthy uncertainty quantification, robustness to covariate shift and model misspecification, and reproducible integration within scientific pipelines. This white paper surveys the current landscape of AI/ML across DESC's primary cosmological probes and cross-cutting analyses, revealing that the same core methodologies and fundamental challenges recur across disparate science cases. Since progress on these cross-cutting challenges would benefit multiple probes simultaneously, we identify key methodological research priorities, including Bayesian inference at scale, physics-informed methods, validation frameworks, and active learning for discovery. With an eye on emerging techniques, we also explore the potential of the latest foundation model methodologies and LLM-driven agentic AI systems to reshape DESC workflows, provided their deployment is coupled with rigorous evaluation and governance. Finally, we discuss critical software, computing, data infrastructure, and human capital requirements for the successful deployment of these new methodologies, and consider associated risks and opportunities for broader coordination with external actors.

Aubourg, Eric [APC, Paris] (ORCID:000000025592023X↗

Uncertainty analysis of Raman spectra for measuring ortho-parahydrogen compositions

Determining ortho-parahydrogen compositions at cryogenic temperatures is an important quality control for liquid hydrogen custody exchange. High orthohydrogen compositions lead to an exothermic reaction resulting in increased boil off and increased venting losses of liquid hydrogen product by either the supplier or consumer. Traditional methods for measuring ortho-parahydrogen compositions such as hot-wire anemometry, nuclear magnetic resonance, and infrared spectroscopy are typically inconvenient at best. As the number of cryogenic hydrogen systems continue to increase, there is an increased need for flexible, standardized techniques for post processing composition measurements. A Raman spectrometer implemented in the Cryo-Catalysis Hydrogen Experimental Facility (CHEF) has demonstrated incredible flexibility for composition measurements of cryogenic hydrogen flows through a wide range of temperatures and pressures with the Raman probe in-situ at cryogenic temperatures. This article analyzes the uncertainty of Raman spectroscopy data of cryogenic hydrogen equilibrated at temperatures in IONEX catalyst. Spectra from hydrogen catalyzed to the equilibrium ortho-parahydrogen composition are processed and compared to expected statistical distributions for method verification.

08 HYDROGEN↗

First 3D vector tracking of helium recoils for fast neutron measurements at Super$\mathrm{KEKB}$

Here we present results from the first deployment of novel, high definition, compact gas Time Projection Chambers (TPCs) with pixel chip readout as part of the BEAST II beam background measurement project at SuperKEKB. The TPCs provide detailed 3D imaging of ionization from neutron-induced nuclear recoils in a helium and carbon dioxide target gas mixture at standard temperature and pressure. We present the TPC performance and the neutron backgrounds observed during the initial stage of collider commissioning. We find excellent electron background rejection, leading to background-free nuclear recoil measurements above 50 keV ee , despite the extreme high-background environment. We measure an angular resolution better than 20° for recoil tracks longer than 1.7 mm, corresponding to an average ionization energy of approximately 100 keV ee . We also obtain the full 3D vector direction of helium recoils by utilizing charge profile measurements along the recoil axis, with a correct head/tail assignment efficiency of approximately 80%. With this performance, we present comparisons between measured and simulated event rates, recoil energy spectra, and directional distributions originating from beam–gas and Touschek beam losses at SuperKEKB. We utilize head/tail recognition to distinguish neutron components traveling with positive radial velocity in the Belle II coordinate system from those traveling in the opposite direction. Finally, we present a novel method of discriminating beam–gas interactions from Touschek beam losses that can eliminate the need for dedicated accelerator runs for background measurements. This method is still statistics-limited. However, future studies should be able to verify this method, which in turn could lead to neutron background analysis runs symbiotic with normal Belle II operation. The capabilities demonstrated here also suggest that high definition recoil imaging in gas TPCs is applicable to low energy, low-background experiments, such as directional dark matter searches.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Estimating Code Biases for Criticality Safety Applications with Few Relevant Benchmarks

Criticality safety analyses rely on the availability of relevant benchmark experiments to determine justifiable margins of subcriticality. When a target application lacks neutronically similar benchmark experiments, validation studies must provide justification to the regulator that the impact of modeling and simulation limitations is well understood for the application and often must provide additional subcritical margin to ensure safe operating conditions. This study estimated the computational bias in the critical eigenvalue for several criticality safety applications supported by only a few relevant benchmark experiments. The accuracy of the following three methods for predicting computational biases was evaluated: the Upper Subcritical Limit STATisticS (USLSTATS) trending analysis method; the Whisper nonparametric method; and TSURFER, which is based on the generalized linear least-squares technique. Here, these methods were also applied to estimate computational biases and recommended upper subcriticality limits for several critical experiments with known biases and for several cases from a blind benchmark study. The methods are evaluated based on both the accuracy of their predicted computation bias and upper subcriticality limit estimates, as well as on the consistency of the methods’ estimates, as the model parameters, covariance data libraries, and set of available benchmark data were varied. Data assimilation methods typically have not been used for criticality safety licensing activities, and this study explores a methodology to address concerns regarding the reliability of such methods in criticality safety bias prediction applications.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗