Search NASA⌕ Search

SEARCH · Search NASA

Results for “DYNAMICS”

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 145 records · Page 8

Application of a New Ensemble Conserving Quantum Dynamics Simulation Algorithm to Liquid para-Hydrogen and ortho-Deuterium

Here, we apply the Feynman-Kleinert Quasi-Classical Wigner (FK-QCW) method developed in our previous work [Smith et al., J. Chem. Phys. 142, 244112 (2015)] for the determination of the dynamic structure factor of liquid para-hydrogen and ortho-deuterium at state points of (T = 20.0 K, n = 21.24 nm -3 ) and (T = 23.0 K, n = 24.61 nm -3 ), respectively. When applied to this challenging system, it is shown that this new FK-QCW method consistently reproduces the experimental dynamic structure factor reported by Smith et al. [J. Chem. Phys. 140, 034501 (2014)] for all momentum transfers considered. Moreover, this shows that FK-QCW provides a substantial improvement over the Feynman-Kleinert linearized path-integral method, in which purely classical dynamics are used. Furthermore, for small momentum transfers, it is shown that FK-QCW provides nearly the same results as ring-polymer molecular dynamics (RPMD), thus suggesting that FK-QCW provides a potentially more appealing algorithm than RPMD since it is not formally limited to correlation functions involving linear operators.

36 MATERIALS SCIENCE↗

Quantum dynamics of hydrogen atoms on graphene. I. System-bath modeling

An accurate system-bath model to investigate the quantum dynamics of hydrogen atoms chemisorbed on graphene is presented. The system comprises a hydrogen atom and the carbon atom from graphene that forms the covalent bond, and it is described by a previously developed 4D potential energy surface based on density functional theory ab initio data. The bath describes the rest of the carbon lattice and is obtained from an empirical force field through inversion of a classical equilibrium correlation function describing the hydrogen motion. By construction, model building easily accommodates improvements coming from the use of higher level electronic structure theory for the system. Further, it is well suited to a determination of the system-environment coupling by means of ab initio molecular dynamics. This paper details the system-bath modeling and shows its application to the quantum dynamics of vibrational relaxation of a chemisorbed hydrogen atom, which is here investigated at T = 0 K with the help of the multi-configuration time-dependent Hartree method. Paper II deals with the sticking dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio molecular dynamics with nuclear quantum effects at classical cost: Ring polymer contraction for density functional theory

Path integral molecular dynamics simulations, combined with an ab initio evaluation of interactions using electronic structure theory, incorporate the quantum mechanical nature of both the electrons and nuclei, which are essential to accurately describe systems containing light nuclei. However, path integral simulations have traditionally required a computational cost around two orders of magnitude greater than treating the nuclei classically, making them prohibitively costly for most applications. Here we show that the cost of path integral simulations can be dramatically reduced by extending our ring polymer contraction approach to ab initio molecular dynamics simulations. By using density functional tight binding as a reference system, we show that our ring polymer contraction scheme gives rapid and systematic convergence to the full path integral density functional theory result. We demonstrate the efficiency of this approach in ab initio simulations of liquid water and the reactive protonated and deprotonated water dimer systems. We find that the vast majority of the nuclear quantum effects are accurately captured using contraction to just the ring polymer centroid, which requires the same number of density functional theory calculations as a classical simulation. Combined with a multiple time step scheme using the same reference system, which allows the time step to be increased, this approach is as fast as a typical classical ab initio molecular dynamics simulation and 35× faster than a full path integral calculation, while still exactly including the quantum sampling of nuclei. This development thus offers a route to routinely include nuclear quantum effects in ab initio molecular dynamics simulations at negligible computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toward more realistic projections of soil carbon dynamics by Earth system models

Soil carbon (C) is a critical component of Earth system models (ESMs), and its diverse representations are a major source of the large spread across models in the terrestrial C sink from the third to fifth assessment reports of the Intergovernmental Panel on Climate Change (IPCC). Improving soil C projections is of a high priority for Earth system modeling in the future IPCC and other assessments. To achieve this goal, we suggest that (1) model structures should reflect real-world processes, (2) parameters should be calibrated to match model outputs with observations, and (3) external forcing variables should accurately prescribe the environmental conditions that soils experience. First, most soil C cycle models simulate C input from litter production and C release through decomposition. The latter process has traditionally been represented by first-order decay functions, regulated primarily by temperature, moisture, litter quality, and soil texture. While this formulation well captures macroscopic soil organic C (SOC) dynamics, better understanding is needed of their underlying mechanisms as related to microbial processes, depth-dependent environmental controls, and other processes that strongly affect soil C dynamics. Second, incomplete use of observations in model parameterization is a major cause of bias in soil C projections from ESMs. Optimal parameter calibration with both pool-and flux-based data sets through data assimilation is among the highest priorities for near-term research to reduce biases among ESMs. Third, external variables are represented inconsistently among ESMs, leading to differences in modeled soil C dynamics. Furthermore, we recommend the implementation of traceability analyses to identify how external variables and model parameterizations influence SOC dynamics in different ESMs. Overall, projections of the terrestrial C sink can be substantially improved when reliable data sets are available to select the most representative model structure, constrain parameters, and prescribe forcing fields.

54 ENVIRONMENTAL SCIENCES↗

Dynamic Data Reconciliation and Validation of a Dynamic Model for Solvent-Based CO 2 Capture Using Pilot-Plant Data

This paper is focused on the development and validation of a dynamic model for an MEA-based CO 2 capture unit. Starting with a rigorous steady-state process model, the dynamic model is developed in Aspen Plus Dynamics ® using a modified Murphree-efficiency approach and is validated using dynamic data collected from the National Carbon Capture Center (NCCC) in Wilsonville, Alabama. The dynamic test runs were designed with due consideration of process excitation and process nonlinearities constrained by the available time to conduct the test runs and implement the experimental designs in the plant control system. The experimental data include solvent composition and loading. Since the experimental data were found to violate mass and energy balances, filtering algorithms coupled with dynamic data reconciliation techniques were used. Without adjusting any model parameters, the dynamic model satisfactorily predicted dynamic response of CO 2 capture due to step changes in the solvent flowrate, flue gas flowrate, and steam flowrates. Transient studies show that the process gain and time constants can considerably change depending on the direction of the disturbance or manipulated variables, indicating process nonlinearities. Due to a storage tank between the absorber and stripper, the time constant of the full plant can be considerably longer than the time constants of the absorber or stripper individually. Furthermore, impact of the control system on the process efficiency while following an optimal CO 2 capture schedule was also studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Heterogeneous nature of relaxation dynamics of room-temperature ionic liquids (EMIm) 2 [Co(NCS) 4 ] and (BMIm) 2 [Co(NCS) 4 ]

Dynamic crossover above T g has been recognized as a characteristic feature of molecular dynamics of liquids approaching glass transition. Experimentally, it is manifested as a change in Vogel Fulcher Tammann dependence or a breakdown of the Stokes Einstein and related relations. In this paper, we report the exception from this rather general pattern of behavior. By means of dielectric, ultrasonic, rheological, and calorimetric methods, dynamics of two good ionic conductors (BMIm) 2 [Co(NCS) 4 ] and (EMIm) 2 [Co(NCS) 4 ] of less common stoichiometry (2:1) was studied in a very broad temperature range. However, none of the mentioned dynamic changes was observed in the entire studied temperature range. On the contrary, the single VFT and the same fractional Walden coefficient were found for conductivity and viscosity changes over 12 decades. Furthermore, ultrasonic studies revealed that the data at temperatures which cover the normal liquid region cannot be fitted by a single exponential decay, and the Cole Cole function should be used instead.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Graph dynamical networks for unsupervised learning of atomic scale dynamics in materials

Understanding the dynamical processes that govern the performance of functional materials is essential for the design of next generation materials to tackle global energy and environmental challenges. Many of these processes involve the dynamics of individual atoms or small molecules in condensed phases, e.g. lithium ions in electrolytes, water molecules in membranes, molten atoms at interfaces, etc., which are difficult to understand due to the complexity of local environments. In this work, we develop graph dynamical networks, an unsupervised learning approach for understanding atomic scale dynamics in arbitrary phases and environments from molecular dynamics simulations. We show that important dynamical information, which would be difficult to obtain otherwise, can be learned for various multi-component amorphous material systems. With the large amounts of molecular dynamics data generated every day in nearly every aspect of materials design, this approach provides a broadly applicable, automated tool to understand atomic scale dynamics in material systems.

36 MATERIALS SCIENCE↗

Chaotic dynamics of flexible Euler-Bernoulli beams

Mathematical modeling and analysis of spatio-temporal chaotic dynamics of flexible simple and curved Euler-Bernoulli beams are carried out. The Kármán-type geometric non-linearity is considered. Algorithms reducing partial differential equations which govern the dynamics of studied objects and associated boundary value problems are reduced to the Cauchy problem through both Finite Difference Method with the approximation of O(c{sup 2}) and Finite Element Method. The obtained Cauchy problem is solved via the fourth and sixth-order Runge-Kutta methods. Validity and reliability of the results are rigorously discussed. Analysis of the chaotic dynamics of flexible Euler-Bernoulli beams for a series of boundary conditions is carried out with the help of the qualitative theory of differential equations. We analyze time histories, phase and modal portraits, autocorrelation functions, the Poincaré and pseudo-Poincaré maps, signs of the first four Lyapunov exponents, as well as the compression factor of the phase volume of an attractor. A novel scenario of transition from periodicity to chaos is obtained, and a transition from chaos to hyper-chaos is illustrated. In particular, we study and explain the phenomenon of transition from symmetric to asymmetric vibrations. Vibration-type charts are given regarding two control parameters: amplitude q{sub 0} and frequency ω{sub p} of the uniformly distributed periodic excitation. Furthermore, we detected and illustrated how the so called temporal-space chaos is developed following the transition from regular to chaotic system dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Vibrational dynamics of aqueous hydroxide solutions probed using broadband 2DIR spectroscopy

We employed ultrafast transient absorption and broadband 2DIR spectroscopy to study the vibrational dynamics of aqueous hydroxide solutions by exciting the O–H stretch vibrations of the strongly hydrogen-bonded hydroxide solvation shell water and probing the continuum absorption of the solvated ion between 1500 and 3800 cm{sup −1}. We observe rapid vibrational relaxation processes on 150–250 fs time scales across the entire probed spectral region as well as slower vibrational dynamics on 1–2 ps time scales. Furthermore, the O–H stretch excitation loses its frequency memory in 180 fs, and vibrational energy exchange between bulk-like water vibrations and hydroxide-associated water vibrations occurs in ∼200 fs. The fast dynamics in this system originate in strong nonlinear coupling between intra- and intermolecular vibrations and are explained in terms of non-adiabatic vibrational relaxation. These measurements indicate that the vibrational dynamics of the aqueous hydroxide complex are faster than the time scales reported for long-range transport of protons in aqueous hydroxide solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Photonic Doppler velocimetry probe used to measure grain boundaries of dynamic shocked materials

Material scientists have developed computational modeling to predict the dynamic response of materials undergoing stress, but there is still a need to make precision measurements of surfaces undergoing shock compression. Miniature photonic Doppler velocimetry (PDV) probes have been developed to measure the velocity distribution from a moving surface traveling tens of kilometers per second. These probes use hundreds of optical fibers imaged by optical relays onto different regions of this moving surface. While previous work examined large surface areas, we have now developed a PDV microscope that can interrogate 37 different spots within a field of view of <1 mm, with a standoff distance of 17 mm, to analyze the motion differences across grain boundaries of the material undergoing dynamic stress. Each PDV fiber interrogates a 10 μm spot size on the moving surface. A separate imaging system using a coherent bundle records the location of the PDV spots relative to the grain boundaries prior to the dynamic event. Designing the mounting structures for the lenses, fibers, and coherent bundle was a challenge. To minimize back reflections, the fibers are index matched onto the first relay lens, which is made of fused silica. The PDV fibers are aligned normal to the moving surface. The imaging probe views the surface at an 18° angle. The coherent bundle is tilted 11° to its optical relay. All components are assembled into a single probe head assembly. The coherent bundle is removed from the probe head to be used for the next dynamic event. Alignment issues will also be discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A computational fluid dynamics model to estimate local quantities in firebrand char oxidation

Firebrand burning is a complex phenomenon that is influenced by several parameters which are difficult to fully explore experimentally. Computational fluid dynamics models capable of predicting local quantities are essential for accurate prediction of char oxidation in firebrands. This article presents a computational fluid dynamics model to estimate firebrand mass loss, diameter change, and surface temperature during char oxidation. The model was validated using previously conducted wind tunnel experiments. These experiments were conducted for firebrands of two different aspect ratios, which were arranged in three different configurations (single, horizontal array, and vertical array), and for four different wind speeds (0.5, 1, 1.5, and 2 m/s). The computational fluid dynamics results were compared with a previous 1 D model. In all the test cases, the computational fluid dynamics model predicted the physical phenomena with significantly improved accuracy compared to a 1 D model. The char oxidation model presented in this article can be coupled with other models to study firebrand generation and trajectory, biomass pyrolysis, fluidized bed reactors, and coal combustion.

Engineering↗

Important powertrain dynamics for developing models for control of connected and automated electrified vehicles

Connected and Automated Vehicles (CAV) technology presents significant opportunities for energy saving in the transportation sector. CAV technology forecasts vehicle and powertrain power needs under various terrain, ambient, and traffic conditions. Integration of the CAV technology in Hybrid Electric Vehicles (HEVs) provides the opportunity for optimal vehicle operation. Indeed, Hybrid Electric Vehicle powertrains present high degrees of flexibility and possibility for choosing optimum powertrain modes based on the predicted traction power needs. In modeling complex CAV powertrain dynamics, the modeler needs to consider short-time scale powertrain dynamics, such as engine transients, and hysteresis of mode-switching for a multi-mode HEV. Therefore, the powertrain dynamics essential for developing powertrain controllers for a class of connected HEVs is presented. To this end, control-oriented powertrain dynamic models for a test vehicle consisting of full electric, hybrid, and conventional engine operating modes are developed. The resulting powertrain model can forecast vehicle traction torque and energy consumption for the specified prediction horizon of the test vehicle. The model considers different operating modes and associated energy penalty terms for mode switching. Thus, the vehicle controller can determine the optimum powertrain mode, torque, and speed for forecasted vehicle operation via utilizing connectivity data. The powertrain model is validated against the experimental data and shows prediction error of less than 5% for predicting vehicle energy consumption. The model is used to create energy penalty maps that can be used for CAV control, for example fuel penalty map for engine torque changes (10–40 Nm) at each engine speed. The results of model-based optimization show optimum switching delays ranging from 0.4 to 1.4 s to avoid hysteresis in mode switching.

Engineering↗

On the characteristics of aerosol indirect effect based on dynamic regimes in global climate models

Aerosol-cloud interactions continue to constitute a major source of uncertainty for the estimate of climate radiative forcing. The variation of aerosol indirect effects (AIE) in climate models is investigated across different dynamical regimes, determined by monthly mean 500 hPa vertical pressure velocity (ω500), lower-tropospheric stability (LTS) and large-scale surface precipitation rate derived from several global climate models (GCMs), with a focus on liquid water path (LWP) response to cloud condensation nuclei (CCN) concentrations. The LWP sensitivity to aerosol perturbation within dynamic regimes is found to exhibit a large spread among these GCMs. It is in regimes of strong large-scale ascend (ω500 < -25 hPa/d) and low clouds (stratocumulus and trade wind cumulus) where the models differ most. Shortwave aerosol indirect forcing is also found to differ significantly among different regimes. Shortwave aerosol indirect forcing in ascending regimes is as large as that in stratocumulus regimes, which indicates that regimes with strong large-scale ascend are as important as stratocumulus regimes in studying AIE. 42" It is further shown that shortwave aerosol indirect forcing over regions with high monthly large-scale surface precipitation rate (> 0.1 mm/d) contributes the most to the total aerosol indirect forcing (from 64% to nearly 100%). Results show that the uncertainty in AIE is even larger within specific dynamical regimes than that globally, pointing to the need to reduce the uncertainty in AIE in different dynamical regimes.

Zhang, S.↗

Mapping the conformational landscape of a dynamic enzyme by multitemperature and XFEL crystallography

Determining the interconverting conformations of dynamic proteins in atomic detail is a major challenge for structural biology. Conformational heterogeneity in the active site of the dynamic enzyme cyclophilin A (CypA) has been previously linked to its catalytic function, but the extent to which the different conformations of these residues are correlated is unclear. Here we compare the conformational ensembles of CypA by multitemperature synchrotron crystallography and fixed-target X-ray free-electron laser (XFEL) crystallography. The diffraction-before-destruction nature of XFEL experiments provides a radiation-damage-free view of the functionally important alternative conformations of CypA, confirming earlier synchrotron-based results. We monitored the temperature dependences of these alternative conformations with eight synchrotron datasets spanning 100-310 K. Multiconformer models show that many alternative conformations in CypA are populated only at 240 K and above, yet others remain populated or become populated at 180 K and below. These results point to a complex evolution of conformational heterogeneity between 180-–240 K that involves both thermal deactivation and solvent-driven arrest of protein motions in the crystal. The lack of a single shared conformational response to temperature within the dynamic active-site network provides evidence for a conformation shuffling model, in which exchange between rotamer states of a large aromatic ring in the middle of the network shifts the conformational ensemble for the other residues in the network. Together, our multitemperature analyses and XFEL data motivate a new generation of temperature- and time-resolved experiments to structurally characterize the dynamic underpinnings of protein function.

Keedy, Daniel A. (ORCID:0000000291847586)↗

Rotational Dynamics and Transition Mechanisms of Surface-Adsorbed Proteins

Assembly of biomolecules at solid-water interfaces requires molecules to traverse complex orientation-dependent energy landscapes through processes that are poorly understood, largely due to the dearth of in-situ single molecule measurements and statistical analyses of the rotational dynamics that define directional selection. Emerging capabilities in high-speed atomic force microscopy and machine learning have allowed us to directly determine the orientational energy landscape and observe and quantify the rotational dynamics for protein nanorods on the surface of muscovite mica under a variety of conditions. Comparisons with kinetic Monte Carlo simulations show that the transition rates between adjacent orientation-specific energetic minima can largely be understood through traditional models of in-plane Brownian rotation across a biased energy landscape, with resulting transition rates that are exponential in the energy-barriers between states. However, transitions between more distant angular states are decoupled from barrier height, with jump-size distributions showing a power-law decay that is characteristic of a non-classical Levy-flight random walk, indicating that large jumps are enabled by alternative modes of motion via activated states. The findings provide new insights into the dynamics of biomolecules at solid-liquid interfaces that lead to self-assembly, epitaxial matching and other orientationally anisotropic outcomes and define a general procedure for exploring such dynamics with implications for hybrid biomolecular-inorganic materials design.

Orientational energy landscapes, Rotational dynami↗

Abnormal Slow Phonon Dynamics Toward Prolonging Excited States Dynamics Enabled by Crystalline‐Assembling Donor–Acceptor Molecules

Phonon dynamics are a critical factor to control the optical properties of excited states in light-emitting materials. Here, we report an extremely slow relaxation of photoexcited lattice vibrations enabled by assembling the donor-acceptor (D–A) molecules [2-(9,9-dimethylacridin-10(9H)-yl)-9,9-dimethyl-9H-thioxanthene 10,10-dioxide], namely AC molecules, into dipolar crystal. By using photoexcitation-modulated Raman spectroscopy, we find that the crystalline-lattice vibrations monitored by Raman-scattering laser beam of 785 nm demonstrate an un-usual slow relaxation in the time scale of seconds after ceasing photoexcitation beam of 343 nm in such dipolar crystal. This presents extremely slow phonon dynamics enabled by crystalline-assembling the D–A molecules into a dipolar crystal. Simultaneously, the photoluminescence (PL) exhibits a prolonged behavior, lasting 10 ms after ceasing photoexcitation in dipolar AC crystal. This phenomenon provides an experimental hypothesis that the slow phonon dynamics function as an important mechanism to unusually prolong excited states dynamics upon crystalline-assembling the D–A molecules into dipolar crystal. Here, this hypothesis can be verified by directly suppressing the phonon dynamics through freezing D–A molecular liquid into dipolar crystalline solid at 77 K to largely prolong the PL to 1 s- after removing photoexcitation. Clearly, crystalline-assembling D–A molecules provide the necessary conditions to enable slow phonon dynamics toward prolonging excited states dynamics.

36 MATERIALS SCIENCE↗

ReaxFF molecular dynamics simulations on the structure and dynamics of electrolyte water systems at ambient temperature

Here, ReaxFF molecular dynamics simulations have been performed to study the effect of cations Li + , Na+ and K + and anion Cl - on the structural and dynamical properties of water, using the force field recently developed by Fedkin and co-workers. The structural relationship of ion and water has been analyzed from the radial distribution function and angular distribution. Comparisons of ReaxFF angle variation of ions and water within the first solvation shell were made and found to be in good agreement with literature. The disruption of hydrogen bond network of water by ions is elucidated by ion-water residence times, water-water hydrogen bond dynamics and reorientational dynamics. ReaxFF diffusion coefficient and residence times of electrolyte water system were compared with ab initio and non-reactive potentials to analyze the difference in dynamics. We gained insight into the ion interaction with water and how it can accelerate or decelerate water dynamics. ReaxFF outlines the formation and dissolution of metal hydroxides and metal chlorides over the course of simulation to explain the diffusion dynamics of water in salt solutions, allowing us to elucidate the impact of concentration on the self-diffusivity of water and ions in solutions, and to reveal that this effect always decreases the mobility and is not at all ion-specific. The obtained findings have opened new opportunities to extend the ReaxFF methodology towards systems involving electrolytes.

42 ENGINEERING↗

Non-adiabatic molecular dynamics by accelerated semiclassical Monte Carlo

Non-adiabatic dynamics, where systems non-radiatively transition between electronic states, plays a crucial role in many photo-physical processes, such as fluorescence, phosphorescence, and photoisomerization. Methods for the simulation of non-adiabatic dynamics are typically either numerically impractical, highly complex, or based on approximations which can result in failure for even simple systems. Recently, the Semiclassical Monte Carlo (SCMC) approach was developed in an attempt to combine the accuracy of rigorous semiclassical methods with the efficiency and simplicity of widely used surface hopping methods. However, while SCMC was found to be more efficient than other semiclassical methods, it is not yet as efficient as is needed to be used for large molecular systems. Here, we have developed two new methods: the accelerated-SCMC and the accelerated-SCMC with re-Gaussianization, which reduce the cost of the SCMC algorithm up to two orders of magnitude for certain systems. In most cases shown here, the new procedures are nearly as efficient as the commonly used surface hopping schemes, with little to no loss of accuracy. This implies that these modified SCMC algorithms will be of practical numerical solutions for simulating non-adiabatic dynamics in realistic molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗