Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic 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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Isothermal Compressor Computational Fluid Dynamics Simulations (Final Report)

Carnot Compression is a startup company developing an innovative technology for air and gas compression. This technology is inherently oil-free and isothermal due to the use of water to simultaneously compress and cool the gas throughout the process. Isothermal gas compression eliminates the need to cool the compressed gas so it may lead to significant energy savings. The development of Carnot’s isothermal compressor is limited by the lack of insight into the flow and detailed behavior of the fluids (water and air/gas) inside the air end. Previous and current attempts at Computational Fluid Dynamics (CFD) simulations by Carnot have been unable to provide the level of accuracy required to use simulations for technology development. In this Cooperative Research and Development Agreement (CRADA) project, Oak Ridge National Laboratory (ORNL) used the CFD package Star-CCM+ to successfully develop a CFD simulation, providing the much-needed insight required to speed up the development. The work is intended to enable Carnot Compression to unlock its technology potential to a level sufficient for commercialization of the intended first product. The work will also prepare Carnot to scale the technology to much larger and more energy intensive applications, positioning it to broaden the product applications. Successful development of the technology has the potential to result in 20% or more efficiency gains for compression processes, or about $3 billion in annual energy savings potential for the U.S. alone.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Simulating Dynamical Phases of Chiral p + i p Superconductors with a Trapped ion Magnet

Two-dimensional p + ip superconductors and superfluids are systems that feature chiral behavior emerging from the Cooper pairing of electrons or neutral fermionic atoms with nonzero angular momentum. Their realization has been a longstanding goal because they offer great potential utility for quantum computation and memory. However, they have so far eluded experimental observation both in solid-state systems as well as in ultracold quantum gases. Here, we propose to leverage the tremendous control offered by rotating two-dimensional trapped-ion crystals in a Penning trap to simulate the dynamical phases of two-dimensional p + ip superfluids. This is accomplished by mapping the presence or absence of a Cooper pair into an effective spin-1/2 system encoded in the ions’ electronic levels. We show how to infer the topological properties of the dynamical phases, and discuss the role of beyond mean-field corrections. More broadly, our work opens the door to use trapped-ion systems to explore exotic models of topological superconductivity and also paves the way to generate and manipulate skyrmionic spin textures in these platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A Method for Obtaining Liquid–Solid Adsorption Rates from Molecular Dynamics Simulations: Applied to Methanol on Pt(111) in H 2 O

Adsorption is an important step in heterogeneous catalysis as it predetermines how many reactant molecules can participate in a surface reaction per unit time. While the rate of adsorption processes is well studied in gas–solid adsorption in both theory and experiment, such rates are still not well studied for liquid–solid adsorption. This is partly because the ever-changing configurations of liquid-phase solvent molecules impede the ability to study a molecule approaching a surface from a liquid phase by either experiment or theory. In this work, we develop a method using molecular dynamics (MD) simulations to study the rate of adsorption in liquid–solid adsorption processes. Specifically, we use MD to model the diffusion of a methanol molecule in aqueous solvent and its adsorption to a Pt(111) surface. We find that by approximating the solute motion as following the same displacement rates as a random walk model, the adsorbed and non-adsorbed states of the methanol molecule near the Pt(111) surface can be discerned and quantified. In particular, this methodology enables extracting a sticking coefficient and a macroscopically relatable adsorption rate. This method can be applied to arbitrary types of reactants and surfaces, as well as different liquid environments, thus providing a general tool for predicting quantitative adsorption rates of liquid–solid adsorption systems.

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

Molecular dynamics simulations of the secondary-binding site in disaccharide-modified glycopeptide antibiotics

Oritavancin is a semisynthetic glycopeptide antibiotic used to treat severe infections by multidrug-resistant Gram-positive pathogens. Oritavancin is known to be a thousand times more potent than vancomycin against Gram-positive bacteria due to the additional interactions with bacterial peptidoglycan (PG) facilitated by a secondary-binding site. The presence of this secondary-binding site is evident in desleucyl-oritavancin, an Edman degradation product of oritavancin, still retaining its potency against Gram-positive bacteria, whereas desleucyl-vancomycin is devoid of any antimicrobial activities. Herein, using explicit solvent molecular dynamics (MD) simulations, steered MD simulations, and umbrella sampling, we show evidence of a secondary-binding site mediated by the disaccharide-modified hydrophobic sidechain of oritavancin interactions with the pentaglycyl-bridge segment of the PG. The interactions were characterized through comparison to the interaction of PG with chloroeremomycin, vancomycin, and the desleucyl analogs of the glycopeptides. Our results show that the enhanced binding of oritavancin to PG over the binding of the other complexes studied is due to an increase in the hydrophobic effect, electrostatic and van der Waals interactions, and not the average number of hydrogen bonds. Our ranking of the binding interactions of the biomolecular complexes directly correlates with the order based on their experimental minimum inhibitory concentrations. The results of our simulations provide insight into the modification of glycopeptides to increase their antimicrobial activities or the design of novel antibiotics against pathogenic Gram-positive bacteria.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Molecular dynamics simulation of the effect of cooling rate on the structure and properties of lithium disilicate glass

The effects of cooling rate on the structure and properties of lithium disilicate (LS2) glass are investigated using molecular dynamics (MD) computer simulations. The evolution of structural features such as pair distribution function, bond angle distribution, and Li coordination distribution are determined, and correlated with dynamic and static properties to elucidate the effects of cooling rate. The density, elastic moduli, and diffusion coefficient are found to be highly sensitive to cooling rate, whereas the Si pair distribution function and bond angle distribution are weakly affected by the cooling rate. Additionally, the changes of Si-O-Si bond angle and Li coordination number suggest the formation of Li cluster at lower cooling rates. Furthermore, by comparing results from other simulations and reported experiments, we confirm that the increase of cooling rate leads to an increase of conductivity and a decrease of density. Finally, at very high cooling rates, we find that all the atoms in LS2-glass do not relax simultaneously, but do so in two distinct configurations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular dynamics simulations of displacement cascades in LiAlO2 and LiAl5O8 ceramics

Abstract Molecular dynamics was employed to investigate the radiation damage due to collision cascades in LiAlO 2 and LiAl 5 O 8 , the latter being a secondary phase formed in the former during irradiation. Atomic displacement cascades were simulated by initiating primary knock-on atoms (PKA) with energy values = 5, 10 and 15 keV and the damage was quantified by the number of Frenkel pairs formed for each species: Li, Al and O. The primary challenges of modeling an ionic system with and without a core–shell model for oxygen atoms were addressed and new findings on the radiation resistance of these ceramics are presented. The working of a variable timestep function and the kinetics in the background of the simulations have been elaborated to highlight the novelty of the simulation approach. More importantly, the key results indicated that LiAlO 2 experiences much more radiation damage than LiAl 5 O 8 , where the number of Li Frenkel pairs in LiAlO 2 was 3–5 times higher than in LiAl 5 O 8 while the number of Frenkel pairs for Al and O in LiAlO 2 are ~ 2 times higher than in LiAl 5 O 8 . The primary reason is high displacement threshold energies (E d ) in LiAl 5 O 8 for Li cations. The greater E d for Li imparts higher resistance to damage during the collision cascade and thus inhibits amorphization in LiAl 5 O 8 . The presented results suggest that LiAl 5 O 8 is likely to maintain structural integrity better than LiAlO 2 in the irradiation conditions studied in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Characterizing Porous and Nonporous Phenolic Resins from Molecular Dynamics Simulations

Phenolic resins are an important component of many ablative heat shield materials, which protect spacecrafts from the extreme temperatures reached during atmospheric entry. Examples include the high-density Heritage Carbon Phenolic (HCP) used in the Pioneer-Venus and Galileo missions, as well as the low-density Phenolic Impregnated Carbon Ablator (PICA) used in the Mars Science Laboratory and Mars 2020 missions. Additionally, recent developments within NASA have produced the mid-density Heatshield for Extreme Entry Environment Technology (HEEET) and its derivative 3D Woven Mid-Density Carbon Phenolic (3MDCP). Unlike the nonporous phenolic in HCP, PICA and HEEET/3MDCP are fabricated by infusing preforms with diluted phenolic formulations to obtain a lower density porous matrix. Despite the importance of the phenolic phase to the material response during entry, the variation in properties of porous and nonporous phenolic is not well understood. Here, we present an investigation of porous and nonporous phenolic resins using molecular dynamics (MD) simulations. Resin cure is mimicked in the simulations through the inclusion of representative reaction templates to generate accurate models of the complex crosslinked structures. To create porous models, explicit solvent molecules are included during the cure simulations. We observe nanoscale separation of the phenolic and solvent phases, which results in significant differences in the final structures of porous and nonporous models. In addition to a quantitative assessment of the network structure and porosity, we elucidate the effects of the phenolic formulation on the final material properties. These results are compared with experimental data as appropriate.

phenolic↗
DOE OSTI2023

Applying molecular dynamics simulations to the study of micro cold spray of yttria

Micro cold spray (MCS) is a novel manufacturing process with proven success depositing dense, thick films of technical ceramics without the need for elevated processing temperatures or chemical precursors. The process transforms a loose sub-micrometer powder into an adherent film with mechanical properties approaching that of the bulk material through a phenomenon known as room temperature impact consolidation. While this phenomenon has been demonstrated for a wide variety of ceramic materials and film substrates, yttria stands out as a material that is both well-characterized in experiments and of significant interest for industrial applications due to its excellent refractory properties. This dissertation seeks to investigate the fundamental mechanisms that enable room temperature impact consolidation by using molecular dynamics (MD) simulations of yttria particle impacts and characterizing the deformation behavior of the system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational Fluid Dynamics Simulations to Predict Oxidation in Heat Recovery Steam Generator Tubes

Heat Recovery Steam Generators (HRSGs) are widely used across United States in combined cycle power plants to recover waste heat from the gas turbine (GT). HRSGs are used either to generate electricity or to produce process steam for industrial applications. The primary components of HRSG consists of a duct and a heat exchanger (HX). It is known across the power industry that the high temperature oxidation and ensuing exfoliation problem is a major cause for the damage of HX materials of HRSG. Alloys and/or coatings that can prevent or mitigate oxidation are very expensive, therefore they must be used or applied on the select regions of the HX tubes where the tendency of oxide formation is the highest. The main goal of this project is to identify such regions through Computational Fluid Dynamics (CFD) simulations. Therefore, in this work, we developed a CFD framework using commercial code StarCCM+ for the prediction of the fluid flow and heat transfer in a HRSG and associated oxidation inside the tubes of the HX. The developed CFD framework was verified and validated with experimental data before deployment. We also developed an innovative method to model the effect of the fins on the heat transfer and the pressure drop using a porous media model (PMM) approach to keep the mesh size within reasonable limits. After validation, we performed high-fidelity CFD simulations of a real-scale HRSG using the PMM, with High Performance Computing (HPC) resources of ORNL. From the simulation results, we acquired oxide thickness maps for all the tubes of the select HX sections of HRSG prone to oxidation. These oxide maps can inform regions of the HX tubes that requires oxide-resistant coatings, thereby guiding engineers for cost-efficient manufacturing of the HX that can combat oxidation in HRSGs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

CHARMM-GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field

Explicit treatment of electronic polarizability in empirical force fields (FFs) represents an extension over a traditional additive or pairwise FF and provides a more realistic model of the variations in electronic structure in condensed phase, macromolecular simulations. To facilitate utilization of the polarizable FF based on the classical Drude oscillator model, Drude Prepper has been developed in CHARMM-GUI. Drude Prepper ingests additive CHARMM protein structures file (PSF) and pre-equilibrated coordinates in CHARMM, PDB, or NAMD format, from which the molecular components of the system are identified. These include all residues and patches connecting those residues along with water, ions, and other solute molecules. This information is then used to construct the Drude FF-based PSF using molecular generation capabilities in CHARMM, followed by minimization and equilibration. In addition, inputs are generated for molecular dynamics (MD) simulations using CHARMM, GROMACS, NAMD, and OpenMM. Validation of the Drude Prepper protocol and inputs is performed through conversion and MD simulations of various heterogeneous systems that include proteins, nucleic acids, lipids, polysaccharides, and atomic ions using the aforementioned simulation packages. Stable simulations are obtained in all studied systems, including 5 μs simulation of ubiquitin, verifying the integrity of the generated Drude PSFs. Additionally, the ability of the Drude FF to model variations in electronic structure is shown through dipole moment analysis in selected systems. Finally, the capabilities and availability of Drude Prepper in CHARMM-GUI is anticipated to greatly facilitate the application of the Drude FF to a range of condensed phase, macromolecular systems.

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

Evaluation of 6-OxP-CD, an Oxime-based cyclodextrin as a viable medical countermeasure against nerve agent poisoning: Experimental and molecular dynamic simulation studies on its inclusion complexes with cyclosarin, soman and VX

The ability of the cyclodextrin-oxime construct 6-OxP-CD to bind and degrade the nerve agents Cyclosarin (GF), Soman (GD) and S -[2-[Di(propan-2-yl)amino]ethyl] O -ethyl methylphosphonothioate (VX) has been studied using 31 P-nuclear magnetic resonance (NMR) under physiological conditions. While 6-OxP-CD was found to degrade GF instantaneously under these conditions, it was found to form an inclusion complex with GD and significantly improve its degradation (t 1/2 ~ 2 hrs) relative over background (t 1/2 ~ 22 hrs). Consequently, effective formation of the 6-OxP-CD:GD inclusion complex results in the immediate neutralization of GD and thus preventing it from inhibiting its biological target. In contrast, NMR experiments did not find evidence for an inclusion complex between 6-OxP-CD and VX, and the agent’s degradation profile was identical to that of background degradation (t 1/2 ~ 24 hrs). As a complement to this experimental work, molecular dynamics (MD) simulations coupled with Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) calculations have been applied to the study of inclusion complexes between 6-OxP-CD and the three nerve agents. These studies provide data that informs the understanding of the different degradative interactions exhibited by 6-OxP-CD with each nerve agent as it is introduced in the CD cavity in two different orientations (up and down). For its complex with GF, it was found that the oxime in 6-OxP-CD lies in very close proximity (P GF …O Oxime ~ 4–5 Å) to the phosphorus center of GF in the ‘down GF ’ orientation for most of the simulation accurately describing the ability of 6-OxP-CD to degrade this nerve agent rapidly and efficiently. Further computational studies involving the center of masses (COMs) for both components (GF and 6-OxP-CD) also provided some insight on the nature of this inclusion complex. Distances between the COMs (ΔCOM) lie closer in space in the ‘down GF ’ orientation than in the ‘up GF ’ orientation; a correlation that seems to hold true not only for GF but also for its congener, GD. In the case of GD, calculations for the ‘down GD ’ orientation showed that the oxime functional group in 6-OxP-CD although lying in close proximity (P GD …O Oxime ~ 4–5 Å) to the phosphorus center of the nerve agent for most of the simulation, adopts another stable conformation that increase this distance to ~ 12–14 Å, thus explaining the ability of 6-OxP-CD to bind and degrade GD but with less efficiency as observed experimentally (t 1/2 ~ 4 hr. vs. immediate). Lastly, studies on the VX:6-OxP-CD system demonstrated that VX does not form a stable inclusion complex with the oxime-bearing cyclodextrin and as such does not interact in a way that is conducive to an accelerated degradation scenario. Collectively, these studies serve as a basic platform from which the development of new cyclodextrin scaffolds based on 6-OxP-CD can be designed in the development of medical countermeasures against these highly toxic chemical warfare agents.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Machine learning enabled quantification of the hydrogen bonds inside the polyelectrolyte brush layer probed using all-atom molecular dynamics simulations

The configuration of densely grafted charged polyelectrolyte (PE) brushes is strongly dictated by the properties and behavior of the counterions that screen the PE brush charges and the solvent molecules (typically water) that solvate the brush molecules and these screening counterions. Only recently, efforts have been made to study the PE brushes atomistically, thereby shedding light on the properties of brush-supported ions and water molecules. However, even for such efforts, there are limitations associated with using a generic definition to estimate certain properties of water and ions inside the brush layer. For example, water–water hydrogen bonds (HBs) will behave differently for locations outside and inside the brush layer, given the fact that the densely closely grafted PE brush molecules create a soft nanoconfinement where the water connectivity becomes highly disrupted: therefore, using the same definition to quantify the HBs inside and outside the brush layer will be unwise. In this paper, we address this limitation by employing an unsupervised machine learning (ML) approach to predict the water–water hydrogen bonding inside a cationic PE brush layer modeled using all-atom molecular dynamics (MD) simulations. Here, the ML method, which relies on a clustering approach and uses the equilibrium coordinates of the water molecules (obtained from the all-atom MD simulations) as the input, is capable of identifying the structural modification of water–water HBs (revealed through appropriate clustering of the data) inside the PE brush layer induced soft nanoconfinement. Such capabilities would not have been possible by using a generic definition of the HBs. Our calculations lead to four key findings: (1) the clusters formed inside and outside the brush layer are structurally similar; (2) the margin of the cluster is shorter inside the PE brush layer confirming the possible disruption of the HBs inside the PE brush layer; (3) the average “hydrogen–acceptor-oxygen–donor-oxygen” angle that defines the HB is reduced for the HBs formed inside the brush layer; (4) the use of the generic definition (definition usable for characterizing the HBs in brush-free bulk) leads to an overprediction of the number of HBs formed inside the PE brush layer.

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

Computational strategies in the dynamic simulation of constrained flexible MBS

This research focuses on the computational dynamics of flexible constrained multibody systems. At first a recursive mapping formulation of the kinematical expressions in a minimum dimension as well as the matrix representation of the equations of motion are presented. The method employs Kane's equation, FEM, and concepts of continuum mechanics. The generalized active forces are extended to include the effects of high temperature conditions, such as creep, thermal stress, and elastic-plastic deformation. The time variant constraint relations for rolling/contact conditions between two flexible bodies are also studied. The constraints for validation of MBS simulation of gear meshing contact using a modified Timoshenko beam theory are also presented. The last part deals with minimization of vibration/deformation of the elastic beam in multibody systems making use of time variant boundary conditions. The above methodologies and computational procedures developed are being implemented in a program called DYAMUS.

Amirouche, F. M. L.↗
NASA NTRS1998

Test of a flexible spacecraft dynamics simulator

There are a number of approaches one can take to modeling the dynamics of a flexible body. While one can attempt to capture the full dynamical behavior subject to disturbances from actuators and environmental torques, such a detailed description often is unnecessary. Simplification is possible either by limiting the amplitude of motion to permit linearization of the dynamics equations or by restricting the types of allowed motion. In this work, we study the nonlinear dynamics of bending deformations of wire booms on spinning spacecraft. The theory allows for large amplitude excursions from equilibrium while enforcing constraints on the dynamics to prohibit those modes that are physically less relevant or are expected to damp out fast. These constraints explicitly remove the acoustic modes (i.e., longitudinal sound waves and shear waves) while allowing for arbitrary bending and twisting, motions which typically are of lower frequency. As a test case, a spin axis reorientation maneuver by the Polar Plasma Laboratory (POLAR) spacecraft has been simulated. POLAR was chosen as a representative spacecraft because it has flexible wire antennas that extend to a length of 65 meters. Bending deformations in these antennas could be quite large and have a significant effect on the attitude dynamics of the spacecraft body. Summary results from the simulation are presented along, with a comparison with POLAR flight data.

Dichmann, Donald↗
DOE OSTI2021

Local structure analysis of low-temperature neutron pair distribution function coupled with molecular dynamics simulations of CH 4 and CO 2 hydrates from 2 to 210 K

CH 4 hydrates occur naturally and are an abundant potential fuel source with a corresponding risk of potent greenhouse gas release, due to their stability conditions at low temperature and high pressure. Byproduct CO 2 can be exchanged with CH 4 in these natural deposits and can potentially stabilize the hydrates at higher temperatures. CH 4 , CO 2 , and mixed CH 4 -CO 2 hydrates are studied with in situ neutron pair distribution function experiments from 2 to 210 K to investigate the impact of varying the CH 4 -CO 2 guest composition in the gas hydrate structure. These experiments combined with Reverse Monte Carlo analysis allow for the characterization of intermolecular CH 4 and CO 2 interactions with the water molecules which form the hydrate lattice and how they impact the local structure of the lattice itself. Results indicate that when CH 4 and CO 2 co-occupy the hydrate, the host is more strongly distorted than in the pure CH 4 and CO 2 hydrates, but this becomes less defined with increasing temperature. The presence of CO 2 in mixed hydrate increases the stability range and creates a barrier for CH 4 to completely leave the structure.

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

Towards computational polar-topotronics: Multiscale neural-network quantum molecular dynamics simulations of polar vortex states in SrTiO3/PbTiO3 nanowires

Recent discoveries of polar topological structures ( e.g ., skyrmions and merons) in ferroelectric/paraelectric heterostructures have opened a new field of polar topotronics. However, how complex interplay of photoexcitation, electric field and mechanical strain controls these topological structures remains elusive. To address this challenge, we have developed a computational approach at the nexus of machine learning and first-principles simulations. Our multiscale neural-network quantum molecular dynamics molecular mechanics approach achieves orders-of-magnitude faster computation, while maintaining quantum-mechanical accuracy for atoms within the region of interest. This approach has enabled us to investigate the dynamics of vortex states formed in PbTiO 3 nanowires embedded in SrTiO 3 . We find topological switching of these vortex states to topologically trivial, uniformly polarized states using electric field and trivial domain-wall states using shear strain. These results, along with our earlier results on optical control of polar topology, suggest an exciting new avenue toward opto-electro-mechanical control of ultrafast, ultralow-power polar topotronic devices.

Linker, Thomas↗
DOE OSTI2021

Accelerated molecular dynamics simulations of dislocation climb in nickel

We report the mechanical behavior of materials operating under high temperatures is strongly influenced by creep mechanisms such as dislocation climb, which is controlled by the diffusion of vacancies. However, atomistic simulations of these mechanisms have traditionally been impractical due to the long time scales required. To overcome these time scale challenges, we use Parallel Trajectory Splicing (ParSplice), an accelerated molecular dynamics method, to simulate dislocation climb in nickel. We focus on modeling the activity of a vacancy near a jog on an edge dislocation in order to observe vacancy pipe diffusion and vacancy absorption at the jog. From rigorously constructed trajectories encompassing more than 2000 vacancy absorption events over a simulation time of more than 4μs at 900 K, a comprehensive sampling of available atomistic mechanisms is collated and analyzed further with molecular statics calculations. We estimate average rates for pipe diffusion and vacancy absorption into the jog using data from the dynamic and static calculations, finding very good agreement. Our results strongly suggest that the dominant mechanism for vacancy absorption by jogs is via biased diffusion to the dislocation core followed by fast pipe diffusion to the jog.

36 MATERIALS SCIENCE↗