Search NASA⌕ Search

SEARCH · Search NASA

Results for “randomized sampling”

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 415 records · Page 23

Performance of plane wavefront Fizeau interferometers in power spectral density measurements with tilted plane optics

Binary pseudo-random array (BPRA) test samples are useful devices for calibrating the instrument transfer function (ITF) of Fizeau interferometers, interferometric microscopes, and other optical and non-optical surface and wavefront metrology instruments. The intrinsic white noise character of the power spectral density (PSD) function of the BPRA pattern simplifies the extraction of the ITF from the measured PSD. The ITF determined in a dedicated calibration experiment can be used to reconstruct the surface height profile from the measured data, effectively enhancing the instrument’s spatial resolution. For a high confidence reconstruction procedure, a reliable analytical model of the IFT is desirable. Usually, the model accounts for the contributions to the ITF related the imperfections of the instrument’s optical and detector systems. Here, we experimentally demonstrate that in the case of surface height metrology with Fizeau interferometers, the PSD measurements and, therefore, the efficacy of the ITF calibration of the tool, are strongly affected by the instrument data acquisition and processing procedures, as well as by the shape of the optic under test and its alignment with respect to the interferometer.

Yashchuk, V V↗

Characteristic properties of macroporous silicon sintering in an argon atmosphere

The temperature and time dependences of the sintering of macroporous silicon in Ar or Ar + 3% H{sub 2} are studied. The contribution of various mechanisms governing this process is determined. The specific features of the sintering of macroporous silicon are examined by means of isochronous and isothermal annealing of the samples with ordered and random macropores in the temperature range 1000–1225°C. It is found that the sintering of macroporous silicon under atmospheric pressure in an inert gas flow containing 2 × 10–4% O{sub 2} is greatly affected by thermal etching. Thermal etching competes with the substance-transfer processes characteristic of sintering and hinders the formation of a defect-free surface crust. The reason for etching consists in that gaseous silicon monoxide is generated and then carried away by the gas flow. The etching effect is dominant in the low-temperature range and is independent of whether H{sub 2} is added. The values obtained for the activation energy of the silicon diffusion coefficient, E{sub a} = 2.57 eV, and for the exponent n = 3.31–3.74 in the time dependence of the pore radius, r ~ t{sup 1/n} are indicative of a mixed substance-transfer mechanism via the surface and volume diffusion of silicon atoms.

36 MATERIALS SCIENCE↗

Learning Planar Ising Models Software

Learning Planar Ising Models is a software package written in Matlab for learning relationships among variable in a dataset using graphical models. The software package implements a generally-applicable algorithm for learning planar Ising models from any multivariate dataset. The code provides an algorithm for learning the best planar Ising model to approximate an arbitrary collection of binary random variables (possibly from sample data). Given the set of all pairwise correlations among variables, we select a planar graph and optimal planar Ising model defined on this graph to best approximate that set of correlations. The software includes demonstrations of the algorithm in simulations and for applications on publicly available datasets. Details of the algorithm, demonstration simulations, and applications are given in Johnson, et al; 2016. Reference: Johnson, J. K., Oyen, D., Chertkov, M., and Netrapalli, P. (2016). Learning planar Ising models. Journal of Machine Learning Research.

Oyen, Diane↗

Learning planar Ising models

Inference and learning of graphical models are both well-studied problems in statistics and machine learning that have found many applications in science and engineering. However, exact inference is intractable in general graphical models, which suggests the problem of seeking the best approximation to a collection of random variables within some tractable family of graphical models. In this paper, we focus on the class of planar Ising models, for which exact inference is tractable using techniques of statistical physics. Based on these techniques and recent methods for planarity testing and planar embedding, we propose a greedy algorithm for learning the best planar Ising model to approximate an arbitrary collection of binary random variables (possibly from sample data). Given the set of all pairwise correlations among variables, we select a planar graph and optimal planar Ising model defined on this graph to best approximate that set of correlations. Finally, we demonstrate our method in simulations and for two applications: modeling senate voting records and identifying geo-chemical depth trends from Mars rover data.

97 MATHEMATICS AND COMPUTING↗

Identification of RR Lyrae stars in multiband, sparsely-sampled data from the Dark Energy Survey using template fitting and Random Forest classification

Many studies have shown that RR Lyrae variable stars (RRL) are powerful stellar tracers of Galactic halo structure and satellite galaxies. The Dark Energy Survey (DES), with its deep and wide coverage (g ~ 23.5 mag in a single exposure; over 5000 deg2) provides a rich opportunity to search for substructures out to the edge of the Milky Way halo. However, the sparse and unevenly sampled multiband light curves from the DES wide-field survey (a median of four observations in each of grizY over the first three years) pose a challenge for traditional techniques used to detect RRL. We present an empirically motivated and computationally efficient template-fitting method to identify these variable stars using three years of DES data. When tested on DES light curves of previously classified objects in SDSS stripe 82, our algorithm recovers 89% of RRL periods to within 1% of their true value with 85% purity and 76% completeness. Using this method, we identify 5783 RRL candidates, ~28% of which are previously undiscovered. This method will be useful for identifying RRL in other sparse multiband data sets.

79 ASTRONOMY AND ASTROPHYSICS↗

Extraction efficiency and implications for absolute quantitation of propranolol in mouse brain, liver and kidney thin tissue sections using droplet-based liquid microjunction surface sampling-HPLC ESI-MS/MS

Currently, absolute quantitation aspects of droplet-based surface sampling for thin tissue analysis using a fully automated autosampler/HPLC-ESI-MS/MS system are not fully evaluated. Knowledge of extraction efficiency and its reproducibility is required to judge the potential of the method for absolute quantitation of analytes from thin tissue sections. Methods: Adjacent thin tissue sections of propranolol dosed mouse brain (10- μm-thick), kidney (10- μm-thick) and liver (8-, 10-, 16- and 24- μm-thick) were obtained. Absolute concentration of propranolol was determined in tissue punches from serial sections using standard bulk tissue extraction protocols and subsequent HPLC separations and tandem mass spectrometric analysis. These values were used to determine propranolol extraction efficiency from the tissues with the droplet-based surface sampling approach. Results: Extraction efficiency of propranolol using 10- μm-thick brain, kidney and liver thin tissues using droplet-based surface sampling varied between ~45-63%. Extraction efficiency decreased from ~65% to ~36% with liver thickness increasing from 8 μm to 24 μm. Randomly selecting half of the samples as standards, precision and accuracy of propranolol concentrations obtained for the other half of samples as quality control metrics were determined. Resulting precision ( ±15%) and accuracy ( ±3%) values, respectively, were within acceptable limits. In conclusion, comparative quantitation of adjacent mouse thin tissue sections of different organs and of various thicknesses by droplet-based surface sampling and by bulk extraction of tissue punches showed that extraction efficiency was incomplete using the former method, and that it depended on the organ and tissue thickness. However, once extraction efficiency was determined and applied, the droplet-based approach provided the required quantitation accuracy and precision for assay validations. Furthermore, this means that once the extraction efficiency was calibrated for a given tissue type and drug, the droplet-based approach provides a non-labor intensive and high-throughput means to acquire spatially resolved quantitative analysis of multiple samples of the same type.

59 BASIC BIOLOGICAL SCIENCES↗

Neutron Next-Event Estimators Kinematics (Rev.2)

This paper reviews the kinematics of neutron elastic and inelastic scattering with moving and stationary targets for contributions to Neutron Next-Event Estimators (NEEs). NEEs are often used for simulating the response of detectors in locations that have few random walk particles. Contributions to NEEs from collisions differ from sampling the outgoing particle state in standard random walk collisions as the location of the estimator is fixed, thus fixing the scattering angle between the incoming and outgoing directions. To calculate contributions to NEEs, the outgoing energies and probabilities of scatter toward the estimator must be calculated. This paper presents the most general case that encompasses elastic and inelastic scattering with both moving and stationary targets. The moving target case is required for elastic scattering with thermal motion due to the free gas thermal treatment approximation. The stationary target equations are presented for comparison with equations presented in canonical Monte Carlo texts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The 3D Lyman- α forest power spectrum from eBOSS DR16

We measure the three-dimensional power spectrum (P3D) of the transmitted flux in the Lyman-α (Ly α) forest using the complete extended Baryon Oscillation Spectroscopic Survey data release 16 (eBOSS DR16). This sample consists of ~205 000 quasar spectra in the redshift range 2 ≤ z ≤ 4 at an effective redshift z = 2.334. We propose a pair-count spectral estimator in configuration space, weighting each pair by exp( i k ∙ r), for wave vector k and pixel pair separation r, effectively measuring the anisotropic power spectrum without the need for fast Fourier transforms. This accounts for the window matrix in a tractable way, avoiding artefacts found in Fourier-transform based power spectrum estimators due to the sparse sampling transverse to the line of sight of Ly α skewers. We extensively test our pipeline on two sets of mocks: (i) idealized Gaussian random fields with a sparse sampling of Ly α skewers, and (ii) log-normal LyaCoLoRe mocks including realistic noise levels, the eBOSS survey geometry and contaminants. On eBOSS DR16 data, the Kaiser formula with a non-linear correction term obtained from hydrodynamic simulations yields a good fit to the power spectrum data in the range $(0.02 ≤ k ≤ 0.35)$ h Mpc -1 at the 1–2σ level with a covariance matrix derived from LyaCoLoRe mocks. We demonstrate a promising new approach for full-shape cosmological analyses of Ly α forest data from cosmological surveys such as eBOSS, the currently observing Dark Energy Spectroscopic Instrument and future surveys such as the Prime Focus Spectrograph, WEAVE-QSO, and 4MOST.

79 ASTRONOMY AND ASTROPHYSICS↗

Interferometric measurements of refractive index and dispersion at high pressure

Abstract We describe a high precision interferometer system to measure the pressure dependence of the refractive index and its dispersion in the diamond anvil cell (DAC). The reflective Fabry–Perot fringe patterns created by both a white light and a monochromatic beam are recorded to determine both the sample thickness and its index at the laser wavelength and to characterize the dispersion in the visible range. Advances in sample preparation, optical setup, and data analysis enable us to achieve $$10^{-4}$$ 10 - 4 random uncertainty, demonstrated with an air sample, a factor of five improvement over the best previous DAC measurement. New data on $${\text{H}}_{2}\text{O}$$ H 2 O liquid water and ice VI up to 2.21 GPa at room temperature illustrate how higher precision measurements of the index and its optical dispersion open up new opportunities to reveal subtle changes in the electronic structure of water at high pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An adaptive local reduced basis method for solving PDEs with uncertain inputs and evaluating risk

Many physical systems are modeled using partial differential equations (PDEs) with uncertain or random inputs. For such systems, naively propagating a fixed number of samples of the input probability law (or an approximation thereof) through the PDE is often inadequate to accurately quantify the "risk" associated with critical system responses. In this paper, we develop a goal-oriented, adaptive sampling and local reduced basis approximation for PDEs with random inputs. Our method determines a set of samples and an associated (implicit) Voronoi partition of the parameter domain on which we build local reduced basis approximations of the PDE solution. The samples are selected in an adaptive manner using an a posteriori error indicator. A notable advantage of the proposed approach is that the computational cost of the approximation during the adaptive process remains constant. We provide theoretical error bounds for our approximation and numerically demonstrate the performance of our method when compared to widely used adaptive sparse grid techniques. In addition, we tailor our approach to accurately quantify the risk of quantities of interest that depend on the PDE solution. We demonstrate our method on an advection-diffusion example and a Helmholtz example. Keywords: stochastic partial differential equation, adaptive sampling, reduced basis, a posteriori error estimation, coherent risk measure, conditional value-at-risk.

42 ENGINEERING↗

Reproducible emission from nonlinear random lasers

Multiple scattering of light serves as a mechanism for feedback in random lasers. Consequently, internal spatial mode patterns, lasing wavelengths, and output directionality can all be random. Strong mode interaction can occur in such devices due to spatially overlapping modes resulting in nonlinearity with respect to the pump input power. Nevertheless, temporal coherence and lasing mode amplitude can be fixed at a constant pumping rate. This is a property desirable for applications where unique randomness is exploited but expected to be reliable over time, such as physical unclonable functions. Random lasers can also be cheaply and easily fabricated, exhibit relatively low lasing thresholds and high emission intensity. However, the precise scattering properties of such structures and fluctuations in the pump field can make device emission irreproducible, thereby limiting random laser applications. Here, in this work, we directly compare the random lasing spectra from zinc oxide samples fabricated in four distinct ways: spin-coating, sputtering, solgel deposition, and atomic layer deposition. The particular method of fabrication has a strong impact. Samples made through atomic layer deposition here exhibit both reproducibility and strong nonlinearity desirable for applications. Randomness in emission spectra persists across hundreds of repeated and averaged measurements irrespective of spatial location and is demonstrably nonlinear with respect to input signal intensity.

47 OTHER INSTRUMENTATION↗

Deep Learning without Global Optimization by Random Fourier Neural Networks

Here we introduce a new training algorithm for deep neural networks that utilize random complex exponential activation functions. Our approach employs a Markov chain Monte Carlo sampling procedure to iteratively train network layers, avoiding global and gradient-based optimization while maintaining error control. It consistently attains the theoretical approximation rate for residual networks with complex exponential activation functions, determined by network complexity. Additionally, it enables efficient learning of multiscale and high-frequency features, producing interpretable parameter distributions. Despite using sinusoidal basis functions, we do not observe Gibbs phenomena in approximating discontinuous target functions.

97 MATHEMATICS AND COMPUTING↗

Method for conducting in situ high-temperature digital image correlation with simultaneous synchrotron measurements under thermomechanical conditions

This work presents a novel method of obtaining in situ strain measurements at high temperature by simultaneous digital image correlation (DIC), which provides the total strain on the specimen surface, and synchrotron x-ray diffraction (XRD), which provides lattice strains of crystalline materials. DIC at high temperature requires specialized techniques to overcome the effects of increased blackbody radiation that would otherwise overexpose the images. The technique presented herein is unique in that it can be used with a sample enclosed in an infrared heater, remotely and simultaneously with synchrotron XRD measurements. The heater included a window for camera access, and the light of the heater lamps is used as illumination. High-temperature paint is used to apply a random speckle pattern to the sample to allow the tracking of displacements and the calculation of the DIC strains. An inexpensive blue theatrical gel filter is used to block interfering visible and infrared light at high temperatures. This technique successfully produces properly exposed images at 870 °C and is expected to perform similarly at higher temperatures. The average strains measured by DIC were validated by an analytical calculation of the theoretical strain. Simultaneous DIC and XRD strain measurements of Inconel 718 (IN718) tensile test specimens were performed under thermal and mechanical loads and evaluated. Finally, this approach uses the fact that with DIC, the total strain is measured, including plastic strain, while with XRD, only elastic strain is captured. The observed differences were discussed with respect to the effective deformation mechanisms.

47 OTHER INSTRUMENTATION↗

Do Epistemic Uncertainties Allow for Replacing Microstructural Experiments with Reconstruction Algorithms?

A novel problem in computational materials modeling is addressed: “Are the computational microstructure reconstruction techniques reliable enough to replace experiments?” Here, “reliable” computations are associated with producing “expected” reconstructions that are adequately close to the experimental data. The output of computational techniques can deviate from the experimental measurements because of the epistemic uncertainties in the algorithms. In this work, an analytical formulation for quantification of epistemic uncertainties in a microstructure reconstruction algorithm based on Markov random field is presented. The method is used to predict the large-scale spatial distribution of a microstructure given an experimental input measured over a small spatial domain. However, small variations are observed on the microstructural features of the synthesized samples due to the Markov random field algorithm. The proposed analytical technique aims to quantify these uncertainties and estimate their propagation to the macroscale material properties to provide a significant understanding on how reliable it is to replace the experiments with the Markov random field model to predict microstructural maps over large spatial domains.

97 MATHEMATICS AND COMPUTING↗

Southwest Pacific tropical cyclone development classification utilizing machine learning and synoptic composites

This study evaluates the ability of machine learning algorithms to classify tropical depressions (TDs) and tropical storms (TSs) in the western region of the southwest Pacific Ocean (SWPO). Decision rules are generated to predict the environment required for a depression to fully develop into a mature storm, and the most influential predictors in the classification decision are ranked. TD and TS are discriminated based on a maximum sustained wind speed threshold (≥17 ms -1 ). Various aerosol, thermodynamic, and dynamic parameters are extracted closest to the initiation point of each non-developing and developing sample. The covariates associated with each labelled sample are used to train a decision tree and random forest model. Results using a testing dataset suggest the random forest approach more accurately distinguishes between non-developing and developing samples. The classification accuracy of the decision tree and random forest are 72% and 91%, respectively. Random forest outperformed the decision tree by providing higher accuracy in test data. The most important variables for binary classification are sea salt aerosol optical depth (AOD), 1,000 mb relative humidity, and sea surface temperature. AOD is a quantitative estimate of the aerosols presents in the air through the extinction of a ray of light as it passes through the atmosphere. Mean composite maps constructed in an unsupervised manner have been created for the most important variables identified by the random forest classifier during TD and TS events to highlight the difference in geophysical and aerosol variables' climatology during the two different classifications. This work will advance the risk management strategies for northeastern Australia and other SWPO basin islands to control their tropical cyclone related losses through prioritizing forecasting variables that are the strongest predictors of the strengthening of tropical depressions into tropical cyclones.

54 ENVIRONMENTAL SCIENCES↗

Characterizing quantum supremacy in near-term devices

A critical question for quantum computing in the near future is whether quantum devices without error correction can perform a well-defined computational task beyond the capabilities of supercomputers. Such a demonstration of what is referred to as quantum supremacy requires a reliable evaluation of the resources required to solve tasks with classical approaches. Here, we propose the task of sampling from the output distribution of random quantum circuits as a demonstration of quantum supremacy. We extend previous results in computational complexity to argue that this sampling task must take exponential time in a classical computer. We introduce cross-entropy benchmarking to obtain the experimental fidelity of complex multiqubit dynamics. This can be estimated and extrapolated to give a success metric for a quantum supremacy demonstration. We study the computational cost of relevant classical algorithms and conclude that quantum supremacy can be achieved with circuits in a two-dimensional lattice of 7 × 7 qubits and around 40 clock cycles. This requires an error rate of around 0.5% for two-qubit gates (0.05% for one-qubit gates), and it would demonstrate the basic building blocks for a fault-tolerant quantum computer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Synchrotron X-ray microbeam diffraction measurements of full elastic long range internal strain and stress tensors in commercial-purity aluminum processed by multiple passes of equal-channel angular pressing

Synchrotron X-ray microbeam diffraction was used to measure the full elastic long range internal strain and stress tensors of low dislocation density regions within the submicrometer grain/subgrain structure of equal-channel angular pressed (ECAP) aluminum alloy AA1050 after 1, 2, and 8 passes using route B C . This is the first time that full tensors were measured in plastically deformed metals at this length scale. The maximum (most tensile or least compressive) principal elastic strain directions for the unloaded 1 pass sample for the grain/subgrain interiors align well with the pressing direction, and are more random for the 2 and 8 pass samples. The measurements reported here indicate that the local stresses and strains become increasingly isotropic (homogenized) with increasing ECAP passes using route BC. The average maximum (in magnitude) LRISs are -0.43 σ a for 1 pass, -0.44 σ a for 2 pass, and 0.14 σ a for the 8 pass sample. Furthermore, these LRISs are larger than those reported previously because those earlier measurements were unable to measure the full stress tensor. Significantly, the measured stresses are inconsistent with the two-component composite model.

36 MATERIALS SCIENCE↗

Random Forest Regressor-Based Approach for Detecting Fault Location and Duration in Power Systems

Power system failures or outages due to short-circuits or “faults” can result in long service interruptions leading to significant socio-economic consequences. It is critical for electrical utilities to quickly ascertain fault characteristics, including location, type, and duration, to reduce the service time of an outage. Existing fault detection mechanisms (relays and digital fault recorders) are slow to communicate the fault characteristics upstream to the substations and control centers for action to be taken quickly. Fortunately, due to availability of high-resolution phasor measurement units (PMUs), more event-driven solutions can be captured in real time. In this paper, we propose a data-driven approach for determining fault characteristics using samples of fault trajectories. A random forest regressor (RFR)-based model is used to detect real-time fault location and its duration simultaneously. This model is based on combining multiple uncorrelated trees with state-of-the-art boosting and aggregating techniques in order to obtain robust generalizations and greater accuracy without overfitting or underfitting. Four cases were studied to evaluate the performance of RFR: 1. Detecting fault location (case 1), 2. Predicting fault duration (case 2), 3. Handling missing data (case 3), and 4. Identifying fault location and length in a real-time streaming environment (case 4). A comparative analysis was conducted between the RFR algorithm and state-of-the-art models, including deep neural network, Hoeffding tree, neural network, support vector machine, decision tree, naive Bayesian, and K-nearest neighborhood. Experiments revealed that RFR consistently outperformed the other models in detection accuracy, prediction error, and processing time.

24 POWER TRANSMISSION AND DISTRIBUTION↗