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 OSTI2021

Effect of Helium Flux on Near-surface Helium Accumulation in Plasma-Exposed Tungsten

We report results of object kinetic Monte Carlo (OKMC) simulations aimed at understanding the effect of helium flux on the near-surface helium accumulation in plasma-facing tungsten, which is initially defect-free and has a W(100) surface orientation. These OKMC simulations are performed at 933 K for fluxes ranging from 10 22 to 4 × 10 25 He/m 2 s, with 100 eV helium atoms impinging on a W(100) surface up to a maximum fluence of 4×10 19 He/m 2 . In the near-surface region, helium clusters interact elastically with the free surface. The interaction is attractive and results in the drift of mobile helium clusters towards the surface as well as increased trap mutation rates. The associated kinetics and energetics of the above-mentioned processes obtained from molecular dynamics simulations are also considered. The OKMC simulations indicate that as the flux decreases, the retention of implanted helium decreases, and its depth distribution shifts to deeper below the surface in initially pristine tungsten. Furthermore, the fraction of retained helium diffusing into the bulk increases as well, so much so that for 10 22 He/m 2 s, almost all of the retained helium diffused into the bulk with minimal/negligible near-surface helium accumulation. At a given flux, with increasing fluence, the fraction of retained helium initially decreases and then starts to increase after reaching a minimum. The occurrence of the retention minimum shifts to higher fluences as the flux decreases. Although the near-surface helium accumulation spreads deeper into the material with decreasing flux and increasing fluence, the spread appears to saturate at depths between 80 and 100 nm. Finally, we present a detailed analysis of the influence of helium flux on the size and depth distribution of total helium and helium bubbles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Asymmetric Composition of Ionic Aggregates and the Origin of High Correlated Transference Number in Water-in-Salt Electrolytes

"Water-in-salt" electrolytes open up exciting new avenues for expanding the electrochemical window of aqueous electrolytes. We investigated the solvation structure and dynamics of highly concentrated lithium bis(trifluoromethane)sulfonimide aqueous electrolyte using experimentally corroborated molecular dynamics simulations. The simulations revealed that the heterogeneous structure of the electrolyte comprises percolating networks of ion and water domains/aggregates. Interestingly, the ionic regions are composed of more TFSI - ions than Li + ions. The Li + -ion transport mechanism was further explored. Li + ions can hop along the coordinated TFSI- ions in the ionic aggregates. The calculated correlated transference number of the 20 m electrolyte is similar to 0.32, which is reasonably high for the high concentration due to a weak negative correlation between the motion of cations and anions within the heterogeneous microscopic domains. These molecular dynamics results connect the heterogeneous structure of the electrolyte to the correlated dynamics of the Li + ion and provide a new understanding of the Li + -ion transport mechanism in this novel electrolyte.

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

Cold-cap structure in a slurry-fed electric melter

We report as glass batch is charged into an electric melter, a cold cap forms on the glass melt surface. Heat transfer to the cold cap from the molten glass below and the melter atmosphere above determines the melting rate. A mathematical model of the cold cap and the experimental kinetic data of the feed-to-glass conversion that were collected for several simulated low-activity and high-level waste melter feeds allowed us to develop relationships between the internal structure of the cold cap, its properties, its thickness, and the internal heat transfer. This contribution shows the distribution of major crystalline phases and the cumulative evolution of gases within the cold cap. It also examines the temperature, conversion degree, and heating rate the melter feed is experiencing during the passage through the cold cap and their effects on the cold-cap bottom temperature and morphology, which are important for the computational fluid dynamics simulations of melters.

36 MATERIALS SCIENCE↗
DOE OSTI2020

MoorDyn V2: New Capabilities in Mooring System Components and Load Cases

MoorDyn, an open-source mooring dynamics model, is being expanded with capabilities for additional mooring system features and load cases. As floating wind turbine technology matures, mooring systems are becoming more sophisticated and more complex scenarios need to be considered in the design process. Mooring systems may have synthetic line materials, ballast/buoyancy bodies along the lines, or interconnections between platforms. Failure modes may involve multiple cascading line failures that depend on mooring system dynamics. Features recently added to MoorDyn aim to address these emerging needs. MoorDyn’s linear elasticity model has been supplemented to support user-defined stress-strain curves, which can be adjusted to represent synthetic mooring materials. Rigid six-degree-of-freedom bodies in the mooring system can now be modeled using two new model objects. “Rod” objects provide an option for rigid cylindrical bodies. They use the existing Morison equation-based hydrodynamics model and can be connected to mooring lines at either end. “Body” objects provide a generic six-degree-of-freedom rigid-body representation based on a lumped-parameter model of translational and rotational properties. Rod objects can be added to Body objects and mooring lines can be attached at any location, allowing a wide variety of submerged structures to be integrated into the mooring system. Lastly, a means of dynamically simulating mooring line failures has been implemented. These new features, currently in the C++ version of MoorDyn, are described and then demonstrated on a two-turbine shared-mooring array. A qualitative view of the results suggests the new features are functioning as expected.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2020

MoorDyn V2: New Capabilities in Mooring System Components and Load Cases: Preprint

MoorDyn, an open-source mooring dynamics model, is being expanded with capabilities for additional mooring system features and load cases. As floating wind turbine technology matures, mooring systems are becoming more sophisticated and more complex scenarios need to be considered in the design process. Mooring systems may have synthetic line materials, ballast/buoyancy bodies along the lines, or interconnections between platforms. Failure modes may involve multiple cascading line failures that depend on mooring system dynamics. Features recently added to MoorDyn aim to address these emerging needs. MoorDyn’s linear elasticity model has been supplemented to support user-defined stress-strain curves, which can be adjusted to represent synthetic mooring materials. Rigid six-degree-of-freedom bodies in the mooring system can now be modeled using two new model objects. “Rod” objects provide an option for rigid cylindrical bodies. They use the existing Morison equation-based hydrodynamics model and can be connected to mooring lines at either end. “Body” objects provide a generic six-degree-of-freedom rigid-body representation based on a lumped-parameter model of translational and rotational properties. Rod objects can be added to Body objects and mooring lines can be attached at any location, allowing a wide variety of submerged structures to be integrated into the mooring system. Lastly, a means of dynamically simulating mooring line failures has been implemented. These new features, currently in the C++ version of MoorDyn, are described and then demonstrated on a two-turbine shared-mooring array. A qualitative view of the results suggests the new features are functioning as expected.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2022

Dynamic aspects of graphene deformation and fracture from approximate density functional theory

Graphene is one of the most intriguing two-dimensional carbon materials. Its mechanical strength and failure are key concerns for materials engineering and applications. Despite the success of fracture mechanics, the mechanism of how pristine materials fail remains an elusive problem. While many theoretical studies based on molecular dynamics using empirical forcefields have tried to address this question, atomic-scale mechanics are not clearly understood. Especially, a widely employed bond-breaking approach based on the critical bond length has not been rigorously tested. Here, utilizing molecular dynamics simulations with density functional based tight binding, we investigate how the failure of the pristine material initiates. The Wiberg bond order (W BO ) to estimate the change of chemical bonds shows a transition from sp 2 (W BO ~ 1.33) to sp 3 (WBO < 1.0) during the deformation. However, it reveals that a single threshold value for either the WBO or bond length is insufficient to decide failure of pristine material without free edges or defects. Instead, collective behaviors of the local atomic group govern the fracture initiation of pristine graphene. Overall, our study provides dynamic mechanical responses based on quantum mechanics, which have not been captured by empirical forcefields, opening opportunities to design properties by precisely coupling the mechanics and quantum mechanics.

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

Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth

Abstract. Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT) as well as in the PBL at low relative humidity and low temperature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Dynamics of collision-induced energy transfer

Collisions are of fundamental importance in many relevant physico-chemical phenomena in the gas phase and at the gas-surface interface. Collision-induced energy transfer not only offers a sensitive probe of interatomic interactions, but also impacts an array of processes ranging from transport to reactivity. This perspective highlights some recent advances on collision-induced energy transfer dynamics. In the gas phase, cold collisions revealed remarkable quantum effects and stereodynamics of bimolecular inelastic scattering, which motivated development of new state-to-state quantum scattering methods. For surface processes, experiments of hydrogen atom scattering from various surfaces combined with first-principles molecular dynamics simulations shed valuable light on the relative importance of adiabatic energy transfer to surface phonons and nonadiabatic energy transfer to surface electrons. We outline further challenges in these fields and expect more fruitful interplay between experiment and theory.

Cold collision↗
DOE OSTI2022

Rapid water dynamics structures the OH-stretching spectra of solitary water in ionic liquids and dipolar solvents

In a recent study [J. Phys. Chem. B 126, 4584–4598 (2022)], we have used infrared spectroscopy to investigate the solvation and dynamics of solitary water in ionic liquids and dipolar solvents. Complex shapes observed for water OH-stretching bands, common to all high-polarity solvents, were assigned to water in several solvation states. Here, in the present study, classical molecular dynamics simulations of a single water molecule in four ionic liquids and three dipolar solvents were used to test and refine this interpretation. Consistent with past assignments, simulations show solitary water usually donates two hydrogen bonds to distinct solvent molecules. Such symmetrically solvated water produces the primary pair of peaks identified in the OH spectra of water in nearly all solvents. We had further proposed that additional features flanking this main peak are due to asymmetric solvation states, states in which only one OH group makes a hydrogen bond to solvent. Such states were found in significant concentrations in all of the systems simulated. Simulations of the OH stretching spectra using a semiclassical description and the vibrational map developed by Auer and Skinner [J. Chem. Phys. 128, 224511–224512 (2008)] provided semi-quantitative agreement with experiment. Analysis of species-specific spectra confirmed assignment of the additional features in the experimental spectra to asymmetrically solvated water. The simulations also showed that rapid water motions cause a marked motional narrowing compared with the inhomogeneous limit. This narrowing is largely responsible for making the additional features due to minority solvation states manifest in the spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗