Search NASA⌕ Search

SEARCH · Search NASA

Results for “manifold 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 19 records

Structure-aware methods for expensive derivative-free nonsmooth composite optimization

We present new methods for solving a broad class of bound-constrained nonsmooth composite minimization problems. These methods are specially designed for objectives that are some known mapping of outputs from a computationally expensive function. We provide accompanying implementations of these methods: in particular, a novel manifold sampling algorithm (MS-P) with subproblems that are in a sense primal versions of the dual problems solved by previous manifold sampling methods and a method (GOOMBAH) that employs more difficult optimization subproblems. For these two methods, we provide rigorous convergence analysis and guarantees. We demonstrate extensive testing of these methods. Open-source implementations of the methods developed in this manuscript can be found at https://github.com/POptUS/ IBCDFO/.

97 MATHEMATICS AND COMPUTING↗

Airborne laser induced fluorescence system for measuring OH and other trace gases in the parts-per-quadrillion to parts-per-trillion range

Described in detail is a laser induced fluorescence system which has been successfully interfaced with two aircraft sampling platforms (i.e., Sabreliner jet and an L-188C Electra). This system, which has been under development for four years, presently consists of the following major components: (1) a Nd-Yag laser driven oscillator-amplifier dye laser; (2) a sampling manifold with associated fluorescence detection optics; (3) an OH calibration chamber; (4) a laser beam steering assembly; and (5) sampling electronics and data processing hardware. During the last three years, this system has been flown some 50,000 air miles making tropospheric OH radical measurements over the latitude range of 70 N to 57 S. OH concentrations measured during these flights have ranged from 30 parts-per-quadrillion (3.7 x 10 to the 5th molecules/sq cm) at altitudes of 6 km to 0.8 parts-per-trillion (2.0 x 10 to the 7th molecules/sq cm) at 0.5 km. Computations have been completed which indicate that the existing aircraft system with modest modifications should also be capable of detecting natural tropospheric levels of NO, SO2, CH2O, NO2, HNO2, NO3, H2O2 and CS2 by using both conventional laser-induced fluorescence methodology and multiphoton techniques.

Davis, D. D.↗

Optimization on Manifolds via Graph Gaussian Processes

This paper integrates manifold learning techniques within a Gaussian process upper confidence bound algorithm to optimize an objective function on a manifold. Our approach is motivated by applications where a full representation of the manifold is not available and querying the objective is expensive. We rely on a point cloud of manifold samples to define a graph Gaussian process surrogate model for the objective. Query points are sequentially chosen using the posterior distribution of the surrogate model given all previous queries. We establish regret bounds in terms of the number of queries and the size of the point cloud. Several numerical examples complement the theory and illustrate the performance of our method.

Bayesian optimization↗

Development of an Atmospheric Pressure Ionization Mass Spectrometer

A commercial atmospheric pressure ionization mass spectrometer (APIMS) was purchased from EXTREL Mass Spectrometry, Inc. (Pittsburgh, PA). Our research objectives were to adapt this instrument and develop techniques for real-time determinations of the concentrations of trace species in the atmosphere. The prototype instrument is capable of making high frequency measurements with no sample preconcentrations. Isotopically labeled standards are used as an internal standard to obtain high precision and to compensate for changes in instrument sensitivity and analyte losses in the sampling manifold as described by Bandy and coworkers. The prototype instrument is capable of being deployed on NASA C130, Electra, P3, and DC8 aircraft. After purchasing and taking delivery by June 1994, we assembled the mass spectrometer, data acquisition, and manifold flow control instrumentation in electronic racks and performed tests.

Source record↗

Generative learning for slow manifolds and bifurcation diagrams

In dynamical systems characterized by separation of time scales, the approximation of so called “slow manifolds”, on which the long term dynamics lie, is a useful step for model reduction. Initializing on such slow manifolds is a useful step in modeling, since it circumvents fast transients, and is crucial in multiscale algorithms (like the equation-free approach) alternating between fine scale (fast) and coarser scale (slow) simulations. In a similar spirit, when one studies the infinite time dynamics of systems depending on parameters, the system attractors (e.g., its steady states) lie on bifurcation diagrams (curves for one-parameter continuation, and more generally, on manifolds in state parameter space. Sampling these manifolds gives us representative attractors (here, steady states of ODEs or PDEs) at different parameter values. Algorithms for the systematic construction of these manifolds (slow manifolds, bifurcation diagrams) are required parts of the “traditional” numerical nonlinear dynamics toolkit. In more recent years, as the field of Machine Learning develops, conditional score-based generative models (cSGMs) have been demonstrated to exhibit remarkable capabilities in generating plausible data from target distributions that are conditioned on some given label. It is tempting to exploit such generative models to produce samples of data distributions (points on a slow manifold, steady states on a bifurcation surface) conditioned on (consistent with) some quantity of interest (QoI, observable). In this work, we present a framework for using cSGMs to quickly (a) initialize on a low-dimensional (reduced-order) slow manifold of a multi-time-scale system consistent with desired value(s) of a QoI (a “label”) on the manifold, and (b) approximate steady states in a bifurcation diagram consistent with a (new, out-of-sample) parameter value. This conditional sampling can help uncover the geometry of the reduced slow-manifold and/or approximately “fill in” missing segments of steady states in a bifurcation diagram. Finally, the quantity of interest, which determines how the sampling is conditioned, is either known a priori or identified using manifold learning-based dimensionality reduction techniques applied to the training data.

Dynamical systems↗

Sample Processor for Life on Icy Worlds (SPLIce): Monolithic Manifold-Based System to Recover, Prepare, and Deliver Samples and Standards to Instrumentation Suites for Ocean World Life-Search Missions

A claim of life detection on one of the solar system’s icy ocean worlds would necessitate extraordinarily convincing evidence. Limited energy availability in the oceans of such bodies as Europa and Enceladus argues for microbes as most probable among possible life forms, but evidence of their existence in the surface layers of an icy moon or in a frozen plume ejected into space could take various forms, pointing to instrumentation suites as a preferred means to detect diverse molecular and morphological life indicators. Multiple disparate categories of positive detections could provide truly convincing evidence from samples that may be only a few micro-liters. SPLIce’s Foundation. Teams led by NASA’s Ames Research Center have developed and operated numerous small, live-biology and astrobiology science payloads in space over two decades. Since 2016, we have adapted and augmented their biological sample-handling systems to create compact, robust search-for-life fluidic processors designed to function after a decade or more in transit, in environments with very little gravity and lots of radiation: up to 100’s of kilorads.

Instrumentation suite↗

Beagle 2: The Next Exobiology Mission to Mars

Beagle 2 is a 60 kg probe (with a 30 kg lander) developed in the United Kingdom for inclusion on the European Space Agency's 2003 Mars Express. Beagle 2 will deliver to the Martian surface a payload which consists of a high percentage of science instruments to landed spacecraft mass. Beagle 2 will be launched in June, 2003 with Mars Express on a Soyuz-Fregat rocket from the Baikonur Cosmodrome in Kazakhstan. Beagle 2 will land on Mars on December 26, 2003 in the Isidis Planitia basin (approximately 10 degrees N and 275 degrees W), a large sedimentary basin that overlies the boundary between ancient highlands and northern plains. Isidis Planitia, the third largest basin on Mars, which is possibly filled with sediment deposited at the bottom of long-standing lakes or seas, offers an ideal environment for preserving traces of life. Beagle 2 was developed to search for organic material and other volatiles on and below the surface of Mars in addition to the study of the inorganic chemistry and mineralogy. Beagle 2 will utilize a mechanical mole and grinder to obtain samples from below the surface, under rocks and inside rocks. A pair of stereo cameras will image the landing site along with a microscope for examination of surface and rock samples. Analyses will include both rock and soil samples at various wavelengths, X-ray spectrometer and Mossbauer spectrometer as well as a search for organics and other light element species (e.g. carbonates and water) and measurement of their isotopic compositions. Beagle 2 has as its focus the goal of establishing whether evidence for life existed in the past on Mars at the Isidis Planitia site or at least establishing if the conditions were ever suitable. Carbonates and organic components were first recognized as existing on Mars when they were found in the Martian meteorite Nakhla. Romanek et al showed the carbonates in ALH84001 were formed at low temperatures. McKay et al noted possible evidence of early life on Mars within the ALH84001 meteorite. Thomas-Keprta et al showed the magnetite biomarkers in ALH84001's carbonates are indistinguishable from those formed by magnetotactic bacteria found on Earth. Gibson et al showed there was significant evidence for liquid water and biogenic products present on Mars across a 3.9 billion year period. A mechanical arm (PAW) operates from the lander and is used for science operations along with sample acquisition. Instruments attached to the PAW include: stereo cameras, Mossbauer instrument, X-ray fluorescence instrument, microscope, environmental sensors, rock corer/grinder, a spoon, mirror, brushes, a mole attachment for acquisition of subsurface to depths of 1 to 2 meters and an illumination device. Each camera has 14 filters which have been optimized for mineralogy composition, dust and water vapor detection. The microscope's camera is designed for viewing the size and shape of dust particles, rock surfaces, microfossils, and characteristics of the samples prior to introduction into the gas analysis package (GAP). The camera has a resolution of 4 microns/pixel, features four color capability (red, green, blue and UV (ultraviolet) fluoresence), a depth of focus of 40 micrometers and translation stage of +3 millimeters. The heart of the Beagle 2's life detection package is the gas analysis package (GAP), which consists of a mass spectrometer with collectors at fixed masses for precise isotopic ratio measurements and voltage scanning for spectral analysis. Primary aim of the GAP is to search for the presence of bulk constituents, individual species, and isotopic fractionations for both extinct and extent life along with studying the low-temperature geochemistry of the hydrogen, carbon, nitrogen and oxygen components on Mars from both the surface and atmosphere. GAP is a magnetic sector mass spectrometer with the range of 1 to 140 amu which can be operated in both the static and dynamic modes. A triple Faraday collector array will be used for C, N and O ratios along with a double Faraday array for H/D. Pulse counting electron multiplier will be utilized for noble gases and selected organics. Anticipated detection limits are at the picomole level for operation in the static mode of operation and high precision isotopic measurements will be made in the dynamic mode. Sample processing and preparation system consists of reaction vessels along with references. Sample ovens capable of being heated are attached to the manifold for sample combustion. Surface, subsurface materials and interior rock specimens will be combusted in pure oxygen gas at various temperature intervals to release organic matter and volatiles. Combustion process will permit detection of all forms and all atoms of carbon present in the samples. A chemical processing system is capable of a variety of conversion reactions. Gases are manipulated either by cryogenic or chemical reactions and passed through the gas handling portion of the vacuum system. There are two modes of operation: quantitative analysis and precise isotopic measurements. Additional information is contained in the original extended abstract.

Gibson, Everett K., Jr.↗

Generative learning of densities on manifolds

A generative modeling framework is proposed that combines diffusion models and manifold learning to efficiently sample data densities on manifolds. The approach utilizes Diffusion Maps to uncover possible low-dimensional underlying (latent) spaces in the high-dimensional data (ambient) space. Two approaches for sampling from the latent data density are described. The first is a score-based diffusion model, which is trained to map a standard normal distribution to the latent data distribution using a neural network. The second one involves solving an Itô stochastic differential equation in the latent space. Additional realizations of the data are generated by lifting the samples back to the ambient space using Double Diffusion Maps , a recently introduced technique typically employed in studying dynamical system reduction; here the focus lies in sampling densities rather than system dynamics. The proposed approaches enable sampling high dimensional data densities restricted to low-dimensional, a priori unknown manifolds. The efficacy of the proposed framework is demonstrated through a benchmark problem and a material with multiscale structure.

Double diffusion maps↗

Observations of the J = 10 manifold of the pure rotational band of phosphine on Saturn

Saturn was observed in the vicinity of the J = 10 manifold of the pure rotational band of phosphine on 1984 July 10 and 12 from NASA's Kuiper Airborne Observatory with the facility far-infrared cooled grating spectrometer. On each night observations of the full disk plus rings were made at 4 to 6 discrete wavelengths which selectively sampled the manifold and the adjacent continuum. The previously reported detection of this manifold is confirmed. After subtraction of the flux due to the rings, the data are compared with disk-averaged models of Saturn. It is found that PH3 must be strongly depleted above the thermal inversion (approx. 70 mbar). The best fitting models consistent with other observational constaints indicate that PH3 is significantly depleted at even deeper atmospheric levels ( or = 500 mbar), implying an eddy diffusion coefficient for Saturn of 10 to the 4 cm sq/sec.

Haas, M. R.↗

Observations of the J = 10 manifold of the pure rotational band of phosphine on Saturn

Saturn was observed in the vicinity of the J = 10 manifold of the pure rotational band of phosphine on 1984 July 10 and 12 from NASA's Kuiper Airborne Observatory with the facility far-infrared cooled grating spectrometer. On each night observations of the full disk plus rings were made at 4 to 6 discrete wavelengths which selectively sampled the manifold and the adjacent continuum. The previously reported detection of this manifold is confirmed. After subtraction of the flux due to the rings, the data are compared with disk-averaged models of Saturn. It is found that PH3 must be strongly depleted above the thermal inversion (approx. 70 mbar). The best fitting models consistent with other observational constraints indicate that PH3 is significantly depleted at even deeper atmospheric levels (or = 500 mbar), implying an eddy diffusion coefficient for Saturn of 10 to the 4 cm sq/sec.

Haas, M. R.↗

Economic gas chromatograph system for subambient pressure gas sampling

Gas chromatograph sampling system consists of a manifold with a gas-sample valve, a minimum-volume pressure transducer with a portable monitor, a vacuum-source valve, and a sample inlet valve. Increased accuracy of analysis is obtained by better control of sample size.

Mitchell, S. M.↗

Physics-inspired spatiotemporal-graph AI ensemble for the detection of higher order wave mode signals of spinning binary black hole mergers

We present a new class of AI models for the detection of quasi-circular, spinning, non-precessing binary black hole mergers whose waveforms include the higher order gravitational wave modes ($\ell$, |m|) = {(2,2), (2,1), (3,3), (3,2), (4,4)}, and mode mixing effects in the $\ell$ = 3, |m| = 2 harmonics. These AI models combine hybrid dilated convolution neural networks to accurately model both short- and long-range temporal sequential information of gravitational waves; and graph neural networks to capture spatial correlations among gravitational wave observatories to consistently describe and identify the presence of a signal in a three detector network encompassing the Advanced LIGO and Virgo detectors. We first trained these spatiotemporal-graph AI models using synthetic noise, using 1.2 million modeled waveforms to densely sample this signal manifold, within 1.7 h using 256 NVIDIA A100 GPUs in the Polaris supercomputer at the Argonne Leadership Computing Facility. This distributed training approach exhibited optimal classification performance, and strong scaling up to 512 NVIDIA A100 GPUs. With these AI ensembles we processed data from a three detector network, and found that an ensemble of 4 AI models achieves state-of-the-art performance for signal detection, and reports two misclassifications for every decade of searched data. We distributed AI inference over 128 GPUs in the Polaris supercomputer and 128 nodes in the Theta supercomputer, and completed the processing of a decade of gravitational wave data from a three detector network within 3.5 h. Finally, we fine-tuned these AI ensembles to process the entire month of February 2020, which is part of the O3b LIGO/Virgo observation run, and found 6 gravitational waves, concurrently identified in Advanced LIGO and Advanced Virgo data, and zero false positives. This analysis was completed in one hour using one NVIDIA A100 GPU.

79 ASTRONOMY AND ASTROPHYSICS↗

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗

An intercomparison of nitric oxide measurement techniques

Results from an intercomparison of techniques to measure tropospheric levels of nitric oxide (NO) are discussed. The intercomparison was part of the National Aeronautics and Space Administration's Global Tropospheric Experiment and was conducted at Wallops Island, VA, in July 1983. Instruments intercompared included a laser-induced fluorescence system and two chemiluminescence instruments. The intercomparisons were performed with ambient air at NO mixing ratios ranging from 10 to 60 pptv and NO-enriched ambient air at mixing ratios from 20 to 170 pptv. All instruments sampled from a common manifold. The techniques exhibited a high degree of correlation among themselves and with changes in the NO mixing ratio. Agreement among the three techniques was placed at approximately + or - 30 percent. Within this level of agreement, no artifacts or species interferences were identified.

Hoell, J. M., Jr.↗

Generative models on phase space

Deep generative models such as diffusion and flow matching are powerful machine learning tools capable of learning and sampling from high-dimensional distributions. They are particularly useful when the training data appears to be concentrated on a submanifold of the data embedding space. For high-energy physics data, consisting of collections of relativistic energy-momentum 4-vectors, this submanifold can enforce extremely strong physically-motivated priors, such as energy and momentum conservation. If these constraints are learned only approximately, rather than exactly, this can inhibit the interpretability and reliability of such generative models. To remedy this deficiency, we introduce generative models which are, by construction, confined at every step of their sampling trajectory to the manifold of massless N-particle Lorentz-invariant phase space in the center-of-momentum frame. In the case of diffusion models, the "pure noise" forward process endpoint corresponds to the uniform distribution on phase space, which provides a clear starting point from which to identify how correlations among the particles emerge during the reverse (de-noising) process. We demonstrate that our models are able to learn both few-particle and many-particle distributions with various singularity structures, paving the way for future interpretability studies using generative models trained on simulated jet data.

Bogorad, Zachary [Fermilab]↗

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation) [SWR-26-095]

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation): Multifidelity aerodynamic polar data generation for hydrofoil/tidal-turbine airfoil sections. Foilpolars ties together three pieces: *AeroSandbox supplies the baseline airfoil coordinates (UIUC database). *G2Aero parameterizes those shapes on a Grassmannian manifold (Karcher mean + PGA basis) and samples new perturbed shapes around that basis. *XFoil (panel method) and NeuralFoil (neural-network surrogate, shipped with AeroSandbox) each solve the resulting shapes for lift, drag, moment, and pressure at the swept angles of attack, Reynolds numbers, and n_crit values. Design optimization of foil shapes in a computationally efficient way requires polars data across many candidate shapes, not just a handful of baseline foils. However, high-fidelity CFD at that scale is too costly, and naive shape perturbation strays from realistic geometries. FOILPOLARS addresses this by loading baseline airfoils (via AeroSandbox) and mapping them onto a Grassmannian manifold (via G2Aero), computing a Karcher mean and principal geodesic analysis (PGA) basis. New shapes are sampled by perturbing PGA coefficients, keeping them close to the manifold of realistic foils. Each sampled shape is evaluated across a configurable sweep of angle of attack, Reynolds number, and critical amplification factor using two solvers: XFoil (panel method) and NeuralFoil (neural-network surrogate), producing a paired dataset of lift, drag, moment, pressure, convergence, and confidence, indexed alongside each shape's PGA coefficients and shared Grassmannian basis in a single xarray dataset. From this, FOILPOLARS produces convergence summaries and comparison plots per shape, Reynolds number, and n_crit. A command-line interface exposes each pipeline stage independently, supporting data-driven design, optimization, and machine-learning workflows for foils.

Sandhu, Rimple [National Laboratory of the Rockies↗

Sampling functions for geophysics

A set of spherical sampling functions is defined such that they are related to spherical-harmonic functions in the same way that the sampling functions of information theory are related to sine and cosine functions. An orderly distribution of (N + 1) squared sampling points on a sphere is given, for which the (N + 1) squared spherical sampling functions span the same linear manifold as do the spherical-harmonic functions through degree N. The transformations between the spherical sampling functions and the spherical-harmonic functions are given by recurrence relations. The spherical sampling functions of two arguments are extended to three arguments and to nonspherical reference surfaces. Typical applications of this formalism to geophysical topics are sketched.

Giacaglia, G. E. O.↗

Deep nonparametric estimation of intrinsic data structures by chart autoencoders: Generalization error and robustness

Autoencoders have demonstrated remarkable success in learning low-dimensional latent features of high-dimensional data across various applications. Assuming that data are sampled near a low-dimensional manifold, we employ chart autoencoders, which encode data into low-dimensional latent features on a collection of charts, preserving the topology and geometry of the data manifold. Our paper establishes statistical guarantees on the generalization error of chart autoencoders, and we demonstrate their denoising capabilities by considering n noisy training samples, along with their noise-free counterparts, on a d-dimensional manifold. By training autoencoders, we show that chart autoencoders can effectively denoise the input data with normal noise. We prove that, under proper network architectures, chart autoencoders achieve a squared generalization error in the order of n–$\frac{2}{d+2}$log 4 n, which depends on the intrinsic dimension of the manifold and only weakly depends on the ambient dimension and noise level. We further extend our theory on data with noise containing both normal and tangential components, where chart autoencoders still exhibit a denoising effect for the normal component. As a special case, our theory also applies to classical autoencoders, as long as the data manifold has a global parametrization. Furthermore, our results provide a solid theoretical foundation for the effectiveness of autoencoders, which is further validated through several numerical experiments.

97 MATHEMATICS AND COMPUTING↗