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,729 records · Page 96

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↗

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↗

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↗

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↗

Mitigating keyhole pore formation by nanoparticles during laser powder bed fusion additive manufacturing

Keyhole pore formation is one of the most detrimental subsurface defects in the laser metal additive manufacturing process. However, effective ways to mitigate keyhole pore formation beyond tuning laser processing conditions during keyhole mode laser melting are still lacking. Here we report a novel approach to mitigate keyhole pore formation during laser powder bed fusion (LPBF) process by using stable nanoparticles. The critical keyhole depth for keyhole pore generation (i.e., the largest keyhole depth without keyhole pore formation) during LPBF of Al6061 increases from 246 µm to 454 µm (85% increase) after adding TiC nanoparticles. In-depth x-ray imaging studies and thermo-fluid dynamics simulation enable us to identify that two mechanisms work together to mitigate keyhole pore generation: (1) adding nanoparticles prevents the keyhole from collapsing by increasing the liquid viscosity to impede the protrusion development; (2) adding nanoparticles slows down the keyhole pore movement by increasing the liquid viscosity, resulting in the recapturing of the pore by the keyhole. We further demonstrate that adding TiC nanoparticles can also eliminate the keyhole fluctuation induced keyhole pore during LPBF of Al6061. Our research provides a potential way to mitigate keyhole pore formation for defect lean metal additive manufacturing.

36 MATERIALS SCIENCE↗

Which way does the dendrite grow? Competition among epitaxy, preferred growth direction, and thermal gradients in powder bed fusion additive manufacturing

The as-processed microstructure of metal alloy parts manufactured through laser powder bed fusion (LPBF) is heavily derived from the cellular dendritic solidification. The growth direction of dendrites within the melt pool is determined through competition among epitaxial growth, preferred growth directions, and maximum thermal gradients. However, the dominant factor and the specific role of each in developing melt pool microstructures remain unknown. Here, in this study, we performed single laser track scans on an SS316L single crystal substrate and combined experimental characterization of microstructure and crystal orientations with Computational Fluid Dynamics simulations of thermal gradients to evaluate the role of each factor in determining dendritic growth direction and evolution. Our results reveal that epitaxial growth dominates microstructure development by preferentially growing along a single 〈100〉 variant of the single crystal substrate adjacent to the melt pool boundary. Under LPBF’s highly curved and rapidly evolving thermal field, this preferential dendrite variant selection and its continued growth from the melt pool boundary to the centerline are governed by the local temperature gradient magnitude at the solid-liquid interface, rather than by the instantaneous maximum temperature gradient direction alone. Using these findings, we successfully predict changes in the dendrite growth direction with changing laser scan direction on a single crystal substrate, and show that the geometric melt pool centerline can deviate from the microstructural centerline because asymmetric local temperature gradient magnitudes transiently limit growth, resulting in different dendrite travel distances on each side of the melt pool.

36 MATERIALS SCIENCE↗

Neutron thermalization in nuclear graphite: A modern story of a classic moderator

Graphite has been of great interest since the dawn of the nuclear era. In this work, we will focus on the thermal neutron scattering of two types of nuclear graphite that have different microstructures, namely, the historic PGA which has a very anisotropic medium-grained needle coke, and the new G347A which is a near-isotropic fine-grained graphite. We report the results of room temperature inelastic neutron scattering measurements for the scattering function, double differential scattering cross-section, and phonon density of states. The obtained results are discussed in conjunction with results extracted from first-principles lattice dynamics simulations of crystalline graphite, using the incoherent approximation and by also including the coherent-one phonon contribution. Excellent agreements have been found between the measured and calculated scattering functions and phonon densities of states. The inclusion of the coherent-one phonon scattering shows a significant improvement of the calculated scattering functions, double differentials, and inelastic scattering cross-sections. We also highlight the physical concerns about the current thermal neutron scattering ENDF/B-VIII of nuclear graphite library and its corresponding phonon density of states. This work aims to provide a new presentation of neutron thermalization by graphite and offers a new set of high quality experimental and computational data for consideration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Macroscale Control of Reactivity using 3D Printed Materials with Intrinsic Catalytic Properties

The morphology of heterogeneous catalysts can impact their performance. However, standard manufacturing methods like extrusion or pelleting offer little options for tailoring catalyst shape. Herein, stereolithographic 3D printing is used to produce catalysts with controlled topologies to enhance their performance. A series of magnetic stir-bar compartments (SBC) were 3D printed and tested as catalysts for sucrose hydrolysis. The SBC were printed using acrylic acid (AA) and 1,6-hexanediol diacrylate (HDDA) as acid sites and hydrophobic crosslinking domains, respectively. Variations in the number and tilt direction of the SBC blades produced significant changes in their apparent catalytic activities. These changes resulted from differences in the fraction of active surface effectively interacting with the reactants in solution, as revealed by computational fluid dynamics simulations. Moreover, varying HDDA:AA ratios in SBC regulated reactant-surface interactions to control catalytic activity. Overall, 3D printing catalysts enables quick performance optimization by simultaneously controlling macroscopic structure and molecular composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring NaCl-PuCl 3 molten salts with machine learning interatomic potentials and graph theory

Actinide molten salts are the basis of the liquid fuels used in molten salt reactors. Due to the inherent difficulties associated with high temperature and hazardous conditions, experimental investigations of fundamental properties of these materials are usually challenging. In this work, we describe the structure and transport of NaCl-PuCl 3 mixtures using computational techniques. Three compositions were considered (16, 25, and 36 mol% PuCl 3 ) over a temperature range (730 – 1257K) using ab initio molecular dynamics, which provided the necessary data sets for training machine learned interatomic potentials. Further, molecular dynamics simulations based on these potentials were then used to determine structure and transport properties. A substantial change was noted in the structure factor when increasing the PuCl 3 content from 25 to 36 mol%. This change is linked to the aggregation of larger Pu 3+ clusters. In addition, the similarity of the atomic environments of metal cations in molten salt systems to their solid states counterparts was investigated using an unsupervised learning technique. Finally, graph theory was employed to explore the structure and size of actinide networks. Consistent with the structure factor, a dense Pu 3+ intermolecular structure is observed within the 36 mol% PuCl 3 mixture. The structure of cation-cation inter-junctions is also discussed. In all cases, the diffusion of Pu 3+ is significantly lower than that of Na + and Cl - .

36 MATERIALS SCIENCE↗

Deformation behavior and phase transformation of nanotwinned Al/Ti multilayers

We report nanotwinned Al/Ti multilayers have exhibited size-dependent microstructure evolution and high strength. However, their deformation mechanisms are less well understood. In this work, we investigated the deformation mechanisms of nanotwinned Al/Ti multilayers with FCC/HCP layer interfaces by using in situ micropillar compression tests. Nanotwinned Al/Ti multilayers exhibit compressive strength up to 2.4 GPa and good work hardening capability. Post-compression TEM analyses reveal high-density stacking faults and the HCP-to-FCC phase transformations in Ti. Molecular dynamics simulations elucidate the mechanisms of deformation induced phase transformation in Ti and the influence of collective movement of partial dislocations on the deformability of Al/Ti multilayers.

36 MATERIALS SCIENCE↗

Traction-separation response of bilayer graphene interfaces: The role of an intercalated single atomic layer of water molecules and hydroxyl groups

Molecular dynamics simulations of the normal and shear traction-separation response of dry and wet bilayer graphene interfaces were performed. The influence of a single, atomic layer of water molecules sandwiched between graphene layers and hydroxyl functionalization on one or both layers were investigated. Hydroxyl groups provided a shielding effect during normal separation while favoring stick-slip friction under sliding conditions. The single, atomic layer of water molecules confined between hydrophilic and asymmetric hydrophobic/hydrophilic surfaces was found to ease sliding friction. The mechanisms arose from the balance between the attractive and repulsive graphene/water interactions through bridges of single water molecules between the two surfaces and the structure of the water monolayer under confinement for which an out-of-plane water/water hydrogen bond network could not be developed. The single, atomic layer of water molecules reduced the sliding friction the most when both surfaces were hydrophilic and lowered the sliding friction of the asymmetric hydrophobic/hydrophilic interfaces to that when both surfaces were hydrophobic. An optimal hydroxyl surface coverage of less than 10% may exist for efficient water lubrication, with diminishing effect with increasing number of hydroxyl groups. Finally, the findings offer a design strategy for assembling bilayer graphene interfaces with tailored friction and sliding contact properties.

36 MATERIALS SCIENCE↗

Formation of pyrophosphates across grain boundaries induces the formation of mismatched but oriented interfaces in silver phosphate polypods

The interfaces and their misfit defects determine the materials properties for a wide range of applications, such as electronic devices, photocatalysis, and mechanical engineering, etc. However, current understanding of atomic interfacial structures is limited. Here we discover a special interfacial structure, mismatched but oriented interface via two distinct facets. Transmission and scanning electron microscopy results suggest that, in Ag 3 PO 4 polypods structures, interfaces of {1 0 0} and {1 1 0}, ({1 0 0}/(1 1 0}), {1 0 0}/{1 1 1}, {1 1 0}/{1 1 1}, and {1 0 0}/{1 0 0}, etc., have a certain orientation relationship, corresponding to the energy minima and coincident site lattice of interfacial atoms as demonstrated by molecular dynamics simulations. Density functional theory demonstrates that the formation of pyrophosphate and/or phosphates rotation to bond across the interface compensate the lattice mismatch at the interfaces, as well as deformations of Ag-O bonds. Furthermore, our work opens up a new avenue for a much wider range of interfacial structures, allow for a higher diversity of structures, and shine light on tailoring crystal structures, morphologies, and the resulting properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal energy evolution and mechanical deformation of monocrystalline yttria-stabilized zirconia nanoparticles in aerosol deposition processes

Aerosol deposition (AD) is a coating process wherein aerosol particles are impacted on a target substrate. There are fundamental differences between the AD process (cold impact), where particle translational kinetic energy is high and thermal energy is low, and thermal spray deposition (thermal impact), where translational energy is lower but thermal energy is high. To better compare cold and thermal impact effects on particles, we carried out molecular dynamics simulations using yttria-stabilized zirconia (YSZ) nanoparticles on YSZ substrates as a model system. We performed cold impact simulations at 300 K with variable impact velocity in the 500 ms –1 –1500 ms –1 range to understand how increasing translational kinetic energy affects thermal energy and mechanical evolution. We then performed thermal impact simulations at variable temperature and impact velocity, but where the total kinetic energy of the nanoparticle was equivalent to that of a 300 K, 1000 impact. In cold impact, the temperature increases in YSZ nanoparticles at a rate of 10 13 – 10 14 K s –1 , and large temperature gradients result. Conversely, in thermal impact, nanoparticle temperatures remain uniform. Furthermore, the temperature gradients during cold impact coincide with plastic deformation in nanoparticles, while with larger thermal energies, plastic deformation is reduced.

42 ENGINEERING↗

The relation between lignin sequence and its 3D structure

Background: Lignin, the second most abundant biopolymer on earth, plays a major structural role in plants, conferring mechanical strength and regulating water conduction. Understanding the three-dimensional structure of lignin is important for fundamental reasons as well as engineering plants towards lignin valorization. Lignin lacks a specific primary sequence, making its average chemical composition the focus of most recent studies. However, it remains unclear whether the 3D structure of lignin molecules depends on their sequence. Methods: In this work, we performed all-atom molecular dynamics simulation of three S/G-lignin molecules with the same average composition but different sequence. Results: A detailed statistical analysis of the radius of gyration and relative shape anisotropy reveals that the lignin sequence has no statistically significant effect on the global three-dimensional structure. We found however, that homopolymers of C-lignin with the same molecular weight have smaller radii of gyration than S/G-lignin. We attribute this to lower hydroxyl content of C-lignin, which makes it more compact and rigid. Conclusions: The 3D structure of lignin is influenced by the overall content of monomeric units and interunit linkages and not by its precise primary sequence. General Significance: Lignin is assumed to not have a well-defined primary structure. The results presented here demonstrate there are no significant differences in the global 3D structure of lignin molecules with the same average composition but different primary sequence.

59 BASIC BIOLOGICAL SCIENCES↗

Conformational space exploration of cryo-EM structures by variability refinement

Cryo-EM observation of biological samples enables visualization of sample heterogeneity, in the form of discrete states that are separable, or continuous heterogeneity as a result of local protein motion before flash freezing. Variability analysis of this continuous heterogeneity describes the variance between a particle stack and a volume, and results in a map series describing the various steps undertaken by the sample in the particle stack. While this observation is absolutely stunning, it is very hard to pinpoint structural details to elements of the maps. Here, in order to bridge the gap between observation and explanation, we designed a tool that refines an ensemble of structures into all the maps from variability analysis. Using this bundle of structures, it is easy to spot variable parts of the structure, as well as the parts that are not moving. Comparison with molecular dynamics simulations highlights the fact that the movements follow the same directions, albeit with different amplitudes. Ligand can also be investigated using this method. Variability refinement is available in the Phenix software suite, accessible under the program name phenix.varref.

59 BASIC BIOLOGICAL SCIENCES↗

Lpp positions peptidoglycan at the AcrA-TolC interface in the AcrAB-TolC multidrug efflux pump

The multidrug efflux pumps of Gram-negative bacteria are a class of complexes that span the periplasm, coupling both the inner and outer membranes to expel toxic molecules. The best-characterized example of these tripartite pumps is the AcrAB-TolC complex of Escherichia coli. However, how the complex interacts with the peptidoglycan (PG) cell wall, which is anchored to the outer membrane (OM) by Braun’s lipoprotein (Lpp), is still largely unknown. In this work, we present molecular dynamics simulations of a complete, atomistic model of the AcrAB-TolC complex with the inner membrane, OM, and PG layers all present. We find that the PG localizes to the junction of AcrA and TolC, in agreement with recent cryo-tomography data. Free-energy calculations reveal that the positioning of PG is determined by the length and conformation of multiple Lpp copies anchoring it to the OM. The distance between the PG and OM measured in cryo-electron microscopy images of wild-type E. coli also agrees with the simulation-derived spacing. Finally, sequence analysis of AcrA suggests a conserved role for interactions with PG in the assembly and stabilization of efflux pumps, one that may extend to other trans-envelope complexes as well.

59 BASIC BIOLOGICAL SCIENCES↗

Steric effects of central dogma processes on the compaction and segregation of bacterial nucleoids

The bacterial cytoplasm is characterized by a distinctive membrane-less organelle, the nucleoid, which harbors the chromosomal DNA. Here, we investigate the steric effects of dynamic processes associated with transcription and translation on the structure of this organelle using coarse-grained molecular dynamics simulations that incorporate out-of-equilibrium reactions. Our model captures the scale of the entire cell and incorporates a reaction-diffusion system for ribosomes and polyribosomes, coupling their nonequilibrium kinetics to DNA excluded volume interactions. Our findings demonstrate that out-of-equilibrium reactions increase the size of the nucleoid and the number of ribosomes in this subcellular region. In addition, we show that nucleoid size is proportional to transcriptional activity. Our model reproduces the time-dependent change in nucleoid size observed in rifampicin treatment experiments, where the pool of polyribosomes is depleted. Furthermore, we find that these active processes are essential for complete sister chromosome separation and correct nucleoid positioning within the cell. Overall, our study reveals the effects of the central dogma processes on the internal organization and localization of bacterial nucleoids.

Chang, Mu-Hung [Univ. of Tennessee, Knoxville, TN ↗

Converting PBO fibers into carbon fibers by ultrafast carbonization

Here, we report an ultrafast transformation of poly(p-phenylene-2,6-benzobisoxazole) (PBO) fibers into car- bon fibers (CFs) by direct carbonization that eliminates high-cost oxidation in traditional CF conversion. Ultrafast heating and cooling render PBO-derived CFs superlative mechanical properties. ReaxFF atomistic-scale reactive molecular dynamics simulations unveil the PBO to CF conversion mechanism in which ultrafast heating suppresses the O-contained gas release, facilitating all-carbon rings alignment during the carbonization of PBO. These results offer new guidelines for identifying and evaluating alternative CF precursors and fresh perspectives on smart manufacturing of CFs.

42 ENGINEERING↗