Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistically designed experiments”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Muon Neutrino Disappearance in NOvA with a Deep Convolutional Neural Network Classifier

The NuMI Off-axis Neutrino Appearance Experiment (NOvA) is designed to study neutrino oscillation in the NuMI (Neutrinos at the Main Injector) beam. NOvA observes neutrino oscillation using two detectors separated by a baseline of 810 km; a 14 kt Far Detector in Ash River, MN and a functionally identical 0.3 kt Near Detector at Fermilab. The experiment aims to provide new measurements of Δm 2 and θ23 and has potential to determine the neutrino mass hierarchy as well as observe CP violation in the neutrino sector. Essential to these analyses is the classification of neutrino interaction events in NOvA detectors. Raw detector output from NOvA is interpretable as a pair of images which provide orthogonal views of particle interactions. A recent advance in the field of computer vision is the advent of convolutional neural networks, which have delivered top results in the latest image recognition contests. This work presents an approach novel to particle physics analysis in which a convolutional neural network is used for classification of particle interactions. The approach has been demonstrated to improve the signal efficiency and purity of the event selection, and thus physics sensitivity. Early NOvA data has been analyzed (2.74×10 20 POT, 14 kt equivalent) to provide new best- fit measurements of sin 2 (θ23) = 0.43 (with a statistically-degenerate compliment near 0.60) and |Δm2 | = 2.48 × 10 -3 eV 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Comparison Of 252 Cf Time Correlated Induced Fisssion With AmLi Induced Fission On Fresh MTR Research Reactor Fuel

The effective application of international safeguards to research reactors requires verification of spent fuel as well as fresh fuel. To accomplish this goal various nondestructive and destructive assay techniques have been developed in the US and around the world. The Advanced Experimental Fuel Counter (AEFC) is a nondestructive assay (NDA) system developed at Los Alamos National Laboratory (LANL) combining both neutron and gamma measurement capabilities. Since spent fuel assemblies are stored in water, the system was designed to be watertight to facilitate underwater measurements by inspectors. The AEFC is comprised of six 3 He detectors as well as a shielded and collimated ion chamber. The 3 He detectors are used for active and passive neutron coincidence counting while the ion chamber is used for gross gamma counting. Active coincidence measurement data is used to measure residual fissile mass, whereas the passive coincidence measurement data along with passive gamma measurement can provide information about burnup, cooling time, and initial enrichment. In the past, most of the active interrogation systems along with the AEFC used an AmLi neutron interrogation source. Owing to the difficulty in obtaining an AmLi source, a 252 Cf spontaneous fission (SF) source was used during a 2014 field trail in Uzbekistan as an alternative. In this study, experiments were performed to calibrate the AEFC instrument and compare use of the 252 Cf spontaneous fission source and the AmLi (α,n) neutron emission source. The 252 Cf source spontaneously emits bursts of time-correlated prompt fission neutrons that thermalize in the water and induce fission in the fuel assembly. The induced fission (IF) neutrons are also time correlated resulting in more correlated neutron detections inside the 3 He detector, which helps reduce the statistical errors in doubles when using the 252 Cf interrogation source instead of the AmLi source. In this work, two MTR fuel assemblies varying both in size and number of fuel plates were measured using 252 Cf and AmLi active interrogation sources. This paper analyzes time correlated induced fission (TCIF) from fresh MTR fuel assemblies due to 252 Cf and AmLi active interrogation sources.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2015

Results of Lunar Impact Observations During Geminid Meteor Shower Events

Meteoroids are natural particles with origins from comets, asteroids, and planets from within the solar system. On average, 33 metric tons (73,000 lb) of meteoroids hit Earth everyday with velocities ranging between 20 and 72 km/s. However, the vast majority of these meteoroids disintegrate in the atmosphere and never make it to the ground. The Moon also encounters the same meteoroid flux, but has no atmosphere to stop them from striking the surface. At such speeds even a small meteoroid has incredible energy. A meteoroid with a mass of only 5 kg can excavate a crater over 9 m across, hurling 75 metric tons (165,000 lb) of lunar soil and rock on ballistic trajectories above the lunar surface. Meteoroids with particle sizes as small as 100 micrometer (1 Microgram) can do considerable damage to spacecraft in Earth's orbit and beyond. Impacts can damage thermal protection systems, radiators, windows, and pressurized containers. Secondary effects might include partial penetration or pitting, local deformation, and surface degradation that can cause a failure upon reentry. The speed, mass, density, and flux of meteoroids are important factors for design considerations and mitigation during operations. Lunar operations (unmanned and manned) are also adversely affected by the meteoroid flux. Ejecta from meteoroid impacts is also part of the lunar environment and must be characterized. Understanding meteoroid fluxes and the associated risk of meteoroids impacting spacecraft traveling in and beyond Earth's orbit is the objective of the Meteoroid Environment Office (MEO) located at Marshall Space Flight Center (MSFC). One of the MEO's programs is meteoroid impact monitoring of the Moon. The large collecting area of the night side of the lunar disk provides statistically significant counts of meteoroids that can provide useful information about the flux of meteoroids in the hundreds of grams to kilograms size range. This information is not only important for characterizing the lunar environment associated with larger lunar impactors, but also provides statistical data for verification and improving meteoroid prediction models. Current meteoroid models indicate that the Moon is struck by a sporadic meteoroid with a mass greater than 1 kg over 260 times per year. This number is very uncertain since observations for objects in this mass range are few. Factors of several times, higher or lower, are easily possible. Meteor showers are also present to varying degrees at certain times of the year. The Earth experiences meteor showers when encountering the debris left behind by comets, which is also the case with 2 the Moon. During such times, the rate of shower meteoroids can greatly exceed that of the sporadic background rate for larger meteoroids. Looking for meteor shower impacts on the Moon at about the same time as they occur on Earth will yield important data that can be fed into meteor shower forecasting models, which can then be used to predict times of greater meteoroid hazard on the Moon. The Geminids are one such meteor shower of interest. The Geminids are a major meteor shower that occur in December with a peak intensity occurring usually during the 13th and 14th of the month and appearing to come from a radiant in the constellation Gemini. The Geminids are interesting in that the parent body of the debris stream is an asteroid, which along with the Quadrantids, are the only major meteor showers not originating from a comet. The Geminids parent body, 3200 Phaethon, is about 5 km in diameter and has an orbit that has a 22deg inclination which intersects the main asteroid belt and has a perihelion less than half of Mercury's perihelion distance. Thus, its orbit crosses those of Mars, Earth, Venus, and Mercury. The Geminid debris stream is by far the most massive as compared to the others. When the Earth passes through the stream in mid-December, a peak intensity of approx. equal 120 meteors per hour can be seen. Because of the Geminids' relatively large intensity and unique origin, it is important to monitor and gain information about the Geminids so as to improve their forecasts and understand their contribution to the meteoroid environment in Earth's orbit and at the Moon. It is the purpose of this Technical Memorandum (TM) to document two lunar observing periods coinciding with the Geminid meteor showers that occurred in 2006 and 2010.

Suggs, R. J.↗
DOE OSTI2022

Worldwide Physics-Based Lifetime Prediction of c-Si Modules Due to Solder-Bond Failure

Lifetime prediction of the fielded c-Si solar modules due to location-specific weather conditions has been an important topic of photovoltaic research and the economic viability of solar energy. Data analytic techniques such as the performance ratio method, Statistical clear sky model, and Suns-Vmp methods quantify the degradation from measured data of a solar farm, however, the nonlinear time-dependence and correlated degradations make it difficult to use the empirical degradation rates for ultimate lifetime projection. In this article, we propose a complementary physics-based model to predict the solder bond failure caused by mechanical stress associated with the variations of the temperature. Integrating the worldwide weather information from NASA/NSRDB databases, the model predicts the location-specific output-power degradation and the lifetime of a module due to solder bond failure. The model parameters are calibrated against qualification tests involving thermal cycling of specific batches of modules from a specific technology/manufacturer. The results may be summarized as: 1) Modules installed at higher latitudes show a longer lifetime due to reduced damage accumulation. 2) The reduction of temperature fluctuation close to large bodies of water, such as seashores, increases solder bond lifetime significantly. 3) Relatively speaking, modules installed close to the Tropic of Cancer/Capricorn (23.5 degrees North/South) suffer from a higher solder bond damage and have a shorter lifetime, suggesting a conservative design. This model should serve as a building block of a comprehensive reliability framework that can predict the lifetime of a module that experiences simultaneous and correlated degradation mechanisms involving yellowing, corrosion, and potential-induced degradation.

14 SOLAR ENERGY↗
DOE OSTI2016

Evaluation of two Vaisala RS92 radiosonde solar radiative dry bias correction algorithms

Solar heating of the relative humidity (RH) probe on Vaisala RS92 radiosondes results in a large dry bias in the upper troposphere. Two different algorithms (Miloshevich et al., 2009, MILO hereafter; and Wang et al., 2013, WANG hereafter) have been designed to account for this solar radiative dry bias (SRDB). These corrections are markedly different with MILO adding up to 40 % more moisture to the original radiosonde profile than WANG; however, the impact of the two algorithms varies with height. The accuracy of these two algorithms is evaluated using three different approaches: a comparison of precipitable water vapor (PWV), downwelling radiative closure with a surface-based microwave radiometer at a high-altitude site (5.3 km m.s.l.), and upwelling radiative closure with the space-based Atmospheric Infrared Sounder (AIRS). The PWV computed from the uncorrected and corrected RH data is compared against PWV retrieved from ground-based microwave radiometers at tropical, midlatitude, and arctic sites. Although MILO generally adds more moisture to the original radiosonde profile in the upper troposphere compared to WANG, both corrections yield similar changes to the PWV, and the corrected data agree well with the ground-based retrievals. The two closure activities – done for clear-sky scenes – use the radiative transfer models MonoRTM and LBLRTM to compute radiance from the radiosonde profiles to compare against spectral observations. Both WANG- and MILO-corrected RHs are statistically better than original RH in all cases except for the driest 30 % of cases in the downwelling experiment, where both algorithms add too much water vapor to the original profile. In the upwelling experiment, the RH correction applied by the WANG vs. MILO algorithm is statistically different above 10 km for the driest 30 % of cases and above 8 km for the moistest 30 % of cases, suggesting that the MILO correction performs better than the WANG in clear-sky scenes. Lastly, the cause of this statistical significance is likely explained by the fact the WANG correction also accounts for cloud cover – a condition not accounted for in the radiance closure experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

An adaptive adversarial domain adaptation approach for corn yield prediction

Recently, statistical machine learning and deep learning methods have been widely explored for corn yield prediction. Though successful, machine learning models generated within a specific spatial domain often lose their validity when directly applied to new regions. To address this issue, we designed an unsupervised adaptive domain adversarial neural network (ADANN). Specifically, through domain adversarial training, the ADANN model reduced the impact of domain shift by projecting data from different domains into the same subspace. Also, the ADANN model was designed to be trained in an adaptive way, which guaranteed the model can learn the domain-invariant features and perform accurate yield prediction simultaneously. Informative variables including time-series vegetation indices and sequential weather observations were first collected from multiple data sources and aggregated to the county level. Then, we trained the ADANN model with the extracted features and corresponding reported county-level corn yield from the U.S. Department of Agriculture (USDA). Finally, the trained model was evaluated in four testing years 2016–2019. The U.S. corn belt was used as the study area and counties under study were grouped into two diverse ecological regions. Overall, the experimental results showed that the developed ADANN model had better performance than three other state-of-the-art machine learning models in both local experiments (train and test in the same region) and transfer experiments (train and test in different regions). As the first study using adversarial learning for crop yield prediction, this research demonstrates a novel solution for improving model transferability on crop yield prediction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Conceptual Design of Neutron Sensor Qualification Device

The Advanced Sensors and Instrumentation Program at Idaho National Laboratory has been formulating strategies to qualify sensors for use in nuclear environments, particularly in irradiation experiments and advanced reactors. When qualifying neutron sensors for use in high-temperature environments, the wide range of neutron flux levels and representative energy spectra presents significant challenges. This paper discusses the development of the Neutron Sensor Qualification Device (NQD), which is designed to test neutron sensors in high temperature controlled environments with known neutron spectra, addressing the spatial and spectral complexities of neutron fluxes in reactor cores. The proposed NQD will be situated in the exposure room at the Armed Forces Radiobiology Research Institute, thus affording a unique capability to expose sensors to high neutron and gamma fluxes. To achieve thermal control, the device will utilize a radiation-hardened tube furnace, accommodating multiple sensors and neutron activation dosimetry wires. Titanium, iron, and cobalt dosimeter wires are chosen from the American Society for Testing and Materials and International Reactor Dosimetry and Fusion File libraries as references for providing energy-dependent fluence measurements. The design ensures precise sensor positioning to minimize mutual shielding and flux perturbation, which are evaluated via Monte Carlo N particle Transport Code (MCNP) simulations. These simulations have informed the development of guidelines on sensor placement within the NQD. The NQD is essential to the qualification of neutron sensors for advanced reactor technologies. It enables controlled testing of a statistically significant number of sensors, thereby supporting assessments of sensor performance across various neutron flux levels and temperatures. This paper highlights the detailed planning for the NQD prototype, along with its inaugural irradiation (scheduled for fiscal year [FY] 2025). The results from this initial testing will be fundamental in evaluating the device’s performance and establishing measurement uncertainty for in-pile neutron sensor measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Comparative Analysis of Standard and Advanced USL Methodologies for Nuclear Criticality Safety

The American National Standards Institute/American Nuclear Society national standards 8.1 and 8.24 provide guidance on the requirements and recommendations for establishing confidence in the results of the computerized models used to support operation with fissionable materials. By design, the guidance is not prescriptive, leaving freedom to the analysts to determine how the various sources of uncertainties are to be statistically aggregated. Due to the involved use of statistics entangled with heuristic recipes, the resulting safety margins are often difficult to interpret. Also, these technical margins are augmented by additional administrative margins, which are required to ensure compliance with safety standards or regulations, eliminating the incentive to understand their differences. With the new resurgent wave of advanced nuclear systems, e.g., advanced reactors, fuel cycles, and fuel concepts, focused on economizing operation, there is a strong need to develop a clear understanding of the uncertainties and their consolidation methods to reduce them in manners that can be scientifically defended. In response, the current studies compare the analyses behind four notable methodologies for upper subcriticality limit estimation that have been documented in the nuclear criticality safety literature: the parametric, nonparametric, Whisper, and TSURFER methodologies. Specifically, the work offers a deep dive into the various assumptions of the noted methodologies, their adequacies, and their limitations to provide guidance on developing confidence for the emergent nuclear systems that are expected to be challenged by the scarcity of experimental data. Here, to limit the scope, the current work focuses on the application of these methodologies to criticality safety experiments, where the goal is to calculate a bias, a bias uncertainty, and a tolerance limit for k eff in support of determining an upper subcriticality limit for nuclear criticality safety.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
NASA NTRS1996

Sampling of prenatal and postnatal offspring from individual rat dams enhances animal use without compromising development

To assess prenatal and postnatal developmental status in the offspring of a group of animals, it is typical to examine fetuses from some of the dams as well as infants born to the remaining dams. Statistical limitations often arise, particularly when the animals are rare or especially precious, because all offspring of the dam represent only a single statistical observation; littermates are not independent observations (biologically or statistically). We describe a study in which pregnant laboratory rats were laparotomized on day 7 of gestation (GD7) to ascertain the number and distribution of uterine implantation sites and were subjected to a simulated experience on a 10-day space shuttle flight. After the simulated landing on GD18, rats were unilaterally hysterectomized, thus providing a sample of fetuses from 10 independent uteruses, followed by successful vaginal delivery on GD22, yielding postnatal samples from 10 uteruses. A broad profile of maternal and offspring morphologic and physiologic measures indicated that these novel sampling procedures did not compromise maternal well-being and maintained normal offspring development and function. Measures included maternal organ weights and hormone concentrations, offspring body size, growth, organ weights, sexual differentiation, and catecholamine concentrations.

NASA Discipline Developmental Biology↗
DOE OSTI2023

Tensor-to-scalar ratio forecasts for extended LiteBIRD frequency configurations

LiteBIRD is a planned JAXA-led cosmic microwave background (CMB) B-mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertainty on the tensor-to-scalar ratio, δr, down to δr < 0.001. A key aspect of this performance is accurate astrophysical component separation, and the ability to remove polarized thermal dust emission is particularly important. In this paper we note that the CMB frequency spectrum falls off nearly exponentially above 300 GHz relative to the thermal dust spectral energy distribution, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compared the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the High-Frequency Telescope (HFT) frequency range was shifted logarithmically toward higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measured the tensor-to-scalar ratio r uncertainty and bias using both parametric and minimum-variance component-separation algorithms. When the thermal dust sky model includes a spatially varying spectral index and temperature, we find that the statistical uncertainty on r after foreground cleaning may be reduced by as much as 30–50% by extending the upper limit of the frequency range from 400 to 600 GHz, with most of the improvement already gained at 500 GHz. We also note that a broader frequency range leads to higher residuals when fitting an incorrect dust model, but also it is easier to discriminate between models through higher χ 2 sensitivity. Even in the case in which the fitting procedure does not correspond to the underlying dust model in the sky, and when the highest frequency data cannot be modeled with sufficient fidelity and must be excluded from the analysis, the uncertainty on r increases by only about 5% for a 500 GHz configuration compared to the baseline.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2023

Tensor-to-scalar ratio forecasts for extended LiteBIRD frequency configurations

LiteBIRD is a planned JAXA-led cosmic microwave background (CMB) B -mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertainty on the tensor-to-scalar ratio, δr , down to δr < 0.001. A key aspect of this performance is accurate astrophysical component separation, and the ability to remove polarized thermal dust emission is particularly important. In this paper we note that the CMB frequency spectrum falls off nearly exponentially above 300 GHz relative to the thermal dust spectral energy distribution, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compared the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the High-Frequency Telescope (HFT) frequency range was shifted logarithmically toward higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measured the tensor-to-scalar ratio r uncertainty and bias using both parametric and minimum-variance component-separation algorithms. When the thermal dust sky model includes a spatially varying spectral index and temperature, we find that the statistical uncertainty on r after foreground cleaning may be reduced by as much as 30–50% by extending the upper limit of the frequency range from 400 to 600 GHz, with most of the improvement already gained at 500 GHz. We also note that a broader frequency range leads to higher residuals when fitting an incorrect dust model, but also it is easier to discriminate between models through higher χ 2 sensitivity. Even in the case in which the fitting procedure does not correspond to the underlying dust model in the sky, and when the highest frequency data cannot be modeled with sufficient fidelity and must be excluded from the analysis, the uncertainty on r increases by only about 5% for a 500 GHz configuration compared to the baseline.

U. Fuskeland↗
DOE OSTI2021

Thermonuclear neutron emission from a sheared-flow stabilized Z-pinch

We report the fusion Z-pinch experiment (FuZE) is a sheared-flow stabilized Z-pinch designed to study the effects of flow stabilization on deuterium plasmas with densities and temperatures high enough to drive nuclear fusion reactions. Results from FuZE show high pinch currents and neutron emission durations thousands of times longer than instability growth times. While these results are consistent with thermonuclear neutron emission, energetically resolved neutron measurements are a stronger constraint on the origin of the fusion production. This stems from the strong anisotropy in energy created in beam-target fusion, compared to the relatively isotropic emission in thermonuclear fusion. In dense Z-pinch plasmas, a potential and undesirable cause of beam-target fusion reactions is the presence of fast-growing, “sausage” instabilities. This work introduces a new method for characterizing beam instabilities by recording individual neutron interactions in plastic scintillator detectors positioned at two different angles around the device chamber. Histograms of the pulse-integral spectra from the two locations are compared using detailed Monte Carlo simulations. These models infer the deuteron beam energy based on differences in the measured neutron spectra at the two angles, thereby discriminating beam-target from thermonuclear production. An analysis of neutron emission profiles from FuZE precludes the presence of deuteron beams with energies greater than 4.65 keV with a statistical uncertainty of 4.15 keV and a systematic uncertainty of 0.53 keV. This analysis demonstrates that axial, beam-target fusion reactions are not the dominant source of neutron emission from FuZE. These data are promising for scaling FuZE up to fusion reactor conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Apparent Monte Carlo Source Convergence Problem with BWR Fuel Depleted with Partial Control Blade Insertion

Oak Ridge National Laboratory (ORNL) and the US Nuclear Regulatory Commission (NRC) have initiated a multiyear project to investigate application of burnup credit (BUC) to boiling-water reactor (BWR) fuel in storage and transportation casks and spent fuel pools (SFPs). Ideally, a single combination of all operating parameters including moderator density and axial burnup profiles, control blade exposure, and others can be identified to allow a simple, bounding approach to BWR BUC. One aspect of this investigation involves examining the effect of control blade insertion during operation on discharged assembly reactivity. In the course of these calculations, an apparent source convergence problem was discovered for fuel that has experienced partial blade insertion for the entire lifetime of the assembly. The source convergence problem appears not to impact k{sub eff}, but it manifests in predictions of the axial fission density profile. This phenomenon is of interest in identifying areas of concern with BWR BUC since it has not been studied as thoroughly as BUC for pressurized-water (PWR) reactor fuel. There may also be interest in this configuration for further study of a realistic source Monte Carlo convergence problem as opposed to some earlier manufactured source convergence problems through careful selection of initial conditions. This summary presents the depletion conditions used to generate the problem, a description of the models, a background based on well-behaved calculations, the identification of the problem, and attempts made to resolve the problem. Only a brief description of the depletion conditions and models is possible within the limited scope of this summary. Finally, conclusions are presented. This summary identifies a source convergence problem in the process of investigating BUC for BWR spent fuel. The depletion scenario which created the depleted fuel and burnable absorber compositions included a control blade inserted in the bottom 14 nodes throughout the entire depletion. This insertion is slightly over half the length of the assembly, and it appears to create two regions of nearly equal reactivity. As expected, one region is near the top of the assembly, and the other region is near the middle of the assembly. This node includes significant {sup 239}Pu generation as a result of rodded operation, but it is also the top node of the full lattice zone. While the calculated k{sub eff} values appear to be statistically equivalent, the axial neutron production profile varies significantly among independent simulations of the same model. A range of standard approaches has been used to attempt to resolve the source convergence problem and determine the true axial neutron production profile. This profile is of interest in cask criticality safety analysis because it influences the interaction among assemblies loaded in the cask. The forthcoming SCALE 6.2 release will include source convergence testing based on Shannon entropy, which may assist in identifying these types of problems in the future. It should be noted that this source convergence problem results from the depletion conditions and the final distribution of fissile and absorber material within the discharged fuel assembly. The results, while entirely generated within the SCALE code system, are representative of any Monte Carlo program attempting to calculate k{sub eff} for such a system. The most important conclusion of this work is that the BUC criticality safety community must be aware of the potential for this type of source convergence problem. Analysts must approach such partially rodded scenarios with caution and critically examine the results. Fortunately, it is extremely unlikely that BWR assemblies would experience the bladed operation necessary to create this particular source convergence problem. These results may also have some bearing on the analysis of PWR BUC for new plant designs such as AP1000 that incorporate significant rodded operations. In a broader context, the source convergence problem identified here could have been magnified if it had occurred in a Monte Carlo depletion simulation. It may be possible to identify this type of problem in a depletion calculation by comparing the apparently converged fluxes between steps to ensure that no oscillations or other serious shifts are occurring. The depletion calculations performed in this work were performed with a two-dimensional deterministic code and are therefore not subject to this problem. These results may also raise questions about the validity of 2D depletion models with extruded geometry in the 3D k{sub eff} calculation. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Understanding Dynamical Systems in High-Dimensional Parameter Spaces (Final Report)

In the last few decades, computational models, especially of the complex dynamical systems of greatest interest, i.e., climate, high energy density science, combustion, etc. have become increasingly accurate and sophisticated. Consequently, they are now routinely used to predict outcomes of experiments or to inform design decisions. However, many such simulations depend on a (large) number of parameters in order to, for example, approximate physical phenom whose exact solution is unknown or describe material properties in situations with little or no experimental data. Choosing these parameters and understanding how the outcomes of such a simulation changes with respect to the inputs are rapidly becoming some of the primary challenges in many application areas. Simulations have long been too complex for humans to fully comprehend and thus are often treated as black-box systems. The inputs to the systems are typically a collection of parameters with expert-selected ranges, and the outputs are various quantities of interest (QoIs) describing the performance of the system. Therefore, the fundamental challenge can be described as one of understanding a high-dimensional (many input parameters), multi-variate (many QoIs) function. Unfortunately, existing analysis tools are only of limited use for this task as most techniques can neither deal with many input dimensions nor offer opportunities to jointly study several QoIs. Furthermore, data analysis has long relied on an exploratory process in which scientist first visualize the available data, for example, using scatter plots or 3D renderings, form new hypotheses on some phenomena of interest, and subsequently test these hypotheses using dedicated analysis or statistics approaches. However, for high-dimensional functions this exploration is difficult if not impossible using existing tools as there are too many QoIs to simultaneously observe, and visualizations of high-dimensional spaces are abstract and difficult to understand if they exist at all. The goal of this project is to build new analysis techniques and tool that offer high level insights for exploring high-dimensional scientific data. And to extend existing approach beyond single scalar functions and instead consider multiple QoIs simultaneously.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Physics‐based iterative reconstruction for dual‐source and flying focal spot computed tomography

Purpose For single‐source helical Computed Tomography (CT), both Filtered‐Back Projection (FBP) and statistical iterative reconstruction have been investigated. However, for dual‐source CT with flying focal spot (DS‐FFS CT), a statistical iterative reconstruction that accurately models the scanner geometry and acquisition physics remains unknown to researchers. Therefore, our purpose is to present a novel physics‐based iterative reconstruction method for DS‐FFS CT and assess its image quality. Methods Our algorithm uses precise physics models to reconstruct from the native cone‐beam geometry and interleaved dual‐source helical trajectory of a DS‐FFS CT. To do so, we construct a noise physics model to represent data acquisition noise and a prior image model to represent image noise and texture. In addition, we design forward system models to compute the locations of deflected focal spots, the dimension, and sensitivity of voxels and detector units, as well as the length of intersection between x‐rays and voxels. The forward system models further represent the coordinated movement between the dual sources by computing their x‐ray coverage gaps and overlaps at an arbitrary helical pitch. With the above models, we reconstruct images by an advanced Consensus Equilibrium (CE) numerical method to compute the maximum a posteriori estimate to a joint optimization problem that simultaneously fits all models. Results We compared our reconstruction with Siemens ADMIRE, which is the clinical standard hybrid iterative reconstruction (IR) method for DS‐FFS CT, in terms of spatial resolution, noise profile, and image artifacts through both phantoms and clinical scan datasets. Experiments show that our reconstruction has a higher spatial resolution, with a Task‐Based Modulation Transfer Function (MTF task ) consistently higher than the clinical standard hybrid IR. In addition, our reconstruction shows a reduced magnitude of image undersampling artifacts than the clinical standard. Conclusions By modeling a precise geometry and avoiding data rebinning or interpolation, our physics‐based reconstruction achieves a higher spatial resolution and fewer image artifacts with smaller magnitude than the clinical standard hybrid IR.

Wang, Xiao↗
DOE OSTI2023

Manganese effects on plant residue decomposition and carbon distribution in soil fractions depend on soil nitrogen availability

Recent studies have highlighted the critical role of manganese (Mn) in plant litter decomposition and soil organic carbon (C) cycling in forest ecosystems. Long term nitrogen (N) deposition and N fertilization can increase soil acidity and mobilize bioavailable Mn (Mn 2+ ) in soil. However, no studies have examined the interactive effect of N and Mn fertilization on litter decomposition and carbon distribution in agricultural soils, despite agroecosystems being subject to both N and Mn management. We hypothesized that increased soil N and Mn availability would accelerate plant residue decomposition and transfer of its C to mineral-associated organic matter (MAOM), and that the combined effect of Mn and N enrichment would be greater than the individual effect. Here, we conducted a laboratory incubation experiment by adding 13 C-labeled residue of perennial grass Glyceria striata (Lam.) to agricultural soils that had received 225 kg N ha –1 yr –1 for 27 years (N 1 ) and comparable soils that received no N (N 0 ). Before the experiment, these soils also received three levels of dissolved Mn 2+ , designated M 0 (no additional Mn), M 1 (50 mg kg –1 ), or M 2 (250 mg kg –1 ). We measured total CO 2 production as well as distribution of 13 C from the residue into CO 2 , particulate organic matter (POM), MAOM, and dissolved organic carbon (DOC) over a 1-year period. Manganese amendments significantly increased CO 2 production from residue decomposition in the N 1 soil, but no such effect was observed in the N 0 soil. Manganese also accelerated the loss of residue-derived C from POM and DOC, but increased its recovery in MAOM. However, the positive effect of added Mn in decomposition and recovery in MAOM in the presence of N fertilization occurred only during the initial 30-day decomposition period, where M 2 showed a 12% increase in cumulative CO 2 production from residue, 8% increase in POM loss, and 43% increase in recovery of residue C in MAOM compared to M 0 . For M 1 , only CO 2 emission from residue was significantly higher than Mo during this period. At 365 days M 2 showed 8% increase in CO 2 production, 1% increase in POM loss, and 16% increase in recovery of residue C in MAOM compared to M 0 , but none of these were statistically significant (p < 0.05). This study adds to the growing evidence that increasing Mn availability enhances plant litter decomposition. However, the occurrence and magnitude of Mn-induced stimulation of decomposition is context specific. Further investigation with greater temporal resolution, involving a multitude of litter and soil types and including microbial compositional and functional characterization, is recommended to fully elucidate the interactive role of Mn and N on C cycling.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2005

Case Study of the Space Shuttle Cockpit Avionics Upgrade Software

The purpose of the Space Shuttle Cockpit Avionics Upgrade project was to reduce crew workload and improve situational awareness. The upgrade was to augment the Shuttle avionics system with new hardware and software. An early version of this system was used to gather human factor statistics in the Space Shuttle Motion Simulator of the Johnson Space Center for one month by multiple teams of astronauts. The results were compiled by NASA Ames Research Center and it was was determined that the system provided a better than expected increase in situational awareness and reduction in crew workload. Even with all of the benefits nf the system, NASA cancelled the project towards the end of the development cycle. A major success of this project was the validation of the hardware architecture and software design. This was significant because the project incorporated new technology and approaches for the development of human rated space software. This paper serves as a case study to document knowledge gained and techniques that can be applied for future space avionics development efforts. The major technological advances were the use of reflective memory concepts for data acquisition and the incorporation of Commercial off the Shelf (COTS) products in a human rated space avionics system. The infused COTS products included a real time operating system, a resident linker and loader, a display generation tool set, and a network data manager. Some of the successful design concepts were the engineering of identical outputs in multiple avionics boxes using an event driven approach and inter-computer communication, a reconfigurable data acquisition engine, the use of a dynamic bus bandwidth allocation algorithm. Other significant experiences captured were the use of prototyping to reduce risk, and the correct balance between Object Oriented and Functional based programming.

Ferguson, Roscoe C.↗
DOE OSTI2020

A Bayesian Approach for Statistical–Physical Bulk Parameterization of Rain Microphysics. Part I: Scheme Description

A novel framework is proposed for the bulk parameterization of rain microphysics: the Bayesian Observationally Constrained Statistical–Physical Scheme (BOSS). It is designed to facilitate direct constraint by observations using Bayesian inference. BOSS combines existing process-level microphysical knowledge with flexible process rate formulations and parameters constrained by observations within a Bayesian framework. Furthermore, using a raindrop size distribution (DSD) normalization method that relates DSD moments to one another via generalized power series, generalized multivariate power expressions are derived for the microphysical process rates as functions of a set of prognostic DSD moments. The scheme is flexible and can utilize any number and combination of prognostic moments and any number of terms in the process rate formulations. This means that both uncertainty in parameter values and structural uncertainty associated with the process rate formulations can be investigated systematically, which is not possible using traditional schemes. In this paper, BOSS is compared to two- and three-moment versions of a traditional bulk rain microphysics scheme (denoted as MORR). It is demonstrated that some process formulations in MORR are analytically equivalent to the generalized power expressions in BOSS using one or two terms, while others are not. BOSS is able to replicate the behavior of MORR in idealized one-dimensional rainshaft tests, but with a much more flexible and systematic design. Part II of this study describes the application of BOSS to derive rain microphysical process rates and posterior parameter distributions in Bayesian experiments using Markov chain Monte Carlo sampling constrained by synthetic observations.

58 GEOSCIENCES↗