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,711 records · Page 95

Ab Initio Bulk Free Energy Surface of Proper Ferroelectrics

We report a systematic and accurate approach for deriving the bulk free energy surface (FES), a function of temperature, polarization, and strain, from the first-principles density functional theory (DFT) of proper ferroelectrics. The core of our approach is the metadynamics algorithm that extracts the polarization dependence of the FES from all-atom molecular dynamics simulations without an a priori ansatz. The rest of the FES is derived from the metadynamics trajectories that span the relevant phase space. We demonstrate our approach in the case of lead titanate. The errors across the phase transition, due to DFT numerics, all-atom molecular dynamics, and free energy evaluation by enhanced sampling, can be systematically controlled and are of the order of 1 meV/atom. The accuracy of the resulting ab initio FES is only limited by the adopted functional approximation of DFT.

Xie, Pinchen [Lawrence Berkeley National Laborator↗

Patient-based radiographic exposure factor selection: a systematic review

Digital technology has wider exposure latitude and post-processing algorithms which can mask the evidence of underexposure and overexposure. Underexposure produces noisy, grainy images which can impede diagnosis and overexposure results in a greater radiation dose to the patient. These exposure errors can result from inaccurate adjustment of exposure factors in response to changes in patient thickness. This study aims to identify all published radiographic exposure adaptation systems which have been, or are being, used in general radiography and discuss their applicability to digital systems. Studies in EMBASE, MEDLINE, CINAHL and SCOPUS were systematically reviewed. Some of the search terms used were exposure adaptation, exposure selection, exposure technique, 25% rule, 15% rule, DuPont™ Bit System and radiography. A manual journal-specific search was also conducted in The Radiographer and Radiologic Technology. Studies were included if they demonstrated a system of altering exposure factors to compensate for variations in patients for general radiography. Studies were excluded if they focused on finding optimal exposures for an ‘average’ patient or focused on the relationship between exposure factors and dose. The database search uncovered 11 articles and the journal-specific search uncovered 13 articles discussing systems of exposure adaptation. They can be categorised as simple one-step guidelines, comprehensive charts and computer programs. Only two papers assessed the efficacy of exposure adjustment systems. No literature compares the efficacy of exposure adaptations system for film/screen radiography with digital radiography technology nor is there literature on a digital specific exposure adaptation system.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Evaluating Cloud Feedback Components in Observations and Their Representation in Climate Models

This study quantifies the contribution of individual cloud feedbacks to the total short-term cloud feedback in satellite observations over the period 2002–2014 and evaluates how they are represented in climate models. The observed positive total cloud feedback is primarily due to positive high-cloud altitude, extratropical high- and low-cloud optical depth, and land cloud amount feedbacks partially offset by negative tropical marine low-cloud feedback. Seventeen models from the Atmosphere Model Intercomparison Project of the sixth Coupled Model Intercomparison Project are analyzed. The models generally reproduce the observed moderate positive short-term cloud feedback. However, compared to satellite estimates, the models are systematically high-biased in tropical marine low-cloud and land cloud amount feedbacks and systematically low-biased in high-cloud altitude and extratropical high- and low-cloud optical depth feedbacks. Errors in modeled short-term cloud feedback components identified in this analysis highlight the need for improvements in model simulations of the response of high clouds and tropical marine low clouds. Our results suggest that skill in simulating interannual cloud feedback components may not indicate skill in simulating long-term cloud feedback components.

58 GEOSCIENCES↗

Status of RENO Experiment

The RENO (Reactor Experiment for Neutrino Oscillation) experiment is under construction to measure the smallest neutrino mixing angle {theta}{sub 13} using anti-neutrinos emitted from the Yonggwang nuclear power plant in Korea. The experiment is planning to start data-taking in early 2011 with two identical 16-ton Gadolinium loaded liquid scintillator detectors located near and far from the center of the reactor array. The estimated systematic uncertainty associated with the measurement is less than 0.6%. Based on three years of data, the expected statistical error is about 0.3% and it would be sensitive to measure the neutrino mixing angle in the range, sin{sup 2}(2{theta}{sub 13})> 0.02. In this talk, the overview and status of RENO experiment are described.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Markov state models from short non-equilibrium simulations—Analysis and correction of estimation bias

To begin, many state-of-the-art methods for the thermodynamic and kinetic characterization of large and complex biomolecular systems by simulation rely on ensemble approaches, where data from large numbers of relatively short trajectories are integrated. In this context, Markov state models (MSMs) are extremely popular because they can be used to compute stationary quantities and long-time kinetics from ensembles of short simulations, provided that these short simulations are in “local equilibrium” within the MSM states. However, over the last 15 years since the inception of MSMs, it has been controversially discussed and not yet been answered how deviations from local equilibrium can be detected, whether these deviations induce a practical bias in MSM estimation, and how to correct for them. In this paper, we address these issues: We systematically analyze the estimation of MSMs from short non-equilibrium simulations, and we provide an expression for the error between unbiased transition probabilities and the expected estimate from many short simulations. We show that the unbiased MSM estimate can be obtained even from relatively short non-equilibrium simulations in the limit of long lag times and good discretization. Further, we exploit observable operator model (OOM) theory to derive an unbiased estimator for the MSM transition matrix that corrects for the effect of starting out of equilibrium, even when short lag times are used. Finally, we show how the OOM framework can be used to estimate the exact eigenvalues or relaxation time scales of the system without estimating an MSM transition matrix, which allows us to practically assess the discretization quality of the MSM. Applications to model systems and molecular dynamics simulation data of alanine dipeptide are included for illustration. The improved MSM estimator is implemented in PyEMMA of version 2.3.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coherent population trapping in Raman-pulse atom interferometry

Raman pulse atom interferometry is an important modality for precision measurements of inertial forces and tests of fundamental physics. Typical Raman atom optics use two coherent laser fields applied at gigahertz-scale detunings from optical resonance, so that spontaneous emission produces a minor or negligible source of decoherence. An additional consequence of spontaneous emission is coherent population trapping (CPT). We show that CPT produces coherences and population differences which induce systematic effects in Raman pulse atom interferometers. We do not believe that CPT has been previously identified as an error mechanism in Raman pulse atom interferometry. We present an experimental characterization of CPT coherences and population differences induced in laser-cooled cesium atoms by application of Raman pulses at detunings near 1 GHz, commensurate with detunings used in several precision measurement experiments. We are not aware of previous demonstrations of CPT-induced population difference. We argue that CPT effects could induce phase shifts of several milliradians in magnitude for typical experimental parameters and stipulate that these errors can be suppressed by propagation direction reversal in Raman interferometer-based precision measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum-accurate spectral neighbor analysis potential models for Ni-Mo binary alloys and fcc metals

In recent years, efficient interatomic potentials approaching the accuracy of density functional theory (DFT) calculations have been developed using rigorous atomic descriptors satisfying strict invariances, for example, for translation, rotation, permutation of homonuclear atoms, among others. In this work, we generalize the spectral neighbor analysis potential (SNAP) model to bcc-fcc binary alloy systems. We demonstrate that machine-learned SNAP models can yield significant improvements even over the well-established high-performing embedded atom method (EAM) and modified EAM potentials for fcc Cu and Ni. We also report on the development of a SNAP model for the fcc Ni-bcc Mo binary system by machine learning a carefully constructed large computed data set of elemental and intermetallic compounds. We demonstrate that this binary Ni-Mo SNAP model can achieve excellent agreement with experiments in the prediction of a Ni-Mo phase diagram as well as near-DFT accuracy in the prediction of many key properties, such as elastic constants, formation energies, melting points, etc., across the entire binary composition range. In contrast, the existing Ni-Mo EAM has significant errors in the prediction of the phase diagram and completely fails in binary compounds. Lastly, this paper provides a systematic model development process for multicomponent alloy systems, including an efficient procedure to optimize the hyperparameters in the model fitting, and paves the way for long-time large-scale simulations of such systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Probing the top-quark width using the charge identification of b jets

We propose a new method for measuring the top-quark width based on the on-/off-shell ratio of b -charge asymmetry in pp → Wbj production at the LHC. The charge asymmetry removes virtually all backgrounds and related uncertainties, while remaining systematic and theoretical uncertainties can be taken under control by the ratio of cross sections. Limited only by statistical error, in an optimistic scenario, we find that our approach leads to good precision at high integrated luminosity, at a few hundred MeV assuming 300 – 3000 fb -1 at the LHC. The approach directly probes the total width, in such a way that model-dependence can be minimized. It is complementary to existing cross section measurements which always leave a degeneracy between the total rate and the branching ratio, and provides valuable information about the properties of the top quark. Here, the proposal opens up new opportunities for precision top measurements using a b-charge identification algorithm.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A machine-learning approach to measuring the escape of ionizing radiation from galaxies in the reionization epoch

Recent observations of galaxies at z≳7, along with the low value of the electron scattering optical depth measured by the Planck mission, make galaxies plausible as dominant sources of ionizing photons during the epoch of reionization. However, scenarios of galaxy-driven reionization hinge on the assumption that the average escape fraction of ionizing photons is significantly higher for galaxies in the reionization epoch than in the local universe. The NIRSpec instrument on the James Webb Space Telescope (JWST) will enable spectroscopic observations of large samples of reionization-epoch galaxies. While the leakage of ionizing photons will not be directly measurable from these spectra, the leakage is predicted to have an indirect effect on the spectral slope and the strength of nebular emission lines in the rest-frame ultraviolet and optical. Here, we apply a machine learning technique known as lasso regression on mock JWST/NIRSpec observations of simulated z = 7 galaxies in order to obtain a model that can predict the escape fraction from JWST/NIRSpec data. Barring systematic biases in the simulated spectra, our method is able to retrieve the escape fraction with a mean absolute error of Δf{sub esc}≈0.12 for spectra with signal-to-noise ratio ≈ 5 at a rest-frame wavelength of 1500 Å for our fiducial simulation. This prediction accuracy represents a significant improvement over previous similar approaches.

79 ASTRONOMY AND ASTROPHYSICS↗

A COMPARISON OF STELLAR ELEMENTAL ABUNDANCE TECHNIQUES AND MEASUREMENTS

Stellar elemental abundances are important for understanding the fundamental properties of a star or stellar group, such as age and evolutionary history, as well as the composition of an orbiting planet. However, as abundance measurement techniques have progressed, there has been little standardization between individual methods and their comparisons. As a result, different stellar abundance procedures determine measurements that vary beyond the quoted error for the same elements within the same stars. The purpose of this paper is to better understand the systematic variations between methods and offer recommendations for producing more accurate results in the future. We invited a number of participants from around the world (Australia, Portugal, Sweden, Switzerland, and the United States) to calculate 10 element abundances (C, O, Na, Mg, Al, Si, Fe, Ni, Ba, and Eu) using the same stellar spectra for four stars (HD 361, HD 10700, HD 121504, and HD 202206). Each group produced measurements for each star using (1) their own autonomous techniques, (2) standardized stellar parameters, (3) a standardized line list, and (4) both standardized parameters and a line list. We present the resulting stellar parameters, absolute abundances, and a metric of data similarity that quantifies the homogeneity of the data. We conclude that standardization of some kind, particularly stellar parameters, improves the consistency between methods. However, because results did not converge as more free parameters were standardized, it is clear there are inherent issues within the techniques that need to be reconciled. Therefore, we encourage more conversation and transparency within the community such that stellar abundance determinations can be reproducible as well as accurate and precise.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimizing Long-term Monitoring of Radiation Air Dose Rates - 19169

Radiation air dose rates near the Fukushima Daiichi Nuclear Power Plant (NPP) have been steadily decreasing over the past seven years since the accident in March 2011. After seven years, the radiation monitoring program is expected to transition from short-term emergency responses to long-term monitoring. The objectives of long-term monitoring are often different from early emergency responses, since extensive data accumulation has led to an understanding of contaminant distributions and mobility. The main long-term monitoring objectives are (1) to confirm system stability and the continuing reduction of contaminant and hazard levels, (2) to provide assurance for the public, (3) to accumulate the basic datasets for scientific knowledge and future preparation, and (4) to detect changes or anomalies in contaminant mobility (if they occur), or any unexpected processes or events. In this work, we have developed a methodology for optimizing the sampling locations of radiation air dose-rate monitoring. Three steps are taken in order to determine sampling locations in a systematic manner: (1) prioritizing the critical locations, such as schools or regulatory requirement locations, (2) diversifying locations across the key environmental controls that are known to influence contaminant mobility and distributions, and (3) capturing the heterogeneity of radiation air dose rates across the domain. In the second step, we use a Gaussian mixture model to identify the representative locations in multiple environmental variables, such as elevation and land-cover types. In the third step, we use a Gaussian process model to capture and estimate the heterogeneity of air dose rates across the domain. We use the integrated dose-rate map from the Bayesian geostatistical method as the reference map and distribute the sampling in such a way as to capture the heterogeneity of the reference map. Although we assume hypothetical sets of measurement locations to evaluate the performance of the methodology and algorithms, we have demonstrated this approach using real, spatially extensive datasets, including the integrated air-dose-rate map of 2016. Our results have shown that this approach allows us to add sampling locations in a systematic manner. Increasing the number of sampling locations can better capture the heterogeneity of dose rates, although the estimation error does not decrease further after a certain number of samples. We have also found that when there are restrictions such as pre-existing monitoring locations or the ones along roads, the spatial estimation becomes poor even with the same number of monitoring locations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Two-step fabrication technique of gold tips for use in point-contact spectroscopy

For a successful point-contact spectroscopy (PCS) measurement, metallic tips of proper shape and smoothness are essential to ensure the ballistic nature of a point-contact junction. Until recently, the fabrication of Au tips suitable for use in point-contact spectroscopy has remained more of an art involving a trial and error method rather than an automated scientific process. To address these issues, we have developed a technique with which one can prepare high quality Au tips reproducibly and systematically. It involves an electronic control of the driving voltages used for an electrochemical etching of a gold wire in a HCl-glycerol mixture or a HCl solution. We find that a stopping current, below which the circuit is set to shut off, is a single very important parameter to produce an Au tip of desired shape. We present detailed descriptions for a two-step etching process for Au tips and also test results from PCS measurements using them.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

BeyondPlanck: I. Global Bayesian analysis of the Planck Low Frequency Instrument data

We describe the BEYONDPLANCK project in terms of our motivation, methodology, and main products, and provide a guide to a set of companion papers that describe each result in more detail. Building directly on experience from ESA’s Planck mission, we implemented a complete end-to-end Bayesian analysis framework for the Planck Low Frequency Instrument (LFI) observations. The primary product is a full joint posterior distribution P(ω | d), where ω represents the set of all free instrumental (gain, correlated noise, bandpass, etc.), astrophysical (synchrotron, free-free, thermal dust emission, etc.), and cosmological (cosmic microwave background – CMB – map, power spectrum, etc.) parameters. Some notable advantages of this approach compared to a traditional pipeline procedure are seamless end-to-end propagation of uncertainties; accurate modeling of both astrophysical and instrumental effects in the most natural basis for each uncertain quantity; optimized computational costs with little or no need for intermediate human interaction between various analysis steps; and a complete overview of the entire analysis process within one single framework. As a practical demonstration of this framework, we focus in particular on low-ℓ CMB polarization reconstruction with Planck LFI. In this process, we identify several important new effects that have not been accounted for in previous pipelines, including gain over-smoothing and time-variable and non-1/f correlated noise in the 30 and 44 GHz channels. Modeling and mitigating both previously known and newly discovered systematic effects, we find that all results are consistent with the ΛCDM model, and we constrained the reionization optical depth to τ = 0.066 ± 0.013, with a low-resolution CMB-based χ 2 probability to exceed of 32%. This uncertainty is about 30% larger than the official pipelines, arising from taking a more complete instrumental model into account. The marginal CMB solar dipole amplitude is 3362.7 ± 1.4 μK, where the error bar was derived directly from the posterior distribution without the need of any ad hoc instrumental corrections. We are currently not aware of any significant unmodeled systematic effects remaining in the Planck LFI data, and, for the first time, the 44 GHz channel is fully exploited in the current analysis. We argue that this framework can play a central role in the analysis of many current and future high-sensitivity CMB experiments, including LiteBIRD, and it will serve as the computational foundation of the emerging community-wide COSMOGLOBE effort, which aims to combine state-of-the-art radio, microwave, and submillimeter data sets into one global astrophysical model.

79 ASTRONOMY AND ASTROPHYSICS↗

Enhanced reaction rates in NDP analysis with neutron scattering

Neutron depth profiling (NDP) makes accessible quantitative information on a few isotopic concentration profiles ranging from the surface into the sample a few micrometers. Because the candidate analytes for NDP are few, there is little interference encountered. Furthermore, neutrons have no charge so mixed chemical states in the sample are of no direct concern. There are a few nuclides that exhibit large probabilities for neutron scattering. The effect of neutron scattering on NDP measurements has not previously been evaluated as a basis for either enhancing the reaction rates or as a source of measurement error. Hydrogen is a common element exhibiting large neutron scattering probability found in or around sample volumes being analyzed by NDP. A systematic study was conducted to determine the degree of signal change when neutron scattering occurs during analysis. The relative signal perturbation was evaluated for materials of varied neutron scattering probability, concentration, total mass, and geometry. Signal enhancements up to 50% are observed when the hydrogen density is high and in close proximity to the region of analysis with neutron beams of sub thermal energies. Greater signal enhancements for the same neutron number density are reported for thermal neutron beams. Even adhesive tape used to position the sample produces a measureable signal enhancement. Because of the shallow volume, negligible distortion of the NDP measured profile shape is encountered from neutron scattering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rapid systematic assessment of the detection and attribution of regional anthropogenic climate change

Despite being a well-established research field, the detection and attribution of observed climate change to anthropogenic forcing is not yet provided as a climate service. One reason for this is the lack of a methodology for performing tailored detection and attribution assessments on a rapid time scale. Here we develop such an approach, based on the translation of quantitative analysis into the “confidence” language employed in recent Assessment Reports of the Intergovernmental Panel on Climate Change. While its systematic nature necessarily ignores some nuances examined in detailed expert assessments, the approach nevertheless goes beyond most detection and attribution studies in considering contributors to building confidence such as errors in observational data products arising from sparse monitoring networks. When compared against recent expert assessments, the results of this approach closely match those of the existing assessments. Where there are small discrepancies, these variously reflect ambiguities in the details of what is being assessed, reveal nuances or limitations of the expert assessments, or indicate limitations of the accuracy of the sort of systematic approach employed here. Deployment of the method on 116 regional assessments of recent temperature and precipitation changes indicates that existing rules of thumb concerning the detectability of climate change ignore the full range of sources of uncertainty, most particularly the importance of adequate observational monitoring.

54 ENVIRONMENTAL SCIENCES↗

High-throughput cancer hypothesis testing with an integrated PhysiCell-EMEWS workflow

BackgroundCancer is a complex, multiscale dynamical system, with interactions between tumor cells and non-cancerous host systems. Therapies act on this combined cancer-host system, sometimes with unexpected results. Systematic investigation of mechanistic computational models can augment traditional laboratory and clinical studies, helping identify the factors driving a treatment's success or failure. However, given the uncertainties regarding the underlying biology, these multiscale computational models can take many potential forms, in addition to encompassing high-dimensional parameter spaces. Therefore, the exploration of these models is computationally challenging. We propose that integrating two existing technologiesone to aid the construction of multiscale agent-based models, the other developed to enhance model exploration and optimizationcan provide a computational means for high-throughput hypothesis testing, and eventually, optimization.ResultsIn this paper, we introduce a high throughput computing (HTC) framework that integrates a mechanistic 3-D multicellular simulator (PhysiCell) with an extreme-scale model exploration platform (EMEWS) to investigate high-dimensional parameter spaces. We show early results in applying PhysiCell-EMEWS to 3-D cancer immunotherapy and show insights on therapeutic failure. We describe a generalized PhysiCell-EMEWS workflow for high-throughput cancer hypothesis testing, where hundreds or thousands of mechanistic simulations are compared against data-driven error metrics to perform hypothesis optimization.ConclusionsWhile key notational and computational challenges remain, mechanistic agent-based models and high-throughput model exploration environments can be combined to systematically and rapidly explore key problems in cancer. These high-throughput computational experiments can improve our understanding of the underlying biology, drive future experiments, and ultimately inform clinical practice.

Agent-based model↗

Uncertainty of Rotating Shadowband Irradiometers and Si-Pyranometers Including the Spectral Irradiance Error

Concentrating solar power projects require accurate direct normal irradiance (DNI) data including uncertainty specifications for plant layout and cost calculations. Ground measured data are necessary to obtain the required level of accuracy and are often obtained with Rotating Shadowband Irradiometers (RSI) that use photodiode pyranometers and correction functions to account for systematic effects. The uncertainty of Si-pyranometers has been investigated, but so far basically empirical studies were published or decisive uncertainty influences had to be estimated based on experience in analytical studies. One of the most crucial estimated influences is the spectral irradiance error because Si-photodiode-pyranometers only detect visible and color infrared radiation and have a spectral response that varies strongly within this wavelength interval. Furthermore, analytic studies did not discuss the role of correction functions and the uncertainty introduced by imperfect shading. In order to further improve the bankability of RSI and Si-pyranometer data, a detailed uncertainty analysis following the Guide to the Expression of Uncertainty in Measurement (GUM) has been carried out. The study defines a method for the derivation of the spectral error and spectral uncertainties and presents quantitative values of the spectral and overall uncertainties. Data from the PSA station in southern Spain was selected for the analysis. Average standard uncertainties for corrected 10 min data of 2% for global horizontal irradiance (GHI), and 2.9% for DNI (for GHI and DNI over 300 W/m2) were found for the 2012 yearly dataset when separate GHI and DHI calibration constants were used. Also the uncertainty in 1 min resolution was analyzed. The effect of correction functions is significant. The uncertainties found in this study are consistent with results of previous empirical studies.

14 SOLAR ENERGY↗

Performance Analysis of OpenSHMEM Applications with TAU Commander

The TAU Performance System® (TAU) is a powerful and highly versatile profiling and tracing tool ecosystem for performance engineering of parallel programs. Developed over the last twenty years, TAU has evolved with each new generation of HPC systems and scales efficiently to hundreds of thousands of cores. TAU’s organic growth has resulted in a loosely coupled software toolbox such that novice users first encountering TAU’s complexity and vast array of features are often intimidated and easily frustrated. To lower the barrier to entry for novice TAU users, ParaTools and the US Department of Energy have developed “TAU Commander,” a performance engineering workflow manager that facilitates a systematic approach to performance engineering, guides users through common profiling and tracing workflows, and offers constructive feedback in case of error. This work compares TAU and TAU Commander workflows for common performance engineering tasks in OpenSHMEM applications and demonstrates workflows targeting two different SHMEM implementations, Intel Xeon “Haswell” and “Knights Landing” processors, direct and indirect measurement methods, callsite, profiles, and traces.

97 MATHEMATICS AND COMPUTING↗