Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic 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 127 records · Page 7

Electron-ion equilibration in superheated gold

Electron-ion equilibration dynamics in samples driven into superheated states with multi-eV electron temperatures represent a fundamental process in nonequilibrium physics. A clear picture of the evolving ion temperature is essential for understanding the strength of the electron-ion coupling. However, direct, model-independent measurements of the ion temperature in laser-irradiated samples have remained experimentally elusive. Using inelastic X-ray scattering with meV-resolution in gold samples, the ion response to ultrafast heating by the hot electrons can be resolved, quantifying the electron-ion equilibration dynamics. We report measurements revealing a strongly enhanced energy transfer rate compared to that in weakly excited gold. Moreover, we obtain a quasi-constant electron-ion coupling at the highly elevated electron temperatures. These results place new constraints on electron-ion energy transfer in warm dense matter and establish a path toward quantitative benchmarking of models for nonequilibrium dynamics.

Griffin, Travis D. [Univ. of Nevada, Reno, NV (Uni↗

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

Practical aspects of modeling aircraft dynamics from flight data

The purpose of parameter estimation, a subset of system identification, is to estimate the coefficients (such as stability and control derivatives) of the aircraft differential equations of motion from sampled measured dynamic responses. In the past, the primary reason for estimating stability and control derivatives from flight tests was to make comparisons with wind tunnel estimates. As aircraft became more complex, and as flight envelopes were expanded to include flight regimes that were not well understood, new requirements for the derivative estimates evolved. For many years, the flight determined derivatives were used in simulations to aid in flight planning and in pilot training. The simulations were particularly important in research flight test programs in which an envelope expansion into new flight regimes was required. Parameter estimation techniques for estimating stability and control derivatives from flight data became more sophisticated to support the flight test programs. As knowledge of these new flight regimes increased, more complex aircraft were flown. Much of this increased complexity was in sophisticated flight control systems. The design and refinement of the control system required higher fidelity simulations than were previously required.

Iliff, K. W.↗

Water Imbibition and Oil Recovery in Shale: Dynamics and Mechanisms Using Integrated Centimeter-to-Nanometer-Scale Imaging

Water imbibition, and the associated oil displacement, is an important process in shale oil reservoirs after hydraulic fracturing and in water-based enhanced oil recovery (EOR). Current techniques for water imbibition measurement are mostly “black-box”-type methods. A more explicit understanding of the water imbibition/oil recovery dynamics and geological controls is in demand. In this paper, a multiscale imaging technique that covers centimeter to nanometer scale (i.e., core to pore scale), integrating neutron radiography, microcomputed tomography (micro-CT), and scanning electron microscope (SEM) is applied to investigate the water imbibition depth and rate and the cause of heterogeneity of imbibition in shale samples. The dynamic processes of water imbibition in the 1-in. (25.4-mm) core sample were explicitly demonstrated, and the imbibition along the matrix and imbibition through microfractures are distinguished through neutron radiography image analysis. The causes of observed imbibition heterogeneity were further investigated through micro-CT and SEM image analysis for 1.5-mm diameter miniplug samples from different laminas of the 1-in. core samples. Imbibition depth and rate were calculated on the basis of image analysis as well. Estimation of oil recovery through water imbibition in shale matrix was performed for an example shale field. Here this innovative and integrated multiscale imaging technique provides a “white/gray-box” method to understand water imbibition and water-oil displacement in shale. The wide span of the length scale (from centimeter to nanometer) of this technique enables a more comprehensive, accurate, and specific understanding of both the core-scale dynamics and pore-scale mechanisms of water imbibition, oil recovery, and matrix-fracture interaction.

04 OIL SHALES AND TAR SANDS↗

A compact x-ray diffraction system for dynamic compression experiments on pulsed-power generators

Pulsed-power generators can produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies using the magnetic loading technique. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength–strain rate dependence, structural phase transitions, and density of crystal defects, such as dislocations. Here, we present a cost-effective, compact, pulsed x-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically driven ramp compression of materials with a single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials’ community to investigate in situ dynamic phase transitions critical to equation of states. Finally, we present results using this new diagnostic to evaluate lattice compression in Zr and Al and to capture signatures of phase transitions in CdS.

47 OTHER INSTRUMENTATION↗

A Platform-Independent X-ray Diffraction Diagnostic for Phase Transition Kinetics in Traditional and Synthetic Microstructure Materials (LDRD Project 213088 Final Report)

Pulsed-power generators using the magnetic loading technique are able to produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength-strain rate dependence, structural phase transitions, and density of crystal defects such as dislocations. Here, we present a cost effective, compact X-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically-driven ramp compression of materials with single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials community. The success in fielding this new XRD diagnostic dramatically improves our predictive capability and understanding of rate-dependent behavior at or near phase transition. As Sandia plans the next-generation pulse-power driver platform, a key element needed to deliver new state-of-the-art experiments will be having the necessary diagnostic tools to probe new regimes and phenomena. These diagnostics need to be as versatile, compact, and portable as they are powerful. The development of a platform-independent XRD diagnostic gives Sandia researchers a new window to study the microstructure and phase dynamics of materials under load. This project has paved the way for phase transition research in a variety of materials with mission interest.

36 MATERIALS SCIENCE↗

Simulating water dynamics related to pedogenesis across space and time: Implications for four-dimensional digital soil mapping

Digital soil mapping (DSM) relies on machine-learning and geostatistics to represent soil property observations across space. DSM techniques are powerful but often empirical, being limited to the quality and density of point samples. Water dynamics are closely related to soil variability, and the physics that govern water movement are well known. Hydrological properties can hence be simulated by physical models through space and time, unveiling key characteristics about soils. We propose the use of hydrologic models to map soils across the surface (2D), depth (1D), and time (1D)–which provides a 4D approach to digital soil mapping (4DSM). The Distributed Hydrology Soil Vegetation Model (DHSVM) was applied to a watershed currently under pasture. Moisture sensors and wells were installed at different depths in the watershed on summit, sideslope and toeslope positions to validate the model. DHSVM simulations of soil moisture distribution and depth to saturation were performed during the hydrological year (October 2008-September 2009). Clusters of similar pixels based on soil moisture values were determined using Dynamic Time Warping (DTW) to align temporal data and K-means. Clustering was performed both seasonally and for the entire year. Temporal patterns simulated by DHSVM matched measurements given by moisture sensors and wells. Seasonal clusters differed from the annual cluster. Distinct clusters were observed for each season and with depth, showing that spatiotemporal soil variability is lost when statically assessing soils. Spatiotemporal clusters corroborated field observations of fragipan occurrence not explicitly spatially mapped by Soil Survey Geographic Database (SSURGO). If a connection can be made between water and soils, static and dynamic soil variability can be predicted using physically based hydrologic models. Hydrologic models can benefit soil mapping by enabling reliable 4D simulation of water dynamics, which are fundamental to soil variability and soil classification and directly relate to biological, physical and chemical soil processes not captured by typical soil sampling protocols.

54 ENVIRONMENTAL SCIENCES↗

Sampling Low-Dimensional Markovian Dynamics for Preasymptotically Recovering Reduced Models from Data with Operator Inference

This work introduces a method for learning low-dimensional models from data of high-dimensional black-box dynamical systems. The novelty is that the learned models are exactly the reduced models that are traditionally constructed with classical projection-based model reduction techniques. Thus, the proposed approach learns models that are guaranteed to have the well-studied properties of reduced models known from model reduction, without requiring full knowledge of the governing equations and without requiring the operators of the high-dimensional systems. The key ingredient is a new data sampling scheme to obtain re-projected trajectories of high-dimensional systems that correspond to Markovian dynamics in low-dimensional subspaces. The exact recovery of reduced models from these re-projected trajectories is guaranteed pre-asymptotically under certain conditions for finite amounts of data and for a large class of systems with polynomial nonlinear terms. Numerical results demonstrate that the low-dimensional models learned with the proposed approach match reduced models from traditional model reduction up to numerical errors in practice. In conclusion, the numerical results further indicate that low-dimensional models fitted to re-projected trajectories are predictive even in situations where models fitted to trajectories without re-projection are inaccurate and unstable.

97 MATHEMATICS AND COMPUTING↗

Interacting spin and charge density waves in the kagome metal FeGe

Unveiling the interplay between spin density wave (SDW) and charge density wave (CDW) orders in correlated electron materials is important in obtaining a comprehensive understanding of their electronic, structural, and magnetic properties. Kagome lattice materials are interesting because their flat electronic bands, Dirac points, and Van Hove singularities can enable a variety of exotic electronic and magnetic phenomena. The kagome metal FeGe (the B35 phase), which exhibits a CDW order deep within an A-type antiferromagnetic (AFM) phase, was found to respond dramatically to postgrowth annealing—with the ability to tune the CDW repeatedly from long-range order to negligible order. Additionally, neutron scattering studies suggest that incommensurate magnetic peaks that onset at 𝑇 Canting = 𝑇 SDW ≈ 60 K in the system arise from a SDW order instead of the AFM double-cone structure. Here, in this study, we use inelastic neutron scattering to show that two distinct spin excitations exist below 𝑇 Canting corresponding to two coexisting magnetic orders in the system in both sets of annealed samples with and without CDW. While CDW order or negligible order can dramatically affect the onset temperature of 𝑇 Canting and elastic incommensurate magnetic scattering, its impact on low-energy spin fluctuations is more limited. In both samples, a pair of gapless incommensurate spin excitations arising from the SDW order wave vector coexist with gapped commensurate spin waves from the A-type AFM order across 𝑇 Canting . The low-energy spin excitations for both samples couple dynamically to the lattice through enhanced magnetic scattering intensity on cooling below 𝑇 CDW , regardless of the status of the static long-range CDW order. The incommensurate SDW order in the long-range CDW ordered sample also induces a tiny in-plane lattice distortion of the kagome lattice that is absent in the negligible CDW ordered sample, in a way that is different from the previously known SDW and CDW ordering materials.

charge density waves↗

Dynamic and quasi-static strength of additively repaired aluminum

Additive manufacturing has the potential to repair high value components, saving significant time and resources; however, the level of reliability and performance of additive repairs is still relatively unknown. In this work, the structure–property and performance of laser wire additive manufacturing repairs in 1100 aluminum are investigated. Two types of intentional damage are inflicted on the samples and subsequently repaired with pulsed laser deposition additive manufacturing. Quasi-static (10 –3 s –1 ) and high strain-rate (10 –3 s –1 ) mechanical testing is carried out with in situ diagnostics and post-mortem imaging. The results show that while the quasi-static strength and ductility of samples with a repaired region are lower than a pristine sample, the dynamic strength under shock loading is comparable. This work highlights both the potential utility of additive manufacturing for repair purposes, the significant risk of compromised performance of additive parts under specific conditions, and the need to test at varying strain rates to fully characterize material performance.

36 MATERIALS SCIENCE↗

Integrated transmission electron microscope

An integrated transmission electron microscope comprising multiple electron sources for tuned beams of ultrafast, scanning probe, and parallel illumination in varied beam energies can be alternated within sub-microseconds onto a sample with dynamic ‘transient state’ processes to acquire atomic-scale structural/chemical data with site specificity. The various electron sources and condenser optics enable high-resolution imaging, high-temporal resolution imaging, and chemical imaging, using fast-switching magnets to direct the different electron beams onto a single maneuverable objective pole piece where the sample resides. Such multimodal in situ characterization tools housed in a single microscope have the potential to revolutionize materials science.

Jungjohann, Katherine L.↗

Saddle point search with dynamic active volume

Sampling potential energy surface (PES) is critical for many problems in materials science, chemistry, physics, and biology and requires highly efficient saddle point searches (SPS). In the study, we introduce the concept of dynamic active volume (DAV) in addition to the active volume in self-evolving atomistic kinetic Monte Carlo (SEAKMC). The DAV method has further reduced the dimensionality of the PES at the elevation stage of a SPS. At the subsequent converging stage, the dynamic boundary in DAV is lifted to allow the system to converge to the right location of a saddle point. Coupled with the dimer method, the DAV method not only significantly reduces the time cost for a given search attempt, but also dramatically increases the probability of finding relevant saddle points for PES sampling. Further, a Python software package within the framework SEAKMC (SEAKMC_py) with the DAV method has been developed.

36 MATERIALS SCIENCE↗

Protonation Dynamics of Confined Ethanol–Water Mixtures in H-ZSM-5 from Machine Learning-Driven Metadynamics

Zeolites are indispensable heterogeneous catalysts in industrial chemical processes, valued for their strong Brønsted acidity, well-defined microporous frameworks, and tunable pore structures. Their catalytic activity arises primarily from Brønsted acid sites (BAS), typically present as bridging hydroxyl groups (Si–OH–Al). Under aqueous reaction conditions, these protons interact dynamically with water and alcohol molecules, leading to complex solvation and protonation behavior within confined pores. In this study, we investigate the protonation equilibrium occurring between ethanol and water at the BAS of acidic zeolites under varying hydration levels, i.e., C2H5OH–(H2O)n, n=1–4. Local structure was analyzed through an adaptive-learning global optimization algorithm, while enhanced sampling molecular dynamics simulations with Well-Tempered Metadynamics (WMetaD) and machine learning interatomic potentials (MLPs) provide free-energy surfaces (FES) at variable hydration levels. The results reveal a strong dependence of proton localization on the degree of hydration. At low hydration (1 water molecule), the proton resides predominantly on ethanol; with 2 water molecules, it shifts toward water, and at higher hydration (3 or more water molecules), it becomes extensively delocalized over the water cluster. These findings underscore the critical role of solvation in modulating acid site behavior and suggest that a minimum of three water molecules is necessary to fully stabilize the proton on water within the zeolite framework. This solvation threshold has significant implications for catalytic processes, particularly in biomass conversion reactions where alcohol protonation is a key step in dehydration mechanisms.

machine learning↗

Exceptionally high spallation strength for a high-entropy alloy demonstrated by experiments and simulations

High-entropy alloys are materials with an increasing number of technological applications. Amongst them, the Cantor alloy, FeMnCoCrNi, shows desirable mechanical properties at normal loading conditions. In this study we focus on the performance of the Cantor alloy at the ultra-high deformation rates of shock waves. We study shock-induced spallation using both experiments and atomistic simulations. Experimental loading is achieved using high power laser, with VISAR to obtain velocity profiles and spall strength, followed by transmission electron microscopy of the recovered samples. Molecular Dynamics (MD) simulations of shock-induced spallation are compared with experiments. Both experiments and simulations show a high spall strength which would be beneficial for certain applications, with experiments giving ~8 GPa at ~10 7 s -1 and MD giving almost ~30 GPa at ~10 9 s -1 . The difference between experiments and simulations can be explained by the difference in strain rate. Post-mortem analysis of the experimental samples shows nanotwins near the spall plane, while MD simulations show a highly disordered region giving rise to void nucleation and spall during loading.

36 MATERIALS SCIENCE↗

Adaptive Sampling-Based Bi-Fidelity Stochastic Trust Region Method for Stochastic Derivative-Free Optimization

Bi-fidelity stochastic optimization has gained increasing attention as an efficient approach to reduce computational costs by leveraging a low-fidelity (LF) model to optimize an expensive high-fidelity (HF) objective. In this paper, we propose ASTRO-BFDF, an adaptive sampling trust-region method specifically designed for unconstrained bi-fidelity stochastic derivative-free optimization problems. In ASTRO-BFDF, the LF function serves two purposes: (i) to identify better iterates for the HF function when the optimization process indicates a high correlation between them and (ii) to reduce the variance of the HF function estimates using bi-fidelity Monte Carlo (BFMC). The algorithm dynamically determines sample sizes while adaptively choosing between crude Monte Carlo and BFMC to balance the trade-off between optimization and sampling errors. We prove that the iterates generated by ASTRO-BFDF converge to a first-order stationary point almost surely. Additionally, we demonstrate the effectiveness of the proposed algorithm through numerical experiments on synthetic benchmarks and simulation optimization problems involving discrete event systems.

97 MATHEMATICS AND COMPUTING↗

Scalable molecular dynamics on CPU and GPU architectures with NAMD

NAMD is a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

high-performance computing↗

Solution-phase sample-averaged single-particle spectroscopy of quantum emitters with femtosecond resolution

Here, the development of many optical quantum technologies depends on the availability of solid-state single quantum emitters with near-perfect optical coherence. However, a standing issue that limits systematic improvement is the significant sample heterogeneity and lack of mechanistic understanding of microscopic energy flow at the single emitter level and ultrafast timescales. Here we develop solution-phase single-particle pump-probe spectroscopy with photon correlation detection that captures sample-averaged dynamics in single molecules and/or defect states with unprecedented clarity at femtosecond resolution. We apply this technique to single quantum emitters in two-dimensional hexagonal boron nitride, which suffers from significant heterogeneity and low quantum efficiency. From millisecond to nanosecond timescales, the translation diffusion, metastable-state-related bunching shoulders, rotational dynamics, and antibunching features are disentangled by their distinct photon-correlation timescales, which collectively quantify the normalized two-photon emission quantum yield. Leveraging its femtosecond resolution, spectral selectivity and ultralow noise (two orders of magnitude improvement over solid-state methods), we visualize electron-phonon coupling in the time domain at the single defect level, and discover the acceleration of polaronic formation driven by multi-electron excitation. Corroborated with results from a theoretical polaron model, we show how this translates to sample-averaged photon fidelity characterization of cascaded emission efficiency and optical decoherence time. Our work provides a framework for ultrafast spectroscopy in single emitters, molecules, or defects prone to photoluminescence intermittency and heterogeneity, opening new avenues of extreme-scale characterization and synthetic improvements for quantum information applications.

36 MATERIALS SCIENCE↗