Mesoscale modeling of shock waves in Ni+Al multilayers poses significant challenges that are due, in part, to shock-induced chemical reactions. Current modeling approaches utilize reactive molecular dynamics (MD), but they are limited to resolving domains of only a few hundred nanometers. In contrast, actual multilayer superlattices can be tens of micrometers thick, and they exhibit non-ideal (i.e., wavy) interfaces. The second part of our research builds upon previous work developing physically based, thermodynamically complete equations of state for various Ni and Al intermetallic compositions. Here, we introduce a novel workflow for high-fidelity mesoscale simulations of Ni+Al multilayers using a continuum hydrocode. By increasing the simulation domain size beyond MD limitations (e.g., 2 × 6 μm2) and incorporating explicit interfacial roughness, we investigate the shock response of Ni+Al multilayers at previously unexplored scales. Our experimental design encompasses nine multilayer geometries with varying roughness amplitudes and tilt angles (θ = 15°, 30°, and 45°), alongside 19 flyer impact velocities ranging from 0.3 to 3.0 km/s, resulting in a total of 171 high-fidelity simulations. The bulk shock state from inert 2D mesoscale simulations aligns with the law of mixtures, while temperature and pressure fluctuations strongly correlate with multilayer geometry types. A new metric dubbed the “hot spot probability integral” shows a greater dependence on a tilt angle than interfacial roughness.
Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni xAl y, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni xAl y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.
The mechanical response of a component is affected by defects, such as porosity, arising from the laser powder bed fusion (LPBF) fabrication process. Thus, it is important to develop accurate and efficient inspection methods for identifying porosity. In this work, porosity identified in an X-ray computed tomography (XCT) volume of a Ti-5553 coupon was compared to pores identified in a serial sectioned volume that represented the ground truth. The porosity of the XCT scan was identified using contrast-based, ISO-based, and machine learning (ML) methods for segmentation. Large inherent porosity was easy to identify, but the ISO thresholding still struggled due to the intensity gradient resulting from both the beam hardening in XCT and the uneven lighting of the serial sectioning panels. Further, the results show that ML-based methods were better suited for identifying small pores and reducing the amount of false positives. Additionally, high strain-rate impact testing was done on some of the XCT samples as well as post-mortem XCT inspection, and the same suite of segmentation and quantification tools were used to identify the large spallation cavities. The comparison of porosity pre- and post-mortem provides insight on the influence of the LPBF porosity on the formation of spall cavities.
This report summarizes the result of a one year seedling project to investigate unusual twinning behavior in shock loaded additive Ti5552. The twinning behavior only occurs when the β phase of Ti5553 is metastable, and it appears to be a type of double twin involving two different twin variants, first a {332}⟨113⟩ twin forms before being consumed by a specific {112}⟨111⟩ twin variant to create a 20o. This behavior has only been detected during shock loading around incipient spall damage. The twinning is investigated by performing postmortem EBSD and PED analysis of gas-gun loaded specimens and preliminary molecular dynamics simulations.
We measured the austenitic FeCr18Ni12.5 stainless steel Hugoniot as a function of crystallographic direction to approximately 60 GPa. We shock-compressed FeCr18Ni12.5 samples oriented along ⟨ 100 ⟩ , ⟨ 110 ⟩ , and ⟨ 111 ⟩ to mean stresses ranging 30.5-58.1 GPa via Ta plate impact in a large-bore powder gun and measured the free-surface velocities with laser interferometry. We unambiguously observed the largest post-shock free-surface velocity along ⟨ 100 ⟩ in each experiment, which consequently produced the lowest shock velocity along that orientation. However, the propagation of experimental uncertainties through the impedance matching scheme used to compute the shock velocity produced sufficient uncertainty overlap to preclude definitive conclusion of Hugoniot anisotropy.
Both shock and shockless compression experiments were performed on laser powder bed fusion (LPBF) Ti-5Al-5V-5Mo-3Cr (Ti-5553) to peak compressive stresses near 15 GPa. Experiments were performed on the as-built material, containing a purely β (body centered cubic) microstructure, and two differing heat treatments resulting in a dual phase α (hexagonal close packed) and β microstructure. The Hugoniot, Hugoniot elastic limit (HEL), and spallation strength were measured and compared to wrought Ti-6Al-4V (Ti-64). The results indicate the LPBF Ti-5553 Hugoniot response is similar between heat treatments and to Ti-64. The HEL stress observed in the LPBF Ti-5553 was considerably higher than Ti-64, with the as-built, fully β alloy exhibiting the largest values. The spallation strength of the LPBF Ti-5553 was also similar to Ti-64. Clear evidence of initial porosity serving as initiation sites for spallation damage was observed when comparing computed tomography measurements before and after loading. Post-mortem scanning electron microscopy images of the recovered spallation samples showed no evidence of retained phase changes near the spall plane. The spall plane was found to have kinks aligned with the loading direction near areas with large concentrations of twin-like, crystallographic defects in the as-built condition. For the heat-treated samples, the concentrations of twin-like, crystallographic defects were absent, and no preference for failure at the interface between the α and β phases was observed.
Pressure-shear plate impact experiments were performed to quantify flow strength of wrought, as-built additively manufactured (AM), and heat-treated and recrystallized AM 304 L stainless steel (SS304L) under combined loading. Impact velocities spanned between 0.03 and 0.24 mm/μs, resulting in corresponding pressures of 0.62–5.93 GPa. Flow strength measurements are comparable for the sample variants across the studied loading conditions; however, shear wave structures significantly differ between sample type. Microstructurally aware simulations indicate local strain differences attributed to anisotropic elastic constants of large grains (~1 mm) in the as-built and heat-treated AM may impede the ability to uniformly transmit a shear wave.
Diffusion bonding of two immiscible, binary metallic systems, Cu-Ta and Cu-W was employed to make repeatable and predictable dual-layer impactors for shock-reshock experiments. The diffusion bonded impactors were characterized using ultrasonic imaging and optical microscopy to ensure bonding and the absence of excessive Cu grain coarsening. The diffusion bonded impactors were launched via a two-stage gas gun at [100] LiF windows instrumented with multiple interferometry probes spanning nearly the entire impactor area. Consistent interferometry data was obtained from all experiments with no evidence of release prior to recompression, indicating a uniform bond. Comparisons to hydrocode simulations show excellent agreement for all experiments, facilitating easy application of these impactors to future experiments.
Additive manufactured Ti-5Al-5V-5Mo-3Cr (Ti-5553) is being considered as an AM repair material for engineering applications because of its superior strength properties compared to other titanium alloys. Here, we describe the failure mechanisms observed through computed tomography, electron backscatter diffraction (EBSD), and scanning electron microscopy (SEM) of spall damage as a result of tensile failure in as-built and annealed Ti-5553. We also investigate the phase stability in native powder, as-built and annealed Ti-5553 through diamond anvil cell (DAC) and ramp compression experiments. We then explore the effect of tensile loading on a sample containing an interface between a Ti-6Al-V4 (Ti-64) baseplate and additively manufactured Ti-5553 layer. Post-mortem materials characterization showed spallation occurred in regions of initial porosity and the interface provides a nucleation site for spall damage below the spall strength of Ti-5553. Preliminary peridynamics modeling of the dynamic experiments is described. Finally, we discuss further development of Stochastic Parallel PARticle Kinteic Simulator (SPPARKS) Monte Carlo (MC) capabilities to include the integration of alpha (α)-phase and microstructural simulations for this multiphase titanium alloy.
