Update on multi-megabar shockless compression at Sandia’s Z machine (2024)
Abstract not provided.
Abstract not provided.
Journal of Applied Physics
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. 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 .
Journal of Applied Physics
Dynamic shockless compression experiments provide the ability to explore material behavior at extreme pressures but relatively low temperatures. Typically, the data from these types of experiments are interpreted through an analytic method called Lagrangian analysis. In this work, alternative analysis methods are explored using modern statistical methods. Specifically, Bayesian model calibration is applied to a new set of platinum data shocklessly compressed to 570 GPa. Several platinum equation-of-state models are evaluated, including traditional parametric forms as well as a novel non-parametric model concept. The results are compared to those in Paper I obtained by inverse Lagrangian analysis. The comparisons suggest that Bayesian calibration is not only a viable framework for precise quantification of the compression path, but also reveals insights pertaining to trade-offs surrounding model form selection, sensitivities of the relevant experimental uncertainties, and assumptions and limitations within Lagrangian analysis. The non-parametric model method, in particular, is found to give precise unbiased results and is expected to be useful over a wide range of applications. The calibration results in estimates of the platinum principal isentrope over the full range of experimental pressures to a standard error of 1.6%, which extends the results from Paper I while maintaining the high precision required for the platinum pressure standard.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Journal of Dynamic Behavior of Materials
Magnetic loading was used to shocklessly compress four different metals to extreme pressures. Velocimetry monitored the behavior of the material as it was loaded to a desired peak state and then decompressed back down to lower pressures. Two distinct analysis methods, including a wave profile analysis and a novel Bayesian calibration approach, were employed to estimate quantitative strength metrics associated with the loading reversal. Specifically, we report for the first time on strength estimates for tantalum, gold, platinum, and iridium under shockless compression at strain rates of ∼ 5 × 10 5/s in the pressure range of ∼ 100–400 GPa. The magnitude of the shear stresses supported by the different metals under these extreme conditions are surprisingly similar, representing a dramatic departure from ambient conditions.
Abstract not provided.
Journal of Geophysical Research. Solid Earth
The outer core of the Earth is composed primarily of liquid iron, and the inner core boundary is governed by the intersection of the melt line and the geotherm. While there are many studies on the thermodynamic equation of state for solid iron, the equation of state of liquid iron is relatively unexplored. In this work, we use dynamic compression to diagnose the high-pressure liquid equation of state of iron by utilizing the shock-ramp capability at Sandia National Laboratories’ Z-Machine. This technique enables measurements of material states off the Hugoniot by initially shocking samples and subsequently driving a further, shockless compression. Planetary studies benefit greatly from isentropic, off-Hugoniot experiments since they can cover pressure-temperature (P-T) conditions that are close to adiabatic profiles found in planetary interiors. We used this method to drive iron to P-T conditions similar to those of the Earth’s outer-inner core boundary, along an elevated-temperature isentrope in the liquid from 275 GPa to 400 GPa. We derive the equation of state using a hybrid backward integration – forward Lagrangian technique on particle velocity traces to determine the pressure-density history of the sample. Our results are in excellent agreement with SESAME 92141, a previously published equation of state table. With our data and previous experimental data on liquid iron we provide new information on the iron melting line and derive new parameters for a Vinet-based equation of state. The table and our parameterized equation of state are applied to provide an updated means of modeling the pressure, mass, and density of liquid iron cores in exoplanetary interiors.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Physical Review Letters
Nanosecond in situ x-ray diffraction and simultaneous velocimetry measurements were used to determine the crystal structure and pressure, respectively, of ramp-compressed aluminum at stress states between 111 and 475 GPa. The solid-solid Al phase transformations, fcc-hcp and hcp-bcc, are observed at 216±9 and 321±12 GPa, respectively, with the bcc phase persisting to 475 GPa. The high-pressure crystallographic texture of the hcp and bcc phases suggests close-packed or nearly close-packed lattice planes remain parallel through both transformations.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
The CHEDS researchers are engaged in a collaborative research project to study the properties of iron and iron alloys under Earth’s core conditions. The Earth’s core, inner and outer, is composed primarily of iron, thus studying iron and iron alloys at high pressure and temperature conditions will give the best estimate of its properties. Also, comparing studies of iron alloys with known properties of the core can constrain the potential light element compositions found within the core, such as fitting sound speeds and densities of iron alloys to established inner- Earth models. One of the lesser established properties of the core is the thermal conductivity, where current estimates vary by a factor of three. Therefore, one of the primary goals of this collaboration is to make relevant measurements to elucidate this conductivity.
AIP Conference Proceedings
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.
Journal of Applied Physics
Single crystal lithium fluoride (LiF), oriented [100], was shock loaded and subsequently shocklessly compressed in two experiments at the Z Machine. Velocimetry measurements were employed to obtain an impactor velocity, shock transit times, and in-situ particle velocities for LiF samples up to ∼1.8 mm thick. A dual thickness Lagrangian analysis was performed on the in-situ velocimetry data to obtain the mechanical response along the loading path of these experiments. An elastic response was observed on one experiment during initial shockless compression from 100 GPa before yielding. The relatively large thickness differences utilized for the dual sample analyses (up to ∼1.8 mm) combined with a relative timing accuracy of ∼0.2 ns resulted in an uncertainty of less than 1% on density and stress at ∼200 GPa peak loading on one experiment and <4% on peak loading at ∼330 GPa for another. The stress-density analyses from these experiments compare favorably with recent equation of state models for LiF.
Journal of Applied Physics
An understanding of the mechanical and optical properties of lithium fluoride (LiF) is essential to its use as a transparent tamper and window for dynamic materials experiments. In order to improve models for this material, we applied iterative Lagrangian analysis to ten independent sets of data from magnetically driven planar shockless compression experiments on single crystal [100] LiF to pressures as high as 350 GPa. We found that the compression response disagreed with a prevalent tabular equation of state for LiF that is commonly used to interpret shockless compression experiments. We also present complementary data from ab initio calculations performed using the diffusion quantum Monte Carlo method. The agreement between these two data sets lends confidence to our interpretation. In order to aid in future experimental analysis, we have modified the tabular equation of state to match the new data. We have also extended knowledge of the optical properties of LiF via shock-compression and shockless compression experiments, refining the transmissibility limit, measuring the refractive index to ∼300 GPa, and confirming the nonlinear dependence of the refractive index on density. We present a new model for the refractive index of LiF that includes temperature dependence and describe a procedure for correcting apparent velocity to true velocity for dynamic compression experiments.
Abstract not provided.
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.
Abstract not provided.
Physical Review Accelerators and Beams
In this study, we have developed a conceptual design of a next-generation pulsed-power accelerator that is optmized for driving megajoule-class dynamic-material-physics experiments at pressures as high as 1 TPa. The design is based on an accelerator architecture that is founded on three concepts: single-stage electrical-pulse compression, impedance matching, and transit-time-isolated drive circuits. Since much of the accelerator is water insulated, we refer to this machine as Neptune. The prime power source of Neptune consists of 600 independent impedance-matched Marx generators. As much as 0.8 MJ and 20 MA can be delivered in a 300-ns pulse to a 16-mΩ physics load; hence Neptune is a megajoule-class 20-MA arbitrary waveform generator. Neptune will allow the international scientific community to conduct dynamic equation-of-state, phase-transition, mechanical-property, and other material-physics experiments with a wide variety of well-defined drive-pressure time histories. Because Neptune can deliver on the order of a megajoule to a load, such experiments can be conducted on centimeter-scale samples at terapascal pressures with time histories as long as 1 μs.
Abstract not provided.
Abstract not provided.
Physical Review B
An absolute stress-density path for shocklessly compressed copper is obtained to over 450 GPa. A magnetic pressure drive is temporally tailored to generate shockless compression waves through over 2.5-mm-thick copper samples. The free-surface velocity data is analyzed for Lagrangian sound velocity using the iterative Lagrangian analysis (ILA) technique, which relies upon the method of characteristics. We correct for the effects of strength and plastic work heating to determine an isentropic compression path. By assuming a Debye model for the heat capacity, we can further correct the isentrope to an isotherm. Our determination of the isentrope and isotherm of copper represents a highly accurate pressure standard for copper to over 450 GPa.
Abstract not provided.
Physical Review Special Topics - Accelerators and Beams
We have developed the design of Thor: a pulsed power accelerator that delivers a precisely shaped current pulse with a peak value as high as 7 MA to a strip-line load. The peak magnetic pressure achieved within a 1-cm-wide load is as high as 100 GPa. Thor is powered by as many as 288 decoupled and transit-time isolated bricks. Each brick consists of a single switch and two capacitors connected electrically in series. The bricks can be individually triggered to achieve a high degree of current pulse tailoring. Because the accelerator is impedance matched throughout, capacitor energy is delivered to the strip-line load with an efficiency as high as 50%. We used an iterative finite element method (FEM), circuit, and magnetohydrodynamic simulations to develop an optimized accelerator design. When powered by 96 bricks, Thor delivers as much as 4.1 MA to a load, and achieves peak magnetic pressures as high as 65 GPa. When powered by 288 bricks, Thor delivers as much as 6.9 MA to a load, and achieves magnetic pressures as high as 170 GPa. We have developed an algebraic calculational procedure that uses the single brick basis function to determine the brick-triggering sequence necessary to generate a highly tailored current pulse time history for shockless loading of samples. Thor will drive a wide variety of magnetically driven shockless ramp compression, shockless flyer plate, shock-ramp, equation of state, material strength, phase transition, and other advanced material physics experiments.
Abstract not provided.
Abstract not provided.
Physical Review Letters
We report on experiments where cylindrical beryllium liners filled with liquid deuterium were compressed to extreme pressure and density with current pulse shaping. In one set of experiments the pressure at stagnation is inferred to be & 100 Mbar using penetrating radiography. A peak liner convergence ratio (initial radius over final radius) of 7.6 was measured resulting in an average deuterium density of 10 g=cm3 and areal density of 0.45 g=cm2. The stagnation shock propagating radially outward through the liner wall was directly measured with a strength of ≈ 120 Mbar. In a second set of experiments the liner was imploded to a peak convergence of 19 resulting in a density of 55 g=cm3 and areal density of 0.5 g=cm2. The pressure at stagnation in this experiment is estimated to be 2 Gbar. This platform enables the study of high-pressure, high-density, implosion deceleration and stagnation dynamics at spatial scales that are readily diagnosable (R ~ 0.1 -- 0.4 mm). Thus, these experiments are directly relevant to both Inertial Con nement Fusion and the study of material properties under extreme conditions.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Journal of Applied Physics
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 330GPa 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.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
AIP Conference Proceedings
We describe a technique for measuring the pressure and density of a metallic solid, shocklessly compressed to multi-megabar pressure, through x-ray radiography of a magnetically driven, cylindrical liner implosion. Shockless compression of the liner produces material states that correspond approximately to the principal compression isentrope (quasi-isentrope). This technique is used to determine the principal quasi-isentrope of solid beryllium to a peak pressure of 2.4 Mbar from x-ray images of a high current (20 MA), fast (∼100 ns) liner implosion. © 2012 American Institute of Physics.
Abstract not provided.
Physics of Plasmas
Current pulse shaping techniques, originally developed for planar dynamic material experiments on the Z-machine [M. K. Matzen, Phys. Plasmas 12, 055503 (2005)], are adapted to the design of controlled cylindrical liner implosions. By driving these targets with a current pulse shape that prevents shock formation inside the liner, shock heating is avoided along with the corresponding decrease in electrical conductivity ahead of the magnetic diffusion wave penetrating the liner. This results in an imploding liner with a significant amount of its mass in the solid phase and at multi-megabar pressures. Pressures in the solid region of a shaped pulse driven beryllium liner fielded on the Z-machine are inferred to 5.5 Mbar, while simulations suggest implosion velocities greater than 50 kms-1. These solid liner experiments are diagnosed with multi-frame monochromatic x-ray backlighting which is used to infer the material density and pressure. This work has led to a new platform on the Z-machine that can be used to perform off-Hugoniot measurements at higher pressures than are accessible through magnetically driven planar geometries. © 2012 American Institute of Physics.
Abstract not provided.
Abstract not provided.
Digest of Technical Papers-IEEE International Pulsed Power Conference
The success of dynamic materials properties research at Sandia National Laboratories has led to research into ultra-low impedance, compact pulsed power systems capable of multi-MA shaped current pulses with rise times ranging from 220-500 ns. The Genesis design consists of two hundred and forty 200 kV, 80 kA modules connected in parallel to a solid dielectric disk transmission line and is capable of producing 280 kbar of magnetic pressure (>500 kbar pressure in high Z materials) in a 1.75 nH, 20 mm wide stripline load. Stripline loads operating under these conditions expand during the experiment resulting in a time-varying load that can impact the performance and lifetime of the system. This paper provides analysis of time-varying stripline loads and the impact of these loads on system performance. Further, an approach to reduce dielectric stress levels through active damping is presented as a means to increase system reliability and lifetime. © 2011 IEEE.
Digest of Technical Papers-IEEE International Pulsed Power Conference
Genesis is a compact pulsed power platform designed by Sandia National Laboratories to generate precision shaped multi-MA current waves with a rise time of 200-500 ns. In this system, two hundred and forty, 200 kV, 80 kA modules are selectively triggered to produce 280 kbar of magnetic pressure (>500 kbar pressure in high Z materials) in a stripline load for dynamic materials properties research. This new capability incorporates the use of solid dielectrics to reduce system inductance and size, programmable current shaping, and gas switches that must perform over a large range of operating conditions. Research has continued on this technology base with a focus on demonstrating the integrated performance of key concepts into a Genesis-like prototype called Protogen. Protogen measures approximately 1.4 m by 1.4 m and is designed to hold twelve Genesis modules. A fixed inductance load will allow rep-rate operation for component reliability and system lifetime experiments at the extreme electric field operating conditions expected in Genesis. © 2011 IEEE.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.