Engineering Papers⌕ Search

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.

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

At least 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Molecular Dynamics Simulation of a Multi-Walled Carbon Nanotube Based Gear

We used molecular dynamics to investigate the properties of a multi-walled carbon nanotube based gear. Previous work computationally suggested that molecular gears fashioned from (14,0) single-walled carbon nanotubes operate well at 50-100 gigahertz. The gears were formed from nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. The gear in this study was based on the smallest multi-walled nanotube supported by some experimental evidence. Each gear was a (52,0) nanotube surrounding a (37,10) nanotube with approximate 20.4 and 16,8 A radii respectively. These sizes were chosen to be consistent with inter-tube spacing observed by and were slightly larger than graphite inter-layer spacings. The benzyne teeth were attached via 2+4 cycloaddition to exterior of the (52,0) tube. 2+4 bonds were used rather than the 2+2 bonds observed by Hoke since 2+4 bonds are preferred by naphthalene and quantum calculations by Jaffe suggest that 2+4 bonds are preferred on carbon nanotubes of sufficient diameter. One gear was 'powered' by forcing the atoms near the end of the outside buckytube to rotate to simulate a motor. A second gear was allowed to rotate by keeping the atoms near the end of its outside buckytube on a cylinder. The ends of both gears were constrained to stay in an approximately constant position relative to each other, simulating a casing, to insure that the gear teeth meshed. The stiff meshing aromatic gear teeth transferred angular momentum from the powered gear to the driven gear. The simulation was performed in a vacuum and with a software thermostat. Preliminary results suggest that the powered gear had trouble turning the driven gear without slip. The larger radius and greater mass of these gears relative to the (14,0) gears previously studied requires a smaller rotation rate and multiple rows of teeth to avoid excessive force on the gear teeth resulting, in slip and failure of the driven gear to turn. We hope that studies such as these will eventually lead to synthesis of components that can be assembled into atomically precise fullerene machines. These machines, in turn, may someday be used in machine-phase fullerene materials with remarkable properties.

Han, Jie↗
NASA NTRS1990

Staggered solution procedures for multibody dynamics simulation

A computational procedure is presented for the direct integration of the multibody dynamical (MBD) equations with constraints. Numerical experiments conducted with this MBD procedure indicate that it yields robust solutions in cases where the step size furnishes more than 20 samples for the period of highest apparent response frequency in a given multibody system. The MBD solution procedure is implemented in two separate modules: the generalized coordinated solver, 'CINT', and the constraint Lagrange multiplier solver, 'LINT'.

Park, K. C.↗
DOE OSTI2025

First Principles Molecular Dynamics Simulations of Ammonia Adsorption onto MFI Zeolite Nanosheets

MFI zeolite nanosheets membranes are promising candidates for ammonia separation from nitrogen and hydrogen, yet questions remain on the origin of their high selectivity. Silanols, Si-OH, are present in high concentration at the surface of zeolite nanosheets, and force-field-based simulations indicate that surface adsorption at the silanols contributes to selectivity. Silanols can chemically react with ammonia, which may further contribute to the ability of zeolite nanosheet membranes to separate it from other gases. In this work, we use first-principles molecular dynamics techniques to simulate ammonia’s behavior within stacked MFI zeolite nanosheets. We find that at 523 K and a loading corresponding to 35 bar, conditions desired for the ammonia separation, about 30% of ammonia reacts with surface silanols. Our work explores H-bonding and proton transfer within this system.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Orientational Pair Correlations and Local Structure of Benzonitrile from Molecular Dynamics Simulations with Comparisons to Experiments [plus supplemental information]

Here, we present an experimentally parametrized molecular dynamics study of single-molecule and collective orientational relaxation in neat benzonitrile through the analysis of the reorientational anisotropy and polarizability anisotropy time-correlation function (PA-TCF). The simulations show that the PA-TCF is dominated by collective reorientation after 20 ps. Collective reorientation is found to be slower than single-molecule reorientation by a factor of 1.67, consistent with recent experiments. The simulations provide direct evidence of local antiparallel benzonitrile configurations. These structures, which have been the center of some debate, are responsible for the slower rate of collective vs. single-molecule reorientation in the liquid. Further structural analysis indicates that significant Coulombic interactions between the nitrile group and hydrogen atoms on adjacent molecules play a role in the formation of the antiparallel structures. The single-molecule dynamics reflected in the anisotropy are complex and consist of a ballistic regime, restricted angular diffusion, and spatially anisotropic free diffusion. The principal components of the rotational diffusion tensor are independently obtained and shown to reproduce the free diffusion regime of the anisotropy for each principal axis according to the predictions of a previous theory.

14 SOLAR ENERGY↗
DOE OSTI2018

The Kinetics of the Liquid Phase Nucleation in a Stretched FCC Crystal: A Molecular Dynamics Simulation

The kinetics of a spontaneous formation of liquid phase in a stretched (superheated) Lennard-Jones fcc crystal is studied. Molecular dynamics experiments are used to determine the main parameters of the nucleation process: nucleation frequency J, diffusion coefficient of nuclei D{sub *}, nonequilibrium Zel’dovich factor Z, and critical nucleus size R{sub *}. The calculations are performed at negative pressures from the endpoint of the melting line and at positive pressures that are higher by a factor of eight than the critical pressure. The simulation results are compared to the classical homogeneous nucleation theory. It is found that the theory qualitatively correctly reproduces the dynamics of developing the process. The theory and the simulation demonstrate good quantitative agreement for the transition rate of the liquid phase nucleus through the critical size, but there is large difference in the numbers of critical nuclei in the unit volume of the metastable phase. In the case of significant superheatings and negative pressures, the contribution of the energy of elastic stresses to the moving force of the phase transformation is small and it can be neglected in a first approximation. The mismatch between the theory and the simulation results can be eliminated taking that the surface free energy of a curved “crystal–liquid droplet” interface is smaller than that of a plane interface by 30–35%.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

First-Principles Molecular Dynamics Simulations of Ammonia Adsorption onto MFI Zeolite Nanosheets

MFI zeolite nanosheet membranes are promising candidates for ammonia separation from nitrogen and hydrogen gases, yet questions remain on the origin of their high selectivity. Silanols, Si–OH, are present in high concentration at the surface of zeolite nanosheets, and force-field-based simulations indicate that surface adsorption at the silanols contributes to selectivity. Acidic silanol groups can chemically react with ammonia via transfer of a proton to form ammonium ions, which may further contribute to the ability of zeolite nanosheet membranes to separate ammonia from other gases. In this work, we used first-principles molecular dynamics techniques to simulate ammonia’s behavior within stacked MFI zeolite nanosheets. We found that at 523 K and a loading corresponding to 35 bar, conditions desired for the ammonia separation, about 30% of ammonia reacts with surface silanols. Our work explored H-bonding and proton transfer within this system.

ammonia↗
DOE OSTI2016

Computational Fluid Dynamics Simulation Study of Active Power Control in Wind Plants

This paper presents an analysis performed on a wind plant's ability to provide active power control services using a high-fidelity computational fluid dynamics-based wind plant simulator. This approach allows examination of the impact on wind turbine wake interactions within a wind plant on performance of the wind plant controller. The paper investigates several control methods for improving performance in waked conditions. One method uses wind plant wake controls, an active field of research in which wind turbine control systems are coordinated to account for their wakes, to improve the overall performance. Results demonstrate the challenge of providing active power control in waked conditions but also the potential methods for improving this performance.

17 WIND ENERGY↗
DOE OSTI2020

Nonadiabatic Dynamics Simulation of the Wavelength-Dependent Photochemistry of Azobenzene Excited to the nπ* and ππ* Excited States

Azobenzene is one of the most ubiquitous photoswitches in photochemistry and a prototypical model for photoisomerizing systems. Despite this, its wavelength-dependent photochemistry has puzzled researchers for decades. Upon excitation to the higher energy ππ* excited state instead of the dipole-forbidden nπ* state, the quantum yield of isomerization from trans- to cis-azobenzene is halved. The difficulties associated with unambiguously resolving this effect both experimentally and theoretically have contributed to lasting controversies regarding the photochemistry of azobenzene. Here, we systematically characterize the dynamic photoreaction pathways of azobenzene by performing first-principles simulations of the nonadiabatic dynamics following excitation to both the ππ* and the nπ* states. In this work, we demonstrate that ground-state recovery is mediated by two distinct S1 decay pathways: a reactive twisting pathway and an unreactive planar pathway. Increased preference for the unreactive pathway upon ππ* excitation largely accounts for the wavelength-dependent behavior observed in azobenzene.

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

Efficacy of the radial pair potential approximation for molecular dynamics simulations of dense plasmas

Macroscopic simulations of dense plasmas rely on detailed microscopic information that can be computationally expensive and is difficult to verify experimentally. In this work, we delineate the accuracy boundary between microscale simulation methods by comparing Kohn–Sham density functional theory molecular dynamics (KS-MD) and radial pair potential molecular dynamics (RPP-MD) for a range of elements, temperature, and density. By extracting the optimal RPP from KS-MD data using force matching, we constrain its functional form and dismiss classes of potentials that assume a constant power law for small interparticle distances. Our results show excellent agreement between RPP-MD and KS-MD for multiple metrics of accuracy at temperatures of only a few electron volts. The use of RPPs offers orders of magnitude decrease in computational cost and indicates that three-body potentials are not required beyond temperatures of a few eV. Due to its efficiency, the validated RPP-MD provides an avenue for reducing errors due to finite-size effects that can be on the order of ∼20%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Study of iron structure stability in high temperature molten lead-bismuth eutectic with oxygen injection using molecular dynamics simulation

Corrosion of structural materials in high temperature molten lead-bismuth eutectic is a major problem for design of PbBi cooled reactor. One technique to inhibit corrosion process is to inject oxygen into coolant. In this paper we study and focus on a way of inhibiting the corrosion of iron using molecular dynamics method. For the simulation results we concluded that effective corrosion inhibition of iron may be achieved by injection 0.0532 wt% to 0.1156 wt% oxygen into liquid lead-bismuth. At this oxygen concentration the structure of iron material will be maintained at about 70% in bcc crystal structure during interaction with liquid metal.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRSDate not supplied

Validation of Cryogenic Propellant Tank Self-Pressurization by Leveraging Reduced Order Modeling within Computational Fluid Dynamics Simulation

Validation of cryogenic propellant tank self-pressurization was performed using a hybrid Computational Fluid Dynamics (CFD) and reduced order modeling methodology. Data from a liquid hydrogen ground test conducted at the K-site facility at the National Aeronautics and Space Administration (NASA) Glenn Research Center was used for the validation effort. Liquid phase dynamics were explicitly resolved with a CFD tool. Vapor phase dynamics were modeled as a point mass that communicated heat from the tank wall to the liquid phase via a boundary condition used at the gas-liquid interface. The method proved to be more accurate, robust, and efficient than explicit resolution of the dynamics using a standard Volume of Fluid (VOF) methodology. The subject pressurization process was found to be heavily dependent upon both the relatively high liquid temperature gradient near the gas-liquid interface and the natural convection flow path. Modeling the gas-liquid interface as an immovable surface eliminated temperature gradient destroying gas-liquid interface velocities observed in VOF simulations, and correspondingly enabled more rapid simulation since interface advection was not allowed. The single phase computational domain also facilitated the ability to demonstrate spatial resolution convergence of natural convection cells within the liquid which significantly impacted the tank pressurization rate. This work was used to demonstrate the critical physics for tank self-pressurization and numerical methodologies that may be used to best resolve those physics. The findings informed development and operation of production level CFD tools used in the Fluid Dynamics Branch at NASA Marshall Space Flight Center.

J. M. Brodnick↗
DOE OSTI2021

AENET–LAMMPS and AENET–TINKER : Interfaces for accurate and efficient molecular dynamics simulations with machine learning potentials

Machine-learning potentials (MLPs) trained on data from quantum-mechanics based first-principles methods can approach the accuracy of the reference method at a fraction of the computational cost. To facilitate efficient MLP-based molecular dynamics and Monte Carlo simulations, an integration of the MLPs with sampling software is needed. Here, we develop two interfaces that link the atomic energy network (ænet) MLP package with the popular sampling packages TINKER and LAMMPS. The three packages, ænet, TINKER, and LAMMPS, are free and open-source software that enable, in combination, accurate simulations of large and complex systems with low computational cost that scales linearly with the number of atoms. Scaling tests show that the parallel efficiency of the ænet–TINKER interface is nearly optimal but is limited to shared-memory systems. The ænet–LAMMPS interface achieves excellent parallel efficiency on highly parallel distributed memory systems and benefits from the highly optimized neighbor list implemented in LAMMPS. We demonstrate the utility of the two MLP interfaces for two relevant example applications: the investigation of diffusion phenomena in liquid water and the equilibration of nanostructured amorphous battery materials.

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

Thermal Protection System Cavity Heating for Simplified and Actual Geometries Using Computational Fluid Dynamics Simulations with Unstructured Grids

Thermal Protection System (TPS) Cavity Heating is predicted using Computational Fluid Dynamics (CFD) on unstructured grids for both simplified cavities and actual cavity geometries. Validation was performed using comparisons to wind tunnel experimental results and CFD predictions using structured grids. Full-scale predictions were made for simplified and actual geometry configurations on the Space Shuttle Orbiter in a mission support timeframe.

McCloud, Peter L.↗
DOE OSTI2016

Nonintrusive Uncertainty Quantification of Computational Fluid Dynamics Simulations of a Bench-Scale Fluidized-Bed Gasifier

Uncertainty quantification (UQ) analysis is increasingly becoming one of the major requirements of simulation-based engineering to assess the confidence in the results and make better-informed decisions based on the insight derived from the simulations. In an earlier study, Bayesian UQ analysis was applied to existing bench-scale fluidized-bed gasifier experiment results. In the current study, a series of simulations were carried over with the open-source computational fluid dynamics software MFiX to reproduce the experimental conditions, where three operating factors, i.e., coal flow rate, coal particle diameter, and steam-to-oxygen ratio, were systematically varied to understand their effect on the syngas composition. Bayesian UQ analysis was this time performed on the numerical results for comparison purposes. This is part of ongoing research efforts to explore the applicability of advanced UQ methods and processes such as Bayesian methods for large-scale complex multiphase flow simulations. As part of Bayesian UQ analysis, a global sensitivity analysis was performed based on the simulation results, which shows that the predicted syngas composition is strongly affected not only by the steam-to-oxygen ratio (which was observed in experiments as well) but also by variation in the coal flow rate and particle diameter (which was not observed in experiments). The carbon monoxide mole fraction is underpredicted at lower steam-to-oxygen ratios and overpredicted at higher steam-to-oxygen ratios. The opposite trend is observed for the carbon dioxide mole fraction. These discrepancies are attributed to either excessive segregation of the phases that leads to the fuel-rich or -lean regions or alternatively the selection of reaction models, where different reaction models and kinetics can lead to different syngas compositions throughout the gasifier.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Micro-beam bending of FCC bicrystals: A comparison between defect dynamics simulations and experiments

To understand the role of the grain boundary (GB) in plasticity at small scale, a concurrently coupled mesoscale plasticity model was developed to simulate micro-bending of bicrystalline micron-sized beams. By coupling dislocation dynamics (DD) with a finite element model (FEM), a novel defect dynamics model provides the means to investigate intricate interactions between dislocations and GBs under various loading conditions. Our simulations of micro-bending agree well with corresponding micro-bending experiments, and they show that mechanical response of bicrystals could have not only hardening but also softening depending on the characters of the GB. In addition, changing the location of the GB in the microbeams results in different mechanical responses; GBs located at the neutral plane show softening compared to single crystals, while inclined GBs located halfway along the length of the beam show little effect. In conclusion, simulation results could provide a clear picture on detailed dislocation-GB interactions, and quantitative resolved shear stress analysis supplemented by dislocation density distribution is used to analyze the mechanical response of bicrystalline samples.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Effects of Frothers and Oil at Saltwater–Air Interfaces for Oil Separation: Molecular Dynamics Simulations and Experimental Measurements

Separating oil from saltwater is a process relevant to some industries and may be aided by bubble and froth generation. Simulating saltwater–air interfaces adsorbed with surfactants and oil molecules can assist in understanding froth stability to improve separation. Here, combining with surface tension experimental measurements, in this work we employ molecular dynamics with a united-atom force field to linear alkane oil and three surfactant frothers, methyl isobutyl carbinol (MIBC), terpineol, and ethyl glycol butyl ether (EGBE), to investigate their synergistic behaviors for oil separation. The interfacial phenomena were measured for a range of frother surface coverages on saltwater. Density profiles of the hydrophilic and hydrophobic portions of the frothers show an expected orientation of alcohol groups adsorbing to the polar water. A decrease in surface tension with increasing surface coverage of MIBC and terpineol was observed and reflected in experiments where the frother concentration increased. Relations between surface coverage and bulk concentration were observed by comparing the surface tension decreases. Additionally, a range of oil surface coverages was explored when the interface has a thin layer of adsorbed frother molecules. Finally, the obtained results indicate that an increase in surface coverage of oil molecules led to an increase in surface tension for all frother types and the pair correlation functions depicted MIBC and terpineol as having higher distributions with water at closer distances than with oil.

54 ENVIRONMENTAL SCIENCES↗