Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical 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,693 records · Page 94

Deeply uncertain: comparing methods of uncertainty quantification in deep learning algorithms

We present a comparison of methods for uncertainty quantification (UQ) in deep learning algorithms in the context of a simple physical system. Three of the most common uncertainty quantification methods - Bayesian Neural Networks (BNN), Concrete Dropout (CD), and Deep Ensembles (DE) - are compared to the standard analytic error propagation. We discuss this comparison in terms endemic to both machine learning ("epistemic" and "aleatoric") and the physical sciences ("statistical" and "systematic"). The comparisons are presented in terms of simulated experimental measurements of a single pendulum - a prototypical physical system for studying measurement and analysis techniques. Our results highlight some pitfalls that may occur when using these UQ methods. For example, when the variation of noise in the training set is small, all methods predicted the same relative uncertainty independently of the inputs. This issue is particularly hard to avoid in BNN. On the other hand, when the test set contains samples far from the training distribution, we found that no methods sufficiently increased the uncertainties associated to their predictions. This problem was particularly clear for CD. In light of these results, we make some recommendations for usage and interpretation of UQ methods.

59 BASIC BIOLOGICAL SCIENCES↗

Machine learning for neutron reflectometry data analysis of two-layer thin films

Neutron reflectometry (NR) is a powerful tool for probing thin films at length scales down to nanometers. We investigated the use of a neural network to predict a two-layer thin film structure to model a given measured reflectivity curve. Application of this neural network to predict a thin film structure revealed that it was accurate and could provide an excellent starting point for traditional fitting methods. Employing prediction-guided fitting has considerable potential for more rapidly producing a result compared to the labor-intensive but commonly-used approach of trial and error searches prior to refinement. A deeper look at the stability of the predictive power of the neural network against statistical fluctuations of measured reflectivity profiles showed that the predictions are stable. We conclude that the approach presented here can provide valuable assistance to users of NR and should be further extended for use in studies of more complex n-layer thin film systems. This result also opens up the possibility of developing adaptive measurement systems in the future.

36 MATERIALS SCIENCE↗

Block Lanczos algorithm for lattice QCD spectroscopy and matrix elements

Recent work introduced a new framework for analyzing correlation functions with improved convergence and signal-to-noise properties, as well as rigorous quantification of excited-state effects, based on the Lanczos algorithm and spurious eigenvalue filtering with the Cullum-Willoughby test. Here, we extend this framework to the analysis of correlation-function matrices built from multiple interpolating operators in lattice quantum chromodynamics (QCD) by constructing an oblique generalization of the block Lanczos algorithm, as well as a new physically motivated reformulation of the Cullum-Willoughby test that generalizes to block Lanczos straightforwardly. The resulting block Lanczos method directly extends generalized eigenvalue problem (GEVP) methods, which can be viewed as applying a single iteration of block Lanczos. Block Lanczos provides qualitative and quantitative advantages over GEVP methods analogous to the benefits of Lanczos over the standard effective mass, including faster convergence to ground- and excited-state energies, explicitly computable two-sided error bounds, straightforward extraction of matrix elements of external currents, and asymptotically constant signal-to-noise. No fits or statistical inference are required. Proof-of-principle calculations are performed for noiseless mock-data examples as well as two-by-two proton correlation-function matrices in lattice QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Robust Medium-Voltage Distribution System State Estimation using Multi-Source Data

Due to the lack of sufficient online measurements for distribution system observability, pseudo-measurements from short-term load or distributed renewable energy resources (DERs) forecasting are used. However, the accuracy of them is low and thus significantly limits the performance of distribution system state estimation (DSSE). In this paper, a robust DSSE that integrates multi-source measurement data is proposed. Specifically, the historical low-voltage (LV) side smart meters are used to forecast load and DERs injections via the support vector machine (SVM) with optimally tuned parameters. By contrast, the online smart meters at LV side are utilized to derive equivalent power injections at the MV/LV transformers, yielding more accurate pseudo-measurements compared to the forecasted injections. Furthermore, to deal with bad data caused by communication loss, instrumental errors and cyber attacks, robust DSSE that relies on generalized maximum-likelihood (GM)-estimation criterion is developed. The projection statistics are developed to adjust the weights of each measurement, leading to better balance between pseudo- and real-time measurements. Numerical results conducted on modified IEEE 33-bus system with DG integration demonstrate the effectiveness and robustness of the proposed method.

distribution system state estimation↗

SU-E-T-550: Modulation Index for VMAT

Purpose: To present modulation indices (MIs) for volumetric modulated arc therapy (VMAT). Methods: A total of 40 VMAT plans were retrospectively selected. To investigate the delivery accuracy of each VMAT plan, gamma passing rates, differences in modulating parameters between plans and log files, and differences between the original plans and the plans reconstructed with the log files were acquired. A modulation index (MIt) was designed by multiplications of the weighted quantifications of MLC speeds, MLC accelerations, gantry accelerations and dose-rate variations. Textural features including angular second moment, inverse difference moment, contrast, variance, correlation and entropy were calculated from the fluences of each VMAT plan. To test the performance of suggested MIs, Spearman’s rank correlation coefficients (r) with the plan delivery accuracy were calculated. Conventional modulation indices for VMAT including the modulation complexity score for VMAT (MCSv), leaf travel modulation complexity score (LTMCS) and MI by Li & Xing were calculated, and their correlations were also analyzed in the same way. Results: The r values of contrast (particular displacement distance, d = 1), variance (d = 1), MIt, MCSv, LTMCS and MI by Li&Xing to the local gamma passing rates with 2%/2 mm were 0.547 (p < 0.001), 0.519 (p < 0.001), −0.658 (p < 0.001), 0.186 (p = 0.251), 0.312 (p = 0.05) and −0.455 (p = 0.003), respectively. The r values of those to the MLC errors were −0.863, −0.828, 0.917, −0.635, − 0.857 and 0.795, respectively (p < 0.001). For dose-volumetric parameters, MIt showed higher statistically significant correlations than did the conventional modulation indices. Conclusion: The MIt, contrast (d = 1) and variance (d = 1) showed good performance to predict the VMAT delivery accuracy showing higher correlations to the results of various types of verification methods for VMAT. This work was in part supported by the National Research Foundation of Korea (NRF) grant (no. 490-20140029 and no. 490-20130047) funded by the Korea government.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-C-BRB-05: Determining the Adequacy of Auto-Contouring Via Probabilistic Assessment of Ensuing Treatment Plan Metrics in Comparison with Manual Contours

Purpose: To determine if auto-contour and manual-contour—based plans differ when evaluated with respect to probabilistic coverage metrics and biological model endpoints for prostate IMRT. Methods: Manual and auto-contours were created for 149 CT image sets acquired from 16 unique prostate patients. A single physician manually contoured all images. Auto-contouring was completed utilizing Pinnacle’s Smart Probabilistic Image Contouring Engine (SPICE). For each CT, three different 78 Gy/39 fraction 7-beam IMRT plans are created; PD with drawn ROIs, PAS with auto-contoured ROIs, and PM with auto-contoured OARs with the manually drawn target. For each plan, 1000 virtual treatment simulations with different sampled systematic errors for each simulation and a different sampled random error for each fraction were performed using our in-house GPU-accelerated robustness analyzer tool which reports the statistical probability of achieving dose-volume metrics, NTCP, TCP, and the probability of achieving the optimization criteria for both auto-contoured (AS) and manually drawn (D) ROIs. Metrics are reported for all possible cross-evaluation pairs of ROI types (AS,D) and planning scenarios (PD,PAS,PM). Bhattacharyya coefficient (BC) is calculated to measure the PDF similarities for the dose-volume metric, NTCP, TCP, and objectives with respect to the manually drawn contour evaluated on base plan (D-PD). Results: We observe high BC values (BC≥0.94) for all OAR objectives. BC values of max dose objective on CTV also signify high resemblance (BC≥0.93) between the distributions. On the other hand, BC values for CTV’s D95 and Dmin objectives are small for AS-PM, AS-PD. NTCP distributions are similar across all evaluation pairs, while TCP distributions of AS-PM, AS-PD sustain variations up to %6 compared to other evaluated pairs. Conclusion: No significant probabilistic differences are observed in the metrics when auto-contoured OARs are used. The prostate auto-contour needs improvement to achieve clinically equivalent plans.

60 APPLIED LIFE SCIENCES↗

Efficient design of random metasurfaces

Random media introduce large degrees of freedom in device design and can thus address challenges in manipulating optical waves. Wave shaping with metasurfaces has mainly utilized periodic or quasi-periodic grids, and the potential of random arrangement of particles for devices has only come under investigation recently. The main difficulty in pursuing random metasurfaces is the identification of the degrees of freedom that optimize their efficiencies and functions. They can also encode information using the statistics of particle distribution. We propose a phase-map that accounts for the statistical nature of random media. The method takes into account effects of random near-field couplings that introduce phase errors by affecting the phase shift of elements. The proposed approach increases the efficiency of our random metasurface devices by up to ~20%. Furthermore, this work paves the way toward the efficient design of random metasurfaces with potential applications in highly secure optical cryptography and information encoding.

Nasari, Hadiseh↗

High-resolution hypernuclear decay pion spectroscopy at MAMI and future

Hypernuclear decay pion spectroscopy was established in 2012 at MAMI as a mass spectroscopy method for light hypernuclei. A monochromatic pion peak from $^4_Λ$H was successfully observed, and the Λ binding energy was determined to be B Λ = 2.157±0.005(stat.)±0.077(syst.) MeV in the 2014 run. In 2022, an upgrade experiment for $^3_Λ$H spectroscopy was conducted using a newly developed Li target. The absolute electron beam energy will be measured by the synchrotron radiation interferometry, which will be applied with the spectrometer calibration to improve the systematic error. The decay pion spectroscopy is planned to be performed at JLab Hall-C, which would significantly increase the statistics thanks to the higher energy beam, the better K + identification, and the faster data acquisition system. This experiment has been submitted as a Letter of Intent in JLab PAC51. The upgraded decay pion spectroscopy method is expected to provide new, accurate hypernuclear data, which will contribute to the advancement of our understanding of hypernuclear physics

47 OTHER INSTRUMENTATION↗

Analysis and Synthesis of Load Forecasting Data for Renewable Integration Studies: Preprint

As renewable energy constitutes greater portions of the generation fleet, the importance of modeling uncertainty as part of integration studies also increases. In pursuit of optimal system operations, it is important to capture not only the definitive behavior of power plants, but also the risks associated with systemwide interactions. This research examines the dependence of load forecast errors on external predictor variables such as temperature, day type, and time of day. The analysis was utilized to create statistically relevant instances of sequential load forecasts with only a time series of historic, measured load available. The creation of such load forecasts relies on Bayesian techniques for informing and updating the model, thus providing a basis for networked and adaptive load forecast models in future operational applications.

17 WIND ENERGY↗

Model-Form Epistemic Uncertainty Quantification for Modeling with Differential Equations: Application to Epidemiology

Modeling real-world phenomena to any degree of accuracy is a challenge that the scientific research community has navigated since its foundation. Lack of information and limited computational and observational resources necessitate modeling assumptions which, when invalid, lead to model-form error (MFE). The work reported herein explored a novel method to represent model-form uncertainty (MFU) that combines Bayesian statistics with the emerging field of universal differential equations (UDEs). The fundamental principle behind UDEs is simple: use known equational forms that govern a dynamical system when you have them; then incorporate data-driven approaches – in this case neural networks (NNs) – embedded within the governing equations to learn the interacting terms that were underrepresented. Utilizing epidemiology as our motivating exemplar, this report will highlight the challenges of modeling novel infectious diseases while introducing ways to incorporate NN approximations to MFE. Prior to embarking on a Bayesian calibration, we first explored methods to augment the standard (non-Bayesian) UDE training procedure to account for uncertainty and increase robustness of training. In addition, it is often the case that uncertainty in observations is significant; this may be due to randomness or lack of precision in the measurement process. This uncertainty typically manifests as “noisy” observations which deviate from a true underlying signal. To account for such variability, the NN approximation to MFE is endowed with a probabilistic representation and is updated using available observational data in a Bayesian framework. By representing the MFU explicitly and deploying an embedded, data-driven model, this approach enables an agile, expressive, and interpretable method for representing MFU. In this report we will provide evidence that Bayesian UDEs show promise as a novel framework for any science-based, data-driven MFU representation; while emphasizing that significant advances must be made in the calibration of Bayesian NNs to ensure a robust calibration procedure.

97 MATHEMATICS AND COMPUTING↗

The Backup-Gilbert method and its application to the electrical conductivity problem

The theory of Backus and Gilbert gives a technique for solving the general linear inverse problem. Observational error and lack of data are shown to reduce the reliability of the solution in different ways: the former introduces statistical uncertainties in the model, while the latter smooths out the detail. Precision can be improved by sacrificing resolving power, and vice versa, so that some compromise may be made between the two in choosing the best model. Nonlinear inverse problems can be brought into the domain of the theory by linearizing about a typical solution. The inverse problem of electrical conductivity in the mantle is used to illustrate the Backus-Gilbert technique; an example of the tradeoff diagram is given.

Parker, R. L.↗

Estimation of terrain iso-gradients from a stochastic range data measurement matrix

The problem of estimating terrain iso-gradients for an autonomous roving vehicle is complicated by the measurement error inherent in the range data matrix. In this paper, the in-path and cross-path slopes at each data point in the multivariable range matrix are expressed in terms of parameters at that data point and those at its surrounding points. The sensitivity in the change of these slopes - due to the errors in measurements in range, elevation angle and azimuth angle are formulated. These sensitivities are used to determine the appropriate statistics of the gradient, which are then used to describe the terrain with a known probability of accuracy.

Shen, C. N.↗

The HEAO A-2 survey of Abell clusters and the X-ray luminosity function

The HEAO A-2 all sky data base was surveyed for 2-10KeV X-rau emission from the 225 Abell clusters og galaxies listed in Abell's (1958) catalog which are of distance class four or less, and are within the fraction of the sky surveyed completely by Abell. Thirty-two identifications of clusters with X-ray sources were made, for which 2-10 KeV fluxes and 90% error boxes are presented. Twelve of these identification are new. The X-ray luminosity function was derived for this statistically complete sample and the best exponential fit was found to be f(L) = 20.2 x 10 to the minus 8 power exp (-l(44)/1.9) per cu Mpc 2-10KeV. The relationship between X-ray luminosity and richness was examined and a correlation was found for richness classes 0, 1, and 2. The relationship of X-rau luminosity, Bautz-Morgan type, and Rood-Sastry type was examined. It was found that BM type I's and RS type cD and B have the greatest average luminosity. The contribution of clusters to the X-ray background was calculated from the luminosity function and was found to be 5%, and with 90% certainty, less than 8% in the 2-10 KeV band pass.

Mckee, J. D.↗

A comparison of head and manual control for a position-control pursuit tracking task

Head control was compared with manual control in a pursuit tracking task involving proportional controlled-element dynamics. An integrated control/display system was used to explore tracking effectiveness in horizontal and vertical axes tracked singly and concurrently. Compared with manual tracking, head tracking resulted in a 50 percent greater rms error score, lower pilot gain, greater high-frequency phase lag and greater low-frequency remnant. These differences were statistically significant, but differences between horizontal- and vertical-axis tracking and between 1- and 2-axis tracking were generally small and not highly significant. Manual tracking results were matched with the optimal control model using pilot-related parameters typical of those found in previous manual control studies. Head tracking performance was predicted with good accuracy using the manual tracking model plus a model for head/neck response dynamics obtained from the literature.

Levison, W. H.↗

A variational assimilation method for satellite and conventional data: Development of basic model for diagnosis of cyclone systems

A three-dimensional diagnostic model for the assimilation of satellite and conventional meteorological data is developed with the variational method of undetermined multipliers. Gridded fields of data from different type, quality, location, and measurement source are weighted according to measurement accuracy and merged using least squares criteria so that the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation are satisfied. The model is used to compare multivariate variational objective analyses with and without satellite data with initial analyses and the observations through criteria that were determined by the dynamical constraints, the observations, and pattern recognition. It is also shown that the diagnoses of local tendencies of the horizontal velocity components are in good comparison with the observed patterns and tendencies calculated with unadjusted data. In addition, it is found that the day-night difference in TOVS biases are statistically different (95% confidence) at most levels. Also developed is a hybrid nonlinear sigma vertical coordinate that eliminates hydrostatic truncation error in the middle and upper troposphere and reduces truncation error in the lower troposphere. Finally, it is found that the technique used to grid the initial data causes boundary effects to intrude into the interior of the analysis a distance equal to the average separation between observations.

Achtemeier, G. L.↗

Mars Rover/Sample Return landing strategy

This paper describes the analysis and results of an investigation of the Mars Rover/Sample Return mission's landing strategy, together with the trade-offs of different landing strategies. The percentile points and the descriptive statistics of the probability distribution of traverse distances are calculated using a simple model formulated on the basis of landing error characteristics. The results show that variations in the landing stratgegy can significantly affect the traverse distance requirements, which range from 20 to 200 km.

Friedlander, Alan L.↗

Data Assimilation Into a Coupled Ocean Atmosphere Model: Application to the 1997-1998 El Nino

As part of JPL's ocean data assimilation effort to study ocean circulation and seasonal-interannual climate variability, sea level anomaly observed by TOPEX altimeter, together with sea surface temperature and wind stress data, are assimilated into a simple coupled ocean atmosphere model of the tropical Pacific. Model-data consistency is examined. Impact of the assimilation (as initialization) on El Nino Southern Oscillation (ENSO) forecasts is evaluated. The coupled model consists of a shallow water component with two baroclinic modes, an Ekman shear layer, a simplified mixed-layer temperature equation, and a statistical atmosphere based on dominant correlations between historical surface temperature and wind stress anomaly data. The adjoins method is used to fit the coupled model to the data over various six-month periods from late 1996 to early 1998 by optimally adjusting the initial state, model parameters, and basis functions of the statistical atmosphere. On average, the coupled model can be fitted to the data to approximately within the data and representation errors (5 cm, 0.5 C, and 10 sq m/sq m for sea level, surface temperature, and pseudo wind stress anomalies, respectively). The estimated fields resemble observed spatio-temporal structure reasonably well. Hindcasts/forecasts of the 1997/1998 El Nino initialized from forced estimated ocean states and parameters are much more realistic than those simply initialized from ocean states (see figure below). In particular, the ability of the model to produce significant warming beyond the initial state is dramatically improved. Parameter estimation, which compensates for some model errors, is found to be important to obtaining better fits of the model to data and to improving forecasts.

Lee, Tong↗

(abstract) A VLBI Test of Tropospheric Delay Calibration with WVRs

Dual frequency (S/X band) very long baseline interferometry (VLBI) observations were used to test troposphere calibration by water vapor radiometers (WVRs). Comparison of the VLBI and WVR measurements show a statistical agreement (specifically, their structure functions agree) on time scales less than 700 seconds. On longer time scales, VLBI instrumental errors become important. The improvement in VLBI residual delays from WVR calibration was consistent with the measured level of tropospheric fluctuations.

VLBI water vapor radiometers troposphere interfero↗