SEARCH · Engineering Papers

Results for “dynamics 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Molecular Dynamics Simulation of Resin Adsorption at Kaolinite Edge Sites: Effect of Surface Deprotonation on Interfacial Structure

Low-salinity water flooding, a method of enhanced oil recovery, consists of injecting low ionic strength fluids into an oil reservoir in order to detach oil from mineral surfaces in the underlying formation. Although highly successful in practice, the approach is not completely understood at the molecular scale. Molecular dynamics simulations have been used to investigate the effect of surface protonation on the adsorption of an anionic crude oil component on clay mineral edge surfaces. A set of interatomic potentials appropriate for edge simulations has been applied to the kaolinite (010) surface in contact with an aqueous nanopore. Decahydro-2-napthoic acid in its deprotonated form (DHNA – ) was used as a representative resin component of crude oil, with monovalent and divalent counterions, to test the observed trends in low-salinity water flooding experiments. Surface models include fully protonated (neutral) and deprotonated (negative) edge sites, which require implementation of a new deprotonation scheme. The surface adsorptive properties of the kaolinite edge under neutral and deprotonated conditions have been investigated for low and high DHNA – concentrations with Na + and Ca 2+ as counterions. The tendency of DHNA – ions to coordinate with divalent (Ca 2+ ) rather than monovalent (Na + ) ions greatly influences adsorption tendencies of the anion. Additionally, the formation of net positively charged surface sites due to Ca 2+ at deprotonated sites results in increased DHNA – adsorption. Divalent cations such as Ca 2+ are able to efficiently bridge surface sites and organic anions. Replacing those cations with monovalent cations such as Na + diminishes the bridging mechanism, resulting in reduced adsorption of the organic species. As a result, a clear trend of decreased DHNA – adsorption is observed in the simulations as Ca 2+ is replaced by Na + for deprotonated surfaces, as would be expected for oil detachment from reservoir formations following a low-salinity flooding event.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2014

Molecular dynamics simulations of uranyl adsorption and structure on the basal surface of muscovite

Anthropogenic activities have led to an increased concentration of uranium on the Earth’s surface and potentially in the subsurface with the development of nuclear waste repositories. Uranium is soluble in groundwater, and its mobility is strongly affected by the presence of clay minerals in soils and in subsurface sediments. We use molecular dynamics simulations to probe the adsorption of aqueous uranyl (UO 2 2+ ) ions onto the basal surface of muscovite, a suitable proxy for typically ultrafine-grained clay phases. Model systems include the competitive adsorption between potassium counterions and aqueous ions (0.1 M and 1.0 M UO 2 Cl 2 , 0.1 M NaCl). We find that for systems with potassium and uranyl ions present, potassium ions dominate the adsorption phenomenon. Potassium ions adsorb entirely as inner-sphere complexes associated with the ditrigonal cavity of the basal surface. Uranyl ions adsorb in two configurations when it is the only ion species present, and in a single configuration in the presence of potassium. Finally, the majority of adsorbed uranyl ions are tilted less than 45° relative to the muscovite surface, and are associated with the Si 4 Al 2 rings near aluminum substitution sites.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2018

A fast scheme to calculate electronic couplings between P3HT polymer units using diabatic orbitals for charge transfer dynamics simulations

We propose a fast and accurate calculation method to compute the electronic couplings between molecular units in a thiophene-ring-based polymer chain mimicking a real organic semiconducting polymer, poly(3-hexylthiophene). Through a unit block diabatization scheme, the method employed minimal number of diabatic orbitals to compute the site energies and electronic couplings, which were validated by comparing with benchmark density functional theory calculations. In addition, by using the obtained electronic couplings, a quantum dynamics simulation was carried out to propagate a hole initially localized in a thiophene-ring unit of the polymer chain. Here, this work establishes a simple, efficient, and robust means for the simulation of electron or hole transfer processes in organic semiconducting materials, an important capability for study and understanding of the class of organic optoelectronic and photovoltaic materials.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Formation and dissociation of shear-induced high-energy dislocations: insight from molecular dynamics simulations

Solid-phase processing (SPP) allows one to create complex microstructures, not achievable via thermal processing alone. The resulting structures exhibit a rich palette of defects, both thermal and non-thermal, including defect substructures, such as dislocation networks. It is essential to understand the mechanisms of deformation and defect structure formation to guide SPP towards achieving desired microstructures and material properties. In this study, large-scale molecular dynamics simulations are used to investigate the effects of inhomogeneous strain distribution, that mimics deformation conditions of tribological tests, on the evolution of defects under severe shear deformation in polycrystalline Al. Analysis of defect nucleation and reaction pathways reveals that strong geometric constraints suppress the nucleation and slide of low energy dislocation 1/2$\langle110\rangle${111} but promote the nucleation and slide of high energy dislocations, such as [$1\bar1$0] (001) and 1/2 [$1\bar1$$\bar2$]($1\bar1$1). Further, a rough contact surface, characteristic to tribological tests, imposes an inhomogeneous stress field leading to inhomogeneous defect substructures due to location-dependent activation of slip systems. The results suggest that high-energy dislocations can dominate the evolution of grain structures in highly constrained environments, which should be considered in modeling plastic deformation and grain refinement during SPP.

42 ENGINEERING↗

DOE OSTI2016

Effect of sequence dispersity on morphology of tapered diblock copolymers from molecular dynamics simulations

Tapered diblock copolymers are similar to typical AB diblock copolymers but have an added transition region between the two blocks which changes gradually in composition from pure A to pure B. This tapered region can be varied from 0% (true diblock) to 100% (gradient copolymer) of the polymer length, and this allows some control over the microphase separated domain spacing and other material properties. We perform molecular dynamics simulations of linearly tapered block copolymers with tapers of various lengths, initialized from fluids density functional theory predictions. To investigate the effect of sequence dispersity, we compare systems composed of identical polymers, whose taper has a fixed sequence that most closely approximates a linear gradient, with sequentially disperse polymers, whose sequences are created statistically to yield the appropriate ensemble average linear gradient. Especially at high segregation strength, we find clear differences in polymer conformations and microstructures between these systems. Importantly, the statistical polymers are able to find more favorable conformations given their sequence, for instance, a statistical polymer with a larger fraction of A than the median will tend towards the A lamellae. The conformations of the statistically different polymers can thus be less stretched, and these systems have higher overall density. Consequently, the lamellae formed by statistical polymers have smaller domain spacing with sharper interfaces

36 MATERIALS SCIENCE↗

DOE OSTI2022

Predicting the Mechanical Response of Polyhydroxyalkanoate Biopolymers Using Molecular Dynamics Simulations

Polyhydroxyalkanoates (PHAs) have emerged as a promising class of biosynthesizable, biocompatible, and biodegradable polymers to replace petroleum-based plastics for addressing the global plastic pollution problem. Although PHAs offer a wide range of chemical diversity, the structure–property relationships in this class of polymers remain poorly established. In particular, the available experimental data on the mechanical properties is scarce. In this contribution, we have used molecular dynamics simulations employing a recently developed forcefield to predict chemical trends in mechanical properties of PHAs. Specifically, we make predictions for Young’s modulus, and yield stress for a wide range of PHAs that exhibit varying lengths of backbone and side chains as well as different side chain functional groups. Deformation simulations were performed at six different strain rates and six different temperatures to elucidate their influence on the mechanical properties. Our results indicate that Young’s modulus and yield stress decrease systematically with increase in the number of carbon atoms in the side chain as well as in the polymer backbone. In addition, we find that the mechanical properties were strongly correlated with the chemical nature of the functional group. The functional groups that enhance the interchain interactions lead to an enhancement in both the Young’s modulus and yield stress. Finally, we applied the developed methodology to study composition-dependence of the mechanical properties for a selected set of binary and ternary copolymers. Overall, our work not only provides insights into rational design rules for tailoring mechanical properties in PHAs, but also opens up avenues for future high throughput atomistic simulation studies geared towards identifying functional PHA polymer candidates for targeted applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Computing virtual dark-field X-ray microscopy images of complex discrete dislocation structures from large-scale molecular dynamics simulations

Dark-field X-ray microscopy (DFXM) is a novel diffraction-based imaging technique that non-destructively maps the local deformation from crystalline defects in bulk materials. While studies have demonstrated that DFXM can spatially map 3D defect geometries, it is still challenging to interpret DFXM images of the high-dislocation-density systems relevant to macroscopic crystal plasticity. This work develops a scalable forward model to calculate virtual DFXM images for complex discrete dislocation structure(s) (DDS) obtained from atomistic simulations. Our new DDS-DFXM model integrates a non-singular formulation for calculating the local strain from the DDS and an efficient geometrical optics algorithm for computing the DFXM image from the strain field. We apply the model to complex DDS obtained from a large-scale molecular dynamics simulation of compressive loading on single-crystal silicon. Simulated DFXM images exhibit prominent contrast for dislocation features between the multiple slip systems, demonstrating the potential of DFXM to resolve features from dislocation multiplication. In conclusion, the integrated DDS-DFXM model provides a toolbox for DFXM experimental design and image interpretation in the context of bulk crystal plasticity for a range of measurements across shock plasticity and the broader materials science community.

X-ray imaging↗

DOE OSTI2019

Structural features of sodium silicate glasses from reactive force field-based molecular dynamics simulations

Atomistic computer simulations can offer insights into silicate glass-environment interactions with the recent development of reactive potentials. Yet, the accuracy of generated glass structures with these potential was usually not fully examined. In this paper, the capability of the reactive force field (ReaxFF) to describe the short and medium range structure features of sodium silicate glasses in molecular dynamics simulations is investigated by comparing a widely used partial charge pairwise potential and available experimental data. Glass structure information such as pair distribution function (PDF), coordination number, Qn species, neutron broadened structure factor, and X-ray broadened structure factor of the glass structures from ReaxFF simulations were calculated and compared to evaluate the generated glass structure. Advantages and limitations of the potentials and glass forming procedures, as well as areas of further improvement, were discussed. The results show that the recently refined ReaxFF parameters through the proposed procedure enable the simulations of sodium silicate glass structures with minimal defects, which paves the way to investigate water-glass interaction mechanisms with the reactive enabled potentials.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Effect of Tethering Anions in Block Copolymer Electrolytes via Molecular Dynamics Simulations

Ion-containing block copolymers (BCPs) are promising solid-state electrolytes as they can combine ion conduction and mechanical robustness via multiple microphases. However, salt-doped BCPs can have low ion conduction and transference number. One strategy to increase transference number is to tether anions to the polymer backbone (creating a single-ion conductor), though this may impact cation mobility. Leveraging coarse-grained molecular dynamics simulations, we compare the effects of ion concentration and polymer dielectric constant on BCPs with and without bonding anions to the chain (without changing any other chemical properties). We find that single-ion BCPs have higher cation conductivity than analogous salt-doped BCPs, especially at moderate ion concentrations. Furthermore, we also compare several structural features of these systems; notably, tethering anions causes the ions to spread more evenly throughout the conducting phase, which is related to improved cation conduction.

25 ENERGY STORAGE↗

NASA NTRS2010

Implementing and Simulating Dynamic Traffic Assignment with Intelligent Transportation Systems in Cube Avenue

As urban populations and traffic congestion levels increase, effective use of information and communication tools and intelligent transportation systems as becoming increasingly important in order to maximize the efficiency of transportation networks. The appropriate placement and employment of these tools within a network is critical to their effectiveness. This presentation proposes and demonstrates the use of a commercial transportation simulation tool to simulate dynamic traffic assignment and rerouting to model route modifications as a result of traffic information.

Foytik, Peter↗

DOE OSTI2021

Molecular dynamics simulations of radiation cascade evolution near cellular dislocation structures in additively manufactured stainless steels

The dislocation dense cellular structures observed in additively manufactured (AM) stainless steels can allow these materials to have enhanced strength and ductility compared to conventionally manufactured materials. Effective design of radiation resistant materials often necessitates creating microstructures containing of a high density of dislocations that act as sinks for radiation induced defects. This work uses molecular dynamics simulations to study the impact of pre-existing defects on radiation damage in stainless steel 316L fabricated by the laser powder bed fusion process. The evolution of synthetically generated dislocations in a crystal system in response to multiple radiation collision cascades suggest that dislocation dense regions reduce the probability of surviving point defects forming clusters. This finding seems to result from the dislocation cores’ ability to absorb point defects, thereby limiting the formation of clusters during the recombination phase. However, after successive cascades, the ability of the dislocation entanglement to prevent defect formation is diminished and defect cluster formation trends to an equilibrium. Accordingly, it is hypothesized that in an experimental setting the dislocation cells will initially act as neutral sinks for point defects, potentially delaying the onset of radiation damage effects; but, these structures will be degraded due to radiation-enhanced diffusion and lose their effectiveness at doses above 1.5–2 displacements per atom. Nevertheless, the ability to produce complex and feature-specific microstructures suggests promise for the design of radiation tolerant materials using AM methods.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Ab initio molecular dynamics simulations of short-range order in Zr 50 Cu 45 Al 5 and Cu 50 Zr 45 Al 5 metallic glasses

Comparative analysis between Zr-rich Zr 50 Cu 45 Al 5 and Cu-rich Cu 50 Zr 45 Al 5 metallic glasses (MGs) is extensively performed to locate the key structural motifs accounting for their difference of glass forming ability. Here we adopt ab initio molecular dynamics simulations to investigate the local atomic structures of Zr 50 Cu 45 Al 5 and Cu 50 Zr 45 Al 5 MGs. A high content of icosahedral-related (full and distorted) orders was found in both samples, while in the Zr-rich MG full icosahedrons < 0,0,12,0 > is dominant, and in the Cu-rich one the distorted icosahedral orders, especially < 0,2,8,2 > and < 0,2,8,1 >, are prominent. And the < 0,2,8,2 > polyhedra in Cu 50 Zr 45 Al 5 MG mainly originate from Al-centered clusters, while the < 0,0,12,0 > in Zr 50 Cu 45 Al 5 derives from both Cu-centered clusters and Al-centered clusters. These difference may be ascribed to the atomic size difference and chemical property between Cu and Zr atoms. Lastly, the relatively large size of Zr and large negative heat of mixing between Zr and Al atoms, enhancing the packing density and stability of metallic glass system, may be responsible for the higher glass forming ability of Zr 50 Cu 45 Al 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Fast neutrino cooling of nuclear pasta in neutron stars: Molecular dynamics simulations

The direct Urca process of rapid neutrino emission can occur in nonuniform nuclear pasta phases that are expected in the inner crusts of neutron stars. Here, the periodic potential for a nucleon in nuclear pasta allows momentum conservation to be satisfied for direct Urca reactions. We improve on earlier work by modeling a rich variety of pasta phases (gnocchi, waffle, lasagna, and anti-spaghetti) with large-scale molecular dynamics simulations. We find that the neutrino luminosity due to direct Urca reactions in nuclear pasta can be 3 to 4 orders of magnitude larger than that from the modified Urca process in the NS core. Thus neutrino radiation from pasta could dominate radiation from the core and this could significantly impact the cooling of neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Initial beam dynamics simulations of a high-average-current field-emission electron source in a superconducting radiofrequency gun

High-power electron beams are sought-after tools in support to a wide array of societal applications. This paper explores the production of high-power electron beams by combining a high-current field-emission electron source to a superconducting radio-frequency (SRF) cavity. We especially carry out beam-dynamics simulations that demonstrate the viability of the scheme to form ~ 300 kW average-power electron beam using a 1+1/2-cell SRF gun.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Molecular Modeling and Molecular Dynamics Simulation of a Packed and Intact Bacterial Microcompartment

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria which encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semipermeable shell proteins and are packed full of enzyme cargoes and metabolites as they fulfill their function. Coupling together recent SAXS and proteomics work, it is possible to develop molecular models for these microcompartments and interrogate enzyme and metabolite dynamics within. Our primary goal of this study is to quantify the permeability of metabolite glyceraldehyde-3-phosphate (G3P) and dihydroxyacetone phosphate (DHAP) across the BMC shell through classical molecular dynamics simulation. The Haliangium ochraceum model of BMC shell (PDB: 6MZX) was used to model an intact BMC of approximately 10 million atoms. Working at this scale presented its own challenges in managing large data sets, with multiple challenges and hardware advances discussed that facilitated this work. Over approximately 750 ns of aggregate simulation, we see multiple permeation events for these metabolites that were added at high concentration through the pores present within BMC shell tiles. When compared to independent permeability estimates for the same metabolites determined through replica exchange umbrella sampling simulations, the permeabilities varied by approximately 3 orders of magnitude. Regardless, the permeability coefficients for both G3P and DHAP are highly similar and very high, such that only very small concentration gradients can be maintained across the BMC shell between the cytosol and BMC interior. The large simulation systems also facilitated comparisons for molecular diffusivity in the crowded environment within the BMC shell. By our estimates, the viscosity within a packed BMC shell is at least 10-fold higher than it would be in neat solution and is the real driver for varying permeability estimates we obtained through simulation. These findings will be used as design inputs for future bioengineering efforts to make products from BMCs, highlighting how permeable BMC shells can be.

Diffusion↗

DOE OSTI2025

Molecular Origin of Nanoscale Anion Ordering of LiTFSI Electrolytes Revealed through SAXS/WAXS and Molecular Dynamics Simulations

Recent developments in “water-in-salt” electrolytes have precipitated a renewed effort to study imide-based electrolytes. While previous small-/wide-angle X-ray scattering (SAXS/WAXS) studies have attributed the emergence of a low-Q peak in the SAXS profile of aqueous LiTFSI electrolytes to nanometer-scale anion clustering, a molecular-level understanding of the root of these clusters remains unclear. Here, in this study, we combined molecular dynamics simulations and SAXS/WAXS to study the solvation structures of LiTFSI in acetonitrile, methanol, and water. We concluded that hydrogen bonding in water and MeOH stabilizes anion clusters, while nonpolar methyl groups on methanol and acetonitrile interrupt the nanoscale ordering of TFSI anions. This causes LiTFSI in water and MeOH electrolytes to exhibit two low-Q SAXS profile peaks while LiTFSI in acetonitrile exhibits only a single peak below Q = 1 Å –1 . These findings shed light on the underlying molecular origins of nanoscale anion clusters, which may help in the design of the next generation of electrolyte chemistries.

Trojanowski, Lucas [Univ. of Michigan, Ann Arbor, ↗

DOE OSTI2018

Reactive molecular dynamics simulations of an excess proton in polyethylene glycol-water solutions

We investigate the effects of both triethylene glycol (TEG) and polyethylene glycol (PEG) on the structure and transport of an excess proton (H + ) in aqueous solution. A new reactive force field (ReaxFF) parameterisation is carried out to better describe proton transfer of H + with PEG polymers in aqueous solution via training against density functional theory (DFT) calculations. Using this ReaxFF parameterisation, reactive molecular dynamics simulations are carried out for three different systems containing H + : bulk water, aqueous TEG, and aqueous PEG with a single, oriented polymer chain. The H + diffusivity is suppressed due to the presence of TEG and enhanced in the oriented chain PEG system, particularly along the axis parallel to the PEG chain, relative to bulk water. The enhancement in proton mobility along the chain axis is mainly due to the structural, rather than vehicular, component. From free energy calculations for H + transfer, the increase in proton mobility is not due to a reduction in activation energy. Finally, from pair distribution function analysis of the local structure for water and hydronium we hypothesise the increased mobility is due to an entropic effect, specifically a favourable hydrogen bond network for long range proton transport due to the oriented PEG chain.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Integrated atomistic chemical imaging and reactive force field molecular dynamic simulations on silicon oxidation

In this paper, we quantitatively investigate with atom probe tomography, the effect of temperature on the interfacial transition layer suboxide species due to the thermal oxidation of silicon. The chemistry at the interface was measured with atomic scale resolution, and the changes in chemistry and intermixing at the interface were identified on a nanometer scale. We find an increase of suboxide (SiOx) concentration relative to SiO{sub 2} and increased oxygen ingress with elevated temperatures. Our experimental findings are in agreement with reactive force field molecular dynamics simulations. This work demonstrates the direct comparison between atom probe derived chemical profiles and atomistic-scale simulations for transitional interfacial layer of suboxides as a function of temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗