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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Design of partially supervised classifiers for multispectral image data

A partially supervised classification problem is addressed, especially when the class definition and corresponding training samples are provided a priori only for just one particular class. In practical applications of pattern classification techniques, a frequently observed characteristic is the heavy, often nearly impossible requirements on representative prior statistical class characteristics of all classes in a given data set. Considering the effort in both time and man-power required to have a well-defined, exhaustive list of classes with a corresponding representative set of training samples, this 'partially' supervised capability would be very desirable, assuming adequate classifier performance can be obtained. Two different classification algorithms are developed to achieve simplicity in classifier design by reducing the requirement of prior statistical information without sacrificing significant classifying capability. The first one is based on optimal significance testing, where the optimal acceptance probability is estimated directly from the data set. In the second approach, the partially supervised classification is considered as a problem of unsupervised clustering with initially one known cluster or class. A weighted unsupervised clustering procedure is developed to automatically define other classes and estimate their class statistics. The operational simplicity thus realized should make these partially supervised classification schemes very viable tools in pattern classification.

Jeon, Byeungwoo↗
DOE OSTI2026

Multidimensional scaling informed by F -statistic: Visualizing grouped microbiome data with inference

Multidimensional scaling (MDS) is a widely used dimensionality reduction technique in microbial ecology data analysis that captures the multivariate structure of the data while preserving pairwise distances between samples. While improvements in MDS have enhanced the ability to reveal group-specific data patterns, these MDS-based methods require prior assumptions for inference, limiting their application in general microbiome analysis. Here, in this study, we introduce a new MDS-based ordination method, “F-informed MDS,” which configures the data distribution based on the F-statistic, the ratio of dispersion between groups sharing common and different characteristics. Using semisynthetic datasets, we demonstrate that the proposed method is robust to hyperparameter selection while maintaining statistical significance throughout the ordination process. Various quality metrics for evaluating dimensionality reduction confirm that F-informed MDS is comparable to state-of-the-art methods in preserving both local and global data structures. Its application to a diatom-associated bacterial community suggests the role of this new method in interpreting the community’s response to the host. Our approach offers a well-founded refinement of MDS that aligns with statistical test results, which can be beneficial for broader multidimensional data analyses in microbiology and ecology. This new visualization tool can be incorporated into standard microbiome data analyses.

Biological and medical sciences↗
NASA NTRS1988

Quasar microlensing and dark matter

The amplification of quasar brightness due to gravitational lensing by foreground objects is discussed. It is shown that a recently published sample of X-ray-selected quasars behind foreground galaxies shows a statistically significant brightening compared to a control sample. Correlations with galaxy redshift and impact parameter predicted by microlensing are also demonstrated. A technique is described to measure the mean density of the lenses from a small number of identified cases of microlensing. It is shown that, in this sample, amplification bias is important in determining the mean intensity enhancement and must be included in the density estimate. Assuming that at least two of the four intrinsically brightest quasars behind galaxies are indeed microlensed, the present data yield a formal lower limit on the mean density parameter of lenses Omega(l) greater than 0.25 at 95 percent confidence. These data also imply that a considerable quantity of dark matter exists in macroscopic objects outside the visible parts of galaxies but is still highly correlated with them.

Rix, Hans-Walter↗
NASA NTRS1990

Potential, velocity, and density fields from sparse and noisy redshift-distance samples - Method

A method for recovering the three-dimensional potential, velocity, and density fields from large-scale redshift-distance samples is described. Galaxies are taken as tracers of the velocity field, not of the mass. The density field and the initial conditions are calculated using an iterative procedure that applies the no-vorticity assumption at an initial time and uses the Zel'dovich approximation to relate initial and final positions of particles on a grid. The method is tested using a cosmological N-body simulation 'observed' at the positions of real galaxies in a redshift-distance sample, taking into account their distance measurement errors. Malmquist bias and other systematic and statistical errors are extensively explored using both analytical techniques and Monte Carlo simulations.

Dekel, Avishai↗
NASA NTRS2004

Physical Validation of TRMM TMI and PR Monthly Rain Products Over Oklahoma

The Tropical Rainfall Measuring Mission (TRMM) provides monthly rainfall estimates using data collected by the TRMM satellite. These estimates cover a substantial fraction of the earth's surface. The physical validation of TRMM estimates involves corroborating the accuracy of spaceborne estimates of areal rainfall by inferring errors and biases from ground-based rain estimates. The TRMM error budget consists of two major sources of error: retrieval and sampling. Sampling errors are intrinsic to the process of estimating monthly rainfall and occur because the satellite extrapolates monthly rainfall from a small subset of measurements collected only during satellite overpasses. Retrieval errors, on the other hand, are related to the process of collecting measurements while the satellite is overhead. One of the big challenges confronting the TRMM validation effort is how to best estimate these two main components of the TRMM error budget, which are not easily decoupled. This four-year study computed bulk sampling and retrieval errors for the TRMM microwave imager (TMI) and the precipitation radar (PR) by applying a technique that sub-samples gauge data at TRMM overpass times. Gridded monthly rain estimates are then computed from the monthly bulk statistics of the collected samples, providing a sensor-dependent gauge rain estimate that is assumed to include a TRMM equivalent sampling error. The sub-sampled gauge rain estimates are then used in conjunction with the monthly satellite and gauge (without sub- sampling) estimates to decouple retrieval and sampling errors. The computed mean sampling errors for the TMI and PR were 5.9% and 7.796, respectively, in good agreement with theoretical predictions. The PR year-to-year retrieval biases exceeded corresponding TMI biases, but it was found that these differences were partially due to negative TMI biases during cold months and positive TMI biases during warm months.

Fisher, Brad L.↗
NASA NTRS1976

Monte Carlo investigation of thrust imbalance of solid rocket motor pairs

The Monte Carlo method of statistical analysis is used to investigate the theoretical thrust imbalance of pairs of solid rocket motors (SRMs) firing in parallel. Sets of the significant variables are selected using a random sampling technique and the imbalance calculated for a large number of motor pairs using a simplified, but comprehensive, model of the internal ballistics. The treatment of burning surface geometry allows for the variations in the ovality and alignment of the motor case and mandrel as well as those arising from differences in the basic size dimensions and propellant properties. The analysis is used to predict the thrust-time characteristics of 130 randomly selected pairs of Titan IIIC SRMs. A statistical comparison of the results with test data for 20 pairs shows the theory underpredicts the standard deviation in maximum thrust imbalance by 20% with variability in burning times matched within 2%. The range in thrust imbalance of Space Shuttle type SRM pairs is also estimated using applicable tolerances and variabilities and a correction factor based on the Titan IIIC analysis.

Sforzini, R. H.↗
DOE OSTI2026

CADIS and FW-CADIS Variance Reduction in Gamma Transport for Predicting Prompt Forensics Signatures

The goal of prompt nuclear forensics is to determine the characteristics of a nuclear detonation based on the signatures available almost immediately after the explosion. An important characteristic is the reaction time history (RTH), a measure of the device’s rate of neutron multiplication. The RTH can be estimated by observation of the gamma radiation emitted from the detonation, which can be detected directly or observed indirectly as Teller light. Gamma transport simulations used to predict these radiation fields are often modeled stochastically using the Monte Carlo N-Particle (MCNP) code, which can be a computationally demanding task due to the number of particle histories needed to achieve statistical convergence. In an attempt to improve the efficiency of these calculations, we evaluate two variance reduction techniques: Consistent Adjoint-Driven Importance Sampling (CADIS) and Forward-Weighted Consistent Adjoint-Driven Importance Sampling (FW-CADIS). These methods use a deterministically calculated adjoint flux to create weight windows and source biasing that guide MCNP sampling. We study the utility of CADIS and FW-CADIS for their use in MCNP gamma transport for nuclear forensics prediction simulations. Furthermore, the results demonstrate that both CADIS and FW-CADIS improve the accuracy for forensics-focused simulations, with CADIS being most beneficial in direct detection and FW-CADIS being ideal for computing a global Teller light source.

CADIS↗
NASA NTRS1974

On evaluating compliance with air pollution levels 'not to be exceeded more than once per year'

The point of view taken is that the Environmental Protection Agency (EPA) Air Quality Standards (AQS) represent conditions which must be made to exist in the ambient environment. The statistical techniques developed should serve as tools for measuring the closeness to achieving the desired quality of air. It is shown that the sampling frequency recommended by EPA is inadequate to meet these objectives when the standard is expressed as a level not to be exceeded more than once per year and sampling frequency is once every three days or less frequent.

Neustadter, H. E.↗
NASA NTRS1990

Probing the tides in interacting galaxy pairs

Detailed spectroscopic and imaging observations of colliding elliptical galaxies revealed unmistakable diagnostic signatures of the tidal interactions. It is possible to compare both the distorted luminosity distributions and the disturbed internal rotation profiles with numerical simulations in order to model the strength of the tidal gravitational field acting within a given pair of galaxies. Using the best-fit numerical model, one can then measure directly the mass of a specific interacting binary system. This technique applies to individual pairs and therefore complements the classical methods of measuring the masses of galaxy pairs in well-defined statistical samples. The 'personalized' modeling of galaxy pairs also permits the derivation of each binary's orbit, spatial orientation, and interaction timescale. Similarly, one can probe the tides in less-detailed observations of disturbed galaxies in order to estimate some of the physical parameters for larger samples of interacting galaxy pairs. These parameters are useful inputs to the more universal problems of (1) the galaxy merger rate, (2) the strength and duration of the driving forces behind tidally stimulated phenomena (e.g., starbursts and maybe quasi steller objects), and (3) the identification of long-lived signatures of interaction/merger events.

Borne, Kirk D.↗
DOE OSTI2022

The DESI N -body Simulation Project – II. Suppressing sample variance with fast simulations

ABSTRACT Dark Energy Spectroscopic Instrument (DESI) will construct a large and precise three-dimensional map of our Universe. The survey effective volume reaches $\sim 20\, h^{-3}\, \mathrm{Gpc}^{3}$. It is a great challenge to prepare high-resolution simulations with a much larger volume for validating the DESI analysis pipelines. AbacusSummit is a suite of high-resolution dark-matter-only simulations designed for this purpose, with $200\, h^{-3}\, \mathrm{Gpc}^{3}$ (10 times DESI volume) for the base cosmology. However, further efforts need to be done to provide a more precise analysis of the data and to cover also other cosmologies. Recently, the CARPool method was proposed to use paired accurate and approximate simulations to achieve high statistical precision with a limited number of high-resolution simulations. Relying on this technique, we propose to use fast quasi-N-body solvers combined with accurate simulations to produce accurate summary statistics. This enables us to obtain 100 times smaller variance than the expected DESI statistical variance at the scales we are interested in, e.g. $k \lt 0.3\, h\, \mathrm{Mpc}^{-1}$ for the halo power spectrum. In addition, it can significantly suppress the sample variance of the halo bispectrum. We further generalize the method for other cosmologies with only one realization in AbacusSummit suite to extend the effective volume ∼20 times. In summary, our proposed strategy of combining high-fidelity simulations with fast approximate gravity solvers and a series of variance suppression techniques sets the path for a robust cosmological analysis of galaxy survey data.

cosmology: theory↗
DOE OSTI2024

Cross sections for the formation of Rb84m,g, Rb83, and Rb82m in Sr86(d,x) reactions up to deuteron energies of 49 MeV: Competition between α-particle and multinucleon emission processes

Cross sections of Sr86(d,x) reactions leading to the products Rb84m,g, Rb83, and Rb82m were measured by the stacked-sample activation technique up to deuteron energies of 49 MeV. Nuclear model calculations were performed using the codes talys and empire, which combine the statistical, precompound, and direct interaction components. In all cases, the empire results were much higher than the talys calculation. Fairly good agreement was obtained between measured data and the talys calculation after some optimization of the input model parameters. Insight into competition between α-particle and multinucleon emission in the Y88 compound-nucleus system was also gained.

59 ≤ A ≤ 89↗
NASA NTRS2006

Proteomic Retrieval from Nucleic Acid Depleted Space-Flown Human Cells

Compared to experiments utilizing humans in microgravity, cell-based approaches to questions about subsystems of the human system afford multiple advantages, such as crew safety and the ability to achieve statistical significance. To maximize the science return from flight samples, an optimized method was developed to recover protein from samples depleted of nucleic acid. This technique allows multiple analyses on a single cellular sample and when applied to future cellular investigations could accelerate solutions to significant biomedical barriers to human space exploration. Cell cultures grown in American Fluoroseal bags were treated with an RNA stabilizing agent (RNAlater - Ambion), which enabled both RNA and immunoreactive protein analyses. RNA was purified using an RNAqueous(registered TradeMark) kit (Ambion) and the remaining RNA free supernatant was precipitated with 5% trichloroacetic acid. The precipitate was dissolved in SDS running buffer and tested for protein content using a bicinchoninic acid assay (1) (Sigma). Equal loads of protein were placed on SDS-PAGE gels and either stained with CyproOrange (Amersham) or transferred using Western Blotting techniques (2,3,4). Protein recovered from RNAlater-treated cells and stained with protein stain, was measured using Imagequant volume measurements for rectangles of equal size. BSA treated in this way gave quantitative data over the protein range used (Fig 1). Human renal cortical epithelial (HRCE) cells (5,6,7) grown onboard the International Space Station (ISS) during Increment 3 and in ground control cultures exhibited similar immunoreactivity profiles for antibodies to the Vitamin D receptor (VDR) (Fig 2), the beta isoform of protein kinase C (PKC ) (Fig 3), and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (Fig 4). Parallel immunohistochemical studies on formalin-fixed flight and ground control cultures also showed positive immunostaining for VDR and other biomarkers (Fig 5). These results are consistent with data from additional antigenic recovery experiments performed on human Mullerian tumor cells cultured in microgravity (8).

Hammond, D. K.↗
DOE OSTI2023

Structure optimization with stochastic density functional theory

Linear-scaling techniques for Kohn–Sham density functional theory are essential to describe the ground state properties of extended systems. Still, these techniques often rely on the localization of the density matrix or accurate embedding approaches, limiting their applicability. In contrast, stochastic density functional theory (sDFT) achieves linear- and sub-linear scaling by statistically sampling the ground state density without relying on embedding or imposing localization. In return, ground state observables, such as the forces on the nuclei, fluctuate in sDFT, making optimizing the nuclear structure a highly non-trivial problem. In this work, we combine the most recent noise-reduction schemes for sDFT with stochastic optimization algorithms to perform structure optimization within sDFT. We compare the performance of the stochastic gradient descent approach and its variations (stochastic gradient descent with momentum) with stochastic optimization techniques that rely on the Hessian, such as the stochastic Broyden–Fletcher–Goldfarb–Shanno algorithm. In conclusion, we further provide a detailed assessment of the computational efficiency and its dependence on the optimization parameters of each method for determining the ground state structure of bulk silicon with varying supercell dimensions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2007

GEO Population Estimates using Optical Survey Data

Optical survey data taken using the NASA Michigan Orbital Debris Survey Telescope (MODEST) gives us an opportunity to statistically sample faint object population in the Geosynchronous (GEO) and near-GEO environment. This paper will summarize the MODEST survey work that has been conducted by NASA since 2002, and will outline the techniques employed to arrive at the current population estimates in the GEO environment for dim objects difficult to detect and track using current systems in the Space Surveillance Network (SSN). Some types of orbits have a higher detection rate based on what parts of the GEO belt is being observed, a straightforward statistical technique is used to debias these observations to arrive at an estimate of the total population potentially visible to the telescope. The size and magnitude distributions of these fainter debris objects are markedly different from the catalogued population. GEO debris consists of at least two different populations, one which follows the standard breakup power law and one which has anomalously high Area-to-Mass Ratios (1 to approx. 30 square meters per kilogram; a sheet of paper = approx. 13 square meters per kilogram). The Inter-Agency Space Debris Coordination Committee (IADC) is investigating objects in GEO orbits with anomalously high Area-to-Mass Ratios (AMRs). The ESA Space Debris Telescope discovered this population and has and its properties of inclinations (0 to 30 degrees), changing eccentricities (0 and 0.6), and mean motions (approx. 1 rev), will be presented. The accepted interpretation of this orbital behavior is that solar radiation pressure drives the perturbations causing time varying inclinations and eccentricities. The orbital parameters are unstable for this population and thus difficult to predict. Their dim visual magnitudes and photometric variability make observations a challenge. The IADC has enlisted a series of observatories (participating institutions: University of Michigan/CTIO, Astronomical Institute University of Bern, Boeing LTS / AMOS, Keldysh Institute of Applied Mathematics) at different longitudes. Complete observational coverage over periods of days to months will provide a better understanding of the properties, such as solar radiation pressure effects on orbital elements, size, shape, attitude, color variations, and spectral characteristics. Results from recent observational programs will be summarized, and includes a description of the orbit elements prediction processes, a summary of the metric tracking performance, and some photometric characteristics of this class of debris.

Barker, Edwin S.↗
DOE OSTI2023

Discovering strongly lensed quasar candidates with catalogue-based methods from DESI Legacy Surveys

The Hubble tension, revealed by a ~5σ discrepancy between measurements of the Hubble-Lemaitre constant among observations of the early and local Universe, is one of the most significant problems in modern cosmology. In order to better understand the origin of this mismatch, independent techniques to measure H 0 , such as strong lensing time delays, are required. Notably, the sample size of such systems is key to minimising the statistical uncertainties and cosmic variance, which can be improved by exploring the datasets of large-scale sky surveys such as Dark Energy Spectroscopic Instrument (DESI). We identify possible strong lensing time-delay systems within DESI by selecting candidate multiply imaged lensed quasars from a catalogue of 24 440 816 candidate QSOs contained in the ninth data release of the DESI Legacy Imaging Surveys (DESI-LS). Using a friend-of-friends-like algorithm on spatial co-ordinates, our method generates an initial list of compact quasar groups. This list is subsequently filtered using a measure of the similarity of colours among a group’s members and the likelihood that they are quasars. A visual inspection finally selects candidate strong lensing systems based on the spatial configuration of the group members. We identified 620 new candidate multiply imaged lensed quasars (101 grade-A, 214 grade-B, 305 grade-C). This number excludes 53 known spectroscopically confirmed systems and existing candidate systems identified in other similar catalogues. When available, these new candidates will be further checked by combining the spectroscopic and photometric data from DESI.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

The OSU self-phased array for propagation measurements using the 11.7 GHz CTS beacon

A self phased array was developed for propagation measurements on an earth-space path. The 11.7 GHz CTS beacon was used as the signal source. The self phased array was used to measure angle of arrival as well as attenuation and scintillation statistics. The performance of the array is described, and sample data are presented. The tracking capability of the self phased array was also studied. This technique permits fully electronic, nonmechanical satellite tracking, thus simplifying unmanned operation and eliminating severe weather tracking constraints.

Theobold, D. M.↗
NASA NTRS2014

Optimizing Integrated Terminal Airspace Operations Under Uncertainty

In the terminal airspace, integrated departures and arrivals have the potential to increase operations efficiency. Recent research has developed geneticalgorithm- based schedulers for integrated arrival and departure operations under uncertainty. This paper presents an alternate method using a machine jobshop scheduling formulation to model the integrated airspace operations. A multistage stochastic programming approach is chosen to formulate the problem and candidate solutions are obtained by solving sample average approximation problems with finite sample size. Because approximate solutions are computed, the proposed algorithm incorporates the computation of statistical bounds to estimate the optimality of the candidate solutions. A proof-ofconcept study is conducted on a baseline implementation of a simple problem considering a fleet mix of 14 aircraft evolving in a model of the Los Angeles terminal airspace. A more thorough statistical analysis is also performed to evaluate the impact of the number of scenarios considered in the sampled problem. To handle extensive sampling computations, a multithreading technique is introduced.

air traffic optimization↗