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 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simultaneous optimal system and controller design for multibody systems with joint friction using direct sensitivities

Abstract Real-world multibody systems are often subject to phenomena like friction, joint clearances, and external events. These phenomena can significantly impact the optimal design of the system and its controller. This work addresses the gradient-based optimization methodology for multibody dynamic systems with joint friction using a direct sensitivity approach. The Brown–McPhee model has been used to characterize the joint friction in the system. This model is suitable for the study due to its accuracy for dynamic simulation and its compatibility with sensitivity analysis. This novel methodology supports codesign of the multibody system and its controller, which is especially relevant for applications like robotics and servo-mechanical systems, where the actuation and design are highly dependent on each other. Numerical results are obtained using a software package written in Julia with state-of-the-art libraries for automatic differentiation and differential equations. Three case studies are provided to demonstrate the attractive properties of simultaneous optimal design and control approach for certain applications.

Verulkar, Adwait↗
DOE OSTI2022

Computational Investigations of the Reactivity of Metalloporphyrins for Ammonia Oxidation

Density functional theory and molecular dynamics simulations were used to assess the ability of tetraphenylporphyrin (TPP)M complexes (where M = Cr, Mn, Fe, Co, Ni, Mo, Ru, W, and Os) to coordinate and weaken the N-H bonds of ammonia, as well as their reactivity towards N-N bond formation for N 2 generation. Compared to other metalloporphyrins, bis-ammonia complexes (TPP)Mo(NH 3 ) 2 and (TPP)W(NH 3 ) 2 exhibit low and level N-H BDFEs due to a stabilized (TPP)M(NH 3 )(NH) intermediate by multiple metal-ligand bonding. These results resemble those previously obtained for polypyridyl metal complexes, suggesting that broad trends in reactivity towards N-H bond cleavage are more metal-dependent rather than ligand-dependent for a metal in a nitrogen pseudo-octahedral environment. We investigated N-N bond formation via NH 3 nucleophilic attack on M-NH and M-N intermediates, compared to bimolecular coupling of M-NH x intermediates. We evaluated the reactivity of (TPP)Fe(NH 3 ) 2 towards N-N bond formation via a hydrazine pathway, and found amide-amide coupling to form a bridged hydrazido complex to be the most favorable pathway for N-N bond formation. Further investigation of possible N-N bond formation pathways by reaction with NH 3 led us to identify a possible Fe III -•NH species with significant aminyl character that bypasses the nucleophilic attack of NH 3 and that promotes homolytic N-H bond cleavage of ammonia. This reaction forms a Fe-NH 2 moiety and a transient •NH 2 radical that subsequently forms an N-N bond with the Fe-NH 2 moiety to form a (TPP)Fe(NH 3 )(N 2 H 4 ) species. Furthermore, these results indicate the need to evaluate the radical character of imido species and their reactivity towards N-H bond cleavage of ammonia.

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

Structurally Driven Environmental Degradation of Friction in MoS2 Films

Abstract We report an investigation of the friction mechanisms of MoS 2 thin films under changing environments and contact conditions using a variety of computational and experimental techniques. Molecular dynamics simulations were used to study the effects of water and molecular oxygen on friction and bonding of MoS 2 lamellae during initial sliding. Characterization via photoelectron emission microscopy (PEEM) and Kelvin probe force microscopy (KPFM) were used to determine work function changes in shear modified material within the top few nanometers of MoS 2 wear scars. The work function was shown to change with contact conditions and environment, and shown by density functional theory (DFT) calculations and literature reports to be correlated with lamellae size and thickness of the basally oriented surface layer. Results from nanoscale simulations and macroscale experiments suggest that the evolution of the friction behavior of MoS 2 is linked primarily to the formation or inhibition of a basally oriented, molecularly thin surface film with long-range order.

Curry, John F.↗
DOE OSTI2022

Quantifying Equiaxed vs Epitaxial Solidification in Laser Melting of CMSX-4 Single Crystal Superalloy

The competition between epitaxial vs. equiaxed solidification has been investigated in CMSX-4 single crystal superalloy during laser melting as practiced in additive manufacturing. Single-track laser scans were performed on a powder-free surface of directionally solidified CMSX-4 alloy with several combinations of laser power and scanning velocity. Electron backscattered diffraction (EBSD) mapping facilitated identification of new orientations, i.e., “stray grains” that nucleated within the fusion zone along with their area fraction and spatial distribution. Using high-fidelity computational fluid dynamics simulations, both the temperature and fluid velocity fields within the melt pool were estimated. This information was combined with a nucleation model to determine locations where nucleation has the highest probability to occur in melt pools. In conformance with general experience in metals additive manufacturing, the as-solidified microstructure of the laser-melted tracks is dominated by epitaxial grain growth; nevertheless, stray grains were evident in elongated melt pools. It was found that, though a higher laser scanning velocity and lower power are generally helpful in the reduction of stray grains, the combination of a stable keyhole and minimal fluid velocity further mitigates stray grains in laser single tracks.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Support of Adhesion Mechanisms in Al 2 O 3 Aerosol Deposition Through Laser-Induced Particle Impact Testing

Aerosol deposition (AD) is a kinetic spray process capable of depositing ceramic coatings at room temperature, but AD process development is generally a laborious exploration of a large process parameter space. Here, this paper presents a case study investigating whether laser-induced particle impact testing (LIPIT) could be applied to expedite development of an alumina (Al 2 O 3 ) coating on nickel (Ni): Specifically, whether LIPIT measurements could predict critical velocities of adhesion on Ni and Al 2 O 3 , and the effect of ball milling the Al 2 O 3 powder. Because LIPIT has a diffraction-limited lower bound on imageable particle size, the usefulness of Al 2 O 3 powder agglomerates as a proxy for single particles was additionally studied. Overall, LIPIT measurements and AD sprays agreed that ball milling dramatically improves adhesion. Additionally, LIPIT measurements of critical velocity of adhesion of Al 2 O 3 powder agglomerates on Ni and Al 2 O 3 substrates (150 meters per second [m/s] and 250 m/s, respectively) quantitatively agreed with predictions from a previously published model based on picoindentation and molecular dynamics simulations. Together, these findings support the established hypothesis that Al 2 O 3 adheres via a dislocation-mediated mechanism in AD, that Al 2 O 3 powder agglomerates adhere as individual constituent particles rather than collectively, and that, for this case study, LIPIT measurements were predictive of AD process parameters.

Al2O3↗
DOE OSTI2021

Observation of Fundamental Mechanisms in Compression-Induced Phase Transformations Using Ultrafast X-ray Diffraction

As theoretically hypothesized for several decades in group IV transition metals, we have discovered a dynamically stabilized body-centered cubic (bcc) intermediate state in Zr under uniaxial loading at sub-nanosecond timescales. Under ultrafast shock wave compression, rather than the transformation from α-Zr to the more disordered hex-3 equilibrium ω-Zr phase, in its place we find the formation of a previously unobserved nonequilibrium bcc metastable intermediate. We probe the compression-induced phase transition pathway in zirconium using time-resolved sub-picosecond x-ray diffraction analysis at the Linac Coherent Light Source. We also present molecular dynamics simulations using a potential derived from first-principles methods which independently predict this intermediate phase under ultrafast shock conditions. In contrast with experiments on longer timescale (> 10 ns) where the phase diagram alone is an adequate predictor of the crystalline structure of a material, our recent study highlights the importance of metastability and time dependence in the kinetics of phase transformations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Data-Driven Insights into the Structural Essence of Plasticity in High-Entropy Alloys

The heterogeneous mechanical response of a crystalline alloy with multiple principal elements was investigated using molecular dynamics simulations. The local configuration of the alloy in its quiescent state was characterized by the variables derived from the gyration tensor and the atomic electronegativity. A multivariate analysis identified the geometric and chemical factors that influenced the atomic packing variations. Further, upon straining, the non-affine displacement exhibited spatial heterogeneity. A statistical correlation was established between the local yield events and the specific features of the local configuration. Our findings, validated by the performance metrics analysis, provided a structural criterion for the instability mechanisms in high-entropy alloys (HEAs) and enhanced the understanding of their plasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Towards an Understanding of Oxidative Damage in an a-L-Arabinofuranosidase of Trichoderma reesei : a Molecular Dynamics Approach

Trichoderma reesei is a ‘workhorse’ fungus that produces glycosyl hydrolases (e.g., cellulases) at high titers for use in industrial bioprocessing. In this study, we focused on two side-chain hemicellulases, a-L-arabinofuranosidases and a-D-galactosidases, which are important for the treatment of lignocellulosic biomass, but susceptible to oxidative damage that can occur during industrial processing. The molecular details that render these enzymes inactive have not yet been identified. To approach this issue, we used proteomics to identify amino acid residues that were oxidized after a relevant oxidative treatment (Fenton reaction). These oxidative modifications were included in the 3D-protein structures, and using molecular dynamics simulations, we then studied the behaviors of non-modified and oxidized enzymes. These simulations showed significant alterations of the conformational stability of proteins when oxidized, as evidenced by changes in root mean square deviation (RMSD) and principal component analyses (PCA) trajectories. Likewise, enzyme-ligand interactions such as hydrogen bonds were greatly reduced in quantity and quality in the oxidized proteins. Finally, free energy landscape plots showed that there was a more rugged energy surface in the oxidized proteins, implying a less favorable reaction pathway. These results reveal the basis for loss of function in these carbohydrate active enzymes (CAZYs) in T. reesei.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Differences in Sb 2 Te 3 growth by pulsed laser and sputter deposition

High quality van der Waals chalcogenides are important for phase change data storage, thermoelectrics, and spintronics. Using a combination of statistical design of experiments and density functional theory, we clarify how out-of-equilibrium van der Waals epitaxial deposition methods can improve the crystal quality of Sb 2 Te 3 films. We compare films grown by radio frequency sputtering and pulsed laser deposition (PLD). The growth factors that influence the crystal quality for each method are different. For PLD grown films a thin amorphous Sb 2 Te 3 seed layer most significantly influences the crystal quality. In contrast, the crystalline quality of films grown by sputtering is rather sensitive to the deposition temperature and less affected by the presence of a seed layer. This difference is somewhat surprising as both methods are out-of-thermal-equilibrium plasma-based methods. Non-adiabatic quantum molecular dynamics simulations show that this difference originates from the density of excited atoms in the plasma. The PLD plasma is more intense and with higher energy than that used in sputtering, and this increases the electronic temperature of the deposited atoms, which concomitantly increases the adatom diffusion lengths in PLD. In contrast, the adatom diffusivity is dominated by the thermal temperature for sputter grown films. These results explain the wide range of Sb 2 Te 3 and superlattice crystal qualities observed in the literature. Finally, these results indicate that, contrary to popular belief, plasma-based deposition methods are suitable for growing high quality crystalline chalcogenides.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Short- and medium-range orders in Al90Tb10 glass and their relation to the structures of competing crystalline phases

Molecular dynamics simulations using an interatomic potential developed by artificial neural network deep machine learning are performed to study the local structural order in Al90Tb10 metallic glass. We show that more than 80% of the Tb-centered clusters in Al90Tb10 glass have short-range order (SRO) with their 17 first coordination shell atoms stacked in a ‘3661’ or ‘15551’ sequence. Medium-range order (MRO) in Bergman-type packing extended out to the second and third coordination shells is also clearly observed. Analysis of the network formed by the ‘3661’ and ‘15551’ clusters show that ~82% of such SRO units share their faces or vertexes, while only ~6% of neighboring SRO pairs are interpenetrating. Such a network topology is consistent with the Bergman-type MRO around the Tb-centers. Moreover, crystal structure searches using genetic algorithm and the neural network interatomic potential reveal several low-energy metastable crystalline structures in the composition range close to Al 90 Tb 10 . Some of these crystalline structures have the ‘3661’ SRO while others have the ‘15551’ SRO. While the crystalline structures with the ‘3661’ SRO also exhibit the MRO very similar to that observed in the glass, the ones with the ‘15551’ SRO have very different atomic packing in the second and third shells around the Tb centers from that of the Bergman-type MRO observed in the glassy phase.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Segregation competition and complexion coexistence within a polycrystalline grain boundary network

Interfacial segregation can stabilize grain structures and even lead to grain boundary complexion transitions. However, understanding of the complexity of such phenomena in polycrystalline materials is limited, as most studies focus on bicrystal geometries. In this work, we investigate interfacial segregation and subsequent complexion transitions in polycrystalline Cu-Zr alloys using hybrid Monte Carlo/molecular dynamics simulations. No significant change in the grain size or structure is observed upon Zr dopant addition to a pure Cu polycrystal at moderate temperature, where grain boundary segregation is the dominant behavior. Segregation within the boundary network is inhomogeneous, with some boundaries having local concentrations that are an order of magnitude larger than the global value and others having almost no segregation, and changes to physical parameters such as boundary free volume and energy are found to correlate with dopant concentration. Further, another alloy sample is investigated at a higher temperature to probe the occurrence of widespread transitions in interfacial structure, where a significant fraction of the originally ordered boundaries transition to amorphous complexions, demonstrating the coexistence of multiple complexion types, each with their own distribution of boundary chemical composition. Overall, this work highlights that interfacial segregation and complexion structure can be diverse in a polycrystalline network. Here, the findings shown here complement existing computational and experimental studies of individual interfaces and help pave the way for unraveling the complexity of interfacial structure in realistic microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hall–Petch strengthening limit through partially active segregation in nanocrystalline Ag-Cu alloys

The breakdown from grain-size strengthening to softening mechanisms is generally well understood for high-purity nanocrystalline materials when the mean grain size decreases to the nanometer range. In nanocrystalline alloys, however, the stabilization of nanosized grains by grain-boundary solute segregation complicates the above mechanisms. Moreover, current segregation models have little predictive power for determining the optimal solute content that maximizes Hall-Petch strengthening effects. In this article, using large-scale hybrid Monte-Carlo/molecular dynamic simulations, we present a systematic study of the Hall–Petch breakdown in Cu-segregated Ag alloys with grain sizes ranging from 8 nm to 59 nm, where three concentration-dependent regimes of plasticity are described: (1) Classical segregation strengthening behavior at low solute contents, (2) shear band-induced softening at high solute contents, and (3) a previously unknown, but extended plateau of maximum strengths for intermediate solute contents from 4 to 15 at.%, which we term as nanocrystalline Sterling silver. We find that flow strengths in nanocrystalline Sterling alloys naturally exhibit a zero-slope limit at the smallest grain sizes that is well below the ideal Hall–Petch strengthening trend. This phenomenon results from partially active grain-boundary segregation that acts to influence interfacial plasticity in some, but not all, grain boundary regions. Our findings amplify the atomic nature of solute segregation and interaction at grain boundaries and its complex roles on grain boundary-mediated plasticity mechanisms in nanocrystalline alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nanotwin assisted reversible formation of low angle grain boundary upon reciprocating shear load

Severe plastic deformation of metals is known to lead to superior properties that cannot be achieved by any traditional metallurgic process. Origin of the superior properties is perceived to be closely associated with grain refinement, a fundamental process during the severe plastic deformation, which is essentially the formation of new grain boundaries. However, the atomistic mechanism of grain boundary formation remains largely obscure. Here, by using in-situ transmission electron microscopy and molecular dynamic simulation, we reveal, for the first time at atomic level, the shear-induced low-angle grain boundary (LAGB) formation processes in Au nanocrystal. In this work, we discover the LAGB formation is accomplished through inward propagation of nanotwins accompanied by dislocations gliding on the twin boundaries, a nanotwin-mediated dislocation slip mechanism, which shows reversible characteristic under reciprocating shear load and is affected by the nanocrystal microstructure and orientation. Our result unveils unprecedented atomistic insights on shear driven grain refinement towards nanostructure of superior properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic modeling of surface and grain boundary dislocation nucleation in FCC metals

Dislocation nucleation plays a critical role in the plastic deformation of crystalline materials. However, it is challenging to predict the active mode and associated rate of dislocation nucleation under typical experimental loading conditions through molecular dynamics simulation due to timescale limitations. In this report we use the free-end nudged elastic band method to determine the activation energies and activation volumes of dislocation nucleation in four typical face-centered cubic metals of Au, Al, Cu and Ni. We focus on the representative processes of surface and grain boundary dislocation nucleation. The atomistically determined activation volumes of these dislocation nucleation processes are larger than 10b3 (with b being the Burgers vector length) under typical experimental loading conditions. These results are compared with experimentally measured activation volumes in ultrafine-grained and nanocrystalline metals, thereby providing mechanistic insight into their rate-controlling deformation mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mechanical properties, failure mechanisms, and scaling laws of bicontinuous nanoporous metallic glasses

Molecular dynamics simulations are employed to study the mechanical properties of nanoporous Cu x Zr 1-x metallic glasses (MGs) with five different compositions, x = 0.28, 0.36, 0.50, 0.64, and 0.72, and porosity in the range 0.1 < Φ < 0.7. Results from tensile loading simulations indicate a strong dependence of Young's modulus, E, and Ultimate Tensile Strength (UTS) on porosity and composition. By increasing the porosity from Φ = 0.1 to Φ = 0.7, the topology of the nanoporous MG shifts from closed cell to open-cell bicontinuous. The change in nanoporous topology enables a brittle-to-ductile transition in deformation and failure mechanisms from a single critical shear band to necking and rupture of ligaments. Genetic Programming (GP) is employed to find scaling laws for E and UTS as a function of porosity and composition. Here, a comparison of the GP-derived scaling laws against existing relationships shows that the GP method is able to uncover expressions that can predict accurately both the values of E and UTS in the whole range of porosity and compositions considered.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Detecting Thermally-Induced Spinodal Decomposition with Picosecond Ultrasonics in Cast Austenitic Stainless Steels

Destructive techniques to monitor nuclear reactor component health may not always be available during service, as they are time-consuming and often require pre-installed inspection coupons. Non-destructive evaluation (NDE) techniques can bridge this gap by rapidly identifying the state of mission-critical reactor components, via inference between NDE-measurable material properties and those of ultimate interest, such as ductility and toughness. Here, we demonstrate one such inference about the health of thermally aged cast austenitic stainless steels. Observations of surface acoustic wave peak (SAW) splitting correlate with spinodal decomposition-induced embrittlement as destructively measured by Charpy impact energy. Elastodynamic calculations and molecular dynamics simulations of the effects of spinodal decomposition on elastic moduli support that the new acoustic modes present are due to stiffening in the δ-ferrite domains. Finally, this discovery enables one to probe structure-property relationships in materials in a greatly accelerated manner, suggesting that similar inference methods can be used to determine material fitness-for-service, or to quickly uncover new structure-property relationships.

304-type stainless steel↗
DOE OSTI2023

Athermal annealing of pre-existing defects in crystalline silicon

Systematic investigations of electronic energy loss (S e ) effects on pre-existing defects in crystalline silicon (Si) are crucial to provide reliance on the use of ionizing irradiation to anneal pre-existing defects, leading to successful implementation of this technology in the fabrication of Si-based devices. In this regard, the S e effects on nonequilibrium defect evolution in pre-damaged Si single crystals at 300 K has been investigated using intermediate-energy ions (12 MeV O and Si ions) that interact with the pre-damaged surface layers of Si mainly by ionization, except at the end of their range where the nuclear energy loss (S n ) is no longer negligible. Furthermore, under these irradiation conditions, experimental results and molecular dynamics simulations have revealed that pre-existing disorder in Si can be almost fully annealed by subsequent irradiation with intermediate-energy incident ions with S e values as low as 1.5–3.0 keV/nm. Selective annealing of pre-existing defect levels in Si at room temperature can be considered as an effective strategy to mediate the transient enhanced diffusion of dopants in Si. This approach is more desirable than the regular thermal annealing, which is not compatible with the processing requirements that fall below the typical thermal budget.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mobility of twinning dislocations in copper up to supersonic speeds

Understanding the mobility of twinning dislocations is important for multiscale modeling of crystal plasticity, especially at high strain rates, where such dislocations may reach transonic or supersonic speeds. Here, we used molecular dynamics simulations to investigate the relationship between dislocation velocity and the applied resolved shear stress of an edge twinning dislocation in copper up to supersonic speeds. The twinning dislocation mobility relation is composed of two branches separated by a band of forbidden velocities. The lower velocity branch is limited by the first transverse sound speed ~2000 m/s while the upper branch stretches from ~3500 m/s in the transonic regime to supersonic velocities. Twinning dislocations cannot undergo uniform steady-state motion at velocities within the forbidden band. Our simulation results also reveal that edge twinning dislocation motion in copper is kink-mediated. We discuss the implications of our findings for the motion of twins, twinning dislocations, and twinning dislocation kinks in copper.

36 MATERIALS SCIENCE↗