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 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Pseudoelastic deformation in Mo-based refractory multi-principal element alloys

Phase diagrams supported by density functional theory methods can be crucial for designing high-entropy alloys that are subset of multi-principal-element alloys. We present phase and property analysis of quinary (MoW) x Zr y (TaTi) 1-x-y refractory high-entropy alloys from combined Calculation of Phase Diagram (CALPHAD) and density-functional theory results, supplemented by molecular dynamics simulations. Both CALPHAD and density-functional theory analysis of phase stability indicates a Mo-W-rich region of this quinary has a stable single-phase body-centered-cubic structure. We report first quinary composition from Mo-W-Ta-Ti-Zr family of alloy with pseudo-elastic behavior, i.e., hysteresis in stress-strain. Our analysis shows that only Mo-W-rich compositions of Mo-W-Ta-Ti-Zr, i.e., Mo+W, show reproducible hysteresis in stress-strain responsible for pseudo-elastic behavior. The (MoW) 85 Zr 7.5 (TaTi) 7.5 was down-selected based on temperature-dependent phase diagram analysis and molecular dynamics simulations predicted elastic behavior that reveals twinning-assisted pseudoelastic behavior. While mostly unexplored in body-centered-cubic crystals, twinning is a fundamental deformation mechanism that competes against dislocation slip in crystalline solids. Furthermore, this alloy shows identical cyclic deformation characteristics during uniaxial < 100 > loading, i.e., the pseudoelasticity is isotropic in loading direction. Additionally, a temperature increase from 77 to 1,500 K enhances the elastic strain recovery in load-unload cycles, offering possibly control to tune the pseudoelastic behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2015

A Serially-Connected Compensator for Eliminating the Unbalanced Three-Phase Voltage Impact on Wind Turbine Generators

Untransposed transmission lines, unbalanced tap changer operations, and unbalanced loading in weak distribution lines can cause unbalanced-voltage conditions. The resulting unbalanced voltage at the point of interconnection affects proper gird integration and reduces the lifetime of wind turbines due to power oscillations, torque pulsations, mechanical stresses, energy losses, and uneven and overheating of the generator stator winding. This work investigates the dynamic impact of unbalanced voltage on the mechanical and electrical components of integrated Fatigue, Aerodynamics, Structures, and Turbulence (FAST) wind turbine generation systems (WTGs) of Type 1 (squirrel-cage induction generator) and Type 3 (doubly-fed induction generator). To alleviate this impact, a serially-connected compensator for a three-phase power line is proposed to balance the wind turbine-side voltage. Dynamic simulation studies are conducted in MATLAB/Simulink to compare the responses of these two types of wind turbine models under normal and unbalanced-voltage operation conditions and demonstrate the effectiveness of the proposed compensator.

17 WIND ENERGY↗
DOE OSTI2018

Experimental characterization and modeling of optical tweezer particle handling dynamics [Modeling and Simulation of Optical Tweezer Particle Handling Dynamics]

Here, we report a new framework for a quantitative understanding of optical trapping (OT) particle handling dynamics. We present a novel three-dimensional particle-based model that includes optical, hydrodynamic, and inter-particle forces. This semi-empirical colloid model is based on an open-source simulation code known as LAMMPS (large-scale atomic/molecular massively parallel simulator) and properly recapitulates the full OT force profile beyond the typical linear approximations valid near the trap center. Simulations are carried out with typical system parameters relevant for our experimental holographic optical trapping (HOT) system, including varied particle sizes, trap movement speeds, and beam powers. Furthermore, we present a new experimental method for measuring both the stable and metastable boundaries of the optical force profile to inform or validate the model’s underlying force profile. We show that our framework is a powerful tool for accurately predicting particle behavior in a practical experimental OT setup and can be used to characterize and predict particle handling dynamics within any arbitrary OT force profile.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Simulation of injector dynamics during steady inductive helicity injection current drive in the HIT-SI experiment

We present simulations of inductive helicity injection in the Helicity Injected Torus with Steady Inductive helicity injection (HIT-SI) device that treats the entire plasma volume in a single dynamic MHD model. A new fully 3D numerical tool, the PSI-center TETrahedral mesh code, was developed that provides the geometric flexibility required for this investigation. Implementation of a zero-β Hall MHD model using PSI-TET will be presented including formulation of a new self-consistent magnetic boundary condition for the wall of the HIT-SI device. Results from simulations of HIT-SI are presented focusing on injector dynamics that are investigated numerically for the first time. Asymmetries in the plasma loading between the two helicity injectors and progression of field reversal in each injector are observed. Analysis indicates cross-coupling between injectors through confinement volume structures. Injector impedance is found to scale with toroidal current at fixed density, consistent with experimental observation. Comparison to experimental data with an injector drive frequency of 14.5 kHz shows good agreement with magnetic diagnostics. Global mode structures from Bi-Orthogonal decomposition agree well with experimental data for the first four modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Spatial description of dislocation nucleation in the shock response of single-crystal aluminum

Nonequilibrium molecular dynamics simulations of shock loaded single-crystal Al in the $\langle$100$\rangle$, $\langle$110$\rangle$, $\langle$111$\rangle$, and $\langle$123$\rangle$ orientations are conducted to study elastic and plastic shockwave formation and details associated with dislocation activity. A computer vision-based approach is implemented to capture the presence of dislocations and describe their spatial characteristics in the zone of nucleation behind the propagating shockwave. The methodology developed relies on the sequences of images extracted during shock loading that show dislocation activity within a cross section of the sample. Results reveal that the spacing between activated slip systems is orientation dependent and exhibits a modest reduction for the $\langle$100$\rangle$ and $\langle$111$\rangle$ orientations as shock pressure increases. Comparisons are made to existing theoretical models. Such relationships between shock pressure and dislocation activity, extracted from molecular dynamics simulations, can be used to inform higher length scale simulations or modeling of dislocation-based plasticity during shock.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Low-Resolution Models for the Interaction Dynamics of Coated Gold Nanoparticles with β2-microglobulin

A large number of low-resolution models have been proposed in the last decades to reduce the computational cost of molecular dynamics simulations for bio-nano systems, such as those involving the interactions of proteins with functionalized nanoparticles (NPs). For the proteins, “minimalist” models at the one-bead-per residue (Cα-based) level and with implicit solvent are well established. For the gold NPs, widely explored for biotechnological applications, mesoscale (MS) models treating the NP core with a single spheroidal object are commonly proposed. In this representation, the surface details (coating, roughness, etc.) are lost. These, however, and the specificity of the functionalization, have been shown to have fundamental roles for the interaction with proteins. We presented a mixed-resolution coarse-grained (CG) model for gold NPs in which the surface chemistry is reintroduced as superficial smaller beads. We compared molecular dynamics simulations of the amyloid β2-microglobulin represented at the minimalist level interacting with NPs represented with this model or at the MS level. Our finding highlights the importance of describing the surface of the NP at a finer level as the chemical-physical properties of the surface of the NP are crucial to correctly understand the protein-nanoparticle association.

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

The mobility of small vacancy/helium complexes in tungsten and its impact on retention in fusion-relevant conditions

Tungsten is a promising plasma facing material for fusion reactors. Despite many favorable properties, helium ions incoming from the plasma are known to dramatically affect the microstructure of tungsten, leading to bubble growth, blistering, and/or to the formation of fuzz. In order to develop mitigation strategies, it is essential to understand the atomistic processes that lead to bubble formation and subsequent microstructural changes. In this work, we use large-scale Accelerated Molecular Dynamics simulations to investigate small (N=1,2) V N He M vacancy/helium complexes, which serve as the nuclei for larger helium bubble growth, over timescales reaching into the milliseconds under conditions typical of the operation of fusion reactors. These complexes can interconvert between different I L V N+L He M variants via Frenkel pair nucleation (leading to the creation of a additional vacancy/interstitial pair) and annihilation events; sequences of these events can lead to net migration of these embryonic bubbles. The competition between nucleation and annihilation produces a very complex dependence of the diffusivity on the number of heliums. Finally, through cluster dynamics simulations, we show that diffusion of these complexes provides an efficient pathway for helium release at fluxes expected in fusion reactors, and hence that accounting for the mobility of these complexes is crucial.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Connecting Structural and Transport Properties of Ionic Liquids with Cationic Oligoether Chains

We used X-ray diffraction and molecular dynamics simulations to probe the structures of two families of ionic liquids containing oligoether tails on the cations. Imidazolium and pyrrolidinium bis(trifluoromethylsulfonyl)amide ILs with side chains ranging from 4 to 10 atoms in length, including both linear alkyl and oligo-ethylene oxide tails, were prepared. Furthermore, their physical properties, such as viscosity, conductivity and thermal profile, were measured and compared for systematic trends. Consistent with earlier literature, a single ether substituent substantially decreases the viscosity of pyrrolidinium and imidazolium ILs compared to their alkyl congeners. Remarkably, as the number of ether units in the pyrrolidinium ILs increases there is hardly any increase in the viscosity, in contrast to alkylpyrrolidinium ILs where the viscosity increases steadily with chain length. Viscosities of imidazolium ether ILs increase with chain length but always remain well below their alkyl congeners. To complement the experimentally determined properties, molecular dynamics simulations were run on the two ILs with the longest ether chains. Our results point to specific aspects that could be useful for researchers designing ILs for specific applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Validation of a Mid-Fidelity Approach for Aircraft Stability and Control Characterization

This paper describes an efficient computational approach for aircraft aerodynamic model development intended for use in flight dynamics simulations. A commercial surface-vorticity flow solver called FlightStream is employed to predict performance, stability, and control characteristics for a NASA subscale electric vertical takeoff and landing aircraft in its isolated-airframe configuration. A framework is described to rapidly convert OpenVSP geometry into a form compatible with the FlightStream software. FlightStream predictions for the variation of aerodynamic force and moment coefficients with airflow angles and control surface deflection angles are compared to static wind-tunnel data to illustrate the process and assess the accuracy of FlightStream solutions. The results and approach presented in this paper are encouraging for use in rapid aerodynamic modeling and flight dynamics simulation development early in the aircraft design process.

aerodynamic modeling↗
DOE OSTI2023

Design of high polarization low switching barrier hybrid improper ferroelectric perovskite oxide superlattices

Hybrid improper ferroelectricity is a useful tool to design ABO 3 /A'BO 3 polar superlattices from non polar building blocks. In this study, we have designed high polarization-low switching barrier hybrid improper ferroelectric superlattices with efficient polarization, and polarization-magnetization switching properties above room temperature, using density functional theory and ab initio molecular dynamics simulations. Superlattices with a chemical formula of (AAlO 3 ) m /(A'AlO 3 ) n , where m/n = 1/1, 1/3, 3/1, 1/5 and 5/1, A, A' = Lanthanide and Y cations are considered to outline the design principles behind polarization switching and (LaFeO 3 ) 3 /(CeFeO 3 ) 1 is investigated for polarization-magnetization switching. We find that the unconventional switching paths via out-of-phase rotation Q R- (a 0 a 0 c - ) and tilt precession QTP always yield lower switching barrier compared to those via in-phase rotation Q R+ (a 0 a 0 c + ) and tilt QT (a - a - c 0 ) of BO 6 octahedra. Results from ab initio molecular dynamics simulations estimate the temperature at which the lowest energy barrier can be overcome. Furthermore, it is possible to tune the polarization switching barrier by tuning the tolerance factor, A,A' cation radius mismatch and super lattice periodicity. For switching via Q R- , the switching barrier varies exponentially with rotation angle, indicating how high switching barrier is expected for systems, away from cubic symmetry. We provide a recipe to overcome such a bottleneck by tuning superlattice periodicity. Finally, we have proposed the multiferroic device application concept through a proposed polarization-temperature hysteresis loop and magnetization switching.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Anomalous entropy-driven kinetics of dislocation nucleation

The kinetics of dislocation reactions, such as dislocation multiplication, controls the plastic deformation in crystals beyond their elastic limit, therefore critical mechanisms in a number of applications in materials science. We present a series of large-scale molecular dynamics simulations that shows that one such type of reactions, the nucleation of dislocation at free surfaces, exhibit unconventional kinetics, including unexpectedly large nucleation rates under compression, very strong entropic stabilization under tension, as well as strong non-Arrhenius behavior. These unusual kinetics are quantitatively rationalized using a variational transition state theory approach coupled with an efficient numerical scheme for the estimation of vibrational entropy changes. These results highlight the need for a variational treatment of the kinetics to quantitatively capture dislocation reaction kinetics, especially at low-to-moderate strains where large deformations are required to activate reactions. These observations suggest possible explanations to previously observed unconventional deformation kinetics in both molecular dynamics simulations and experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Tunable quadruple-well ferroelectric van der Waals crystals

The family of layered thio- and seleno-phosphates has gained attention as potential control dielectrics for the rapidly growing family of two-dimensional and quasi-two-dimensional electronic materials. Here we report a combination of density functional theory calculations, quantum molecular dynamics simulations and variable-temperature, -pressure and -bias piezoresponse force microscopy data to predict and verify the existence of an unusual ferroelectric property—a uniaxial quadruple potential well for Cu displacements—enabled by the van der Waals gap in copper indium thiophosphate (CuInP 2 S 6 ). The calculated potential energy landscape for Cu displacements is strongly influenced by strain, accounting for the origin of the negative piezoelectric coefficient and rendering CuInP 2 S 6 a rare example of a uniaxial multi-well ferroelectric. Experimental data verify the coexistence of four polarization states and explore the temperature-, pressure- and bias-dependent piezoelectric and ferroelectric properties, which are supported by bias-dependent molecular dynamics simulations. As a result, these phenomena offer new opportunities for both fundamental studies and applications in data storage and electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Estimated Fluid Force and Damping Characteristics of a Thin Film Damper Comparison Between Closed-Form Solutions and Numerical Analysis

Analytical solutions of thin film dampers are useful for determining critical speeds and stability of rotor systems. Most thin film dampers in use are of short axial length, and closed-form solutions to the Reynolds equations exist for estimating pressure, forces, and damping for these types of dampers. This article compares the fluid film forces and damping estimated by the short film bearing model form of the Reynolds equations to the calculated forces and damping of a transient computational fluid dynamic simulation. For this comparison, the fluid was assumed to be incompressible, laminar, and isoviscous. The fluid film forces and damping are calculated from integrating the pressure distribution over the surface of the damper due to small amplitude motions about a steady state static off-center circular orbit. In this case, no cavitation is assumed, and the journal has no angular velocity, so direct stiffness cannot be calculated from the closed-form solution. Radial clearance, journal length, and journal eccentricity have a significant effect on fluid force and damping within a thin film damper. Fluid density does not affect fluid force or damping substantially, while fluid viscosity does. Both the closed-form solutions and computational fluid dynamics simulation compare well with each other and reflect these trends.

Cook, Jason↗
NASA NTRS2015

Fracture of Carbon Nanotube - Amorphous Carbon Composites: Molecular Modeling

Carbon nanotubes (CNTs) are promising candidates for use as reinforcements in next generation structural composite materials because of their extremely high specific stiffness and strength. They cannot, however, be viewed as simple replacements for carbon fibers because there are key differences between these materials in areas such as handling, processing, and matrix design. It is impossible to know for certain that CNT composites will represent a significant advance over carbon fiber composites before these various factors have been optimized, which is an extremely costly and time intensive process. This work attempts to place an upper bound on CNT composite mechanical properties by performing molecular dynamics simulations on idealized model systems with a reactive forcefield that permits modeling of both elastic deformations and fracture. Amorphous carbon (AC) was chosen for the matrix material in this work because of its structural simplicity and physical compatibility with the CNT fillers. It is also much stiffer and stronger than typical engineering polymer matrices. Three different arrangements of CNTs in the simulation cell have been investigated: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. The SWNT and MWNT array systems are clearly idealizations, but the SWNT bundle system is a step closer to real systems in which individual tubes aggregate into large assemblies. The effect of chemical crosslinking on composite properties is modeled by adding bonds between the CNTs and AC. The balance between weakening the CNTs and improving fiber-matrix load transfer is explored by systematically varying the extent of crosslinking. It is, of course, impossible to capture the full range of deformation and fracture processes that occur in real materials with even the largest atomistic molecular dynamics simulations. With this limitation in mind, the simulation results reported here provide a plausible upper limit on achievable CNT composite properties and yield some insight on the influence of processing conditions on the mechanical properties of CNT composites.

Jensen, Benjamin D.↗
DOE OSTI2026

Dataset_for_Conserved_macromolecular_architecture_of_Poplar_secondary_cell_walls_revealed_by_ssNMR_and_atomistic_modeling

This dataset contains solid-state 13C NMR data and atomistic molecular dynamics simulation files supporting the study of nanoscale secondary cell wall architecture across 13 genetically diverse Populus trichocarpa genotypes grown under uniform greenhouse conditions in 13C-enriched CO2 atmospheres (~89% 13C enrichment).The dataset contains two collections of solid-state 13C NMR data. (1) 200 MHz data (Bruker Avance III HD, 4 mm HX probe, 10 kHz MAS): raw Bruker TopSpin experiment folders and DMFIT-exported ascii spectra for selective and non-selective 1D 13C-13C spin diffusion experiments (3000 ms mixing) used to quantify inter-polymer spatial proximities, and short-mixing (1 ms) reference spectra used for polymeric abundance quantification by spectral deconvolution. (2) 600 MHz data (Bruker Avance III, 1.6 mm PhoenixNMR HXY probe, 30 kHz MAS): raw Bruker TopSpin experiment folders containing 2D CORD, 2D CP-INADEQUATE, and 13C/1H relaxation (T1, T1rho) experiments for all 13 genotypes, with processed Excel workbooks per experiment type. Molecular dynamics simulation code, coordinate files, and analysis scripts (NAMD/CHARMM/Python) for six atomistic cell wall models are included. Summarized ssNMR data are compiled into a single excel file and subjected to statistical analysis. Multivariate analysis code (PCA, Pearson correlation) and summary data are provided as excel worksheets and Jupyter notebooks (Python 3).

09 BIOMASS FUELS↗
DOE OSTI2019

Density functional theory based neural network force fields from energy decompositions

In order to develop force fields (FF) for molecular dynamics simulations that retain the accuracy of ab initio density functional theory (DFT), in this work we developed a machine learning protocol based on an energy decomposition scheme that extracts atomic energies from DFT calculations. Our DFT to FF (DFT2FF) approach provides almost hundreds of times more data for the DFT energies, which dramatically improves accuracy with less DFT calculations. In addition, we use piecewise cosine basis functions to systematically construct symmetry invariant features into the neural network model. We illustrate this DFT2FF approach for amorphous silicon where only 800 DFT configurations are sufficient to achieve an accuracy of 1 meV/atom for energy and 0.1 eV/A for forces. We then use the resulting FF model to calculate the thermal conductivity of amorphous Si based on long molecular dynamics simulations. The dramatic speedup in training in our DFT2FF protocol allows the adoption of a simulation paradigm where an accurate and problem specific FF for a given physics phenomenon is trained on-the-spot through a quick DFT precalculation and FF training.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗