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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Ferroelectric molecular field-switch based on double proton transfer process: Static and dynamical simulations

In this work, we present a reversible ferroelectric molecular switch controlled by an external electric field. The studied (2Z)-1-(6-((Z)-2-hydroxy-2-phenylvinyl)pyridin-3-yl)-2-(pyridin-2(1H) -ylidene)ethanone (DSA) molecule is polarized by two uniaxial intramolecular hydrogen bonds. Two protons can be transferred along hydrogen bonds upon an electric field applied along the main molecular axis. The process results in reversion of the dipole moment of the system. Static ab initio and on-the-fly dynamical simulations of the DSA molecule placed in an external electric field give insight into the mechanism of the double proton transfer (DPT) in the system and allow for estimation of the time scale of this process. The results indicate that with increasing strength of the electric field, the step-wise mechanism of DPT changes into the downhill barrierless process in which the synchronous and asynchronous DPTs compete with each other.

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

Beam dynamics simulations of post low energy beam transport section in RAON heavy ion accelerator

RAON (Rare isotope Accelerator Of Newness) heavy ion accelerator of the rare isotope science project in Daejeon, Korea, has been designed to accelerate multiple-charge-state beams to be used for various science programs. In the RAON accelerator, the rare isotope beams which are generated by an isotope separation on-line system with a wide range of nuclei and charges will be transported through the post Low Energy Beam Transport (LEBT) section to the Radio Frequency Quadrupole (RFQ). In order to transport many kinds of rare isotope beams stably to the RFQ, the post LEBT should be devised to satisfy the requirement of the RFQ at the end of post LEBT, simultaneously with the twiss parameters small. We will present the recent lattice design of the post LEBT in the RAON accelerator and the results of the beam dynamics simulations from it. In addition, the error analysis and correction in the post LEBT will be also described.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Homogeneous states in driven granular mixtures: Enskog kinetic theory versus molecular dynamics simulations

The homogeneous state of a binary mixture of smooth inelastic hard disks or spheres is analyzed. The mixture is driven by a thermostat composed by two terms: a stochastic force and a drag force proportional to the particle velocity. The combined action of both forces attempts to model the interaction of the mixture with a bath or surrounding fluid. The problem is studied by means of two independent and complementary routes. First, the Enskog kinetic equation with a Fokker-Planck term describing interactions of particles with thermostat is derived. Then, a scaling solution to the Enskog kinetic equation is proposed where the dependence of the scaled distributions φ{sub i} of each species on the granular temperature occurs not only through the dimensionless velocity c = v/v{sub 0} (v{sub 0} being the thermal velocity) but also through the dimensionless driving force parameters. Approximate forms for φ{sub i} are constructed by considering the leading order in a Sonine polynomial expansion. The ratio of kinetic temperatures T{sub 1}/T{sub 2} and the fourth-degree velocity moments λ{sub 1} and λ{sub 2} (which measure non-Gaussian properties of φ{sub 1} and φ{sub 2}, respectively) are explicitly determined as a function of the mass ratio, size ratio, composition, density, and coefficients of restitution. Second, to assess the reliability of the theoretical results, molecular dynamics simulations of a binary granular mixture of spheres are performed for two values of the coefficient of restitution (α = 0.9 and 0.8) and three different solid volume fractions (ϕ = 0.00785, 0.1, and 0.2). Comparison between kinetic theory and computer simulations for the temperature ratio shows excellent agreement, even for moderate densities and strong dissipation. In the case of the cumulants λ{sub 1} and λ{sub 2}, good agreement is found for the lower densities although significant discrepancies between theory and simulation are observed with increasing density.

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

Investigating Alkylsilane Monolayer Tribology at a Single-Asperity Contact with Molecular Dynamics Simulation

Chemisorbed monolayer films are known to possess favorable characteristics for nanoscale lubrication of micro- and nanoelectromechanical systems (MEMS/NEMS). Prior studies have shown that the friction observed for monolayer-coated surfaces features a strong dependence on the geometry of contact. Specifically, tip-like geometries have been shown to penetrate into monolayer films, inducing defects in the monolayer chains and leading to plowing mechanisms during shear, which result in higher coefficients of friction (COF) than those observed for planar geometries. In this work, we use molecular dynamics simulations to examine the tribology of model silica single-asperity contacts under shear with monolayer-coated substrates featuring various film densities. It is observed that lower monolayer densities lead to reduced COFs, in contrast to results for planar systems where COF is found to be nearly independent of monolayer density. This is attributed to a liquid-like response to shear, whereby fewer defects are imparted in monolayer chains from the asperity, and chains are easily displaced by the tip as a result of the higher free volume. This transition in the mechanism of molecular plowing suggests that liquid-like films should provide favorable lubrication at single-asperity contacts.

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

Method of Obtaining High Resolution Intrinsic Wire Boom Damping Parameters for Multi-Body Dynamics Simulations

The goal of NASA's Magnetospheric MultiScale (MMS) mission is to understand magnetic reconnection with sensor measurements from four spinning satellites flown in a tight tetrahedron formation. Four of the six electric field sensors on each satellite are located at the end of 60- meter wire booms to increase measurement sensitivity in the spin plane and to minimize motion coupling from perturbations on the main body. A propulsion burn however, might induce boom oscillations that could impact science measurements if oscillations do not damp to values on the order of 0.1 degree in a timely fashion. Large damping time constants could also adversely affect flight dynamics and attitude control performance. In this paper, we will discuss the implementation of a high resolution method for calculating the boom's intrinsic damping, which was used in multi-body dynamics simulations. In summary, experimental data was obtained with a scaled-down boom, which was suspended as a pendulum in vacuum. Optical techniques were designed to accurately measure the natural decay of angular position and subsequently, data processing algorithms resulted in excellent spatial and temporal resolutions. This method was repeated in a parametric study for various lengths, root tensions and vacuum levels. For all data sets, regression models for damping were applied, including: nonlinear viscous, frequency-independent hysteretic, coulomb and some combination of them. Our data analysis and dynamics models have shown that the intrinsic damping for the baseline boom is insufficient, thereby forcing project management to explore mitigation strategies.

Yew, Alvin G.↗
DOE OSTI2021

Cavity molecular dynamics simulations of vibrational polariton-enhanced molecular nonlinear absorption

We report recent experiments have observed that the chemical and photophysical properties of molecules can be modified inside an optical Fabry– Pérot microcavity under collective vibrational strong coupling (VSC) conditions, and such modification is currently not well understood by theory. In an effort to understand the origin of such cavity-induced phenomena, some recent studies have focused on the effect of the cavity environment on the nonlinear optical response of the molecular subsystem. Here, we use a recently proposed protocol for classical cavity molecular dynamics simulations to numerically investigate the linear and the nonlinear response of liquid carbon dioxide under such VSC conditions following an optical pulse excitation. We find that applying a strong pulse of excitation to the lower hybrid light–matter state, i.e., the lower polariton (LP), can lead to an overall molecular nonlinear absorption that is enhanced by up to two orders of magnitude relative to the excitation outside the cavity. This polariton-enhanced multiphoton absorption also causes an ultrashort LP lifetime (0.2 ps) under strong illumination. Unlike usual polariton relaxation processes—whereby polaritonic energy transfers directly to the manifold of singly excited vibrational dark states—under the present mechanism, the LP transfers energy directly to the manifold of higher vibrationally excited dark states; these highly excited dark states subsequently relax to the manifold of singly excited states with a lifetime of tens of ps. Because the present mechanism is generic in nature, we expect these numerical predictions to be experimentally observed in different molecular systems and in cavities with different volumes.

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

Conformational and oligomeric states of SPOP from small-angle X-ray scattering and molecular dynamics simulations

Speckle-type POZ protein (SPOP) is a substrate adaptor in the ubiquitin proteasome system, and plays important roles in cell-cycle control, development, and cancer pathogenesis. SPOP forms linear higher-order oligomers following an isodesmic self-association model. Oligomerization is essential for SPOP’s multivalent interactions with substrates, which facilitate phase separation and localization to biomolecular condensates. Structural characterization of SPOP in its oligomeric state and in solution is, however, challenging due to the inherent conformational and compositional heterogeneity of the oligomeric species. Here, we develop an approach to simultaneously and self-consistently characterize the conformational ensemble and the distribution of oligomeric states of SPOP by combining small-angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations. We build initial conformational ensembles of SPOP oligomers using coarse-grained molecular dynamics simulations, and use a Bayesian/maximum entropy approach to refine the ensembles, along with the distribution of oligomeric states, against a concentration series of SAXS experiments. Our results suggest that SPOP oligomers behave as rigid, helical structures in solution, and that a flexible linker region allows SPOP’s substrate-binding domains to extend away from the core of the oligomers. Additionally, our results are in good agreement with previous characterization of the isodesmic self-association of SPOP. In the future, the approach presented here can be extended to other systems to simultaneously characterize structural heterogeneity and self-assembly.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

On the equivalence of nonequilibrium and equilibrium measurements of slip in molecular dynamics simulations

Here, we show that nonequilibrium and equilibrium measurements of slip are consistent, provided the hydrodynamic wall location associated with the equilibrium measurement is properly taken into account. The latter is a strong function of the fluid state and wall-fluid interaction and cannot be neglected as it typically has been. Our results are based on an alternative approach for calculating the hydrodynamic wall location which alleviates most of the difficulties associated with its calculation via a Green-Kubo integral that appear to have contributed to its neglect. Extensive molecular dynamics simulations are used to validate our approach including a model for calculating the slip length that does not involve a Green-Kubo integral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Comparing theoretical predictions of radiation-free velocities of edge dislocations to molecular dynamics simulations

Transonic defect motion is of interest for high strain-rate plastic deformation as well as for crack propagation. Ever since Eshelby's 1949 prediction [J. D. Eshelby, Proc. Phys. Soc. A 62, 307 (1949)] in the isotropic limit of a “radiation-free” transonic velocity $v$ RF = $\sqrt{2}$c T , where shock waves are absent, there has been speculation about the significance of radiation-free velocities (if they truly exist) for defect mobility. Here, in this work, we argue that they do not play any significant role in dislocation dynamics in metals, based on comparing theoretical predictions of radiation-free velocities for transonic edge dislocations with molecular dynamics simulations for two face-centered cubic metals: Ag, where theory predicts radiation-free states, and Cu, where it does not.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Thermal transport in semicrystalline polyethylene by molecular dynamics simulation

Recent research has highlighted the potential to achieve high-thermal-conductivity polymers by aligning their molecular chains. Combined with other merits, such as low-cost, corrosion resistance, and light weight, such polymers are attractive for heat transfer applications. Due to their quasi-one-dimensional structural nature, the understanding on the thermal transport in those ultra-drawn semicrystalline polymer fibers or films is still lacking. Here, we built the ideal repeating units of semicrystalline polyethylene and studied their dependence of thermal conductivity on different crystallinity and interlamellar topology using the molecular dynamics simulations. We found that the conventional models, such as the Choy-Young's model, the series model, and Takayanagi's model, cannot accurately predict the thermal conductivity of the quasi-one-dimensional semicrystalline polyethylene. A modified Takayanagi's model was proposed to explain the dependence of thermal conductivity on the bridge number at intermediate and high crystallinity. We also analyzed the heat transfer pathways and demonstrated the substantial role of interlamellar bridges in the thermal transport in the semicrystalline polyethylene. Lastly, our work could contribute to the understanding of the structure–property relationship in semicrystalline polymers and shed some light on the development of plastic heat sinks and thermal management in flexible electronics.

36 MATERIALS SCIENCE↗
NASA NTRS1993

A fast recursive algorithm for molecular dynamics simulation

The present recursive algorithm for solving molecular systems' dynamical equations of motion employs internal variable models that reduce such simulations' computation time by an order of magnitude, relative to Cartesian models. Extensive use is made of spatial operator methods recently developed for analysis and simulation of the dynamics of multibody systems. A factor-of-450 speedup over the conventional O(N-cubed) algorithm is demonstrated for the case of a polypeptide molecule with 400 residues.

Jain, A.↗
DOE OSTI2024

Effect of Zwitterionic Additives on Solvation and Transport of Sodium and Potassium Cations in (Ethylene Oxide)10: A Molecular Dynamics Simulation Study

Sodium- (Na + ) and potassium- (K + ) ion batteries are cost-effective alternatives to lithium-ion (Li + ) batteries due to the abundant sodium and potassium resources. Solid polymer electrolytes (SPEs) are essential for safer and more efficient Na + and K + batteries because they often exhibit low ionic conductivity at room temperature. While zwitterionic (ZW) materials enhance Li+ battery conductivity, their potential for Na + and K + transport in batteries remains unexplored. In this study, we investigated the effect of three ZW molecules (ChoPO4, i.e., 2-methacryloyloxyethyl phosphorylcholine, ImSO3, i.e., sulfobetaine ethylimidazole, and ImCO2, i.e., carboxybetaine ethylimidazole) on the dissociation of Na + and K + coordination with ethylene oxide (EO) chains in EO-based electrolytes through molecular dynamics simulations. Our results showed that ChoPO4 possessed the highest cation–EO10 dissociation ability, while ImSO3 exhibited the lowest. Such dissociation ability correlated with the cation–ZW molecule coordination strength: ChoPO4 and ImSO3 showed the strongest and the weakest coordination with cations. However, the cation–ZW molecule coordination could slow the cationic diffusion. The competition of these effects resulted in accelerating or decelerating cationic diffusion. Our simulated results showed that ImCO2 enhanced Na + diffusion by 20%, while ChoPO4 and ImSO3 led to a 10% reduction. For K + , ChoPO4 reduced its diffusion by 40%, while ImCO2 and ImSO3 caused a similar decrease of 15%. These findings suggest that the ZW structure and the cationic size play an important role in the ionic dissociation effect of ZW materials.

25 ENERGY STORAGE↗
DOE OSTI2015

Frequency-dependent local field factors in dielectric liquids by a polarizable force field and molecular dynamics simulations

A force field model for calculating local field factors, i.e. the linear response of the local electric field for example at a nucleus in a molecule with respect to an applied electric field, is discussed. It is based on a combined charge-transfer and point-dipole interaction model for the polarizability, and thereby it includes two physically distinct terms for describing electronic polarization: changes in atomic charges arising from transfer of charge between the atoms and atomic induced dipole moments. A time dependence is included both for the atomic charges and the atomic dipole moments and if they are assumed to oscillate with the same frequency as the applied electric field, a model for frequency-dependent properties are obtained. Furthermore, if a life-time of excited states are included, a model for the complex frequency-dependent polariability is obtained including also information about excited states and the absorption spectrum. We thus present a model for the frequency-dependent local field factors through the first molecular excitation energy. It is combined with molecular dynamics simulations of liquids where a large set of configurations are sampled and for which local field factors are calculated. We are normally not interested in the average of the local field factor but rather in configurations where it is as high as possible. In electrical insulation, we would like to avoid high local field factors to reduce the risk for electrical breakdown, whereas for example in surface-enhanced Raman spectroscopy, high local field factors are desired to give dramatically increased intensities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Advanced analysis of laser-driven pulsed magnetic diffusion based on quantum molecular dynamics simulation

Magnetic diffusion plays an important role in inertial confinement fusion with strong magnetic fields. In this paper, we improve a previous analysis of the generation and diffusion of the magnetic field [Morita et al., Phys. Plasmas 25, 094505 (2018)]. For the generation process, we calculate the temporal evolution of the coil current using a self-consistent circuit model. The results show that the peak of the calculated magnetic field is delayed by 1.2 ns compared with that of the incident laser pulse. For the diffusion process, we evaluate the electrical conductivity of warm dense gold over a wide temperature range (300 K–100 eV) by combining the Kubo–Greenwood formula based on a quantum molecular dynamics simulation with the modified Spitzer model. Our simulation shows that the maximum magnetic field (530 T) that penetrates the cone is delayed by 2.5 ns compared with the laser peak. This result is consistent with experiments [Sakata et al., Nat. Commun. 9, 3937 (2018)] that showed that applying a strong magnetic field improved the heating efficiency of fusion fuel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Nucleation of Salt Crystals in Clay Minerals: Molecular Dynamics Simulation

Nucleation of salt crystals in confined media occurs in many processes of high importance, such as injection of CO 2 in geological formations for its sequestration. Specifically, salt precipitation in clays, a main component of sedimentary rock, is an important phenomenon. The crystals precipitate on the pores’ surface, modify the pore space morphology, and reduce its flow and transport properties. Despite numerous efforts to understand the mechanisms of nucleation of salt crystals in confined media, the effect of the clay’s chemistry on the growth, distribution, and properties of the crystals is not well understood. We discuss the results of extensive molecular dynamics simulation of nucleation and growth of NaCl crystals in a clay pore using molecular models of two types of clay minerals, Na-montmorillonite and kaolinite. Clear evidence is introduced for the nucleation of the salt crystals that indicates that the molecular structure of clay minerals affects their spatial distribution, although the nucleation mechanism is the same in both types of clays.

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

Effect of different manufacturing methods on polyamide reverse-osmosis membranes for desalination: Insights from molecular dynamics simulations

Membranes are a key technology platform for a broad application of energy-efficient separations. To best serve the separation demands of industry, the manufacturing processes for these membranes are garnering increasing attention. In particular, for water desalination, the industry leading polyamide (PA) reverse osmosis (RO) membranes can be manufactured via molecular layer-by-layer (mLBL) deposition, interfacial polymerization (IP), and 3D-printing technique. However, the influence of different manufacturing methods on PA membrane’s properties is far from understood. Here, in this study, we present the high-pressure transport behavior of water and salt ions for PA membranes formed with IP, mLBL, and 3D-printing through non-equilibrium molecular dynamics simulations. Studies show that membranes fabricated with 3D-printing have similar performances to those manufactured using mLBL, quantified by water permeability, rejection of salt ions, structural integrity, and porosity features. However, the membranes formed with IP exhibit faster water transport, lower rejection, worse structural integrity, and more inhomogeneous network pores than those constructed using mLBL and 3D-printing. The unconnected water-accessible space governs water transport for PA membranes formed with mLBL and 3D-printing, which offers the impermanent open-closed pores that enable water to jump through PA membranes. In contrast, the permeated water-enterable space plays a prominent role in water movement across the PA membrane formed with IP, providing a continuous transport channel at high pressure. Importantly, we observe the more significant compaction features at high pressure for PA membranes formed with IP than mLBL and 3D printing. In short, these findings provide a comprehensive understanding of existing membrane preparation technologies. It also provides a guide for developing the new membrane preparation process at the molecular level.

3D-printing↗
DOE OSTI2019

Studies of NMR, molecular docking, and molecular dynamics simulation of new promising inhibitors of cruzaine from the parasite Trypanosoma cruzi

Cruzaine is the major cysteine protease of Trypanosoma cruzi. Cruzaine is involved throughout the parasite’s life cycle in host cells, and is a promising target in the search for new antichagasic agents. Quantum chemical calculations based on density functional theory (DFT B3LYP/cc-pVDZ) were performed to obtain nuclear magnetic resonance data and to optimize the geometry of four dihydrochalcones. The results showed good agreement with the experimental data and were used to suggest the relative stereochemistry of one of the four dihydrochalcones studied. In addition, we evaluated the interaction of cruzaine with these new inhibitors. We used molecular dynamics simulations, free energy calculations, and a per-residue energy decomposition method. It was observed that these molecules are capable of interacting with residues important for enzymatic activity, like Cys25, His161, and Asp160. The ranking of the inhibitors obtained from the binding free energy calculations is in agreement with that experimentally reported. The evaluation of the energy components involved in these calculations demonstrated that the van der Waals term is the major contributor to the drug–receptor stabilizing interactions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Active learning of a crystal plasticity flow rule from discrete dislocation dynamics simulations

Continuum-scale material deformation models, such as crystal plasticity (CP), can significantly enhance their predictive accuracy by incorporating input from lower-scale (i.e. mesoscale) models. The procedure to generate and extract the relevant information is however typically complex and ad hoc, involving decision and intervention by domain experts, leading to long development times. In this study, we develop a principled approach for calibration of continuum-scale models using lower scale information by representing a CP flow rule as a Gaussian process model. This representation allows for efficient parameter space exploration, guided by the uncertainty embedded in the model through a process known as Bayesian optimization (BO). We demonstrate a semi-autonomous BO loop which instantiates discrete dislocation dynamics simulations whose initial conditions are automatically chosen to optimize the uncertainty of a model CP flow rule. Our self-guided computational pipeline efficiently generated a dataset and corresponding model whose error, uncertainty, and physical feature sensitivities were validated with comparison to an independent dataset four times larger, demonstrating a valuable and efficient active learning implementation readily transferable to similar material systems.

36 MATERIALS SCIENCE↗