Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Systematic error”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

Discretization error estimation and exact solution generation using the method of nearby problems.

The Method of Nearby Problems (MNP), a form of defect correction, is examined as a method for generating exact solutions to partial differential equations and as a discretization error estimator. For generating exact solutions, four-dimensional spline fitting procedures were developed and implemented into a MATLAB code for generating spline fits on structured domains with arbitrary levels of continuity between spline zones. For discretization error estimation, MNP/defect correction only requires a single additional numerical solution on the same grid (as compared to Richardson extrapolation which requires additional numerical solutions on systematically-refined grids). When used for error estimation, it was found that continuity between spline zones was not required. A number of cases were examined including 1D and 2D Burgers equation, the 2D compressible Euler equations, and the 2D incompressible Navier-Stokes equations. The discretization error estimation results compared favorably to Richardson extrapolation and had the advantage of only requiring a single grid to be generated.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Isospin-breaking corrections to light-meson leptonic decays from lattice simulations at physical quark masses

The decreasing uncertainties in theoretical predictions and experimental measurements of several hadronic observables related to weak processes, which in many cases are now smaller than O(1%), require theoretical calculations to include subleading corrections that were neglected so far. Precise determinations of leptonic and semi-leptonic decay rates, including QED and strong isospin-breaking effects, can play a central role in solving the current tensions in the first-row unitarity of the CKM matrix. In this work we present the first RBC/UKQCD lattice calculation of the isospin-breaking corrections to the ratio of leptonic decay rates of kaons and pions into muons and neutrinos. The calculation is performed at fixed lattice spacing (a –1 ≃ 1.730 GeV) on a 48 3 × 96 volume with N f = 2 + 1 dynamical quarks close to the physical point and domain wall fermions in the Möbius formulation are employed. Long-distance QED interactions are included according to the QED L prescription and the crucial role of finite-volume electromagnetic corrections in the determination of leptonic decay rates, which produce a large systematic uncertainty, is extensively discussed. Finally, we study the different sources of uncertainty on |V us |/|V ud | and observe that, if finite-volume systematics can be reduced, the error from isospin-breaking corrections is potentially sub-dominant in the final precision of the ratio of the CKM matrix elements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Data-driven high-dimensional statistical inference with generative models

Crucial to many measurements at the LHC is the use of correlated multi-dimensional information to distinguish rare processes from large backgrounds, which is complicated by the poor modeling of many of the crucial backgrounds in Monte Carlo simulations. In this work, we introduce HI-SIGMA, a method to perform unbinned high-dimensional statistical inference with data-driven background distributions. In contradistinction to many applications of Simulation Based Inference in High Energy Physics, HI-SIGMA relies on generative ML models, rather than classifiers, to learn the signal and background distributions in the high-dimensional space. These ML models allow for interpretable inference while also incorporating model errors and other sources of systematic uncertainties. We showcase this methodology on a simplified version of a di-Higgs measurement in the bbγγ final state, where the di-photon resonance allows for background interpolation from sidebands into the signal region. We demonstrate that HI-SIGMA provides improved sensitivity as compared to standard classifier-based methods, and that systematic uncertainties can be straightforwardly incorporated by extending methods which have been used for histogram based analyses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm

The origin of the Finite Grid Instability (FGI) is studied by resolving the dynamics in the 1D electrostatic Particle-In-Cell (PIC) model in the spectral domain at the single particle level and at the collective motion level. The spectral fidelity of the PIC model is contrasted with the underlying physical system or the gridless model. The systematic spectral phase and amplitude errors from the charge deposition and field interpolation are quantified for common particle shapes used in the PIC models. Lastly, it is shown through such analysis and in simulations that the lack of spectral fidelity relative to the physical system due to the existence of aliased spatial modes is the major cause of the FGI in the PIC model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Inadequacy of current approaches for characterizing membrane transport properties at high salinities

Cost optimal design of osmotic membrane processes requires an accurate estimate of membrane transport parameters across their full operational range. However, standard approaches for estimating these parameters rely on empirical methods, the accuracy of which remains unquantified as a function of temperature, salinity, and measurement error. Herein, we present a systematic accuracy analysis of previously developed methods for estimation of membrane transport properties in reverse osmosis, high-pressure reverse osmosis, forward osmosis, pressure retarded reverse osmosis, and osmotically assisted reverse osmosis. We use a Monte Carlo approach to sample the full range of feasible membrane water permeabilities, salt permeabilities, structural parameters, and operating conditions for these processes. These material and process parameters are then incorporated into a physical transport model for each process. Our analysis shows that the statistical uncertainty of current empirical methods for estimating membrane parameters increases by 5 times from low-salinity to high-salinity conditions. Here, the result of this work demonstrates that empirical methods are inadequate for precisely estimating membrane transport properties at high salinity and highlight a critical need for the development of statistically validated higher accuracy methods.

42 ENGINEERING↗
DOE OSTI2025

Disentangling the Impacts of Microtopography and Shrub Distribution on Snow Depth in a Subarctic Watershed: Toward a Predictive Understanding of Snow Spatial Variability

Snow plays a critical role in carbon cycling, vegetation dynamics, and permafrost hydrology at high latitudes by influencing surface energy exchange. Predicting snow distribution patterns is essential for understanding the evolution of Arctic ecosystems, yet scaling process-level knowledge to landscape predictions remains challenging. Here, we analyze snow depth (2019 and 2022), terrain elevation, and vegetation height from a watershed on the Seward Peninsula, Alaska, to examine how topography and shrubs shape snow redistribution across spatial scales. We find that snow depth is strongly coupled to terrain at scales below ∼60 m but becomes increasingly decoupled at larger scales. The topographic model of snow depth variation, which transforms terrain data to align with these scale-dependent snow patterns, is well correlated with local snow depth variations (linear fit R 2 > 0.5 for 85% of 100-m patches). A machine learning reconstruction of shrub canopy snow trapping reveals a simple exponential relationship between canopy structure and snow accumulation ( R 2 = 0.59), highlighting the combined influence of topography and vegetation on snow distribution. Together, these empirical relationships capture much of the observed snow variability in the watershed ( R 2 = 0.49, root mean square error (RMSE) = 30 cm), though systematic limitations persist in areas of strong scour and at coarser scales where wind-terrain interactions are more complex. These findings provide a framework for more efficient snow depth prediction and offer insights to improve snow-vegetation feedback representation in Earth System Models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

An experimental limit on the charge of antihydrogen

The properties of antihydrogen are expected to be identical to those of hydrogen, and any differences would constitute a profound challenge to the fundamental theories of physics. The most commonly discussed antiatom-based tests of these theories are searches for antihydrogen-hydrogen spectral differences (tests of CPT (charge-parity-time) invariance) or gravitational differences (tests of the weak equivalence principle). Here we, the ALPHA Collaboration, report a different and somewhat unusual test of CPT and of quantum anomaly cancellation. A retrospective analysis of the influence of electric fields on antihydrogen atoms released from the ALPHA trap finds a mean axial deflection of 4.1±3.4 mm for an average axial electric field of 0.51 V mm –1 . Combined with extensive numerical modelling, this measurement leads to a bound on the charge Qe of antihydrogen of Q=(–1.3±1.1±0.4) × 10 –8 . Here, e is the unit charge, and the errors are from statistics and systematic effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Equilibrium reconstruction based on core magnetic measurement and its applications on equilibrium transition in Joint-TEXT tokamak

Evaluation and reconstruction of plasma equilibrium, especially to resolve the safety factor profile, is imperative for advanced tokamak operation and physics study. Based on core magnetic measurement by the high resolution laser polarimeter-interferometer system (POLARIS), the equilibrium of Joint-TEXT (J-TEXT) plasma is reconstructed and profiles of safety factor, current density, and electron density are, therefore, obtained with high accuracy and temporal resolution. The equilibrium reconstruction procedure determines the equilibrium flux surfaces essentially from the data of POLARIS. Refraction of laser probe beam, a major error source of the reconstruction, has been considered and corrected, which leads to improvement of accuracy more than 10%. The error of reconstruction has been systematically assessed with consideration of realistic diagnostic performance and scrape-off layer region of plasma, and its accuracy has been verified. Fast equilibrium transitions both within a single sawtooth cycle and during the penetration of resonant magnetic perturbation have been investigated.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Active helium target: Neutron scalar polarizability extraction via Compton scattering

Precise measurement of the neutron scalar polarizabilities has been a lasting challenge because of the lack of a free-neutron target. Led by the University of Glasgow and the Mount Allison University groups of the A2 collaboration in Mainz, Germany, preparations have begun to test a recent theoretical model with an active helium target with the hope of determining these elusive quantities with small statistical, systematic, and model-dependent errors. Apparatus testing and background-event simulations have been carried out, with the full experiment projected to run in 2015. Once determined, these values can be applied to help understand quantum chromodynamics in the nonperturbative region.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Designing ITER motional Stark effect line shift (MSE-LS) spectrometers

As a part of ITER beam aided diagnostics, the design of Motional Stark Effect (MSE) diagnostic observing the emission from the Balmer-α line is underway. The physics of Stark splitting shows that the Stark manifold is polarization dependent, and the energy splitting results in a line shift proportional to the electric field. Due to the challenges of maintaining the calibration of the plasma facing mirrors in ITER, the conventional MSE polarimetry measurement technique is replaced with a spectral approach that is deemed more favorable in the ITER environment. The MSE line shift (LS) diagnostic is designed to quantify the Lorentz electric field magnitude by measuring the Stark manifold using visible spectroscopy. In the presence of large magnetic fields and high energy heating beams of 1 MeV, the expected Stark splitting is much larger than in typical devices. The MSE-LS design has unique challenges requiring careful consideration and modeling of its viewing geometry and photon budget. The MSE-LS approach on ITER is promising but has stringent demands on the allowable errors for the statistical and systematic fitting uncertainties. In this study, a full system model and numerical simulations of data for each sightline are completed. For a range of optical transmission fractions, photon noise analysis is conducted to determine the statistical uncertainties. This provides guidance on the spectrometer throughput, dispersion at the detector, optics, and other design choices. A conceptual design of a high throughput spectrometer with a volume phase transmission grating is presented.

Instruments & Instrumentation↗
DOE OSTI2012

COSMOLOGICAL CONSTRAINTS FROM GALAXY CLUSTERING AND THE MASS-TO-NUMBER RATIO OF GALAXY CLUSTERS

We place constraints on the average density ({Omega}{sub m}) and clustering amplitude ({sigma}{sub 8}) of matter using a combination of two measurements from the Sloan Digital Sky Survey: the galaxy two-point correlation function, w{sub p} (r{sub p} ), and the mass-to-galaxy-number ratio within galaxy clusters, M/N, analogous to cluster M/L ratios. Our w{sub p} (r{sub p} ) measurements are obtained from DR7 while the sample of clusters is the maxBCG sample, with cluster masses derived from weak gravitational lensing. We construct nonlinear galaxy bias models using the Halo Occupation Distribution (HOD) to fit both w{sub p} (r{sub p} ) and M/N for different cosmological parameters. HOD models that match the same two-point clustering predict different numbers of galaxies in massive halos when {Omega}{sub m} or {sigma}{sub 8} is varied, thereby breaking the degeneracy between cosmology and bias. We demonstrate that this technique yields constraints that are consistent and competitive with current results from cluster abundance studies, without the use of abundance information. Using w{sub p} (r{sub p} ) and M/N alone, we find {Omega}{sup 0.5}{sub m}{sigma}{sub 8} = 0.465 {+-} 0.026, with individual constraints of {Omega}{sub m} = 0.29 {+-} 0.03 and {sigma}{sub 8} = 0.85 {+-} 0.06. Combined with current cosmic microwave background data, these constraints are {Omega}{sub m} = 0.290 {+-} 0.016 and {sigma}{sub 8} = 0.826 {+-} 0.020. All errors are 1{sigma}. The systematic uncertainties that the M/N technique are most sensitive to are the amplitude of the bias function of dark matter halos and the possibility of redshift evolution between the SDSS Main sample and the maxBCG cluster sample. Our derived constraints are insensitive to the current level of uncertainties in the halo mass function and in the mass-richness relation of clusters and its scatter, making the M/N technique complementary to cluster abundances as a method for constraining cosmology with future galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

ON THE MASS-LOSS RATE OF MASSIVE STARS IN THE LOW-METALLICITY GALAXIES IC 1613, WLM, AND NGC 3109

We present a spectroscopic analysis of Very Large Telescope/X-Shooter observations of six O-type stars in the low-metallicity (Z {approx} 1/7 Z {sub sun}) galaxies IC 1613, WLM, and NGC 3109. The stellar and wind parameters of these sources allow us, for the first time, to probe the mass loss versus metallicity dependence of stellar winds at metallicities below that of the Small Magellanic Cloud (at Z {approx} 1/5 Z {sub sun}) by means of a modified wind momentum versus luminosity diagram. The wind strengths that we obtain for the objects in WLM and NGC 3109 are unexpectedly high and do not agree with theoretical predictions. The objects in IC 1613 tend toward a higher than expected mass-loss rate, but remain consistent with predictions within their error bars. We discuss potential systematic uncertainties in the mass-loss determinations to explain our results. However, if further scrutinization of these findings point towards an intrinsic cause for this unexpected sub-SMC mass-loss behavior, implications would include a higher than anticipated number of Wolf-Rayet stars and Ib/Ic supernovae in low-metallicity environments, but a reduced number of long-duration gamma-ray bursts produced through a single-star evolutionary channel.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey quasar sample: testing observational systematics on the Baryon Acoustic Oscillation measurement

ABSTRACT Baryon Acoustic Oscillations (BAO) are considered to be a very robust standard ruler against various systematics. This premise has been tested against observational systematics, but not to the level required for the next generation of galaxy surveys such as the Dark Energy Spectroscopic Instrument (DESI) and Euclid. In this paper, we investigate the effect of observational systematics on the BAO measurement of the final sample of quasars from the extended Baryon Oscillation Spectroscopic Survey Data Release 16 in order to prepare and hone a similar analysis for upcoming surveys. We employ catalogues with various treatments of imaging systematic effects using linear and neural network-based non-linear approaches and consider how the BAO measurement changes. We also test how the variations to the BAO fitting model respond to the observational systematics. As expected, we confirm that the BAO measurements obtained from the DR16 quasar sample are robust against imaging systematics well within the statistical error, while reporting slightly modified constraints that shift the line-of-sight BAO signal by less than 1.1 per cent. We use realistic simulations with similar redshift and angular distributions as the DR16 sample to conduct statistical tests for validating the pipeline, quantifying the significance of differences, and estimating the expected bias on the BAO scale in future high-precision data sets. Although we find a marginal impact for the eBOSS QSO data, the work presented here is of vital importance for constraining the nature of dark energy with the BAO feature in the new era of big data cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Consistent lensing and clustering in a low- S 8 Universe with BOSS, DES Year 3, HSC Year 1, and KiDS-1000

ABSTRACT We evaluate the consistency between lensing and clustering based on measurements from Baryon Oscillation Spectroscopic Survey combined with galaxy–galaxy lensing from Dark Energy Survey (DES) Year 3, Hyper Suprime-Cam Subaru Strategic Program (HSC) Year 1, and Kilo-Degree Survey (KiDS)-1000. We find good agreement between these lensing data sets. We model the observations using the Dark Emulator and fit the data at two fixed cosmologies: Planck (S8 = 0.83), and a Lensing cosmology (S8 = 0.76). For a joint analysis limited to large scales, we find that both cosmologies provide an acceptable fit to the data. Full utilization of the higher signal-to-noise small-scale measurements is hindered by uncertainty in the impact of baryon feedback and assembly bias, which we account for with a reasoned theoretical error budget. We incorporate a systematic inconsistency parameter for each redshift bin, A, that decouples the lensing and clustering. With a wide range of scales, we find different results for the consistency between the two cosmologies. Limiting the analysis to the bins for which the impact of the lens sample selection is expected to be minimal, for the Lensing cosmology, the measurements are consistent with A = 1; A = 0.91 ± 0.04 (A = 0.97 ± 0.06) using DES+KiDS (HSC). For the Planck case, we find a discrepancy: A = 0.79 ± 0.03 (A = 0.84 ± 0.05) using DES+KiDS (HSC). We demonstrate that a kinematic Sunyaev–Zeldovich-based estimate for baryonic effects alleviates some of the discrepancy in the Planck cosmology. This analysis demonstrates the statistical power of small-scale measurements; however, caution is still warranted given modelling uncertainties and foreground sample selection effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Analysis of DESI × DES using the Lagrangian effective theory of LSS

In this work we use Lagrangian perturbation theory to analyze the harmonic space galaxy clustering signal of the Bright Galaxy Survey (BGS) and luminous red galaxies (LRGs) targeted by the dark energy spectroscopic instrument (DESI), combined with the galaxy-galaxy lensing signal measured around these galaxies using Dark Energy Survey Year 3 source galaxies. The BGS and LRG galaxies are extremely well characterized by DESI spectroscopy and, as a result, lens galaxy redshift uncertainty and photometric systematics contribute negligibly to the error budget of our “2 × 2-point” analysis. On the modeling side, this work represents the first application of the SPINOSAURUS code, implementing an effective field theory model for galaxy intrinsic alignments, and we additionally introduce a new scheme ( MAIAR ) for marginalizing over the large uncertainties in the redshift evolution of the intrinsic alignment signal. Furthermore, this is the first application of a hybrid effective field theory model for galaxy bias based on the A emulus 𝜈 simulations. Our main result is a measurement of the amplitude of the lensing signal, 𝑆 8 = 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 = 0.85⁢0$^{+0.042}_{−0.050}$, consistent with values of this parameter derived from the primary cosmic microwave background. This constraint is artificially improved by a factor of 51% if we assume a more standard, but restrictive parametrization for the redshift evolution and sample dependence of the intrinsic alignment signal, and 63% if we additionally assume the nonlinear alignment model. Furthermore, we show that when fixing the cosmological model to the best-fit values from Planck PR4 there is > 5⁢𝜎 evidence for a deviation of the evolution of the intrinsic alignment signal from the functional form that is usually assumed in cosmic shear and galaxy-galaxy lensing studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

All-fiber interferometry: design and analysis

This paper presents a design procedure for a fiber interferometer, the optical system, and its associated electronic control. Analog and digital circuits were optimized to achieve an inexpensive and compact system. The lock-in amplifier required for phase control was designed using a field programmable gate array that was also configured to carry out the required phase stepping. The interferometer was built into two stages. The first stage used only one wavelength to measure samples with step heights in the hundreds of nanometers, with improvements in accuracy through averaging of data. The next stage included two measurement wavelengths to extend the measurement range to a few hundreds of micrometers. The paper conducts a systematic study of sources of errors. A least squares method that minimizes measurement time but retains accuracy has been used to measure height.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Evaluating a Transition-jump System for the Fermilab Main Injector Using Xsuite

We describe an Xsuite simulation framework for the Fermilab Main Injector (MI) along with an evaluation of transition-crossing behaviors in the accelerator. In particular, we studied the introduction of quadrupole magnets into the lattice as part of a transition-jump system that will be implemented through the $2^{nd}$ Proton Improvement Plan (PIP-II). Simulated beam losses spurred by transition-induced instabilities were assessed under several systematic effects, including MI quad errors, magnet-to-magnet variability in the jump magnets, power-supply errors, and timing jitter.

Schreckenberger, A. P. [Fermilab] (ORCID:000000016↗
DOE OSTI2017

Model Uncertainties for Valencia RPA Effect for MINERvA

This technical note describes the application of the Valencia RPA multi-nucleon effect and its uncertainty to QE reactions from the GENIE neutrino event generator. The analysis of MINERvA neutrino data in Rodrigues et al. PRL 116 071802 (2016) paper makes clear the need for an RPA suppression, especially at very low momentum and energy transfer. That published analysis does not constrain the magnitude of the effect; it only tests models with and without the effect against the data. Other MINERvA analyses need an expression of the model uncertainty in the RPA effect. A well-described uncertainty can be used for systematics for unfolding, for model errors in the analysis of non-QE samples, and as input for fitting exercises for model testing or constraining backgrounds. This prescription takes uncertainties on the parameters in the Valencia RPA model and adds a (not-as-tight) constraint from muon capture data. For MINERvA we apply it as a 2D ($q_0$,$q_3$) weight to GENIE events, in lieu of generating a full beyond-Fermi-gas quasielastic events. Because it is a weight, it can be applied to the generated and fully Geant4 simulated events used in analysis without a special GENIE sample. For some limited uses, it could be cast as a 1D $Q^2$ weight without much trouble. This procedure is a suitable starting point for NOvA and DUNE where the energy dependence is modest, but probably not adequate for T2K or MicroBooNE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗