Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling statistics”

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,459 records · Page 81

Photometric study on stellar magnetic activity. I. Flare variability of red dwarf stars in the open cluster M37

Based on one-month long MMT time-series observations of the open cluster M37, we monitored light variations of nearly 2500 red dwarfs and successfully identified 420 flare events from 312 cluster M dwarf stars. For each flare light curve, we derived observational and physical parameters, such as flare shape, peak amplitude, duration, energy, and peak luminosity. We show that cool stars produce serendipitous flares energetic enough to be observed in the r-band, and their temporal and peak characteristics are almost the same as those in traditional U-band observations. We also found many large-amplitude flares with inferred Δu>6 mag in the cluster sample which had been rarely reported in previous ground-based observations. Following the ergodic hypothesis, we investigate in detail statistical properties of flare parameters over a range of energy (E{sub r} ≃ 10{sup 31}–10{sup 34} erg). As expected, there are no statistical differences in the distributions of flare timescales, energies, and frequencies among stars of the same age and mass group. We note that our sample tend to have longer rise and decay timescales compared to those seen in field flare stars of the same spectral type and be more energetic. Flare frequency distributions follow power-law distributions with slopes β∼0.62–1.21 for all flare stars and β∼0.52–0.97 for stars with membership information (P{sub mem}⩾0.2). These are in general agreement with previous works on flare statistics of young open clusters and nearby field stars. Our results give further support to the classical age–activity relations.

79 ASTRONOMY AND ASTROPHYSICS↗

SpotSDC: Revealing the Silent Data Corruption Propagation in High-Performance Computing Systems

We report the trend of rapid technology scaling is expected to make the hardware of high-performance computing (HPC) systems more susceptible to computational errors due to random bit flips. Some bit flips may cause a program to crash or have a minimal effect on the output, but others may lead to silent data corruption (SDC), i.e., undetected yet significant output errors. Classical fault injection analysis methods employ uniform sampling of random bit flips during program execution to derive a statistical resiliency profile. However, summarizing such fault injection result with sufficient detail is difficult, and understanding the behavior of the fault-corrupted program is still a challenge. In this article, we introduce SpotSDC, a visualization system to facilitate the analysis of a program's resilience to SDC. SpotSDC provides multiple perspectives at various levels of detail of the impact on the output relative to where in the source code the flipped bit occurs, which bit is flipped, and when during the execution it happens. SpotSDC also enables users to study the code protection and provide new insights to understand the behavior of a fault-injected program. Based on lessons learned, we demonstrate how what we found can improve the fault injection campaign method.

97 MATHEMATICS AND COMPUTING↗

Recombination smooths the time-signal disrupted by latency in within-host HIV phylogenies

Within-host HIV evolution involves several features that may disrupt standard phylogenetic reconstruction. One important feature is re-activation of latently integrated provirus, which has the potential to disrupt the temporal signal, leading to variation in the branch lengths and apparent evolutionary rates in a tree. Yet, real within-host HIV phylogenies tend to show clear, ladder-like trees structured by the time of sampling. Another important feature is recombination, which violates the fundamental assumption that evolutionary history can be represented by a single bifurcating tree. Thus, recombination complicates the within-host HIV dynamic by mixing genomes and creating evolutionary loop structures that cannot be represented in bifurcating trees. In this paper, we develop a coalescent-based simulator of within-host HIV evolution that includes latency, recombination, and effective population size dynamics that allows us to study the relationship between the true, complex genealogy of within-host HIV evolution, encoded as an Ancestral Recombination Graph (ARG), and the observed phylogenetic tree. To compare our ARG results to the familiar phylogeny format, we calculate the expected bifurcating tree after decomposing the ARG into all unique site trees, their combined distance matrix, and the overall corresponding bifurcating tree. While latency and recombination separately disrupt the phylogenetic signal, remarkably, we find that recombination recovers the temporal signal of within-host HIV evolution caused by latency by mixing fragments of old, latent genomes into the contemporary population. In effect, recombination averages over extant heterogeneity, whether it stems from mixed time-signals or population bottlenecks. Further, we establish that the signals of latency and recombination can be observed in phylogenetic trees despite being an incorrect representation of the true evolutionary history. Using an Approximate Bayesian Computation method, we develop a set of statistical probes to tune our simulation model to nine longitudinally-sampled within-host HIV phylogenies. Because ARGs are exceedingly difficult to infer from real HIV data, our simulation system allows investigating effects of latency, recombination, and population size bottlenecks by matching decomposed ARGs to real data as observed in standard phylogenies.

59 BASIC BIOLOGICAL SCIENCES↗

Statistical Analysis of Supernova Remnants in the Large Magellanic Cloud

We construct the most complete sample of supernova remnants (SNRs) in any galaxy—the Large Magellanic Cloud (LMC) SNR sample. We study their various properties such as spectral index (α), size, and surface brightness. We suggest an association between the spatial distribution and environment density of LMC SNRs, and their tendency to be located around supergiant shells. We find evidence that the 16 known type Ia LMC SNRs are expanding in a lower density environment compared to the Core-Collapse (CC) type. The mean diameter of our entire population (74) is 41 pc, which is comparable to nearby galaxies. We did not find any correlation between the type of SN explosion, ovality, or age. The N(<D) relationship with slope a = 0.96 implies that the randomized diameters are readily mimicking such an exponent. The rate of SNe occurring in the LMC is estimated to be ∼1 per 200 yr. The mean α of the entire LMC SNR population is −0.52, which is typical of most SNRs. However, our estimates show a clear flattening of the synchrotron α as the remnants age. As predicted, the CC SNRs in our sample are significantly brighter radio emitters than type Ia remnants. We also estimate the Σ--D relation for the LMC to have a slope ∼3.8, which is comparable with other nearby galaxies. We also find the residency time of electrons in the galaxy (4.0–14.3 Myr), implying that SNRs should be the dominant mechanism for the production and acceleration of CRs.

79 ASTRONOMY AND ASTROPHYSICS↗

Random Dopant Induced Threshold Voltage Lowering and Fluctuations in Sub 50 nm MOSFETs: a Statistical 3D 'Atomistic' Simulation Study

A 3D 'atomistic' simulation study of random dopant induced threshold voltage fluctuations and lowering in sub 50 nm MOSFETs is presented. The attention is focused mainly on devices with 30 nm effective channel length which represent the expected level of scaling at the end of the Silicon Roadmap. An efficient algorithm, based on a single 3D ap solution of the Poisson equation and a simplified current continuity equation, is used in the simulations. Large samples of microscopically different devices (typically 200) arc used in order to obtain statistically reliable results. The influence of different aspects of the conventional MOSFET design on the threshold voltage fluctuations and lowering are investigated. Results for fluctuation resistant device architectures based on low-doped epitaxial channel MOSFETs are also presented.

Asenov, Asen↗

Improved Sampling Algorithms in the Risk-Informed Safety Margin Characterization Toolkit

The RISMC approach is developing advanced set of methodologies and algorithms in order to perform Probabilistic Risk Analyses (PRAs). In contrast to classical PRA methods, which are based on Event-Tree and Fault-Tree methods, the RISMC approach largely employs system simulator codes applied to stochastic analysis tools. The basic idea is to randomly perturb (by employing sampling algorithms) timing and sequencing of events and internal parameters of the system codes (i.e., uncertain parameters) in order to estimate stochastic parameters such as core damage probability. This approach applied to complex systems such as nuclear power plants requires to perform a series of computationally expensive simulation runs given a large set of uncertain parameters. These types of analysis are affected by two issues. Firstly, the space of the possible solutions (a.k.a., the issue space or the response surface) can be sampled only very sparsely, and this precludes the ability to fully analyze the impact of uncertainties on the system dynamics. Secondly, large amounts of data are generated and tools to generate knowledge from such data sets are not yet available. This report focuses on the first issue and in particular employs novel methods that optimize the information generated by the sampling process by sampling unexplored and risk-significant regions of the issue space: adaptive (smart) sampling algorithms. They infer system response from surrogate models constructed from existing samples and predict the most relevant location of the next sample. It is therefore possible to understand features of the issue space with a small number of carefully selected samples. In this report, we will present how it is possible to perform adaptive sampling using the RISMC toolkit and highlight the advantages compared to more classical sampling approaches such Monte-Carlo. We will employ RAVEN to perform such statistical analyses using both analytical cases but also another RISMC code: RELAP-7.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Southeast Carbon Sequestration Partnership (SECARB) Phase III Anthropogenic Test at Citronelle Field, Mobile Co., AL Surface and Shallow Subsurface Monitoring

The Southeast Carbon Sequestration Partnership Phase III Anthropogenic Test at Citronelle Field is one of the United States’ largest demonstrations of CO 2 capture, transportation, injection, storage and monitoring using CO 2 derived from a coal-fired electric generating unit. This report summarizes the background and initial injection phase groundwater and soil flux monitoring results. The purpose of the groundwater monitoring program, as prescribed in the Underground Injection Control (UIC) permits is to assure that the project’s activities do not harm Underground Sources of Drinking Water (USDWs). In particular, the objective of the program is to monitor for the release of injected CO2 into groundwater aquifers at the project site. Three shallow USDW monitoring wells were installed as part of the groundwater monitoring program. In addition, a water supply well located near the project’s deep characterization well was utilized for monitoring. Seventeen metals, general chemistry and field parameters (e.g. pH, specific conductance, turbidity) were sampled three times prior to injection and quarterly during the injection period to date. Changes in groundwater geochemistry potentially attributable to a release of carbon dioxide into the USDW from the injection activities are decreasing pH and increasing TIC. The change in pH may also accompany the mobilization of certain metals within the groundwater, but the changes are likely site-specific, based on the natural aquifer mineralogy (Wilkin and Diguilio, 2010). At this time, no downward trends for pH and increasing trends for TIC have been statistically determined in the monitor wells. Valid detections exceeding both the UIC discharge limit and the range of background concentrations observed as of the N=11 event include: Exceedance of the UIC permit discharge limit of manganese at D-9-7 MW-2S and D-9-9 MW-1, and exceedance of the UIC permit discharge limit of iron at D-9-9 MW-1. In addition, a statistically significant increase of barium and manganese has been observed at D-9-9 MW-1. However, Additional lines of evidence (namely, expected shifts in TIC or pH, per the research cited above) to interpret the cause (e.g., potential carbon dioxide influence) of the statistically determined exceedances and/or increases are absent at this time. Quarterly groundwater sampling and analysis will continue throughout the project. The UIC permit requires that soil CO 2 flux shall be monitored at the land surface at the injection site and at nearby locations. Twelve locations were chosen for surveys including a mix of near wellbore and offset locations. To date, 24 surveys have taken place with 12 occurring prior to the onset of injection operations in August 2012. Soil CO 2 flux appears to positively covary with soil temperature. A notably high soil CO 2 flux value (21 µmol/m 2 /sec) was observed from the D-9-9 #2 chamber location during sampling event #24 (9/24/2014). As a result of this finding, the frequency of soil CO 2 flux monitoring events will be increased to twice-monthly until this anomaly can be confirmed or determined to be an outlier. Soil CO 2 flux monitoring will continue for the duration of the project.

01 COAL, LIGNITE, AND PEAT↗

Using yttrium as an indicator to estimate total rare earth element concentration: a case study of anthracite-associated clays from northeastern Pennsylvania

Abstract This study demonstrated using yttrium (Y) as an indicator to estimate the total rare earth element and Y contents (REY) in coal-associated samples and to facilitate selection of samples with high REY assays in a fast and inexpensive manner. More than 10 anthracite-associated samples were collected from each of three Pennsylvanian sites (sites B, J and C) based on Thorium gamma ray logging suggesting high REY content. Several samples from each site were analyzed by ICP-MS to determine the rare earth distribution patterns and to establish the site-specific linear equations of Y and REY. The Y contents of the remaining samples were measured by a portable X-ray fluorescence analyzer, and the REY values were estimated based on the site-specific linear equation developed earlier. R-squared values above 0.70 were obtained for all the estimation equations from all three sites on both a whole sample basis and an ash basis. Previously, ash content has been widely used as an indicator of high REY content. This may not be applicable for a specific site. Site B in this study is an example where ash contents could not be statistically correlated with REY, so using Y for estimation is more applicable. The demonstrated sample screening process is suitable for samples from sites that share more similar distribution patterns (either MREY or LREY or HREY) as well as for samples from sites that share multiple distribution patterns (LREY/MREY/HREY) depending on the desirable accuracy. The demonstrated process lowers the analytical cost from $70 to 80 dollars per sample to $10–15 per sample while significantly reducing the processing time and acid consumption for ICP digestion. This is particularly true when a relatively large sample size is involved, for example, 100 samples from one site analyzed by ICP-MS/OES.

Yang, Xiaojing↗

A novel variant of torque teno virus 7 identified in patients with Kawasaki disease

Kawasaki disease (KD), first identified in 1967, is a pediatric vasculitis of unknown etiology that has an increasing incidence in Japan and many other countries. KD can cause coronary artery aneurysms. Its epidemiological characteristics, such as seasonality and clinical picture of acute systemic inflammation with prodromal intestinal/respiratory symptoms, suggest an infectious etiology for KD. Interestingly, multiple host genotypes have been identified as predisposing factors for KD. To explore experimental methodology for identifying etiological agent(s) for KD and to optimize epidemiological study design (particularly the sample size) for future studies, we conducted a pilot study. For a 1-year period, we prospectively enrolled 11 patients with KD. To each KD patient, we assigned two control individuals (one with diarrhea and the other with respiratory infections), matched for age, sex, and season of diagnosis. During the acute phase of disease, we collected peripheral blood, nasopharyngeal aspirate, and feces. We also determined genotypes, to identify those that confer susceptibility to KD. There was no statistically significant difference in the frequency of the risk genotypes between KD patients and control subjects. We also used unbiased metagenomic sequencing to analyze these samples. Metagenomic sequencing and PCR detected torque teno virus 7 (TTV7) in two patients with KD (18%), but not in control subjects (P = 0.111). Sanger sequencing revealed that the TTV7 found in the two KD patients contained almost identical variants in nucleotide and identical changes in resulting amino acid, relative to the reference sequence. Additionally, we estimated the sample size that would be required to demonstrate a statistical correlation between TTV7 and KD. Future larger scale studies with carefully optimized metagenomic sequencing experiments and adequate sample size are warranted to further examine the association between KD and potential pathogens, including TTV7.

59 BASIC BIOLOGICAL SCIENCES↗

TURBULENCE AND STAR FORMATION IN A SAMPLE OF SPIRAL GALAXIES

We investigate turbulent gas motions in spiral galaxies and their importance to star formation in far outer disks, where the column density is typically far below the critical value for spontaneous gravitational collapse. Following the methods of Burkhart et al. on the Small Magellanic Cloud, we use the third and fourth statistical moments, as indicators of structures caused by turbulence, to examine the neutral hydrogen (H i) column density of a sample of spiral galaxies selected from The H i Nearby Galaxy Survey. We apply the statistical moments in three different methods—the galaxy as a whole, divided into a function of radii and then into grids. We create individual grid maps of kurtosis for each galaxy. To investigate the relation between these moments and star formation, we compare these maps with their far-ultraviolet images taken by the Galaxy Evolution Explorer satellite.We find that the moments are largely uniform across the galaxies, in which the variation does not appear to trace any star-forming regions. This may, however, be due to the spatial resolution of our analysis, which could potentially limit the scale of turbulent motions that we are sensitive to greater than ∼700 pc. From comparison between the moments themselves, we find that the gas motions in our sampled galaxies are largely supersonic. This analysis also shows that the Burkhart et al. methods may be applied not just to dwarf galaxies but also to normal spiral galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Batch Effect Correction Methods for NASA GeneLab Transcriptomic Datasets

RNA sequencing (RNA-seq) data from space biology experiments promise to yield invaluable insights into the effects of spaceflight on terrestrial biology. However, sample numbers from each study are low due to limited crew availability, hardware, and space. To increase statistical power, spaceflight RNA-seq datasets from different missions are often aggregated together. However, this can introduce technical variation or "batch effects", often due to differences in sample handling, sample processing, and sequencing platforms. Several computational methods have been developed to correct for technical batch effects, thereby reducing their impact on true biological signals. In this study, we combined 7 mouse liver RNA-seq datasets from NASA GeneLab (part of the NASA Open Science Data Repository) to evaluate several common batch effect correction methods (ComBat and ComBat-seq from the sva R package, and Median Polish, Empirical Bayes, and ANOVA from the MBatch R package). We quantitatively evaluated the ability of these methods to correct for technical batch variables in space biology RNA-seq data using the following criteria: BatchQC, principal component analysis, dispersion separability criterion, log fold change correlation, and differential gene expression analysis. Each batch variable / correction method combination was then assessed using a custom scoring approach to identify the optimal correction method for the combined dataset, by geometrically probing the space of all allowable scoring functions to yield an aggregate volume-based scoring measure. Finally, we describe the way in which the GeneLab multi-study analysis and visualization portal will allow users to examine the presence or absence of batch effects using multiple metrics. If the user chooses to perform batch effect correction, the scoring approach described here can be implemented to identify the optimal correction method to use for their specific combined dataset prior to analysis.

Lauren M. Sanders↗

Ceci N'est Pas a globular cluster: the metallicity distribution of the stellar system Terzan 5

We present new determinations of the iron abundance for 220 stars belonging to the stellar system Terzan 5 in the Galactic bulge. The spectra have been acquired with FLAMES at the Very Large Telescope of the European Southern Observatory and DEIMOS at the Keck II Telescope. This is by far the largest spectroscopic sample of stars ever observed in this stellar system. From this data set, a subsample of targets with spectra unaffected by TiO bands was extracted and statistically decontaminated from field stars. Once combined with 34 additional stars previously published by our group, a total sample of 135 member stars covering the entire radial extent of the system has been used to determine the metallicity distribution function of Terzan 5. The iron distribution clearly shows three peaks: a super-solar component at [Fe/H] ≅ 0.25 dex, accounting for ∼29% of the sample, a dominant sub-solar population at [Fe/H] ≅ –0.30 dex, corresponding to ∼62% of the total, and a minor (6%) metal-poor component at [Fe/H] ≅ –0.8 dex. Such a broad, multi-modal metallicity distribution demonstrates that Terzan 5 is not a genuine globular cluster but the remnant of a much more complex stellar system.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistical fracture behavior of doped UO 2 using a ball-on-ring equibiaxial flexure test method

Metal oxide dopants, such as titanium and chromium oxides, have garnered considerable attention for their potential to increase grain size (≥ 30 µm) in UO 2 fuel, purportedly enhancing fission gas retention during reactor operation. Fuel performance is significantly impacted by fuel fracture behavior, so it is important to understand the effects of enhanced grain size and dopant content on UO 2 fuel fracture. UO 2 pellets were doped with 0.1 wt% TiO 2 and 0.3 wt% Cr 2 O 3 to alter density and grain size. Inductively coupled plasma mass spectroscopy measured dopant levels pre- and post-sintering. X-ray diffraction revealed lattice changes and microstrain via Rietveld refinement. Field emission scanning electron microscopy determined grain sizes of approximately 30 µm for TiO 2 doping and 7 µm for Cr 2 O 3 doping. Transverse rupture strength tests were performed on over 30 samples per dataset to obtain characteristic strength and Weibull modulus. Results indicate no statistical difference in fracture strength between 0.1 wt% TiO 2 doped UO 2 and undoped UO 2 , while 0.3 wt% Cr 2 O 3 doped UO 2 exhibited a 20% decrease in fracture strength. Doped UO 2 samples also showed reduced Weibull modulus compared to undoped UO 2 , suggesting increased scatter in fracture strength. This study's findings suggest that titanium and chromium oxide doping in UO 2 , regardless of grain size, induce residual stresses, decreasing fracture strength and increasing variability in fracture behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Arm structure in normal spiral galaxies, 1: Multivariate data for 492 galaxies

Multivariate data have been collected as part of an effort to develop a new classification system for spiral galaxies, one which is not necessarily based on subjective morphological properties. A sample of 492 moderately bright northern Sa and Sc spirals was chosen for future statistical analysis. New observations were made at 20 and 21 cm; the latter data are described in detail here. Infrared Astronomy Satellite (IRAS) fluxes were obtained from archival data. Finally, new estimates of arm pattern radomness and of local environmental harshness were compiled for most sample objects.

Magri, Christopher↗

Evaluating Optical Classification for Fermi Blazar Candidates with a Statistical Method Using Broadband Spectral Indices

We aim to test whether a blazar candidate of uncertain type (BCU) in the third Fermi catalog of active galactic nuclei (3LAC) can be potentially classified as a BL Lac object or a flat-spectrum radio quasar (FSRQ) by performing a statistical analysis of its broadband spectral properties. We find that 34% of the radio-selected BCUs (583 BCUs) are BL Lac-like and 20% are FSRQ-like, at a 90% level of confidence. Similarly, 77.3% of the X-ray-selected BCUs (176 BCUs) are evaluated as BL Lac-like and 6.8% may be FSRQ-like sources. And 88.7% of the BL Lac-like BCUs that have synchrotron peak frequencies available are high synchrotron peaked BL Lacs in the X-ray-selected BCUs. The percentages are accordingly 62% and 7.3% in the sample of 124 optically selected BCUs. The high ratio of the number of BL Lac-like sources to the number of FSRQ-like BCUs in the X-ray-selected and optically selected BCU samples is due to selection effects. Examining the consistency between our evaluation and spectroscopic identification case by case with a sample of 78 radio-selected BCUs, it is found that the statistical analysis and its resulting classifications agree with the results of the optical follow-up spectroscopic observations. Our observation campaign for high-|ρ{sub s}| BCUs selected with our method, i.e., |ρ{sub s}|>0.8, is ongoing.

79 ASTRONOMY AND ASTROPHYSICS↗

Determining the tensile strength of fuel surrogate TRISO-coated particle buffer, IPyC, and buffer-IPyC interlayer regions

A novel micro tensile sample fabrication technique for determining the tensile strength of the buffer, IPyC, and buffer -IPyC interlayer regions of surrogate (ZrO2) TRISO fuel particle layers was refined and implemented. Copper micro tensile samples served as baseline materials to verify the methods used. Tensile tests performed in this study, while limited in number, were analyzed using standard and Weibull statistics. As expected, the buffer layer was weakest, with an average ultimate tensile strength of 138.70 MPa, and the IPyC layer samples, were strongest, with an average ultimate tensile strength of 189.74 MPa. In the buffer -IPyC interface samples, all breaks occurred in the buffer region, though the average ultimate tensile strength of the samples, 159.80 MPa, was between the pure buffer and IPyC samples. These results suggest the interlayer region has unique properties, perhaps associated with pyrocarbon infiltration into the buffer layer during particle coating. All interlayer samples fractured within the buffer side; however, the stress strain behavior of some of these samples resembled the behavior of the IPyC layer samples. Here, the buffer and IPyC layer strengths had a normal distribution under Weibull analysis, while the interlayer region had a Rayleigh distribution. Further testing is needed to clarify both the standard and Weibull statistical results

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Inspection error and its adverse effects - A model with implications for practitioners

Inspection error has clearly been shown to have adverse effects upon the results desired from a quality assurance sampling plan. These effects upon performance measures have been well documented from a statistical point of view. However, little work has been presented to convince the QC manager of the unfavorable cost consequences resulting from inspection error. This paper develops a very general, yet easily used, mathematical cost model. The basic format of the well-known Guthrie-Johns model is used. However, it is modified as required to assess the effects of attributes sampling errors of the first and second kind. The economic results, under different yet realistic conditions, will no doubt be of interest to QC practitioners who face similar problems daily. Sampling inspection plans are optimized to minimize economic losses due to inspection error. Unfortunately, any error at all results in some economic loss which cannot be compensated for by sampling plan design; however, improvements over plans which neglect the presence of inspection error are possible. Implications for human performance betterment programs are apparent, as are trade-offs between sampling plan modification and inspection and training improvements economics.

Collins, R. D., Jr.↗