SEARCH · Engineering Papers
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.
This poster provides an overview of the application of the Solar Uncertainty Integrator (SUNI) software to atmospheric radiation measurement (ARM) solar radiation measurements.
We present a multi-redshift Baryon Acoustic Oscillations (BAO) analysis of the DESI Data Release 2 (DR2) Lyman-$α$ (Ly$α$) forest, splitting the forest auto-correlation and its cross-correlation with quasars into three redshift bins. We obtain BAO measurements at effective redshifts $z_{\rm eff} = 2.13$, $2.40$, and $2.81$ with $\sim2.0$--$2.5\%$ precision per bin in the radial and transverse directions, corresponding to $\sim1.1$--$1.2\%$ precision for the isotropic BAO measurement. Using the same data products and modeling framework as the DESI DR2 Ly$α$ BAO analysis, we validate the pipeline on $400$ synthetic datasets and find unbiased BAO recovery with well-calibrated uncertainties. The measurements show an increase in the isotropic dilation parameter $D_V/r_d$ from $30.26\pm0.39$ to $32.22\pm0.47$ and in the Alcock-Paczyński parameter $D_M/D_H$ from $3.96\pm0.15$ to $5.63^{+0.22}_{-0.24}$. The Hubble distance $D_H/r_d$ decreases from $9.40\pm0.20$ to $7.22\pm0.17$, providing a direct measurement of the expansion history consistent with $Λ$CDM and the expected matter-dominated scaling, with $H(z)\propto(1+z)^n$ giving $n=1.34\pm0.16$. The redshift split also provides a self-consistent measurement of clustering evolution: the Ly$α$ forest bias evolves as $(1+z)^γ$ with $γ_α=3.05\pm0.16$, the RSD parameter has a redshift evolution described by $γ_β=-0.97\pm0.26$, and the quasar bias evolves with $γ_Q=1.56\pm0.23$, consistent with independent quasar clustering measurements. Combining these three-bin BAO measurements with DESI DR2 galaxy and quasar BAO measurements yields cosmological constraints consistent with the single-bin Ly$α$ BAO analysis in flat $Λ$CDM and $w_0w_a$CDM and improves curvature constraints by $\sim12\%$ in $Λ$CDM$+Ω_\mathrm{K}$.
This paper discusses the on-orbit performance and performance degradation of the International Space Station P6 solar array wings (SAWs) from the period of December 2000 through February 2003. Data selection considerations and data reduction methods are reviewed along with the approach for calculating array performance degradation based on measured string shunt current levels. Measured degradation rates are compared with those predicted by the computational tool SPACE and prior degradation rates measured with the same SAW technology on the Mir space station. Initial results show that the measured SAW short-circuit current is degrading 0.2 to 0.5 percent per year. This degradation rate is below the predicted rate of 0.8 percent per year and is well within the 3 percent estimated uncertainty in measured SAW current levels. General contributors to SAW degradation are briefly discussed.
An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. Here in this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.
The NOvA Test Beam Program was commissioned to complement the full-scale NOvA experiment in its goal of measuring neutrino oscillations. The NOvA Experiment consists of a 300-ton Near Detector and a 14-kton Far Detector, both exposed to the same beam, which can be operated in neutrino or antineutrino mode from Fermilab at an off-axis angle of 14.6~mrad. These detectors can discriminate between $\nu_\mu$-CC and $\nu_e$-CC neutrino (or corresponding antineutrino) events, enabling appearance and disappearance measurements of these flavors across the \SI{810}{\kilo\meter} near-to-far detector baseline. From these measurements, the underlying parameters of the three-flavor PMNS oscillations framework are measured. The NOvA Test Beam uses identical technology to the full-scale detectors, but scaled down to 30~tons. It is exposed to a beam of charged particles closely matched to the daughter products of NOvA (anti)neutrino interactions in both energy and species. An instrumented beamline determines both particle momentum and species upstream of the NOvA detector, allowing for a crisp measurement of the NOvA detector response to these known inputs. This dissertation presents an analysis of the electromagnetic response of the NOvA hardware and establishes a framework for incorporating these results to enhance the physics reach of the experiment. Prior to this work, the electromagnetic energy response of the NOvA detector relied primarily on simulation-driven calibration together with a globally applied energy scale uncertainty. The measurements presented here provide direct test beam constraints on the detector response, enabling particle-dependent and potentially energy-dependent uncertainty models derived from controlled electromagnetic inputs.
This report documents the initial algorithm that could be used by the Waste Treatment and Immobilization Plant (WTP) in batching high-level waste (HLW) and glass-forming chemicals (GFCs) in the HLW melter feed preparation vessel (MFPV) (HFP-VSL-00001 and -00005). Not all Hanford tank waste can be accommodated by the models developed for this report and significant expansion of the model boundaries could be achievable to reduce the WTP mission life and total canister production count. The immobilized HLW (IHLW) must meet a series of constraints to be acceptable for disposal in the Monitored Geologic Repository, which are contained in the Specification 1 of the Contract (DOE 2000), the Waste Acceptance Product Specifications (WAPS, DOE 1996), and the Waste Acceptance System Requirements Document (WASRD, DOE 2007). The IHLW Waste Form Compliance Plan (WCP, 24590-HLW-PL-RT-07-0001, Rev 3) specifies that the formulation algorithm will be developed and used to comply with the constraints associated with glass composition and properties. This report is not an engineering calculation, does not provide design input, and is not an engineering study. Algorithm inputs include the chemical analyses of the blended HLW in the HLW blend vessel (HBV) (HLP-VSL-00028, the volume and composition of the MFPV heel, the volume and composition of the MFPV after waste addition, the volume and composition of MFPV batch after GFC addition, the compositions of individual GFCs, and the mass of glass in each canister. In addition to these inputs, uncertainties in the HLW composition and processing parameters are included in the algorithm. Using the above inputs, the algorithm calculates the following outputs: 1) the volume of HLW to be transferred from the HBV to the MFPV, 2) the mass of each GFC for addition to the MFPV, 3) the composition of the glass that will be produced along with uncertainties, and 4) the predicted properties, with associated uncertainties, of the resulting IHLW. The algorithm uses the property-composition models to calculate properties with associated uncertainties and compares them with various constraints to ensure that a processable feed is formulated and a compliant IHLW is produced. The GFC additions are determined using an optimization approach to provide high confidence that the HLW glass will meet all product quality requirements and key processing constraints. For most HLW batches there are many possible glass compositions that meet all constraints. In these cases, the glass composition is optimized for a series of target component concentrations and target property values. The algorithm also incorporates process measurement and product quality uncertainties, based on the work of Piepel et al. (2005). Estimates of the various process and measurement uncertainties that affect glass compositions and predicted glass properties have been previously reported (Piepel et al. 2005, 2006) and the impacts of these estimated uncertainties on the IHLW composition envelope that meets product quality and processing-related properties with sufficient confidence were evaluated. The details of work performed to date to develop this initial GFC addition and batching algorithm are summarized in Sections 4 and 5. An example data set is used to illustrate the calculations of the algorithm summarized in Section 6. Finally, in Section 7, there is a statement of the required work to achieve a final operational IHLW formulation control algorithm. This report is not an engineering calculation, does not provide design input, and is not an engineering study.
Experimenters at Los Alamos National Laboratory (LANL) measure the prompt neutron decay constant for many experiments at the National Criticality Experiments Research Center (NCERC) to infer reactivity and the effective neutron multiplication factor. These quantities are very important for nuclear criticality safety and validating nuclear data. Uncertainty in measures of criticality of an experimental configuration can be determined prior to physically performing the experiment by applying first order perturbation theory to Monte Carlo codes, such as MCNP®. The first order perturbation theory produces first derivatives of some nuclear parameter to nuclear data (e.g., cross section data). This first derivative is commonly referred to as a sensitivity coefficient. Currently, the MCNP® software has the capability of computing effective neutron multiplication factor sensitivity coefficients to cross section data. This work builds off of this MCNP® capability and the first order perturbation theory to provide a method of calculating sensitivity coefficients for the prompt neutron decay constant at or near delayed critical to cross section data. The prompt neutron decay constant sensitivity coefficient calculated in this work does not depend on any modification of the MCNP® source code. Prompt neutron decay constant sensitivity coefficient calculations can be used to infer reactivity and effective neutron multiplication factor sensitivity coefficient values as well. By investigating the trends of prompt neutron decay constant sensitivity coefficients for nuclide-reaction pairs across energy spectra, experiments can be designed to maximize or minimize the uncertainty in the prompt neutron decay constant in a particular energy region, which can lead to further optimization studies. Prompt neutron decay constant sensitivity coefficients will be calculated for the Jezebel benchmark. Subsequently, these sensitivity coefficients will be used in a data assimilation process to determine if there is or are optimal experiments that can be performed to provide insight into adjustments of uncertain/inaccurate cross section data. Specifically, the effect of the prompt neutron decay constant sensitivity coefficients on the nuclear data-induced uncertainty in the effective neutron multiplication factor will be examined
To advance MSR MC&A practices, Argonne National Laboratory developed and optimized robust flow enhanced electrochemical sensors (FEES) and multielectrode array voltammetry sensors (MAVS) which provide accurate near-real time measurements of actinides in molten salts. Flow-enhanced electrochemical sensors are installed directly into MSR flow conduits for inline salt chemistry measurement, while MAVS are deployed in quiescent salt conditions enabling online species concentration determination in stationary salt vessels. This report summarizes efforts to improve electrochemical sensor technological readiness through (1) sensor testing in complex, high concentration molten salt systems to demonstrate low uncertainty actinide measurements in MSR-representative solutions and (2) demonstrations of sensor performance in challenging real-world conditions in collaboration with partner MSR institutions.
Data-reduction procedures for determining free stream and post-normal shock kinetic and thermodynamic quantities are derived. These procedures are applicable to imperfect real air flows in thermochemical equilibrium for temperatures to 15 000 K and a range of pressures from 0.25 N/sq m to 1 GN/sq m. Although derived primarily to meet the immediate needs of the 6-inch expansion tube, these procedures are applicable to any supersonic or hypersonic test facility where combinations of three of the following flow parameters are measured in the test section: (1) Stagnation pressure behind normal shock; (2) freestream static pressure; (3) stagnation point heat transfer rate; (4) free stream velocity; (5) stagnation density behind normal shock; and (6) free stream density. Limitations of the nine procedures and uncertainties in calculated flow quantities corresponding to uncertainties in measured input data are discussed. A listing of the computer program is presented, along with a description of the inputs required and a sample of the data printout.
The transition to a more circular economy (CE) is complex and hard to predict, including in sustainable energy technologies. Many sources of uncertainty make it challenging to model CE scenarios and their potential benefits. As an example, the high variability in costs and revenues of different recycling options makes future wind turbine blade recycling highly uncertain. To better understand this challenge, the circular economy, life cycle assessment and visualization (CELAVI) framework - a discrete event simulation and life cycle assessment framework - is modified to incorporate uncertainty analysis capabilities. Moreover, a 3-step procedure that covers different aspects of uncertainty in CE studies and includes a Monte-Carlo analysis is proposed. The procedure is tested in a case study on wind turbine blade recycling using CELAVI. Results highlight that grinding and landfilling costs are the most influential parameters for wind power circularity. The model's output coefficients of variation (when input parameter uncertainties are propagated) are between 92% and 384% depending on the indicator. The approach developed in this study may help researchers and decision-makers who study circularity prioritize their data collection effort. Finally, our method contributes to a mounting yet critical body of research: the measurement of uncertainty in circularity transitions.
Efficient analysis of computer models is essential for the validation and uncertainty quantification of those models. Surrogate models, and Gaussian processes in particular, are a common and powerful approach to analyzing computer models that treat computer models as a black-box function. Gaussian processes form a Bayesian model over a space of functions that gives a measure of uncertainty about the computer model output at unobserved locations and a framework for sequential sampling. This tutorial will show how to fit Gaussian processes on a series of test functions and apply sequential design techniques to estimate extrema, level sets, and reliabilities of those test functions.
The GMn experiment (E12-09-019) was conducted at Jefferson Laboratory from late 2021 into early 2022. The goal was to make high-precision measurement of the neutron’s magnetic form factor (GMn) at multiple kinematic points, including Q2 = 3.5, 4.5, 6.5, 8.5, 10, 12, 13.5, 16, and 18 (GeV/c)2. Limited data exist for GMn in the region up to about Q2 = 10 (GeV/c)2, with existing data having large systematic uncertainty. In this experiment, the ratio method was employed to reduce systematic uncertainty by measuring the ratio of neutron and proton yields. The experiment took place at Jefferson Laboratory in Hall A, where the BigBite spectrometer was used to detect the scattered electrons, while the HCal in the SuperBigbite spectrometer was used to detect both neutrons and protons. The protons were deflected slightly upwards with the use of a large-aperture dipole magnet named BigBen, allowing for enhanced particle identification. Extraction of GMn requires taking the ratio of proton and neutron yields. Analysis efforts are still currently underway to refine corrections to the data, including detector efficiencies, radiative corrections, neutron and proton mass identification, and charge exchange.
Distribution System State Estimation (DSSE) is frequently constrained by limited real-time measurements, the uncertainties introduced by distributed energy resources, and the presence of bad data. To address them, this paper proposes an enhanced Multi-Fidelity Bayesian Neural Network (MFBNN) DSSE approach. A low-fidelity layer based on a Deep Neural Network (DNN) is first pre-trained on pseudo-measurement data to learn fundamental state features. Subsequently, a high-fidelity Bayesian Neural Network (BNN) layer leverages limited but high-quality real-time measurements to refine these features, thereby achieving accurate DSSE. Additionally, the deep SHapley Additive exPlanation (SHAP) is developed to quantify the influence of measurement data on DSSE through dual perspectives of global feature importance and local nodal contributions, establishing a hierarchical explainability framework for machine learning-based DSSE. Comparative studies conducted on the IEEE 13-bus system and a real-world 2135-node system from Dominion Energy demonstrate that the proposed method excels in estimation accuracy, even under situations of high noise levels, bad data, and missing data. Further comparisons with Weighted Least Squares (WLS) and other machine learning-based DSSE approaches verify that the proposed framework offers higher accuracy, improved interpretability, and enhanced robustness.
The Electron-Ion Collider will be the first collider to use both polarized electron beams and polarized protons and light ions . This will offer unique opportunities to study the structure of protons and nuclei and answer fundamental questions in QCD. The uncertainties on the polarization measurement translate directly into the uncertainties of final physics observables. Hence, a precise measurement of the hadron beam polarization and a good control of the uncertainties are critical for the success of the spin program at the EIC. The requirements for beam polarimetry are non-destructive with uncertainty less than 1%. At the Electron Storage Ring (ESR) and the Rapid Cycling Synchrotron (RCS), the electron beam polarization will be measured using wellestablished Compton polarimetry techniques. However, the EIC Compton polarimeter will face unique challenges, demanding further developments. A Mott polarimeter will also be employed at the source for initial polarization measurements. Both longitudinal and transverse polarization will be measured, with the capability to monitor polarization on a bunch-by-bunch basis. Achieving these precise measurements will be critical to controlling systematic uncertainties and ensuring the overall success of the EIC’s physics objectives.
In spite of increased interest in thermal Hall measurements for the analysis of insulating quantum materials, there remains large uncertainty in such measurements due to contact misalignment. In this paper, we propose that sample geometry and uncertain boundary conditions provide an additional source of uncertainty in the measurement of $D_{xy}$ or $κ_{xy}$. By running simple simulations in an open source finite-element solver, we demonstrate that measured $D_{xy}$ can be changed by a non-negligible fraction in samples with similar width and length. This geometric corrective factor depends on the distinction between a uniform heat flow and constant temperature boundary couplings to a bath. Sample geometry and boundary conditions can be accounted for through simulation or by using rectangular samples to make thermal Hall measurements more reliable and reproducible. Finally, we detail a contactless optical method for measuring $D_{xy}$ based on an existing photothermal microscope technique. This method is insensitive to the longitudinal diffusivity pollution caused by contact misalignment.
The complexity of modern military systems, as well as the cost and difficulty associated with experimentally verifying system and subsystem design makes the use of high-fidelity based simulation a future alternative for design and development. The predictive ability of such simulations such as computational fluid dynamics (CFD) and computational structural mechanics (CSM) have matured significantly. However, for numerical simulations to be used with confidence in design and development, quantitative measures of uncertainty must be available. The AVT 147 Symposium has been established to compile state-of-the art methods of assessing computational uncertainty, to identify future research and development needs associated with these methods, and to present examples of how these needs are being addressed and how the methods are being applied. Papers were solicited that address uncertainty estimation associated with high fidelity, physics-based simulations. The solicitation included papers that identify sources of error and uncertainty in numerical simulation from either the industry perspective or from the disciplinary or cross-disciplinary research perspective. Examples of the industry perspective were to include how computational uncertainty methods are used to reduce system risk in various stages of design or development.
The goal of this program is to develop an electrical resistance strain gage system which will accurately measure the static strains of superalloy blades and vanes in gas turbine engines running on a test stand. Accurate knowledge of these strains is essential to reaching the goals of the HOST program in the selection and experimental verification of the various theoretical models developed to understand and improve the performance of these engines. The specific objective is to develop a complete system capable of making strain measurements of up to + or - 10 percent of full scale during a 50 hour period at temperatures as high as 1250 K. In addition to survival and stability, attaining a low temperature coefficient of resistance, of the order of 20 ppm/K or less, was a major goal. This requirement arises from the presently unavoidable uncertainties in measurement of the exact temperatures inside gas turbines for use in making corrections for apparent strain due to temperature.
In this study, shock initiation experiments have been carried out on the polymer-bonded explosive PBX 9012 [nominally 90% octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazine (HMX), 10% vinylidene–hexafluoropropylene copolymer (Viton A) by weight] in order to provide calibration data for the explosive’s unreacted equation-of-state (EOS) and shock-to-detonation transition for reactive burn rate calibration. The input pressures covered the range of 1.86–4.43 GPa. This provided run-to-detonation depths ranging from > 22.3 mm for the lowest pressure shot to 4.91 mm at the highest pressure. The relative shock sensitivity of PBX 9012 is compared to other HMX-explosives in terms of the Pop-plot, showing that the studied explosive is more sensitive than other similar HMX-based counterparts (with notable exceptions). The uncertainty in the shock velocity determinations from the shock tracker measurements are also investigated, yielding new uncertainty measures in the generated Hugoniot data and run-to-detonation coordinates. Finally, the unreacted equation-of-state is determined using a linear U s – u p Mie–Grüneisen relation and the Davis reactants EOS analytical form, the latter being more suited for reactive burn modeling.