Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational structural 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Capturing Coupled Structural and Electronic Motions During Excited-State Intramolecular Proton Transfer via Computational Multiedge Resonant Inelastic X-ray Scattering

Proton-transfer processes build the foundation of many chemical processes. In Excited State Intramolecular Proton Transfer (ESIPT) processes the proton transfer process is impulsively started through light. Here, in this study, we explore the changes in coupled atomic and electronic motions during and following ESIPT through computational time-resolved Resonant Inelastic X-ray Scattering (trRIXS). Excited-state Ab Initio Molecular Dynamics (AIMD) simulations combined with Time-Dependent Density Functional Theory (TDDFT) calculations were performed for 10-hydroxybenzo[h]quinoline to obtain trRIXS signatures. RIXS at both the nitrogen and oxygen K-edges were computed to resolve the dynamics electronic structure from both the proton donor and acceptor perspective. The results show how RIXS can reveal the local electronic structure, the coupling between different electronic states and how electronic structure and coupling change during the proton transfer process. Additionally, we observe a strong correlation between spectral changes and structural changes during ESIPT.

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

Fluid-structure interaction solver for transient dynamics of fracturing media

Computer implemented methods include providing a numerical simulation solid mesh comprising a plurality of solid nodes and a numerical simulation fluid mesh comprising a plurality of fluid nodes, and performing a numerical simulation with the solid mesh and fluid mesh including: determining a fluid-solid interaction between the fluid mesh and the solid mesh at a time step by (i) immersing boundaries of the solid and fluid meshes, (ii) computing current corrected fluid and solid velocities using an interaction parameter applied to an intermediate fluid-solid relative velocity, and (iii) computing a fluid-solid interaction force based on the current corrected fluid and solid velocities.

Rougier, Esteban↗
DOE OSTI2021

Full-field imaging learning machine (FILM)

A method of determining dynamic properties of a structure (linear or nonlinear) includes receiving spatio-temporal inputs, generating mode shapes and modal components corresponding to the spatio-temporal inputs using a trained deep complexity coding artificial neural network, and subsequently generating the dynamic properties by analyzing each modal component using a trained learning machine. A computing system for non-contact determination of dynamic properties of a structure includes a camera, a processor, and a memory including computer-executable instructions. When the instructions are executed, the system is caused to receive spatio-temporal image data, decompose the spatio-temporal image data into constituent manifold components using an autoencoder, and analyze the constituent manifold components using a trained learning machine to determine the dynamic properties.

Yang, Yongchao↗
DOE OSTI2025

Phonon spectrum in the spin-Peierls phase of CuGeO 3

CuGeO 3 has long been studied as a prototypical example of the spin-Peierls transition in a 𝑆 = 1/2 Heisenberg chain. Despite intensive investigation of this quasi-one-dimensional material, systematic measurements and calculations of the phonon excitations in the dimerized phase have not to date been possible, leaving certain aspects of the spin-Peierls phenomenon unresolved. We perform state-of-the-art density functional theory (DFT) calculations to compute the electronic structure and phonon dynamics in the low-temperature dimerized phase. We also perform high-resolution neutron spectroscopy to measure the full phonon spectrum over multiple Brillouin zones. We find excellent agreement between our numerical and experimental results that extend to all measurement temperatures. Notable features of our phonon spectra include a number of steeply dispersive modes, nonmonotonic dispersion features, and specific phonon anticrossings, which we relate to the mode eigenvectors. By calculating the magnetic interactions within DFT and studying the effects of different phonon modes on the superexchange paths, we discuss the possibility of observing spin-phonon hybridization effects in experiments performed both in and out of equilibrium.

density functional theory↗
DOE OSTI2024

Aqueous Titania Interfaces

Water–metal oxide interfaces are central to many phenomena and applications, ranging from material corrosion and dissolution to photoelectrochemistry and bioengineering. In particular, the discovery of photocatalytic water splitting on TiO 2 has motivated intensive studies of water-TiO 2 interfaces for decades. So far, a broad understanding of the interaction of water vapor with several TiO 2 surfaces has been obtained. However, much less is known about liquid water–TiO 2 interfaces, which are more relevant to many practical applications. Probing these complex systems at the molecular level is experimentally challenging and is sometimes possible only through computational studies. This review summarizes recent advances in the atomistic understanding, mostly through computational simulations, of the structure and dynamics of interfacial water on TiO 2 surfaces. The main focus is on the nature, molecular or dissociated, of water in direct contact with low-index defect-free crystalline surfaces. The hydroxyls resulting from water dissociation are essential in the photooxidation of water and critically affect the surface chemistry of TiO 2.

density functional theory↗
DOE OSTI2023

4P (Plastic Parallel Pathways Platform) [SWR 23-84]

The Plastic Parallel Pathways Platform (4P) combines life cycle assessment, agent-based modeling within a dynamic material flow analysis structure to compute the environmental impacts of different recycling options under various behavioral interventions.

Walzberg, Julien↗
DOE OSTI2025

Emergent Nanostructure and Ion Transport in Polyzwitterion/Polyanion Blends

We investigated blends of poly(1-(3-sulfonatopropyl)-2-vinylpyridinium) (P2VPPS) and poly(lithium (trifluoromethane)sulfonimide methacrylate) (poly(MTFSI)Li) at varying molar ratios to gain a mechanistic understanding of ionic conductivity in a miscible polyzwitterion/polyanion system. This dataset contains the raw numerical data corresponding to the figures in the manuscript. The data files include the following information: (1) Experimental Data – includes X-ray and neutron scattering measurements, broadband dielectric spectroscopy (BDS) data, extracted DC conductivity values, differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA) results, extracted glass transition temperatures, etc. (2) CGMD Data – includes molecular dynamics (MD) trajectory files and computed structural correlations. All data files are organized/named according to the figure numbers in the manuscript.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Literary review and modeling outlook for radiation damage in fused-silica optical fibers

The purpose of this document is to distill, from literature, the effects of radiation damage in fused-silica optical fibers. This information will be used to design a theoretical framework capable of predictive modeling that will be published elsewhere. By inspection of the relevant literature, the devised method is calculation of defect amorphous structures by molecular dynamic techniques, and computation of the macroscopic dielectric function using density functional theory. The dielectric function is then used to directly compute the transmittance of the core material before and after irradiation via solution of Maxwell’s equations in a homogeneous medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Modeling and Analysis for Spent Nuclear Fuel Seismic Testing

Note: This is the final public release version of PNNL-31671 Draft, which was previously released to the sponsor for review. The intent is the public release version will be PNNL-31671. There are no significant changes from the earlier draft version. This unlimited distribution report is the deliverable for M3SF-21PN010202014: Modeling and Analysis for Spent Nuclear Fuel Seismic Testing. This report summarizes the modeling, analysis, and test plan support completed by Pacific Northwest National Laboratory for the spent nuclear fuel dry storage system seismic test plan through May of 2021. Test plan preparation is planned to continue until the seismic test is completed in July of 2022. This report covers preliminary structural dynamic model development and computer-aided design of test hardware. Based on preliminary modeling, the strongest earthquakes under consideration in this test program provide mechanical loading on the spent nuclear fuel that is comparable to the 30 cm cask drop scenario, although the loads do not appear to be strong enough to cause significant permanent deformation of the fuel assembly spacer grids. The potential for grid deformation during the test will be assessed when the proposed shake table motion becomes available. The weakest earthquakes considered in this test are expected to be comparable to the mechanical loads witnessed in the multimodal transportation test of 2017. The horizontal canister system is predicted to provide nearly-uniform loading condition on the fuel assemblies it contains, while the vertical cask system is predicted to cause non-uniform dynamic loads on the fuel assemblies. In the horizontal case, gravity keeps the fuel assemblies in contact with one basket wall surface unless the loads are strong enough to cause a separation. In the vertical case, the fuel assemblies are relatively long and slender, and seismic motion in the anticipated test range is predicted to cause the assemblies to lean, tilt, and impact the basket walls throughout the seismic event. These gap closures are a nonlinear force transmission condition, so the vertical system is expected to have more variation and variability than the horizontal system. While the horizontal system is expected to have a relatively more linear response than the vertical system, there is still the potential for nonlinear behavior in the horizontal system because the fuel assemblies are free to slide, bounce, and impact the basket walls if the seismic loads are strong enough. One major conclusion of this study is that the use of mixed fuel assemblies in the canister will be acceptable. There are differences in overall system response when the mass, center of gravity, or gaps are changed, but the changes in response are within the bounds of a system that contains completely homogeneous fuel assemblies. One important observation is that the loading for each individual fuel assembly within the vertical canister is expected to be different from that of the others because of nonlinearities in the system. The horizontal canister system is expected to have a more uniform and more predictable response than the vertical canister, making variations in fuel assembly characteristics easier to account for. One important recommendation that comes from this modeling work is to repeat some of the strongest shake tests at different angles, particularly for the vertical cask system. Modeling predicts that strong seismic motion will cause the fuel assemblies to close initial gaps and impact the fuel basket walls when the canister is in the vertical cask configuration. Each shake test will impose a pre-defined three-dimensional time history on the cask system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Using Dissipative Particle Dynamics to Model Effects of Chemical Reactions Occurring within Hydrogels

Computational models that reveal the structural response of polymer gels to changing, dissolved reactive chemical species would provide useful information about dynamically evolving environments. However, it remains challenging to devise one computational approach that can capture all the interconnected chemical events and responsive structural changes involved in this multi-stage, multi-component process. Here, we augment the dissipative particle dynamics (DPD) method to simulate the reaction of a gel with diffusing, dissolved chemicals to form kinetically stable complexes, which in turn cause concentration-dependent deformation of the gel. Using this model, we also examine how the addition of new chemical stimuli and subsequent reactions cause the gel to exhibit additional concentration-dependent structural changes. Through these DPD simulations, we show that the gel forms multiple latent states (not just the “on/off”) that indicate changes in the chemical composition of the fluidic environment. Hence, the gel can actuate a range of motion within the system, not just movements corresponding to the equilibrated swollen or collapsed states. Moreover, the system can be used as a sensor, since the structure of the layer effectively indicates the presence of chemical stimuli.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Atomistic Insights into Lithium–Glyme Solvate Ionic Liquids: Effects of Chain Length and Anion Coordination

For this study, mixtures of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in diglyme (G2), triglyme (G3) and tetraglyme (G4) at solvate ionic liquid (SIL) concentrations were investigated using classical molecular dynamics (cMD) simulations with a physically motivated force-field specifically developed for modeling these systems. The structural and dynamical properties of the mixtures were computed and analyzed. Lithium solvation shells, radial distribution functions, and X-ray structure factors were studied across the different SIL systems. Translational diffusion and rotational relaxation times were also evaluated, exhibiting similar trends with increasing glyme chain length. The results are consistent with experimental data and in good agreement with previous computational studies on G3 and G4. These findings validate the accuracy of the force field in modeling glyme systems and its use for describing the [Li(G2) 4/3 ][TFSI] mixture. Additionally, the thermal and electrochemical stability of these electrolytes were systematically examined. The thermal stability appears to be governed by cooperative interactions among glyme molecules, while the electrochemical stability is primarily influenced by Li + -anion interactions, which vary significantly with glyme chain length. Overall, the study sheds light on the crucial role of the anion in these glyme-based SILs and offers valuable insights into Li + -glyme systems at SIL concentrations, highlighting their promise as potential Li-ion battery electrolytes.

anions↗
DOE OSTI2021

Temperature Homogeneity under Selective and Localized Microwave Heating in Structured Flow Reactors

Selective heating of different phases of multiphase systems via microwaves can result in energy savings and suppression of side reactions. However, materials properties and operating conditions that maximize temperature gradients are poorly understood. Here we utilize computational fluid dynamics (CFD) computations and temperature measurements in structured flow reactors (monoliths) in a monomodal microwave cavity to assess the temperature difference between the walls and the fluid and develop a simple lumped model to estimate when temperature gradients exist. We also explore the material’s thermal and electrical properties of structured reactors for isothermal catalyst conditions. We propose that CFD simulations can be used as a nonintrusive, predictive tool of temperature homogeneity. Importantly, we demonstrate that localized heating in the bed under several conditions rather than selective heating is responsible for the selectivity enhancement. Furthermore, our results indicate that structured beds made of high thermal conductivity materials avoid arcing and enable temperature homogeneity and low electrical conductivity materials allow microwaves to penetrate the domain.

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

Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases

A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

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

What is (quantitative) system dynamics modeling? Defining characteristics and the opportunities they create

A clear definition of system dynamics modeling can provide shared understanding and clarify the impact of the field. We introduce a set of characteristics that define quantitative system dynamics, selected to capture core philosophy, describe theoretical and practical principles, and apply to historical work but be flexible enough to remain relevant as the field progresses. The defining characteristics are: (1) models are based on causal feedback structure, (2) accumulations and delays are foundational, (3) models are equation-based, (4) concept of time is continuous, and (5) analysis focuses on feedback dynamics. We discuss the implications of these principles and use them to identify research opportunities in which the system dynamics field can advance. These research opportunities include causality, disaggregation, data science and AI, and contributing to scientific advancement. Progress in these areas has the potential to improve both the science and practice of system dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Structural evolution in Au- and Pd-based metallic glass forming liquids and the case for improved molecular dynamics force fields

The results of a combined experimental and computational investigation of the structural evolution of Au 81 Si 19 , Pd 82 Si 18 , and Pd 77 Cu 6 Si 17 metallic glass forming liquids are presented. Electrostatically levitated metallic liquids are prepared, and synchrotron x-ray scattering studies are combined with embedded atom method molecular dynamics simulations to probe the distribution of relevant structural units. Metal–metalloid based metallic glass forming systems are an extremely important class of materials with varied glass forming ability and mechanical processibility. High quality experimental x-ray scattering data are in poor agreement with the data from the molecular dynamics simulations, demonstrating the need for improved interatomic potentials. The first peak in the x-ray static structure factor in Pd 77 Cu 6 Si 17 displays evidence for a Curie–Weiss type behavior but also a peak in the effective Curie temperature. In conclusion, a proposed order parameter distinguishing glass forming ability, 1/(S(T,q 1 )−1, shows a peak in the effective Curie temperature near a crossover temperature established by the behavior of the viscosity, T A .

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

Molecular Origins of Temperature-Dependent X-ray Absorption in YAG:Dy

The absorption and emission of X-rays in dysprosium-doped yttrium aluminum garnet (YAG:Dy) has produced unexpected thermographic behavior, which is investigated using a combination of finite temperature ab initio molecular dynamic simulations, structural characterization, and electronic structure calculations of X-ray characteristics. Calculated average peak X-ray absorption spectra (XAS) from simulations between 300 and 600 K result in peak intensity loss due to thermalization effects, matching experimentally measured behavior of YAG:Dy. Here, investigation of atomic snapshots indicates structural factors that correlated with the X-ray behavior, with the first Y–O coordination sphere identified as the primary structural feature unique to high XAS intensity as calculated by radial and pair distribution functions.

Ab initio molecular dynamics↗
DOE OSTI2023

Machine learning assisted derivation of minimal low-energy models for metallic magnets

Abstract We consider the problem of extracting a low-energy spin Hamiltonian from a triangular Kondo Lattice Model (KLM). The non-analytic dependence of the effective spin-spin interactions on the Kondo exchange excludes the use of perturbation theory beyond the second order. We then introduce a Machine Learning (ML) assisted protocol to extract effective two- and four-spin interactions. The resulting spin model reproduces the phase diagram of the original KLM as a function of magnetic field and single-ion anisotropy and reveals the effective four-spin interactions that stabilize the field-induced skyrmion crystal phase. Moreover, this model enables the computation of static and dynamical properties with a much lower numerical cost relative to the original KLM. A comparison of the dynamical spin structure factor in the fully polarized phase computed with both models reveals a good agreement for the magnon dispersion even though this information was not included in the training data set.

Sharma, Vikram (ORCID:0000000189105519)↗