Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analysis of Variance”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Global Maps of Lunar Neutron Fluxes from the LEND Instrument

The latest neutron spectrometer measurements with the Lunar Exploration Neutron Detector (LEND) onboard the Lunar Reconnaissance Orbiter (LRO) are presented. It covers more than 1 year of mapping phase starting on 15 September 2009. In our analyses we have created global maps showing regional variations in the flux of thermal (energy range < 0.015 eV) and fast neutrons (>0.5 MeV), and compared these fluxes to variances in soil elemental composition, and with previous results obtained by the Lunar Prospector Neutron Spectrometer (LPNS). We also processed data from LEND collimated detectors and derived a value for the collimated signal of epithermal neutrons based on the comparative analysis with the LEND omnidirectional detectors. Finally, we have compared our final (after the data reduction) global epithermal neutron map with LPNS data.

Epithermal neutrons↗
DOE OSTI2018

Graphite Characterization: Baseline Variability Analysis Report

The Idaho National INL’s Baseline Material Property program provides many measurements of different material properties on a variety of nuclear graphite grades being investigated for irradiation-induced changes in material properties. The baseline material properties provide an unirradiated reference for comparison with measurements of those properties after irradiation. This study describes how uncertainty, or variance, in estimates of important material properties, like compressive strength, may depend on sample size, sample distribution and parameter estimation methods. Compressive strength data from two large collections of specimens from a billet of PCEA graphite (sample size = 230) and a billet of IG-110 graphite (sample size = 48) are used to illustrate these dependencies. The grades of graphite used represent approximate end members in terms of relative variability in strength measurements. Monte Carlo simulation and bootstrapping methods were used to examine variance in Weibull-parameter estimates as a function of sample size. Results were consistent with what previous statistical analyses have shown about the relationship between sample size and width of the confidence interval for maximum likelihood (ML) estimated parameters when samples are repeatedly drawn from a Weibull distribution. As a function of sample size, the precision of the modulus estimate increases faster than that of the characteristic strength. The modulus maximum likelihood estimates (MLEs) are biased for low sample sizes. The American Society for Testing and Materials document ASTM D7846-16 provides equations that describe the dependence of the variance of parameter MLEs on sample size. Using the same methodology, we show that the polynomials provided for variance dependence on sample size are inaccurate for sample sizes greater than ~120. We provide approximations to that relationship that illustrate relative constancy for sample size > ~10, in the logspace relationship of the modulus variance and the log transformed characteristic strength variance. Bootstrapping provided a means of comparing Weibull parameters obtained using ASME’s guidance for specimen collection to those obtained using different weighting in a case where strength is related to position in a billet and to grain orientation and a single distribution is used to characterize the billet. ASME specifies that samples be collected with equal representation in all locations and grain orientations. Parameter values obtained using ASME’s method were different from those obtained using random sampling of the complete data sets, and this was particularly true for the PCEA data. In addition, while the ASME approach yielded more conservative parameters for the PCEA data, the opposite was true for the IG-110 data. The difference observed for the PCEA data set likely reflects the larger number of against-grain specimens collected, which had generally higher strength than the with-grain specimens. This emphasizes the importance of the dependence of strength on location or grain orientation and suggests that it may be prudent to weight specimen collection toward the weaker subgroups if a conservative estimate as desired, as the equal weighting method recommended by ASME would will not always provide the most conservative Weibull parameters. Because spatial correlations commonly exist in certain grades of graphite, future work may involve developing a methodology for parameter inference that relaxes the assumption of independent specimens. ASTM and ASME codes recommend different parameter estimation methods and models of the Weibull distribution in different situations. Descriptions of the effects on parameter estimates for the example billets used in this study provide an example of how the method may affect the parameter estimate, and how—in this case—that compares to differences due to other factors, such as sample size and dependence on grain orientation or sample location.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Comparative Analysis of Thrombin Calibration Algorithms and Correction for Thrombin-α2macroglobulin Activity

Background: The thrombin generation (TG) test is useful for characterizing global hemostasis potential, but fluorescence substrate artifacts, such as thrombin-α2macroglobulin (T-α2MG) signal, inner filter effect (IFE), substrate consumption, and calibration algorithms have been suggested as sources of intra- and inter-laboratory variance, which may limit its clinical utility. Methods: Effects of internal vs. external normalization, IFE and T-α2MG on TG curves in normal plasma supplemented with coagulation factors, thrombomodulin, and tissue factor were studied using the Calibrated Automated Thrombinography (CAT; Diagnostica Stago, Parsippany, NJ, USA) and in-house software. Results: The various calibration methods demonstrated no significant difference in producing TG curves, nor increased the robustness of the TG assay. Several TG parameters, including thrombin peak height (TPH), produced from internal linear calibration did not differ significantly from uncalibrated TG parameters. Further, TPH values from internal linear and nonlinear calibration with or without T-α2MG correction correlated well with TPH from external calibration. Higher coefficients of variation (CVs) for TPH values were observed in both platelet-free and platelet-rich plasma with added thrombomodulin. Conclusions: Our work suggests minimal differences between distinct computational approaches toward calibrating and correcting fluorescence signals into TG levels, with most samples returning similar or equivalent TPH results.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1975

M-071 critical data analysis

A prototype balance study was conducted on earth prior to the balance studies conducted in Skylab itself. Collected were daily dietary intake data of 6 minerals and nitrogen, and fecal and urinary outputs on each of three astronauts. Essential statistical issues show what quantities need to be estimated and establish the scope of inference associated with alternative variance estimates. The procedures for obtaining the final variability due both to errors of measurement and total error (total = measurement and biological variability) are exhibited.

Hegsted, D. M.↗
NASA NTRS1991

Analysis of upper and lower bounds of the frame noise in linear detector arrays

This paper estimates the upper and lower bounds of the frame noise of a linear detector array that uses a one-dimensional scan pattern. Using chi-square distribution, it is analytically shown why it is necessary to use the average of the variances and not the average of the standard deviations to estimate these bounds. Also, a criteria for determining whether any excessively noisy lines exist among the detectors is derived from these bounds. Using a Gaussian standard random variable generator, these bounds are demonstrated to be accurate within the specified confidence interval. A silicon detector array is then used for actual dark current measurements. The criterion developed for determination of noisy detectors is checked on the experimentally obtained data.

Jaggi, S.↗
DOE OSTI2014

A stochastic analysis of steady and transient heat conduction in random media using a homogenization approach

We present a new stochastic analysis for steady and transient one-dimensional heat conduction problem based on the homogenization approach. Thermal conductivity is assumed to be a random field K consisting of random variables of a total number N. Both steady and transient solutions T are expressed in terms of the homogenized solution (symbol) and its spatial derivatives (equation), where homogenized solution (symbol) is obtained by solving the homogenized equation with effective thermal conductivity. Both mean and variance of stochastic solutions can be obtained analytically for K field consisting of independent identically distributed (i.i.d) random variables. The mean and variance of T are shown to be dependent only on the mean and variance of these i.i.d variables, not the particular form of probability distribution function of i.i.d variables. Variance of temperature field T can be separated into two contributions: the ensemble contribution (through the homogenized temperature (symbol)); and the configurational contribution (through the random variable Ln(x)Ln(x)). The configurational contribution is shown to be proportional to the local gradient of (symbol). Large uncertainty of T field was found at locations with large gradient of (symbol) due to the significant configurational contributions at these locations. Numerical simulations were implemented based on a direct Monte Carlo method and good agreement is obtained between numerical Monte Carlo results and the proposed stochastic analysis.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2016

Mechanical vulnerability and resistance to snapping and uprooting for Central Amazon tree species

We report high descending winds generated by convective storms are a frequent and a major source of tree mortality disturbance events in the Amazon, affecting forest structure and diversity across a variety of scales, and more frequently observed in western and central portions of the basin. Soil texture in the Central Amazon also varies significantly with elevation along a topographic gradient, with decreasing clay content on plateaus, slopes and valleys respectively. In this study we investigated the critical turning moments ( M crit - rotational force at the moment of tree failure, an indicator of tree stability or wind resistance) of 60 trees, ranging from 19.0 to 41.1 cm in diameter at breast height ( DBH ) and located in different topographic positions, and for different species, using a cable-winch load-cell system. Our approach used torque as a measure of tree failure to the point of snapping or uprooting. This approach provides a better understanding of the mechanical forces required to topple trees in tropical forests, and will inform models of wind throw disturbance. Across the topographic positions, size controlled variation in M crit was quantified for cardeiro ( Scleronema mincranthum (Ducke) Ducke), mata-matá ( Eschweilera spp.), and a random selection of trees from 19 other species. Our analysis of M crit revealed that tree resistance to failure increased with size (DBH and ABG) and differed among species. No effects of topography or failure mode were found for the species either separately or pooled. For the random species, total variance in M crit explained by tree size metrics increased from an R 2 of 0.49 for DBH alone, to 0.68 when both DBH and stem fresh wood density (SWD) were included in a multiple regression model. This mechanistic approach allows the comparison of tree vulnerability induced by wind damage across ecosystems, and facilitates the use of forest structural information in ecosystem models that include variable resistance of trees to mortality inducing factors. Furthermore, our results indicate that observed topographic differences in wind throw vulnerability are likely due to elevational differences in wind velocities, rather than by differences in soil-related factors that might effect M crit .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Graph Sparsification by Approximate matrix Multiplication

Graphs arising in statistical problems, signal processing, large networks, combinatorial optimization, and data analysis are often dense, which causes both computational and storage bottlenecks. One way of sparsifying a weighted graph, while sharing the same vertices as the original graph but reducing the number of edges, is through spectral sparsification. We study this problem through the perspective of RandNLA. Specifically, we utilize randomized matrix multiplication to give a clean and simple analysis of how sampling according to edge weights gives a spectral approximation to graph Laplacians, without requiring spectral information. Through the CR–MM algorithm, we attain a simple and computationally efficient sparsifier whose resulting Laplacian estimate is unbiased and of minimum variance. Here, we define a new notion of additive spectral sparsifiers, which has not been considered in the literature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Uncertainty analysis for infrasound waveform inversion: Application to explosion yield estimation

While the acoustic waveform inversion method is increasingly used in geophysical acoustics to constrain source parameters, the inversion results are often provided without any uncertainty analysis. Here, this study presents a probabilistic representation for acoustic waveform inversion and method to evaluate the inversion uncertainty using ground-truth data. A posteriori probability distribution of source estimate is described by a priori waveform misfit covariance and the variance of acoustic source model. The probabilistic framework is applied to local explosion infrasound to estimate the yields of explosions and uncertainty. Estimated yields showed overall good agreement with the true yields (less than 25% errors). The uncertainty of the estimated yield is represented by the sum of the waveform inversion uncertainty and source model uncertainty. It is shown that the yield uncertainty attributed to local infrasound inversion (within 10 km) is as small as the uncertainty caused by 10% prediction errors in the acoustic source model. Finally, these results indicate that the acoustic source model uncertainty should also be considered for accurate yield estimation and that local infrasound can be a valuable tool to understand the magnitude of the source uncertainty.

58 GEOSCIENCES↗
NASA NTRS1981

The current status of REH theory

A response is made to the evaluation of Fitch (1980) of REH (random evolutionary hits) theory for the evolutionary divergence of proteins and nucleic acids. Correct calculations for the beta hemoglobin mRNAs of the human, mouse and rabbit in the absence and presence of selective constraints are summarized, and it is shown that the alternative evolutionary analysis of Fitch underestimates the total fixed mutations. It is further shown that the model used by Fitch to test for the completeness of the count of total base substitutions is in fact a variant of REH theory. Considerations of the variance inherent in evolutionary estimations are also presented which show the REH model to produce no more variance than other evolutionary models. In the reply, it is argued that, despite the objections raised, REH theory applied to proteins gives inaccurate estimates of total gene substitutions. It is further contended that REH theory developed for nucleic sequences suffers from problems relating to the frequency of nucleotide substitutions, the identity of the codons accepting silent and amino acid-changing substitutions, and estimate uncertainties.

Holmquist, R.↗
NASA NTRS1986

Rendezvous BET Program

Computes relative positions of two vehicles in concentric orbits. LRBET3 program best-estimate-of-trajectory (BET) calculation for postflight trajectory analysis of Shuttle orbital rendezvous maneuvers. LRBET3 produces estimated measurements for reconstructing relative positions of two vehicles. Kalman filter and smoothing filter applied to relative measurement input data to estimate state vector, reduce noise, and produce BET output. BET calculation minimizes variances of all trajectory estimation errors. LRBET3 written in FORTRAN IV for batch execution.

Lear, W. M.↗
NASA NTRS1986

Magnetic field directional discontinuities - Characteristics between 0.46 and 1.0 AU

Based on Mariner 10 data, a statistical survey and an application of the Sonnerup-Cahill variance procedure to a visual identification with 1.2-s averages for time intervals corresponding to the equally spaced heliocentric distances of 1.0, 0.72 and 0.46 AU, are employed to study the characteristics of directional discontinuities (DDs) in the interplanetary magnetic field. Analysis using two methods demonstrated that the ratio of tangential discontinuities (TDs) to rotational discontinuities (RDs) decreased with decreasing radial distance. Decreases in average discontinuity thickness of 41 percent between 1.0 and 0.72 AU, and 56 percent between 1.0 and 0.46 AU, were found for both TDs and RDs, in agreement with Pioneer 10 data between 1 and 5 AU. Normalization of the individual DD thicknesses with respect to the estimated local proton gyroradius (R sub L) gave a nearly constant average thickness at the three locations, 36 + or - 5 R sub L, for both RDs and TDs.

Lepping, R. P.↗
NASA NTRS2000

Comparison of On-Orbit Performance of Rate Sensing Gyroscopes

This work presents results from the study of a large volume of spacecraft flight data pertaining to gyroscope performance. We have examined long and short term trends of gyroscope biase of first, an exact solution for the time-dependence of the attitude part of the state error covariance, averaged over the three spacecraft axes, is given. This solution is more complete than the usual cubic polynomial for the variance in that it includes coupling among the three gyroscope axes and includes an important term arising from the initial correlations. Second, several continuous 24-hour spans of gyroscope data are examined to verify the short-term statistical model. This analysis demonstrates that in-flight data can be used to determine the strength of the white noise driving the random walk of the gyroscope bias. This may be useful for postlaunch improvement to the noise model and for diagnosing the health of the gyroscope. Third, the long-term trends in gyroscope biases show a nearly linear systematic variation over time scales of years. This has been found on three different missions. While the random walk model is adequate as a basis for onboard Kalman filters or for state estimation using relatively short time spans, these trends indicate that some applications could benefit by accounting for the secular changes in the biases. One example is a new gyroscope calibration method that is under development that allows for multi-epoch bias solutions.

Sedlak, Joseph↗
NASA NTRS2005

Global Precipitation Variations and Long-term Changes Derived from the GPCP Monthly Product

Global and large regional rainfall variations and possible long-term changes are examined using the 25-year (1979-2004) monthly dataset from the Global Precipitation Climatology Project (GPCP). The emphasis is to discriminate among the variations due to ENSO, volcanic events and possible long-term changes. Although the global change of precipitation in the data set is near zero, the data set does indicate an upward trend (0.13 mm/day/25yr) and a downward trend (-0.06 mm/day/25yr) over tropical oceans and lands (25S-25N), respectively. This corresponds to a 4% increase (ocean) and 2% decrease (land) during this time period. Techniques are applied to attempt to eliminate variations due to ENSO and major volcanic eruptions. The impact of the two major volcanic eruptions over the past 25 years is estimated to be about a 5% reduction in tropical rainfall. The modified data set (with ENSO and volcano effect removed) retains the same approximate change slopes, but with reduced variance leading to significance tests with results in the 90-95% range. Inter-comparisons between the GPCP, SSWI (1988-2004), and TRMM (1998-2004) rainfall products are made to increase or decrease confidence in the changes seen in the GPCP analysis.

Adler, Robert F.↗
NASA NTRS2020

MSL Telecom Automated Anomaly Detection

The Mars Science Laboratory (MSL) Telecom Operations Team at the Jet Propulsion Laboratory (JPL) has implemented a machine learning system in order to automate the anomaly detection process as a part of daily operations. Machine learning enables reliable detection of anomalies in Telecom-related telemetry and automated reporting of Telecom subsystem status, resulting in an 90% reduction in team workload and improved anomaly detection reliability. At present, machine learning methods are used to detect: 1. Anomalous long-term trends in telemetry data 2. Anomalous time-domain evolution of telemetry values Both types of anomalies pose their own unique challenges that are addressed in different ways. In the first case, long term trending of daily minima, maximum, and mean telemetry values in temperatures, currents, voltages, and radio frequency (RF) power levels is used in addition to hard threshold safety checks to look for changes in long-term equipment health and performance. Long-term trending methods allow for ordinary seasonal variations in these quantities caused by temperature changes over the course of the Martian year while allowing operators to determine whether current performance remains in line with historical values from previous years. Changes in long-term trends can provide important insights into the health and status of the rover's on-board systems as well as valuable early warning if subtle degradation begins to take hold. But while trending of daily statistics is valuable, it does not detect anomalies in the short-term time evolution of data over the course of minutes or hours during a day, and this task is handled with short-term shape analysis. Principal components analysis (PCA) has been found to provide robust detection of short-term anomalies, and several examples of the use of PCA to detect actual anomalous events will be provided here. In using PCA, we use both the percentage of explained variance and also a log likelihood test on the PCA expansion coefficients to flag telemetry data for human review. Previous work in the field of spacecraft anomaly detection includes [1] for MSL and [2] for some other JPL missions.

Mukai, Ryan↗
DOE OSTI2015

Latent uncertainties of the precalculated track Monte Carlo method

Purpose: While significant progress has been made in speeding up Monte Carlo (MC) dose calculation methods, they remain too time-consuming for the purpose of inverse planning. To achieve clinically usable calculation speeds, a precalculated Monte Carlo (PMC) algorithm for proton and electron transport was developed to run on graphics processing units (GPUs). The algorithm utilizes pregenerated particle track data from conventional MC codes for different materials such as water, bone, and lung to produce dose distributions in voxelized phantoms. While PMC methods have been described in the past, an explicit quantification of the latent uncertainty arising from the limited number of unique tracks in the pregenerated track bank is missing from the paper. With a proper uncertainty analysis, an optimal number of tracks in the pregenerated track bank can be selected for a desired dose calculation uncertainty. Methods: Particle tracks were pregenerated for electrons and protons using EGSnrc and GEANT4 and saved in a database. The PMC algorithm for track selection, rotation, and transport was implemented on the Compute Unified Device Architecture (CUDA) 4.0 programming framework. PMC dose distributions were calculated in a variety of media and compared to benchmark dose distributions simulated from the corresponding general-purpose MC codes in the same conditions. A latent uncertainty metric was defined and analysis was performed by varying the pregenerated track bank size and the number of simulated primary particle histories and comparing dose values to a “ground truth” benchmark dose distribution calculated to 0.04% average uncertainty in voxels with dose greater than 20% of D{sub max}. Efficiency metrics were calculated against benchmark MC codes on a single CPU core with no variance reduction. Results: Dose distributions generated using PMC and benchmark MC codes were compared and found to be within 2% of each other in voxels with dose values greater than 20% of the maximum dose. In proton calculations, a small (≤1 mm) distance-to-agreement error was observed at the Bragg peak. Latent uncertainty was characterized for electrons and found to follow a Poisson distribution with the number of unique tracks per energy. A track bank of 12 energies and 60000 unique tracks per pregenerated energy in water had a size of 2.4 GB and achieved a latent uncertainty of approximately 1% at an optimal efficiency gain over DOSXYZnrc. Larger track banks produced a lower latent uncertainty at the cost of increased memory consumption. Using an NVIDIA GTX 590, efficiency analysis showed a 807 × efficiency increase over DOSXYZnrc for 16 MeV electrons in water and 508 × for 16 MeV electrons in bone. Conclusions: The PMC method can calculate dose distributions for electrons and protons to a statistical uncertainty of 1% with a large efficiency gain over conventional MC codes. Before performing clinical dose calculations, models to calculate dose contributions from uncharged particles must be implemented. Following the successful implementation of these models, the PMC method will be evaluated as a candidate for inverse planning of modulated electron radiation therapy and scanned proton beams.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Measurements of the E-mode polarization and temperature-E-mode correlation of the CMB from SPT-3G 2018 data

We present measurements of the E -mode ( E E ) polarization power spectrum and temperature- E -mode ( T E ) cross-power spectrum of the cosmic microwave background using data collected by SPT-3G, the latest instrument installed on the South Pole Telescope. This analysis uses observations of a 1500 deg 2 region at 95, 150, and 220 GHz taken over a four-month period in 2018. We report binned values of the E E and T E power spectra over the angular multipole range 300 ≤ ℓ < 3000 , using the multifrequency data to construct six semi-independent estimates of each power spectrum and their minimum-variance combination. These measurements improve upon the previous results of SPTpol across the multipole ranges 300 ≤ ℓ ≤ 1400 for E E and 300 ≤ ℓ ≤ 1700 for T E , resulting in constraints on cosmological parameters comparable to those from other current leading ground-based experiments. We find that the SPT-3G data set is well fit by a Λ CDM cosmological model with parameter constraints consistent with those from Planck and SPTpol data. From SPT-3G data alone, we find H 0 = 68.8 ± 1.5 km s - 1 Mpc - 1 and σ 8 = 0.789 ± 0.016 , with a gravitational lensing amplitude consistent with the Λ CDM prediction ( A L = 0.98 ± 0.12 ). We combine the SPT-3G and the Planck data sets and obtain joint constraints on the Λ CDM model. The volume of the 68% confidence region in six-dimensional Λ CDM parameter space is reduced by a factor of 1.5 compared to Planck-only constraints, with no significant shifts in central values. We note that the results presented here are obtained from data collected during just half of a typical observing season with only part of the focal plane operable, and that the active detector count has since nearly doubled for observations made with SPT-3G after 2018.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Analysis of core asymmetries in inertial confinement fusion implosions using three-dimensional hot-spot reconstruction

Three-dimensional effects play a crucial role during the hot-spot formation in inertial confinement fusion (ICF) implosions. A data analysis technique for 3D hot-spot reconstruction from experimental observables has been developed to characterize the effects of low modes on 3D hot-spot formations. In nuclear measurements, the effective flow direction, governed by the maximum eigenvalue in the velocity variance of apparent ion temperatures, has been found to agree with the measured hot-spot flows for implosions dominated by mode ℓ = 1. Asymmetries in areal-density (ρR) measurements were found to be characterized by a unique cosine variation along the hot-spot flow axis. In x-ray images, a 3D hot-spot x-ray emission tomography method was developed to reconstruct the 3D hot-spot plasma emissivity using a generalized spherical-harmonic Gaussian function. The gradient-descent algorithm was used to optimize the mapping between the projections from the 3D hot-spot emission model and the measured x-ray images along multiple views. Furthermore, this work establishes a platform to analyze 3D low-mode core asymmetries in ICF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗