Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error analysis”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Crosstalk in microwave SQUID multiplexers

Low-temperature detector technologies provide extraordinary sensitivity for applications ranging from precision measurements of the cosmic microwave background to high-resolution, high-rate x-ray, and γ-ray spectroscopy. To utilize this sensitivity, new instruments are being built, and new instruments are imagined, with ever greater pixel counts, but the scale of these instruments is limited by the capability of the readout electronics. Microwave SQUID multiplexing addresses the needs of these future instruments, exploiting gigahertz of bandwidths of coaxial cables and broadband components to combine hundreds to thousands of signals on a single readout line. A key feature of any multiplexer is the level of crosstalk between input channels. This crosstalk can degrade the sensitivity of the instrument, introduce systematic error, or simply confound data analysis. In this letter, we explain the primary mechanisms of crosstalk in a microwave SQUID multiplexer, calculate and measure their magnitude, and consider their effect and methods of mitigation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Recent Studies on the Asymptotic Convergence of the Spatial Discretization for Two-Dimensional Discrete Ordinates Solutions

In this paper, four types of quadrature schemes are used to define discrete directions in the solution of a two-dimensional fixed-source discrete ordinates problem in Cartesian geometry. Such schemes enable generating numerical results for averaged scalar fluxes over specified regions of the domain with high number (up to 105) of directions per octant. Two different nodal approaches, the ADO and AHOT-N0 methods, are utilized to obtain the numerical results of interest. The AHOT-N0 solutions on a sequence of refined meshes are then used to develop an asymptotic analysis of the spatial discretization error in order to derive a reference solution. It was more clearly observed that the spatial discretization error converges asymptotically with second order for the source region with all four quadratures employed, while for the other regions refined meshes along with tighter convergence criterion must be applied to evidence the same behavior. Although, in that case, some differences among the four quadrature schemes results were found.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Material identification in x-ray microscopy and micro CT using multi-layer, multi-color scintillation detectors

We demonstrate that a dual-layer, dual-color scintillator construct for microscopic CT, originally proposed to increase sensitivity in synchrotron imaging, can also be used to perform material quantification and classification when coupled with polychromatic illumination. We consider two different approaches to data handling: (1) a data-domain material decomposition whose estimation performance can be characterized by the Cramer-Rao lower bound formalism but which requires careful calibration and (2) an image-domain material classification approach that is more robust to calibration errors. The data-domain analysis indicates that useful levels of SNR (>5) could be achieved in one second or less at typical bending magnet fluxes for relatively large amounts of contrast (several mm path length, such as in a fluid flow experiment) and at typical undulator fluxes for small amount of contrast (tens of microns path length, such as an angiography experiment). The tools introduced could of course be used to study and optimize parameters for a wider range of potential applications. The image domain approach was analyzed in terms of its ability to distinguish different elemental stains by characterizing the angle between the lines traced out in a two-dimensional space of effective attenuation coefficient in the front and back layer images. As a result, this approach was implemented at a synchrotron and the results were consistent with simulation predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2015

New axion and hidden photon constraints from a solar data global fit

We present a new statistical analysis that combines helioseismology (sound speed, surface helium and convective radius) and solar neutrino observations (the {sup 8}B and {sup 7}Be fluxes) to place upper limits to the properties of non standard weakly interacting particles. Our analysis includes theoretical and observational errors, accounts for tensions between input parameters of solar models and can be easily extended to include other observational constraints. We present two applications to test the method: the well studied case of axions and axion-like particles and the more novel case of low mass hidden photons. For axions we obtain an upper limit at 3σ for the axion-photon coupling constant of g{sub aγ} < 4.1 · 10{sup −10} GeV{sup −1}. For hidden photons we obtain the most restrictive upper limit available accross a wide range of masses for the product of the kinetic mixing and mass of χ m < 1.8 ⋅ 10{sup −12} eV at 3σ. Both cases improve the previous solar constraints based on the Standard Solar Models showing the power of using a global statistical approach.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Limiting masses and radii of neutron stars and their implications

We combine the equation of state of dense matter up to twice nuclear saturation density 𝑛 sat obtained using chiral effective field theory (𝜒⁢EFT) and recent observations of neutron stars to gain insights about the high-density matter encountered in their cores. A key element in our study is the recent Bayesian analysis of correlated EFT truncation errors based on order-by-order calculations up to next-to-next-to-next-to-leading order in the 𝜒⁢EFT expansion. We refine the bounds on the maximum mass imposed by causality at high densities and provide stringent limits on the maximum and minimum radii of ∼1.4M ⊙ and ∼2.0M ⊙ stars. Including 𝜒⁢EFT predictions from 𝑛sat to 2𝑛 sat reduces the permitted ranges of the radius of a ~1.4M ⊙ star, 𝑅 1.4 , by ∼3.5km. If observations indicate 𝑅 1.4 < 11.2km, then our study implies that either the squared speed of sound 𝑐$^{2}_{𝑠}$ > 1/2 for densities above 2𝑛 sat or that 𝜒⁢EFT breaks down below 2𝑛 sat . Here, we also comment on the nature of the secondary compact object in GW190814 with mass ≃ 2.6M ⊙ and discuss the implications of massive neutron stars > 2.1⁢M ⊙ ⁡(2.6M ⊙ ) in future radio and gravitational-wave searches. Some form of strongly interacting matter with 𝑐$^{2}_{𝑠}$ > 0.35⁢(0.55) must be realized in the cores of such massive neutron stars. In the absence of phase transitions below 2𝑛 sat , the small tidal deformability inferred from GW170817 lends support for the relatively small pressure predicted by 𝜒⁢EFT for the baryon density 𝑛 B in the range 1–2𝑛 sat . Together they imply that the rapid stiffening required to support a high maximum mass should occur only when 𝑛 B ≳ 1.5–1.8 𝑛 sat .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

3D-printed sheet jet for stable megahertz liquid sample delivery at X-ray free-electron lasers

X-ray free-electron lasers (XFELs) can probe chemical and biological reactions as they unfold with unprecedented spatial and temporal resolution. A principal challenge in this pursuit involves the delivery of samples to the X-ray interaction point in such a way that produces data of the highest possible quality and with maximal efficiency. This is hampered by intrinsic constraints posed by the light source and operation within a beamline environment. For liquid samples, the solution typically involves some form of high-speed liquid jet, capable of keeping up with the rate of X-ray pulses. However, conventional jets are not ideal because of radiation-induced explosions of the jet, as well as their cylindrical geometry combined with the X-ray pointing instability of many beamlines which causes the interaction volume to differ for every pulse. This complicates data analysis and contributes to measurement errors. An alternative geometry is a liquid sheet jet which, with its constant thickness over large areas, eliminates the problems related to X-ray pointing. Since liquid sheets can be made very thin, the radiation-induced explosion is reduced, boosting their stability. These are especially attractive for experiments which benefit from small interaction volumes such as fluctuation X-ray scattering and several types of spectroscopy. Although their use has increased for soft X-ray applications in recent years, there has not yet been wide-scale adoption at XFELs. Here, gas-accelerated liquid sheet jet sample injection is demonstrated at the European XFEL SPB/SFX nano focus beamline. Its performance relative to a conventional liquid jet is evaluated and superior performance across several key factors has been found. This includes a thickness profile ranging from hundreds of nanometres to 60 nm, a fourfold increase in background stability and favorable radiation-induced explosion dynamics at high repetition rates up to 1.13 MHz. Its minute thickness also suggests that ultrafast single-particle solution scattering is a possibility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

A Bloom Filter Based Scalable Data Integrity Check Tool for Large-Scale Dataset

Large scale HPC applications are becoming increasingly data intensive. At Oak Ridge Leadership Computing Facility (OLCF), we are observing the number of files curated under individual project are reaching as high as 200 millions and project data size is exceeding petabytes. These simulation datasets, once validated, often needs to be transferred to archival system for long term storage or shared with the rest of the research community. Ensuring the data integrity of the full dataset at this scale is paramount important but also a daunting task. This is especially true considering that most conventional tools are serial and file-based, unwieldy to use and/or can't scale to meet user's demand.To tackle this particular challenge, this paper presents the design, implementation and evaluation of a scalable parallel checksumming tool, fsum, which we developed at OLCF. It is built upon the principle of parallel tree walk and work-stealing pattern to maximize parallelism and is capable of generating a single, consistent signature for the entire dataset at extreme scale. We also applied a novel bloom-filter based technique in aggregating signatures to overcome the signature ordering requirement. Given the probabilistic nature of bloom filter, we provided a detailed error and trade-off analysis. Using multiple datasets from production environment, we demonstrated that our tool can efficiently handle both very large files as well as many small-file based datasets. Our preliminary test showed that on the same hardware, it outperforms conventional tool by as much as 4×. It also exhibited near-linear scaling properties when provisioned with more compute resources.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Verification assessment of piston boundary conditions for Lagrangian simulation of compressible flow similarity solutions

This work is concerned with the use of similarity solutions of the compressible flow equations as benchmarks or verification test problems for finite-volume compressible flow simulation software. In practice, this effort can be complicated by the infinite spatial/temporal extent of many candidate solutions or “test problems.” Methods can be devised with the intention of ameliorating this inconsistency with the finite nature of computational simulation; the exact strategy will depend on the code and problem archetypes under investigation. For example, self-similar shock wave propagation can be represented in Lagrangian compressible flow simulations as rigid boundary-driven flow, even if no such “piston” is present in the counterpart mathematical similarity solution. The purpose of this work is to investigate in detail the methodology of representing self-similar shock wave propagation as a piston-driven flow in the context of various test problems featuring simple closed-form solutions of infinite spatial/temporal extent. The closed-form solutions allow for the derivation of similarly closed-form piston boundary conditions (BCs) for use in Lagrangian compressible flow solvers. Finally, the consequences of utilizing these BCs (as opposed to directly initializing the self-similar solution in a computational spatial grid) are investigated in terms of common code verification analysis metrics (e.g., shock strength/position errors and global convergence rates).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

SU-E-T-52: A New Device for Quality Assurance of a Single Isocenter Technique for the Simultaneous Treatment of Multiple Brain Metastases

Purpose: Comprehensive quality assurance (QA) of a single isocenter technique for the simultaneous treatment of multiple brain metastases is presently impractical due to the time consuming nature of measuring each lesion’s dose on film or with a micro-chamber. Three dimensional diode array and full field film measurements are sometimes used to evaluate these plans, but gamma analysis may not reveal local errors that have significant effects on one or a few of several targets. This work aimed to design, build and test a phantom to simplify comprehensive measurement and evaluation. Methods: A phantom was designed with 28 stackable slabs. The top and bottom slabs are 1.5 centimeters (cm) in thickness, and central 26 slabs are 0.5 cm thick. When assembled with radiochromic film in all 27 gaps, the phantom measures 16.5 x 15 x 19 cm. Etchings were designed to aide in identification of specific film planes on computed tomography (CT) images and correlation of individual PTVs with closest bisecting planes. Patient verification plans with a total of 16 PTVs were calculated on the phantom CT, and test deliveries both with and without couch kicks were performed to test the ability to identify correct film placements and subsequent PTV specific dose distributions on the films. Results: Bisecting planes corresponding to PTV locations were easily identified, and PTV specific dose distributions were clear for all 16 targets. For deliveries with couch kicks, the phantom PTV dose distributions closely approximated those calculated on the patient’s CT. For deliveries without couch kicks, PTV specific dosimetry was also possible, although the distributions had ‘ghosts’ equaling the number of couch kicks, with distance between ghosts increasing with distance from the isocenter. Conclusion: A new phantom facilitates fast comprehensive commissioning validation and PTV specific dosimetry for a single isocenter technique for treating multiple brain metastases. This work was partially funded by CIRS, Inc.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Error Bounds for Dynamical Spectral Estimation

Dynamical spectral estimation is a well-established numerical approach for estimating eigenvalues and eigenfunctions of the Markov transition operator from trajectory data. Although the approach has been widely applied in biomolecular simulations, its error properties remain poorly understood. Here we analyze the error of a dynamical spectral estimation method called “the variational approach to conformational dynamics" (VAC). We bound the approximation error and estimation error for VAC estimates. Our analysis establishes VAC's convergence properties and suggests new strategies for tuning VAC to improve accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Evaluating radiation impact on transmon qubits in above and underground facilities

Superconducting qubits can be sensitive to energy deposits caused by cosmic rays and ambient radioactivity. While previous studies have explored correlated effects in time and space due to cosmic ray interactions, we present the first direct comparison of a transmon qubit’s performance measured at two distinct sites: the above-ground SQMS facility (Fermilab, US) and the deep-underground Gran Sasso Laboratory (Italy). Despite the stark difference in radiation levels, we observe a similar average qubit relaxation time of approximately 80 microseconds at both locations. To investigate radiation-induced events, we employ a fast decay detection protocol, comparing the relative rates of events between the two environments. Although intrinsic noise remains the dominant source of errors in superconducting qubits, our analysis revealed a significant excess of radiation-induced events for high-coherence transmon qubits operated above-ground. Finally, using γ-ray sources with increasing activity levels, we evaluate the qubit response in a controlled low-background environment.

De Dominicis, Francesco [GSSI, Aquila; INFN, LNS]↗
DOE OSTI2017

Managing the Performance/Error Tradeoff of Floating-point Intensive Applications

Modern embedded systems are becoming more reliant on real-valued arithmetic as they employ mathematically complex vision algorithms and sensor signal processing. Double-precision floating point is the most commonly used precision in computer vision algorithm implementations. A single-precision floating point can provide a performance boost due to less memory transfers, less cache occupancy, and relatively faster mathematical operations on some architectures. However, adopting it can result in loss of accuracy. Identifying which parts of the program can run in single-precision floating point with low impact on error is a manual and tedious process. In this paper, we propose an automatic approach to identify parts of the program that have a low impact on error using shadow-value analysis. Our approach provides the user with a performance/error tradeoff, using which the user can decide how much accuracy can be sacrificed in return for performance improvement. We illustrate the impact of the approach using a well known implementation of Apriltag detection used in robotics vision. In conclusion, we demonstrate that an average 1.3x speedup can be achieved with no impact on tag detection, and a 1.7x speedup with only 4% false negatives.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Saturn Facility Oil Transfer Automation System

The Saturn accelerator, owned by Sandia National Laboratories, has been in operation since the early 1980s and still has many of the original systems. A critical legacy system is the oil transfer system which transfers 250,000 gallons of transformer oil from outside storage tanks to the Saturn facility. The oil transfer system was identified for upgrade to current technology standards. Using the existing valves, pumps, and relay controls, the system was automated using the National Instruments cRIO FGPA platform. Engineered safety practices, including a failure mode effects analysis, were used to develop error handling requirements. The uniqueness of the Saturn Oil Automated Transfer System (SOATS) is in the graphical user interface. The SOATS uses an HTML interface to communicate to the cRIO, creating a platform independent control system. The SOATS was commissioned in April 2013.

02 PETROLEUM↗
DOE OSTI2025

Reassessment of Dose Contributions from the M-Area Glass Special Waste Form Buried in the E-Area Low-Level Waste Facility

During a recent revision to the groundwater (GW) and inadvertent human intruder (IHI) screening analysis report (Aleman and Hamm, 2023), weaknesses within the screening logic were observed resulting in U-235 being brought back into the list of radionuclides requiring inventory limits within the M-Area Glass SWF. In addition, while reviewing the earlier vadose zone (VZ) transport analysis efforts, an input deck error was observed in the U-234G VZ transport simulation.

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

Gate Set Tomography

Gate set tomography (GST) is a protocol for detailed, predictive characterization of logic operations (gates) on quantum computing processors. Early versions of GST emerged around 2012-13, and since then it has been refined, demonstrated, and used in a large number of experiments. This paper presents the foundations of GST in comprehensive detail. The most important feature of GST, compared to older state and process tomography protocols, is that it is calibration-free . GST does not rely on pre-calibrated state preparations and measurements. Instead, it characterizes all the operations in a gate set simultaneously and self-consistently, relative to each other. Long sequence GST can estimate gates with very high precision and efficiency, achieving Heisenberg scaling in regimes of practical interest. In this paper, we cover GST's intellectual history, the techniques and experiments used to achieve its intended purpose, data analysis, gauge freedom and fixing, error bars, and the interpretation of gauge-fixed estimates of gate sets. Our focus is fundamental mathematical aspects of GST, rather than implementation details, but we touch on some of the foundational algorithmic tricks used in the pyGSTi implementation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Mitigation of LEO Satellite Brightness and Trail Effects on the Rubin Observatory LSST

In this work, we report studies on the mitigation of optical effects of bright low-Earth-orbit (LEO) satellites on Vera C. Rubin Observatory and its Legacy Survey of Space and Time (LSST). These include options for pointing the telescope to avoid satellites, laboratory investigations of bright trails on the Rubin Observatory LSST camera sensors, algorithms for correcting image artifacts caused by bright trails, experiments on darkening SpaceX Starlink satellites, and ground-based follow-up observations. The original Starlink v0.9 satellites are g ~ 4.5 mag, and the initial experiment "DarkSat" is g ~ 6.1 mag. Future Starlink darkening plans may reach g ~ 7 mag, a brightness level that enables nonlinear image artifact correction to well below background noise. However, the satellite trails will still exist at a signal-to-noise ratio ~ 100, generating systematic errors that may impact data analysis and limit some science. For the Rubin Observatory 8.4 m mirror and a satellite at 550 km, the full width at half maximum of the trail is about 3'' as the result of an out-of-focus effect, which helps avoid saturation by decreasing the peak surface brightness of the trail. For 48,000 LEOsats of apparent magnitude 4.5, about 1% of pixels in LSST nautical twilight images would need to be masked.

79 ASTRONOMY AND ASTROPHYSICS↗