Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials dynamics”

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 271 records · Page 15

Ultrafast X-ray Diffraction Study of a Shock-Compressed Iron Meteorite above 100 GPa

Natural kamacite samples (Fe92.5Ni7.5) from a fragment of the Gibeon meteorite were studied as a proxy material for terrestrial cores to examine phase transition kinetics under shock compression for a range of different pressures up to 140 GPa. In situ time-resolved X-ray diffraction (XRD) data were collected of a body-centered cubic (bcc) kamacite section that transforms to the high-pressure hexagonal close-packed (hcp) phase with sub-nanosecond temporal resolution. The coarse-grained crystal of kamacite rapidly transformed to highly oriented crystallites of the hcp phase at maximum compression. The hcp phase persisted for as long as 9.5 ns following shock release. Comparing the c/a ratio with previous static and dynamic work on Fe and Fe-rich Fe-Ni alloys, it was found that some shots exhibit a larger than ideal c/a ratio, up to nearly 1.65. This work represents the first time-resolved laser shock compression structural study of a natural iron meteorite, relevant for understanding the dynamic material properties of metallic planetary bodies during impact events and Earth’s core elasticity.

36 MATERIALS SCIENCE↗

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗

Constitutive Behavior and Processing Map of T2 Pure Copper Deformed from 293 to 1073 K

The deformation behavior of T2 pure copper compressed from 293 to 1073 K with strain rates from 0.01 to 10 s{sup −1} was investigated. The constitutive equations were established by the Arrhenius constitutive model, which can be expressed as a piecewise function of temperature with two sections, in the ranges 293-723 K and 723-1073 K. The processing maps were established according to the dynamic material model for strains of 0.2, 0.4, 0.6, and 0.8, and the optimal processing parameters of T2 copper were determined accordingly. In order to obtain a better understanding of the deformation behavior, the microstructures of the compressed samples were studied by electron back-scattered diffraction. The grains tend to be more refined with decreases in temperature and increases in strain rate.

36 MATERIALS SCIENCE↗

A New Method for Optimizing Hot Processing Parameters of Mg-6.0Zn-0.5Mn-0.5Er Alloy Based on Kinetic Model of Dynamic Recrystallization and Processing Map

Hot deformation behavior of the cast-homogenized Mg-6.0Zn-0.5Mn-0.5Er alloy was studied using dynamic recrystallization (DRX) kinetic model and processing map. The compressing tests were conducted in the temperature range of 250-450 °C and strain rates of 0.001-1 s{sup −1}. According to the evolution of microstructures, under lower strain rates, the main DRX mechanism of Mg-6.0Zn-0.5Mn-0.5Er alloy is continuous DRX (CDRX); twinning induced DRX (TDRX) and CDRX both become the main DRX mechanism under higher strain rates. The DRX kinetics model of Mg-6.0Zn-0.5Mn-0.5Er alloy is calculated as X{sub DRX}=1−exp[−1.9463(((ε−ε{sub c})/(ε{sup ∗}))){sup 1.4608}], which is corresponding to the microstructure evolutions of DRX under different deformation conditions. The contour map of DRX was proposed based on the calculation results of DRX kinetics model. The processing maps are constructed to predict processing parameters of the alloy, and the predictability was evaluated combining with the contour map of DRX and dynamic materials model (DMM) processing map. It is deduced from the microstructures evolution and processing map that the optimum processing domain is mainly at 380-450 °C and 0.01-0.001 s{sup −1}, and 410-420 °C and 0.01-1 s{sup −1}.

36 MATERIALS SCIENCE↗

Hot Deformation Behavior and Processing Map of a Fe-11Mn-10Al-0.9C Duplex Low-Density Steel Susceptible to κ-Carbides

The hot deformation behavior of a Fe-11Mn-10Al-0.9C duplex low-density steel has been investigated, based on a series of isothermal compression tests at temperatures of 800-1100 °C and strain rates of 0.001-10 s{sup −1}. At relatively lower temperatures, with austenite decomposing, continuous dynamic recrystallization of ferrite and the formation of inter-granular κ-carbide were responsible for significant softening. At high deformation temperature where austenite was a dominant phase, dynamic recrystallization (DRX) is the restoration mechanism. Processing maps were developed at different plastic strains employing dynamic materials model and further verified through microstructural characterization. The optimal hot deformation condition at a large strain (0.7) was identified as deformation temperature of 950-1100 °C and strain rate of 0.01-1.0 s{sup −1}. In this domain, the original coarse grains could be replaced by the fine and uniform recrystallized grains, indicating that the high efficiency was dissipated by DRX. Meanwhile, two instable regions resulted from the formation of inter-granular κ-carbides and necklace structure should be avoided during hot processing.

36 MATERIALS SCIENCE↗

A Systematic Assessment of Various Stability/Instability Criteria in Predicting the Hot Deformation-Related Instabilities in Super-304H Stainless Steels

In this study, different existing stability/instability criteria (viz Prasad, Murty, Gegel, Alexander, Jonas, and Semiatin) are used to characterize the hot deformation instabilities of super-304H stainless steel. The instability maps are developed based on flow stress data acquired from uniaxial hot compression tests performed at various temperatures (1173-1423 K) and strain rates (0.001-10 s{sup −1}). The deformation instabilities are further validated through detailed microstructural observations. Flow localization and cracking have been identified as the signature of deformation instabilities in this alloy. The microcracks have been found to be associated with the grain boundary precipitation of coarse (~ 1- to 5-µm-sized) Nb-rich particles (mainly NbC). The instability criteria developed based on dynamic material modeling (i.e., Prasad’s and Murty’s criteria) could accurately predict the formation of microcracks in the entire hot deformation range employed in this study. Albeit the Gegel’s and Alexander’s stability criteria could correctly predict the flow localization, it has overestimated the deformation instabilities at high temperatures (> 1250 K) and low strain rates (≤ 0.1 s{sup −1}). In contrast, Jonas’ and Semiatin’s instability criteria have under-predicted the unstable domains in the studied alloy.

36 MATERIALS SCIENCE↗

Constitutive Behavior and Hot Workability of a Hot Isostatic Pressed Ti-22Al-25Nb Alloy during Hot Compression

A Ti-22Al-25Nb alloy was fabricated from prealloyed powders by hot isostatic pressing for 2 h at a temperature of 1050 °C and pressure of 100 MPa. The hot deformation behavior of the Ti-22Al-25Nb alloy was characterized by isothermal compression testing at deformation temperatures between 900 and 1060 °C and strain rates between 0.001 and 1 s{sup −1}. Based on the true stress–strain curves, a constitutive equation was constructed to describe the flow stress as a function of the strain rate and deformation temperature. Three-dimensional (3D) processing maps were developed based on dynamic material model theory using the stress flow data to identify the instability and optimization regions of the hot processing parameters. The results show that the flow stress decreases with increasing deformation temperature and decreasing strain rate, and the softening mechanisms are different under different deformation conditions. The apparent activation energies in the (α{sub 2} + β/B2 + O) and (α{sub 2} + B2) phase regions are calculated to be 865.177 and 590.661 kJ/mol, respectively, which suggests that the hot isostatic pressed (HIPed) Ti-22Al-25Nb alloy requires a high hot deformation activation energy. Combined with microstructural observations, the optimal processing domains are determined to be a temperature range of 1030-1060 °C and strain rate range of 0.01-0.1 s{sup −1} in the α{sub 2} + B2 phase region. Moreover, the results indicate that adiabatic shear bands and severe inhomogeneous deformation cause flow instability at lower temperatures (900-1005 °C) and higher strain rates (> 0.1 s{sup −1}).

36 MATERIALS SCIENCE↗

Assessing the role of surface carbon on the surface stability and reactivity of β-Mo 2 C catalysts

Transition-metal carbides, e.g., Mo 2 C, promote industrially important catalytic reactions with activities comparable to precious-metal catalysts. Yet, the nature of the catalytically active sites remains unclear. Herein, density functional theory calculations are applied to systematically assess the stability of low-Miller index facets of β-Mo 2 C as a function of carbon chemical potential (uC) and surface-carbon coverage (nC). When reconstruction of surface carbon is considered to characterize β-Mo 2 C under varying environments, the (111) and (101) surfaces are predicted to be most stable, dominating the Wulff particle surface area across a broad range of uC. The surface-carbon coverage is predicted to steadily increase for all low-index facets under more carburizing reaction conditions, leading to distinct effects on the adsorption energetics for atomic hydrogen and oxygen. Whereas H binding is only slightly affected by surface carbon, O is systematically destabilized. The findings reported indicate that consideration of only bulk-terminated surfaces, and the common focus of theoretical studies on the bulk-truncated, close-packed (100) surface, is likely insufficient to capture the reactivity of the working ..beta..-Mo 2 C catalyst. Similar to metal oxides, β-Mo2C catalysts are predicted to be dynamic materials with surface compositions and reactivities that evolve in response to the environment.

36 MATERIALS SCIENCE↗

A fourth-order Lagrangian discontinuous Galerkin method using a hierarchical orthogonal basis on curvilinear grids

The existing high-order Lagrangian discontinuous Galerkin (DG) hydrodynamic methods are restricted to using quadratic meshes with quadratic polynomials (P2), which in turn, yield up to third-order accuracy. Additionally, these existing DG hydrodynamic schemes, when extended to work with cubic meshes and cubic polynomials (P3), can be unstable on strong-shock problems. Therefore, this paper presents a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method to simulate compressible material dynamics ( e.g., gasses, fluids, and solids) with strong-shocks using cubic meshes and cubic polynomials, and delivers up to fourth-order accuracy on smooth flows. The stability on shock problems is achieved using new hierarchical orthogonal basis functions and a new subcell mesh stabilization (SMS) scheme for cubic meshes. The accuracy and robustness of the new high-order accurate Lagrangian DG hydrodynamic method is demonstrated by simulating a diverse suite of challenging test problems covering gas and solid dynamic problems on curvilinear meshes.

97 MATHEMATICS AND COMPUTING↗

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.

42 ENGINEERING↗

Fatigue crack growth mechanisms at the microstructure scale in as-fabricated and heat treated Ti-6Al-4V ELI manufactured by electron beam melting (EBM)

Electron beam melting (EBM) is a metal powder bed fusion additive manufacturing (AM) technology that fabricates parts by selectively scanning consecutive powder layers with an electron beam. Additive manufacturing technologies are increasing in importance for aerospace and medical applications, where the demand for a fundamental understanding and predictability of static and dynamic material properties are high. Ti-6Al-4V is the most widely used and studied alloy for this technology, and is the focus of this work in its Extra Low Interstitial (ELI) variation. The layered manufacturing of metallic components by EBM creates a unique directional microstructure, and consequently, anisotropic properties. Microstructure evolution, and its influence on mechanical properties of the alloy in the as-fabricated condition, has been documented by various researchers. However, fatigue crack propagation and the effects of the directional structure have not been sufficiently studied, imposing a barrier for this technology’s potential extension to high-integrity applications. Here in this study, fatigue crack growth (FCG) both parallel and perpendicular to the build directions was studied for different stress ratios and crack growth stages. The interaction between the directional as-fabricated EBM microstructure and FCG was investigated and compared to that of the equiaxed β annealed microstructure obtained by annealing above the β transus temperature. Finally, the FCG threshold, ΔK th , was analytically modeled for the two relative crack propagation directions at different stress ratios, and FCG microstructural mechanisms were established for all three regions of crack propagation.

36 MATERIALS SCIENCE↗

Strength of tantalum shocked at ultrahigh pressures

High purity polycrystalline tantalum (Ta) was shocked through 1–3.5 Mbar pressures creating Richtmyer-Meshkov unstable interfaces that were used to determine the dynamic material strength. The experiments were performed on the Omega laser at the University of Rochester Laboratory for Laser Energetics. Prior to shock, the driven surfaces of the tantalum targets where coined with a sinusoidal pattern. The targets were recovered post-shock, and the growth of the sinusoid amplitudes was used to characterize the relative extent of plastic deformation as a function of laser energy. Additionally, analogous data were extracted prior to the experiments from phenomenological Von Mises plasticity simulations that considered equation of state tables for Ta. The simulations showed the best agreement with the experiments (less than 5% difference between mean ripple growth measures) for shock pressures ranging from 1.2 to 2.7 Mbar. A fluid model studied as a function of viscosity was also used to qualitatively indicate the sensitivity of the experiments to strength. These results verify the ability to use a phenomenological, equation of state based model to simulate very high-strain rate, high-pressure deformation of tantalum.

36 MATERIALS SCIENCE↗

Circular economy pathways for decarbonizing aluminum and steel automotive body sheet components in the United States

Decarbonizing vehicle production is essential to reducing automotive sector emissions. This study quantifies greenhouse gas (GHG) emissions from aluminum and steel auto-body sheet components produced in the US. It evaluates the effectiveness of circular economy (CE) strategies (greater closed-loop recycling of pre-consumer scrap, post-consumer scrap, and increased manufacturing yields) to reduce supply chain emissions across different process technology and electricity grid decarbonization pathways. We combine dynamic material flow analysis (2025–2050) with cradle-to-gate life-cycle modeling to assess production emissions and the potential reductions associated with the CE strategies under frozen, moderate, and aggressive technology and grid decarbonization scenarios. Current emissions intensities are estimated at approximately 12.3 kg.CO₂eq/kg of aluminum and 4.3 kg.CO₂eq/kg of steel sheet embedded in the vehicle. Under the frozen decarbonization scenario and current levels of circularity, annual emissions attributable to US aluminum and steel auto-body sheet supply chains could rise by 54 % and 18 % respectively by 2050. Rapid deployment of the CE strategies can cut these annual emissions in 2050 by 52 % for aluminum and 44 % for steel. However, scrap quality constraints lead to saturation points, limiting these benefits unless addressed. Aggressive deployment of low-carbon production technologies and a low-carbon grid reduces the relative benefit of implementing the CE strategies; however, even under the aggressive technology and grid decarbonization scenario, the CE strategies reduce annual emissions by a further 23 %-54 % by 2050. These findings highlight the urgent need to integrate CE strategies into the sheet metal supply chain to support decarbonization efforts and help meet climate targets.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

ExactPack: A python library of exact analytic solutions

Verification of multi-physics simulation software against problems with known analytic or semi-analytic solutions is an important aspect of research into a wide variety of fields involving the motion of fluids, shock physics and other dynamic material properties. Previous work comparing simulation results against analytic solutions has been ad-hoc, with developers frequently writing their own analytic solvers. This has resulted in a large amount of duplicated effort. The python library ExactPack has been developed as a collection of analytic and semi-analytic solvers to a variety of multi-physics problems, providing a consistent API to a set of well-tested solver implementations.

97 MATHEMATICS AND COMPUTING↗

Protein Assembly by Design

Proteins are nature’s primary building blocks for the construction of sophisticated molecular machines and dynamic materials, ranging from protein complexes such as photosystem II and nitrogenase that drive biogeochemical cycles to cytoskeletal assemblies and muscle fibers for motion. Such natural systems have inspired extensive efforts in the rational design of artificial protein assemblies in the last two decades. As molecular building blocks, proteins are highly complex, in terms of both their three-dimensional structures and chemical compositions. To enable control over the self-assembly of such complex molecules, scientists have devised many creative strategies by combining tools and principles of experimental and computational biophysics, supramolecular chemistry, inorganic chemistry, materials science, and polymer chemistry, among others. Owing to these innovative strategies, what started as a purely structure-building exercise two decades ago has, in short order, led to artificial protein assemblies with unprecedented structures and functions and protein-based materials with unusual properties. Furthermore, our goal in this review is to give an overview of this exciting and highly interdisciplinary area of research, first outlining the design strategies and tools that have been devised for controlling protein self-assembly, then describing the diverse structures of artificial protein assemblies, and finally highlighting the emergent properties and functions of these assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Davis Computational Spectroscopy Workflow—From Structure to Spectra

Here, we describe an automated workflow that connects a series of atomic simulation tools to investigate the relationship between atomic structure, lattice dynamics, materials properties, and inelastic neutron scattering (INS) spectra. Starting from the atomic simulation environment (ASE) as an interface, we demonstrate the use of a selection of calculators, including density functional theory (DFT) and density functional tight binding (DFTB), to optimize the structures and calculate interatomic force constants. We present the use of our workflow to compute the phonon frequencies and eigenvectors, which are required to accurately simulate the INS spectra in crystalline solids like diamond and graphite as well as molecular solids like rubrene. We have also implemented a machine-learning force field based on Chebyshev polynomials called the Chebyshev interaction model for efficient simulation (ChIMES) to improve the accuracy of the DFTB simulations. We then explore the transferability of our DFTB/ChIMES models by comparing simulations derived from different training sets. We show that DFTB/ChIMES demonstrates ~100× reduction in computational expense while retaining most of the accuracy of DFT as well as yielding high accuracy for different materials outside of our training sets. The DFTB/ChIMES method within the workflow expands the possibilities to use simulations to accurately predict materials properties of increasingly complex structures that would be unfeasible with ab initio methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tracking the Dynamics of a Ag-MnO x Oxygen Reduction Catalyst Using In Situ and Operando X-ray Absorption Near-Edge Spectroscopy

Sustainable electricity generation via hydrogen fuel cells requires the development of efficient oxygen reduction reaction (ORR) catalysts. In situ/operando experiments are necessary to uncover the extent of dynamic material changes during catalysis. Herein, we use in situ/operando X-ray absorption near-edge spectroscopy to track Mn valence changes of a promising, ultrathin, porous MnO x layer on a Ag thin film. Mn–K-edge measurements as a function of electrochemical environment and ORR conditions reveal that, interestingly, when driving the ORR at 0.8 VRHE, the Mn is distinctively more reduced and the MnO x redox is nonreversible in contrast to measurements in N 2 -saturated electrolyte. Ex situ inductively coupled plasma mass spectrometry, atomic force microscopy, and X-ray photoelectron spectroscopy indicate that these phenomena do not correlate to metal dissolution but might be associated with morphological surface reconstruction related to Ag valency. Furthermore, this study highlights how the microenvironment and catalysis play a key role in the in situ/operando surface structure and chemical state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Schools of skyrmions with electrically tunable elastic interactions

Coexistence of order and fluidity in soft matter often mimics that in biology, allowing for complex dynamics and applications-like displays. In active soft matter, emergent order can arise because of such dynamics. Powered by local energy conversion, this behavior resembles motions in living systems, like schooling of fish. Similar dynamics at cellular levels drive biological processes and generate macroscopic work. Inanimate particles capable of such emergent behavior could power nanomachines, but most active systems have biological origins. Here we show that thousands-to-millions of topological solitons, dubbed “skyrmions”, while each converting macroscopically-supplied electric energy, exhibit collective motions along spontaneously-chosen directions uncorrelated with the direction of electric field. Within these “schools” of skyrmions, we uncover polar ordering, reconfigurable multi-skyrmion clustering and large-scale cohesion mediated by out-of-equilibrium elastic interactions. Remarkably, this behavior arises under conditions similar to those in liquid crystal displays and may enable dynamic materials with strong emergent electro-optic responses.

Science & Technology - Other Topics↗