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 1,423 records · Page 79

Vapor–liquid equilibrium of water with the MB-pol many-body potential

Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. Furthermore, these results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational Fluid Dynamics Demonstration of Rigid Bodies in Motion

The Design Analysis Branch (NE-Ml) at the Kennedy Space Center has not had the ability to accurately couple Rigid Body Dynamics (RBD) and Computational Fluid Dynamics (CFD). OVERFLOW-D is a flow solver that has been developed by NASA to have the capability to analyze and simulate dynamic motions with up to six Degrees of Freedom (6-DOF). Two simulations were prepared over the course of the internship to demonstrate 6DOF motion of rigid bodies under aerodynamic loading. The geometries in the simulations were based on a conceptual Space Launch System (SLS). The first simulation that was prepared and computed was the motion of a Solid Rocket Booster (SRB) as it separates from its core stage. To reduce computational time during the development of the simulation, only half of the physical domain with respect to the symmetry plane was simulated. Then a full solution was prepared and computed. The second simulation was a model of the SLS as it departs from a launch pad under a 20 knot crosswind. This simulation was reduced to Two Dimensions (2D) to reduce both preparation and computation time. By allowing 2-DOF for translations and 1-DOF for rotation, the simulation predicted unrealistic rotation. The simulation was then constrained to only allow translations.

Camarena, Ernesto↗

Self-Healing Ceramics and Mechanical Metamaterials

We performed large scale molecular dynamics (MD) simulations and reactive & quantum molecular dynamics simulations to investigate: (1) faceting, grain growth, and crack healing in ceramics; (2) stress-release mechanism and healing in fibers; and (3) mechanical properties of two-dimensional metamaterials. We also designed highly accurate interatomic potentials and forces using artificial neural networks (ANN) and first-principles molecular dynamics simulation data for use in our MD and reactive MD simulations. We demonstrated the efficacy of modeling quantum materials atomistic simulation data using unsupervised deep learning approaches on a quantum computer. These achievements are elaborated in the following sections.

36 MATERIALS SCIENCE↗

Avionics Tether Operations Control

The activities described in this Final Report were authorized and performed under Purchase Order Number H32835D, issued as part of NASA contract number NAS8-00114. The period of performance of this PO was from March 1 to September 30, 2001. The primary work activity was the continued development and updating of the tether dynamic simulation tools GTOSS (Generalized Tethered Object System Simulation) and TSSIM (Tethered Satellite System) and use of these and other tools in the analysis of various tether dynamics problems. Several updated versions of GTOSS were delivered during the period of performance by the author of the simulation, Lang Associates' David Lang. These updates had mainly to do with updated documentation and an updated coordinate system definition to the J2000 standards. This Final Report is organized by the months in which the activities described were performed. The following sections review the Statement of Work (SOW) and activities performed to satisfy it.

Glaese, John R.↗

Simulation of polymerization induced phase separation in model thermosets

Polymerization induced phase separation (PIPS) in a three component thermoset is studied using molecular dynamics simulations of a new coarse-grained thermoset model. The system includes two crosslinker molecules, which differ in their glass transition temperatures (T g ) and chain length and thus have the potential for phase separation. Here, one crosslinker has a high T g corresponding to a rubbery behavior, and simulations were performed for a short length (4 beads) and a long length (33 beads). The resin and other crosslinker have low T g . A coarse-grained model is developed with these features and with interaction parameters determined so that for either rubbery crosslinker length, the system is in the liquid state at the cure temperature. For sufficiently slow reaction rates, the long rubbery molecule exhibits PIPS into a bicontinuous array of nanoscale domains, but the short one does not, reproducing recent experimental results. The simulations demonstrate that the reaction rates must be slow enough to allow diffusion to yield phase separation. Particularly, the reaction rate corresponding to the secondary amine must be very slow, else the structure of crosslinked clusters and the substantially increased diffusion time will prevent PIPS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Anisotropic Mach Cone Aligned Mesh Adaptation for Low Boom Simulations

An automated off-body Mach cone aligned structured curvilinear grid generation procedure is presented for near-field computational fluid dynamics simulations. This procedure combines output-based indicators and mesh redistribution to perform anisotropic mesh adaptation while maintaining Mach cone alignment. Automation is achieved through a novel direction-based adaptation indicator formulation. The adaptation procedure is demonstrated on the JAXA Wing Body geometry and X-59 C608 demonstrator model from the Second and Third AIAA Sonic Boom Prediction Workshops, respectively. It is demonstrated that anisotropic mesh adaptation may result in a greater than fifty percent reduction in resource usage required to achieve the same level of accuracy as uniform and user constructed Mach cone aligned grids for near-field pressure signatures, ground-level overpressure signatures, and loudness metrics.

CST↗

Anisotropic Mach Cone Aligned Mesh Adaptation for Low Boom Simulations

An automated off-body Mach cone aligned structured curvilinear grid generation procedure is presented for near-field computational fluid dynamics simulations. This procedure combines output-based indicators and mesh redistribution to perform anisotropic mesh adaptation while maintaining Mach cone alignment. Automation is achieved through a novel direction-based adaptation indicator formulation. The adaptation procedure is demonstrated on the JAXA Wing Body geometry and X-59 C608 demonstrator model from the Second and Third AIAA Sonic Boom Prediction Workshops, respectively. It is demonstrated that anisotropic mesh adaptation may result in a greater than fifty percent reduction in resource usage required to achieve the same level of accuracy as uniform and user constructed Mach cone aligned grids for near-field pressure signatures, ground-level overpressure signatures, and loudness metrics.

CST↗

Acclimation of Phytoplankton Fe:C Ratios Dampens the Biogeochemical Response to Varying Atmospheric Deposition of Soluble Iron

Abstract Dissolved iron (dFe) plays an important role in regulating marine productivity. In high nutrient, low chlorophyll regions (>33% of the global ocean), iron is the primary growth limiting nutrient, and elsewhere iron can regulate nitrogen fixation by diazotrophs. The link between iron availability and carbon export is strongly dependent on the phytoplankton iron quotas or cellular Fe:C ratios. This ratio varies by more than an order of magnitude in the open ocean and is positively correlated with ambient dFe concentrations in field observations. Representing Fe:C ratios within models is necessary to investigate how ocean carbon cycling will interact with perturbations to iron cycling in a changing climate. The Community Earth System Model ocean component was modified to simulate dynamic, group‐specific, phytoplankton Fe:C that varies as a function of ambient iron concentration. The simulated Fe:C ratios improve the representation of the spatial trends in the observed Fe:C ratios. The acclimation of phytoplankton Fe:C ratios dampens the biogeochemical response to varying atmospheric deposition of soluble iron, compared to a fixed Fe:C ratio. However, varying atmospheric soluble iron supply has first order impacts on global carbon and nitrogen fluxes and on nutrient limitation spatial patterns. Our results suggest that pyrogenic Fe is a significant dFe source that rivals mineral dust inputs in some regions. Changes in dust flux and iron combustion sources (anthropogenic and wildfires) will modify atmospheric Fe inputs in the future. Accounting for dynamic phytoplankton iron quotas is critical for understanding ocean biogeochemistry and projecting its response to variations in atmospheric deposition.

Environmental Sciences & Ecology↗

X-composer: enabling cross-environments in-situ workflows between HPC and cloud

As large-scale scientific simulations and big data analyses become more popular, it is increasingly more expensive to store huge amounts of raw simulation results to perform post-analysis. To minimize the expensive data I/O, "in-situ" analysis is a promising approach, where data analysis applications analyze the simulation generated data on the fly without storing it first. However, it is challenging to organize, transform, and transport data at scales between two semantically different ecosystems due to the distinct software and hardware difference. To tackle these challenges, we design and implement the X-Composer framework. X-Composer connects cross-ecosystem applications to form an "in-situ" scientific workflow, and provides a unified approach and recipe for supporting such hybrid in-situ workflows on distributed heterogeneous resources. X-Composer reorganizes simulation data as continuous data streams and feeds them seamlessly into the Cloud-based stream processing services to minimize I/O overheads. For evaluation, we use X-Composer to set up and execute a cross-ecosystem workflow, which consists of a parallel Computational Fluid Dynamics simulation running on HPC, and a distributed Dynamic Mode Decomposition analysis application running on Cloud. Our experimental results show that X-Composer can seamlessly couple HPC and Big Data jobs in their own native environments, achieve good scalability, and provide high-fidelity analytics for ongoing simulations in real-time.

Wang, Dali↗

Berni Julian Alder, theoretical physicist and inventor of molecular dynamics, 1925–2020

Berni Julian Alder, one of the leading figures in the invention of molecular dynamics simulations used for a wide array of problems in physics and chemistry, died on September 7th, 2020. His career, spanning more than 65 years, transformed statistical mechanics, many body physics, the study of chemistry and the microscopic dynamics of fluids, by making atomistic computational simulation (in parallel with traditional theory and experiment) a new pathway to unexpected discoveries. Among his many honors, the CECAM prize, recognizing exceptional contributions to the simulation of the microscopic properties of matter is named for him. He was awarded the National Medal of Science by President Obama in 2008.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Amplifying Nanoparticle Reinforcement through Low Volume Topologically Controlled Chemical Coupling

Here, we present a streamlined method to covalently bond hydroxylated carbon nanotubes (CNOH) within a polyphenol matrix, all achieved through a direct, solvent-free process. Employing an extremely small concentration of CNOH (0.01% w/w) along with topologically contrasting linkers led to a maximum of 5-fold increase in modulus and a 25% enhancement in tensile strength compared to the unaltered matrix, an order of magnitude greater reinforcement (w/w) compared to state-of-the-art melt-processed nanocomposites. Through dynamic mechanical analysis, low field solid-state nuclear magnetic resonance spectroscopy, and molecular dynamics simulations, we uncovered the profound influence of linker’s conformational degrees of freedom on the segmental dynamics and therefore the material’s properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ligand modulation of the conformational dynamics of the A 2A adenosine receptor revealed by single-molecule fluorescence

G protein-coupled receptors (GPCRs) are the largest class of transmembrane proteins, making them an important target for therapeutics. Activation of these receptors is modulated by orthosteric ligands, which stabilize one or several states within a complex conformational ensemble. The intra- and inter-state dynamics, however, is not well documented. Here, we used single-molecule fluorescence to measure ligand-modulated conformational dynamics of the adenosine A 2A receptor (A 2A R) on nanosecond to millisecond timescales. Experiments were performed on detergent-purified A 2R in either the ligand-free (apo) state, or when bound to an inverse, partial or full agonist ligand. Single-molecule Förster resonance energy transfer (smFRET) was performed on detergent-solubilized A 2A R to resolve active and inactive states via the separation between transmembrane (TM) helices 4 and 6. The ligand-dependent changes of the smFRET distributions are consistent with conformational selection and with inter-state exchange lifetimes ≥ 3 ms. Local conformational dynamics around residue 229 6.31 on TM6 was measured using fluorescence correlation spectroscopy (FCS), which captures dynamic quenching due to photoinduced electron transfer (PET) between a covalently-attached dye and proximal aromatic residues. Global analysis of PET-FCS data revealed fast (150–350 ns), intermediate (50–60 μs) and slow (200–300 μs) conformational dynamics in A 2A R, with lifetimes and amplitudes modulated by ligands and a G-protein mimetic (mini-G s ). Most notably, the agonist binding and the coupling to mini-G s accelerates and increases the relative contribution of the sub-microsecond phase. Molecular dynamics simulations identified three tyrosine residues (Y112, Y288 7.53 , and Y290 7.55 ) as being responsible for the dynamic quenching observed by PET-FCS and revealed associated helical motions around residue 229 6.31 on TM6. As a result, this study provides a quantitative description of conformational dynamics in A 2A R and supports the idea that ligands bias not only GPCR conformations but also the dynamics within and between distinct conformational states of the receptor.

59 BASIC BIOLOGICAL SCIENCES↗

Fast Na diffusion and anharmonic phonon dynamics in superionic Na 3 PS 4

The design of new solid electrolytes (SEs) hinges on identifying and tuning relevant descriptors. Phonons describe the atomic dynamics in crystalline materials and provide a basis to encode possible minimum energy pathways for ion migration, but anharmonic effects can be large in SEs. Identifying and controlling the pertinent phonon modes coupled most strongly with ionic conductivity, and assessing the role of anharmonicity, could therefore pave the way for discovering and designing new SEs via phonon engineering. In this work, we investigate phonons in Na 3 PS 4 and their coupling to fast Na diffusion, using a combination of neutron scattering, ab initio molecular dynamics (AIMD), and extended molecular dynamics based on machine-learned potentials. We identify that anharmonic soft-modes at the Brillouin zone boundary of the anharmonically stabilized cubic phase constitute key phonon modes that control the Na diffusion process in Na 3 PS 4 . We demonstrate how these strongly anharmonic phonon modes enable Na-ions to hop along the minimum energy pathways. Further, the quasi-elastic neutron scattering (QENS) measurements, supplemented with ab initio and machine-learned molecular dynamics simulations, probe the Na diffusivity and diffusion mechanism. These results offer detailed microscopic insights into the dynamic mechanism of fast Na diffusion and provide an avenue to search for further Na solid electrolytes.

42 ENGINEERING↗

Effects of interaction strength of associating groups on linear and star polymer dynamics

A small number of associating groups incorporated onto a polymer backbone have dramatic effects on the mobility and viscoelastic response of the macromolecules in melts. These associating groups assemble, driving the formation of clusters, whose lifetime affects the properties of the polymers. In this study, we probe the effects of the interaction strength on the structure and dynamics of two topologies, linear and star polymer melts, and further investigate blends of associative and non-associating polymers using molecular dynamics simulations. Polymer chains of approximately one entanglement length are described by a bead–spring model, and the associating groups are incorporated in the form of interacting beads with an interaction strength between them that is varied from 1 to 20 k B T. We find that, for all melts and blends, interaction of a few k B T between the associating groups drives cluster formation, where the size of the clusters increases with increasing interaction strength. These clusters act as physical crosslinkers, which slow the chain mobility. Blends of chains with and without associating groups macroscopically phase separate for interaction strength between the associating groups of a few k B T and above. For weakly interacting associating groups, the static structure function S(q) is well fit by functional form predicted by the random phase approximation where a clear deviation occurs as phase segregation takes place, providing a quantitative assessment of phase segregation.

36 MATERIALS SCIENCE↗

Improving the electrostatic design of the Jefferson Lab 300 kV DC photogun

The 300 kV DC high voltage photogun at Jefferson Lab was redesigned to deliver electron beams with a much higher bunch charge and improved beam properties. The original design provided only a modest longitudinal electric field ( E z ) at the photocathode, which limited the achievable extracted bunch charge. To reach the bunch charge goal of approximately few nC with 75 ps full-width at half-maximum Gaussian laser pulse width, the existing DC high voltage photogun electrodes and anode–cathode gap were modified to increase E z at the photocathode. In addition, the anode aperture was spatially shifted with respect to the beamline longitudinal axis to minimize the beam deflection introduced by the non-symmetric nature of the inverted insulator photogun design. We present the electrostatic design of the original photogun and the modified photogun and beam dynamics simulations that predict vastly improved performance. We also quantify the impact of the photocathode recess on beam quality, where recess describes the actual location of the photocathode inside the photogun cathode electrode relative to the intended location. A photocathode unintentionally recessed/misplaced by sub-millimeter distance can significantly impact the downstream beam size.

47 OTHER INSTRUMENTATION↗