Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic simulation”

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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Bond clusters control rupture force limit in shear loaded histidine-Ni 2+ metal-coordinated proteins

Dynamic noncovalent interactions are pivotal to the structure and function of biological proteins and have been used in bioinspired materials for similar roles. Metal-coordination bonds, in particular, are especially tunable and enable control over static and dynamic properties when incorporated into synthetic materials. Despite growing efforts to engineer metal-coordination bonds to produce strong, tough, and self-healing materials, the systematic characterization of the exact contribution of these bonds towards mechanical strength and the effect of geometric arrangements is missing, limiting the full design potential of these bonds. In this work, we engineer the cooperative rupture of metal-coordination bonds to increase the rupture strength of metal-coordinated peptide dimers. Utilizing all-atom steered molecular dynamics simulations on idealized bidentate histidine-Ni 2+ coordinated peptides, we show that histidine-Ni 2+ bonds can rupture cooperatively in groups of two to three bonds. We find that there is a strength limit, where adding additional coordination bonds does not contribute to the additional increase in the protein rupture strength, likely due to the highly heterogeneous rupture behavior exhibited by the coordination bonds. Further, we show that this coordination bond limit is also found natural metal-coordinated biological proteins. Using these insights, we quantitatively suggest how other proteins can be rationally designed with dynamic noncovalent interactions to exhibit cooperative bond breaking behavior. Altogether, this work provides a quantitative analysis of the cooperativity and intrinsic strength limit for metal-coordination bonds with the aim of advancing clear guiding molecular principles for the mechanical design of metal-coordinated materials.

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

Electrostatic Potential of Functional Cations as a Predictor of Hydroxide Diffusion Pathways in Nanoconfined Environments of Anion Exchange Membranes

Nanoconfined anion exchange membranes (AEMs) play a vital role in emerging electrochemical technologies. The ability to control dominant hydroxide diffusion pathways is an important goal in the design of nanoconfined AEMs. Such control can shorten hydroxide transport pathways between electrodes, reduce transport resistance, and enhance device performance. In this work, we propose an electrostatic potential (ESP) approach to explore the effect of the polymer electrolyte cation spacing on hydroxide diffusion pathways from a molecular perspective. By exploring cation ESP energy surfaces and validating outcomes through prior ab initio molecular dynamics simulations of nanoconfined AEMs, we find that we can achieve control over preferred hydroxide diffusion pathways by adjusting the cation spacing. The results presented in this work provide a unique and straightforward approach to predict preferential hydroxide diffusion pathways, enabling efficient design of highly conductive nanoconfined AEM materials for electrochemical technologies.

25 ENERGY STORAGE↗
DOE OSTI2022

Deformation behavior of crystalline/amorphous Al-Si nanocomposites with nanolaminate or nanofibrous microstructures

Deformation mechanisms in sputter-deposited crystalline Al/amorphous Si nanocomposites with nanolaminate or nanofibrous morphology are characterized by nanoindentation, micropillar compression testing and transmission electron microscopy (TEM). The nanofibrous composite having crystalline Al nanofibers with ~40–50 nm in length and 15–20 nm in diameter embedded in amorphous Si exhibits strain hardening to a maximum flow stress of 2.9 GPa and no shear band (SB) formation in compression up to plastic strain exceeding 24%. On the other hand, nanolaminate composite that is composed of 80 nm crystalline Al layers and 20 nm amorphous Si layers exhibits catastrophic SBs starting at plastic strains in the range of 5–10%. Cross-sectional TEM of the deformed samples reveals a high density of stacking faults and twin boundaries in Al nanofibers and no microshear bands in the nanofibrous composite, suggesting plastic deformation in amorphous Si phase and crystalline Al nanofibers. Further, molecular dynamics simulations revealed that the plastic deformation in amorphous Si phase in the cosputtered films could be favored by the decrease in flow strength of amorphous Si with increasing Al solute concentration trapped in Si.

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

A Molecular Description of Hydrogel Forming Polymers for Cement-Based Printing Paste Applications

This research endeavors to link the physical and chemical characteristics of select polymer hydrogels to differences in printability when used as printing aids in cement-based printing pastes. A variety of experimental probes including differential scanning calorimetry (DSC), NMR-diffusion ordered spectroscopy (DOSY), quasi-elastic neutron scattering (QENS) using neutron backscattering spectroscopy, and X-ray powder diffraction (XRD), along with molecular dynamic simulations, were used. Conjectures based on objective measures of printability and physical and chemical-molecular characteristics of the polymer gels are emerging that should help target printing aid selection and design, and mix formulation. Molecular simulations were shown to link higher hydrogen bond probability and larger radius of gyration to higher viscosity gels. Furthermore, the higher viscosity gels also produced higher elastic properties, as measured by neutron backscattering spectroscopy.

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

Evidence for nitrogen binding to surface defects for topological insulator Bi 2 Se 3

Using scanning tunneling spectroscopy and theoretical simulations we have studied the effects of nitrogen gas exposure on the electronic density of states of Bi 2 Se 3 , a well-studied topological insulator. In carefully controlled measurements, Bi 2 Se 3 crystals were initially cleaved in a helium gas environment and then exposed to a 22 SCFH flow of ultra-high purity N 2 gas. We observe a resulting change in the spectral curves, with the exposure effect saturating after approximately 50 min, ultimately bringing the Dirac point about 50 meV closer to the Fermi level. Furthermore, these results are compared to density functional theory calculations, which support a picture of N 2 molecules physisorbing near Se vacancies. Furthermore, ab initio molecular dynamics simulations aided by a Blue Moon ensemble method reveal the dissociative adsorption of N 2 molecules which then bind strongly to Se vacancies at the surface. In this scenario, the binding of the N atom to a Se vacancy site removes the surface defect state created by the vacancy and changes the position of the Fermi energy with respect to the Dirac point.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ultrafast visualization of incipient plasticity in dynamically compressed matter

Plasticity is ubiquitous and plays a critical role in material deformation and damage; it inherently involves the atomistic length scale and picosecond time scale. A fundamental understanding of the elastic-plastic deformation transition, in particular, incipient plasticity, has been a grand challenge in high-pressure and high-strain-rate environments, impeded largely by experimental limitations on spatial and temporal resolution. Here, we report femtosecond MeV electron diffraction measurements visualizing the three-dimensional (3D) response of single-crystal aluminum to the ultrafast laser-induced compression. We capture lattice transitioning from a purely elastic to a plastically relaxed state within 5 ps, after reaching an elastic limit of ~25 GPa. Our results allow the direct determination of dislocation nucleation and transport that constitute the underlying defect kinetics of incipient plasticity. Large-scale molecular dynamics simulations show good agreement with the experiment and provide an atomic-level description of the dislocation-mediated plasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ab-initio simulations of coherent phonon-induced pumping of carriers in zirconium pentatelluride

Abstract Laser-driven coherent phonons can act as modulated strain fields and modify the adiabatic ground state topology of quantum materials. Here we use time-dependent first-principles and effective model calculations to simulate the effect of the coherent phonon induced by strong terahertz electric field on electronic carriers in the topological insulator ZrTe 5 . We show that a coherent A 1g Raman mode modulation can effectively pump carriers across the band gap, even though the phonon energy is about an order of magnitude smaller than the equilibrium band gap. We reveal the microscopic mechanism of this effect which occurs via Landau–Zener-Stückelberg tunneling of Bloch electrons in a narrow region in the Brillouin zone center where the transient energy gap closes when the system switches from strong to weak topological insulator. The quantum dynamics simulation results are in excellent agreement with recent pump-probe experiments in ZrTe 5 at low temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Defects vibrations engineering for enhancing interfacial thermal transport in polymer composites

To push upper boundaries of thermal conductivity in polymer composites, understanding of thermal transport mechanisms is crucial. Despite extensive simulations, systematic experimental investigation on thermal transport in polymer composites is limited. To better understand thermal transport processes, we design polymer composites with perfect fillers (graphite) and defective fillers (graphite oxide), using polyvinyl alcohol (PVA) as a matrix model. Measured thermal conductivities of ~1.38 ± 0.22 W m -1 K -1 in PVA/defective filler composites is higher than those of ~0.86 ± 0.21 W m -1 K -1 in PVA/perfect filler composites, while measured thermal conductivities in defective fillers are lower than those of perfect fillers. We identify how thermal transport occurs across heterogeneous interfaces. Thermal transport measurements, neutron scattering, quantum mechanical modeling, and molecular dynamics simulations reveal that vibrational coupling between PVA and defective fillers at PVA/filler interfaces enhances thermal conductivity, suggesting that defects in polymer composites improve thermal transport by promoting this vibrational coupling.

42 ENGINEERING↗
DOE OSTI2023

Switching of Hybrid Improper Ferroelectricity in Oxide Double Perovskites

In ABO 3 -type perovskite oxides with Pnma symmetry, rotation (Q R+ , a 0 a 0 c + ) and tilt (Q T , a – a – c 0 ) of BO 6 octahedra are the two primary order parameters. These order parameters establish an inherent trilinear coupling with anti-ferroelectric A-site displacement (Q AFE ) to form the low-symmetry phase. The symmetry is further lowered in double perovskite oxides (DPOs) due to A/A' cation ordering. It in turn makes these systems polar via hybrid improper ferroelectric mechanism, primarily driven by Q R+ and Q T . Naturally, it has been believed that functionalities such as polarization can also be switched by tuning these primary order parameters. However, mystery around finding switching mechanism still remains. Our study based on density functional theory calculations combined with finite-temperature molecular dynamics simulations shows that the polarization switching is a two-step process, driven by out-of-phase rotation (Q R– , a 0 a 0 c – when Q T = 0 or, a – a – b – when Q T ≠ 0). A series of polar DPOs such as KLnFeOsO 6 [Ln = Sm, Gd, Dy, Tm (lanthanides) and Y (rare earth)], all belonging to P2 1 symmetry, are considered in this investigation. The polarization switching P ($\overrightarrow{P}$) occurs at a very high temperature of ~1150 K through a phase transition, from a polar (P2 1 ) phase with $\overrightarrow{P}$(+) to $\overrightarrow{P}$(-) via a non-polar P4/n phase. The switching itself is metastable in nature. The switching (both polarization and spin state) is only observed for a very short period of time (~23 ps) that poses limitation on using such a mechanism in memory device realization. We demonstrate a concurrent heating–cooling procedure to overcome such shortcoming. In conclusion, simulations conducted at 600 K further imply that long lasting switching can be achieved, at least for 1.2 ns for 600 K, and ideally for an infinite time, if the material is heated just above the T c followed by rapid cooling to a temperature below T c .

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

Surface Segregation of Liquid Metal Plasma-Facing Component Alloys: A ReaxFF Investigation

Engineering liquid metal alloys offers a transformative pathway for plasma-facing components by enabling chemically tailored surfaces that can simultaneously optimize plasma-material interactions, reduce divertor heat flux, and enhance core plasma confinement, thereby advancing the commercial viability of nuclear fusion power plants. This study, employing an atomistic simulation approach, provides direct evidence that incorporating nonmetal surface-active agents (such as O and H, or their combination) enables strong surface segregation. This capability makes tin−aluminum (Sn−Al) and tin−lithium (Sn−Li) alloys, with suitable compositions, good candidates for PFC applications. Specifically, the presence of low-Z solutes (Li, Al) leads to preferential surface enrichment, which imparts low-Z sputtering characteristics, while the Sn solvent maintains thermophysical stability. To systematically examine this behavior, we developed a ReaxFF force field spanning the full Sn/Al/Li/O/H chemistry, validated it against formation energies and elastic constants, and applied it in reactive molecular dynamics simulations at fusion-relevant temperatures. We also introduced an overlapbased segregation index that captures interfacial compositional separation directly from atomistic density distributions. Here, this metric reveals a clear hierarchy of segregation regimes and provides a unified view across all systems studied. Together, these findings establish a mechanistic link between nonmetal chemistry and interfacial structure, providing a predictive framework for designing self-adaptive, low-sputtering liquid metal alloys for fusion applications.

Alloys↗
DOE OSTI2023

Cation–polymer interactions and local heterogeneity determine the relative order of alkali cation diffusion coefficients in PEGDA hydrogels

Current research efforts are focused on endowing polymer membranes with ion–ion selectivity by incorporating ion–polymer interactions into materials to bias the selective partitioning and or diffusivity of one species over another. However, little is known about the impact of such interactions on the mechanisms of ion transport. In this study, we probe the influence of cation–polymer interactions on cation, anion, and salt diffusivity in a model membrane material, poly(ethylene glycol) diacrylate (PEGDA) by modeling concentrated polyethylene oxide solutions via molecular dynamics simulations. These results are compared to published experimental data for LiCl, NaCl, and KCl diffusion in PEGDA. Experimentally, the order of salt and cation diffusion coefficients for LiCl, NaCl, and KCl deviate from the order in aqueous solutions. Here, simulations identify these deviations to arise from cation–polymer coordination in the membrane. Both the fraction of bound cations and the average binding lifetime increases with decreasing cation hydration free energy (moving down the alkali series), leading to different diffusivity trends in the membrane compared to solution. However, to recover the experimentally observed order of diffusivities cations and salt in our simulations, we needed to incorporate membrane heterogeneity explicitly via a polymer charge scaling procedure. Together, our results indicate that cation–polymer interactions, as well as spatial heterogeneity within the membrane, play a critical role in dictating the observed order of alkali cation and salt diffusion coefficients in membranes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cryo-EM, Protein Engineering, and Simulation Enable the Development of Peptide Therapeutics against Acute Myeloid Leukemia

Cryogenic electron microscopy (cryo-EM) has emerged as a viable structural tool for molecular therapeutics development against human diseases. However, it remains a challenge to determine structures of proteins that are flexible and smaller than 30 kDa. The 11 kDa KIX domain of CREB-binding protein (CBP), a potential therapeutic target for acute myeloid leukemia and other cancers, is a protein which has defied structure-based inhibitor design. Here, we develop an experimental approach to overcome the size limitation by engineering a protein double-shell to sandwich the KIX domain between apoferritin as the inner shell and maltose-binding protein as the outer shell. To assist homogeneous orientations of the target, disulfide bonds are introduced at the target–apoferritin interface, resulting in a cryo-EM structure at 2.6 Å resolution. We used molecular dynamics simulations to design peptides that block the interaction of the KIX domain of CBP with the intrinsically disordered pKID domain of CREB. The double-shell design allows for fluorescence polarization assays confirming the binding between the KIX domain in the double-shell and these interacting peptides. Further cryo-EM analysis reveals a helix–helix interaction between a single KIX helix and the best peptide, providing a possible strategy for developments of next-generation inhibitors.

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

The nano- and meso-scale structure of amorphous calcium carbonate

Abstract Understanding the underlying processes of biomineralization is crucial to a range of disciplines allowing us to quantify the effects of climate change on marine organisms, decipher the details of paleoclimate records and advance the development of biomimetic materials. Many biological minerals form via intermediate amorphous phases, which are hard to characterize due to their transient nature and a lack of long-range order. Here, using Monte Carlo simulations constrained by X-ray and neutron scattering data together with model building, we demonstrate a method for determining the structure of these intermediates with a study of amorphous calcium carbonate (ACC) which is a precursor in the bio-formation of crystalline calcium carbonates. We find that ACC consists of highly ordered anhydrous nano-domains of approx. 2 nm that can be described as nanocrystalline. These nano-domains are held together by an interstitial net-like matrix of water molecules which generate, on the mesoscale, a heterogeneous and gel-like structure of ACC. We probed the structural stability and dynamics of our model on the nanosecond timescale by molecular dynamics simulations. These simulations revealed a gel-like and glassy nature of ACC due to the water molecules and carbonate ions in the interstitial matrix featuring pronounced orientational and translational flexibility. This allows for viscous mobility with diffusion constants four to five orders of magnitude lower than those observed in solutions. Small and ultra-small angle neutron scattering indicates a hierarchically-ordered organization of ACC across length scales that allow us, based on our nano-domain model, to build a comprehensive picture of ACC formation by cluster assembly from solution. This contribution provides a new atomic-scale understanding of ACC and provides a framework for the general exploration of biomineralization and biomimetic processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Probing the electronic properties of the electrified silicon/water interface by combining simulations and experiments.

Silicon (Si) is broadly used in electrochemical and photoelectrochemical devices, where the capacitive and Faradaic reactions at the Si/water interfaces are critical for signal transduction or noise generation. However, probing the electrified Si/water interface at the microscopic level remains a challenging task. Here we focus on hydrogenated Si surfaces in contact with water, relevant to transient electronics and photoelectrochemical modulation of biological cells and tissues. We show that by carrying out first principles molecular dynamics simulations of the Si(100)/water interface in the presence of an electric field we can realistically correlate the computed flat-band potential and tunneling current images at the interface with experimentally measured capacitive and Faradaic currents. Specifically, we validate our simulations in the presence of bias by performing pulsed chronoamperometry measurements on Si wafers in solution. Consistent with prior experiments, our measurements and simulations indicate the presence of voltage-dependent capacitive currents at the interface. We also find that Faradaic currents are weakly dependent on the applied bias, which we relate to surface defects present in newly prepared samples.

Ye, Zifan↗
DOE OSTI2024

The nonexistence of a paddlewheel effect in superionic conductors

Since the 1980s, the paddlewheel effect has been suggested as a mechanism to boost lithium-ion diffusion in inorganic materials via the rotation of rotor-like anion groups. However, it remains unclear whether the paddlewheel effect, defined as large-angle anion group rotations assisting Li hopping, indeed exists; furthermore, the physical mechanism by which the anion-group dynamics affect lithium-ion diffusion has not yet been established. In this work, we differentiate various types of rotational motions of anion groups and develop quaternion-based algorithms to detect, quantify, and relate them to lithium-ion motion in ab initio molecular dynamics simulations. Our analysis demonstrates that, in fact, the paddlewheel effect, where an anion group makes a large angle rotation to assist a lithium-ion hop, does not exist and thus is not responsible for the fast lithium-ion diffusion in superionic conductors, as historically claimed. Instead, we find that materials with topologically isolated anion groups can enhance lithium-ion diffusivity via a more classic nondynamic soft-cradle mechanism, where the anion groups tilt to provide optimal coordination to a lithium ion throughout the hopping process to lower the migration barrier. This anion-group disorder is static in nature, rather than dynamic and can explain most of the experimental observations. Our work substantiates the nonexistence of the long-debated paddlewheel effect and clarifies any correlation that may exist between anion-group rotations and fast ionic diffusion in inorganic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Anisotropic X-Ray Scattering of Transiently Oriented Water

We study the structural dynamics of liquid water by time-resolved anisotropic x-ray scattering under the optical Kerr effect condition. In this way, we can separate the anisotropic scattering decay of 160 fs from the delayed temperature increase of ~ 0.1 K occurring at 1 ps and quantify transient changes in the O-O pair distribution function. Polarizable molecular dynamics simulations reproduce well the experiment, indicating transient alignment of molecules along the electric field, which shortens the nearest-neighbor distances. In addition, analysis of the simulated water local structure provides evidence that two hypothesized fluctuating water configurations exhibit different polarizability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Reinforcement Learning for Intentional Islanding in Resilient Power Transmission Systems

Intentional islanding is the process of identifying and deliberately decomposing the transmission network to form self-sustained islands from an endangered network during disruptions to improve resilience and security. Most existing intentional islanding models are offline resilience decision tools and hence do not provide outage responses in a timely manner. In this paper, a reinforcement learning (RL) based model for intentional islanding is developed, which offers real-time switching control, online deployability, and adaptability to varying system conditions. The intentional islanding process is formulated as a Markov decision process, where the optimal transmission switching policy is learned using the RL approach. The control policy is learned over an environment that encompasses a Power System Simulator for Engineering (PSS/E) model of the transmission network, facilitated by an interface to the standard openAI Gym framework. The proposed RL-based methodology aims to form stable and self-sustainable islands by ensuring voltage stability while reducing the power mismatch in the formed islands. A proximal policy optimization algorithm is designed, which is suitable for controlling the on/off status of the switches with multi-layer perceptron as value and actor networks. The effectiveness of the proposed framework in the self-recovery of the grid by island formation is applied on the modified IEEE 39-bus test network and validated by dynamic simulations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Demonstrate Improved Ag Diffusion and Describe the Basis for Pd Penetration Modeling in SiC

In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗