Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical techniques”

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,081 records · Page 60

A performance model for GPUs with caches

To exploit the abundant computational power of the world's fastest supercomputers, an even workload distribution to the typically heterogeneous compute devices is necessary. While relatively accurate performance models exist for conventional CPUs, accurate performance estimation models for modern GPUs do not exist. This paper presents two accurate models for modern GPUs: a sampling-based linear model, and a model based on machine-learning (ML) techniques which improves the accuracy of the linear model and is applicable to modern GPUs with and without caches. We first construct the sampling-based linear model to predict the runtime of an arbitrary OpenCL kernel. Based on an analysis of NVIDIA GPUs' scheduling policies we determine the earliest sampling points that allow an accurate estimation. The linear model cannot capture well the significant effects that memory coalescing or caching as implemented in modern GPUs have on performance. We therefore propose a model based on ML techniques that takes several compiler-generated statistics about the kernel as well as the GPU's hardware performance counters as additional inputs to obtain a more accurate runtime performance estimation for modern GPUs. We demonstrate the effectiveness and broad applicability of the model by applying it to three different NVIDIA GPU architectures and one AMD GPU architecture. On an extensive set of OpenCL benchmarks, on average, the proposed model estimates the runtime performance with less than 7 percent error for a second-generation GTX 280 with no on-chip caches and less than 5 percent for the Fermi-based GTX 580 with hardware caches. On the Kepler-based GTX 680, the linear model has an error of less than 10 percent. On an AMD GPU architecture, Radeon HD 6970, the model estimates with 8 percent of error rates. As a result, the proposed technique outperforms existing models by a factor of 5 to 6 in terms of accuracy.

97 MATHEMATICS AND COMPUTING↗

The Discovery of λ Bootis Stars: The Southern Survey I

The λ Boo stars are a class of chemically peculiar Population I A-type stars characterized by under-abundances of the refractory elements, but near-solar abundances of carbon, nitrogen, oxygen, and sulfur. There is some evidence that λ Boo stars have higher frequencies of “bright” debris disks than normal A-type stars. The discovery of four exoplanets orbiting HR 8799, a λ Boo star with a resolved debris disk, suggests that the λ Boo phenomenon may be related to the presence of a dynamic debris disk, perhaps perturbed by migrating planets. However, only 64 λ Boo stars are known, and those stars were discovered with different techniques, making it problematic to use that sample for statistical purposes, including determining the frequency of debris disks. The purpose of this paper is to derive a new sample of λ Boo stars using a technique that does not lead to biases with respect to the presence of infrared excesses. Through spectroscopic observations in the southern hemisphere, we have discovered 33 λ Boo stars and have confirmed 12 others. As a step toward determining the proportion of λ Boo stars with infrared excesses, we have used WISE data to examine the infrared properties of this sample out to 22 μ m. On this basis, we cannot conclude that λ Boo stars have a greater tendency than normal A-type stars to show infrared excesses. However, observing this sample at longer wavelengths may change that conclusion, as many λ Boo debris disks are cool and do not radiate strongly at 22 μ m.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistics of Delta v magnitude for a trajectory correction maneuver containing deterministic and random components

A number of interplanetary missions now being planned involve placing deterministic maneuvers along the flight path to alter the trajectory. Lee and Boain (1973) examined the statistics of trajectory correction maneuver (TCM) magnitude with no deterministic ('bias') component. The Delta v vector magnitude statistics were generated for several values of random Delta v standard deviations using expansions in terms of infinite hypergeometric series. The present investigation uses a different technique (Monte Carlo simulation) to generate Delta v magnitude statistics for a wider selection of random Delta v standard deviations and also extends the analysis to the case of nonzero deterministic Delta v's. These Delta v magnitude statistics are plotted parametrically. The plots are useful in assisting the analyst in quickly answering questions about the statistics of Delta v magnitude for single TCM's consisting of both a deterministic and a random component. The plots provide quick insight into the nature of the Delta v magnitude distribution for the TCM.

Bollman, W. E.↗

Gaussian process based optimization of molecular geometries using statistically sampled energy surfaces from quantum Monte Carlo

Optimization of atomic coordinates and lattice parameters remains a significant challenge to the wide use of stochastic electronic structure methods such as quantum Monte Carlo (QMC). Measurements of the forces and stress tensor by these methods contain statistical errors, challenging conventional gradient-based numerical optimization methods that assume deterministic results. Additionally, forces are not yet available for some methods, wavefunctions, and basis sets and when available may be expensive to compute to sufficiently high statistical accuracy near energy minima, where the energy surfaces are flat. Here, we explore the use of Gaussian process based techniques to sample the energy surfaces and reduce sensitivity to the statistical nature of the problem. We utilize Latin hypercube sampling, with the number of sampled energy points scaling quadratically with the number of optimized parameters. Furthermore, we show these techniques may be successfully applied to systems consisting of tens of parameters, demonstrating QMC optimization of a benzene molecule starting from a randomly perturbed, broken symmetry geometry.

74 ATOMIC AND MOLECULAR PHYSICS↗

Sci-Thur PM - Colourful Interactions: Highlights 08: ARC TBI using Single-Step Optimized VMAT Fields

Purpose: This work outlines a new TBI delivery technique to replace a lateral POP full bolus technique. The new technique is done with VMAT arc delivery, without bolus, treating the patient prone and supine. The benefits of the arc technique include: increased patient experience and safety, better dose conformity, better organ at risk sparing, decreased therapist time and reduction of therapist injuries. Methods: In this work we build on a technique developed by Jahnke et al. We use standard arc fields with gantry speeds corrected for varying distance to the patient followed by a single step VMAT optimization on a patient CT to increase dose inhomogeneity and to reduce dose to the lungs (vs. blocks). To compare the arc TBI technique to our full bolus technique, we produced plans on patient CTs for both techniques and evaluated several dosimetric parameters using an ANOVA test. Results and Conclusions: The arc technique is able reduce both the hot areas to the body (D2% reduced from 122.2% to 111.8% p<0.01) and the lungs (mean lung dose reduced from 107.5% to 99.1%, p<0.01), both statistically significant, while maintaining coverage (D98% = 97.8% vs. 94.6%, p=0.313, not statistically significant). We developed a more patient and therapist-friendly TBI treatment technique that utilizes single-step optimized VMAT plans. It was found that this technique was dosimetrically equivalent to our previous lateral technique in terms of coverage and statistically superior in terms of reduced lung dose.

60 APPLIED LIFE SCIENCES↗

Design, Detection, and Countermeasure of Frequency Spectrum Attack and Its Impact on Long Short-Term Memory Load Forecasting and Microgrid Energy Management

This paper introduces a frequency-domain false data injection attack called Frequency Spectrum Attack (FSA) and explores its effects on load forecasting and the energy management system (EMS) in a microgrid. The FSA analyzes time-series signals in the frequency domain to identify patterns in their frequency spectrum. It learns the distribution of dominant frequencies in a dataset of healthy signals. Subsequently, it manipulates the amplitudes of dominant frequencies within this healthy distribution, ensuring a stealthy attack against statistical analysis of the signal spectrum. We evaluated the performance of FSA on LSTM, a state-of-the-art network for load forecasting. The results show that FSA can triple the Mean Absolute Error (MAE) of predictions compared to the normal case and increase it by 70% compared to noise injection attacks. Furthermore, FSA indirectly enhances battery utilization in the EMS by 45%. We then proposed a detection method that combines statistical analysis and machine-learning-based classification techniques with features. The model effectively distinguishes FSA from healthy and noisy signals, achieving an accuracy of 98.7% and an F1-score of 98.1% on a load dataset, covering healthy, FSA, and noisy load data. Finally, a countermeasure was introduced based on the statistical analysis of the frequency spectrum of healthy signals to mitigate the impact of FSA. This countermeasure successfully reduces the MAE of the attacked model from 0.135 to 0.053, validating its effectiveness in mitigating FSA.

Nazeri, Amirhossein↗

Sandia Fracture Challenge 3: detailing the Sandia Team Q failure prediction strategy

The third Sandia Fracture Challenge illustrated the geometric and material uncertainties presented by modern additive manufacturing techniques. Tasked with the challenge of predicting failure of a complex additively-manufactured geometry made of 316L stainless steel, we combined a rigorous material calibration scheme with a number of statistical assessments of problem uncertainties. Particularly, we used optimization techniques to calibrate a rate-dependent and anisotropic Hill plasticity model to represent material deformation coupled with a damage model driven by void growth and nucleation. Through targeted simulation studies we assessed the influence of internal voids and surface flaws on the specimens of interest in the challenge which guided our material modeling choices. Employing the Kolmogorov–Smirnov test statistic, we developed a representative suite of simulations to account for the geometric variability of test specimens and the variability introduced by material parameter uncertainty. This method allowed the team to successfully predict the failure mode of the experimental test population as well as the global response with a high degree of accuracy.

316L stainless steel↗

Improvements to the macroscopic-microscopic approach of nuclear fission

The well-established macroscopic-microscopic (mac-mic) description of nuclear fission enables the prediction of fission-fragment yields for a broad range of fissioning systems. In this work, we present several key enhancements to this approach. We improve upon the microscopic sector of nuclear potential-energy surfaces by magnifying the resolution of the Lipkin-Nogami equations and strengthening the Strutinsky procedure, thus reducing spurious effects from the continuum. We further present a novel deterministic method for calculating fission dynamics under the assumption of strongly damped nucleonic motion. Our technique directly determines the evolution of the scissioned shape distribution according to the number of random-walk steps rather than the statistical accumulation of fission events. We show that our new technique is equivalent to the Metropolis random walk pioneered over the past decade by Randrup and colleagues. It further improves upon it because we remove the need for altering the nuclear landscape via a biased potential. With our final improvement, we calculate fission fragments mass and charge distributions using particle number projection, which affords the simultaneous calculation of both mass- and charge-yield distributions. Fission fragments are thus calculated from the quantum-mechanical A -body states of the potential-energy surface rather than from the collective mass asymmetry variable α g of the finite-range liquid-drop model used in past work. We highlight the success of our enhancements by predicting the odd-even staggering and the charge polarization for the neutron-induced fission of 233 U and 235 U .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Machine Learning Predictions of a Multiresolution Climate Model Ensemble

Statistical models of high-resolution climate models are useful for many purposes, including sensitivity and uncertainty analyses, but building them can be computationally prohibitive. We generated a unique multiresolution perturbed parameter ensemble of a global climate model. We use a novel application of a machine learning technique known as random forests to train a statistical model on the ensemble to make high-resolution model predictions of two important quantities: global mean top-of-atmosphere energy flux and precipitation. The random forests leverage cheaper low-resolution simulations, greatly reducing the number of high-resolution simulations required to train the statistical model. We demonstrate that high-resolution predictions of these quantities can be obtained by training on an ensemble that includes only a small number of high-resolution simulations. Finally, we also find that global annually averaged precipitation is more sensitive to resolution changes than to any of the model parameters considered.

54 ENVIRONMENTAL SCIENCES↗

Statistical fluctuations in Monte Carlo calculations

The time counter and modified Nanbu simulation techniques are analyzed, with emphasis placed on the convergence of the calculations to a steady macroscopic state. Such variables as translational and rotational temperature, and flow velocity, sampled at several points in the flowfield, are considered. Both macroscopic averages and molecular distribution functions are analyzed. The calculation of inelastic collisions, in which transfer of energy between translational and internal energy modes is performed, is achieved through the use of the Larsen-Borgnakke phenomenological model. It is noted that, with reference to translational temperature, the time counter method shows less statistical scatter than that found with the modified Nanbu simulation technique.

Boyd, I. D.↗

Improved profile fitting and quantification of uncertainty in experimental measurements of impurity transport coefficients using Gaussian process regression

The need to fit smooth temperature and density profiles to discrete observations is ubiquitous in plasma physics, yet the prevailing techniques for this have many shortcomings that cast doubt on the statistical validity of the results. This issue is amplified in the context of validation of gyrokinetic transport models (Holland et al 2009 Phys. Plasmas 16 052301), where the strong sensitivity of the code outputs to input gradients means that inadequacies in the profile fitting technique can easily lead to an incorrect assessment of the degree of agreement with experimental measurements. In order to rectify the shortcomings of standard approaches to profile fitting, we have applied Gaussian process regression (GPR), a powerful non-parametric regression technique, to analyse an Alcator C-Mod L-mode discharge used for past gyrokinetic validation work (Howard et al 2012 Nucl. Fusion 52 063002). We demonstrate that the GPR techniques can reproduce the previous results while delivering more statistically rigorous fits and uncertainty estimates for both the value and the gradient of plasma profiles with an improved level of automation. Moreover, we discuss how the use of GPR can allow for dramatic increases in the rate of convergence of uncertainty propagation for any code that takes experimental profiles as inputs. The new GPR techniques for profile fitting and uncertainty propagation are quite useful and general, and we describe the steps to implementation in detail in this paper. These techniques have the potential to substantially improve the quality of uncertainty estimates on profile fits and the rate of convergence of uncertainty propagation, making them of great interest for wider use in fusion experiments and modelling efforts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct Detection Doppler Lidar for Spaceborne Wind Measurement

Aerosol and molecular based versions of the double-edge technique can be used for direct detection Doppler lidar spaceborne wind measurement. The edge technique utilizes the edge of a high spectral resolution filter for high accuracy wind measurement using direct detection lidar. The signal is split between an edge filter channel and a broadband energy monitor channel. The energy monitor channel is used for signal normalization. The edge measurement is made as a differential frequency measurement between the outgoing laser signal and the atmospheric backscattered return for each pulse. As a result the measurement is insensitive to laser and edge filter frequency jitter and drift at a level less than a few parts in 10(exp 10). We have developed double edge versions of the edge technique for aerosol and molecular-based lidar measurement of the wind. Aerosol-based wind measurements have been made at Goddard Space Flight Center and molecular-based wind measurements at the University of Geneva. We have demonstrated atmospheric measurements using these techniques for altitudes from 1 to more than 10 km. Measurement accuracies of better than 1.25 m/s have been obtained with integration times from 5 to 30 seconds. The measurements can be scaled to space and agree, within a factor of two, with satellite-based simulations of performance based on Poisson statistics. The theory of the double edge aerosol technique is described by a generalized formulation which substantially extends the capabilities of the edge technique. It uses two edges with opposite slopes located about the laser frequency at approximately the half-width of each edge filter. This doubles the signal change for a given Doppler shift and yields a factor of 1.6 improvement in the measurement accuracy compared to the single edge technique. The use of two high resolution edge filters substantially reduces the effects of Rayleigh scattering on the measurement, as much as order of magnitude, and allows the signal to noise ratio to be substantially improved in areas of low aerosol backscatter. We describe a method that allows the Rayleigh and aerosol components of the signal to be independently determined using the two edge channels and an energy monitor channel. The effects of Rayleigh scattering may then subtracted from the measurement and we show that the correction process does not significantly increase the measurement noise for Rayleigh to aerosol ratios up to 10. We show that for small Doppler shifts a measurement accuracy of 0.4 m/s can be obtained for 5000 detected photon, 1.2 m/s for 1000 detected photons, and 3.7 m/s for 50 detected photons for a Rayleigh to aerosol ratio of 5. Methods for increasing the dynamic range of the aerosol-based system to more than +/- 100 m/s are given.

Korb, C. Laurence↗

Data analytics for intermodal freight transportation applications

With the growth of intermodal freight transportation, it is important that transportation planners and decision-makers are knowledgeable about freight flow data to make informed decisions. This is particularly true with Intelligent Transportation Systems (ITS) offering new capabilities for intermodal freight transportation. Specifically, ITS enables access to multiple different data sources, but they have different formats, resolutions, and time scales. Thus, knowledge of data science is essential to be successful in future ITS-enabled intermodal freight transportation systems. This chapter discusses the commonly used descriptive and predictive data analytic techniques in intermodal freight transportation applications. These techniques cover the entire spectrum of univariate, bivariate, and multivariate analyses. In addition to illustrating how to apply these techniques manually, this chapter will also show how to apply them using the statistical software R. Additional exercises are provided for those who wish to apply the described techniques to more complex problems.

Huynh, Nathan↗

A statistical approach to combining multisource information in one-class classifiers

A new method is introduced in this paper for combining information from multiple sources to support one-class classification. The contributing sources may represent measurements taken by different sensors of the same physical entity, repeated measurements by a single sensor, or numerous features computed from a single measured image or signal. The approach utilizes the theory of statistical hypothesis testing, and applies Fisher's technique for combining p-values, modified to handle nonindependent sources. Classifier outputs take the form of fused p-values, which may be used to gauge the consistency of unknown entities with one or more class hypotheses. The approach enables rigorous assessment of classification uncertainties, and allows for traceability of classifier decisions back to the constituent sources, both of which are important for high-consequence decision support. Application of the technique is illustrated in two challenge problems, one for skin segmentation and the other for terrain labeling. Finally, the method is seen to be particularly effective for relatively small training samples.

97 MATHEMATICS AND COMPUTING↗

Evaluation of precipitation indices in suites of dynamically and statistically downscaled regional climate models over Florida

Abstract The present work evaluates historical precipitation and its indices defined by the Expert Team on Climate Change Detection and Indices (ETCCDI) in suites of dynamically and statistically downscaled regional climate models (RCMs) against NOAA’s Global Historical Climatology Network Daily (GHCN-Daily) dataset over Florida. The models examined here are: (1) nested RCMs involved in the North American CORDEX (NA-CORDEX) program, (2) variable resolution Community Earth System Models (VR-CESM), (3) Coupled Model Intercomparison Project phase 5 (CMIP5) models statistically downscaled using localized constructed analogs (LOCA) technique. To quantify observational uncertainty, three in situ-based (PRISM, Livneh, CPC) and three reanalysis (ERA5, MERRA2, NARR) datasets are also evaluated against the station data. The reanalyses and dynamically downscaled RCMs generally underestimate the magnitude of the monthly precipitation and the frequency of the extreme rainfall in summer. The models forced with CanESM2 miss the phase of the seasonality of extreme precipitation. All models and reanalyses severely underestimate both the mean and interannual variability of mean wet-day precipitation (SDII), consecutive dry days (CDD), and overestimate consecutive wet days (CWD). Metric analysis suggests large uncertainty across NA-CORDEX models. Both the LOCA and VR-CESM models perform better than the majority of models. Overall, RegCM4 and WRF models perform poorer than the median model performance. The performance uncertainty across models is comparable to that in the reanalyses. Specifically, NARR performs poorer than the median model performance in simulating the mean indices and MERRA2 performs worse than the majority of models in capturing the interannual variability of the indices.

54 ENVIRONMENTAL SCIENCES↗

PYSIMFRAC: A Python library for synthetic fracture generation and analysis

In this paper, we introduce PYSIMFRAC, an open-source python library for generating 3-D synthetic fracture realizations, integrating with fluid simulators, and performing analysis. PYSIMFRAC allows the user to specify one of three fracture generation techniques (Box, Gaussian, or Spectral) and perform statistical analysis including the autocorrelation, moments, and probability density functions of the fracture surfaces and aperture. This analysis and accessibility of a python library allows the user to create realistic fracture realizations and vary properties of interest. In addition, PYSIMFRAC includes integration examples to two different pore-scale simulators and the discrete fracture network simulator, dfnWorks. The capabilities developed in this work provides opportunity for quick and smooth adoption and implementation by the wider scientific community for accurate characterization of fluid transport in geologic media. We present PYSIMFRAC along with integration examples and discuss the ability to extend PYSIMFRAC from a single complex fracture to complex fracture networks.

58 GEOSCIENCES↗

Directed Gas-Phase Formation of the 1-Cyanovinyl Radical (H 2 CCCN, X 2 A′) in the Interstellar Medium

The formation pathways to nitrogen-containing molecules and radicals are crucial to the understanding of the carbon–nitrogen chemistry in interstellar and atmospheric environments. While over 65 nitrogen-containing neutral species have been observed in deep space to date, their formation mechanisms─in particular, those of radical species─remain largely speculative. The crossed molecular beam technique in conjunction with electronic structure and statistical calculations was utilized to offer a detailed overview of the fundamental pathways in the gas-phase bimolecular reaction of ground-state atomic carbon (C, 3 P) with acetonitrile-d 3 (CD 3 CN, X 1 A 1 ) under single-collision conditions leading to the formation of the 1-cyanovinyl radical (D 2 CCCN, X 2 A′) coupled with deuterium atom loss. Here, the indirect reaction was initiated by barrierless carbon-atom addition, with the most probable route involving carbon addition across the carbon–nitrogen nitrile triple bond of acetonitrile, forming a three-membered ring intermediate followed by ring-opening and unimolecular decomposition via atomic deuterium loss from the C3 carbon atom. The reaction was overall exoergic, and intermediates and transition states lie lower in energy than the separated reactants, unlocking the reaction of carbon with acetonitrile in low-temperature environments such as cold molecular clouds, e.g., Taurus Molecular Cloud (TMC-1), and planetary atmospheres, e.g., Saturn’s moon Titan. In these environments, the 1-cyanovinyl radical may act as a building block for cyano-substituted polycyclic aromatic hydrocarbons and N-heterocycles, thus furthering our understanding of the complex carbon–nitrogen chemistry in deep space.

Chemical reactions↗