Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Evaluation of Treadmill Exercise in a Lower Body Negative Pressure Chamber as a Countermeasure for Weightlessness-Induced Bone Loss: a Bed Rest Study with Identical Twins

Counteracting bone loss is required for future space exploration. We evaluated the ability of treadmill exercise in a LBNP chamber to counteract bone loss in a 30-day bed rest study. Eight pairs of identical twins were randomly assigned to sedentary control or exercise groups. Exercise within LBNP decreased the bone resorption caused by bed rest and may provide a countermeasure for spaceflight. INTRODUCTION: Bone loss is one of the greatest physiological challenges for extended-duration space missions. The ability of exercise to counteract weightlessness-induced bone loss has been studied extensively, but to date, it has proven ineffective. We evaluated the effectiveness of a combination of two countermeasures-treadmill exercise while inside a lower body negative pressure (LBNP) chamber-on bone loss during a 30-day bed rest study. MATERIALS AND METHODS: Eight pairs of identical twins were randomized into sedentary (SED) or exercise/LBNP (EX/LBNP) groups. Blood and urine samples were collected before, several times during, and after the 30-day bed rest period. These samples were analyzed for markers of bone and calcium metabolism. Repeated measures ANOVA was used to determine statistical significance. Because identical twins were used, both time and group were treated as repeated variables. RESULTS: Markers of bone resorption were increased during bed rest in samples from sedentary subjects, including the collagen cross-links and serum and urinary calcium concentrations. For N-telopeptide and deoxypyridinoline, there were significant (p < 0.05) interactions between group (SED versus EX/LBNP) and phase of the study (sample collection point). Pyridinium cross-links were increased above pre-bed rest levels in both groups, but the EX/LBNP group had a smaller increase than the SED group. Markers of bone formation were unchanged by bed rest in both groups. CONCLUSIONS: These data show that this weight-bearing exercise combined with LBNP ameliorates some of the negative effects of simulated weightlessness on bone metabolism. This protocol may pave the way to counteracting bone loss during spaceflight and may provide valuable information about normal and abnormal bone physiology here on Earth.

Twin Study↗
DOE OSTI2015

Lens statistics with gravitationally lensed yet morphologically regular images

GRAMORs are GRAvitationally highly magnified yet MORphologically regular images. An example of this phenomenon was discovered in the cluster MACS J1149.5+2223 in 2009. We investigate the lens statistics of GRAMORs in detail. Assuming a NFW profile for a sample of clusters, we calculate the expected number and redshift distribution of GRAMORs using parameters from COSMOS data for the number density of the background galaxy. A model with a cluster placed at z = 0.544 based on WMAP5 cosmology predicts the redshift of a GRAMOR at z≃1.49 which is close to the observed z = 1.4906. These results show that the expected number of GRAMORs is about two per cluster in the most likely case, and thus a large number of GRAMORs would be observed in a systematic survey. The probability distribution function of source redshift for GRAMORs depends strongly on dark energy and may be useful for constraining the nature of dark energy.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1974

Thermal histories and crystal distributions in partly devitrified lunar glasses cooled by radiation

Calculations of the expected crystallization behavior of selected lunar compositions (60095 and 14259) are presented. Knowledge of this behavior combined with measurements of the state of crystallinity in partially devitrified lunar samples permits the determination of the thermal history of the samples. The state of crystallinity of a specimen is described by a statistical crystal distribution function psi. The heat flow problems of a plate of opaque glass quenched on an aluminum block and of a semi-transparent glass sphere cooling by radiation are analyzed. The results are combined with laboratory measurements of the crystal growth velocity and the viscosity, and with theoretical expressions for the nucleation frequency (homogeneous and heterogeneous) to obtain psi. It is found that the calculated psi is much too large in comparison with laboratory experiments on lunar samples. Possible reasons for the discrepancy are discussed, and the most likely cause is that the assumed nucleation barrier is too small. It also appears likely that these materials have relatively few nucleating heterogeneities.

Hopper, R. W.↗
DOE OSTI2026

Drop clustering and drop size correlations from holographic imagery suggest cloud droplet spectral broadening via entrainment-mixing

The question of how droplets rapidly grow large enough to initiate collision-coalescence has persisted for decades. Many theories explaining the production of sufficiently large drops (i.e., those in the “bottleneck” size range; ∼ 25–50 µm diameters) involve drop clustering on millimeter scales. A novel method is introduced to evaluate drop clustering trends particle-by-particle (i.e., the number/proximity of neighboring drops for given droplets; defined as drop clustering fields) which are diagnosed relative to drops within their shared drop environments – in contrast to previous studies which diagnose drop clustering of defined sample volumes, or in terms of absolute length scales. Specifically, this study evaluates the statistical likelihood that drops of a given size are associated with either a significant number of neighboring drops, or are significantly isolated from neighboring drops. Observations are acquired from the HOLODEC during the Cloud System Evolution in the Trades campaign, which sampled subtropical marine clouds. The HOLODEC measures drop size distributions and the 3D spatial coordinates of droplets. Results show drops within the bottleneck size range (diameters of ∼ 25–50 µm) are most likely to be significantly isolated from neighboring drops. This “isolated large drop trend” is primarily observed at subsaturated conditions, suggesting entrainment is the contributing factor. Holograms associated with this trend are more likely to have broader drop size distributions, larger maximum drop sizes and overly regions where precipitation reaches the lowest altitudes from the sampled cloud, suggesting entrainment-mixing drop size distribution broadening is a relevant precipitation-initiation mechanism.

D'Alessandro, John J. [Univ. of Washington, Seattl↗
NASA NTRS2011

An Empirical State Error Covariance Matrix for Batch State Estimation

State estimation techniques serve effectively to provide mean state estimates. However, the state error covariance matrices provided as part of these techniques suffer from some degree of lack of confidence in their ability to adequately describe the uncertainty in the estimated states. A specific problem with the traditional form of state error covariance matrices is that they represent only a mapping of the assumed observation error characteristics into the state space. Any errors that arise from other sources (environment modeling, precision, etc.) are not directly represented in a traditional, theoretical state error covariance matrix. Consider that an actual observation contains only measurement error and that an estimated observation contains all other errors, known and unknown. It then follows that a measurement residual (the difference between expected and observed measurements) contains all errors for that measurement. Therefore, a direct and appropriate inclusion of the actual measurement residuals in the state error covariance matrix will result in an empirical state error covariance matrix. This empirical state error covariance matrix will fully account for the error in the state estimate. By way of a literal reinterpretation of the equations involved in the weighted least squares estimation algorithm, it is possible to arrive at an appropriate, and formally correct, empirical state error covariance matrix. The first specific step of the method is to use the average form of the weighted measurement residual variance performance index rather than its usual total weighted residual form. Next it is helpful to interpret the solution to the normal equations as the average of a collection of sample vectors drawn from a hypothetical parent population. From here, using a standard statistical analysis approach, it directly follows as to how to determine the standard empirical state error covariance matrix. This matrix will contain the total uncertainty in the state estimate, regardless as to the source of the uncertainty. Also, in its most straight forward form, the technique only requires supplemental calculations to be added to existing batch algorithms. The generation of this direct, empirical form of the state error covariance matrix is independent of the dimensionality of the observations. Mixed degrees of freedom for an observation set are allowed. As is the case with any simple, empirical sample variance problems, the presented approach offers an opportunity (at least in the case of weighted least squares) to investigate confidence interval estimates for the error covariance matrix elements. The diagonal or variance terms of the error covariance matrix have a particularly simple form to associate with either a multiple degree of freedom chi-square distribution (more approximate) or with a gamma distribution (less approximate). The off diagonal or covariance terms of the matrix are less clear in their statistical behavior. However, the off diagonal covariance matrix elements still lend themselves to standard confidence interval error analysis. The distributional forms associated with the off diagonal terms are more varied and, perhaps, more approximate than those associated with the diagonal terms. Using a simple weighted least squares sample problem, results obtained through use of the proposed technique are presented. The example consists of a simple, two observer, triangulation problem with range only measurements. Variations of this problem reflect an ideal case (perfect knowledge of the range errors) and a mismodeled case (incorrect knowledge of the range errors).

Frisbee, Joseph H., Jr.↗
DOE OSTI2026

Uncovering hidden bias in neutron diffraction residual strain measurements

When calculating residual strain via neutron or X-ray diffraction, uncertainties propagated from the peak fit are often inadequate to describe the true scatter of measurements about a singular strain state, such as one that should describe a macroscopic continuum. Because diffraction is inherently a selective process, orientation-dependent scatter arises from the sub-sampling of strong microstructure and strain gradients. This paper investigates the appropriateness of propagated uncertainties with reference to their original intention, i.e. noise about a mean value. Thirty-six unique orientations of strain measurements are taken at multiple locations within an additive friction-stir deposition component with fine-scale gradients (∼200 µm) of plastic strain, texture and residual elastic strain. Multiple strain and stress calculation pathways are compared: direct substitution of three measurements into Hooke's law, direct inversion of any six unique orientations into the strain state tensor and thirty-six measurement least-squares estimation. For the last two cases, the appropriateness of the uncertainty interval is statistically evaluated on the basis of a physical constraint: common agreement under the strain transformation law. For this sample, the direct inversion of six measurements retains a conservative estimate of the uncertainty. However, propagated uncertainties in the least-squares solution greatly underestimate the true experimental scatter. A simple pathway to estimate appropriate uncertainty intervals is suggested. These results demonstrate that the interpretation of uncertainty in residual strain is strongly dependent on intrinsic sample-dependent effects, and that oversampling orientations and statistical analysis can give more accurate results with realistic uncertainties.

36 MATERIALS SCIENCE↗
DOE OSTI2015

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↗
DOE OSTI2020

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↗
DOE OSTI2023

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↗
DOE OSTI2017

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↗
NASA NTRS1999

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↗
DOE OSTI2015

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↗
DOE OSTI2015

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↗
DOE OSTI2014

The statistical nature of the brightest group galaxies

We examine the statistical properties of the brightest group galaxies (BGGs) by using a complete spectroscopic sample of groups/clusters of galaxies selected from the Data Release 7 of the Sloan Digital Sky Survey. We test whether BGGs and other bright members of groups are consistent with an ordered population among the total population of group galaxies. We find that the luminosity distributions of BGGs do not follow the predictions from the order statistics (OS). The average luminosities of BGGs are systematically brighter than OS predictions. On the other hand, by properly taking into account the brightening effect of the BGGs, the luminosity distributions of the second brightest galaxies are in excellent agreement with the expectations of OS. The brightening of BGGs relative to the OS expectation is consistent with a scenario that the BGGs on average have overgrown about 20% masses relative to the other member galaxies. The growth (ΔM) is not stochastic but correlated with the magnitude gap (G {sub 1,} {sub 2}) between the brightest and the second brightest galaxy. The growth (ΔM) is larger for the groups having more prominent BGGs (larger G {sub 1,} {sub 2}) and averagely contributes about 30% of the final G {sub 1,} {sub 2} of the groups of galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

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↗
DOE OSTI2018

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↗
DOE OSTI2016

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↗