Search NASA⌕ Search

Engineering topics

Brown, Justin L.

Publications and source records attributed to Brown, Justin L..

At least 19 records

Quantifying uncertainty in analysis of shockless dynamic compression experiments on platinum. I. Inverse Lagrangian analysis

Absolute measurements of solid-material compressibility by magnetically driven shockless dynamic compression experiments to multi-megabar pressures have the potential to greatly improve the accuracy and precision of pressure calibration standards for use in diamond anvil cell experiments. Here, to this end, we apply characteristics-based inverse Lagrangian analysis (ILA) to 11 sets of ramp-compression data on pure platinum (Pt) metal and then reduce the resulting weighted-mean stress–strain curve to the principal isentrope and room-temperature isotherm using simple models for yield stress and Grüneisen parameter. We introduce several improvements to methods for ILA and quasi-isentrope reduction, the latter including calculation of corrections in wave speed instead of stress and pressure to render results largely independent of initial yield stress while enforcing thermodynamic consistency near zero pressure. More importantly, we quantify in detail the propagation of experimental uncertainty through ILA and model uncertainty through quasi-isentrope reduction, considering all potential sources of error except the electrode and window material models used in ILA. Compared to previous approaches, we find larger uncertainty in longitudinal stress. Monte Carlo analysis demonstrates that uncertainty in the yield-stress model constitutes by far the largest contribution to uncertainty in quasi-isentrope reduction corrections. We present a new room-temperature isotherm for Pt up to 444 GPa, with 1-sigma uncertainty at that pressure of just under ±1.2%; the latter is about a factor of three smaller than uncertainty previously reported for multi-megabar ramp-compression experiments on Pt. The result is well represented by a Vinet-form compression curve with (isothermal) bulk modulus K 0 = 270.3 ± 3.8 GPa, pressure derivative K$^{'}_{0}$= 5.66 ± 0.10, and correlation coefficient $R_{K_{0},K^{'}_{0}}$= –0.843.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Inference of strength and phase transition kinetics in dynamically-compressed tin

Dynamic compression experiments in condensed matter are of interest in part because they provide opportunities to examine material response under extreme conditions; however, the inference of material behavior from dynamic experiments is challenging in the presence of phase transitions exhibiting kinetic processes. Demonstrating an approach to quantitative interpretation of such dynamic experiments, here we present a Bayesian model calibration of strength and phase transformation parameters to data drawn from pulsed power and gas gun shot experiments. The posterior predictions of the Bayesian model capture the experimental measurements and account for the various uncertainties in the experimental configurations. This holistic approach to model calibration utilizing multiple types of experimental data identifies important cross correlations among kinetics, strength, and the phase boundary. Improved insight into potential sources of current model form error is provided by comparing the differences between calibrations against different subsets of the experimental data.

42 ENGINEERING↗

Quasi-isentropic compression of an additively manufactured aluminum alloy to 14.8 GPa

We uniaxially and shocklessly compressed an additively manufactured aluminum alloy, laser powder bed fusion (LPBF) AlSi10Mg, to peak stresses ranging from 4.4 to 14.8 GPa at peak strain rates on the order of [Formula: see text] via a series of magnetic loading experiments to measure the principal isentrope, yield strength, and shear modulus as a function of material orientation and applied stress. We did not observe significant anisotropy in any of the measured quantities. We found that the principal isentrope, within the uncertainty and up to our peak stress, overlaps the material’s Hugoniot. We measured yield strengths and shear moduli ranging from 0.28 to 0.81 GPa and 36 to 52 GPa, respectively. Our results indicate that LPBF AlSi10Mg behaves similarly to wrought Al alloys under quasi-isentropic compression.

Brown, Nathan P.↗

A Multi-Phase Modeling Framework Suitable for Dynamic Applications

Under dynamic loading conditions and the associated extreme conditions many metals will undergo phase transformations. The change in crystal structure associated with solid–solid phase transformations can significantly alter the subsequent mechanical response of the material. For the interpretation of experiments involving dynamic loading it is beneficial to have a modeling framework that captures key features of the material response while remaining relatively simple. We introduce a candidate framework and apply it to the metal tin to highlight a range of behaviors that are captured by the model. We also discuss potential extensions to capture additional behaviors that could be important for certain materials and loading scenarios. The model is useful for analysis of results from dynamic experiments and offers a point of departure for more complex model formulations.

36 MATERIALS SCIENCE↗

Model for the Solid–Liquid Interfacial Free Energy at High Pressures

The free energy involved in the formation of an interface between two phases (e.g., a solid–liquid interface) is referred to as the interfacial free energy. For the case of solidification, the interfacial free energy dictates the height of the energy barrier required to nucleate stable clusters of the newly forming solid phase and is essential for producing an accurate solidification kinetics model using classical nucleation theory (CNT)-based methods. While various methods have been proposed for modeling the interfacial free energy for solid–liquid interfaces in prior literature, many of these formulations involve making restrictive assumptions or approximations, such as the system being at or near equilibrium (i.e., the system temperature is approximately equal to the melt temperature) or that the system is at pressures close to atmospheric. However, these approximations and assumptions may break down in highly non-equilibrium situations, such as in dynamic-compression experiments where metastable liquids that are undercooled by hundreds of kelvin or over-pressurized by several gigapascals or more are formed before eventually solidifying. We derive a solid–liquid interfacial free-energy model for such high-pressure conditions by considering the enthalpies of interactions between pairs of atoms or molecules. We also consider the contribution of interface roughness (disordering) by incorporating a multilayer interface model known as the Temkin n-layer model. Our formulation is applicable to a diverse variety of materials, and we demonstrate it by developing models specifically for two different materials: water and gallium. We apply our interfacial free-energy formulation to CNT-based kinetics simulations of several suites of dynamic-compression experiments that cause liquid water to solidify to the high-pressure solid polymorph ice VII and have found good agreement to the observed kinetics with only minor empirical fitting.

36 MATERIALS SCIENCE↗

A compact x-ray diffraction system for dynamic compression experiments on pulsed-power generators

Pulsed-power generators can produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies using the magnetic loading technique. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength–strain rate dependence, structural phase transitions, and density of crystal defects, such as dislocations. Here, we present a cost-effective, compact, pulsed x-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically driven ramp compression of materials with a single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials’ community to investigate in situ dynamic phase transitions critical to equation of states. Finally, we present results using this new diagnostic to evaluate lattice compression in Zr and Al and to capture signatures of phase transitions in CdS.

47 OTHER INSTRUMENTATION↗

A Platform-Independent X-ray Diffraction Diagnostic for Phase Transition Kinetics in Traditional and Synthetic Microstructure Materials (LDRD Project 213088 Final Report)

Pulsed-power generators using the magnetic loading technique are able to produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength-strain rate dependence, structural phase transitions, and density of crystal defects such as dislocations. Here, we present a cost effective, compact X-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically-driven ramp compression of materials with single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials community. The success in fielding this new XRD diagnostic dramatically improves our predictive capability and understanding of rate-dependent behavior at or near phase transition. As Sandia plans the next-generation pulse-power driver platform, a key element needed to deliver new state-of-the-art experiments will be having the necessary diagnostic tools to probe new regimes and phenomena. These diagnostics need to be as versatile, compact, and portable as they are powerful. The development of a platform-independent XRD diagnostic gives Sandia researchers a new window to study the microstructure and phase dynamics of materials under load. This project has paved the way for phase transition research in a variety of materials with mission interest.

36 MATERIALS SCIENCE↗

Verification of experimental dynamic strength methods with atomistic ramp-release simulations

Material strength and moduli can be determined from dynamic high-pressure ramp-release experiments using an indirect method of Lagrangian wave profile analysis of surface velocities. This method, termed self-consistent Lagrangian analysis (SCLA), has been difficult to calibrate and corroborate with other experimental methods. Using nonequilibrium molecular dynamics, we validate the SCLA technique by demonstrating that it accurately predicts the same bulk modulus, shear modulus, and strength as those calculated from the full stress tensor data, especially where strain rate induced relaxation effects and wave attenuation are small. We show here that introducing a hold in the loading profile at peak pressure gives improved accuracy in the shear moduli and relaxation-adjusted strength by reducing the effect of wave attenuation. When rate-dependent effects coupled with wave attenuation are large, we find that Lagrangian analysis overpredicts the maximum unload wavespeed, leading to increased error in the measured dynamic shear modulus. Furthermore, these simulations provide insight into the definition of dynamic strength, as well as a plausible explanation for experimental disagreement in reported dynamic strength values.

36 MATERIALS SCIENCE↗

Results from new multi-megabar shockless compression experiments at the Z machine

Sandia’s Z Machine has been used to magnetically drive shockless compression of materials in a planar configuration to multi-megabar pressure levels, allowing accurate measurements of quasi-isentropic mechanical response at relatively low temperatures in the solid phase. This work details recent improvements to design and analysis of such experiments, including the use of new data on the mechanical and optical response of lithium fluoride windows. Comparison of windowed and free-surface data on copper to 350 GPa lends confidence to the window correction method. Preliminary results are presented on gold to 500 GPa and platinum to 450 GPa; both appear stiffer than existing models.

36 MATERIALS SCIENCE↗

Optimization of Isentropic Compression Loads on Current-Adder Pulsed Power Accelerator Architectures

The Thor pulsed power generator is being developed at Sandia National Laboratories. The design consists of up to 288 decoupled and transit time isolated capacitor-switch units, called “bricks”, that can be individually triggered to achieve a high degree of pulse tailoring for magnetically-driven isentropic compression experiments (ICE). The connecting transmission lines are impedance matched to the bricks, allowing the capacitor energy to be efficiently delivered to an ICE strip-line load with peak pressures of over 100 GPa. Thor will drive experiments to explore equation of state, material strength, and phase transition properties of a wide variety of materials. We present an optimization process for producing tailored current pulses, a requirement for many material studies, on the Thor generator. This technique, which is unique to the novel “current-adder” architecture used by Thor, entirely avoids the iterative use of complex circuit models to converge to the desired electrical pulse. We describe the optimization procedure for the Thor design and show results for various materials of interest. Also, we discuss the extension of these concepts to the megajoule-class Neptune machine design. Given this design, we are able to design shockless ramp-driven experiments in the 1 TPa range of material pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Strain-rate dependence of ramp-wave evolution and strength in tantalum

We have conducted molecular dynamics (MD) simulations of quasi-isentropic ramp-wave compression to very high pressures over a range of strain rates from 10 11 down to 10 8 1/s. Using scaling methods, we collapse wave profiles from various strain rates to a master profile curve, which shows deviations when material response is strain-rate dependent. Thus, we can show with precision where, and how, strain-rate dependence affects the ramp wave. We find that strain rate affects the stress-strain material response most dramatically at strains below 20%, and that above 30% strain the material response is largely independent of strain rate. We show good overall agreement with experimental stress-strain curves up to approximately 30% strain, above which simulated response is somewhat too stiff. We postulate that this could be due to our interatomic potential or to differences in grain structure and/or size between simulation and experiment. Strength is directly measured from per-atom stress tensor and shows significantly enhanced elastic response at the highest strain rates. As a result, this enhanced elastic response is less pronounced at higher pressures and at lower strain rates.

36 MATERIALS SCIENCE↗

Physically-based strength model of tantalum incorporating effects of temperature, strain rate and pressure

In this work, we develop a tantalum strength model that incorporates e ects of temperature, strain rate and pressure. Dislocation kink-pair theory is used to incorporate temperature and strain rate e ects while the pressure dependent yield is obtained through the pressure dependent shear modulus. Material constants used in the model are parameterized from tantalum single crystal tests and polycrystalline ramp compression experiments. It is shown that the proposed strength model agrees well with the temperature and strain rate dependent yield obtained from polycrystalline tantalum experiments. Furthermore, the model accurately reproduces the pressure dependent yield stresses up to 250 GPa. The proposed strength model is then used to conduct simulations of a Taylor cylinder impact test and validated with experiments. This approach provides a physically-based multi-scale strength model that is able to predict the plastic deformation of polycrystalline tantalum through a wide range of temperature, strain and pressure regimes.

36 MATERIALS SCIENCE↗

Analysis of shockless dynamic compression data on solids to multi-megabar pressures: Application to tantalum

Magnetically-driven, planar shockless-compression experiments to multi-megabar pressures were performed on tantalum samples using a stripline target geometry. Free-surface velocity waveforms were measured in 15 cases; nine of these in a dual-sample configuration with two samples of different thicknesses on opposing electrodes, and six in a single-sample configuration with a bare electrode opposite the sample. Details are given on the application of inverse Lagrangian analysis (ILA) to these data, including potential sources of error. The most significant source of systematic error, particularly for single-sample experiments, was found to arise from the pulse-shape dependent free-surface reflected wave interactions with the deviatoric-stress response of tantalum. This could cause local, possibly temporary, unloading of material from a ramp compressed state, and thus multi-value response in wave speed that invalidates the free-surface to in-material velocity mapping step of ILA. By averaging all 15 data sets, a final result for the principal quasi-isentrope of tantalum in stress-strain was obtained to a peak longitudinal stress of 330 GPa with conservative uncertainty bounds of ±4.5% in stress. The result agrees well with a tabular equation of state developed at Los Alamos National Laboratory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

2169 Steel Waveform Experiments

In support of LLNL efforts to develop multiscale models of a variety of materials, we have performed a set of eight gas gun impact experiments on 2169 steel (21% Cr, 6% Ni, 9% Mn, balance predominantly Fe). These experiments provided carefully controlled shock, reshock and release velocimetry data, with initial shock stresses ranging from 10 to 50 GPa (particle velocities from 0.25 to 1.05 km/s). Both windowed and free-surface measurements were included in this experiment set to increase the utility of the data set, as were samples ranging in thickness from 1 to 5 mm. Target physical phenomena included the elastic/plastic transition (Hugoniot elastic limit), the Hugoniot, any phase transition phenomena, and the release path (windowed and free-surface). The Hugoniot was found to be nearly linear, with no indications of the Fe α – ε phase transition. Releases were non-hysteretic, and relatively consistent between 3- and 5-mm–thick samples (the 3 mm samples giving slightly lower wavespeeds on release). Reshock tests with explosively welded impactors produced clean results; those with glue bonds showed transient releases prior to the arrival of the reshock, reducing their usefulness for deriving strength information. The free-surface samples, which were steps on a single piece of steel, showed lower wavespeeds for thin (1 mm) samples than for thicker (2 or 4 mm) samples. A configuration used for the last three shots allows release information to be determined from these free surface samples. The sample strength appears to increase with stress from ~1 GPa to ~ 3 GPa over this range, consistent with other recent work but about 40% above the Steinberg model.

Furnish, Michael David↗