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 181 records · Page 10

Topology optimization of an airfoil fin microchannel heat exchanger using artificial intelligence

High-performance microchannel heat exchangers are needed to supply heat for power conversion for nuclear microreactors. An airfoil fin microchannel design, constructed of Alloy 617 with helium as the working fluid, is analyzed and optimized using a design of experiments with artificial intelligence techniques. The use of airfoil fins offers the potential to reduce pressure drop across the heat exchanger, as compared to other types of channel configurations. A framework for topology optimization of airfoil fin printed circuit heat exchangers (PCHEs) has been developed that can be readily extended to different fin sizes and shapes, as well as different inlet and operating conditions, materials of construction, and working fluids. An optimization procedure is developed that employs computational fluid dynamics for a set of design points identified using Latin hypercube sampling. Computational fluid dynamics is used to analyze a simplified two-channel configuration where five design parameters are varied – inlet angle, fin scale, extent of staggering, transverse and longitudinal pitches. Two methods (a 5D polynomial and a regression neural network) are compared for generating surrogate models and the resulting response surface approximation is input to a genetic algorithm that is used to identify a set of optimal parameters. The optimal geometries are found across six channel Reynolds numbers ranging from 1000 to 5000, since inlet conditions affect flow through the heat exchanger. Additionally, a set of optimal designs that maximizes heat transfer and minimizes pressure drop is identified, and a thermal stress analysis is performed on the optimal design. Correlations for the Nusselt number and Darcy friction factor are developed that can be useful for thermal hydraulic analyses using system codes. Thermal stresses are analyzed and a brief discussion of the status of code cases of PCHEs for nuclear applications is given. Testing and thermomechanical modeling is needed to facilitate future code compliance of PCHEs for high pressure and high temperature applications.

42 ENGINEERING↗

FLAMES─Fast, Low-Storage, Accurate, and Memory-Efficient Adaptive Sampling─Approach to Resolve Spatially Dependent Dynamics of Molecular Liquids

Many critical phenomena in soft matter occur at large length scales, necessitating the resolution of their structure and dynamics at low wavenumbers. However, resolving wavenumber-dependent dynamics computationally via molecular dynamics simulations presents significant challenges, as these phenomena span several orders of magnitude in both time and length scales, resulting in high computational costs and memory demands. Here, this work highlights the computational and memory challenges associated with analyzing molecular trajectories in reciprocal space and demonstrates a method to address them. We introduce FLAMESFast, Low-storage, Accurate, and Memory-Efficient adaptive Sampling, which is a direct method for calculation of structure factors, allowing us to select only the required number of wavevectors for binning. We also use wavenumber-dependent time steps to extract dynamics. Our FLAMES approach effectively mitigates computational and memory/storage bottlenecks. We demonstrate the method using simulations of a model system, liquid octane, at various temperatures. Comparisons with experimental data and real space computation show that the FLAMES technique achieves high accuracy in resolving temperature- and spatially dependent dynamics while being significantly more computationally efficient and requiring less memory and storage than methods based on a uniform wavevector grid and fixed temporal spacing.

Chen, Guang [Argonne National Laboratory (ANL), Ar↗

Insights into the structure and dynamics of K + ions at the muscovite–water interface from machine learning potential simulations

The surfaces of many minerals are covered by naturally occurring cations that become partially hydrated and can be replaced by hydronium or other cations when the surface is exposed to water or an aqueous solution. These ion exchange processes are relevant to various chemical and transport phenomena, yet elucidating their microscopic details is challenging for both experiments and simulations. Here, in this work, we make a first step in this direction by investigating the behavior of the native K + ions at the interface between neat water and the muscovite mica (001) surface with ab-initio -based machine learning molecular dynamics and enhanced sampling simulations. Our results show that the desorption of the surface K + ions in pure ion-free water has a significant free energy barrier irrespective of their local surface arrangement. In contrast, facile K + diffusion between mica’s ditrigonal cavities characterized by different Al/Si orderings is observed. This behavior suggests that the K + ions may favor a dynamic disordered surface arrangement rather than complete desorption when exposed to deionized water.

Ab-initio molecular dynamics↗

Measuring H 0 using X-ray and SZ effect observations of dynamically relaxed galaxy clusters

ABSTRACT We use a sample of 14 massive, dynamically relaxed galaxy clusters to constrain the Hubble constant, H0, by combining X-ray and Sunyaev–Zel’dovich (SZ) effect signals measured with Chandra, Planck, and Bolocam. This is the first such analysis to marginalize over an empirical, data-driven prior on the overall accuracy of X-ray temperature measurements, while our restriction to the most relaxed, massive clusters also minimizes astrophysical systematics. For a cosmological-constant model with Ωm = 0.3 and ΩΛ = 0.7, we find $H_0 = 67.3^{+21.3}_{-13.3}\, \mathrm{km}\, \mathrm{s}^{-1}\, \mathrm{Mpc}^{-1}$, limited by the temperature calibration uncertainty (compared to the statistically limited constraint of $H_0 = 72.3^{+7.6}_{-7.6}\, \mathrm{km}\, \mathrm{s}^{-1}\, \mathrm{Mpc}^{-1}$). The intrinsic scatter in the X-ray/SZ pressure ratio is found to be 13 ± 4 per cent (10 ± 3 per cent when two clusters with significant galactic dust emission are removed from the sample), consistent with being primarily due to triaxiality and projection. We discuss the prospects for reducing the dominant systematic limitation to this analysis, with improved X-ray calibration and/or precise measurements of the relativistic SZ effect providing a plausible route to per cent level constraints on H0.

79 ASTRONOMY AND ASTROPHYSICS↗

Tuning the spin dynamics and magnetic phase transitions of the Cantor alloy via composition and sample processing protocols: A muon spin relaxation study

CrMnFeCoNi, also called the Cantor alloy, is a well-known high-entropy alloy whose magnetic properties have recently become a focus of attention. Here, we present a detailed muon spin relaxation study of the influence of chemical composition and sample processing protocols on the magnetic phase transitions and spin dynamics of several different Cantor alloy samples. Specific samples studied include a pristine equiatomic sample, samples with deficient and excess Mn content, and equiatomic samples magnetized in a field of 9 T or plastically deformed in pressures up to 0.5 GPa. The results confirm the sensitive dependence of the transition temperature on composition and demonstrate that post-synthesis pressure treatments cause the transition to become significantly less homogeneous throughout the sample volume. In addition, we observe critical spin dynamics in the vicinity of the transition in all samples, reminiscent of canonical spin glasses and magnetic materials with ideal continuous phase transitions. Application of an external magnetic field suppresses the critical dynamics in the Mn-deficient sample, while the equiatomic and Mn-rich samples show more robust critical dynamics. The spin-flip thermal activation energy in the paramagnetic phase increases with Mn content, ranging from 3.1⁢(3) × 10 -21 J for 0% Mn to 1.2⁢(2) × 10 -20 J for 30% Mn content. These results shed light on critical magnetic behavior in environments of extreme chemical disorder and demonstrate the tunability of spin dynamics in the Cantor alloy via chemical composition and sample processing.

36 MATERIALS SCIENCE↗

Dynamic Leaching Method for Intact Salt-stone Samples - 18495

At the Savannah River Site (SRS) low activity salt solution is stabilized via encapsulation within a grout termed salt-stone. Salt-stone is emplaced into large (multi-million liter) concrete storage containers referred to as Salt-stone Disposal Units (SDUs). Technetium-99 (Tc-99) and iodine-129 (I-129) are long-lived radionuclides contained in the low activity salt waste and subsequently incorporated into the grout waste form: both radionuclides are contributors to the potential long-term radiation exposure of the public and environment surrounding the SDUs. With respect to Tc-99, in the reducing, high pH environment within the grout, Tc-99 is expected to be relatively immobile since it exists in a reduced Tc(IV) oxidation state in the form of sparingly soluble hydrated oxides (TcO{sub 2}.xH{sub 2}O) or sulfides (TcS{sub x}). However, in the presence of O{sub 2} (associated with the future infiltration of air or oxygenated ground waters into the salt-stone monolith) it is possible for redox-sensitive Tc(IV) to transition into highly soluble (and mobile) Tc(VII) species, such as pertechnetate (TcO{sub 4}{sup -}), which are more readily transported to the surrounding environment. In contrast, I-129 is highly soluble and its leaching behavior is not anticipated to be significantly impacted by potential changes in redox environment or pH within salt-stone. Traditional approaches to quantifying the leaching behavior of Tc-99 and I-29 from cementitious matrices have involved partitioning experiments using size-reduced (crushed/ground) samples, and determination of the radionuclide fraction immobilized by the cementitious solids. Such experiments create artificially high solid-liquid contact areas that likely result in higher Tc-99 and I-129 leachate concentrations than would be expected for intact, monolithic samples. In the current study a new technique, termed the Dynamic Leaching Method (DLM), is being used to initially investigate the Tc-99 leaching behavior from monolithic salt-stone samples; the system will be used at a later date to evaluate the leaching characteristics of I-129. The data derived using this technique is intended to inform the SRS Salt-stone Disposal Facility (SDF) Performance Assessment (PA) which models the long-term transport of radionuclides from the SDUs to the environment. The DLM utilizes a flexible-wall permeameter to achieve saturated leaching under an elevated hydraulic gradient to simulate and accelerate the future transport of groundwater through salt-stone. Leachates are regularly analyzed with respect to radionuclide concentrations, key constituent concentrations (e.g., Al, Ca, Mg, Na, Si, CO{sub 3}{sup 2-}, NO{sub 2}{sup -}, NO{sub 3}{sup -}, SO{sub 4}{sup 2-}), pH, and reduction potential (E{sub h}) in order to develop an understanding of how chemical changes associated with the pore solution influence grout weathering and radionuclide leaching behavior. This paper provides a brief update with respect to data derived from preliminary Tc-99 leaching evaluations as the DLM was being developed, but the paper is predominantly focused on new sample data in which the DLM technique, pre- and post-DLM salt-stone sample characterization, and leachate analysis protocols have been optimized. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ensemble Simulations and Experimental Free Energy Distributions: Evaluation and Characterization of Isoxazole Amides as SMYD3 Inhibitors

Optimization of binding affinities for ligands to their target protein is a primary objective in rational drug discovery. Herein, we report on a collaborative study that evaluates various compounds designed to bind to the SET and MYND domain-containing protein 3 (SMYD3). SMYD3 is a histone methyltransferase and plays an important role in transcriptional regulation in cell proliferation, cell cycle, and human carcinogenesis. Experimental measurements using the scintillation proximity assay show that the distributions of binding free energies from a large number of independent measurements exhibit non-normal properties. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to predict the binding free energies and to provide a detailed chemical insight into the nature of ligand-protein binding. Our results show that the 1-trajectory ESMACS protocol works well for the set of ligands studied here. Although one unexplained outlier exists, we obtain excellent statistical ranking across the set of compounds from the ESMACS protocol and good agreement between calculations and experiments for the relative binding free energies from the TIES protocol. ESMACS and TIES are again found to be powerful protocols for the accurate comparison of the binding free energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The performance of ensemble-based free energy protocols in computing binding affinities to ROS1 kinase

Optimization of binding affinities for compounds to their target protein is a primary objective in drug discovery. Herein we report on a collaborative study that evaluates a set of compounds binding to ROS1 kinase. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to rank the binding free energies. The predicted binding free energies from ESMACS simulations show good correlations with experimental data for subsets of the compounds. Consistent binding free energy differences are generated for TIES and ESMACS. Although an unexplained overestimation exists, we obtain excellent statistical rankings across the set of compounds from the TIES protocol, with a Pearson correlation coefficient of 0.90 between calculated and experimental activities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conformational sampling enhancement of replica exchange molecular dynamics simulations using swarm particle intelligence

In this paper, we present a new method based on swarm particle social intelligence for use in replica exchange molecular dynamics simulations. In this method, the replicas (representing the different system configurations) are allowed communicating with each other through the individual and social knowledge, in additional to considering them as a collection of real particles interacting through the Newtonian forces. The new method is based on the modification of the equations of motion in such way that the replicas are driven towards the global energy minimum. The method was tested for the Lennard-Jones clusters of N = 4, 5, and 6 atoms. Our results showed that the new method is more efficient than the conventional replica exchange method under the same practical conditions. In particular, the new method performed better on optimizing the distribution of the replicas among the thermostats with time and, in addition, ergodic convergence is observed to be faster. We also introduce a weighted histogram analysis method allowing analyzing the data from simulations by combining data from all of the replicas and rigorously removing the inserted bias.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic Event Tree advancements and control logic improvements

The RAVEN code has been under development at the Idaho National Laboratory since 2012. Its main goal is to create a multi-purpose platform for the deploying of all the capabilities needed for Probabilistic Risk Assessment, uncertainty quantification, data mining analysis and optimization studies. RAVEN is currently equipped with three different sampling categories: Forward samplers (Monte Carlo, Latin Hyper Cube, Stratified, Grid Sampler, Factorials, etc.), Adaptive Samplers (Limit Surface search, Adaptive Polynomial Chaos, etc.) and Dynamic Event Tree (DET) samplers (Deterministic and Adaptive Dynamic Event Trees). The main subject of this document is to report the activities that have been done in order to: start the migration of the RAVEN/RELAP-7 control logic system into MOOSE, and develop advanced dynamic sampling capabilities based on the Dynamic Event Tree approach. In order to provide to all MOOSE-based applications a control logic capability, in this Fiscal Year an initial migration activity has been initiated, moving the control logic system, designed for RELAP-7 by the RAVEN team, into the MOOSE framework. In this document, a brief explanation of what has been done is going to be reported. The second and most important subject of this report is about the development of a Dynamic Event Tree (DET) sampler named “Hybrid Dynamic Event Tree” (HDET) and its Adaptive variant “Adaptive Hybrid Dynamic Event Tree” (AHDET). As other authors have already reported, among the different types of uncertainties, it is possible to discern two principle types: aleatory and epistemic uncertainties. The classical Dynamic Event Tree is in charge of treating the first class (aleatory) uncertainties; the dependence of the probabilistic risk assessment and analysis on the epistemic uncertainties are treated by an initial Monte Carlo sampling (MCDET). From each Monte Carlo sample, a DET analysis is run (in total, N trees). The Monte Carlo employs a pre-sampling of the input space characterized by epistemic uncertainties. The consequent Dynamic Event Tree performs the exploration of the aleatory space. In the RAVEN code, a more general approach has been developed, not limiting the exploration of the epistemic space through a Monte Carlo method but using all the forward sampling strategies RAVEN currently employs. The user can combine a Latin Hyper Cube, Grid, Stratified and Monte Carlo sampling in order to explore the epistemic space, without any limitation. From this pre-sampling, the Dynamic Event Tree sampler starts its aleatory space exploration. As reported by the authors, the Dynamic Event Tree is a good fit to develop a goal-oriented sampling strategy. The DET is used to drive a Limit Surface search. The methodology that has been developed by the authors last year, performs a Limit Surface search in the aleatory space only. This report documents how this approach has been extended in order to consider the epistemic space interacting with the Hybrid Dynamic Event Tree methodology.

97 MATHEMATICS AND COMPUTING↗

Gas-Dynamic Research Technique of Indirectly Irradiated Targets in Porous Materials at Iskra-5 Fasility

In this paper, studies of gas-dynamic flows in samples of various materials loaded with an X-ray flux in “Illuminator”-type targets on Iskra-5 laser facility of the ILFI VNIIEF are presented. The obtained results on the velocity of shock waves in the sample under the study are compared with model calculations to match the gas-dynamic parameters to their calculated estimates.

42 ENGINEERING↗

ML-AMD/exa-pd

Exa-pd is a highly parallelizable workflow for constructing multi-element phase diagrams (PDs). It uses standard sampling techniques—molecular dynamics (MD) and Monte Carlo (MC)—as implemented in the LAMMPS package, to simultaneously sample multiple phases on a fine temperature–composition mesh for free-energy calculations. The workflow uses Parsl as a global controller to manage the MD/MC jobs to achieve massive parallelization with almost ideal scalability. The resulting free energies of both liquid and solid phases (including solid solutions) are then fed to CALPHAD modeling using the PYCALPHAD package for the construction of a multi-element PD.

Zhang, Feng [Ames Laboratory (AMES), Ames, IA (Uni↗

Dynamical mean field theory simulations with the adaptive sampling configuration interaction method

In the pursuit of accurate descriptions of strongly correlated quantum many-body systems, dynamical mean field theory (DMFT) has been an invaluable tool for elucidating the spectral properties and quantum phases of both phenomenological models and ab initio descriptions of real materials. Key to the DMFT process is the self-consistent map of the original system into an Anderson impurity model, the ground state of which is computed using an impurity solver. The power of the method is thus limited by the complexity of the impurity model the solver can handle. Simulating realistic systems generally requires many correlated sites. By adapting the recently proposed adaptive sampling configuration interaction (ASCI) method as an impurity solver, we enable much more efficient zero-temperature DMFT simulations. The key feature of the ASCI method is that it selects only the most relevant Hilbert space degrees of freedom to describe the ground state. This reduces the numerical complexity of the calculation, which will allow us to pursue future DMFT simulations with more correlated impurity sites than in previous works. Here, we present the ASCI-DMFT method and example calculations on the one- and two-dimensional Hubbard models that exemplify its efficient convergence and timing properties. We show that the ASCI approach is several orders of magnitude faster than the current best published ground-state DMFT simulations, which allows us to study the bath discretization error in simulations with small clusters, as well as to address cluster sizes beyond the current state of the art. Finally, our approach can also be adapted for other embedding methods such as density matrix embedding theory and self-energy embedding theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nanoscale transient gratings excited and probed by extreme ultraviolet femtosecond pulses

Advances in developing ultrafast coherent sources operating at extreme ultraviolet (EUV) and x-ray wavelengths allow the extension of nonlinear optical techniques to shorter wavelengths. Here, we describe EUV transient grating spectroscopy, in which two crossed femtosecond EUV pulses produce spatially periodic nanoscale excitations in the sample and their dynamics is probed via diffraction of a third time-delayed EUV pulse. The use of radiation with wavelengths down to 13.3 nm allowed us to produce transient gratings with periods as short as 28 nm and observe thermal and coherent phonon dynamics in crystalline silicon and amorphous silicon nitride. This approach allows measurements of thermal transport on the ~10-nm scale, where the two samples show different heat transport regimes, and can be applied to study other phenomena showing nontrivial behaviors at the nanoscale, such as structural relaxations in complex liquids and ultrafast magnetic dynamics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Lipidic cubic phase injector is a viable crystal delivery system for time-resolved serial crystallography

Serial femtosecond crystallography (SFX) using X-ray free-electron laser sources is an emerging method with considerable potential for time-resolved pump-probe experiments. Here we present a lipidic cubic phase SFX structure of the light-driven proton pump bacteriorhodopsin (bR) to 2.3 Å resolution and a method to investigate protein dynamics with modest sample requirement. Time-resolved SFX (TR-SFX) with a pump-probe delay of 1 ms yields difference Fourier maps compatible with the dark to M state transition of bR. Importantly, the method is very sample efficient and reduces sample consumption to about 1 mg per collected time point. Accumulation of M intermediate within the crystal lattice is confirmed by time-resolved visible absorption spectroscopy. Furthermore, this study provides an important step towards characterizing the complete photocycle dynamics of retinal proteins and demonstrates the feasibility of a sample efficient viscous medium jet for TR-SFX.

59 BASIC BIOLOGICAL SCIENCES↗

Reagent-Triggered Isomerization of Fluxional Cluster Catalyst via Dynamic Coupling

Metallic cluster catalysts have many thermodynamically accessible isomers with diverse active sites and low reaction barriers, and lately a strong hypothesis emerged that the many catalyst states collectively drive the catalysis. Furthermore, it remained a hypothesis that catalyst isomerization is actually kinetically feasible under the current reaction conditions. Using high-temperature dynamics simulations and sampling, a range of orientations, and vibrational energy distributions, we probe how thermal effects and molecular events affect cluster catalyst dynamics. We show that even such a delicate affair as the dissociation or scattering of a methane molecule on the heavy and thus slow Pt 13 cluster triggers substantial isomerization of the catalyst, far beyond thermal at 700 K. A kinetic coupling between the methane activity and cluster catalyst dynamics is observed. In return, the thermal dynamics of the cluster affects the methane reaction and scattering probabilities. Hence, molecular events at the surfaces of fluxional cluster catalysts should facilitate the population of an ensemble of catalyst states under the current reaction conditions, with implications for available active sites, reaction mechanisms, and apparent rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗