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

Negative thermal expansion coefficient and amorphization in defective 4H-SiC

This paper presents thermal expansion coefficient (TEC) and amorphization in 4H-SiC containing point defects at different concentrations. We considered vacancy defects, interstitial defects, and Frenkel pair defects and investigated the thermomechanical response of the lattice over a wide range of temperatures using classical molecular dynamics simulations. The results show that 4H-SiC with vacancy defects exhibits a negative TEC above a critical defect density of around 9% (irrespective of the temperature). With interstitial defects, it exhibits a positive TEC (regardless of the defect density), and with Frenkel pair defects it shows a transition from positive TEC to negative TEC for a defect density greater than 8%. The coupling between temperature-induced expansion and defect-introduced stress in the lattice forms the mechanistic basis for the observed variation in TEC. Furthermore, the specific heat decreases rapidly with an increase in defect density at room temperature, with the highest sensitivity of the lattice observed for the Frenkel pair defects followed by interstitial defects and then by vacancy defects. Finally, these findings highlight the critical implications of defects on thermal expansion behavior of 4H-SiC with applications in radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spontaneous rearrangement of acetylated xylan on hydrophilic cellulose surfaces

The interaction of xylan, an abundant plant polysaccharide, with cellulose microfibrils is essential for secondary cell wall strength. A deeper understanding of these interactions is crucial both to improve our understanding of plant cell wall architecture and to design alternate strategies to overcome cellulose recalcitrance for the production of biofuels and sustainable biomaterials. Naturally occurring acetate or glucuronic acid substitutions on xylan have been shown to influence xylan-cellulose interactions. In this work, we use unrestrained molecular dynamics simulations to determine the interactions with the (110) hydrophilic face of cellulose fibers of four different xylans. In the absence of cellulose, all xylans, independent of the substitution pattern, adopt a highly flexible threefold helical screw conformation. However, when xylan is spatially close to a cellulose surface 1,2 linked acetyl xylans (2AcX) adopt rigid twofold helical screw conformations. The 2AcX conformations are primarily stabilized by interactions between the acetylated oxygen and the glycosidic linkage with C-O6 of cellulose. In contrast, the glycosidic oxygens and acetyl decorations for 1,3 linked acetyl groups (3AcX) are oriented away from the cellulose surface and the 3AcX xylans maintain threefold helical screw conformations on the cellulose surface. Our results show that evenly spaced chemical functionalization (with acetyl groups) and the position of substitution (1,2) on xylan backbone play key roles in tuning the xylan-cellulose interactions to stabilize the twofold helical screw conformations of xylan on the cellulose surface. A comparison with previous experimental findings further suggests that 1,2 substitutions induce twofold helical screw conformations of xylan on the cellulose surface irrespective of the chemical nature of the substituent, while 1,3 substitutions primarily bind lignin in threefold helical screw conformations rather than cellulose in plant cell walls.

1,2 Ac↗
DOE OSTI2023

Patterns in interactions of variably acetylated xylans with hydrophobic cellulose surfaces

The recalcitrance of plant cell wall lignocellulosic biomass to deconstruction is a major hurdle to sustainable biofuel/bioproduct economy. A multitude of interactions stabilize lignocellulosic biomass structure. Among these, tight packing of hemicellulose-cellulose is partly responsible for biomass recalcitrance. Here, unrestrained molecular dynamics simulations are employed to understand the influence of the nature and pattern of naturally-occuring acetyl decorations of the xylan backbone on interactions with the (100) hydrophobic cellulose surface. Periodically O2-acetylated xylan (2AcX) assume twofold helical screw conformations that are stabilized by a combination of multiple hydrophobic contacts and hydrogen bonds with the hydrophobic cellulose surface. In contrast, acetylation at the O3 position in xylan obstructs interactions, thereby adopting threefold helical screw conformations that potentially preferably interact with lignin rather than cellulose. Fully acetylated xylan desorbs from the surface implying a minimum number of unsubstituted residues on the xylan backbone is required for interaction with the surface. Here, the substituted residues must form ~ 20% fewer contacts than the unsubstituted residues to sustain stable twofold helical screw xylan conformations on the cellulose surface. Thus, specific roles of macromolecular conformations of cellulose and hemicellulose in influencing the supramolecular interactions and function of plant cell walls have been determined.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Integrated data-driven and experimental approaches to accelerate lead optimization targeting SARS-CoV- 2 main protease

Identification of potential therapeutic candidates can be expedited by integrating computational modeling with domain aware machine learning (ML) approaches followed by experimental validation. Generative deep learning models have been recently developed that can generate thousands of new candidates, but their physiochemical properties are typically not optimized. Using our deep learning models and a scaffold as a starting point, we generated tens of thousands of compounds for SARS-CoV-2 M pro that preserve the core scaffold. Here we utilized and implemented several computational tools such as structural alert and toxicity analysis, high throughput virtual screening, ML-based 3D quantitative structure–activity relationships, multi-parameter optimization, and graph neural networks on libraries of generated candidates to predict biological activity and binding affinity a priori. From these collective computational results, eight promising candidates were identified and tested experimentally using Native Mass Spectrometry (MS) and FRET-based functional assays. Two compounds, with quinazoline-2-thiol and acetylpiperidine core moiety showed IC 50 values in the low micromolar range: 2.95±0.0017 µM and 3.41±0.0015 µM, respectively. The molecular dynamics simulations further highlight that binding of these compounds results in allosteric modulations in the chain B and the interface domains of the M pro . The key fragments from these top hits can be used as input for closed loop lead optimization in the integrated pipeline.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

The heterogeneous nature of mechanically accelerated grain growth

Abstract While grain growth is traditionally viewed as a purely thermally driven process, nanocrystalline metals can undergo grain growth under mechanical loads, even at room temperature. We performed a detailed atomistic study of the heterogeneous nature of mechanically accelerated grain growth in a polycrystalline Pt nanowire. Using molecular dynamics simulations, we compared the grain-growth behavior of individual grains during tensile and shear cyclic loading, for three different equivalent strain levels, and at two temperatures. Pure thermal grain growth with no mechanical loading provided a baseline reference case. On average, grains that were already susceptible to thermal grain growth were stimulated to grow faster with mechanical loading, as expected. However, when analyzed on a grain-by-grain basis, the results were far more complex: grains that grew fastest under one stimuli were less accelerated under other stimuli. Even when the magnitude of loading changed, the relative growth of individual grains was distorted. We interpret this complexity from the perspective of superimposed growth mechanisms.

Chen, Elton Y.↗
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↗