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 1,171 records · Page 65

Abiotic proxies for predictive mapping of nearshore benthic assemblages: implications for marine spatial planning

Marine spatial planning (MSP) should assist managers in guiding human activities toward sustainable practices and in minimizing user conflicts in our oceans. A necessary first step is to quantify spatial patterns of marine assemblages in order to understand the ecosystem's structure, function, and services. However, the large spatial scale, high economic value, and density of human activities in nearshore habitats often makes quantifying this component of marine ecosystems especially daunting. To address this challenge, we developed an assessment method that employs abiotic proxies to rapidly characterize marine assemblages in nearshore benthic environments with relatively high resolution. We evaluated this assessment method along 300 km of the State of Maine's coastal shelf (<100 m depth), a zone where high densities of buoyed lobster traps typically preclude extensive surveys by towed sampling gear (i.e., otter trawls). During the summer months of 2010–2013, we implemented a stratified-random survey using a small remotely operated vehicle that allowed us to work around lobster buoys and to quantify all benthic megafauna to species. Stratifying by substrate, depth, and coastal water masses, we found that abiotic variables explained a significant portion of variance (37–59%) in benthic species composition, diversity, biomass, and economic value. Generally, the density, diversity, and biomass of assemblages significantly increased with the substrate complexity (i.e., from sand-mud to ledge). The diversity, biomass, and economic value of assemblages also decreased significantly with increasing depth. Last, demersal fish densities, sessile invertebrate densities, species diversity, and assemblage biomass increased from east to west, while the abundance of mobile invertebrates and economic value decreased, corresponding mainly to the contrasting water mass characteristics of the Maine Coastal Current system (i.e., summertime current direction, speed, and temperature). Integrating modeled predictions with existing GIS layers for abiotic conditions allowed us to scale up important assemblage attributes to define key foundational ecological principles of MSP and to find priority regions where some bottom-disturbing activities would have minimal impact to benthic assemblages. We conclude that abiotic proxies can be strong forcing functions for the assembly of marine communities and therefore useful tools for spatial extrapolations of marine assemblages in congested (heavily used) nearshore habitats.

54 ENVIRONMENTAL SCIENCES↗

Two-Step Estimation in a Heteroscedastic Linear Regression Model

We study the problem of estimating a parameter in some heteroscedastic linear regression model in the case where the regressors consist of all order statistics based on the sample of identically distributed not necessarily independent observations with finite second moment. It is assumed that the random errors depend on the parameter and distributions of the corresponding regressors. We propose a two-step procedure for finding explicit asymptotically normal estimators.

97 MATHEMATICS AND COMPUTING↗

Improved multifidelity Monte Carlo estimators based on normalizing flows and dimensionality reduction techniques

Here, we study the problem of multifidelity uncertainty propagation for computationally expensive models. In particular, we consider the general setting where the high-fidelity and low-fidelity models have a dissimilar parameterization both in terms of number of random inputs and their probability distributions, which can be either known in closed form or provided through samples. We derive novel multifidelity Monte Carlo estimators which rely on a shared subspace between the high-fidelity and low-fidelity models where the parameters follow the same probability distribution, i.e., a standard Gaussian. We build the shared space employing normalizing flows to map different probability distributions into a common one, together with linear and nonlinear dimensionality reduction techniques, active subspaces and autoencoders, respectively, which capture the subspaces where the models vary the most. We then compose the existing low-fidelity model with these transformations and construct modified models with an increased correlation with the high-fidelity model, which therefore yield multifidelity estimators with reduced variance. A series of numerical experiments illustrate the properties and advantages of our approaches.

97 MATHEMATICS AND COMPUTING↗

Next Steps for Compressively Sensed Video in Transmission Electron Microscopy

Here, compressive sensing (CS) is a fast-emerging approach for improving data throughput by taking advantage of the inherent compressibility of real world data—the same kind of compressibility that allows millions of users to stream full-screen video without breaking the internet. Transmission electron microscope (TEM) users regularly collect multi-GB data sets and then use aggressive data reduction strategies to draw out the tiny residual of Shannon entropy that answers the question they wanted to ask. Much of the raw-byte-count information in an acquisition is random, redundant, or irrelevant, yet its collection still represents substantial cost in terms of instrumentation bandwidth and sample damage.

47 OTHER INSTRUMENTATION↗

Influence of orbital symmetry on diffraction imaging with rescattering electron wave packets

The ability to directly follow and time-resolve the rearrangement of the nuclei within molecules is a frontier of science that requires atomic spatial and few-femtosecond temporal resolutions. While laser-induced electron diffraction can meet these requirements, it was recently concluded that molecules with particular orbital symmetries (such as pg) cannot be imaged using purely backscattering electron wave packets without molecular alignment. Here, we demonstrate, in direct contradiction to these findings, that the orientation and shape of molecular orbitals presents no impediment for retrieving molecular structure with adequate sampling of the momentum transfer space. We overcome previous issues by showcasing retrieval of the structure of randomly oriented O 2 and C 2 H 2 molecules, with π g and π u symmetries, respectively, and where their ionization probabilities do not maximize along their molecular axes. As a result, while this removes a serious bottleneck for laser-induced diffraction imaging, we find unexpectedly strong backscattering contributions from low-Z atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗

Smart darting diffusion Monte Carlo: Applications to lithium ion-Stockmayer clusters

In a recent investigation [K. Roberts et al., J. Chem. Phys. 136, 074104 (2012)], we have shown that, for a sufficiently complex potential, the Diffusion Monte Carlo (DMC) random walk can become quasiergodic, and we have introduced smart darting-like moves to improve the sampling. In this article, we systematically characterize the bias that smart darting moves introduce in the estimate of the ground state energy of a bosonic system. We then test a simple approach to eliminate completely such bias from the results. The approach is applied for the determination of the ground state of lithium ion-n–dipoles clusters in the n = 8–20 range. For these, the smart darting diffusion Monte Carlo simulations find the same ground state energy and mixed-distribution as the traditional approach for n < 14. In larger systems we find that while the ground state energies agree quantitatively with or without smart darting moves, the mixed-distributions can be significantly different. Some evidence is offered to conclude that introducing smart darting-like moves in traditional DMC simulations may produce a more reliable ground state mixed-distribution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental quantum Bernoulli factories via Bell-basis measurements

Randomness processing in the Bernoulli factory framework provides a concrete setting in which quantum resources can outperform classical ones. We experimentally demonstrate quantum randomness processing based on Bell-basis measurements of two identical input quoins prepared on IBM superconducting hardware. Using only the measurement outcomes (and no external classical randomness source), we realize the classically inconstructible Bernoulli doubling primitive f(p)=2p, and, as intermediate outputs from the same Bell-measurement statistics, an exact fair coin f(p)=1/2 and the classically inconstructible function f(p)=4p(1−p). We benchmark the measured output biases against ideal predictions and discuss the impact of device noise. Our results establish a simple resource-efficient experimental primitive for quantum-to-classical randomness processing and support the viability of quantum Bernoulli factories for quantum-enhanced stochastic simulation and sampling tasks.

Roy, Tanay [Fermilab] (ORCID:000000019442862X)↗

Identifications of RR Lyrae Stars and Quasars from the Simulated Data of Mephisto-W Survey

We have investigated the feasibilities and accuracies of the identifications of RR Lyrae stars and quasars from the simulated data of the Multi-channel Photometric Survey Telescope (Mephisto) W Survey. Based on the variable sources light curve libraries from the Sloan Digital Sky Survey (SDSS) Stripe 82 data and the observation history simulation from the Mephisto-W Survey Scheduler, we have simulated the uvgriz multi-band light curves of RR Lyrae stars, quasars and other variable sources for the first-year observation of Mephisto W Survey. We have applied the ensemble machine learning algorithm Random Forest Classifier (RFC) to identify RR Lyrae stars and quasars, respectively. We build training and test samples and extract ~150 features from the simulated light curves and train two RFCs respectively for the RR Lyrae star and quasar classification. We find that, our RFCs are able to select the RR Lyrae stars and quasars with remarkably high precision and completeness, with purity = 95.4% and completeness = 96.9% for the RR Lyrae RFC and purity = 91.4% and completeness = 90.2% for the quasar RFC. In conclusion, we have also derived relative importances of the extracted features utilized to classify RR Lyrae stars and quasars.

(galaxies:) quasars: general↗

Randomized Federated Learning Methods for Nonsmooth, Nonconvex, and Hierarchical Optimization (Final Technical Report)

This final technical report summarizes the outcomes of a DOE-funded project on federated scientific machine learning (FL) under nonsmooth, nonconvex, and hierarchical optimization settings. The project develops new mathematical models, algorithms, and theoretical guarantees for decentralized stochastic, bilevel, and minimax optimization problems arising in DOE mission-relevant applications. A unified framework of randomized and zeroth-order federated optimization methods is introduced, providing provable convergence, communication efficiency, and sample-complexity guarantees. The report documents algorithmic design, theoretical analysis, and empirical validation of the proposed federated learning methods. The project also contributes to workforce development through graduate training and dissemination of results via publications and seminars.

97 MATHEMATICS AND COMPUTING↗

Posiform planting: generating QUBO instances for benchmarking

We are interested in benchmarking both quantum annealing and classical algorithms for minimizing quadratic unconstrained binary optimization (QUBO) problems. Such problems are NP-hard in general, implying that the exact minima of randomly generated instances are hard to find and thus typically unknown. While brute forcing smaller instances is possible, such instances are typically not interesting due to being too easy for both quantum and classical algorithms. In this contribution, we propose a novel method, called posiform planting , for generating random QUBO instances of arbitrary size with known optimal solutions, and use those instances to benchmark the sampling quality of four D-Wave quantum annealers utilizing different interconnection structures (Chimera, Pegasus, and Zephyr hardware graphs) and the simulated annealing algorithm. Posiform planting differs from many existing methods in two key ways. It ensures the uniqueness of the planted optimal solution, thus avoiding groundstate degeneracy, and it enables the generation of QUBOs that are tailored to a given hardware connectivity structure, provided that the connectivity is not too sparse. Posiform planted QUBOs are a type of 2-SAT boolean satisfiability combinatorial optimization problems. Our experiments demonstrate the capability of the D-Wave quantum annealers to sample the optimal planted solution of combinatorial optimization problems with up to 5, 627 qubits.

97 MATHEMATICS AND COMPUTING↗

Deep sea tides determination from GEOS-3

GEOS 3 altimeter data in a 5 degree X 5 degree square centered at 30 deg N, 70 deg W were analyzed to evaluate deep sea tide determination from a spacecraft. The signal to noise ratio of known tidal variability to altimeter measurement of sea level above the ellipsoid was 0.1. A sample was obtained in a 5 deg x 5 deg area approximately once every four days. The randomly spaced time series was analyzed using two independent least squares techniques.

Maul, G. A.↗

The graphics and data acquisition software package

A software package was developed for use with micro and minicomputers, particularly the LSI-11/DPD-11 series. The package has a number of Fortran-callable subroutines which perform a variety of frequently needed tasks for biomedical applications. All routines are well documented, flexible, easy to use and modify, and require minimal programmer knowledge of peripheral hardware. The package is also economical of memory and CPU time. A single subroutine call can perform any one of the following functions: (1) plot an array of integer values from sampled A/D data, (2) plot an array of Y values versus an array of X values; (3) draw horizontal and/or vertical grid lines of selectable type; (4) annotate grid lines with user units; (5) get coordinates of user controlled crosshairs from the terminal for interactive graphics; (6) sample any analog channel with program selectable gain; (7) wait a specified time interval, and (8) perform random access I/O of one or more blocks of a sequential disk file. Several miscellaneous functions are also provided.

Crosier, W. G.↗

Improved importance sampling technique for efficient simulation of digital communication systems

A new, improved importance sampling (IIS) approach to simulation is considered. Some basic concepts of IS are introduced, and detailed evolutions of simulation estimation variances for Monte Carlo (MC) and IS simulations are given. The general results obtained from these evolutions are applied to the specific previously known conventional importance sampling (CIS) technique and the new IIS technique. The derivation for a linear system with no signal random memory is considered in some detail. For the CIS technique, the optimum input scaling parameter is found, while for the IIS technique, the optimum translation parameter is found. The results are generalized to a linear system with memory and signals. Specific numerical and simulation results are given which show the advantages of CIS over MC and IIS over CIS for simulations of digital communications systems.

Lu, Dingqing↗

The computation of pi to 29,360,000 decimal digits using Borweins' quartically convergent algorithm

The quartically convergent numerical algorithm developed by Borwein and Borwein (1987) for 1/pi is implemented via a prime-modulus-transform multiprecision technique on the NASA Ames Cray-2 supercomputer to compute the first 2.936 x 10 to the 7th digits of the decimal expansion of pi. The history of pi computations is briefly recalled; the most recent algorithms are characterized; the implementation procedures are described; and samples of the output listing are presented. Statistical analyses show that the present decimal expansion is completely random, with only acceptable numbers of long repeating strings and single-digit runs.

Bailey, David H.↗

Nonlinear probabilistic finite element models of laminated composite shells

A probabilistic finite element analysis procedure for laminated composite shells has been developed. A total Lagrangian finite element formulation, employing a degenerated 3-D laminated composite shell with the full Green-Lagrange strains and first-order shear deformable kinematics, forms the modeling foundation. The first-order second-moment technique for probabilistic finite element analysis of random fields is employed and results are presented in the form of mean and variance of the structural response. The effects of material nonlinearity are included through the use of a rate-independent anisotropic plasticity formulation with the macroscopic point of view. Both ply-level and micromechanics-level random variables can be selected, the latter by means of the Aboudi micromechanics model. A number of sample problems are solved to verify the accuracy of the procedures developed and to quantify the variability of certain material type/structure combinations. Experimental data is compared in many cases, and the Monte Carlo simulation method is used to check the probabilistic results. In general, the procedure is quite effective in modeling the mean and variance response of the linear and nonlinear behavior of laminated composite shells.

Engelstad, S. P.↗

The massive halos of spiral galaxies

We use a sample of satellite galaxies to demonstrate the existence of extended massive dark halos around spiral galaxies. Isolated spirals with rotation velocities near 250 km/s have a typical halo mass within 200 kpc of 1.5-2.6 x 10(exp 12) solar mass (90% confidence range for H(sub 0) = 75 km/s/Mpc). This result is most easily derived using standard mass estimator techniques, but such techniques do not account for the strong observational selection effects in the sample, nor for the extended mass distributions that the data imply. These complications can be addressed using scale-free models similar to those previously employed to study binary galaxies. When satellite velocities are assumed isotropic, both methods imply massive and extended halos. However, the derived masses depend sensitively on the assumed shape of satellite orbits. Furthermore, both methods ignore the fact that many of the satellites in the sample have orbital periods comparable to the Hubble time. The orbital phases of such satellites cannot be random, and their distribution in radius cannot be freely adjusted; rather these properties reflect ongoing infall onto the outer halos of their primaries. We use detailed dynamical models for halo formation to evaluate these problems, and we devise a maximum likelihood technique for estimating the parameters of such models from the data. The most strongly constrained parameter is the mass within 200-300 kpc, giving the confidence limits quoted above. The eccentricity, e, of satellite orbits is also strongly constrained, 0.50 less than e less than 0.88 at 90% confidence, implying a near-isotropic distribution of satellite velocities. The cosmic density parameter in the vicinity of our isolated halos exceeds 0.13 at 90% confidence, with preferred values exceeding 0.3.

Zaritsky, Dennis↗

Star Tracker Performance Estimate with IMU

A software tool for estimating cross-boresight error of a star tracker combined with an inertial measurement unit (IMU) was developed to support trade studies for the Integrated Radio and Optical Communication project (iROC) at the National Aeronautics and Space Administration Glenn Research Center. Typical laser communication systems, such as the Lunar Laser Communication Demonstration (LLCD) and the Laser Communication Relay Demonstration (LCRD), use a beacon to locate ground stations. iROC is investigating the use of beaconless precision laser pointing to enable laser communication at Mars orbits and beyond. Precision attitude knowledge is essential to the iROC mission to enable high-speed steering of the optical link. The preliminary concept to achieve this precision attitude knowledge is to use star trackers combined with an IMU. The Star Tracker Accuracy (STAcc) software was developed to rapidly assess the capabilities of star tracker and IMU configurations. STAcc determines the overall cross-boresight error of a star tracker with an IMU given the characteristic parameters: quantum efficiency, aperture, apparent star magnitude, exposure time, field of view, photon spread, detector pixels, spacecraft slew rate, maximum stars used for quaternion estimation, and IMU angular random walk. This paper discusses the supporting theory used to construct STAcc, verification of the program and sample results.

Error Budget↗

Mapping Post-fire Conifer Regeneration using Snow-on Imagery

The 2000 Jasper Fire in the Black Hills of South Dakota was the largest wildfire to date in the region, burning over 83,000 acres of ponderosa pine forest. In collaboration with partners from the United States Forest Service (USFS) Black Hills Experimental Forest, USFS Rocky Mountain Research Station, and United States Geological Survey Geosciences and Environmental Change Science Center, we characterized post-fire forest regeneration within high severity burn patches. We accomplished this by implementing novel conifer detection techniques using a snow index mask to create a winter, snow-on image composite from Landsat 8 Operational Land Imager (OLI) and Sentinel-2 Multispectral Instrument (MSI) data. We utilized 2015 USFS stem maps of field-observed regeneration plots and ocularly sampled reforestation sites planted from 2001–2013. The field data and imagery were used to train a Random Forest (RF) model in Google Earth Engine. The RF model classified conifer regeneration density as low, medium, or high across the high-severity burn area with an overall accuracy of 81.3% for 2021. Approximately 45.9% of the high-severity burn area had low or no regeneration (0-40 trees per acre) 20 years post-fire. Given our partners' desire to find easily accessible low conifer regeneration zones, we identified 4,079 acres of priority planting sites that were within 1,500 feet of roads, had not been planted previously, and were larger than 50 acres. This method supports the use of snow-on imagery as a successful technique to identify conifer regeneration in a post-wildfire landscape.

Yeshey Seldon↗