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,387 records · Page 77

Modeling optical systematics for the Taurus CMB experiment

We simulate a variety of optical systematics for Taurus, a balloon-borne cosmic microwave background (CMB) polarisation experiment, to assess their impact on large-scale E-mode polarisation measurements and constraints of the optical depth to reionisation τ. We model a one-month flight of Taurus from Wanaka, New Zealand aboard a super-pressure balloon (SPB). We simulate night-time scans of both the CMB and dust foregrounds in the 150 GHz band, one of Taurus's four observing bands. We consider a variety of possible systematics that may affect Taurus's observations, including non-gaussian beams, pointing reconstruction error, and half-wave plate (HWP) non-idealities. For each of these, we evaluate the residual power in the difference between maps simulated with and without the systematic, and compare this to the expected signal level corresponding to Taurus's science goals. Our results indicate that most of the HWP-related systematics can be mitigated to be smaller than sample variance by calibrating with Planck's TT spectrum and using an achromatic HWP model, with a preference for five layers of sapphire to ensure good systematic control. However, additional beam characterization will be required to mitigate far-sidelobe pickup from dust on larger scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Characterization of the Impact of Soft Errors on Iterative Methods

Soft errors caused by transient bit flips have the potential to significantly impact an application’s behavior. This has motivated the design of an array of techniques to detect, isolate, and correct soft errors using microarchitectural, architectural, compilation-based, or application-level techniques to minimize their impact on the executing application. The first step towards the design of good error detection/correction techniques involves an understanding of an application’s vulnerability to soft errors. In this paper, we present the first comprehensive characterization of the impact of soft errors on the convergence characteristics of six iterative methods using application-level fault injection. In particular, we consider the use of iterative methods to incrementally solve a linear systems of equations, which constitutes the core kernel in many scientific applications. We analyze the impact of soft errors in terms of the type of error (single- vs multi-bit), the distribution and location of bits affected, the data structure and the statement impacted, and variation with time. In addition to understanding the vulnerability of iterative solvers to soft errors, this characterization can aid the design of fault injection campaigns that ensure systematic coverage.

fault injection, iterative methods, Silent data co↗

SPECTROSCOPIC AND INTERFEROMETRIC MEASUREMENTS OF NINE K GIANT STARS

We present spectroscopic and interferometric measurements for a sample of nine K giant stars. These targets are of particular interest because they are slated for stellar oscillation observations. Our improved parameters will directly translate into reduced errors in the final masses for these stars when interferometric radii and asteroseismic densities are combined. Here, we determine each star’s limb-darkened angular diameter, physical radius, luminosity, bolometric flux, effective temperature, surface gravity, metallicity, and mass. When we compare our interferometric and spectroscopic results, we find no systematic offsets in the diameters and the values generally agree within the errors. Our interferometric temperatures for seven of the nine stars are hotter than those determined from spectroscopy with an average difference of about 380 K.

79 ASTRONOMY AND ASTROPHYSICS↗

Application of step-response lambda tuning to proportional-integral controllers in water resource recovery facilities

Proportional-integral-derivative (PID) controllers in water resource recovery facilities (WRRFs) feedback control loops are commonplace. While simple to implement, such control loops are rarely tuned optimally or systematically. Heuristic tuning approaches are commonly applied with varying degrees of success using trial-and-error, ad hoc tuning rules, or duplication of tuning values from a similar system. However, there are effective methods, such as lambda tuning, produce acceptable tuning with limited effort. These are based on the step-response method, where a manual process perturbation is used to define the relationship between the manipulated and controlled variables. Based on such an experiment, a simple process model is constructed and used to determine the controller tuning values. In this work, we used the step-response method and lambda tuning for two control systems in full-scale WRRFs. This led to responsive and stable behavior of the controlled system as defined by the absolute average error of the controlled variable to setpoint and standard deviation of the manipulated variable. Tuning of feedback control loops can be completed successfully through a systematic approach, and this work suggests that tuning tools, like lambda, should be part of all wastewater treatment control engineers' toolbox.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Direct Numerical Simulation of Pore-Scale Flow in a Bead Pack: Comparison with Magnetic Resonance Imaging Observations

A significant body of current research is aimed at developing methods for numerical simulation of flow and transport in porous media that explicitly resolve complex pore and solid geometries, and at utilizing such models to study the relationships between fundamental pore-scale processes and macroscopic manifestations at larger (i.e., Darcy) scales. A number of different numerical methods for pore-scale simulation have been developed, and have been extensively tested and validated for simplified geometries. However, validation of pore-scale simulations of fluid velocity for complex, three-dimensional (3D) pore geometries that are representative of natural porous media is challenging due to our limited ability to measure pore-scale velocity in such systems. Recent advances in magnetic resonance imaging (MRI) offer the opportunity to measure not only the pore geometry, but also local fluid velocities under steady-state flow conditions in 3D and with high spatial resolution. In this paper, we present a 3D velocity field measured at sub-pore resolution (tens of micrometers) over a centimeter-scale 3D domain using MRI methods. We have utilized the measured pore geometry to perform 3D simulations of Navier-Stokes flow over the same domain using direct numerical simulation techniques. We present a comparison of the numerical simulation results with the measured velocity field. It is shown that the numerical results match the observed velocity patterns well overall except for a variance and small systematic scaling which can be attributed to the known experimental error in the MRI measurements. The comparisons presented here provide strong validation of the pore-scale simulation methods and new insights for interpretation of uncertainty in MRI measurements of pore-scale velocity. This study also provides a potential benchmark for future comparison of other pore-scale simulation methods.

computational fluid dynamics↗

Parameterization Strategies for Intermolecular Potentials for Predicting Trajectory-Based Collision Parameters

The accuracy of separable strategies for constructing full-dimensional potential energy surfaces for collisional energy transfer and collision rate calculations is studied systematically for three alcohols (A = methanol, ethanol, and butanol) and three bath gases (M = Ar, N 2 , and H 2 O). The fitting efficiency (defined as the number of ab initio data required to achieve parametrizations of a desired accuracy) is quantified for both pairwise (Buckingham or "exp6") and nonpairwise (permutationally invariant polynomials, PIPs, of Morse variables) functional forms, for four sampling strategies, and as a function of the complexity and anisotropy of the interaction potential. We find that convergence with respect to the number of sampled ab initio data is largely independent of the choice of functional form but instead varies nearly linearly with the number of adjustable parameters and depends strongly on the sampling strategy. Specifically, the use of biased Sobol quasirandom sampling is similar to 7X more efficient than using unbiased pseudorandom sampling, on average, requiring just similar to 3 computed ab initio energies per adjustable fitting parameter. The pairwise exp6 functional form is shown to provide accurate and transferable parametrizations for M = Ar but is unable to accurately describe alcohol interactions with M = N 2 and H 2 O. The nonpairwise PIP functional form, which is systematically improvable, can produce separable parametrizations with arbitrarily small fitting errors. However, these can suffer from overfitting, which is demonstrated using dynamics calculations of collision parameters for a large number of exp6 and PIP parametrizations. The tests described here validate a robust strategy for automatically generating A + M potential energy surfaces with minimal human intervention, including a quantifiable out of sample metric for judging the accuracy of the fitted surface. We further analyze this set of automatically generated potential energy surfaces to identify areas where more sophisticated fitting strategies may be desired, including pruning of the PIP expansions for large systems and improved sampling strategies more closely coupled with the description of the functional form.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimization of foreground moment deprojection for semi-blind CMB polarization reconstruction

Abstract Upcoming Cosmic Microwave Background (CMB) experiments, aimed at measuring primordial CMB polarization B-modes, require exquisite control of instrumental systematics and Galactic foreground contamination. Blind minimum-variance techniques, like the Needlet Internal Linear Combination (NILC), have proven effective in reconstructing the CMB polarization signal and mitigating foregrounds and systematics across diverse sky models without suffering from foreground mismodelling errors. Still, residual foreground contamination from NILC may bias the recovered CMB polarization at large angular scales when confronted with the most complex foreground scenarios.By adding constraints to NILC to deproject statistical moments of the Galactic emission, the Constrained Moment ILC (cMILC) method has been demonstrated to further enhance foreground subtraction, albeit with an associated increase in overall noise variance. Faced with this trade-off between foreground bias reduction and overall variance minimization, there is still no recipe on which moments to deproject and which are better suited for blind variance minimization. To address this, we introduce the optimized cMILC (ocMILC) pipeline, which performs full automated optimization of the required number and set of foreground moments to deproject, pivot parameter values, and deprojection coefficients across the sky and angular scales, depending on the actual sky complexity, available frequency coverage, and experiment sensitivity. The optimal number of moments for deprojection, before paying significant noise penalty, is determined through a data diagnosis inspired by the Generalized NILC (GNILC) method.Validated on B-mode simulations of thePICOspace mission concept with four challenging foreground models, ocMILC exhibits lower Galactic foreground contamination compared to NILC and cMILC at all angular scales, with limited noise penalty. This multi-layer optimization enables the ocMILC pipeline to achieve unbiased posteriors of the tensor-to-scalar ratio, regardless of foreground complexity.

Astronomy & Astrophysics↗

Probing the galaxy–halo connection with total satellite luminosity

ABSTRACT We demonstrate how the total luminosity in satellite galaxies is a powerful probe of dark matter haloes around central galaxies. The method cross-correlates central galaxies in spectroscopic galaxy samples with fainter galaxies detected in photometric surveys. Using models, we show that the total galaxy luminosity, Lsat, scales linearly with host halo mass, making Lsat an excellent proxy for Mh. Lsat is also sensitive to the formation time of the halo. We demonstrate that probes of galaxy large-scale environment can break this degeneracy. Although this is an indirect probe of the halo, it yields a high signal-to-noise ratio measurement for galaxies expected to occupy haloes at <1012 M⊙, where other methods suffer from larger errors. In this paper, we focus on observational and theoretical systematics in the Lsat method. We test the robustness of our method of finding central galaxies and our methods of estimating the number of background galaxies. We implement this method on galaxies in the Sloan Digital Sky Survey (SDSS) data, with satellites identified in fainter imaging data. We find excellent agreement between our theoretical predictions and the observational measurements. Finally, we compare our Lsat measurements to weak lensing estimates of Mh for red and blue subsamples. In the stellar mass range where the measurements overlap, we find consistent results, where red galaxies live in larger haloes. However, the Lsat approach allows us to probe significantly lower mass galaxies. At these masses, the Lsat values are equivalent. This example shows the potential of Lsat as a probe of dark haloes.

79 ASTRONOMY AND ASTROPHYSICS↗

Programmer-guided reliability for extreme-scale applications

We present “programmer-guided reliability” (PGR) as a systematic conceptual approach to address the expected rise in soft errors in coming extreme-scale systems at the application level. The approach involves instrumentation of the application with code to detect data corruption errors. The location and nature of these error detectors are at the discretion of the programmer, who uses their knowledge and experience with the problem domain, the application, the solution algorithms, etc., to determine the most vulnerable areas of the code and the most appropriate ways to detect data corruption. To illustrate the approach, we provide examples of error detectors from four different benchmark-scale applications. We also describe a simple control framework that allows for runtime configuration of the error detectors without recompilation of the application, as well as dynamic reconfiguration during the execution of the application. Lastly, we discuss a number of future directions building on the basic PGR approach, including the incorporation of some general error detectors into the programming environment in order to make them more easily usable by the programmer.

97 MATHEMATICS AND COMPUTING↗

Lessons from Past Nuclear Disarmament: What Worked What Did Not

This report studies the efforts of the international community and the United States to verif)/ the denuclearization of various countries: Ukraine, South Africa, Iraq, Taiwan, and Libya. In doing so, it considers the verification of nuclear warhead destruction and the accounting of nuclear materials. Each case study contributes to the understanding we have of which verification procedures worked and which did not and what factors contributed to that success and which did not. The most important factor contributing to successful verification is the cooperation of the subject country. If a country has made the strategic decision to cooperate, then it is possible that verification can be successful. If the country chooses not to cooperate, verification might rest on random and unpredictable events. This unpredictability of verification is politically unacceptable. Even if verification is judged to be successful by the implementing agencies, outsiders can cast doubts on it by pointing out small and potentially unavoidable errors in material amounts verified. It is vitally important to systematically preserve forensic evidence, especially shipping and receiving records to avoid just such issues. Historically, such records are the most important evidence in verifying denuclearization. Technology is far less important than forensic analysis of records. It does play a supporting role in verifying some declarations but there can be an enormous delay in drawing conclusions caused by the necessity to analyze statistically large samples of materials. The confidence one gets from such analysis can be delayed well past the time it is politically significant. The most important advance in technology would be to accelerate that sample analysis process.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Vector Calibration for Vertical Tests of Superconducting RF Cavities

We propose procedure for vector calibration of Q 0 vs E acc curve measurement for superconducting resonators. The proposed vector calibration procedure and hardware has several advantages over the techniques currently used for vertical cavity tests. It provides high level of automation and reduces manpower and preparation time for the test; it provides in situ, real time control of hardware parameters and allows to detect hardware failure and operator's errors; it allows to control and exclude the most significant systematic uncertainties, measured values of Q 0 (E acc ) can be made traceable to national standards.

43 PARTICLE ACCELERATORS↗

Adaptive Determination of the Optimal Exchange Location in Wall-Modeled Large-Eddy Simulation

Wall-modeled large-eddy simulation introduces a modeling interface (or exchange location) separating the wall-modeled layer from the rest of the domain. The current state-of-the-art is to rely on user expertise when choosing where to place this modeling interface, whether this choice is tied to the grid or not. This paper presents a postprocessing algorithm that determines the exchange location systematically. The algorithm is based on a model for the error in the predicted wall shear stress and a model for the computational cost, and then finds the exchange location that minimizes a combination of the two. Here, the algorithm is tested both a priori and a posteriori using an equilibrium wall model for the flow over a wall-mounted hump, a boundary layer in an adverse pressure gradient, and a shock/boundary-layer interaction. The algorithm produces exchange locations that mostly agree with what an experienced user would suggest, with thinner wall-modeled layers in nonequilibrium flow regions and thicker wall-modeled layers where the boundary layer is closer to equilibrium. This suggests that the algorithm should be useful in simulations of realistic and highly complex geometries.

42 ENGINEERING↗

The Sloan Digital Sky Survey Reverberation Mapping Project: An investigation of biases in C iv emission line properties

We investigate the dependence on data quality of quasar properties measured from the C iv emission line region at high redshifts. Our measurements come from 32 epochs of Sloan Digital Sky Survey Reverberation Mapping Project spectroscopic observations of 482 z>1.46 quasars. We compare the differences between measurements made from the single-epoch (SE) and coadded spectra, focusing on the C iv λ1549 emission line because of its importance for studies of high-redshift quasar demographics and physical properties, including black hole masses. In addition to statistical errors increasing (by factors of ∼2–4), we find increasing systematic offsets with decreasing signal-to-noise ratio (S/N). The systematic difference (measurement uncertainty) in our lowest-S/N (<5) subsample between the SE and coadded spectrum (i) C iv equivalent width is 17 Å (31 Å), (ii) centroid wavelength is <1 Å (2 Å), and fractional velocity widths, ΔV/V, characterized by (iii) the line dispersion, σ {sub l}, is 0.104 (0.12), and (iv) the mean absolute deviation (MAD) is 0.072 (0.11). These remain smaller than the 1σ measurement uncertainties for all subsamples considered. The MAD is found to be the most robust line-width characterization. Offsets in the C iv FWHM velocity width and the C iv profile characterized by FWHM/σ {sub l} are only smaller than the statistical uncertainties when S/N > 10, although offsets in lower-S/N spectra exceed the statistical uncertainties by only a factor of ∼1.5 and may depend on the type of functional fit to the line. Characterizing the C iv line profile by the kurtosis is the least robust property investigated, as the median systematic coadded–SE measurement differences are larger than the statistical uncertainties for all S/N subsamples.

79 ASTRONOMY AND ASTROPHYSICS↗

Systematic study of the uncertainties in fitting the cosmic positron data by AMS-02

The operation of AMS-02 opens a new era for the study of cosmic ray physics with unprecedentedly precise data which are comparable with the laboratory measurements. The high precision data allow a quantitative study on the cosmic ray physics and give strict constraints on the nature of cosmic ray sources. However, the intrinsic errors from the theoretical models to interpret the data become dominant over the errors in the data. In the present work we try to give a systematic study on the uncertainties of the models to explain the AMS-02 positron fraction data, which shows the cosmic ray e{sup +}e{sup −} excesses together with the PAMELA and Fermi-LAT measurements. The excesses can be attributed to contributions from the extra e{sup +}e{sup −} sources, such as pulsars or the dark matter annihilation. The possible systematic uncertainties of the theoretical models considered include the cosmic ray propagation, the treatment of the low energy data, the solar modulation, the pp interaction models, the nuclei injection spectrum and so on. We find that in general a spectral hardening of the primary electron injection spectrum above ∼50–100 GeV is favored by the data. Furthermore, the present model uncertainties may lead to a factor of ∼2 enlargement in the determination of the parameter regions of the extra source, such as the dark matter mass, annihilation rate and so on.

79 ASTRONOMY AND ASTROPHYSICS↗

The impact of light polarization effects on weak lensing systematics

ABSTRACT A fraction of the light observed from edge-on disc galaxies is polarized due to two physical effects: selective extinction by dust grains aligned with the magnetic field and scattering of the anisotropic starlight field. Since the reflection and transmission coefficients of the reflecting and refracting surfaces in an optical system depend on the polarization of incoming rays, this optical polarization produces both (a) a selection bias in favour of galaxies with specific orientations and (b) a polarization-dependent point spread function (PSF). In this work, we build toy models to obtain for the first time an estimate for the impact of polarization on PSF shapes and the impact of the selection bias due to the polarization effect on the measurement of the ellipticity used in shear measurements. In particular, we are interested in determining if this effect will be significant for Wide-Field Infrared Survey Telescope (WFIRST). We show that the systematic uncertainties in the ellipticity components are 8 × 10−5 and 1.1 × 10−4 due to the selection bias and PSF errors respectively. Compared to the overall requirements on knowledge of the WFIRST PSF ellipticity (4.7 × 10−4 per component), both of these systematic uncertainties are sufficiently close to the WFIRST tolerance level that more detailed studies of the polarization effects or more stringent requirements on polarization-sensitive instrumentation for WFIRST are required.

79 ASTRONOMY AND ASTROPHYSICS↗

Interlaboratory Study of Digital Volume Correlation Error Due to X-Ray Computed Tomography Equipment and Scan Parameters: an Update from the DVC Challenge

Background: The quality of Digital Volume Correlation (DVC) full-field displacement measurements depends directly on the characteristics of the X-ray Computed Tomography (XCT) equipment, and scan procedures used to acquire the tomographic images. Objective: In this work, we seek to experimentally study the effects of XCT equipment and tomographic scan procedures on the quality of these images for DVC analysis, and to survey the level of DVC error that may be achieved using standard XCT operating procedures. Methods: Six participants in an interlaboratory study acquired high-quality XCT scans of a syntactic foam before and after rigid body motion. The resulting images were correlated using commercial DVC software to quantify error sources due to random image noise, reconstruction artifacts, as well as systematic spatial or temporal distortion. Results: In the absence of rigid body motion, the standard deviation of the displacement measurements ranged from 0.012 to 0.043 voxels using a moderate subvolume size, indicating that subvoxel measurement resolution could readily be achieved with a variety of XCT equipment and scan recipes. Comparison of consecutive scans without rigid body motion showed transient dilatational displacement gradients due to self-heating of the X-ray source and/or thermal expansion of the foam. Evaluation of the scans after rigid body motion showed significant, machine-specific spatial distortion in the displacement fields of up to 0.5 voxels; new approaches to remove this error need to be developed. Conclusions: Analysis of the scan protocols used in the interlaboratory study, as well as a complementary parametric sensitivity study, showed that the DVC error was strongly influenced by the XCT equipment, but could be mitigated by adjusting the total scan duration.

47 OTHER INSTRUMENTATION↗

Primer terminal ribonucleotide alters the active site dynamics of DNA polymerase η and reduces DNA synthesis fidelity

DNA polymerases catalyze DNA synthesis with high efficiency, which is essential for all life. Extensive kinetic and structural efforts have been executed in exploring mechanisms of DNA polymerases, surrounding their kinetic pathway, catalytic mechanisms, and factors that dictate polymerase fidelity. Recent time-resolved crystallography studies on DNA polymerase η (Pol η) and β have revealed essential transient events during the DNA synthesis reaction, such as mechanisms of primer deprotonation, separated roles of the three metal ions, and conformational changes that disfavor incorporation of the incorrect substrate. DNA-embedded ribonucleotides (rNs) are the most common lesion on DNA and a major threat to genome integrity. While kinetics of rN incorporation has been explored and structural studies have revealed that DNA polymerases have a steric gate that destabilizes ribonucleotide triphosphate binding, the mechanism of extension upon rN addition remains poorly characterized. Using steady-state kinetics, static and time-resolved X-ray crystallography with Pol η as a model system, we showed that the extra hydroxyl group on the primer terminus does alter the dynamics of the polymerase active site as well as the catalysis and fidelity of DNA synthesis. During rN extension, Pol η error incorporation efficiency increases significantly across different sequence contexts. Finally, our systematic structural studies suggest that the rN at the primer end improves primer alignment and reduces barriers in C2'-endo to C3'-endo sugar conformational change. Overall, our work provides further mechanistic insights into the effects of rN incorporation on DNA synthesis.

59 BASIC BIOLOGICAL SCIENCES↗

Permutationally Invariant Polynomial Expansions with Unrestricted Complexity

A general strategy is presented for constructing and validating permutationally invariant polynomial (PIP) expansions for chemical systems of any stoichiometry. Demonstrations are made for three categories of gas-phase dynamics and kinetics: collisional energy-transfer trajectories for predicting pressure-dependent kinetics, three-body collisions for describing transient van der Waals adducts relevant to atmospheric chemistry, and nonthermal reactivity via quasiclassical trajectories. In total, 30 systems are considered with up to 15 atoms and 39 degrees of freedom. Permutational invariance is enforced in PIP expansions with as many as 13 million terms and 13 permutationally distinct atom types by taking advantage of petascale computational resources. The quality of the PIP expansions is demonstrated through the systematic convergence of in-sample and out-of-sample errors with respect to both the number of training data and the order of the expansion, and these errors are shown to predict errors in the dynamics for both reactive and nonreactive applications. Here, the parallelized code distributed as part of this work enables the automation of PIP generation for complex systems with multiple channels and flexible user-defined symmetry constraints and for automatically removing unphysical unconnected terms from the basis set expansions, all of which are required for simulating complex reactive systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗