Search NASA⌕ Search

SEARCH · Search NASA

Results for “adaptive 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Adaptive Metropolis Sampling with Product Distributions

The Metropolis-Hastings (MH) algorithm is a way to sample a provided target distribution pi(z). It works by repeatedly sampling a separate proposal distribution T(x,x') to generate a random walk {x(t)}. We consider a modification of the MH algorithm in which T is dynamically updated during the walk. The update at time t uses the {x(t' less than t)} to estimate the product distribution that has the least Kullback-Leibler distance to pi. That estimate is the information-theoretically optimal mean-field approximation to pi. We demonstrate through computer experiments that our algorithm produces samples that are superior to those of the conventional MH algorithm.

Wolpert, David H.↗
NASA NTRS2018

Reliability-Based Design of Thermal Protection Systems with Support Vector Machines

The primary objective of this work was to develop a computationally efficient and accurate approach to reliability analysis of thermal protection systems using support vector machines. An adaptive sampling approach was introduced informs a iterative support vector machine approximation of the limit state function used for measuring reliability. The proposed sampling approach efficient adds samples along the limit state function until the reliability approximation is converged. This methodology is applied to two samples, mathematical functions to test and demonstrate the applicability. Then, the adaptive sampling-based support vector machine approach is applied to the reliability analysis of a thermal protection system. The results of all three problems highlight the potential capability of the new approach in terms of accuracy and computational saving in determining thermal protection system reliability.

White, Laura M.↗
DOE OSTI2024

Two-Stage Estimation and Variance Modeling for Latency-Constrained Variational Quantum Algorithms

The quantum approximate optimization algorithm (QAOA) has enjoyed increasing attention in noisy, intermediate-scale quantum computing with its application to combinatorial optimization problems. QAOA has the potential to demonstrate a quantum advantage for NP-hard combinatorial optimization problems. As a hybrid quantum-classical algorithm, the classical component of QAOA resembles a simulation optimization problem in which the simulation outcomes are attainable only through a quantum computer. The simulation that derives from QAOA exhibits two unique features that can have a substantial impact on the optimization process: (i) the variance of the stochastic objective values typically decreases in proportion to the optimality gap, and (ii) querying samples from a quantum computer introduces an additional latency overhead. In this paper, we introduce a novel stochastic trust-region method derived from a derivative-free, adaptive sampling trust-region optimization method intended to efficiently solve the classical optimization problem in QAOA by explicitly taking into account the two mentioned characteristics. The key idea behind the proposed algorithm involves constructing two separate local models in each iteration: a model of the objective function and a model of the variance of the objective function. Exploiting the variance model allows us to restrict the number of communications with the quantum computer and also helps navigate the nonconvex objective landscapes typical in QAOA optimization problems. In conclusion, we numerically demonstrate the superiority of our proposed algorithm using the SimOpt library and Qiskit when we consider a metric of computational burden that explicitly accounts for communication costs.

Derivative-free Optimization↗
DOE OSTI2022

Efficient data acquisition and training of collisional-radiative model artificial neural network surrogates through adaptive parameter space sampling

Abstract Effective plasma transport modeling of magnetically confined fusion devices relies on having an accurate understanding of the ion composition and radiative power losses of the plasma. Generally, these quantities can be obtained from solutions of a collisional-radiative (CR) model at each time step within a plasma transport simulation. However, even compact, approximate CR models can be computationally onerous to evaluate, and in-situ evaluation of these models within a larger plasma transport code can lead to a rigid bottleneck. As a way to bypass this bottleneck, we propose deploying artificial neural network (ANN) surrogates to allow rapid evaluation of the necessary plasma quantities. However, one issue with training an accurate ANN surrogate is the reliance on a sufficiently large and representative training and validation data set, which can be time-consuming to generate. In this work we explore a data-driven active learning and training routine to allow autonomous adaptive sampling of the problem parameter space to ensure a sufficiently large and meaningful set of training data is assembled for the network training. As a result, we can demonstrate approximately order-of-magnitude savings in required training data samples to produce an accurate surrogate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Enabling Scale-Up Through Multi-Fidelity Adaptive Computing

We present ideas from our ongoing work in adaptive computing - an optimization framework that allows us to strategically deploy various fidelity level experiments and simulations to guide decision making. The framework aims to enable uncertainty quantified scale-up of simulations and experiments, which causes increased complexity. A key feature of the framework is the integration of user-specified local model trustworthiness estimates. Adaptive sampling strategies allow us to optimally exploit the multiple fidelity level information and trustworthiness measures to arrive at the best decisions within a highly limited budget of objective function evaluations.

adaptive sampling↗
NASA NTRS2013

Elucidating Microbial Adaptation Dynamics via Autonomous Exposure and Sampling

The adaptation of micro-organisms to their environments is a complex process of interaction between the pressures of the environment and of competition. Reducing this multifactorial process to environmental exposure in the laboratory is a common tool for elucidating individual mechanisms of evolution, such as mutation rates. Although such studies inform fundamental questions about the way adaptation and even speciation occur, they are often limited by labor-intensive manual techniques. Current methods for controlled study of microbial adaptation limit the length of time, the depth of collected data, and the breadth of applied environmental conditions. Small idiosyncrasies in manual techniques can have large effects on outcomes; for example, there are significant variations in induced radiation resistances following similar repeated exposure protocols. We describe here a project under development to allow rapid cycling of multiple types of microbial environmental exposure. The system allows continuous autonomous monitoring and data collection of both single species and sampled communities, independently and concurrently providing multiple types of controlled environmental pressure (temperature, radiation, chemical presence or absence, and so on) to a microbial community in dynamic response to the ecosystem's current status. When combined with DNA sequencing and extraction, such a controlled environment can cast light on microbial functional development, population dynamics, inter- and intra-species competition, and microbe-environment interaction. The project's goal is to allow rapid, repeatable iteration of studies of both natural and artificial microbial adaptation. As an example, the same system can be used both to increase the pH of a wet soil aliquot over time while periodically sampling it for genetic activity analysis, or to repeatedly expose a culture of bacteria to the presence of a toxic metal, automatically adjusting the level of toxicity based on the number or growth rate of surviving cells. We are on our second prototype iteration, with demonstrated functions of microbial growth monitoring and dynamic exposure to UV-C radiation and temperature. We plan to add functionality for general chemical presence or absence by Nov. 2013. By making the project low-cost and open-source, we hope to encourage others to use it as a basis for future development of a common microbial environmental adaptation testbed.

Microbiology↗
DOE OSTI2023

Predicting rare events using neural networks and short-trajectory data

Estimating the likelihood, timing, and nature of events is a major goal of modeling stochastic dynamical systems. When the event is rare in comparison with the timescales of simulation and/or measurement needed to resolve the elemental dynamics, accurate prediction from direct observations becomes challenging. In such cases a more effective approach is to cast statistics of interest as solutions to Feynman-Kac equations (partial differential equations). Here, we develop an approach to solve Feynman-Kac equations by training neural networks on short-trajectory data. Our approach is based on a Markov approximation but otherwise avoids assumptions about the underlying model and dynamics. This makes it applicable to treating complex computational models and observational data. Additionally, we illustrate the advantages of our method using a low-dimensional model that facilitates visualization, and this analysis motivates an adaptive sampling strategy that allows on-the-fly identification of and addition of data to regions important for predicting the statistics of interest. Finally, we demonstrate that we can compute accurate statistics for a 75-dimensional model of sudden stratospheric warming. This system provides a stringent test bed for our method.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2022

Adaptive Force Control for Small Celestial Body Sampling

An adaptive force control algorithm for small celestial body sampling for a variety of surface properties is developed. The control algorithm consists of an adaptive controller combined with feedback linearization. When a spacecraft makes contact with the surface, it must maintain a desired contact force in order to capture a sample. The properties of the surface are unknown or uncertain before contact with the surface is made. The adaptive controller performs system identification online to create an input-output model of the feedback linearized system. From the input-output model a block observable canonical form is realized and the control input is determined by model predictive control (MPC) to maintain a desired contact force in spite of the unknown surface properties. The approach is applied to a variety of surface properties with linear and nonlinear contact models.

Quadrelli, Marco M↗
DOE OSTI2024

Advanced multimaterial shape optimization methods as applied to advanced manufacturing of wind turbine generators

Abstract Currently, many utility‐scale wind turbine generator original equipment manufacturers are dependent on imported rare earth permanent magnets, which are susceptible to market risks from cost instability. To lower the production costs of these generators and stay competitive in the market, several small wind manufacturers are pursuing continuous improvements to both generator design and manufacturing. However, traditional design and manufacturing methods have yielded marginal improvements in wind power performance. This work presents novel methods to redesign a baseline 15‐kW wind turbine generator with reduced rare‐earth permanent magnets by leveraging cutting‐edge three‐dimensional (3D) printed polymer‐bonded permanent magnets and steel. Symmetric, asymmetric, and multimaterial‐magnet parametrization methods are introduced for shape optimization. We extend the symmetric and asymmetric methods to the back iron in the stator to further investigate the impact and opportunities for performance improvements with lesser active materials. We employ a design‐of‐experiments approach with parametric computer‐aided design for shape generation and evaluate different designs by magneto‐thermal modeling and finite‐element analysis. We use adaptive sampling technique to identify better performing designs with lesser magnet mass, higher efficiency, and lower cogging torque when compared with the baseline generator. Asymmetric pole designs resulted in a magnet mass in the range of 4.77–5.37 kg, which was 27%–35% lighter than the baseline generator, suggesting that a new design freedom exists that can be enabled by advanced manufacturing, such as 3D printing. Shaping the back iron in the stator resulted in material savings in electrical steel of up to 14.62 kg, which was 20% lighter than the baseline stator. We conducted a structural analysis to evaluate an optimized asymmetric rotor design from the point of view of mechanical integrity and air‐gap stiffness. The magnetically optimal shape profile was shown as having a positive impact on the radial stiffness, and an optimal solution was discovered to reduce the structural mass by nearly 30 kg, which was 29% lighter than the baseline.

17 WIND ENERGY↗
DOE OSTI2018

Modeling Mesoscale Processes of Scalable Synthesis Final Scientific/Technical Report

This is the final scientific and technical report for this award. Molecular dynamics (MD) is a powerful tool to study biomolecules at the nanoscale. One of its main limitations however is that such simulations typically require a timestep of 1 femtosecond while reactions and conformational changes typically happen on the time scale of 1 microsecond or more. This discrepancy severely restricts the impact of this method. We have developed a series of algorithms to mitigate this problem and allow predicting long-time scale events using much shorter simulations. The method is based on the concurrent adaptive sampling method. It consists in using probabilistic weights for trajectories and biasing the sampling such that rare but important trajectories, e.g., those that go over energy barriers and lead to relevant conformational changes in the molecule, are observed with much higher frequency that in a standard MD simulation. We will demonstrate below that this approach has been very successful. In collaboration with Jay Grate at PNNL we have applied this novel method to the study the triazine polymers, which can be con gured using different types of side chains in order to adopt different 3-dimensional structures. We have found novel stable con- formations and have analyzed the impact of side chain modifcations on the native folded shape of the polymer. Direct MD simulations would have been untractable for this problem. The method of replica exchange was also used successfully to calculate free energy profiles, in particular for cis-trans torsional angle changes with high energy barriers. Studying molecules at the atomistic scale is essential to design new molecules in chemical engineering, such as the triazine polymers developed by Grate, in material science to study the strength of metals for examples (e.g., dislocation dynamics), or in biochemistry to understand how proteins fold and how they interact with ligands (a key step in designing drugs). Although new hardware makes these simulations increasingly faster (e.g., GPU, Anton from DE Shaw Research), new algorithms are key to making these methods more tractable and increase their impact. At this point, because of limitations in modeling, one of the main techniques for advancing science remains to a large extent numerous and costly experiments. For exam- ple, when designing drug, where multiple variants must be screened, the ability to run tests on a computer rather than experimentally testing with high throughput methods thousands of compounds would represent a signi cant advance with broad benefits for the public at large.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A nonlinear data-driven reduced order model for computational homogenization with physics/pattern-guided sampling

Developing an accurate nonlinear reduced order model from simulation data has been an outstanding research topic for many years. For many physical systems, data collection is very expensive and the optimal data distribution is not known in advance. Thus, maximizing the information gain remains a grand challenge. In a recent paper, Bhattacharjee and Matouš (2016) proposed a manifold-based nonlinear reduced order model for multiscale problems in mechanics of materials. Expanding this work here, we develop a novel sampling strategy based on the physics/pattern-guided data distribution. Our adaptive sampling strategy relies on enrichment of sub-manifolds based on the principal stretches and rotational sensitivity analysis. This novel sampling strategy substantially decreases the number of snapshots needed for accurate reduced order model construction (i.e., ~5 x reduction of snapshots over Bhattacharjee and Matouš (2016)). Moreover, we build the nonlinear manifold using the displacement rather than deformation gradient data. We provide rigorous verification and error assessment. We demonstrate both localization and homogenization of the multiscale solution on a large particulate composite unit cell.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Plug-and-play polymer microfluidic chips for hydrated, room temperature, fixed-target serial crystallography

The practice of serial X-ray crystallography (SX) depends on efficient, continuous delivery of hydrated protein crystals while minimizing background scattering. Of the two major types of sample delivery devices, fixed-target devices offer several advantages over widely adopted jet injectors, including: lower sample consumption, clog-free delivery, and the ability to control on-chip crystal density to improve hit rates. In this study we present our development of versatile, inexpensive, and robust polymer microfluidic chips for routine and reliable room temperature serial measurements at both synchrotrons and X-ray free electron lasers (XFELs). Our design includes highly X-ray-transparent enclosing thin film layers tuned to minimize scatter background, adaptable sample flow layers tuned to match crystal size, and a large sample area compatible with both raster scanning and rotation based serial data collection. The optically transparent chips can be used both for in situ protein crystallization (to eliminate crystal handling) or crystal slurry loading, with prepared samples stable for weeks in a humidified environment and for several hours in ambient conditions. Serial oscillation crystallography, using a multi-crystal rotational data collection approach, at a microfocus synchrotron beamline (SSRL, beamline 12-1) was used to benchmark the performance of the chips. Furthermore, high-resolution structures (1.3–2.7 Å) were collected from five different proteins – hen egg white lysozyme, thaumatin, bovine liver catalase, concanavalin-A (type VI), and SARS-CoV-2 nonstructural protein NSP5. Overall, our modular fabrication approach enables precise control over the cross-section of materials in the X-ray beam path and facilitates chip adaption to different sample and beamline requirements for user-friendly, straightforward diffraction measurements at room temperature.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Demonstration of an AI-driven workflow for dynamic x-ray spectroscopy

X-ray absorption near edge structure (XANES) spectroscopy is a powerful technique for characterizing the chemical state and symmetry of individual elements within materials, but requires collecting data at many energy points which can be time-consuming. While adaptive sampling methods exist for efficiently collecting spectroscopic data, they often lack domain-specific knowledge about the structure of XANES spectra. Here we demonstrate a knowledge-injected Bayesian optimization approach for adaptive XANES data collection that incorporates understanding of spectral features like absorption edges and pre-edge peaks. We show this method accurately reconstructs the absorption edge of XANES spectra using only 15–20% of the measurement points typically needed for conventional sampling, while maintaining the ability to determine the x-ray energy of the sharp peak after the absorption edge with errors less than 0.03 eV, the absorption edge with errors less than 0.1 eV; and overall root-mean-square errors less than 0.005 compared to traditionally sampled spectra. Our experiments on battery materials and catalysts demonstrate the method’s effectiveness for both static and dynamic XANES measurements, improving data collection efficiency and enabling better time resolution for tracking chemical changes. This approach advances the degree of automation in XANES experiments, reducing the common errors of under- or over-sampling points near the absorption edge and enabling dynamic experiments that require high temporal resolution or limited measurement time.

Bayesian optimization↗
NASA NTRS1985

Photoconversion from the light-adapted to the dark-adapted state of bacteriorhodopsin

The dark and light adaptation of the bR(trans)570 (bacteriorhodopsin) and bR(cis)550 isomers is analyzed. The equilibrium between the two bR isomers in light-adapted purple membrane films is studied in terms of the wavelength of actinic light and hydration levels. Absorption spectra observed after light adaptations with red and yellow light reveal that red light is less efficient in converting bR(cis)550 to bR(trans)570 than yellow light and the amount of bR(cis)550 in a light-adapted sample increases with decreasing hydration. The rate constants of dark and light adaptation are evaluated; the rate constant of dark adaptation is independent of the hydration level and the rate constant of light adaptation increases with hydration. The acceleration of a dark adaptation by red light is investigated; the dependence of the accelerated dark adaptation on the light intensity is discussed. The action spectrum of light adaptation in a purple membrane suspension is compared with the absorption spectrum of bR(cis)550; correlation between the spectra reveals that cis-to-transconversion is due to excitation of bR(cis)550 and the mechanism of cis-to-trans conversion in film is not affected by humidity levels. It is noted that the light-driven trans-to-cis conversion is a single photon process. The branching at M410 from the all-trans into the 13-cis photocycle is examined.

Kouyama, T.↗
DOE OSTI2018

Investigating the role of non-covalent interactions in conformation and assembly of triazine-based sequence-defined polymers

Grate and co-workers at Pacific Northwest National Laboratory recently developed high information content triazine-based sequence-defined polymers that are robust by not having hydrolyzable bonds and can encode structure and functionality by having various side chains. Through molecular dynamics (MD) simulations, the triazine polymers have been shown to form particular sequential stacks, have stable backbone-backbone interactions through hydrogen bonding and π–π interactions, and conserve their cis/trans conformations throughout the simulation. However, we do not know the effects of having different side chains and backbone structures on the entire conformation and whether the cis or trans conformation is more stable for the triazine polymers. For this reason, we investigate the role of non-covalent interactions for different side chains and backbone structures on the conformation and assembly of triazine polymers in MD simulations. Since there is a high energy barrier associated with the cis-trans isomerization, we use replica exchange molecular dynamics (REMD) to sample various conformations of triazine hexamers. To obtain rates and intermediate conformations, we use the recently developed concurrent adaptive sampling (CAS) algorithm for dimers of triazine trimers. Here, we found that the hydrogen bonding ability of the backbone structure is critical for the triazine polymers to self-assemble into nanorod-like structures, rather than that of the side chains, which can help researchers design more robust materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Foldamer architectures of triazine-based sequence-defined polymers investigated with molecular dynamics simulations and enhanced sampling methods

Triazine-based sequence-defined polymers have recently been developed that are biomimetic and robust (Grate \latin{et al.} Angewandte Chemie 2016). In molecular dynamics (MD) simulations, the triazine polymers were shown to form linear nanorod foldamers through hydrogen bonding and $\pi$-$\pi$ interactions. The nanorod foldamers have motifs resembling those of DNA, $\alpha$-helices, and $\beta$-sheets, and have potential to be useful building blocks for new macromolecules and materials. To understand the formation of nanorod foldamers, we investigate how linker structures in the middle of the triazine polymers lead to folding using MD simulations. We found that a variety of linkers can participate in folding, but that specific linker structures are more favorable than others, depending on the polymer length. Folding of hexamers into well-defined nanorod foldamers was most favorable with pentanediamine and \emph{ortho}-xylenediamine linkers in the center of the polymers. Foldamers with \emph{ortho}-xylenediamine linkers in the center were investigated for longer polymers, i.e., octamers and decamers, using two different enhanced sampling methods, since regular MD simulations had failed to show any folding for these longer polymers. In particular, the recently developed concurrent adaptive sampling (CAS) algorithm (Ahn \latin{et al.} J. Chem. Phys. 2017) and replica exchange molecular dynamics (REMD) were used. We found that the two enhanced sampling methods did lead to the observation of foldamers, and that REMD revealed new foldamer architectures where \emph{cis}-\emph{trans} isomerizations had occurred. Foldamer formation, diversity, and the strengths and limitations of simulation techniques are discussed. These findings provide new insights into the diversity of foldamer architectures for a new type of biomimetic synthetic polymer.

Foldamers, Nanorods, Free energy landscapes, Molec↗
NASA NTRS1990

A new finite element method for solving compressible Navier-Stokes equations based on an operator splitting method and h-p adaptivity

A new finite element method solving compressible Navier-Stokes equations is proposed. The method is based on a version of Strang's operator splitting and an h-p adaptive finite element approximation in space. This paper contains the formulation of the method with a detailed discussion of boundary conditions, a sample adaptive strategy and numerical examples involving compressible viscous flow over a flat plate with Reynolds number Re = 1000 and Re = 10,000.

Demkowicz, L.↗