Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamics Simulations”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
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 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 OSTI2022

Ionic transport in electrostatic Janus Membranes. An explicit solvent molecular dynamic simulation.

Janus –or two-sided, charged membranes offer promise as ionic current rectifiers. In such systems, pores consisting of two regions of opposite charge can be used to generate a current from a gradient in salinity. The efficiency of nanoscale Janus pores increases dramatically as their diameter becomes smaller. However, little is known about the underlying transport processes, particularly under experimentally accessible conditions. In this work, we examine the molecular basis for rectification in Janus nanopores using an applied electric field. Molecular simulations with explicit water and ions are used to examine the structure and dynamics of all molecular species in aqueous electrolyte solutions. For several macroscopic observables, the results of such simulations are consistent with experimental observations on asymmetric membranes. Our analysis reveals a number of previously unknown features, including a pronounced local re-orientation of water molecules in the pores, and a segregation of ionic species that has not been anticipated by previously reported continuum analyses of Janus pores. Using these insights, a model is proposed for ionic current rectification in which electric leakage at pore entrance controls net transport.

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

Evaluation of Methods for Molecular Dynamics Simulation of Ionic Liquid Electric Double Layers

We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.

Haskins, Justin B.↗
DOE OSTI2020

Imaging phonon dynamics with ultrafast electron microscopy: Kinematical and dynamical simulations

Ultrafast x-ray and electron scattering techniques have proven to be useful for probing the transient elastic lattice deformations associated with photoexcited coherent acoustic phonons. Indeed, femtosecond electron imaging using an ultrafast electron microscope (UEM) has been used to directly image the influence of nanoscale structural and morphological discontinuities on the emergence, propagation, dispersion, and decay behaviors in a variety of materials. Here, we describe our progress toward the development of methods ultimately aimed at quantifying acoustic-phonon properties from real-space UEM images via conventional image simulation methods extended to the associated strain-wave lattice deformation symmetries and extents. Using a model system consisting of pristine single-crystal Ge and a single, symmetric Lamb-type guided-wave mode, we calculate the transient strain profiles excited in a wedge specimen and then apply both kinematical- and dynamical-scattering methods to simulate the resulting UEM bright-field images. While measurable contrast strengths arising from the phonon wavetrains are found for optimally oriented specimens using both approaches, incorporation of dynamical scattering effects via a multi-slice method returns better qualitative agreement with experimental observations. Contrast strengths arising solely from phonon-induced local lattice deformations are increased by nearly an order of magnitude when incorporating multiple electron scattering effects. We also explicitly demonstrate the effects of changes in global specimen orientation on the observed contrast strength, and we discuss the implications for increasing the sophistication of the model with respect to quantification of phonon properties from UEM images.

36 MATERIALS SCIENCE↗
NASA NTRS2019

Molecular Dynamics Simulations of Silicon Carbide, Boron Nitride and Silicon for Ceramic Matrix Composite Applications

A comprehensive computational molecular dynamics study is presented for crystalline α-SiC (6H, 4H, and 2H SiC), β-SiC (3C SiC), layered boron nitride, amorphous boron nitride and silicon, the constituent materials for high-temperature SiC/SiC compositions. Large-scale Atomic/Molecular Parallel Simulator software package was used. The Tersoff Potential force field was utilized to evaluate their mechanical characteristics of most of the materials, and the Reax force field was used to model silicon when the Tersoff Potential did not provide accurate results. Their mechanical behaviors were evaluated at a strain rate of 10(exp 7)/s and the results agree with the experimental data in the literature. The results are foundational for linking constituent behavior to composite performance, particularly when test data is unavailable or suspect.

Aluko, Olanrewaju↗
NASA NTRS2006

Dynamics Simulation Model for Space Tethers

This document describes the development of an accurate model for the dynamics of the Momentum Exchange Electrodynamic Reboost (MXER) system. The MXER is a rotating tether about 100-km long in elliptical Earth orbit designed to catch payloads in low Earth orbit and throw them to geosynchronous orbit or to Earth escape. To ensure successful rendezvous between the MXER tip catcher and a payload, a high-fidelity model of the system dynamics is required. The model developed here quantifies the major environmental perturbations, and can predict the MXER tip position to within meters over one orbit.

Levin, E. M.↗
NASA NTRS1986

Design and dynamic simulation of a fixed pitch 56 kW wind turbine drive train with a continuously variable transmission

The dynamic analysis of a horizontal axis fixed pitch wind turbine generator (WTG) rated at 56 kW is discussed. A mechanical Continuously Variable Transmission (CVT) was incorporated in the drive train to provide variable speed operation capability. One goal of the dynamic analysis was to determine if variable speed operation, by means of a mechanical CVT, is capable of capturing the transient power in the WTG/wind environment. Another goal was to determine the extent of power regulation possible with CVT operation.

Gallo, C.↗
NASA NTRSDate not supplied

Computational Fluid Dynamics Simulation of A Hazardous Gas Purge System in an EUS Aft Adapter Umbilical Plate Concept

The Fluid Dynamics Branch at the NASA Marshall Space Flight Center actively supports flight programs through a variety of in-line analyses. One area of expertise is the mitigation of hazardous gas mixtures that may form during nominal vehicle fueling operation. After the successful Artemis I launch, design cycles for future iterations of the vehicle are in progress. The subject of this work concerns the analysis of a hazardous gas mitigation system concept for a future upper stage.

Riley D Bridges↗
NASA NTRSDate not supplied

Computational Fluid Dynamics Simulation of A Hazardous Gas Purge System in an EUS Aft Adapter Umbilical Plate Concept

The Fluid Dynamics Branch at the NASA Marshall Space Flight Center actively supports flight programs through a variety of in-line analyses. One area of expertise is the mitigation of hazardous gas mixtures that may form during nominal vehicle fueling operation. After the successful Artemis I launch, design cycles for future iterations of the vehicle are in progress. The subject of this work concerns the analysis of a hazardous gas mitigation system concept for a future upper stage.

Riley Bridges↗
DOE OSTI2023

$\mathrm{DRAGON}$: A multi-GPU orbital-free density functional theory molecular dynamics simulation package for modeling of warm dense matter

As progress in electronic structure theoretical methods is made, ab initio molecular dynamics (MD) based on orbital-free density functional theory (OF-DFT) is becoming increasingly more successful at substituting the traditional, very accurate but computationally costly Kohn–Sham (KS) approach for simulations of matter at the challenging warm dense matter (WDM) regime. However, despite the significant cost alleviation of eliminating the dependence on the KS orbitals, OF-DFT MD runs require ~10 2 to 10 3 CPU cores running for days, or even weeks, for simulations of systems comprised of 10 2 to 10 3 atoms, depending on thermodynamic conditions. Here we present DRAGON, a multi-GPU OF-DFT MD code for fast and efficient simulations of WDM. With a relatively small allocation of resources (4 to 8 GPU devices) it can provide an order of magnitude speedup for simulations containing $\mathscr{O}$(10 4 ) atoms and target systems composed of $\mathscr{O}$(10 5 ) atoms at conditions within the WDM regime, which is currently outside the capabilities of CPU codes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

Computational Fluid Dynamics Simulation of the Swordfish Oxygen/Methane Rotating Detonation Rocket Engine

Rotating Detonation Rocket Engines (RDREs) provide several theoretical performance benefits over traditional deflagration combustion rocket engines. RDREs can potentially generate higher specific impulse (ISP) for a given propellant combination enabling significant consumable mass savings for a spaceflight vehicle. Conversely, RDREs are a relatively new technology with a shallower analytical modeling base compared to traditional combustion engine architectures. In this work, we utilize time-accurate, Computational Fluid Dynamics (CFD) modeling to simulate the complex fluid environment inside an RDRE combustion chamber and propellant manifolds. The selected engine of study is the Swordfish RDRE which was developed and tested and at the National Aeronautics and Space Administration (NASA) Marshall Space Flight Center (MSFC). CFD simulations were conducted using the as-tested Swordfish geometry and conditions to facilitate comparison between the analytical predictions and experimental results.

Brian R. Richardson↗
NASA NTRS2001

Textbook Multigrid Efficiency for Computational Fluid Dynamics Simulations

Considerable progress over the past thirty years has been made in the development of large-scale computational fluid dynamics (CFD) solvers for the Euler and Navier-Stokes equations. Computations are used routinely to design the cruise shapes of transport aircraft through complex-geometry simulations involving the solution of 25-100 million equations; in this arena the number of wind-tunnel tests for a new design has been substantially reduced. However, simulations of the entire flight envelope of the vehicle, including maximum lift, buffet onset, flutter, and control effectiveness have not been as successful in eliminating the reliance on wind-tunnel testing. These simulations involve unsteady flows with more separation and stronger shock waves than at cruise. The main reasons limiting further inroads of CFD into the design process are: (1) the reliability of turbulence models; and (2) the time and expense of the numerical simulation. Because of the prohibitive resolution requirements of direct simulations at high Reynolds numbers, transition and turbulence modeling is expected to remain an issue for the near term. The focus of this paper addresses the latter problem by attempting to attain optimal efficiencies in solving the governing equations. Typically current CFD codes based on the use of multigrid acceleration techniques and multistage Runge-Kutta time-stepping schemes are able to converge lift and drag values for cruise configurations within approximately 1000 residual evaluations. An optimally convergent method is defined as having textbook multigrid efficiency (TME), meaning the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in the discretized system of equations (residual equations). In this paper, a distributed relaxation approach to achieving TME for Reynolds-averaged Navier-Stokes (RNAS) equations are discussed along with the foundations that form the basis of this approach. Because the governing equations are a set of coupled nonlinear conservation equations with discontinuities (shocks, slip lines, etc.) and singularities (flow- or grid-induced), the difficulties are many. This paper summarizes recent progress towards the attainment of TME in basic CFD simulations.

Brandt, Achi↗
DOE OSTI2020

Protein Molecular Dynamics Simulations with Approximate QM: What Can We Learn?

Classical force fields are essential for computer simulations of proteins and are typically parameterized to reproduce secondary and tertiary structure of isolated proteins. However, while protein-protein interactions are ubiquitous in nature, they are not considered in parameterization efforts and are far less understood than isolated proteins. A better characterization of intermolecular interactions is widely recognized as a key to revolutionizing drug and therapeutic developments with high-throughput computational screening. Urgently needed is a critical assessment of the performance of modern protein force fields against first-principles electronic structure methods and experiments. In a daring step toward this goal, we here describe a comparison of peptide folding dynamics as predicted by a molecular mechanics force field on the one hand and by an approximate electronic structure quantum mechanical (QM) method based on density-functional tight-binding (DFTB) on the other. We further compare the dynamics from straightforward DFTB simulations with a near-linear scaling version of DFTB for massively parallel computation based on the fragment molecular orbital (FMO-DFTB) method. We illustrate differences between the phenomenology of the folding dynamics from these three methods for a small model peptide, as well as charge polarization and dynamic fluctuations, point out possible correlations and implications for force field developers, and discuss the lessons learned that might become applicable to future predictive high-throughput computer screening for personalized neoantigen cancer therapy.

Irle, Stephan↗
DOE OSTI2021

Anisotropic fluid dynamical simulations of heavy-ion collisions

We present VAH, a (3+1)–dimensional simulation that evolves the far-from-equilibrium quark-gluon plasma produced in ultrarelativistic heavy-ion collisions with anisotropic fluid dynamics. We solve the hydrodynamic equations on an Eulerian grid using the Kurganov–Tadmor algorithm in combination with a new adaptive Runge–Kutta method. Our numerical scheme allows us to start the simulation soon after the nuclear collision, largely avoiding the need to integrate it with a separate pre-equilibrium dynamics module. We test the code's performance by simulating on the Eulerian grid conformal and non-conformal Bjorken flow as well as conformal Gubser flow, whose (0+1)–dimensional solutions are precisely known. Finally, we compare non-conformal anisotropic hydrodynamics to second-order viscous hydrodynamics in central Pb+Pb collisions and find that the former's longitudinal flow profile responds more consistently to the fluid's gradients along the spacetime rapidity direction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗