Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical sampling techniques”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Detection Limits of Low-mass, Long-period Exoplanets Using Gaussian Processes Applied to HARPS-N Solar Radial Velocities

Radial velocity (RV) searches for Earth-mass exoplanets in the habitable zone around Sun-like stars are limited by the effects of stellar variability on the host star. In particular, suppression of convective blueshift and brightness inhomogeneities due to photospheric faculae/plage and starspots are the dominant contribution to the variability of such stellar RVs. Gaussian process (GP) regression is a powerful tool for statistically modeling these quasi-periodic variations. We investigate the limits of this technique using 800 days of RVs from the solar telescope on the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N) spectrograph. These data provide a well-sampled time series of stellar RV variations. Into this data set, we inject Keplerian signals with periods between 100 and 500 days and amplitudes between 0.6 and 2.4 m s{sup −1}. We use GP regression to fit the resulting RVs and determine the statistical significance of recovered periods and amplitudes. We then generate synthetic RVs with the same covariance properties as the solar data to determine a lower bound on the observational baseline necessary to detect low-mass planets in Venus-like orbits around a Sun-like star. Our simulations show that discovering planets with a larger mass (∼0.5 m s{sup −1}) using current-generation spectrographs and GP regression will require more than 12 yr of densely sampled RV observations. Furthermore, even with a perfect model of stellar variability, discovering a true exo-Venus (∼0.1 m s{sup −1}) with current instruments would take over 15 yr. Therefore, next-generation spectrographs and better models of stellar variability are required for detection of such planets.

47 OTHER INSTRUMENTATION↗

Is Terzan 5 the remnant of a building block of the Galactic bulge? Evidence from APOGEE

ABSTRACT It has been proposed that the globular cluster-like system Terzan 5 is the surviving remnant of a primordial building block of the Milky Way bulge, mainly due to the age/metallicity spread and the distribution of its stars in the α–Fe plane. We employ Sloan Digital Sky Survey data from the Apache Point Observatory Galactic Evolution Experiment to test this hypothesis. Adopting a random sampling technique, we contrast the abundances of 10 elements in Terzan 5 stars with those of their bulge field counterparts with comparable atmospheric parameters, finding that they differ at statistically significant levels. Abundances between the two groups differ by more than 1σ in Ca, Mn, C, O, and Al, and more than 2σ in Si and Mg. Terzan 5 stars have lower [α/Fe] and higher [Mn/Fe] than their bulge counterparts. Given those differences, we conclude that Terzan 5 is not the remnant of a major building block of the bulge. We also estimate the stellar mass of the Terzan 5 progenitor based on predictions by the Evolution and Assembly of GaLaxies and their Environments suite of cosmological numerical simulations, concluding that it may have been as low as ∼3 × 108 M⊙ so that it was likely unable to significantly influence the mean chemistry of the bulge/inner disc, which is significantly more massive (∼1010 M⊙). We briefly discuss existing scenarios for the nature of Terzan 5 and propose an observational test that may help elucidate its origin.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistical Generation of Ocean Forcing With Spatiotemporal Variability for Ice Sheet Models

Melting of ice at the base of floating ice shelves that fringe the Antarctic ice sheet has been identified as a significant source of uncertainty in sea level rise projections. Part of this uncertainty derives from chaotic internal variability of the coupled ocean-atmosphere system. For numerical ice sheet model projections, this uncertainty has not previously been quantified because of the prohibitive computational expense of running large climate model ensembles. Here, we develop and demonstrate a technique that generates independent realizations of internal climate variability from a single climate model simulation. Building on previous developments in model emulation, this technique uses empirical orthogonal function decomposition and Fourier-phase randomization to generate statistically consistent realizations of spatiotemporal variability fields for the target climate variable. The method facilitates efficient sampling of a wide range of climate trajectories, which can also be incorporated within ice sheet or other physical models to represent feedback processes.

54 ENVIRONMENTAL SCIENCES↗

Collaborative: in situ visual analytics technologies for extreme scale combustion simulations

This project aims to drastically enhance the usability of in situ analysis and visualization for extreme-scale scientific simulations. Current exascale computing capabilities promise to offer greater predictive ability of simulations and to further push the frontiers of science and technology. However, to validate the simulation output at extreme scale, examine the modeled phenomena, and discover previously unknowns from the output data, the output must be reduced or transformed in situ as it is being generated during the simulation such that the amount of data to examine and store is kept to a minimum. Such in situ approaches allow us to process and analyze the data and any embedded geometry to an extent that would be prohibitively expensive, if not impossible, to perform as a post hoc task. While in situ processing has been demonstrated to be a feasible and promising approach, its full potential has not yet been leveraged. In this project, we have developed comprehensive enhancements to in situ technology based on probability distributions in data. Our research focuses on jointly developing new ways of interacting with massive statistical samples while creatively utilizing new state-of-the-art computational resources to push the boundaries of in situ exploration. Moreover, we have developed new time-dependent techniques to enable previously unattainable capabilities in areas such as intelligent simulation steering and precise feature identification. We have experimentally studied our design and implementation at NERSC and OLCF, and are able to leverage existing in situ infrastructures whenever possible. While the exemplar in this project is combustion, many other fields for which turbulent transport is important, e.g., fusion, climate, astrophysics among others, encounter similar issues as simulations scale up to the exascale. This project shows its potential to generate high impact on DOE missions since the resulting technology promises to improve scientists’ ability to rapidly and correctly interpret and tune extreme-scale simulations, leading to new scientific understanding and advancements.

97 MATHEMATICS AND COMPUTING↗

The HEAO A-1 all-sky survey of fast X-ray transients

Ten fast X-ray transient sources were discovered in the first complete pass of the sky by the HEAO A-1 Sky Survey Experiment, including 7 relatively faint sources and 3 bright ones. The survey technique and each of the transients are described. Four of the transients are identified with active coronal sources; the remaining ones are unidentified. The statistical properties of the total sample of 10 sources are discussed, including event rates, distribution on the sky, and the number-flux distribution.

Ambruster, Carol W.↗

Bayesian Approach to the Joint Inversion of Gravity and Magnetic Data, with Application to the Ismenius Area of Mars

This viewgraph presentation reviews a Bayesian approach to the inversion of gravity and magnetic data with specific application to the Ismenius Area of Mars. Many inverse problems encountered in geophysics and planetary science are well known to be non-unique (i.e. inversion of gravity the density structure of a body). In hopes of reducing the non-uniqueness of solutions, there has been interest in the joint analysis of data. An example is the joint inversion of gravity and magnetic data, with the assumption that the same physical anomalies generate both the observed magnetic and gravitational anomalies. In this talk, we formulate the joint analysis of different types of data in a Bayesian framework and apply the formalism to the inference of the density and remanent magnetization structure for a local region in the Ismenius area of Mars. The Bayesian approach allows prior information or constraints in the solutions to be incorporated in the inversion, with the "best" solutions those whose forward predictions most closely match the data while remaining consistent with assumed constraints. The application of this framework to the inversion of gravity and magnetic data on Mars reveals two typical challenges - the forward predictions of the data have a linear dependence on some of the quantities of interest, and non-linear dependence on others (termed the "linear" and "non-linear" variables, respectively). For observations with Gaussian noise, a Bayesian approach to inversion for "linear" variables reduces to a linear filtering problem, with an explicitly computable "error" matrix. However, for models whose forward predictions have non-linear dependencies, inference is no longer given by such a simple linear problem, and moreover, the uncertainty in the solution is no longer completely specified by a computable "error matrix". It is therefore important to develop methods for sampling from the full Bayesian posterior to provide a complete and statistically consistent picture of model uncertainty, and what has been learned from observations. We will discuss advanced numerical techniques, including Monte Carlo Markov

data analysis↗

Rapid measurement of soluble xylo-oligomers using near-infrared spectroscopy (NIRS) and multivariate statistics: calibration model development and practical approaches to model optimization

Rapid monitoring of biomass conversion processes using techniques such as near-infrared (NIR) spectroscopy can be substantially quicker and less labor-, resource-, and energy-intensive than conventional measurement techniques such as gas or liquid chromatography (GC or LC) due to the lack of solvents and preparation methods, as well as removing the need to transfer samples to an external lab for analytical evaluation. The purpose of this study was to determine the feasibility of rapid monitoring of a biomass conversion process using NIR spectroscopy combined with multivariate statistical modeling, and to examine the impact of (1) subsetting the samples in the original dataset by process location and (2) reducing the spectral range used in the calibration model on model performance. We develop multivariate calibration models for the concentrations of soluble xylo-oligosaccharides (XOS), monomeric xylose, and total solids at multiple points in a biomass conversion process which produces and then purifies XOS compounds from sugar cane bagasse. A single model using samples from multiple locations in the process stream showed acceptable performance as measured by standard statistical measures. However, compared to the single model, we show that separate models built by segregating the calibration samples according to process location show improved performance. We also show that combining an understanding of the sample spectra with simple multivariate analysis tools can result in a calibration model with a substantially smaller spectral range that provides essentially equal performance to the full-range model. We demonstrate that real-time monitoring of soluble xylo-oligosaccharides (XOS), monomeric xylose, and total solids concentration at multiple points in a process stream using NIR spectroscopy coupled with multivariate statistics is feasible. Segregation of sample populations by process location improves model performance. Models using a reduced spectral range containing the most relevant spectral signatures show very similar performance to the full-range model, reinforcing the importance of performing robust exploratory data analysis before beginning multivariate modeling.

09 BIOMASS FUELS↗

Long-term 4.6 µm Variability in Brown Dwarfs and a New Technique for Identifying Brown Dwarf Binary Candidates

Using a sample of 361 nearby brown dwarfs, we have searched for 4.6 μm variability, indicative of large-scale rotational modulations or large-scale, long-term changes on timescales of over 10 yr. Our findings show no statistically significant variability in Spitzer's Infrared Array Camera (IRAC) channel 2 (ch2) or Wide-field Infrared Survey Explorer W2 photometry. For Spitzer the ch2 1σ limits are ∼8 mmag for objects at 11.5 mag and ∼22 mmag for objects at 16 mag. This corresponds to no variability above 4.5% at 11.5 mag and 12.5% at 16 mag. We conclude that highly variable brown dwarfs, at least two previously published examples of which have been shown to have 4.6 μm variability above 80 mmag, are very rare. While analyzing the data, we also developed a new technique for identifying brown dwarf binary candidates in Spitzer data. We find that known binaries have IRAC ch2 point response function (PRF) flux measurements that are consistently dimmer than aperture flux measurements. We have identified 59 objects that exhibit such PRF versus aperture flux differences and are thus excellent binary brown dwarf candidates.

Hunter Brooks↗

An automated approach to the design of decision tree classifiers

The classification of large dimensional data sets arising from the merging of remote sensing data with more traditional forms of ancillary data is considered. Decision tree classification, a popular approach to the problem, is characterized by the property that samples are subjected to a sequence of decision rules before they are assigned to a unique class. An automated technique for effective decision tree design which relies only on apriori statistics is presented. This procedure utilizes a set of two dimensional canonical transforms and Bayes table look-up decision rules. An optimal design at each node is derived based on the associated decision table. A procedure for computing the global probability of correct classfication is also provided. An example is given in which class statistics obtained from an actual LANDSAT scene are used as input to the program. The resulting decision tree design has an associated probability of correct classification of .76 compared to the theoretically optimum .79 probability of correct classification associated with a full dimensional Bayes classifier. Recommendations for future research are included.

Argentiero, P.↗

Surface smoothing for laser powder-bed Ti-6Al-4V by a transient liquid phase

Surface roughness is the primary driver of fatigue for additively manufactured metals. To address surface roughness, this work introduces a new method to smooth features beyond line-of-sight without material removal. The method applies a coating that triggers local surface remelting by activating a eutectic reaction during heat treatment. The associated liquid phase then wets and isothermally solidifies into a smoother surface. For Ti-6Al-4V fabricated with laser powder bed fusion, samples with and without TLP smoothing (using a Cu coating) were characterized with a suite of techniques, including mechanical testing, electron backscatter diffraction, synchrotron X-ray tomography, and fractography. TLP smoothing reduced surface roughness by 80% and amplified compressive residual stress at the surface by about 50%. With statistically equivalent virtual microstructures, crystal plasticity scrutinized the roles of phases, porosity, and surface roughness. Although the tensile strain-to-failure was reduced to 1% strain, the TLP smoothing process increased high-cycle fatigue strength by about 20% compared to control samples, pointing to future opportunities to optimize the new process through various coating compositions and heat treatment schedules. Overall, this work establishes a new paradigm for treating surfaces of materials for smoothness and compressive residual stress.

Additive manufacturing↗

An Energy-Integrated Analysis for Measuring the Anomalous Precession Frequenzc for the Muon g - 2 Experiment at Fermilab

In the search for physics beyond the Standard Model, the Muon g - 2 Experiment at Fermilab (E989) will make the most precise measurement of the anomalous magnetic moment of the muon, a μ . Improvements in precision come from both increased statistics and new techniques to significantly reduce previous systematic uncertainties. The muon a μ is determined by extracting both the anomalous spin precession frequency, ω a , and the average magnetic field sampled by the muons, B. Traditionally an energy threshold analysis method which requires reconstruction of decay positrons from the muon decay, μ + → e + ν e ν μ has been used to extract ω a . This thesis will describe a novel, energy-integrated analysis method for ωa extraction that does not require full reconstruction of positron hits and lowers the energy threshold therefore sampling more of the stored muons. This new technique has different sensitivities to effects like gain fluctuations and positron pulse pile-up. The analysis results will focus on data taken in 2019, the second year of running the experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hydrogen partitioning in pure cast aluminum as determined by dynamic evolution rate measurements

Hydrogen in pure aluminum can be found in two different states. One is related to a presence in gas-filled pores, while the other state involves the formation of a solid solution between hydrogen and aluminum. The considered investigation is concerned with the distribution of the hydrogen between various states. A dynamic technique is employed to measure the evolution of hydrogen from commercially available samples of polycrystalline pure aluminum under ultrahigh vacuum conditions. The obtained data are compared with the results of a statistical analysis concerning the porosity in the cast aluminum. It was found that more than 99 pct of the hydrogen in the aluminum is located in large pores. Less than 1 pct of the hydrogen is partitioned between the solid solution and the small pores.

Outlaw, R. A.↗

A new method for simulating atmospheric turbulence for rotorcraft applications

Simulation of atmospheric turbulence as seen by a rotating blade element involves treatment of cyclostationary processes. Conventional filtering techniques do not lend themselves well to the generation of such turbulence sample functions as are required in rotorcraft flight dynamics simulation codes. A method to generate sample functions containing second order statistics of mean and covariance is presented. Compared to ensemble averaging involving excessive computer time, the novelty is to exploit cycloergodicity and thereby, replace ensemble averaging by averaging over a single path sample function of long duration. The method is validated by comparing its covariance results with the analytical and ensemble averaged results for a widely used 1-D turbulence approximation.

J. Riaz↗

BATSE analysis techniques for probing the GRB spatial and luminosity distributions

The Burst And Transient Source Experiment (BATSE) has measured homogeneity and isotropy parameters from an increasingly large sample of observed gamma-ray bursts (GRBs), while also maintaining a summary of the way in which the sky has been sampled. Measurement of both of these are necessary for any study of the BATSE data statistically, as they take into account the most serious observational selection effects known in the study of GRBs: beam-smearing and inhomogeneous, anisotropic sky sampling. Knowledge of these effects is important to analysis of GRB angular and intensity distributions. In addition to determining that the bursts are local, it is hoped that analysis of such distributions will allow boundaries to be placed on the true GRB spatial distribution and luminosity function. The technique for studying GRB spatial and luminosity distributions is direct. Results of BATSE analyses are compared to Monte Carlo models parameterized by a variety of spatial and luminosity characteristics.

Hakkila, Jon↗

Spectral analysis of groove spacing on Ganymede

The technique used to analyze groove spacing on Ganymede is presented. Data from Voyager images are used determine the surface topography and position of the grooves. Power spectal estimates are statistically analyzed and sample data is included.

Grimm, R. E.↗

A new method for simulating atmospheric turbulence for rotorcraft applications

Simulation of atmospheric turbulence as seen by a rotating blade element involves treatment of cyclostationary processes. Conventional filtering techniques do not lend themselves well to the generation of such turbulence sample functions as are required in rotorcraft flight dynamics simulation codes. A method to generate sample functions containing second-order statistics of mean and covariance is presented. Compared to ensemble averaging involving excessive computer time, the novelty is to exploit cycloergodicity and thereby, replace ensemble averaging by averaging over a single-path sample function of long duration. The method is validated by comparing its covariance results with the analytical and ensemble-averaged results for a widely used one-dimensional turbulence approximation.

Prasad, J. V. R.↗

Chemical studies of H chondrites. 4: New data and comparison of Antarctic suites

We report data for the trace elements Au, Co, Sb, Ga, Rb, Ag, Se, Cs, Te, Zn, Cd, Bi, Ti, and In (ordered by putative volatility during nebular condensation and accretion) determined by neutron activation analysis in 13 H5 chondrites from Victoria Land and 20 H4-6 chondrites from Queen Maud Land, Antarctica. These and earlier results provide Antarctic sample suites of 34 chondrites from Victoria Land and 25 from Queen Maud Land. Treatment of data for the most volatile 10 elements (Rb to In) in these studies by multivariate statistical techniques more robust, as well as more conservative, than conventional linear discriminant analysis and logistic regression demonstrates that compositions differ at marginally significant levels. This difference cannot be explained by trivial (terrestrial) causes and becomes more significant, despite the smaller size of the database, when comparisons are limited to data from a single analyst and when all upper limits are eliminated from consideration. The Victoria Land and Queen Maud Land suites have different mean terrestrial ages (approximately 300 kyr and approximately 100 kyr, respectively) and age distributions, suggesting that a time-dependent variation of chondritic sources with different thermal histories is responsible. As a result, these two Antarctic suites are, on average, chemically distinguishable from each other. Since H chondrites serve as a paradigm for other meteorite classes, these results indicate that the near-Earth populations of planetary materials varied with time on the 10(exp 5)-year timescale.

Wolf, Stephen F.↗

An Overview of Manufacturing Controls for Production of High-Consequence, Single-Use Systems

Extended Testing (Crowder, et al. (2025) ) is a reliability demonstration technique that can be used to dramatically reduce sample size requirements. Manufacturing controls are needed to supplement extended testing by identifying production issues that a reduced sample sizes might overlook, especially built-in, or latent, manufacturing defects. This report focuses on some of the most commonly used, yet most impactful, manufacturing control tools that are used to limit production-related defects and efficiently screen any remaining defects at final inspection. These tools include statistical process control (SPC), acceptance sampling, environmental stress screening (HASS and ESS), and mistake proofing. The goal is to minimize the probability that built-in defects ever reach the customer. In terms of nuclear weapons (NW), the goal is to prevent defective units from ever entering the nation’s NW stockpile. Examples of each of the techniques are illustrated with case studies.

42 ENGINEERING↗