Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 235 records · Page 13

Overview of Dynamic Test Techniques for Flight Dynamics Research at NASA LaRC (Invited)

An overview of dynamic test techniques used at NASA Langley Research Center on scale models to obtain a comprehensive flight dynamics characterization of aerospace vehicles is presented. Dynamic test techniques have been used at Langley Research Center since the 1920s. This paper will provide a partial overview of the current techniques available at Langley Research Center. The paper will discuss the dynamic scaling necessary to address the often hard-to-achieve similitude requirements for these techniques. Dynamic test techniques are categorized as captive, wind tunnel single degree-of-freedom and free-flying, and outside free-flying. The test facilities, technique specifications, data reduction, issues and future work are presented for each technique. The battery of tests conducted using the Blended Wing Body aircraft serves to illustrate how the techniques, when used together, are capable of characterizing the flight dynamics of a vehicle over a large range of critical flight conditions.

Owens, D. Bruce↗

Slosh-Induced Ullage Collapse Validation by Computational Fluid Dynamics and Application to In-space Propellant Dynamics

Ullage collapse predictions made using Computational Fluid Dynamics (CFD) were validated then made for an in-space cryogenic propellant unsettling application. Ullage collapse is the decrease in propellant tank pressure due to cooling of the ullage gas, which can be significant for cryogenic propulsion systems. Vehicles may be developed that are robust to temperature and pressure changes caused by heat leaks, e.g., propellant tanks with high maximum operation pressures, technologies yielding active thermal conditioning, and systems with ample pressurant for conditioning. These methods are robust because additional mass, cost, and complexity are used to address cryogenic fluid management (CFM). As vehicles are optimized for increased payload capacity, tracking of events such as ullage collapse and self-pressurization in tanks will become increasingly important for propellant inventory calculations. Some ullage collapse events can be studied through ground experimentation, but propellant unsettling in and around in-space coast mission phases must be studied in the absence of an acceleration field. Computational Fluid Dynamics is a powerful tool to address this need, enabling prediction of complex physics prior to flight. Validation of the chosen CFD tool and modeling methodology was done using data from a previously conducted subscale slosh-induced ullage collapse ground test performed in part by the Japan Aerospace Exploration Agency (JAXA). A small tank was partially filled with liquid nitrogen and pressurized with gaseous nitrogen at various temperatures. The tank was oscillated laterally to induce high amplitude slosh dynamics resulting in wave breakup. Successful qualitative validation of slosh wave dynamics and quantitative validation of ullage collapse promoted confidence in extending ullage collapse predictions to a flight application. Mechanisms of ullage collapse caused by main engine cutoff while in-space were identified. Both slosh dynamics leading to engine shutdown and body dynamics following shutdown were shown to contribute to propellant unsettling and ullage collapse.

J M Brodnick↗

Dynamic-range studies and improvements for multiplexed photonic Doppler velocimetry

This paper presents studies of the dynamic range achievable with multiplexed photonic Doppler velocimetry (MPDV) measurements, and demonstrates techniques to extend the dynamic range. Improved dynamic range for MPDV measurements is needed in order to track the velocity of the free surface behind a cloud of ejecta. Therefore, the authors have undertaken theoretical and experimental studies of factors affecting dynamic range, particularly in cases where the large number of MPDV probe points precludes high illumination power on each channel. To quantify the potential dynamic range of a given MPDV configuration, they introduce a metric called the frequency-domain number of bits, FNOB, which is less stringent than the formally defined equivalent number of bits (ENOB). This new metric is simple to compute in the lab, and it is well suited to conventional PDV analysis, which does not require digitizer phase coherence beyond tens of nanoseconds.

36 MATERIALS SCIENCE↗

Dynamic response of a freely rotating butterfly valve in the advanced test reactor − dynamic coefficients modeling

Here, in evaluating the water hammer issue pertaining to the primary-coolant-regulating butterfly valve in the Advanced Test Reactor, the dynamic fluid body interaction (DFBI) approach was implemented in the analysis covered in Part I. Although DFBI modeling accurately and simultaneously solved the dynamic motion of the valve’s disk along with the flow field of the surrounding fluid, it shed little light on the reason behind such motion. For Part II, the reacting torque of the fluid on the disk was decomposed into representations of the dynamic coefficients in terms of stiffness, damping, and added mass. These were evaluated via simulations with steady-state static (stiffness), constant angular speed (damping), and variable angular speed (added mass) disks. Substituting the dynamic coefficients into Newton’s second law enabled the response trajectories to be obtained. Stable (by average) and unstable equilibrium positions and thrust tendencies of the valve were determined based on the stiffness coefficient (or static torque), the response amplitude was dampened or enlarged by the damping coefficient (minorly affected by added mass), and the response frequency was altered by the damping and added mass coefficients. Although the dynamic coefficient approach renders slightly different trajectories, due to the averaging effect of the torque in comparison to the DFBI method, the overall trend of the response aligns with the DFBI simulation, thus confirming the conclusion in Part I that a fix to the current butterfly valve is necessary.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Trajectory Ensemble Methods Provide Single-Molecule Statistics for Quantum Dynamical Systems

The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. Here, we observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system–bath entanglement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydration dynamics in water clusters via quantum molecular dynamics simulations

We have investigated the hydration dynamics in size selected water clusters with n = 66, 104, 200, 500, and 1000 water molecules using molecular dynamics simulations. To study the most fundamental aspects of relaxation phenomena in clusters, we choose one of the simplest, still realistic, quantum mechanically treated test solute, an excess electron. The project focuses on the time evolution of the clusters following two processes, electron attachment to neutral equilibrated water clusters and electron detachment from an equilibrated water cluster anion. The relaxation dynamics is significantly different in the two processes, most notably restoring the equilibrium final state is less effective after electron attachment. Nevertheless, in both scenarios only minor cluster size dependence is observed. Significantly different relaxation patterns characterize electron detachment for interior and surface state clusters, interior state clusters relaxing significantly faster. This observation may indicate a potential way to distinguish surface state and interior state water cluster anion isomers experimentally. A comparison of equilibrium and non-equilibrium trajectories suggests that linear response theory breaks down for electron attachment at 200 K, but the results converge to reasonable agreement at higher temperatures. Relaxation following electron detachment clearly belongs to the linear regime. Cluster relaxation was also investigated using two different computational models, one preferring cavity type interior states for the excess electron in bulk water, while the other simulating non-cavity structure. While the cavity model predicts appearance of several different hydrated electron isomers in agreement with experiment, the non-cavity model locates only cluster anions with interior excess electron distribution. The present simulations show that surface isomers computed with the cavity predicting potential show similar dynamical behavior to the interior clusters of the non-cavity type model. Relaxation associated with cavity collapse presents, however, unique dynamical signatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Local vibrational dynamics of hematite (α-Fe{sub 2}O{sub 3}) studied by extended x-ray absorption fine structure and molecular dynamics

The local vibrational dynamics of hematite (α-Fe{sub 2}O{sub 3}) has been investigated by temperature-dependent extended x-ray absorption fine structure spectroscopy and molecular dynamics simulations. The local dynamics of both the short and long nearest-neighbor Fe–O distances has been singled out, i.e., their local thermal expansion and the parallel and perpendicular mean-square relative atomic displacements have been determined, obtaining a partial agreement with molecular dynamics. No evidence of the Morin transition has been observed. More importantly, the strong anisotropy of relative thermal vibrations found for the short Fe–O distance has been related to its negative thermal expansion. The differences between the local dynamics of short and long Fe–O distances are discussed in terms of projection and correlation of atomic motion. As a result, we can conclude that the short Fe–O bond is stiffer to stretching and softer to bending than the long Fe–O bond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Semiclassical Monte Carlo: A first principles approach to non-adiabatic molecular dynamics

Modeling the dynamics of photophysical and (photo)chemical reactions in extended molecular systems is a new frontier for quantum chemistry. Many dynamical phenomena, such as intersystem crossing, non-radiative relaxation, and charge and energy transfer, require a non-adiabatic description which incorporate transitions between electronic states. Additionally, these dynamics are often highly sensitive to quantum coherences and interference effects. Several methods exist to simulate non-adiabatic dynamics; however, they are typically either too expensive to be applied to large molecular systems (10's-100's of atoms), or they are based on ad hoc schemes which may include severe approximations due to inconsistencies in classical and quantum mechanics. We present, in detail, an algorithm based on Monte Carlo sampling of the semiclassical time-dependent wavefunction that involves running simple surface hopping dynamics, followed by a post-processing step which adds little cost. The method requires only a few quantities from quantum chemistry calculations, can systematically be improved, and provides excellent agreement with exact quantum mechanical results. Here we show excellent agreement with exact solutions for scattering results of standard test problems. Additionally, we find that convergence of the wavefunction is controlled by complex valued phase factors, the size of the non-adiabatic coupling region, and the choice of sampling function. These results help in determining the range of applicability of the method, and provide a starting point for further improvement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio two-component Ehrenfest dynamics

We present an ab initio two-component Ehrenfest-based mixed quantum/classical molecular dynamics method to describe the effect of nuclear motion on the electron spin dynamics (and vice versa) in molecular systems. The two-component time-dependent non-collinear density functional theory is used for the propagation of spin-polarized electrons while the nuclei are treated classically. We use a three-time-step algorithm for the numerical integration of the coupled equations of motion, namely, the velocity Verlet for nuclear motion, the nuclear-position-dependent midpoint Fock update, and the modified midpoint and unitary transformation method for electronic propagation. As a test case, the method is applied to the dissociation of H{sub 2} and O{sub 2}. In contrast to conventional Ehrenfest dynamics, this two-component approach provides a first principles description of the dynamics of non-collinear (e.g., spin-frustrated) magnetic materials, as well as the proper description of spin-state crossover, spin-rotation, and spin-flip dynamics by relaxing the constraint on spin configuration. This method also holds potential for applications to spin transport in molecular or even nanoscale magnetic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Combining molecular dynamics with mesoscopic Green’s function reaction dynamics simulations

In many reaction-diffusion processes, ranging from biochemical networks, catalysis, to complex self-assembly, the spatial distribution of the reactants and the stochastic character of their interactions are crucial for the macroscopic behavior. The recently developed mesoscopic Green’s Function Reaction Dynamics (GFRD) method enables efficient simulation at the particle level provided the microscopic dynamics can be integrated out. Yet, many processes exhibit non-trivial microscopic dynamics that can qualitatively change the macroscopic behavior, calling for an atomistic, microscopic description. We propose a novel approach that combines GFRD for simulating the system at the mesoscopic scale where particles are far apart, with a microscopic technique such as Langevin dynamics or Molecular Dynamics (MD), for simulating the system at the microscopic scale where reactants are in close proximity. This scheme defines the regions where the particles are close together and simulated with high microscopic resolution and those where they are far apart and simulated with lower mesoscopic resolution, adaptively on the fly. The new multi-scale scheme, called MD-GFRD, is generic and can be used to efficiently simulate reaction-diffusion systems at the particle level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular dynamics studies of fundamental bulk properties of palladium hydrides for hydrogen storage

Solid-state hydrogen storage materials undergo complex phase transformations whose behavior are collectively determined by thermodynamic (e.g., Gibbs free energy), mechanical (e.g., lattice and elastic constants), and mass transport (e.g., diffusivity) properties. These properties depend on the reaction conditions and evolve continuously during (de)hydrogenation. Thus, they are difficult to measure in experiments. Because of this, past progress to improve solid-state hydrogen storage materials has been prolonged. Using PdH x as a representative example for interstitial metal hydride, we have recently applied molecular dynamics simulations to quantify hydrogen diffusion in the entire reaction space of temperature and composition. In this paper, we have further applied molecular dynamics simulations to obtain well-converged expressions for lattice constants, Gibbs free energies, and elastic constants of PdH x at various stages of the reaction. Finally, our studies confirm significant dependence of elastic constants on temperature and composition. Specifically, a new dynamic effect of hydrogen diffusion on elastic constants is discovered and discussed.

08 HYDROGEN↗

Accurate quantification of lattice temperature dynamics from ultrafast electron diffraction of single-crystal films using dynamical scattering simulations

In ultrafast electron diffraction (UED) experiments, accurate retrieval of time-resolved structural parameters, such as atomic coordinates and thermal displacement parameters, requires an accurate scattering model. Unfortunately, kinematical models are often inaccurate even for relativistic electron probes, especially for dense, oriented single crystals where strong channeling and multiple scattering effects are present. This article introduces and demonstrates dynamical scattering models tailored for quantitative analysis of UED experiments performed on single-crystal films. As a case study, we examine ultrafast laser heating of single-crystal gold films. Comparison of kinematical and dynamical models reveals the strong effects of dynamical scattering within nm-scale films and their dependence on sample topography and probe kinetic energy. Applying to UED experiments on an 11 nm thick film using 750 keV electron probe pulses, the dynamical models provide a tenfold improvement over a comparable kinematical model in matching the measured UED patterns. Also, the retrieved lattice temperature rise is in very good agreement with predictions based on previously measured optical constants of gold, whereas fitting the Debye–Waller factor retrieves values that are more than three times lower. Altogether, these results show the importance of a dynamical scattering theory for quantitative analysis of UED and demonstrate models that can be practically applied to single-crystal materials and heterostructures.

36 MATERIALS SCIENCE↗

Dynamic Load Balancing Based on Constrained K-D Tree Decomposition for Parallel Particle Tracing

Particle tracing is a fundamental technique in flow field data visualization. In this work, we present a novel dynamic load balancing method for parallel particle tracing. Specifically, we employ a constrained k-d tree decomposition approach to dynamically redistribute tasks among processes. Each process is initially assigned a regularly partitioned block along with duplicated ghost layer under the memory limit. During particle tracing, the k-d tree decomposition is dynamically performed by constraining the cutting planes in the overlap range of duplicated data. This ensures that each process is reassigned particles as even as possible, and on the other hand the new assigned particles for a process always locate in its block. Result shows good load balance and high efficiency of our method.

Zhang, Jiang↗

Dynamic Catalysis Fundamentals: I. Fast calculation of limit cycles in dynamic catalysis

Dynamic catalysis—the forced oscillation of catalytic reaction coordinate potential energy surfaces (PES)—has recently emerged as a promising method for the acceleration of heterogeneously-catalyzed reactions. Theoretical study of enhancement of rates and supra-equilibrium product yield via dynamic catalysis has, to-date, been severely limited by onerous computational demands of forward integration of stiff, coupled ordinary differential equations (ODEs) that are necessary to quantitatively describe periodic cycling between PESs. Here, we establish a new approach that reduces, by ≳108×, the computational cost of finding the time-averaged rate at dynamic steady state (i.e. the limit cycle for linear and nonlinear systems of kinetic equations). Our developments are motivated by and conceived from physical and mathematical insight drawn from examination of a simple, didactic case study for which closed-form solutions of rate enhancement are derived in explicit terms of periods of oscillation and elementary step rate constants. Generalization of such closed-form solutions to more complex catalytic systems is achieved by introducing a periodic boundary condition requiring the dynamic steady state solution to have the same periodicity as the kinetic oscillations and solving the corresponding differential equations by linear algebra or Newton-Raphson-based approaches. The methodology is well-suited to extension to non-linear systems for which we detail the potential for multiple solutions or solutions with different periodicities. For linear and non-linear systems alike, the acute decrement in computational expense enables rapid optimization of oscillation waveforms and, consequently, accelerates understanding of the key catalyst properties that enable maximization of reaction rates, conversions, and selectivities during dynamic catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unusual dynamics of tetrahedral liquids caused by the competition between dynamic heterogeneity and structural heterogeneity

Tetrahedral liquids exhibit intriguing thermodynamic and transport properties because of the various ways tetrahedra can be packed and connected. Recently, an unusual temperature dependence of the stretching exponent β in a model tetrahedral liquid ZnCl 2 from T m + 85 K to T m + 35 K has been reported using neutron-spin echo spectroscopy. This discovery stands in sharp contrast to other glass-forming liquids. In this study, we conducted neural network force field driven molecular dynamic simulations of ZnCl 2 . We found a non-monotonic temperature dependence of β from liquid to supercooled liquid temperatures. Further structural decomposition and dynamic analysis suggest that this unusual dynamic behavior is a result of the competition between the decrease in the diversity of tetrahedra motifs (structural heterogeneity) and the increase in glassy dynamic heterogeneity. Furthermore, this result may contribute to new understandings of the structural relaxation of other network liquids.

36 MATERIALS SCIENCE↗

Quantitative Mapping of Nanoscale Chemical Dynamics in Sub-sampled Operando (S)TEM Images Using Spatio-Temporal Analytics

One of the main limitations in the use of operando scanning transmission electron microscopes to study dynamic chemical processes is the effect of the electron beam on the kinetics of the reaction being observed. Here we demonstrate that a flexible Gaussian mixture model can be used to extract quantitative information directly from sub-sampled images, i.e.~images where not all the pixels in the image are illuminated with the beam. The use of this method is demonstrated for the charge/discharge cycle of a Li-battery, where the deposition and dissolution of Li at the anode can be accurately quantified for sub-sampling levels down to 1%. Performing operando imaging using a small fraction of the pixels means that the observations can significantly reduce the effect of the beam, automatically increase the imaging speed and decrease the total data transfer rate required. Such new software capabilities offer the potential to significantly widen the application of operando hardware approaches to study nanoscale dynamics in materials.

Stanfill, Bryan A.↗