Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical Testing Approach”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Estimation of Radiation Doses to U.S. Military Test Participants from Nuclear Testing: A Comparison of Historical Film-Badge Measurements, Dose Reconstruction and Retrospective Biodosimetry

Retrospective radiation dose estimations, whether based on physical or biological measurements, or on theoretical dose reconstruction, are limited in their precision and reliability, particularly for exposures that occurred many decades ago. Here, we studied living U.S. military test participants, believed to have received high-dose radiation exposures during nuclear testing-related activities approximately six decades ago, with two primary goals in mind. The first was to compare three different approaches of assessing past radiation exposures: 1. Historical personnel monitoring data alone; 2. Dose reconstruction based on varying levels of completeness of individual information, which can include film badge data; and 3. Retrospective biodosimetry using chromosome aberrations in peripheral blood lymphocytes. The second goal was to use the collected data to make the best possible estimates of bone marrow dose received by a group with the highest military recorded radiation doses of any currently living military test participants. Six nuclear test participants studied had been on Rongerik Atoll during the 1954 CASTLE Bravo nuclear test. Another six were present at the Nevada Test Site (NTS) and/or Pacific Proving Ground (PPG) and were believed to have received relatively high-dose exposures at those locations. All were interviewed, and all provided a blood sample for cytogenetic analysis. Military dose records for each test participant, as recorded in the Defense Threat Reduction Agency's Nuclear Test Review and Information System, were used as the basis for historical film badge records and provided exposure scenario information to estimate dose via dose reconstruction. Dose to bone marrow was also estimated utilizing directional genomic hybridization (dGH) for high-resolution detection of radiation-induced chromosomal translocations and inversions, the latter being demonstrated for the first time for the purpose of retrospective biodosimetry. As the true dose for each test participant is not known these many decades after exposure, this study gauged the congruence of different methods by assessing the degree of correlation and degree of systematic differences. Overall, the best agreement between methods, defined by statistically significant correlations and small systematic differences, was between doses estimated by a dose reconstruction methodology that exploited all the available individual detail and the biodosimetry methodology derived from a weighted average dose determined from chromosomal translocation and inversion rates. Employing such a strategy, we found that the Rongerik veterans who participated in this study appear to have received, on average, bone marrow equivalent doses on the order of 300–400 mSv, while the NTS/ PPG participants appear to have received approximately 250–300 mSv. The results show that even for nuclear events that occurred six decades in the past, biological signatures of exposure are still present, and when taken together, chromosomal translocations and inversions can serve as reliable retrospective biodosimeters, particularly on a group-average basis, when doses received are greater than statistically-determined detection limits for the biological assays used.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Dark energy survey year 3 results: Cosmology with peaks using an emulator approach

ABSTRACT We constrain the matter density Ωm and the amplitude of density fluctuations σ8 within the ΛCDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps created by taking a harmonic space product of the convergence of two tomographic redshift bins. Our analysis follows a forward-modelling scheme to create a likelihood of these statistics using N-body simulations, using a Gaussian process emulator. We take into account the uncertainty from the remaining, largely unconstrained ΛCDM parameters (Ωb, ns, and h). We include the following lensing systematics: multiplicative shear bias, photometric redshift uncertainty, and galaxy intrinsic alignment. Stringent scale cuts are applied to avoid biases from unmodelled baryonic physics. We find that the additional non-Gaussian information leads to a tightening of the constraints on the structure growth parameter yielding $S_8~\equiv ~\sigma _8\sqrt{\Omega _{\mathrm{m}}/0.3}~=~0.797_{-0.013}^{+0.015}$ (68 per cent confidence limits), with a precision of 1.8 per cent, an improvement of 38 per cent compared to the angular power spectra only case. The results obtained with the angular power spectra and peak counts are found to be in agreement with each other and no significant difference in S8 is recorded. We find a mild tension of $1.5 \, \sigma$ between our study and the results from Planck 2018, with our analysis yielding a lower S8. Furthermore, we observe that the combination of angular power spectra and tomographic peak counts breaks the degeneracy between galaxy intrinsic alignment AIA and S8, improving cosmological constraints. We run a suite of tests concluding that our results are robust and consistent with the results from other studies using DES Y3 data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Greedy emulators for nuclear two-body scattering

Applications of reduced basis method emulators are increasing in low-energy nuclear physics because they enable fast and accurate sampling of high-fidelity calculations, enabling robust uncertainty quantification. Here, in this paper, we develop, implement, and test two model-driven emulators based on the (Petrov-)Galerkin projection using the prototypical test case of two-body scattering with the Minnesota potential and a more realistic local chiral potential. The high-fidelity scattering equations are solved with the matrix Numerov method, a reformulation of the popular Numerov recurrence relation for solving special second-order differential equations as a linear system of coupled equations. A novel error estimator based on reduced-space residuals is applied to an active learning approach (a greedy algorithm) to choosing training samples (“snapshots”) for the emulator and contrasted with a proper orthogonal decomposition (POD) approach. Both approaches allow for computationally efficient offline-online decompositions, but the greedy approach requires many fewer snapshot calculations. These developments set the groundwork for emulating scattering observables based on chiral nucleon-nucleon and three-nucleon interactions and optical models, where computational speed-ups are necessary for Bayesian uncertainty quantification. Our emulators and error estimators are widely applicable to linear systems.

Bayesian methods↗
DOE OSTI2016

A mathematical and numerical framework for the analysis of compressible thermal convection in gases at very high temperatures

The relevance of non-equilibrium phenomena, nonlinear behavior, gravitational effects and fluid compressibility in a wide range of problems related to high-temperature gas-dynamics, especially in thermal, mechanical and nuclear engineering, calls for a concerted approach using the tools of the kinetic theory of gases, statistical physics, quantum mechanics, thermodynamics and mathematical modeling in synergy with advanced numerical strategies for the solution of the Navier–Stokes equations. The reason behind such a need is that in many instances of relevance in this field one witnesses a departure from canonical models and the resulting inadequacy of standard CFD approaches, especially those traditionally used to deal with thermal (buoyancy) convection problems. Starting from microscopic considerations and typical concepts of molecular dynamics, passing through the Boltzmann equation and its known solutions, we show how it is possible to remove past assumptions and elaborate an algorithm capable of targeting the broadest range of applications. Moving beyond the Boussinesq approximation, the Sutherland law and the principle of energy equipartition, the resulting method allows most of the fluid properties (density, viscosity, thermal conductivity, heat capacity and diffusivity, etc.) to be derived in a rational and natural way while keeping empirical contamination to the minimum. Special attention is deserved as well to the well-known pressure issue. With the application of the socalled multiple pressure variables concept and a projection-like numerical approach, difficulties with such a term in the momentum equation are circumvented by allowing the hydrodynamic pressure to decouple from its thermodynamic counterpart. The final result is a flexible and modular framework that on the one hand is able to account for all the molecule (translational, rotational and vibrational) degrees of freedom and their effective excitation, and on the other hand can guarantee adequate interplay between molecular and macroscopic-level entities and processes. Performances are demonstrated by computing some incompressible and compressible benchmark test cases for thermal (gravitational) convection, which are then extended to the high-temperature regime taking advantage of the newly developed features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Raman Spectroscopy for the On Line Analysis of Oxidation States of Oxygen Carrier Particles

The need for cleaner energy systems, including CO 2 capture technologies, is driving the current development of new combustion approaches such as chemical looping combustion (CLC). This process is based on using oxygen carrier particles (OCPs) such as iron oxides, copper oxides, calcium sulfates, etc. to provide the oxygen for the combustion process. In order to optimize the overall process performance, it is critical that the properties of the oxygen carriers are well-defined and maintained for their specific purpose during the different stages of the CLC process. One of the critical properties of the oxygen carrier is its oxidation state (e.g., content of Fe 2 O 3 vs. Fe 3 O 4 ) as it affects the fundamental operation of the CLC process. Unfortunately, the ability to make on-line measurements of the oxidation state of oxygen carriers is lacking and new sensors need to be developed. Our goal is to develop a sensor for the on-line analysis of the oxidation state of oxygen carrier particles and to demonstrate its feasibility. Oxygen carrier particles (OCPs) of interest include iron oxides, copper oxides, and calcium sulfate, with particle sizes typically less than 0.5 mm diameter. The expected operating conditions of the sensor include OCP temperatures in the range of 800 °C to 1000 °C, and pressures of about 10 atm. The expected results are to be of a statistical nature, providing relative concentrations of the different oxidation states. These requirements make a strong case for using non-contact, stand-off, pulsed Raman spectroscopy as a sensor for the high-temperature on-line analysis of oxygen carrier particles. In order to achieve our goal, we have set up and tested a pulsed Raman spectroscopy system in combination with a pressurized high-temperature sample chamber, and have taking initial Raman spectra of various materials at high temperature. We are now in the process of optimizing the operating parameters such as laser wavelength, laser intensity, collection optic design, focal spot size, data acquisition, etc. We will also develop an analysis procedure for the interpretation of the Raman spectra. Typical Raman spectra consist of a convolution of background signal and multiple individual spectra. In order to extract relevant information such as the concentration ratio of a material in different oxidation states, we need to process (background subtraction, deconvolution, fitting, etc.) the measured data. In addition, statistical modeling and multivariate calibration need to be performed to correctly analyze the Raman spectra.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2017

An optimal FFT-based anisotropic power spectrum estimator

Measurements of line-of-sight dependent clustering via the galaxy power spectrum's multipole moments constitute a powerful tool for testing theoretical models in large-scale structure. Recent work shows that this measurement, including a moving line-of-sight, can be accelerated using Fast Fourier Transforms (FFTs) by decomposing the Legendre polynomials into products of Cartesian vectors. Here, we present a faster, optimal means of using FFTs for this measurement. We avoid redundancy present in the Cartesian decomposition by using a spherical harmonic decomposition of the Legendre polynomials. With this method, a given multipole of order ℓ requires only 2ℓ+1 FFTs rather than the (ℓ+1)(ℓ+2)/2 FFTs of the Cartesian approach. For the hexadecapole (ℓ = 4), this translates to 40% fewer FFTs, with increased savings for higher ℓ. The reduction in wall-clock time enables the calculation of finely-binned wedges in P ( k ,μ), obtained by computing multipoles up to a large ℓ{sub max} and combining them. This transformation has a number of advantages. We demonstrate that by using non-uniform bins in μ, we can isolate plane-of-sky (angular) systematics to a narrow bin at 0μ ≅ while eliminating the contamination from all other bins. We also show that the covariance matrix of clustering wedges binned uniformly in μ becomes ill-conditioned when combining multipoles up to large values of ℓ{sub max}, but that the problem can be avoided with non-uniform binning. As an example, we present results using ℓ{sub max}=16, for which our procedure requires a factor of 3.4 fewer FFTs than the Cartesian method, while removing the first μ bin leads only to a 7% increase in statistical error on f σ{sub 8}, as compared to a 54% increase with ℓ{sub max}=4.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2004

A Study of Penalty Function Methods for Constraint Handling with Genetic Algorithm

COMETBOARDS (Comparative Evaluation Testbed of Optimization and Analysis Routines for Design of Structures) is a design optimization test bed that can evaluate the performance of several different optimization algorithms. A few of these optimization algorithms are the sequence of unconstrained minimization techniques (SUMT), sequential linear programming (SLP) and the sequential quadratic programming techniques (SQP). A genetic algorithm (GA) is a search technique that is based on the principles of natural selection or "survival of the fittest". Instead of using gradient information, the GA uses the objective function directly in the search. The GA searches the solution space by maintaining a population of potential solutions. Then, using evolving operations such as recombination, mutation and selection, the GA creates successive generations of solutions that will evolve and take on the positive characteristics of their parents and thus gradually approach optimal or near-optimal solutions. By using the objective function directly in the search, genetic algorithms can be effectively applied in non-convex, highly nonlinear, complex problems. The genetic algorithm is not guaranteed to find the global optimum, but it is less likely to get trapped at a local optimum than traditional gradient-based search methods when the objective function is not smooth and generally well behaved. The purpose of this research is to assist in the integration of genetic algorithm (GA) into COMETBOARDS. COMETBOARDS cast the design of structures as a constrained nonlinear optimization problem. One method used to solve constrained optimization problem with a GA to convert the constrained optimization problem into an unconstrained optimization problem by developing a penalty function that penalizes infeasible solutions. There have been several suggested penalty function in the literature each with there own strengths and weaknesses. A statistical analysis of some suggested penalty functions is performed in this study. Also, a response surface approach to robust design is used to develop a new penalty function approach. This new penalty function approach is then compared with the other existing penalty functions.

Ortiz, Francisco↗
NASA NTRS2022

Charging Experiments and Filtration Testing on Facemask Materials: Towards Goals of Full Performance Restoration of Used N95 Masks and Enhancement of Common-Fabric Facemasks

Different techniques have been developed for sanitizing N95 masks while retaining whatever charge remains on them. However, these sanitizing methods do not address the fact that the electrical charge, filtration efficacy, and breathability of such masks may be expected to degrade during long-term use based on mechanisms known in the literature. In this study, approaches to electrically (re)charging facemask filter materials were explored in response to the COVID-19 induced shortage of N95 masks. The effort consisted of charging experiments followed by filtration testing conducted at the coupon level during a month-long period (May 2020). The main focus of the work was to develop a recharging technique which could eventually be adapted and used as part of a broader full-performance restoration process for used N95 masks, e.g., which might employ simple washing by alcohol or soapy water as a preceding step for sanitization and declogging. Most tested coupons were therefore nonwoven meltblown fibrous polypropylene (PP) filter media extracted from decades old N95 masks. Keeping in mind commonly used cloth facemasks, limited trials were also performed with some commercially available fabrics. Charging was done at elevated temperature and was followed by cooling to freeze in the applied charges produced from corona discharges using different electrode configurations and parameter combinations, or from triboelectric rubbing. Charging was also focused on using potentially scalable equipment, e.g., Van de Graaff generators or consumer ionic hair products. Some degree of success was achieved with one of the fabrics and the PP filter media which was first purposefully discharged to have typical filtration efficiency (FE) of ~35%. Numerous attempts at corona charging discharged media resulted in filtration efficiency (FE) ≥ 80 – 85%. Some measure of charge robustness was also established via stable performance during loading testing. However, results where FE ≥ 90% were not reliably repeatable and no results met the FE ≥ 95% goal due to shortcomings of the charging process. The best three trials for PP media statistically met the FE ≥ 90% criteria, amounted to a small percentage of the trials relying on an ionic hair dryer, and are attributed to good embedding of charges via aerodynamic force. Triboelectric charging significantly enhanced the viability of one of the polyester fabrics –measured FE increased from 19% up to 46% while maintaining a reasonably low pressure drop.

N95↗
NASA NTRS2010

Monitoring System for Storm Readiness and Recovery of Test Facilities: Integrated System Health Management (ISHM) Approach

Severe weather events are likely occurrences on the Mississippi Gulf Coast. It is important to rapidly diagnose and mitigate the effects of storms on Stennis Space Center's rocket engine test complex to avoid delays to critical test article programs, reduce costs, and maintain safety. An Integrated Systems Health Management (ISHM) approach and technologies are employed to integrate environmental (weather) monitoring, structural modeling, and the suite of available facility instrumentation to provide information for readiness before storms, rapid initial damage assessment to guide mitigation planning, and then support on-going assurance as repairs are effected and finally support recertification. The system is denominated Katrina Storm Monitoring System (KStorMS). Integrated Systems Health Management (ISHM) describes a comprehensive set of capabilities that provide insight into the behavior the health of a system. Knowing the status of a system allows decision makers to effectively plan and execute their mission. For example, early insight into component degradation and impending failures provides more time to develop work around strategies and more effectively plan for maintenance. Failures of system elements generally occur over time. Information extracted from sensor data, combined with system-wide knowledge bases and methods for information extraction and fusion, inference, and decision making, can be used to detect incipient failures. If failures do occur, it is critical to detect and isolate them, and suggest an appropriate course of action. ISHM enables determining the condition (health) of every element in a complex system-of-systems or SoS (detect anomalies, diagnose causes, predict future anomalies), and provide data, information, and knowledge (DIaK) to control systems for safe and effective operation. ISHM capability is achieved by using a wide range of technologies that enable anomaly detection, diagnostics, prognostics, and advise for control: (1) anomaly detection algorithms and strategies, (2) fusion of DIaK for anomaly detection (model-based, numerical, statistical, empirical, expert-based, qualitative, etc.), (3) diagnostics/prognostics strategies and methods, (4) user interface, (5) advanced control strategies, (6) integration architectures/frameworks, (7) embedding of intelligence. Many of these technologies are mature, and they are being used in the KStorMS. The paper will describe the design, implementation, and operation of the KStorMS; and discuss further evolution to support other needs such as condition-based maintenance (CBM).

Figueroa, Fernando↗
NASA NTRS2014

Reply to Comment by Laprise on 'the Added Value to Global Model Projections of Climate Change by Dynamical Downscaling: a Case Study over the Continental U.S. Using the GISS-ModelE2 and WRF Models'

In his comment, Laprise raises several points that we agree merit consideration. His primary critique is that our study [Racherla et al., 2012] tested the ability of the WRF regional climate model to reproduce historical temperature and precipitation change relative to the driving global climate model (GCM) using only a single simulation rather than an ensemble. He asserts that the observed changes are smaller than the internal variability in the climate system (i.e., not statistically significant) and that thus a single simulation should not necessarily be able to capture the observations. Laprise points out that the statistical signal is reduced for a multi-decadal trend such as the one we analyzed in comparison with mean climatology and cites two studies showing that for particular climate parameters it can take any years for a signal to be discerned over internal variability. He states that The results of theexperiment as designed were strongly influenced by the presence of internal variability and sampling errors,which masked the rather small climate changes that may have occurred as a consequence of changes inforcing during the period considered. While Laprise discusses statistics in general terms at some length, for the actual climate trends examined in our study, he offers no evidence that the forced signal was smallcompared with internal variability. The two studies he cites [de Ela et al., 2013; Maraun, 2013] do not provide convincing evidence as they concern climate variables averaged over different times and areas. One in fact examines extreme precipitation events, which by definition are rare and thus have a lower significance level. We accept the general point that it is important to consider internal variability, and as noted in our paper we agree that an ensemble of simulations is in principle an optimal, though computationally expensive, approach. While we did not present the statistical significance of the observations in our original paper, we have now evaluated those for the regional temperature trends used in our study to evaluate the added value of WRF and thus can analyze data as to the magnitude of the trends with respect to internal variability.

climate↗
NASA NTRS1993

Testing higher-order Lagrangian perturbation theory against numerical simulation. 1: Pancake models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasi-linear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the Zel'dovich approximation (hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of sigma is approximately 2). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-spectrum with power-index n = -1) using cross-correlation statistics employed in previous work. We found that for all statistical methods used the higher-order corrections improve the results obtained for the first-order solution up to the stage when sigma (linear theory) is approximately 1. While this improvement can be seen for all spatial scales, later stages retain this feature only above a certain scale which is increasing with time. However, third-order is not much improvement over second-order at any stage. The total breakdown of the perturbation approach is observed at the stage, where sigma (linear theory) is approximately 2, which corresponds to the onset of hierarchical clustering. This success is found at a considerable higher non-linearity than is usual for perturbation theory. Whether a truncation of the initial power-spectrum in hierarchical models retains this improvement will be analyzed in a forthcoming work.

Buchert, T.↗
NASA NTRS2003

Single-Vector Calibration of Wind-Tunnel Force Balances

An improved method of calibrating a wind-tunnel force balance involves the use of a unique load application system integrated with formal experimental design methodology. The Single-Vector Force Balance Calibration System (SVS) overcomes the productivity and accuracy limitations of prior calibration methods. A force balance is a complex structural spring element instrumented with strain gauges for measuring three orthogonal components of aerodynamic force (normal, axial, and side force) and three orthogonal components of aerodynamic torque (rolling, pitching, and yawing moments). Force balances remain as the state-of-the-art instrument that provide these measurements on a scale model of an aircraft during wind tunnel testing. Ideally, each electrical channel of the balance would respond only to its respective component of load, and it would have no response to other components of load. This is not entirely possible even though balance designs are optimized to minimize these undesirable interaction effects. Ultimately, a calibration experiment is performed to obtain the necessary data to generate a mathematical model and determine the force measurement accuracy. In order to set the independent variables of applied load for the calibration 24 NASA Tech Briefs, October 2003 experiment, a high-precision mechanical system is required. Manual deadweight systems have been in use at Langley Research Center (LaRC) since the 1940s. These simple methodologies produce high confidence results, but the process is mechanically complex and labor-intensive, requiring three to four weeks to complete. Over the past decade, automated balance calibration systems have been developed. In general, these systems were designed to automate the tedious manual calibration process resulting in an even more complex system which deteriorates load application quality. The current calibration approach relies on a one-factor-at-a-time (OFAT) methodology, where each independent variable is incremented individually throughout its full-scale range, while all other variables are held at a constant magnitude. This OFAT approach has been widely accepted because of its inherent simplicity and intuitive appeal to the balance engineer. LaRC has been conducting research in a "modern design of experiments" (MDOE) approach to force balance calibration. Formal experimental design techniques provide an integrated view to the entire calibration process covering all three major aspects of an experiment; the design of the experiment, the execution of the experiment, and the statistical analyses of the data. In order to overcome the weaknesses in the available mechanical systems and to apply formal experimental techniques, a new mechanical system was required. The SVS enables the complete calibration of a six-component force balance with a series of single force vectors.

Parker, P. A.↗
DOE OSTI2022

Searching for Clues for a Matter Dominated Universe in Liquid Argon Time Projection Chambers

Liquid Argon Time Projection Chambers (LArTPCs) represent one of the most widely utilized neutrino detection techniques in neutrino experiments, for instance, in the Short Baseline Neutrino (SBN) program and the future large-scale LArTPC: Deep Underground Neutrino Experiment (DUNE). The high-end technique, facilitating excellent spatial and calorimetric reconstruction resolution, also enables testing exotic Beyond Standard Model (BSM) theories, such as baryon number violation (BNV) processes (e.g., proton-decay, neutron-antineutron oscillation). At the same time, Machine Learning (ML) techniques have demonstrated their ubiquitous use in recent decades; ML techniques have also become some of the most powerful tools in high-energy physics (HEP) analyses. Furthermore, the development of algorithms to cater to the needs of problems in HEP (i.e., triggering, reconstruction, improving sensitivity, etc.) has also become an active area of research. By developing a combined approach using Convolutional Neural Network (CNN) and Boosted Decision Tree (BDT) techniques, the sensitivity of neutron-antineutron oscillation in DUNE is evaluated for a projected exposure of 400kton·years. Additionally, to meet the triggering requirement to select such rare events in DUNE, such a search is only supported with highly efficient self-triggering algorithms. An ML-based self-triggering scheme for large-scale LArTPCs, such as DUNE, is also developed with the intention of implementation on field-programmable gate arrays (FPGAs). The ML-based approach for searching for neutron-antineutron oscillation can be demonstrated and validated on the current LArTPC MicroBooNE. The analysis in MicroBooNE represents the first-ever search for neutron-antineutron oscillation in a LArTPC. DUNE's projected 90% C.L. sensitivity to the neutron antineutron oscillation lifetime is 6.45×10³² years, assuming 1.327×10³⁵ neutron·years, equivalent to 10 years of DUNE far detector exposure (400kton·years). For MicroBooNE, assuming 372 seconds of exposure (equivalent to 3.13×10³⁶ neutron·years), the 90% C.L. lifetime sensitivity is found at 3.07×10²⁵ yrs, after accounting for Monte-Carlo statistical uncertainty and systematic uncertainty from detector effects.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2019

Audibility of Multiple, Low-Frequency Tonal Signals in Noise

The main purpose of this study is to examine the audibility of multiple, low-frequency tones that are placed in distinct auditory channels. Three experiments are described, the goals of which are to determine if the presence of sound in multiple channels results in enhanced audibility and to assess the applicability of the Statistical Summation Model (SSM) to this frequency range. This model predicts that for the case of multiple signals that are in separate auditory channels, implying statistical independence, each with sensitivity value d prime of i, the resulting total sensitivity is given by the square root of the sum of the squares of the individual d prime of i values. In common with previous studies conducted at higher frequencies, the signals are pure tones and the maskers are broadband noise. The requirement that low frequency tones be placed in separate auditory filters limited the number of tones to a maximum of three. The first of the three experiments measured the change in masked thresholds for two- and three-tone signals relative to the level of the equally-detectable single tones. The multiple tone signals were composed of combinations of 55, 120 and 200 Hz tones. The measured changes in thresholds exceeded those predicted by the SSM, although they did not differ statistically from the model predictions. The second experiment employed the same overall approach but acquired more data and concentrated on the three-tone signal. Once again, the measured changes in masked threshold exceeded the model predictions, this time to a statistically-significant degree. Two issues were postulated with the potential to yield inflated changes in masked threshold: interaction between tones resulting in perceptible intermodulation/difference tones, and the assumption that the tones were in distinct auditory filters and statistically independent of one another. The third experiment used two sets of three-tone signals to address these latter concerns. The first set of three tones was composed of harmonically related tone frequencies of 55, 110 and 165 Hz, which was an attempt to reduce effects of intermodulation difference tones. The second set of three tones was chosen to be 110, 220 and 330 Hz, again reducing effects of difference tones, but also providing greater separation between tones. Results for the first set of three tones compared to those of the earlier experiments indicated that intermodulation was not an important effect. The second set of three tones (110, 220, 330 Hz) yielded changes in masked thresholds that, on average, were in good agreement with the SSM, although intersubject variability was large and prohibited a definitive conclusion regarding the concern that tone spacing was inadequate. The results of the three experiments showed that the masked threshold of sounds with multiple (two or three) equally-detectable low frequency tones was lower than those of the single tones. In other words, it is clear that audibility is enhanced by the presence of signals in multiple auditory filters. This finding is consistent with most previous research conducted at higher frequencies. In contrast with previous research, test subjects were, on average, able to detect multitone sounds at lower levels than those predicted using the SSM. Analyses that included Monte Carlo simulations showed that normally distributed errors in the single tone thresholds result in biased estimates of the thresholds of multitone sounds. This phenomenon is likely responsible for at least a substantial fraction of the unexpected deviation of measurements from SSM predictions.

Rafaelof, Menachem↗
DOE OSTI2025

Part distortion monitoring in additive manufacturing using machining

In additive manufacturing, accumulation of residual stresses can result in severe part distortion from the desired preform shape. Current methods for in-situ part distortion monitoring in additive manufacturing typically require expensive sensors, or capital equipment, and require time-consuming post-processing to understand the shape deviation. This paper presents an in-situ method, in the context of hybrid manufacturing, for part distortion detection using machining of additively manufactured parts. As a surrogate, three test artifacts were used to represent different distorted geometries. The tool axis positions from the machine tool controller and the cutting power were monitored during a facing operation. Cutting power data was used to detect the tool entry and exit in the workpiece using a novel approach with power standard deviation metric. The workpiece geometry and distorted configuration was subsequently predicted for positional and rotational deviations to within 2 mm accuracy using synchronized tool position data with cutting power. The proposed method can be used in a hybrid (additive and subtractive) machine tool to periodically check part distortion in the additive build. The method is applicable for any additive process and is low-cost and computationally inexpensive.

36 MATERIALS SCIENCE↗
NASA NTRS2024

Digital Lunar Exploration Sites (DLES) Terrain Crafting

Humans will soon be returning to the surface of the Moon with NASA’s Artemis program. The Artemis program is an international collaboration that will consist of a complex series of space systems and missions to explore the lunar surface and pave the way for the future exploration of Mars. NASA and its partners rely heavily on simulation for lighting and navigation studies as well as training astronauts, flight controllers, and mission support staff. The NASA Exploration Systems Simulations (NExSyS) team in the Simulation and Graphics Branch (ER7) in the Engineering Directorate at NASA’s Johnson Space Center has built up many simulation products to support this effort, one of which is the Digital Lunar Exploration Sites (DLES). DLES is a collection of products used to simulate and render the lunar surface in a digital environment. We discussed and presented an overview of the DLES products at the 2022 IEEE Aerospace Conference in Big Sky, MT with a paper titled "Digital Lunar Exploration Sites". This “DLES Terrain Crafting” paper will expand on the information previously provided in “DLES” paper and dive deeper into the details of the terrain crafting process and the toolsets used to support this task. The best digital data currently available of the lunar surface is provided by the Lunar Reconnaissance Orbiter (LRO). Its Lunar Orbiter Laser Altimeter (LOLA) achieves an impressive resolution of 5m per pixel at the Lunar South Pole (LSP) and can generate datasets covering a large continuous region near the LSP. There are a few additional methods, such as Shape from Shading which can infer higher resolution data (up to 1m per pixel) from the LRO Narrow Angle Camera (NAC) images. However, surface-based simulations require higher-resolution data, and this paper will discuss the process of enhancing the terrain to meet that need. The process begins with capturing statistical data of craters in the regions of interest using images provided by the LRO NAC. This data is then used to scatter artificial features which are not captured in the truth data, resulting in an enhanced DEM with a much higher resolution of 20cm per pixel. Many tools were built up to assist in the creation of these artificial Digital Elevation Models (DEM), which this paper will discuss in detail. DEMs themselves are a very powerful representation of a planetary surface, and many operations and tools can utilize the data they contain. This paper includes a description of the rendering of the lunar surface in a graphics engine, generation of contact patches to simulate tire to ground interaction, and ray tracing utilities to model Line of Sight (LOS) interactions with the terrain. This paper will also explore some new tool sets currently under development which aim to utilize Machine Learning (ML) to assist in the identification of craters from LRO NAC imagery. While this is not a novel idea, the NExSyS team is developing a unique approach which may result in more robust identification of crater characteristics.

Artemis↗
DOE OSTI2018

Characterizing Natural Attenuation of Cr(VI) in a Chromium Plume at Los Alamos National Laboratory - 18594

Characterization of natural attenuation of hexavalent (Cr(VI)) on aquifer sediments within the Cr plume at Los Alamos National Laboratory has been conducted to (1) refine conceptual and numerical models of chromium fate and transport, (2) evaluate the long-term stability of Cr(III)/Cr(VI) in the complex groundwater system along flow paths from the source to the regional aquifer, and (3) evaluate and develop options for successful remediation. The naturally high chromium abundance in aquifer sediments at Los Alamos pose an unique challenge for characterization of chromium reduction, making conventional approaches such as uptake tests and advanced spectroscopic measurements ineffective in identifying relatively small amounts of anthropogenic Cr(VI) in the sediments. To address this challenge and as complimentary to these conventional methods, sequential leaching tests were conducted on core hole sediments obtained from sonic drilled cores taken from both contaminated and non-contaminated horizons among two mineralogically distinct formations (Puye formation and Miocene Pumiceous Unit) within the Cr plume. A series of three and coherent sequential leaching steps were conducted: 1) an alkaline leach to determine the amount of Cr(VI) sorbed on the sediments. 2) a weak acid leach (using 0.01 M and 0.1 M HNO{sub 3}) to determine the amount of anthropogenic chromium (Cr(III) and Cr(VI)) on aquifer sediment surfaces, and 3) a strong acid leach (using 2 M HNO{sub 3}) to characterize the amount of releasable Cr(III) from aquifer minerals. Additionally we tested the efficiency and detection limit of sequential leaching on bentonite, which consists mainly of smectite, used in drill muds and laboratory-grade biotite spiked with known concentrations of Cr(VI). Smectite and biotite are two phases considered to be potentially important host minerals in aquifer sediments within the plume. Thus, by comparing contaminated and uncontaminated sediment leachate results of chromium and other major and trace elements, the signature and quantity of anthropogenic Cr(VI) can be derived. Our results from the weak acid leaching test indicate that the method detection limit for Cr 'uptake' was about 1 μg/g for biotite and about 0.5 μg/g for bentonite. Our results also indicate that all samples taken from known contaminated horizons in comparison to non-contaminated ones within the Puye formation showed no statistical difference in the amount of Cr leached indicative of its insignificant Cr(VI) attenuation capacity. The amount of excess Cr leached from the magnetic fraction of the contaminated Miocene Pumiceous samples was about 2 μg/g of sediment. Considering that the magnetic fraction was about 11% by mass of the bulk sample, the excess amount of Cr leached from the bulk material would be about 0.2 μg/g of sediment, suggesting the potential for some Cr(VI) attenuation capacity in the Miocene pumice unit. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

nautilus : boosting Bayesian importance nested sampling with deep learning

ABSTRACT We introduce a novel approach to boost the efficiency of the importance nested sampling (INS) technique for Bayesian posterior and evidence estimation using deep learning. Unlike rejection-based sampling methods such as vanilla nested sampling (NS) or Markov chain Monte Carlo (MCMC) algorithms, importance sampling techniques can use all likelihood evaluations for posterior and evidence estimation. However, for efficient importance sampling, one needs proposal distributions that closely mimic the posterior distributions. We show how to combine INS with deep learning via neural network regression to accomplish this task. We also introduce nautilus, a reference open-source python implementation of this technique for Bayesian posterior and evidence estimation. We compare nautilus against popular NS and MCMC packages, including emcee, dynesty, ultranest, and pocomc, on a variety of challenging synthetic problems and real-world applications in exoplanet detection, galaxy SED fitting and cosmology. In all applications, the sampling efficiency of nautilus is substantially higher than that of all other samplers, often by more than an order of magnitude. Simultaneously, nautilus delivers highly accurate results and needs fewer likelihood evaluations than all other samplers tested. We also show that nautilus has good scaling with the dimensionality of the likelihood and is easily parallelizable to many CPUs.

97 MATHEMATICS AND COMPUTING↗