Search NASA⌕ Search

SEARCH · Search NASA

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Testing thermal conductivity models with equilibrium molecular dynamics simulations of the one-component plasma

Equilibrium molecular dynamics simulations are used to calculate the thermal conductivity of the onecomponent plasma via the Green-Kubo formalism over a broad range of Coulomb coupling strength, 0.1 ≤ Γ ≤180. These simulations address previous discrepancies between computations using equilibrium versus nonequilibrium methods. Analysis of heat flux autocorrelation functions show that very long (6×10 5 $ω^{-1}_{p}$) time series are needed to reduce the noise level to allow ≲2% accuracy. The simulations provide accurate data for Γ ≲ 1. Furthermore, this enables a test of the traditional Landau-Spitzer theory, which is found to agree with the simulations for Γ ≲ 0.3. It also enables tests of theories to address moderate and strong Coulomb coupling. Two are found to provide accurate extensions to the moderate coupling regime of Γ ≲ 10, but none are accurate in the Γ ≳ 10 regime where potential energy transport and coupling between mass flow and stress dominate thermal conduction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Accelerators for Classical Molecular Dynamics Simulations of Biomolecules

Atomistic Molecular Dynamics (MD) simulations provide researchers the ability to model biomolecular structures such as proteins and their interactions with drug-like small molecules with greater spatiotemporal resolution than is otherwise possible using experimental methods. MD simulations are notoriously expensive computational endeavors that have traditionally required massive investment in specialized hardware to access biologically relevant spatiotemporal scales. Our goal is to summarize the fundamental algorithms that are employed in the literature to then highlight the challenges that have affected accelerator implementations in practice. We consider three broad categories of accelerators: Graphics Processing Units (GPUs), Field-Programmable Gate Arrays (FPGAs), and Application Specific Integrated Circuits (ASICs). These categories are comparatively studied to facilitate discussion of their relative trade-offs and to gain context for the current state of the art. We conclude by providing insights into the potential of emerging hardware platforms and algorithms for MD.

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

The Impact of Protonation on Early Translocation of Anthrax Lethal Factor: Kinetics from Molecular Dynamics Simulations and Milestoning Theory

We report atomically detailed molecular dynamics simulations of the permeation of the lethal factor (LF) N-terminal segment through the anthrax channel. The N-terminal chain is unstructured and leads the permeation process for the LF protein. The simulations were conducted in explicit solvent with milestoning theory, making it possible to extract kinetic information from nanosecond to millisecond time scales. Here, we illustrate that the initial event is strongly influenced by the protonation states of the permeating amino acids. While the N-terminal segment passes easily at high protonation state through the anthrax channel (and the Φ clamp), the initial permeation represents a critical step, which can be irreversible and establishes a hook in the channel mouth.

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

Molecular Dynamics Simulations of a Liquid Metal Wetting Zirconia Suboxides

Molecular dynamics simulations were performed of a liquid metal (Ag) interacting with various zirconia suboxides to calculate the work of adhesion between the liquid metal and the suboxides. These simulations were then used to assess the impact of oxygen content on wetting. The results of this work show a wetting transition as oxygen content is reduced.

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

Integrated Turbine-Based Combined Cycle Dynamic Simulation Model

A Turbine-Based Combined Cycle (TBCC) dynamic simulation model has been developed to demonstrate all modes of operation, including mode transition, for a turbine-based combined cycle propulsion system. The High Mach Transient Engine Cycle Code (HiTECC) is a highly integrated tool comprised of modules for modeling each of the TBCC systems whose interactions and controllability affect the TBCC propulsion system thrust and operability during its modes of operation. By structuring the simulation modeling tools around the major TBCC functional modes of operation (Dry Turbojet, Afterburning Turbojet, Transition, and Dual Mode Scramjet) the TBCC mode transition and all necessary intermediate events over its entire mission may be developed, modeled, and validated. The reported work details the use of the completed model to simulate a TBCC propulsion system as it accelerates from Mach 2.5, through mode transition, to Mach 7. The completion of this model and its subsequent use to simulate TBCC mode transition significantly extends the state-of-the-art for all TBCC modes of operation by providing a numerical simulation of the systems, interactions, and transient responses affecting the ability of the propulsion system to transition from turbine-based to ramjet/scramjet-based propulsion while maintaining constant thrust.

Haid, Daniel A.↗
DOE OSTI2018

Molecular dynamics simulation of the solid-liquid interface migration in terbium

In this work, we developed a Tb embedded atom method potential which properly reproduces the liquid structure obtained from the ab initio molecular dynamics simulation, the hexagonal close packed (hcp)-body-centered cubic (bcc) phase transformation, and melting temperatures. At least three crystal phases [hcp, face-centered cubic (fcc), and bcc] described by this potential can coexist with the liquid phase. Thus, the developed potential provides an excellent test bed for studies of the completive phase nucleation and growth in a single component system. The molecular dynamics simulation showed that all crystal phases can grow from the liquid phase close to their melting temperatures. However, in the cases of the hcp and fcc growth from the liquid phase at very large supercoolings, the bcc phase forms at the solid-liquid interface in the close packed orientations in spite of the fact that both hcp and fcc phases are more stable than the bcc phase at these temperatures. Finally, this bcc phase closes the hcp and fcc phase from the liquid such that the remaining liquid solidifies into the bcc phase. The initial hcp phase then slowly continues growing in expense of the bcc phase.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Impact of Molecular Dynamics Simulations on Research and Development of Semiconductor Materials

Atomic scale defects critically limit performance of semiconductor materials. To improve materials, defect effects and defect formation mechanisms must be understood. In this paper, we demonstrate multiple examples where molecular dynamics simulations have effectively addressed these issues that were not well addressed in prior experiments. In the first case, we report our recent progress on modelling graphene growth, where we found that defects in graphene are created around periphery of islands throughout graphene growth, not just in regions where graphene islands impinge as believed previously. In the second case, we report our recent progress on modelling TlBr, where we discovered that under an electric field, edge dislocations in TlBr migrate in both slip and climb directions. The climb motion ejects extensive vacancies that can cause the rapid aging of the material seen in experiments. In the third case, we discovered that the growth of InGaN films on (0001) surfaces suffers from a serious polymorphism problem that creates enormous amounts of defects. Growth on ( $11\bar{2}0$ ) surfaces, on the other hand, results in single crystalline wurtzite films without any of these defects. In the fourth case, we first used simulations to derive dislocation energies that do not possess any noticeable statistical errors, and then used these error-free methods to discover possible misuse of misfit dislocation theory in past thin film studies. Lastly, we highlight the significance of molecular dynamics simulations in reducing defects in the design space of nanostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fast predictions of liquid-phase acid-catalyzed reaction rates using molecular dynamics simulations and convolutional neural networks

The rates of liquid-phase, acid-catalyzed reactions relevant to the upgrading of biomass into high-value chemicals are highly sensitive to solvent composition and identifying suitable solvent mixtures is theoretically and experimentally challenging. We show that the complex atomistic configurations of reactant–solvent environments generated by classical molecular dynamics simulations can be exploited by 3D convolutional neural networks to enable accurate predictions of Brønsted acid-catalyzed reaction rates for model biomass compounds. We develop a 3D convolutional neural network, which we call SolventNet, and train it to predict acid-catalyzed reaction rates using experimental reaction data and corresponding molecular dynamics simulation data for seven biomass-derived oxygenates in water–cosolvent mixtures. We show that SolventNet can predict reaction rates for additional reactants and solvent systems an order of magnitude faster than prior simulation methods. This combination of machine learning with molecular dynamics enables the rapid, high-throughput screening of solvent systems and identification of improved biomass conversion conditions.

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

Nonlinear dynamic simulation of single- and multi-spool core engines

In this paper a new computational method for accurate simulation of the nonlinear dynamic behavior of single- and multi-spool core engines, turbofan engines, and power generation gas turbine engines is presented. In order to perform the simulation, a modularly structured computer code has been developed which includes individual mathematical modules representing various engine components. The generic structure of the code enables the dynamic simulation of arbitrary engine configurations ranging from single-spool thrust generation to multi-spool thrust/power generation engines under adverse dynamic operating conditions. For precise simulation of turbine and compressor components, row-by-row calculation procedures were implemented that account for the specific turbine and compressor cascade and blade geometry and characteristics. The dynamic behavior of the subject engine is calculated by solving a number of systems of partial differential equations, which describe the unsteady behavior of the individual components. In order to ensure the capability, accuracy, robustness, and reliability of the code, comprehensive critical performance assessment and validation tests were performed. As representatives, three different transient cases with single- and multi-spool thrust and power generation engines were simulated. The transient cases range from operating with a prescribed fuel schedule, to extreme load changes, to generator and turbine shut down.

Schobeiri, T.↗
NASA NTRS1990

Gamma ray observatory dynamics simulator in Ada (GRODY)

This experiment involved the parallel development of dynamics simulators for the Gamma Ray Observatory in both FORTRAN and Ada for the purpose of evaluating the applicability of Ada to the NASA/Goddard Space Flight Center's flight dynamics environment. The experiment successfully demonstrated that Ada is a viable, valuable technology for use in this environment. In addition to building a simulator, the Ada team evaluated training approaches, developed an Ada methodology appropriate to the flight dynamics environment, and established a baseline for evaluating future Ada projects.

Source record↗
DOE OSTI2021

Molecular dynamics simulation of metallic Al–Ce liquids using a neural network machine learning interatomic potential

Al-rich Al-Ce alloys have the possibility of replacing heavier steel and cast-irons for use in high-temperature applications. Knowledge about the structures and properties of Al-Ce alloys at liquid state is vital for optimizing the manufacture process to produce desired allows. However, reliable molecular dynamics simulation of Al-Ce alloy systems remains a great challenge due to the lack of accurate Al-Ce interatomic potential. In this work, an artificial neural network (ANN) deep machine learning (ML) method is used to develop a reliable interatomic potential for Al-Ce alloy. Ab initio molecular dynamics (AIMD) simulation data on Al-Ce liquid with small unit cell (~200 atoms) and on the known Al-Ce crystalline compounds are collected to train the interatomic potential using ANN-ML. The obtained ANN-ML model reproduces well the energies, forces, and atomic structure of Al 90 Ce 10 liquid and crystalline phases of Al-Ce compounds in comparison with ab initio results. The developed ANN-ML potential is applied in molecular dynamics simulations to study the structures and properties of metallic Al 90 Ce 10 liquid, which would provide useful insight for guiding experimental process to produce desired Al-Ce allows.

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

SRG110 Stirling Generator Dynamic Simulator Vibration Test Results and Analysis Correlation

The U.S. Department of Energy (DOE), Lockheed Martin (LM), and NASA Glenn Research Center (GRC) have been developing the Stirling Radioisotope Generator (SRG110) for use as a power system for space science missions. The launch environment enveloping potential missions results in a random input spectrum that is significantly higher than historical RPS launch levels and is a challenge for designers. Analysis presented in prior work predicted that tailoring the compliance at the generator-spacecraft interface reduced the dynamic response of the system thereby allowing higher launch load input levels and expanding the range of potential generator missions. To confirm analytical predictions, a dynamic simulator representing the generator structure, Stirling convertors and heat sources was designed and built for testing with and without a compliant interface. Finite element analysis was performed to guide the generator simulator and compliant interface design so that test modes and frequencies were representative of the SRG110 generator. This paper presents the dynamic simulator design, the test setup and methodology, test article modes and frequencies and dynamic responses, and post-test analysis results. With the compliant interface, component responses to an input environment exceeding the SRG110 qualification level spectrum were all within design allowables. Post-test analysis included finite element model tuning to match test frequencies and random response analysis using the test input spectrum. Analytical results were in good overall agreement with the test results and confirmed previous predictions that the SRG110 power system may be considered for a broad range of potential missions, including those with demanding launch environments.

Suarez, Vicente J.↗
NASA NTRS2007

SRG110 Stirling Generator Dynamic Simulator Vibration Test Results and Analysis Correlation

The U.S. Department of Energy (DOE), Lockheed Martin (LM), and NASA Glenn Research Center (GRC) have been developing the Stirling Radioisotope Generator (SRG110) for use as a power system for space science missions. The launch environment enveloping potential missions results in a random input spectrum that is significantly higher than historical radioisotope power system (RPS) launch levels and is a challenge for designers. Analysis presented in prior work predicted that tailoring the compliance at the generator-spacecraft interface reduced the dynamic response of the system thereby allowing higher launch load input levels and expanding the range of potential generator missions. To confirm analytical predictions, a dynamic simulator representing the generator structure, Stirling convertors and heat sources were designed and built for testing with and without a compliant interface. Finite element analysis was performed to guide the generator simulator and compliant interface design so that test modes and frequencies were representative of the SRG110 generator. This paper presents the dynamic simulator design, the test setup and methodology, test article modes and frequencies and dynamic responses, and post-test analysis results. With the compliant interface, component responses to an input environment exceeding the SRG110 qualification level spectrum were all within design allowables. Post-test analysis included finite element model tuning to match test frequencies and random response analysis using the test input spectrum. Analytical results were in good overall agreement with the test results and confirmed previous predictions that the SRG110 power system may be considered for a broad range of potential missions, including those with demanding launch environments.

Lewandowski, Edward J.↗
DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Crystallographic structure and molecular dynamics simulations of the major endoglucanase from Xanthomonas campestris pv. campestris shed light on its oligosaccharide products release pattern

Cellulases are critical enzymatic components for the transformation of plant biomass into fuels, renewable materials and green chemicals. In this work, we determined the crystal structure, pattern of hydrolysis products release, and conducted molecular dynamics simulations of the major endoglucanase from the Xanthomonas campestris pv. campestris (XccCel5A). XccCel5A has a TIM barrel fold with the catalytic site centrally placed in a binding groove surrounded by aromatic side chains. Molecular dynamics simulations show that productive position of the substrate is secured by a network of hydrogen bonds in the four main subsites, which differ in details from homologous structures. Capillary zone electrophoresis and computational studies reveal XccCel5A can act both as endoglucanase and licheninase, but there are preferable arrangements of substrate regarding β-1,3 and β-1,4 bonds within the binding cleft which are related to the enzymatic efficiency.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Molecular dynamics simulations study of nano bubble attachment at hydrophobic surfaces

Bubble attachment phenomena are examined using Molecular Dynamics Simulations (MDS) for the first time. The simulation involves a nitrogen nano bubble containing 906 nitrogen molecules in a water phase with 74,000 water molecules at molybdenite surfaces. During a simulation period of 1 ns, film rupture and displacement occurs. The attached nanobubble at the hydrophobic molybdenite face surface results in a contact angle of about 90º. This spontaneous attachment is due to a “water exclusion zone” at the molybdenite face surface and can be explained by a van der Waals (vdW) attractive force, as discussed in the literature. In contrast, the film is stable at the hydrophilic quartz (001) surface and the bubble does not attach. Contact angles determined from MD simulations are reported, and these results agree well with experimental and MDS sessile drop results. In this way, film stability and bubble attachment are described with respect to interfacial water structure for surfaces of different polarity. Interfacial water molecules at the hydrophobic molybdenite face surface have relatively weak interactions with the surface when compared to the hydrophilic quartz (001) surface, as revealed by the presence of a 3 Å “water exclusion zone” at the molybdenite/water interface. The molybdenite armchair-edge and zigzag-edge surfaces show a comparably slow process for film rupture and displacement when compared to the molybdenite face surface, which is consistent with their relatively weak hydrophobic character.

bubble attachment↗
DOE OSTI2018

Nonadiabatic dynamics simulations of photoexcited urocanic acid

Urocanic acid (UA) is a UV filter found in human skin, which has been linked to photoimmunosuppression and the formation of skin cancer. Its UV-light-induced photoisomerization and radiationless deactivation mechanisms have been addressed previously by static calculations. In this paper, we present nonadiabatic trajectory-surface-hopping dynamics simulations of photoexcited UA using the semiempirical OM2/MRCI methodology and an adaptive-timestep algorithm. We have simulated almost 6000 trajectories, each for a simulation time of 1.6 ps, covering the entire conformational space of the E and Z isomers of both possible tautomers of the isolated neutral form of UA (overall 32 conformers). Here, the initial conditions for the excited-state dynamics were obtained from 1 ns ground-state dynamics simulations.

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

Temperature dependence of the solid-liquid interface free energy of Ni and Al from molecular dynamics simulation of nucleation

The temperature dependence of the solid-liquid interfacial free energy, γ, is investigated for Al and Ni at the undercooled temperature regime based on a recently developed persistent-embryo method. The atomistic description of the nucleus shape is obtained from molecular dynamics simulations. The computed γ shows a linear dependence on the temperature. The values of γ extrapolated to the melting temperature agree well with previous data obtained by the capillary fluctuation method. Here, using the temperature dependence of γ, we estimate the nucleation free energy barrier in a wide temperature range from the classical nucleation theory. The obtained data agree very well with the results from the brute-force molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗