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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The relation between lignin sequence and its 3D structure

Background: Lignin, the second most abundant biopolymer on earth, plays a major structural role in plants, conferring mechanical strength and regulating water conduction. Understanding the three-dimensional structure of lignin is important for fundamental reasons as well as engineering plants towards lignin valorization. Lignin lacks a specific primary sequence, making its average chemical composition the focus of most recent studies. However, it remains unclear whether the 3D structure of lignin molecules depends on their sequence. Methods: In this work, we performed all-atom molecular dynamics simulation of three S/G-lignin molecules with the same average composition but different sequence. Results: A detailed statistical analysis of the radius of gyration and relative shape anisotropy reveals that the lignin sequence has no statistically significant effect on the global three-dimensional structure. We found however, that homopolymers of C-lignin with the same molecular weight have smaller radii of gyration than S/G-lignin. We attribute this to lower hydroxyl content of C-lignin, which makes it more compact and rigid. Conclusions: The 3D structure of lignin is influenced by the overall content of monomeric units and interunit linkages and not by its precise primary sequence. General Significance: Lignin is assumed to not have a well-defined primary structure. The results presented here demonstrate there are no significant differences in the global 3D structure of lignin molecules with the same average composition but different primary sequence.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A High-Speed Rotational Diamond Anvil Cell for In Situ Analysis of Hierarchical Microstructural Evolution of Metallic Alloys during Extreme Shear Deformation

High speed shear deformation is ubiquitous in engineering applications, ranging from material processing methods such as friction stir processing/extrusion and in tribological contacts. However, analyzing the microstructural evolution of materials while they are undergoing high speed shear deformation have been a long-standing challenge. This led to predominant reliance on ex situ microscopy before and after shear deformation. But ex situ microscopy lacks the ability to analyze dynamic and transient hierarchical microstructural evolution mechanisms that could occur during shear deformation of materials. Therefore, to better understand the dynamic mechanisms of mass and energy transfer in materials under shear deformation, we developed a first of its kind high-speed rotational diamond anvil cell (HS-RDAC) for synchrotron-based in situ high-energy x-ray diffraction (XRD). We studied the time resolved lattice strain evolution, XRD peak broadening and changes in spatial variation of shear deformation induced alloying in pure metal and metal alloy sheets and powder mixture using the HS-RDAC. These in situ results were combined with detailed ex situ microstructural characterization before and after the shear deformation using transmission electron microscopy and atom probe tomography, which revealed the different stages of evolution of a shear deformation induced hierarchical nanostructure. Multiscale computational simulations including computational fluid dynamics, crystal plasticity, molecular dynamic simulation and density functional theory uncovered the mechanisms behind morphological changes, evolution of defect structures and changes in driving force for shear deformation induced intermixing. In conclusion, this in situ HS-RDAC capability, in combination with ex situ microstructural characterization and computational simulations, can provide new insights into the hierarchical microstructural evolution pathway during shear deformation.

36 MATERIALS SCIENCE↗
NASA NTRS2011

Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 are candidate materials for applications in extreme environments because of their high melting point, good mechanical properties and reasonable oxidation resistance. Unlike many ceramics, these materials have high thermal conductivity which can be advantageous, for example, to reduce thermal shock. Recently, we developed Tersoff style interatomic potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current. Results at room temperature and at elevated temperatures will be reported.

Lawson, JOhn W.↗
NASA NTRS2012

Lattice Thermal Conductivity from Atomistic Simulations: ZrB2 and HfB2

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.

Lawson, John W.↗
NASA NTRS2012

Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.

Lawson, John W.↗
DOE OSTI2020

Comparison of fixed charge and polarizable models for predicting the structural, thermodynamic, and transport properties of molten alkali chlorides

Results from extensive molecular dynamics simulations of molten LiCl, NaCl, KCl, and RbCl over a wide range of temperatures are reported. Comparison is made between the “Polarizable Ion Model” (PIM) and the non-polarizable “Rigid Ion Model” (RIM). Densities, self-diffusivities, shear viscosities, ionic conductivities, and thermal conductivities are computed and compared with experimental data. In addition, radial distribution functions are computed from ab initio molecular dynamics simulations and compared with the two sets of classical simulations as well as experimental data. Here, the two classical models perform reasonably well at capturing structural and dynamic properties of the four molten alkali chlorides, both qualitatively and often quantitatively. With the singular exception of liquid density, for which the PIM is more accurate than the RIM, there are few clear trends to suggest that one model is more accurate than the other for the four alkali halide systems studied here.

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

New insights into the mechanism governing the elasticity of calcium silicate hydrate gels exposed to high temperature: A molecular dynamics study

When exposed to fire, the integrity of cement-based materials is governed by thermally-induced changes in the mechanical properties of their binding phase, i.e., the calcium–silicate–hydrate (C–S–H) gel. However, the effect of temperature on the structure, density, and mechanical properties of C–S–H remains only partially known. Here, based on reactive molecular dynamics simulations, we reveal the nature of thermally-induced damage in C–S–H gels. In general, we show that, at the atomic scale, exposure to high temperature results in partial dehydration, volumetric shrinkage, disordering, and stiffening in the C–S–H grains. However, we show that the thermal response of C–S–H strongly depends on its chemical composition, wherein C–S–H systems associated with lower Ca/Si molar ratios are able to undergo higher temperatures before amorphization. Based on these results, we demonstrate that the stiffness of C–S–H gels (i.e., including porosity—as probed by nanoindentation) is governed by a competition between the stiffening of the grains and the decrease in packing density—wherein the latter eventually become predominant.

36 MATERIALS SCIENCE↗
NASA NTRS2000

Plasticity and Kinky Chemistry of Carbon Nanotubes

Since their discovery in 1991, carbon nanotubes have been the subject of intense research interest based on early predictions of their unique mechanical, electronic, and chemical properties. Materials with the predicted unique properties of carbon nanotubes are of great interest for use in future generations of aerospace vehicles. For their structural properties, carbon nanotubes could be used as reinforcing fibers in ultralight multifunctional composites. For their electronic properties, carbon nanotubes offer the potential of very high-speed, low-power computing elements, high-density data storage, and unique sensors. In a continuing effort to model and predict the properties of carbon nanotubes, Ames accomplished three significant results during FY99. First, accurate values of the nanomechanics and plasticity of carbon nanotubes based on quantum molecular dynamics simulations were computed. Second, the concept of mechanical deformation catalyzed-kinky-chemistry as a means to control local chemistry of nanotubes was discovered. Third, the ease of nano-indentation of silicon surfaces with carbon nanotubes was established. The elastic response and plastic failure mechanisms of single-wall nanotubes were investigated by means of quantum molecular dynamics simulations.

Srivastava, Deepak↗
DOE OSTI2020

Molecular dynamics and continuum analyses of the electrokinetic zeta potential in nanostructured slit channels

This work presents a theoretical and numerical study of electrokinetic flow and the zeta potential for the case of slit channels with nanoscale surface features of dimensions comparable to the Debye length, by employing molecular dynamics simulations and continuum-level analyses. Additionally, a simple analytical model for considering the average effect of such surface features producing nanoscale roughness is proposed by employing matched asymptotic solutions for the charge density and fluid flow field, and matching conditions that satisfy electroneutrality and the Onsager reciprocal relation between the electroosmotic flow rate and streaming current. The proposed analytical model quantitatively accounts for results from molecular dynamics simulations that consider the presence of ion solvation shells and surface hydration layers. Our analysis indicates that a simultaneous knowledge of the electroosmotic and pressure-driven flow rate or streaming current can be instrumental to unambiguously determine the zeta potential in the presence of nanoscale surface features.

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

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid–liquid mixtures

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

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

Synergetic Effects of Soil Organic Matter Components During Interactions with Minerals

Mineral-associated soil organic matter (SOM) is critical for stabilizing organic carbon and mitigating climate change. However, mineral-SOM interactions at the molecular scale, particularly synergetic adsorption through organic-organic interaction on the mineral surface known as organic multilayering, remain poorly understood. This study investigates the impact of organic multilayering on mineral-SOM interactions, by integrating macroscale experiments and molecular-scale simulations that assess the individual and sequential adsorption of major SOM compounds – lauric acid (lipid), pentaglycine (amino acid), trehalose (carbohydrate), and lignin onto soil minerals. Ferrihydrite, Al-hydroxide, and calcite are exposed to SOM compounds to determine adsorption affinities and binding energies. Results show that lauric acid has 20-40 times higher K d than pentaglycine, following the order K d (ferrihydrite) > K d (Al-hydroxide) >> K d (calcite). Molecular-scale simulations confirm that lauric acid has a higher binding energy (30.8 kcal/mol) on ferrihydrite than pentaglycine (6.0 kcal/mol), attributed to lipid hydrophobicity. The lower binding energy of pentaglycine results from its hydrophilic amide groups, facilitating partitioning into water. Sequential experiments examine how the first layer of lipid or amino acid affects the adsorption of carbohydrate/lignin, which show little or no individual adsorption affinities. Macroscale results reveal that lipid and amino acid adsorption induce ferrihydrite particle repulsion increasing reactive surface area and enhancing carbohydrate/lignin adsorption independently and synergistically through organic multilayering. Molecular-scale results reveal that amino acid adsorbed on ferrihydrite interacts more readily with lignin macroaggregates (preformed in solution) than with individual lignin units, indicating organic multilayering via H-bonding. Further, these findings reveal the molecular mechanisms of SOM-mineral interactions, crucial for enhancing soil carbon stabilization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Molecular dynamics analysis of particle number fluctuations in the mixed phase of a first-order phase transition

Molecular dynamics simulations are performed for a finite nonrelativistic system of particles with Lennard-Jones potential. We study the effect of liquid-gas mixed phase on particle number fluctuations in coordinate subspace. A metastable region of the mixed phase, the so-called nucleation region, is analyzed in terms of a noninteracting cluster model. Large fluctuations due to spinodal decomposition are observed. They arise due to the interplay between the size of the acceptance region and that of the liquid phase. These effects are studied with a simple geometric model. The model results for the scaled variance of particle number distribution are compared with those obtained from the direct molecular dynamic simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Stockmayer fluid simulations for viscosity and glass transition temperature of ionic liquids

We develop a Stockmayer fluid model for molecular dynamics simulations of ionic liquids that captures molecular polarization, ionic conductivity, viscosity, and glass transition temperature, using ethylammonium nitrate (EAN) as an example. The ions in EAN are treated as spheres interacting via the Lennard-Jones potential with an embedded point charge and a permanent dipole moment. We show that our simulation results for EAN are consistent with experimental data and then explore the effects of the molecular parameters on the viscosity of ionic liquids. Our results indicate that viscosity monotonically increases with ionic charge and dipole moment but non-monotonically changes with ionic diameter (or molar volume). This non-monotonic trend arises from the competition among the electrostatic interactions, molecular packing, and size asymmetry between the cation and anion. In conclusion, our model also shows that long-lived ion pairs result in higher viscosities.

Coarse-grained simulations↗
DOE OSTI2026

Towards the Ultimate Strength of Medium‐Entropy Alloys Through Pulsed Lasers

The tensile strength of metals at extreme strain rates is a key predictor of their performance in ballistic and structural impact applications. An important experimental method to reach these extreme strain rates is the use of high-amplitude, short-duration pulsed lasers. The Jupiter Laser Facility at the Lawrence Livermore National Laboratory enabled probing for the first time the mechanical response of several promising High Entropy Alloys at times on the order of nanoseconds (strain rates of ∼10 7 and ∼10 9 s −1 ). The measured strength is in the range of 6 to 10 GPa, ten times the quasistatic value. The mechanisms of plastic deformation and failure were identified and quantified through analysis and molecular dynamics simulation. The reflected wave amplitudes, obtained by VISAR, were used to determine the tensile (spall) stress. The high tensile strength obtained is due to two factors: the strain-rate dependence of plastic flow and the kinetics of void nucleation, growth, and coalescence. The experimental results are compared with an analytical prediction considering both grain-interior and grain-boundary void initiation. Molecular dynamics simulations, conducted at strain rates of 10 8 and 10 9 s −1 , rationalize the experimental results. In conclusion, they provide valuable information about the process of failure evolution, and reveal that grain boundary separation plays a pivotal role in spalling.

Materials science↗
DOE OSTI2024

Solvation Structure of Methanol-in-Salt Electrolyte Revealed by Small-Angle X-ray Scattering and Simulations

The solvation structure of water-in-salt electrolytes was thoroughly studied, and two competing structures—anion solvated structure and anion network—were well-defined in recent publications. To further reveal the solvation structure in those highly concentrated electrolytes, particularly the influence of solvent, methanol was chosen as the solvent for this proposed study. In this work, small-angle X-ray scattering, small-angle neutron scattering, Fourier-transform infrared spectroscopy, and Raman spectroscopy were utilized to obtain the global and local structural information. With the concentration increment, the anion network formed by TFSI - became the dominant structure. Meanwhile, the hydrogen bonds among methanol were interrupted by the TFSI - anion and formed a new connection with them. Further, molecular dynamic simulations with two different force fields (GAFF and OPLS-AA) are tested, and GAFF agreed with synchrotron small-angle X-ray scattering/wide-angle X-ray scattering (SAXS/WAXS) results well and provided insightful information about molecular/ion scale solvation structure. This article not only deepens the understanding of the solvation structure in highly concentrated solutions, but more importantly, it provides additional strong evidence for utilizing SAXS/WAXS to validate molecular dynamics simulations.

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

Recent developments in 3D-printed membranes for water desalination

The recognition of membrane separations as a vital technology platform for enhancing the efficiency of separation processes has been steadily increasing. Concurrently, 3D printing has emerged as an innovative approach to fabricating reverse osmosis membranes for water desalination and treatment purposes. This method provides a high degree of control over membrane chemistry and structural properties. In particular, when compared to traditional manufacturing techniques, 3D printing holds the potential to expedite customization, a feat that is typically achieved through conventional manufacturing methods but often involves numerous processes and significant costs. This review aims to present the current advancements in membrane manufacturing technology specifically tailored for water desalination purposes, with a particular focus on the development of 3D-printed membranes. A comprehensive analysis of recent progress in 3D-printed membranes is provided. However, conducting experimental work to investigate various influential factors while ensuring consistent results poses a significant challenge. To address this, we explore how membrane manufacturing processes and performance can be effectively pre-designed and guided through the use of molecular dynamics simulations. Finally, this review outlines the challenges faced and presents future perspectives to shed light on research directions for optimizing membrane manufacturing processes and achieving optimal membrane performance.

3D printing membranes↗
NASA NTRS2006

MODELING FUNCTIONALLY GRADED INTERPHASE REGIONS IN CARBON NANOTUBE REINFORCED COMPOSITES

A combination of micromechanics methods and molecular dynamics simulations are used to obtain the effective properties of the carbon nanotube reinforced composites with functionally graded interphase regions. The multilayer composite cylinders method accounts for the effects of non-perfect load transfer in carbon nanotube reinforced polymer matrix composites using a piecewise functionally graded interphase. The functional form of the properties in the interphase region, as well as the interphase thickness, is derived from molecular dynamics simulations of carbon nanotubes in a polymer matrix. Results indicate that the functional form of the interphase can have a significant effect on all the effective elastic constants except for the effective axial modulus for which no noticeable effects are evident.

Seidel, G. D.↗