Publications

Results 1–50 of 65
Skip to search filters

A new equation of state for copper

Carpenter, John H.

A new copper equation of state is developed utilizing the available experimental data in addition to recent theoretical calculations. Semi-empirical models are fit to the data and the results are tabulated in the SNL SESAME format. Comparison to other copper EOS tables are given, along with recommendations of which tables provide the best accuracy.

More Details

Improved reference system for the corrected rigid spheres equation of state model

Journal of Applied Physics

Cowen, Benjamin J.; Carpenter, John H.

The Corrected Rigid Spheres (CRIS) equation of state (EOS) model [Kerley, J. Chem. Phys. 73, 469 (1980); 73, 478 (1980); 73, 487 (1980)], developed from fluid perturbation theory using a hard sphere reference system, has been successfully used to calculate the EOS of many materials, including gases and metals. The radial distribution function (RDF) plays a pivotal role in choosing the sphere diameter, through a variational principle, as well as the thermodynamic response. Despite its success, the CRIS model has some shortcomings in that it predicts too large a temperature for liquid-vapor critical points, can break down at large compression, and is computationally expensive. We first demonstrate that an improved analytic representation of the hard sphere RDF does not alleviate these issues. Relaxing the strict adherence of the RDF to hard spheres allows an accurate fit to the isotherms and vapor dome of the Lennard-Jones fluid using an arbitrary reference system. The second order correction is eliminated, limiting the breakdown at large compression and significantly reducing the computation cost. The transferability of the new model to real systems is demonstrated on argon, with an improved vapor dome compared to the original CRIS model.

More Details

Shock compression of strongly correlated oxides: A liquid-regime equation of state for cerium(IV) oxide

Physical Review B

Weck, Philippe F.; Cochrane, Kyle C.; Root, Seth R.; Lane, J.M.; Shulenburger, Luke N.; Carpenter, John H.; Sjostrom, Travis; Mattsson, Thomas M.; Vogler, Tracy V.

The shock Hugoniot for full-density and porous CeO2 was investigated in the liquid regime using ab initio molecular dynamics (AIMD) simulations with Erpenbeck's approach based on the Rankine-Hugoniot jump conditions. The phase space was sampled by carrying out NVT simulations for isotherms between 6000 and 100 000 K and densities ranging from ρ=2.5 to 20g/cm3. The impact of on-site Coulomb interaction corrections +U on the equation of state (EOS) obtained from AIMD simulations was assessed by direct comparison with results from standard density functional theory simulations. Classical molecular dynamics (CMD) simulations were also performed to model atomic-scale shock compression of larger porous CeO2 models. Results from AIMD and CMD compression simulations compare favorably with Z-machine shock data to 525 GPa and gas-gun data to 109 GPa for porous CeO2 samples. Using results from AIMD simulations, an accurate liquid-regime Mie-Grüneisen EOS was built for CeO2. In addition, a revised multiphase SESAME-Type EOS was constrained using AIMD results and experimental data generated in this work. This study demonstrates the necessity of acquiring data in the porous regime to increase the reliability of existing analytical EOS models.

More Details

Magnetically launched flyer plate technique for probing electrical conductivity of compressed copper

Journal of Applied Physics

Cochrane, Kyle C.; Lemke, Raymond W.; Riford, Z.; Carpenter, John H.

The electrical conductivity of materials under extremes of temperature and pressure is of crucial importance for a wide variety of phenomena, including planetary modeling, inertial confinement fusion, and pulsed power based dynamic materials experiments. There is a dearth of experimental techniques and data for highly compressed materials, even at known states such as along the principal isentrope and Hugoniot, where many pulsed power experiments occur. We present a method for developing, calibrating, and validating material conductivity models as used in magnetohydrodynamic (MHD) simulations. The difficulty in calibrating a conductivity model is in knowing where the model should be modified. Our method isolates those regions that will have an impact. It also quantitatively prioritizes which regions will have the most beneficial impact. Finally, it tracks the quantitative improvements to the conductivity model during each incremental adjustment. In this paper, we use an experiment on Sandia National Laboratories Z-machine to isentropically launch multiple flyer plates and, with the MHD code ALEGRA and the optimization code DAKOTA, calibrated the conductivity such that we matched an experimental figure of merit to +/-1%.

More Details

Probing off-Hugoniot states in Ta, Cu, and Al to 1000 GPa compression with magnetically driven liner implosions

Journal of Applied Physics

Lemke, Raymond W.; Dolan, Daniel H.; Dalton, D.G.; Brown, Justin L.; Tomlinson, K.; Robertson, G.R.; Knudson, Marcus D.; Harding, Eric H.; Mattsson, A.E.; Carpenter, John H.; Drake, Richard R.; Cochrane, Kyle C.; Blue, B.E.; Robinson, Allen C.; Mattsson, Thomas M.

We report on a new technique for obtaining off-Hugoniot pressure vs. density data for solid metals compressed to extreme pressure by a magnetically driven liner implosion on the Z-machine (Z) at Sandia National Laboratories. In our experiments, the liner comprises inner and outer metal tubes. The inner tube is composed of a sample material (e.g., Ta and Cu) whose compressed state is to be inferred. The outer tube is composed of Al and serves as the current carrying cathode. Another aluminum liner at much larger radius serves as the anode. A shaped current pulse quasi-isentropically compresses the sample as it implodes. The iterative method used to infer pressure vs. density requires two velocity measurements. Photonic Doppler velocimetry probes measure the implosion velocity of the free (inner) surface of the sample material and the explosion velocity of the anode free (outer) surface. These two velocities are used in conjunction with magnetohydrodynamic simulation and mathematical optimization to obtain the current driving the liner implosion, and to infer pressure and density in the sample through maximum compression. This new equation of state calibration technique is illustrated using a simulated experiment with a Cu sample. Monte Carlo uncertainty quantification of synthetic data establishes convergence criteria for experiments. Results are presented from experiments with Al/Ta, Al/Cu, and Al liners. Symmetric liner implosion with quasi-isentropic compression to peak pressure ∼1000 GPa is achieved in all cases. These experiments exhibit unexpectedly softer behavior above 200 GPa, which we conjecture is related to differences in the actual and modeled properties of aluminum.

More Details

ALEGRA Update: Modernization and Resilience Progress

Robinson, Allen C.; Petney, Sharon P.; Drake, Richard R.; Weirs, Vincent G.; Adams, Brian M.; Vigil, Dena V.; Carpenter, John H.; Garasi, Christopher J.; Wong, Michael K.; Robbins, Joshua R.; Siefert, Christopher S.; Strack, Otto E.; Wills, Ann E.; Trucano, Timothy G.; Bochev, Pavel B.; Summers, Randall M.; Stewart, James R.; Ober, Curtis C.; Rider, William J.; Haill, Thomas A.; Lemke, Raymond W.; Cochrane, Kyle C.; Desjarlais, Michael P.; Love, Edward L.; Voth, Thomas E.; Mosso, Stewart J.; Niederhaus, John H.

Abstract not provided.

Fundamental issues in the representation and propagation of uncertain equation of state information in shock hydrodynamics

Computers and Fluids

Robinson, Allen C.; Berry, Robert D.; Carpenter, John H.; Debusschere, Bert D.; Drake, Richard R.; Mattsson, A.E.; Rider, William J.

Uncertainty quantification (UQ) deals with providing reasonable estimates of the uncertainties associated with an engineering model and propagating them to final engineering quantities of interest. We present a conceptual UQ framework for the case of shock hydrodynamics with Euler's equations where the uncertainties are assumed to lie principally in the equation of state (EOS). In this paper we consider experimental data as providing both data and an estimate of data uncertainty. We propose a specific Bayesian inference approach for characterizing EOS uncertainty in thermodynamic phase space. We show how this approach provides a natural and efficient methodology for transferring data uncertainty to engineering outputs through an EOS representation that understands and deals consistently with parameter correlations as sensed in the data.Historically, complex multiphase EOSs have been built utilizing tables as the delivery mechanism in order to amortize the cost of creation of the tables over many subsequent continuum scale runs. Once UQ enters into the picture, however, the proper operational paradigm for multiphase tables become much less clear. Using a simple single-phase Mie-Grüneisen model we experiment with several approaches and demonstrate how uncertainty can be represented. We also show how the quality of the tabular representation is of key importance. As a first step, we demonstrate a particular tabular approach for the Mie-Grüneisen model which when extended to multiphase tables should have value for designing a UQ-enabled shock hydrodynamic modeling approach that is not only theoretically sound but also robust, useful, and acceptable to the modeling community. We also propose an approach to separate data uncertainty from modeling error in the EOS. © 2012 Elsevier Ltd.

More Details

General purpose steam table library :

Carpenter, John H.; Belcourt, Kenneth N.

Completion of the CASL L3 milestone THM.CFD.P7.04 provides a general purpose tabular interpolation library for material properties to support, in particular, standardized models for steam properties. The software consists of three parts, implementations of analytic steam models, a code to generate tables from those models, and an interpolation package to interface the tables to CFD codes such as Hydra-TH. Verification of the standard model is maintained through the entire train of routines. The performance of interpolation package exceeds that of freely available analytic implementation of the steam properties by over an order of magnitude.

More Details

Tabular water properties interface for Hydra-TH :

Carpenter, John H.; Belcourt, Kenneth N.

Completion of the CASL L3 milestone THM.CFD.P6.03 provides a tabular material properties capability to the Hydra code. A tabular interpolation package used in Sandia codes was modified to support the needs of multi-phase solvers in Hydra. Use of the interface is described. The package was released to Hydra under a government use license. A dummy physics was created in Hydra to prototype use of the interpolation routines. Finally, a test using the dummy physics verifies the correct behavior of the interpolation for a test water table. 3

More Details
Results 1–50 of 65
Results 1–50 of 65