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 793 records · Page 44

Probing electrolyte–silica interactions through simulations of the infrared spectroscopy of nanoscale pores

The structural and dynamical properties of nanoconfined solutions can differ dramatically from those of the corresponding bulk systems. Understanding the changes induced by confinement is central to controlling the behavior of synthetic nanostructured materials and predicting the characteristics of biological and geochemical systems. A key outstanding issue is how the molecular-level behavior of nanoconfined electrolyte solutions is reflected in different experimental, particularly spectroscopic, measurements. This is addressed here through molecular dynamics simulations of the OH stretching infrared (IR) spectroscopy of NaCl, NaBr, and NaI solutions in isotopically dilute HOD/D 2 O confined in hydroxylated amorphous silica slit pores of width 1–6 nm and pH ~2. In addition, the water reorientation dynamics and spectral diffusion, accessible by pump–probe anisotropy and two-dimensional IR measurements, are investigated. In this work, the aim is to elucidate the effect of salt identity, confinement, and salt concentration on the vibrational spectra. It is found that the IR spectra of the electrolyte solutions are only modestly blue-shifted upon confinement in amorphous silica slit pores, with both the size of the shift and linewidth increasing with the halide size, but these effects are suppressed as the salt concentration is increased. This indicates the limitations of linear IR spectroscopy as a probe of confined water. However, the OH reorientational and spectral diffusion dynamics are significantly slowed by confinement even at the lowest concentrations. The retardation of the dynamics eases with increasing salt concentration and pore width, but it exhibits a more complex behavior as a function of halide.

Mesoporous material↗

Heteroepitaxy of FCC-on-FCC systems of large misfit

Here, to understand the effects of lattice mismatch on heteroepitaxial growth, we have studied the equilibrium structure and orientation relationships (ORs) of FCC films grown epitaxially on FCC substrates, using molecular dynamics simulations in conjunction with embedded atom method potentials. Three film/substrate systems have been investigated, namely: Ag on Cu, Ag on Ni and Pb on Al. These systems cover a significant range of lattice mismatch, from 12.6% for Ag/Cu to 21.8% for Pb/Al. For each system, the ORs of films on six different substrate orientations, namely: (100), (511), (311), (211), (322) and (111), have been investigated. Films on these susbstrates cover a gradual transition from the oct-cube orientation relationship, which occurs only on (100) substrates, to the heterotwin orientation relationship, which often occurs on (111) substrates. It is found that the resulting ORs vary systematically with substrate orientation, but that the pattern of variation is almost identical for all three systems, and therefore largely independent of mismatch. However, the manner in which mismatch is accommodated does depend on the magnitude of mismatch. Simulations point to an important role for edge-to-edge matching and defects such as stacking faults. An analysis of these results in terms of transformation strains highlights the distinction between the ORs, which are largely independent of mismatch, and the local interfacial structure, which changes directly with mismatch.

36 MATERIALS SCIENCE↗

A molecular dynamics survey study of impact-induced amorphization in yttria nanoparticles

The impact of a single yttria (Y 2 O 3 ) nanoparticle onto a yttria substrate was studied as a function of particle velocity (300–1200 m/s) and diameter (5–12 nm) for three particle orientations and two substrate orientations using molecular dynamics simulations. Impacts resulted in minimal dislocation motion but extensive disordering of the crystal lattice without approaching the melting temperature. The resulting amorphous microstructure was both precisely coincident with plastic strain and retained until the end of the simulation, suggesting a deformation mechanism of amorphization-facilitated viscous flow. Further analysis revealed that the degree of amorphization depended primarily on the impact velocity and secondarily on the orientations of the particle and substrate, with significant second-order effects from variable interactions.

74 ATOMIC AND MOLECULAR PHYSICS↗

Polyphenols Weaken Pea Protein Gel by Formation of Large Aggregates with Diminished Noncovalent Interactions

Polyphenols are well-known native cross-linkers and gel strengthening agents for many animal proteins. However, their role in modifying plant protein gels remains unclear. In this study, multiple techniques were applied to unravel the influence of green tea polyphenols (GTP) on pea protein gels and the underlying mechanisms. We found that the elasticity and viscosity of pea protein gels decreased with increased GTP. The protein backbone became less rigid when GTP was present based on shortened T 1ρ H in relaxation solid-state NMR measurements. Electron microscopy and small-angle X-ray scattering showed that gels weakened by GTP possessed disrupted networks with the presence of large protein aggregates. Solvent extraction and molecular dynamic simulation revealed a reduction in hydrophobic interactions and hydrogen bonds among proteins in gels containing GTP. The current findings may be applicable to other plant proteins for greater control of gel structures in the presence of polyphenols, expanding their utilization in food and biomedical applications.

59 BASIC BIOLOGICAL SCIENCES↗

Excess Density as a Descriptor for Electrolyte Solvent Design

Electrolytes mediate interactions between the cathode and anode and determine performance characteristics of batteries. Mixtures of multiple solvents are often used in electrolytes to achieve desired properties, such as viscosity, dielectric constant, boiling point, and melting point. Conventionally, multi-component electrolyte properties are approximated with linear mixing, but in practice, significant deviations are observed. Excess quantities can provide insights into the molecular behavior of the mixture and could form the basis for designing high-performance electrolytes. Here we investigate the excess density of commonly used Li-ion battery solvents such as cyclic carbonates, linear carbonates, ethers, and nitriles with molecular dynamics simulations. We additionally investigate electrolytes consisting of these solvents and a salt. The results smoothly vary with mole percent and are fit to permutation-invariant Redlich-Kister polynomials. Mixtures of similar solvents, such as cyclic-cyclic carbonate mixtures, tend to have excess properties that are lower in magnitude compared to mixtures of dissimilar substances, such as carbonate-nitrile mixtures. We perform experimental testing using our automated test stand, Clio, to provide validation to the observed simulation trends. We quantify the structure similarity using SOAP fingerprints to create a descriptor for excess density, enabling the design of electrolyte properties. To a first approximation, this will allow us to estimate the deviation of a mixture from ideal behavior based solely upon the structural dissimilarity of the components.

74 ATOMIC AND MOLECULAR PHYSICS↗

Stokesian dynamics with odd viscosity

Stokesian dynamics is a well-established computational method for simulating dynamics of many particles suspended in a conventional passive fluid medium. Active fluids composed of self-driving particles with broken time-reversal symmetry permit the emergence of a so-called odd viscosity. Here, in this work, we extended the conventional Stokesian dynamics formalism to incorporate the additional hydrodynamic effects due to odd viscosity, which enables simulating collective behavior of many particles suspended in an active fluid medium with both even viscosity and odd viscosity.

42 ENGINEERING↗

Disclination-dislocation based model for grain boundary stress field evolution due to slip transmission history and influence on subsequent dislocation transmission

This work demonstrates how the structure of a grain boundary (GB) and its evolution due to slip transmission history influences subsequent dislocation transmission. First, a model for the evolution of stress fields within grain boundaries that accounts for the effects of coherent dislocation transmission is introduced. Starting with a disclination-based construct of GBs at minimum energy (equilibrium), the model describes the evolution of the GB stress field to a state characteristic of excess energy (non-equilibrium) due to the incorporation of residual Burgers vector content following sequential slip transmission events. Several essential features of this model are verified via molecular dynamics simulations of lattice dislocation absorption. Second, this model is implemented into a discrete dislocation dynamics (DDD) code and simulations are performed to understand the influence of Non-equilibrium GB stress fields, conditioned by the slip transmission history, on subsequent dislocation transmission. DDD simulations reveal that the critical resolved stress necessary for slip propagation can be reduced by with continued absorption of residual dislocation content. Moreover, DDD simulations prove that a comprehensive consideration of both the binding and driving stresses, and the evolution of the transmission configuration is necessary to quantify the influence of the mechanical state of the GB on slip propagation. In general, this work provides important insights into the role of GB structure evolution, conditioned by prior deformation history, on intergranular plasticity.

42 ENGINEERING↗

Ion-gated carbon molecular sieve gas separation membranes

Membrane technology lies at the heart of many industrial gas separation processes and applications. Molecular sieving membranes that break the Robeson limit are desirable for energy-efficient gas separation. In this paper, we report a facile strategy of directly integrating ionic liquids (ILs) into porous membranes. Particularly, the ILs form an ultra-thin layer on the carbon molecular sieve (CMS) membranes rather than penetrating into the pores, acting as a smart gate for gas entry to boost the selectivity. The hybrid membrane exhibits CO 2 permeability >600 barrer and enhanced CO 2 /N 2 selectivity >50, which surpasses the Robeson limit and shows potential in CO 2 /N 2 separation process. Molecular dynamics simulations confirm the gating effect of the IL layer of molecular thickness. In conclusion, this work demonstrates a universal strategy to improve CMS membrane performance by creating an IL-membrane interface and tuning the ion-pore interaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrokinetic comparison of CaF 2 and AgCl from a MDS perspective

CaF 2 has limited electrokinetic sensitivity to fluoride ions in solution and the positive CaF 2 surface charge is not reversed at high fluoride ion concentration. Similar insensitivity of surface charge to high fluoride ion concentration has been observed for other alkali and alkaline earth fluoride salts. On the contrary, chloride ions act as potential determining ions for AgCl and reverse the surface charge, as expected. FTIR-IRS results indicate that the AgCl/water interface has a water structure-breaking tendency, whereas the CaF 2 /water interface has a structure-making tendency. Here, it is speculated that the strong hydrogen bonding between fluoride ions and water prevents the accommodation of excess fluoride ions at surface lattice sites of fluorite, and therefore explains the low sensitivity of the fluorite surface charge sign and magnitude to fluoride ion concentration. However, this speculation has not been examined at the molecular level and is evaluated by molecular dynamics simulations (MDS) in the research presented in this paper. Higher hydrogen bond density at the CaF 2 mineral/water interface confirms that the water structure has stronger hydrogen bonding at both the CaF 2 mineral surface and with the F - ion in solution. In this way, the interaction of the F - ion with the CaF 2 surface is inhibited.

36 MATERIALS SCIENCE↗

Sienna Modeling Framework [Slides]

NREL's Sienna modeling framework effectively builds, solves, and analyzes the scheduling problems and dynamic simulations of quasi-static infrastructure systems. It uses a modular framework to answer different questions about future energy systems, fundamentally advancing the nation's ability to model individual and integrated infrastructure systems at a range of spatial and temporal scales. This presentation will include NREL power grid researcher Clayton Barrows.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Rational design of antimicrobial peptides targeting Gram-negative bacteria

Membrane-targeting host antimicrobial peptides (AMPs) can kill or inhibit the growth of Gram-negative bacteria. However, the evolution of resistance among microbes poses a substantial barrier to the long-term utility of the host AMPs. Combining experiment and molecular dynamics simulations, we demonstrate that terminal carboxyl capping enhances both membrane insertion and antibacterial activity of an AMP called P1. Furthermore, we show that a bacterial strain with evolved resistance to this peptide becomes susceptible to P1 variants with either backbone capping or lysine-to-arginine substitutions. Our results suggest that cocktails of closely related AMPs may be useful in overcoming evolved resistance.

59 BASIC BIOLOGICAL SCIENCES↗

Thermal coupling-decoupling mechanism of heat transfer across van der Waals interfaces in n-eicosane

Progress toward improving thermal transport properties of n-eicosane, a promising phase change material for thermal energy storage near human body temperature, is hindered by lack of a general theoretical framework to model thermal interfacial conductance (TIC) of the van der Waals (vdW) force bonded molecular interfaces, which have little compositional/structural mismatch. Combining molecular dynamics simulations at temperatures up to the melting point and a time-domain analysis of interfacial heat currents in individual molecules, we unveil that (1) the heat flux across a molecular vdW interface has an alternating coupling-decoupling pattern at the single molecule level due to thermal fluctuations, and (2) the net heat transfer during the thermally coupled periods is orders of magnitude larger than that during the thermally decoupled periods. Furthermore, we demonstrate that a vdW interface's TIC correlates almost linearly with two atomistic-scale coupling-decoupling parameters: αc and σc, that represent the duration and strength of the thermal coupling, respectively. This thermal coupling-decoupling mechanism proposes a new paradigm that explicitly accounts for the effects of structural flexibility and dynamical fluctuation on heat transfer across molecular interfaces, and will likely provide new insights into modeling thermal transport properties in a wide-ranging soft materials, including self-assembled monolayers, paraffins or lipids.

42 ENGINEERING↗

Probing DNA structural heterogeneity by identifying conformational subensembles of a bicovalently bound cyanine dye

DNA is a re-configurable, biological information-storage unit, and much remains to be learned about its heterogeneous structural dynamics. For example, while it is known that molecular dyes templated onto DNA exhibit increased photostability, the mechanism by which the structural dynamics of DNA affect the dye photophysics remains unknown. Here, we use femtosecond, two-dimensional electronic spectroscopy measurements of a cyanine dye, Cy5, to probe local conformations in samples of single-stranded DNA (ssDNA–Cy5), double-stranded DNA (dsDNA–Cy5), and Holliday junction DNA (HJ–DNA–Cy5). A line shape analysis of the 2D spectra reveals a strong excitation–emission correlation present in only the dsDNA–Cy5 complex, which is a signature of inhomogeneous broadening. Molecular dynamics simulations support the conclusion that this inhomogeneous broadening arises from a nearly degenerate conformer found only in the dsDNA–Cy5 complex. These insights will support future studies on DNA’s structural heterogeneity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

An Environmentally Benign Electrolyte for High Energy Lithium Metal Batteries

In this work, a hybrid electrolyte comprising a high content of H 2 O for a lithium metal cell is reported. At high LiFSI salt concentration, the N-methyl-N-propyl-piperidinium bis(fluorosulfonyl) imide (PMpipFSI) electrolyte can tolerate up to 1 M H2O addition without sacrificing its redox stability on both lithium nickel manganese cobalt oxide (NMC) cathode and lithium metal anode. Molecular dynamics simulations revealed the underpinned mechanism that, at high salt concentrations, H 2 O molecules are embedded in the Li + , PMpip + , and FSI - bulk as a structural material with a strong solvation with Li + and are orderly distributed at the surface of both electrodes. This electrolyte eliminates the critical moisture controls required for the state-of-the-art (SOA) carbonate/LiPF 6 electrolyte, electrode, separator and cell assembly, thus significantly reducing the cost of the mass production of the batteries.

molecular dynamics simulation↗

Molecular Dynamics Investigation of Lipid-Specific Interactions with a Fusion Peptide

The HIV-1 fusion peptide, which is a short hydrophobic peptide from the gp41 coat glycoprotein that participates in the infection of a cell, interacts with model lipid bilayer membranes in a concentration-dependent manner. The interaction of the peptide with the bilayer also strongly depends on the lipid composition. Here, molecular dynamics simulations were performed to investigate lipid-specific interactions that arise shortly after the binding of a less-fusogenic variant of the HIV-1 fusion peptide to a lipid bilayer composed of a mixture of dimyristoyl phosphatidylcholine and dimyristoyl phosphatidylglycerol. The impact of peptide concentration was also studied. An improved understanding was gained of the lipid-specific interactions experienced by the FP. New insight was also gained into how the peptide concentration changes these interactions.

60 APPLIED LIFE SCIENCES↗

Effects of high-pressure hydrogen exposure on filler-elastomer adhesion

Elastomers are known to gain enhanced mechanical properties through compounding with nanosized filler particles such as silica or carbon black. Filler dispersion and filler-polymer interfacial strength are key contributing factors to this improvement. The interfacial strength is critical to part lifetime in pressurized gas sealing applications such as O-rings, where weak binding between the filler particle and polymer matrix can lead to internal void structures. With the aim to build a fundamental understanding of precursors to pressurized hydrogen-induced failure in elastomers, we use all-atom molecular dynamics simulations to study the impact of hydrogen oversaturation on filler-polymer interaction strength. We systematically study the interface between a commonly used elastomer, ethylene-propylene-diene monomer (EPDM) and silica by varying gas concentration, crosslink density, and surface chemistry. Our simulations predict that decompression leads to a localization of excess gas near the interface. In conclusion, we demonstrate that this localized gas can weaken interfacial adhesion and quantify the interaction using thermodynamic approaches.

EPDM↗

Structure and dynamics of SARS-CoV-2 proofreading exoribonuclease ExoN

High-fidelity replication of the large RNA genome of coronaviruses (CoVs) is mediated by a 3'-to-5' exoribonuclease (ExoN) in nonstructural protein 14 (nsp14), which excises nucleotides including antiviral drugs misincorporated by the low-fidelity viral RNA-dependent RNA polymerase (RdRp) and has also been implicated in viral RNA recombination and resistance to innate immunity. Here, we determined a 1.6-Å resolution crystal structure of severe acute respiratory syndrome CoV 2 (SARS-CoV-2) ExoN in complex with its essential cofactor, nsp10. The structure shows a highly basic and concave surface flanking the active site, comprising several Lys residues of nsp14 and the N-terminal amino group of nsp10. Modeling suggests that this basic patch binds to the template strand of double-stranded RNA substrates to position the 3' end of the nascent strand in the ExoN active site, which is corroborated by mutational and computational analyses. We also show that the ExoN activity can rescue a stalled RNA primer poisoned with sofosbuvir and allow RdRp to continue its extension in the presence of the chain-terminating drug, biochemically recapitulating proofreading in SARS-CoV-2 replication. Molecular dynamics simulations further show remarkable flexibility of multidomain nsp14 and suggest that nsp10 stabilizes ExoN for substrate RNA binding to support its exonuclease activity. Our high-resolution structure of the SARS-CoV-2 ExoN–nsp10 complex serves as a platform for future development of anticoronaviral drugs or strategies to attenuate the viral virulence.

60 APPLIED LIFE SCIENCES↗

Direct Correlation of the Salt-Reduced Diffusivities of Organic Solvents with the Solvent’s Mole Fraction

We report that lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in organic solvents (especially propylene carbonate) has demonstrated extraordinary pseudocapacitive performance as an electrolyte in the supercapacitor configuration. However, the influence of the solvated ions on the diffusivity of the solvent molecules is yet to be understood. We examine the impact of LiTFSI on the diffusivity in five organic solvents: acetonitrile (ACN), tetrahydrofuran (THF), methanol (MeOH), dimethyl sulfoxide (DMSO), and propylene carbonate (PC) using a combination of neutron scattering, conductivity measurements, and molecular dynamics simulations. The extent of the diffusivity reduction in the concentration regime of ≤1 M directly correlates with the solvent mole fraction at which the solvation shells around Li + ions are of similar size in all the solvents, resulting in a universal ~50% reduction in the solvent diffusivity. These results provide guidance for formulation of the new electrolytes to enhance the performance of energy storage devices.

conductivity↗