Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular 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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The Sluggish Diffusion of Cations in CeO 2 Probed through Molecular Dynamics and Metadynamics Simulations

Cation diffusion in fluorite‐structured CeO 2 , though far slower than anion diffusion, is an important, high‐temperature process because it governs diverse fabrication and degradation phenomena. Herein, cation diffusion is studied by means of classical molecular dynamics and metadynamics simulations. Three different mechanisms are examined: migration involving an isolated cerium vacancy, migration involving a cerium vacancy in a defect associate with an oxygen vacancy, and migration involving a cation divacancy. For each mechanism, defect diffusion coefficients are calculated as a function of temperature, from which the respective activation enthalpy of defect migration is obtained. Through comparisons with experimental cation diffusion data (specifically, of the absolute magnitude of the cation diffusivity as well as its activation enthalpy), it is concluded that cation diffusion takes place predominantly neither by isolated vacancies nor by cation vacancy–oxygen vacancy associates but by cation divacancies.

36 MATERIALS SCIENCE↗
NASA NTRS2002

Computational Nanomechanics of Carbon Nanotubes and Composites

Nanomechanics of individual carbon and boron-nitride nanotubes and their application as reinforcing fibers in polymer composites has been reviewed with interplay of theoretical modeling, computer simulations and experimental observations. The emphasis in this work is on elucidating the multi-length scales of the problems involved, and of different simulation techniques that are needed to address specific characteristics of individual nanotubes and nanotube polymer-matrix interfaces. Classical molecular dynamics simulations are shown to be sufficient to describe the generic behavior such as strength and stiffness modulus but are inadequate to describe elastic limit and nature of plastic buckling at large strength. Quantum molecular dynamics simulations are shown to bring out explicit atomic nature dependent behavior of these nanoscale materials objects that are not accessible either via continuum mechanics based descriptions or through classical molecular dynamics based simulations. As examples, we discus local plastic collapse of carbon nanotubes under axial compression and anisotropic plastic buckling of boron-nitride nanotubes. Dependence of the yield strain on the strain rate is addressed through temperature dependent simulations, a transition-state-theory based model of the strain as a function of strain rate and simulation temperature is presented, and in all cases extensive comparisons are made with experimental observations. Mechanical properties of nanotube-polymer composite materials are simulated with diverse nanotube-polymer interface structures (with van der Waals interaction). The atomistic mechanisms of the interface toughening for optimal load transfer through recycling, high-thermal expansion and diffusion coefficient composite formation above glass transition temperature, and enhancement of Young's modulus on addition of nanotubes to polymer are discussed and compared with experimental observations.

Srivastava, Deepak↗
DOE OSTI2020

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↗
NASA NTRS2001

Chemical and Flowfield Modeling for Enhanced Analysis of Contamination Experiments

This paper describes the application of a new Direct Simulation Monte Carlo (DSMC) code, the Molecular Beam Simulator (MBS), which is designed to analyze laboratory scale molecular beam-surface (and crossed-beam) experiments. The MBS is primarily intended to model experiments associated with spacecraft contamination effects, but it can also be used to simulate a variety of surface chemistry and reactive flow measurements. The MBS code is fully three-dimensional, includes a wide-range of chemical processes, and can model one or multiple pulsed (non-steady) sources. As an example application of the MBS code, a fast, pulsed, oxygen atom-surface experiment which examines the chemistry behind erosion of graphite by oxygen atoms is analyzed. Unsteady DSMC simulations show that experimental observations of excited molecular states after the pulse has hit the surface are consistent with two distinct chemical mechanisms: a direct one where the excited molecules are formed on the surface, and a two-step mechanism where ground state molecules formed on the surface are collisionally excited after they leave the surface by trailing oxygen atoms in the pulse. Further DSMC calculations suggest experiments which can distinguish between these mechanisms.

Braunstein, Matthew↗
DOE OSTI2023

Driving and characterizing nucleation of urea and glycine polymorphs in water

Crystal nucleation is relevant across the domains of fundamental and applied sciences. However, in many cases, its mechanism remains unclear due to a lack of temporal or spatial resolution. To gain insights into the molecular details of nucleation, some form of molecular dynamics simulations is typically performed; these simulations, in turn, are limited by their ability to run long enough to sample the nucleation event thoroughly. To overcome the timescale limits in typical molecular dynamics simulations in a manner free of prior human bias, here, we employ the machine learning-augmented molecular dynamics framework “reweighted autoencoded variational Bayes for enhanced sampling (RAVE).” We study two molecular systems—urea and glycine—in explicit all-atom water, due to their enrichment in polymorphic structures and common utility in commercial applications. From our simulations, we observe multiple back-and-forth nucleation events of different polymorphs from homogeneous solution; from these trajectories, we calculate the relative ranking of finite-sized polymorph crystals embedded in solution, in terms of the free-energy difference between the finite-sized crystal polymorph and the original solution state. We further observe that the obtained reaction coordinates and transitions are highly nonclassical.

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

Optimized utilization of COMB3 reactive potentials in LAMMPS

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid–liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. Additionally, the use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

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

Developing reliable machine learning interatomic potential for Fe–Cr–Ni austenitic alloys

Gaining atomistic understanding of mechanical behavior of heat-resistant structural materials such as Fe–Cr–Ni-based alloys requires an approach with an accuracy close to density functional theory (DFT) that considers the intrinsic properties of the bulk lattice and important defects such as stacking faults, grain boundaries, and surfaces. This work aims to develop reliable machine learning interatomic potential (MLIAP) at cross-scale for Fe–Cr–Ni ternary alloys with a focus on the face-centered-cubic (fcc) solid solution structure. Leveraging the advantages of moment tensor potentials, which typically necessitate a relatively small training dataset and enable rapid calculations using the large-scale atomic/molecular massively parallel simulator package, we ensure the stability and accuracy of the trained potentials. Important defects such as stacking faults, grain boundaries, and surfaces for wide-range compositions are investigated. Structural, thermal, elastic, and defect properties are determined from molecular dynamics simulations comprising several thousand atoms, generated via canonical Monte Carlo simulations guided by the trained potential. The trained potential allows efficient atomic simulations of structural, thermal, and mechanical properties of fcc Fe–Cr–Ni solid solution alloys as a function of composition and temperature. Therefore, the MLIAP approach represents a major advancement from DFT calculations that are limited to small simulation sizes and traditional molecular dynamics simulations using relatively low accuracy potentials. Furthermore, this work outlines a practical foundation for further investigating the structural evolution and mechanical behavior of austenitic stainless steel and nickel-based alloys in a wide array of applications in extreme environments.

Crystal structure↗
DOE OSTI2023

Diffusional creep in UO 2 informed by lower length scale simulations

Using molecular dynamics, we predict information at the atomistic scale used to develop a mechanistic UO 2 creep model for use in higher length-scale fuel performance codes. The ultimate objective of the model is to not only to capture the creep rates of UO 2 but to determine the dominant mechanism in the diffusional regime, which is still debated in the literature. It is important to have a model to capture the correct mechanisms for creep in UO 2 as this can be used as the foundation when applying to other fuels, such as doped UO 2 , and when irradiation is accounted. In last years NEAMS milestone (FY22), we developed a prelimnary model, however there were issues, for example, excessively high values of uranium vacancy concentrations at the grain boundary. This year we have addressed the issues with the previous version of the model, added a new term that accounts for the nucleation of dislocations at stress raisers (e.g., triple junctions) within the microstructure and discussed where there was disagreement in the literature about the underpinning physics (uranium self-diffusion at the grain boundary).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Effect of Solvent on the Local Structure, Dynamics, and Vibrational Density of States in Sn-BEA Zeolite

Lewis acid zeolites are attractive catalysts for epoxidation and biomass valorization, as they are highly active and selective in the liquid phase and can operate at or near ambient conditions. While a rich experimental literature exists on liquid-phase Lewis acid zeolite catalysis, our understanding of the molecular organization and solvent dynamics in the vicinity of Lewis acid sites with differing metal site speciation remains limited. In this work, we investigate the molecular coordination and diffusion of two common solvents (methanol and water) around the closed and open Sn-BEA zeolite active sites using molecular dynamics simulations with a machine-learned interatomic potential trained on ab initio molecular dynamics trajectories. Molecular dynamics simulations reveal that introducing active sites significantly enhances local order in the first and second solvation shells compared to the pure silica case. For methanol, both closed and open active sites are singly coordinated, while more than two water molecules coordinate the open site. In contrast to methanol, we observed that water molecules dissociate, leading to the formation of additional Sn-OH and silanol groups away from the active site. The diffusion coefficients of water and methanol are functions of the solvent population in the pore. Here, our work provides insights into how active site speciation in Lewis acid zeolites affects solvent coordination, diffusion, and vibrational signature. This information is foundational for catalyst design and optimization of liquid-phase catalytic processes in zeolites. It also demonstrates the suitability of machine-learned interatomic potentials for modeling reactive systems, enabling sufficiently long trajectories for appropriate statistical averaging.

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

Effect of slip on vortex shedding from a circular cylinder in a gas flow

Most studies of vortex shedding from a circular cylinder in a gas flow have explicitly or implicitly assumed that the no-slip condition applies on the cylinder surface. To investigate the effect of slip, vortex shedding is simulated using molecular gas dynamics (the direct simulation Monte Carlo method) and computational fluid dynamics (the incompressible Navier-Stokes equations with a slip boundary condition). A Reynolds number of 100, a Mach number of 0.3, and a corresponding Knudsen number of 0.0048 are examined. For these conditions, compressibility effects are small, and periodic laminar vortex shedding is obtained. Slip on the cylinder is varied using combinations of diffuse and specular molecular reflections with accommodation coefficients from zero (maximum slip) to unity (minimum slip). Although unrealistic, bounce-back molecular reflections are also examined because they approximate the no-slip boundary condition (zero slip). The results from both methods are in reasonable agreement. The shedding frequency increases slightly as the accommodation coefficient is decreased, and shedding ceases at low accommodation coefficients (large slip). The streamwise and transverse forces decrease as the accommodation coefficient is decreased. Based on the good agreement between the two methods, computational fluid dynamics is used to determine the critical accommodation coefficient below which vortex shedding ceases for Reynolds numbers of 60–100 at a Mach number of 0.3. Finally, conditions to observe the effect of slip on vortex shedding appear to be experimentally realizable, although challenging.

42 ENGINEERING↗
DOE OSTI2020

Atomistic Characterization of Gramicidin Channel Formation

We investigated gramicidin A (gA) subunit dimerization in lipid bilayers using microsecond-long replica-exchange umbrella sampling simulations, millisecond-long unbiased molecular dynamics simulations, and machine learning. Our simulations led to a dimer structure that is indistinguishable from the experimentally determined gA channel structures, with the two gA subunits joined by six hydrogen bonds (6HB). The simulations also uncovered two additional dimer structures, with different gA–gA stacking orientations that were stabilized by four or two hydrogen bonds (4HB or 2HB). When examining the temporal evolution of the dimerization, we found that two bilayer-inserted gA subunits can form the 6HB dimer directly, with no discernible intermediate states, as well as through paths that involve the 2HB and 4HB dimers.

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

Interfacial Response and Structural Adaptation of Structured Polyelectrolyte Thin Films

Ionizable block copolymers with distinctive block characteristics display the diversity crucial for the design of macromolecules for targeted applications. In contrast to van der Waals copolymers, their interfaces, which are critical to their function, consist of nanodomains, each of a different nature and thus unique interfacial behavior. Here, the interfacial response of a symmetric block copolymer with a sulfonated polystyrene polyelectrolyte center, tethered to polyethylene-r-propylene and terminated by poly(t-butyl styrene) is probed as polymer films are exposed to three polar solvents, water, propanol, and tetrahydrofuran (THF), using molecular dynamics simulations. Each of the solvents captures a distinctive interaction with the individual blocks. We find that at the film boundary, the interfacial response is initially dominated by that of the hydrophobic blocks to all solvents. At later times, the solvent distribution among the blocks, where water molecules associate predominantly with the sulfonated groups and propanol and THF reside at multiple different sites, determines the chemical composition and the polymer conformation at the interface. Altogether, these simulations provide the first direct molecular insight into the interfacial response of ionizable copolymers.

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

On the mechanisms of ion adsorption to aqueous interfaces: air-water vs. oil-water

The adsorption of ions to water-hydrophobe interfaces influences a wide range of phenomena, including chemical reaction rates, ion transport across biological membranes, and electrochemical and many catalytic processes; hence, developing a detailed understanding of the behavior of ions at water-hydrophobe interfaces is of central interest. Here, we characterize the adsorption of the chaotropic thiocyanate anion (SCN − ) to two prototypical liquid hydrophobic surfaces, water-toluene and water-decane, by surface-sensitive nonlinear spectroscopy and compare the results against our previous studies of SCN − adsorption to the air-water interface. For these systems, we observe no spectral shift in the charge transfer to solvent spectrum of SCN − , and the Gibb’s free energies of adsorption for these three different interfaces all agree within error. We employed molecular dynamics simulations to develop a molecular-level understanding of the adsorption mechanism and found that the adsorption for SCN − to both water-toluene and water-decane interfaces is driven by an increase in entropy, with very little enthalpic contribution. This is a qualitatively different mechanism than reported for SCN − adsorption to the air-water and graphene-water interfaces, wherein a favorable enthalpy change was the main driving force, against an unfavorable entropy change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2004

Cross-Linked Nanotube Materials with Variable Stiffness Tethers

The constitutive properties of a cross-linked single-walled carbon nanotube material are predicted with a multi-scale model. The material is modeled as a transversely isotropic solid using concepts from equivalent-continuum modeling. The elastic constants are determined using molecular dynamics simulation. Some parameters of the molecular force field are determined specifically for the cross-linker from ab initio calculations. A demonstration of how the cross-linked nanotubes may affect the properties of a nanotube/polyimide composite is included using a micromechanical analysis.

Frankland, Sarah-Jane V.↗
DOE OSTI2025

Shock-induced chemistry and high strain-rate viscoelastic behavior of a phenolic polymer

We use impact experiments and a finite element model (up to 1.2 GPa), and molecular dynamics simulations (up to 60 GPa), to examine the behavior of a phenolic polymer under shock compression, spanning both nonreactive and reactive regimes. In the nonreactive regime, relaxation following compression at strain rates of ∼105 s−1 can be explained by viscoelasticity observed at ordinary laboratory rates (≲1 s−1) by accounting for the temperature dependence of the phenolic β-transition. Reasonable agreement is found between the measured shock Hugoniot up to 1.2 GPa and molecular dynamics simulation for cross-linked structures of comparable density. We also observed a first-order mechanical transition near 0.36 GPa shock stress and estimated a spall strength of 0.102 GPa and Hugoniot elastic limit of 1–2 GPa. The shock stress is found to vary up to 24% among phenolics made with different resin and/or cure processes. Finally, molecular dynamics simulations are used to identify a reactive regime at shock pressures ≳20 GPa that is characterized by chemically driven, rate-dependent relaxation processes, including dehydrogenation and dehydration reactions that promote the formation of a dense, highly cross-linked carbonaceous solid and the release of light volatiles.

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

Computational study of the electrostatic potential and charges of multivalent ionic liquid molecules

The electrostatic interactions in ionic liquids (ILs) can be difficult to quantify, but this information is very valuable for improving intermolecular potentials for performing molecular-level simulations. Although polarizable forcefields can provide improved accuracy, the computational efficiency of fixed charge molecular models is often preferred for larger systems. Here, we propose a multi-scale screening approach for analyzing the charge scaling behavior and other electrostatic properties of ILs, which is more attainable than previous analysis methods that rely on the IL crystalline structures. Based on isolated molecules and small clusters, we study the electrostatic potential properties and charges of multivalent IL molecules using density-functional theory calculations and ab initio molecular dynamics simulations. The charges of the IL molecules are calculated and compared using several different charge estimation techniques. Although there are some differences among the estimated partial charges on the atoms, the total anion and cation charges are very similar. In conclusion, the charge scaling factors are substantial, and they are only slightly affected by the IL cluster size model.

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

CTRAMER: An open-source software package for correlating interfacial charge transfer rate constants with donor/acceptor geometries in organic photovoltaic materials

In this paper, we present CTRAMER (Charge-Transfer RAtes from Molecular dynamics, Electronic structure, and Rate theory)—an opensource software package for calculating interfacial charge-transfer (CT) rate constants in organic photovoltaic (OPV) materials based on ab initio calculations and molecular dynamics simulations. The software is based on identifying representative donor/acceptor geometries within interfacial structures obtained from molecular dynamics simulation of donor/acceptor blends and calculating the corresponding Fermi’s golden rule CT rate constants within the framework of the linearized-semiclassical approximation. While the methods used are well established, the integration of these state-of-the-art tools originating from different disciplines to study photoinduced CT processes with explicit treatment of the environment, in our opinion, makes this package unique and innovative. The software also provides tools for investigating other observables of interest. After outlining the features and implementation details, the usage and performance of the software are demonstrated with results from an example OPV system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗