Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
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 OSTI2015

Boltzmann-conserving classical dynamics in quantum time-correlation functions: “Matsubara dynamics”

We show that a single change in the derivation of the linearized semiclassical-initial value representation (LSC-IVR or “classical Wigner approximation”) results in a classical dynamics which conserves the quantum Boltzmann distribution. We rederive the (standard) LSC-IVR approach by writing the (exact) quantum time-correlation function in terms of the normal modes of a free ring-polymer (i.e., a discrete imaginary-time Feynman path), taking the limit that the number of polymer beads N → ∞, such that the lowest normal-mode frequencies take their “Matsubara” values. The change we propose is to truncate the quantum Liouvillian, not explicitly in powers of ħ{sup 2} at ħ{sup 0} (which gives back the standard LSC-IVR approximation), but in the normal-mode derivatives corresponding to the lowest Matsubara frequencies. The resulting “Matsubara” dynamics is inherently classical (since all terms O(ħ{sup 2}) disappear from the Matsubara Liouvillian in the limit N → ∞) and conserves the quantum Boltzmann distribution because the Matsubara Hamiltonian is symmetric with respect to imaginary-time translation. Numerical tests show that the Matsubara approximation to the quantum time-correlation function converges with respect to the number of modes and gives better agreement than LSC-IVR with the exact quantum result. Matsubara dynamics is too computationally expensive to be applied to complex systems, but its further approximation may lead to practical methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
DOE OSTI2023

Evaluating the Feasibility of Using Strain Measured During Sinusoidal Rate Pumping Tests to Characterize an Aquifer

Pumping tests with sinusoidal variation in pumping rate have been proposed as a method for improving aquifer characterization. These tests can interrogate a larger aquifer volume than slug tests and they can be more sensitive to small variations in drawdown. Current methods of using sinusoidal variations of rate are based on measuring pressure signals from the reservoir or aquifer, which requires access to monitoring wells. An alternative approach has been developed that measures the strain in the vadose zone instead of pressure in the reservoir. An instrument has been developed at Clemson University that can measure small strains using optical fiber sensors, and we propose that it can be used to measure strain during sinusoidal pumping tests. The objective of this project is to evaluate the feasibility of measuring and interpreting the vertical strain in the vadose zone generated during a sinusoidal rate pumping tests to estimate aquifer and/or overburden properties. Pumping tests were performed at a site with a water table at 9 m and saturated saprolite from 9 to 30m depth. Sinusoidal pumping tests in the saprolite used different pumping amplitudes and periods while pressure was measured at three piezometers and vertical strain was measured in the vadose zone from 1 to 7 m depth. The delay time of the pressure increased linearly with distance, and it increased with the square root of the period, while the amplitude of the pressure decreased with distance and increased with amplitude. Vertical strain in the vadose zone varied with the same period as the pumping rate and the delay times were similar to those for pressures measured at similar radial distances. Delay times for the strain increased with an increasing propagation of the strain field. The hydraulic diffusivity was estimated using delay times for pressure and strain, and the results were compared with the values estimated from a constant-rate pumping test. The hydraulic diffusivity from the constant rate test was 0.039 m2/sec, whereas it was 0.053 (±0.05) m2/sec and 0.035 (±0.03) m2/sec from the pressure data and 0.08 (±0.022) m2/sec and 0.035 (±0.001) m2/sec from the strain data from sinusoidal tests of various periods and amplitudes. There is no statistical difference among these values. This study indicates that strain measurements in the vadose zone can be used to estimate aquifer properties during sinusoidal pumping tests.

Blais, Riley↗
DOE OSTI2022

Narrow-Band Least-Squares Infrasound Array Processing

Infrasound data from arrays can be used to detect, locate, and quantify a variety of natural and anthropogenic sources from local to remote distances. However, many array processing methods use a single broad frequency range to process the data, which can lead to signals of interest being missed due to the choice of frequency limits or simultaneous clutter sources. In this work, we introduce a new open-source Python code that processes infrasound array data in multiple sequential narrow frequency bands using the least-squares approach. We test our algorithm on a few examples of natural sources (volcanic eruptions, mass movements, and bolides) for a variety of array configurations. Our method reduces the need to choose frequency limits for processing, which may result in missed signals, and it is parallelized to decrease the computational burden. Improvements of our narrow-band least-squares algorithm over broad-band least-squares processing include the ability to distinguish between multiple simultaneous sources if distinct in their frequency content (e.g., microbarom or surf vs. volcanic eruption), the ability to track changes in frequency content of a signal through time, and a decreased need to fine-tune frequency limits for processing. We incorporate a measure of planarity of the wavefield across the array (sigma tau, στ) as well as the ability to utilize the robust least trimmed squares algorithm to improve signal processing and insight into array performance. Our implementation allows for more detailed characterization of infrasound signals recorded at arrays that can improve monitoring and enhance research capabilities.

58 GEOSCIENCES↗
DOE OSTI2016

High Accuracy Signal Validation Framework for Sensor Calibration Assessment in NPPs

Effective and efficient health management of sensing channels and instrumentation supports reliable and accurate process monitoring capacities in any industry. Current practice in the nuclear industry relies on periodic calibration assessment and correction activities. This time based maintenance technique consists of three steps: i) removal of sensor to a test environment, ii) testing and recalibration, and iii) reinstallation of sensor into its working environment. This procedure is highly labor intensive, adds to the outage cost and timeline, and contributes to worker radiation exposure and ALARA requirements. Furthermore, the human element contributes to error factors in all stages of the procedure. The effectiveness of this calibration procedure after sensors are returned to working environmental variables also raises concerns, and the practicability of this periodic approach is uncertain for future reactor designs with longer operating cycles between refueling outages. Studies indicate that periodic calibration assessment find only 10% of sensing instrumentations with calibration issues during a single outage (18- 24 months), meaning at least 90% of sensors are undergoing unnecessary maintenance actions. This provides an opportunity to maintain the current levels of safety and functional integrity with fewer, more targeted maintenance actions and lower economic burden. Performance-based maintenance planning shows much promise over current periodic maintenance approaches. Specifically, online monitoring (OLM) utilizes the continually collected process data to non-intrusively monitor the calibration status of process instrumentation and detect the onset of de-calibration, thereby minimizing the number of recalibration procedures required. Ultimately, this can lead to termination of time based recalibration at NPPs. OLM utilizes the process measurements made for process monitoring and control to evaluate the calibration status of the sensors and instrumentation making the measurements. OLM can detect and differentiate sensor faults and process anomalies, thus providing a valuable tool to suit the industry's and regulatory body's requirements. Though supported by more than a decade of research, OLM for sensor calibration assessment is yet to be implemented in the US fleet due to regulatory requirements and constraints. As shown in Figure 1, OLM utilizes process data to train statistical models that predict normal, error-free process behavior. Measured process values are compared with error-corrected predictions to ascertain the presence of faults. The reliability of OLM is dependent on the quality of its predictions, measured in part by the uncertainty in these predictions. This research develops a framework to accurately quantify the prediction uncertainty based on Bayesian inference techniques to provide the required confidence levels for OLM systems. The expected outcomes include signal validation by high-accuracy, high-confidence uncertainty quantification for variable sensor types and operating ranges. Furthermore, the framework will facilitate fault detection capacities and so-called virtual sensors, a capability to replace faulty sensor measurements with error-corrected values when a sensor fault is detected. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

The Power of Many: An Ensemble Approach to Spectral Similarity

Quantifying the similarity between two mass spectra─a known reference mass spectrum and an unidentified sample mass spectrum─is at the heart of compound identification workflows in gas chromatography–mass spectrometry (GC-MS). The reference spectrum most like the sample is assigned as its identification (provided some quantitative similarity threshold is met, e.g., 80%) and thus accurately measuring similarity is essential. Significant research has gone toward developing metrics for this purpose, each of which has attempted to improve upon existing methods by incorporating GC-MS-specific information (e.g., peak ratios or retention times) or adopting various statistical and algorithmic frameworks. While this active development has led to a plethora of similarity metrics with demonstrated value across different contexts, the unfortunate consequence has been confusion surrounding which metric should be used as a global standard. No such metric is currently accepted as the standard method because different metrics have demonstrated optimal performance in different contexts. In this work, we propose an ensemble approach to spectral similarity scoring that combines the collective information from across existing similarity metrics to form an improved, globally representative similarity metric as a step toward establishing a global standard method. In conclusion, the resulting ensemble metrics are evaluated on over 88,000 spectra of varying complexity and demonstrate improved abilities to accurately rank the correct reference spectrum as the top-matching candidate for a sample relative to the rankings generated by individual similarity scores.

Carbohydrates↗
DOE OSTI2014

Adaptation of the modified Bouc–Wen model to compensate for hysteresis in respiratory motion for the list-mode binning of cardiac SPECT and PET acquisitions: Testing using MRI

Purpose: Binning list-mode acquisitions as a function of a surrogate signal related to respiration has been employed to reduce the impact of respiratory motion on image quality in cardiac emission tomography (SPECT and PET). Inherent in amplitude binning is the assumption that there is a monotonic relationship between the amplitude of the surrogate signal and respiratory motion of the heart. This assumption is not valid in the presence of hysteresis when heart motion exhibits a different relationship with the surrogate during inspiration and expiration. The purpose of this study was to investigate the novel approach of using the Bouc–Wen (BW) model to provide a signal accounting for hysteresis when binning list-mode data with the goal of thereby improving motion correction. The study is based on the authors’ previous observations that hysteresis between chest and abdomen markers was indicative of hysteresis between abdomen markers and the internal motion of the heart. Methods: In 19 healthy volunteers, they determined the internal motion of the heart and diaphragm in the superior–inferior direction during free breathing using MRI navigators. A visual tracking system (VTS) synchronized with MRI acquisition tracked the anterior–posterior motions of external markers placed on the chest and abdomen. These data were employed to develop and test the Bouc–Wen model by inputting the VTS derived chest and abdomen motions into it and using the resulting output signals as surrogates for cardiac motion. The data of the volunteers were divided into training and testing sets. The training set was used to obtain initial values for the model parameters for all of the volunteers in the set, and for set members based on whether they were or were not classified as exhibiting hysteresis using a metric derived from the markers. These initial parameters were then employed with the testing set to estimate output signals. Pearson’s linear correlation coefficient between the abdomen, chest, average of chest and abdomen markers, and Bouc–Wen derived signals versus the true internal motion of the heart from MRI was used to judge the signals match to the heart motion. Results: The results show that the Bouc–Wen model generated signals demonstrated strong correlation with the heart motion. This correlation was slightly larger on average than that of the external surrogate signals derived from the abdomen marker, and average of the abdomen and chest markers, but was not statistically significantly different from them. Conclusions: The results suggest that the proposed model has the potential to be a unified framework for modeling hysteresis in respiratory motion in cardiac perfusion studies and beyond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Sensitivity analysis applied to SiC failure probability in TRISO modeled with BISON

Here, a sensitivity analysis on the failure probability of the Silicon Carbide (SiC) layer in tristructural isotropic (TRISO) nuclear fuel during transient conditions predicted by the BISON fuel performance code is performed. The principal goal of the analysis is to understand the most important parameters dictating SiC failure behavior in BISON during Reactivity Initiated Accidents (RIAs). SiC brittle fracture probability is modeled using Weibull statistics. A total of seven inputs related to SiC failure has been selected for the analysis, including the Weibull statistics parameters, elastic moduli for SiC and Pyrolitic Carbon (PyC) and SiC stress-free temperature. A 1D TRISO BISON model has been established for various reactivity insertions performed at the Nuclear Safety Research Reactor (NSRR). The principal advantage associated with the 1D TRISO model developed in this work is its computational efficiency. The Sobol variance decomposition method is used, and the sensitivity indices are presented for eight different values of the energy deposition. The results show that the two most important parameters impacting the predicted SiC failure probability are the Weibull modulus and the characteristic stress, and a co-variance amongst these parameters is obtained for low reactivity insertions. An additional new finding of this work is that the relative importance of Weibull parameters depends on the energy deposition, and thus reactivity, regime. For low energy depositions, two parameters are of influence on SiC failure probability, while for high energy depositions only one parameter impacts failure probability results. Moreover, optimization of the Weibull modulus and characteristic stress is performed by minimizing the RMSE between BISON failure probability predictions and experimental failure fractions for each energy deposition. This work also demonstrates the validity of the NSRR tests BISON simulations and of the respective sensitivity analysis results as conservative, yet indicative for slower transients characterized by lower deposited energies. Such verification is achieved through the partial extension of the analysis to a group Control Rod Withdrawal reproduced from a previous study. Another novel result of this analysis is that no single set of Weibull parameters can reproduce all reactivity insertion experimental failure results, which is related to the intrinsic nature of SiC failure properties, and a new range for the parameters is proposed to produce a failure probability envelope that encompasses the experimental fractions. Additionaly, this work proposes a new approach for future failure analysis with BISON, consisting in the use of Weibull parameters ranges, rather than fixed sets, along with failure envelopes generation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Evaluation of Hybrid FPOG Applications in Regulated and Deregulated Markets Using HERON

Recent changes in the U.S. energy market, such as low natural gas prices and increased electricity production for variable renewable energy (VRE) sources, have led to an economic crisis for existing light-water reactor (LWR) nuclear power plants (NPP). Many owners and operators of LWRs have elected to decommission these plants rather than continue using them as consistent sources of clean baseload power. This has led to exploration of various possibilities to increase the economic viability of these units, including market restructuring to monetize benefits LWRs already provide to the grid through ancillary markets, load following and economic dispatch, and possible integration of secondary systems directly to the NPP for production of additional products through technologies such as hydrogen electrolysis or water desalination. Previous studies have considered the technologies associated with these Integrated Energy Systems (IES) activities, and the analysis of markets for these secondary products. To analyze the economic viability of various system configurations including IES, especially given the uncertainty surrounding load demand, electricity prices, and the availability of VRE resources, the stochastic technoeconomic analysis package HERON (Heuristic Energy Resource Optimization Network) was released earlier this year as an extension of the risk analysis framework RAVEN (Risk Analysis Virtual Environment). HERON focuses foremost on making the complex uncertainty quantification analysis tools approachable for energy systems analysts, also providing general dispatch optimization algorithms for those workflows. HERON continues to be improved and tested as a significant part of the IES viability analyses performed in this work. HERON is not a capacity expansion model. To consider market and grid energy system development in a variety of scenarios, HERON is best used in coupling with modelling tools such as US-REGEN, which sacrifice some of the uncertainty analysis and resolution of HERON's modelling for the ability to efficiently predict the change in the grid energy system's profile due to economic drivers over decades. HERON can then use this information to explore the economic viability of introducing changes to the predicted outcomes, such as the introduction of an IES. In this work, experts at EPRI using US-REGEN provide six projection scenarios for use in HERON stochastic technoeconomic analysis (STEA) in considering the options available for increasing LWR economic viability through introduction of a hydrogen-centric IES using a high-temperature steam electrolysis plant (HTSE), hydrogen storage, and a constant-rate contracted hydrogen consumer. The results obtained are differential in nature; they do not report expected profits for any configuration, but rather report on the possible increase in the NPV of a configuration with respect to a baseline no-IES configuration. Due to the uncertainty captured in the variable net load of the systems, there is likewise uncertainty in the mean values reported. We consider this viability both in terms of a regulated market, where the energy producers and IES are owned and operated by single entity, as well as a deregulated market, where the IES chooses its bid for electricity generation and is then dispatched by the grid system operator. Results indicate that for deregulated markets, the inclusion of the IES is often statistically beneficial. This is especially true in policies that are not favorable towards nuclear, as nuclear is less often dispatched and is forced to deal with frequent idle capacity. In the nominal case as well as the case of carbon tax policies, inclusion of the IES clearly benefited the economic performance of the NPP. In the regulated case, however, there was a trend towards minimizing the IES, likely due to the optimal sizing performed by US-REGEN of the NPP within the system as well as the lack of penalty for idle capacity at the NPP in the regulated market analyses.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Evaluation of a pair-based, joint-likelihood association approach for regional infrasound event identification

SUMMARY A Bayesian framework for the association of infrasonic detections is presented and evaluated for analysis at regional propagation scales. A pair-based, joint-likelihood association approach is developed that identifies events by computing the probability that individual detection pairs are attributable to a hypothetical common source and applying hierarchical clustering to identify events from the pair-based analysis. The framework is based on a Bayesian formulation introduced for infrasonic source localization and utilizes the propagation models developed for that application with modifications to improve the numerical efficiency of the analysis. Clustering analysis is completed using hierarchical analysis via weighted linkage for a non-Euclidean distance matrix defined by the negative log-joint-likelihood values. The method is evaluated using regional synthetic data with propagation distances of hundreds of kilometres in order to study the sensitivity of the method to uncertainties and errors in backazimuth and time of arrival. The method is found to be robust and stable for typical uncertainties, able to effectively distinguish noise detections within the data set from those in events, and can be made numerically efficient due to its ease of parallelization.

58 GEOSCIENCES↗
DOE OSTI2018

Machine Learning Algorithms for Matching Theories, Simulations, and Observations in Cosmology (Final Project)

The search for simple and elegant equations to explain and predict observable natural phenomena has been a driving force in scientific theory for several centuries. The great successes of Newton's laws, Maxwell's equations, the Schrodinger equation, and the theory of General Relativity have encouraged this perspective, and perhaps a beautiful equation comprising a “theory of everything” remains to be discovered. But now scientists are grappling with systems of ever-increasing complexity, fraught with nonlinearity, and with the quantities of interest often only indirectly observed. Nature may resist a simple description, and the most important discoveries of the next century may be complex theories with countless variables and parameters. The era of big data opens up a promising new approach to scientific discovery in this setting. The predictions of modern theories, even if complex or nonlinear, can be examined through detailed and computationally-intensive simulations that take days or weeks to run and fill petabytes of disk space. The upcoming generation of scientific instruments will provide petabytes of observations at unprecedented resolution and depth. The challenge is to compare data to simulations so as to test theories and identify those that best match the observations. The data are numerous and rich but are subject to noise and various systematic biases. The simulations are large and detailed but are so costly that the effective sample size is small. Moreover, the mapping from data to simulation to the parameters of interest is typically complicated and ill-conditioned. Although these challenges, and the solutions we propose, apply across the sciences, we will ground our efforts by focusing on cosmology. Here, we are truly entering the era of big data in several senses. Ever larger experiments survey larger and larger cosmological volumes, resulting in enormous amounts of pixel data to process and extract the information of cosmological interest. These include large imaging surveys, which collect all the light in broad wavelength ranges over significant fractions of the sky, and spectroscopic surveys, which measure the light as a function of wavelength to map out 3D distances. These experiments are designed to shed light on the big unanswered questions in cosmology. Distilling the petabytes of data into observed quantities for comparison with theory is a daunting task. But these large experiments do not only result in large amounts of data to analyze: if we wish to compare to theories of physics, then we need increasingly extensive predictions for the observed quantities as well, which requires ever-larger cosmological simulations that take substantial -- and at times prohibitive -- computational resources to produce. Our goal is to develop statistical and machine learning methods for using observed and simulated data to advance machine learning with applications cosmology. In particular, we focus our research efforts on the following tasks: The challenge for automated science is that it is computationally impossible and statistically dangerous to consider every possible model in order to find the best one. We will develop Bayesian Optimization based active-learning methods that accelerate both the execution of the simulations and the search for best-fitting parameters. The key idea is to make the simulations adaptive -- across resolutions, time, and parameters -- using the data to search as the simulation runs; Most ML methods operate on simple finite dimensional feature vectors. However, many cosmology and other science applications require ML methods that can operate on more complex objects as inputs or outputs such as functions, distributions, or sets and point clouds. Our goal is to develop efficient ML methods for this problem and demonstrate their applicability in Cosmology, Astrophysics, and other science problems.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Flux Comparison of Master-8 and Ordem 3.1 Modelled Space Debris Population

With ESA’s Meteoroid And Space debris Terrestrial Environment Reference (MASTER-8) model and NASA’s Orbital Debris Engineering Model(ORDEM) 3.1, the two premier orbital debris engineering models have been officially released. The two models come with significant enhancements and now represent state-of-the-art orbital debris modelling for their respective agencies. Both models provide the community with estimates of the space debris environment from low Earth orbit (LEO) up to at least geostationary altitude. The MASTER population is an event-based simulation of all known events that generate debris and objects that are part of the U.S. Space Surveillance Network (SSN) catalog, which provides coverage of objects with diameters down to approximately 10cm in LEO and 1 m in geosynchronous Earth orbit (GEO). Different models are used to simulate the artificial objects and their orbit evolution over time. These models are called “sources” since they assign an origin to each individual object and consist of fragments, solid rocket motor (SRM) remainders, sodium-potassium(NaK) droplets, paintflakes, ejecta, and multi-layer insulation (MLI) fragments. The objects from each source are characterized by having individual release mechanisms, as well as orbital distributions, material composition, size, and mass distributions. Dedicated radar and telescope observation data is used to calibrate the model for objects larger than 1 cm in LEO and larger than 10 cm in GEO. For calibrating the small-sized objects, below 1 cm, impact data from returned surfaces are analyzed. Because the >1cm object population is dominated by fragments, the fragmentation event database was updated to include new events, as well as re-evaluate past events. Special attention was drawn to re-evaluating theFengyun-1C anti-satellite test from 2007 and Cosmos-Iridium collision event from 2009 since these events shape the fragment population because of their severity. After 2009, the two largest fragmentations in terms of number of tracked debris are the Briz-M explosion in 2012 and the NOAA-16explosion in 2015. In total, there are 261 confirmed fragmentations in the database up to November 2016.The baseline population for ORDEM 3.1 is based on the U.S. SSN catalog, and observational datasets from radar, in situ, and optical sources provide a foundation from which the model populations are statistically extrapolated to smaller size regions. These regions are not well-covered by the SSN catalog yet may pose the greatest threat to operational spacecraft. The NASA Standard Satellite Breakup Model is used to generate fragments greater than 1mm from collisions and explosions, and these fragment populations are scaled using ground-based radar data. Specific major debris-producing events, including the Fengyun-1C, Iridium 33, and Cosmos2251 debris clouds, and unique populations, such as NaK droplets, were re-examined, modelled, and added to the ORDEM environment separately. Optical measurement data is used to model the GEO population down to 10 cm. The debris environment is propagated using NASA’s LEO-to-GEO Environment Debris model, and future explosions of intact objects and collisions involving objects greater than 10 cm are assessed statistically. The environment from a few millimetres down to 10 𝜇m is modelled using a special degradation model where small particles are generated from intact spacecraft and rocket bodies, then the populations are scaled to fit in situ cratering data from Space Shuttle returned surfaces. Fragments smaller than 10 cm are differentiated based on material density categories, i.e., high-, medium-, and low-density, to better characterize the potential debris risk posed to upper stages and spacecraft. This paper will discuss the MASTER and ORDEM approaches for modelling populations and compare fluxes for specific orbits, including sun-synchronous, ISS-altitude, geosynchronous transfer, and GEO. In the end, a conclusion is drawn towards the importance of having multiple fundamentally different, yet validated, models to estimate the space debris population.

Andre Horstmann↗