Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ramp compression”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Non-iterative characteristics analysis for high-pressure ramp loading

Here in the canonical ramp compression experiment, a smoothly increasing load is applied to the surface of the sample, and the particle velocity history is measured at two or more different distances into the sample, at interfaces where the surface of the sample can be probed. The velocity histories are used to deduce a stress-density relation, usually using iterative Lagrangian analysis to account for the perturbing effect of the impedance mismatch at the interface. In that technique, a stress-density relation is assumed in order to correct for the perturbation and is adjusted until it becomes consistent with the deduced stress-density relation. This process is subject to the usual difficulties of nonlinear optimization, such as the existence of local minima (sensitivity to the initial guess), possible failure to converge, and relatively large computational effort. We show that, by considering the interaction of successive characteristics reaching a free surface, the stress-density relation can be deduced directly by recursion rather than iteration. This calculation is orders of magnitude faster than iterative analysis and does not require an initial guess. Direct recursion may be less suitable for very noisy data, but it was robust when applied to trial data. The stress-density relation deduced was identical to the result from iterative Lagrangian analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Two-phase equation of state for lithium fluoride

Here, we present an equation of state for the solid and liquid phases of lithium fluoride that covers a wide range of conditions from ambient pressure and temperature to the high pressures and temperatures exhibited in shock- and ramp-compression studies. The particular solid phase we have focused on in this work is the B1 phase. We have followed an approach where the pressure and heat-capacity functions of both phases are fit to experimental data and our own quantum molecular dynamics simulations and are then integrated in a thermodynamically consistent way to obtain the corresponding free-energy functions. This approach yields a two-phase equation of state that provides better overall agreement with experimental data than other equations of state for lithium fluoride, such as SESAME 7271v3, LEOS 2240, and the model presented by Smirnov. The last of these is a three-phase equation of state that predicts a B1–B2 transition along the shock Hugoniot at a pressure of about 140 GPa. This solid–solid transition has been a topic of speculation and debate in the literature for over 50 years, culminating in the work of Smirnov, who has developed the only potentially viable equation of state that allows for this transition. We explain why the proposed B1–B2 transition at 140 GPa is not consistent with recent velocimetry data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Calculating shock Hugoniot and isentropes using multiphase equation of state tables and application to shock and release of diamond ablators in inertial confinement fusion implosions

Advances in shock and ramp compression techniques now allow experimental access to unprecedented extreme conditions of pressure and temperature, providing a means to test theoretical models. Here, we describe a simple methodology to compute multi-phase shock Hugoniot and isentropes using multiphase equation of state tables. We treat explicitly the phase coexistence along the phase boundary to reveal the evolution of the sample as it undergoes the phase transformation in adiabatic conditions. We illustrate the method by calculating the predicted shock and shock-and-release behavior of diamond at conditions relevant for the initial stage of inertial confinement fusion implosions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Ramp-release experiments for strength measurements: Strain-rate dependence

This paper presents an enhanced analysis method for investigating material properties at high strain rates, extending the capability of established experimental techniques to gain more information. The ramp-release method has been applied to many experiments reported at high (≈10 5 − 10 6 s −1 ) strain-rates. More recent data gathered at the National Ignition Facility (NIF) has enabled higher (≈10 8 s −1 ) strain-rates to be studied. Here, we present an initial application of ramp-release analysis to NIF ramp-compression data, illustrating both the opportunities and the practical challenges of extending these methods to laser-driven platforms. The higher strain-rates accessed at the NIF mean that there is more strain-rate enhancement to strength, and the experimental configuration means that this enhancement is more readily seen in the data. This is enabled by the capability of avoiding peak-compression attenuation through the sample thickness with a designed hold period made possible by the pulse-shaping capability of NIF. We propose that this combination of experimental conditions and an enhanced analysis method enables the strain-rate enhancement to strength to be studied, and potentially for this to inform physics models at smaller scales than the continuum.

36 MATERIALS SCIENCE↗
DOE OSTI2016

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↗
DOE OSTI2021

Thermodynamic properties of ε -Fe with thermal electronic excitation effects on vibrational spectra

The thermodynamic properties of hexagonal-close-packed iron (ε–Fe) are essential for investigating the internal structure and dynamic properties of planetary cores. Despite their importance to planetary sciences, experimental investigations of ε–Fe at relevant conditions are still challenging. Therefore, ab initio calculations are crucial to elucidating the thermodynamic properties of this system. Here, we use a free energy calculation scheme based on the phonon gas model compatible with temperature-dependent phonon frequencies. We investigate the effects of electronic thermal excitations, which introduces a temperature dependence on phonon frequencies, and the implication for the thermodynamic properties of ε–Fe at extreme pressure (P) and temperature (T) conditions. We disregard phonon-phonon interactions, i.e., anharmonicity and their effect on phonon frequencies. Nevertheless, the current scheme is also applicable to T -dependent anharmonic frequencies. We conclude that the impact of thermal electronic excitations on vibrational properties is not significant up to ~4000 K at 200 GPa but should not be ignored at higher temperatures or pressures. However, the static free energy F st must always include the effect of thermal excitation fully in a continuum of T. Furthermore, our results for isentropic equations of state show good agreement with data from recent ramp compression experiments up to 1400 GPa conducted at the National Ignition Facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Equation of state and strength of diamond in high-pressure ramp loading

Diamond is used extensively as a component in high-energy-density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first-principles theoretical models of equilibrium multiphase EOS, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10 TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Grüneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. Furthermore, the resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Topological electride phase of sodium at high pressures and temperatures

Sodium, a textbook example of a nearly-free-electron metal, exhibits unforeseen pressure-induced behavior including the stabilization of numerous polymorphs—some possessing extremely complex unit cells—as well as a metal-to-insulator transition to the iconic hP4 phase. However, until recently, most of the experimental and theoretical studies on solid sodium have been restricted to the low-temperature regime. In this work, using ab initio evolutionary structure searches coupled with quasiharmonic calculations, we discover seven new phases of sodium that are more stable than the known hP4 phase at pressure-temperature conditions that were recently attained in ramp-compression experiments. From these, our calculations suggest that P6 3 /m Na is the ground state between ≈250 GPa at 710 K and ≈350 GPa at 1900 K. Electronic structure calculations show that this phase is a topological semimetal with a Dirac nodal surface that is protected by a nonsymmorphic symmetry S 2z and an electride owing to its non-nuclear charge localized within one-dimensional honeycomb channels and zero-dimensional cages. Our results highlight the complexity of dense sodium's free-energy landscape and intricate electronic structure that emerges at finite temperatures and conditions where ionic cores overlap.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Measurement of Body-Centered-Cubic Aluminum at 475 GPa [Observation of Body-Centered-Cubic Aluminum at 475 GPa]

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 GPa and 321 ± 12 GPa, respectively, with the bcc phase persisting to 475 GPa. Here, this is the first in situ observation of the high-pressure bcc phase of Al. High-pressure texture of the hcp and bcc phases suggests close-packed or nearly close-packed lattice planes remain parallel through both transformations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Extreme Hardening of Pb at High Pressure and Strain Rate

We study the high-pressure strength of Pb and Pb-4wt%Sb at the National Ignition Facility. We measure Rayleigh-Taylor growth of preformed ripples ramp compressed to ~400 GPa peak pressure, among the highest-pressure strength measurements ever reported on any platform. We find agreement with 2D simulations using the Improved Steinberg-Guinan strength model for body-centered-cubic Pb; the Pb-4wt%Sb alloy behaves similarly within the error bars. Here, the combination of high-rate, pressure-induced hardening and polymorphism yield an average inferred flow stress of ~3.8 GPa at high pressure, a ~250-fold increase, changing Pb from soft to extremely strong.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Limited-Ranged Two-Phase Iron Equation of State Model

This report describes the construction of a two-phase EOS for Fe that includes both (hcp) and liquid phases, and which is meant to be valid over the restricted range: V ϵ [4 : 7] °A 3 /atom, and T ϵ [1000 : 20000] K. It was recently used as the basis for a theoretical study of liquid → solid phase transformation kinetics in Fe which appeared in the Supplemental Material for a paper describing a series of NIF shock + ramp compression experiments aimed at constraining the high-P melt curve of Fe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Understanding Phase and Interfacial Effects of Spall Fracture in Additively Manufactured Ti-5Al-5V-5Mo-3Cr

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.

36 MATERIALS SCIENCE↗
DOE OSTI2023

REYNO: A Reactive Hydrodynamics Modeling Suite

REYNO is a Python library that is primarily designed as a platform for the agile development and testing of novel reactive flow models. Equations of state and multi-step rate laws can be implemented using the provided class hierarchy with different thermodynamic closure conditions. The library includes a one-dimensional Lagrangian hydrodynamic solver to simulate multi-material reactive or inert problems in planar, cylindrical, or spherical coordinates. Custom time-dependent boundary conditions are also available for problems such as ramp compression. This document describes the different numerical methods and algorithms that are included in REYNO.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Design and Characterization of a Lens-Coupled System for Dynamic X-Ray Diffraction

X-ray diffraction (XRD) is a necessary technique for understanding states of materials under static and dynamic loading conditions. The higher-pressure Equation of State (EOS) of many materials can only be explored via shock or ramp compression at temperatures and pressures of interest. While static XRD work has yielded EOS measurements in the 100 - 200 GPa regime, dynamic X-ray diffraction (DXRD) can explore EOS phases in the TPa regime, which closely resembles inner-core planetary conditions. DXRD hinges on the ability to measure the exact phase or phase change of a material while under dynamic loading conditions. Macroscopic diagnostic systems (e.g. velocimetry and pyrometry) can infer a phase change but not identify the specific phase entered by a material. While microscopic (atomic-level) diagnostic systems (e.g. DXRD) have been designed and implemented in Department of Energy’s (DOE) National Laboratories complex, the unique nature of Sandia National Laboratories’ Pulsed Power Facility (Z Machine) prohibits the use of such devices. The destructive nature of Z experiments presents a challenge to data capture and retrieval. Furthermore there are electromagnetic interference, X-ray background, and mechanical constraints to consider. Thus, a multi-part X-ray diagnostic for use on the Z Machine and Z-Beamlet Laser system has been designed and analyzed. Portions of this new DYnamic SCintillator Optic (DYSCO) have been built, tested and fielded. A data analysis software has been written. Finally, the radiance profile of the DYSCO’s scintillator has been characterized through experiments performed at the University of Arizona.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High-Fidelity Simulations of a Mach 9 Shock-Turbulent Boundary Layer Interaction

High-fidelity large-eddy simulations (LES) and wall-modeled large-eddy simulations (WMLES) are performed on the interaction of a Mach 9 turbulent boundary layer with a 34-degree compression ramp. Our primary interest is the sensitivity of mean and unsteady wall quantities to grid density. Additionally within the context of WMLES, we seek to explore solution sensitivity to the wall-model exchange location.

42 ENGINEERING↗
DOE OSTI2012

X-Ray Diffraction on NIF

The National Ignition Facility (NIF) is currently a 192 beam, 1.6 MJ laser. NIF Ramp-Compression Experiments have already made the relevant exo-planet pressure range from 1 to 50 Mbar accessible. We Proposed to Study Carbon Phases by X-Ray Diffraction on NIF. Just a few years ago, ultra-high pressure phase diagrams for materials were very 'simple'. New experiments and theories point out surprising and decidedly complex behavior at the highest pressures considered. High pressures phases of aluminum are also predicted to be complex. Recent metadynamics survey of carbon proposed a dynamic pathway among multiple phases. We need to develop diagnostics and techniques to explore this new regime of highly compressed matter science. X-Ray Diffraction - Understand the phase diagram/EOS/strength/texture of materials to 10's of Mbar. Strategy and physics goals: (1) Powder diffraction; (2) Begin with diamond; (3) Continue with metals etc.; (4) Explore phase diagrams; (5) Develop liquid diffraction; and (6) Reduce background/improve resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

KIVA-Update June 2012

Development of a fractional step, a Predictor-Corrector Split (PCS), or what is often known as a projection method combining hp-adaptive system in a Finite Element Method (FEM) for combustion modeling has been achieved. This model will advance the accuracy and range of applicability of the KIVA combustion model and software used typically for internal combustion engine modeling. This abstract describes a PCS hp-adaptive FEM model for turbulent reactive flow spanning all velocity regimes and fluids that is being developed for the new KIVA combustion algorithm, particularly for internal combustion engines. The method and general solver is applicable to Newtonian and non- Newtonian fluids and also for incompressible solids, porous media, solidification modeling, and fluid structure interaction problems. The fuel injection and injector modeling could easily benefit from the capability of solving the fluid structure interaction problem in an injector, helping to understand cycle to cycle variation and cavitation. This is just one example where the new algorithm differs from the old, in addition to handling Conjugate Heat Transfer (CHT), although there a numerous features that makes the new system more robust and accurate. In these ways, the PCS hp-adaptive algorithm does not compete with commercial software packages, those often used in conjunction with the currently distributed KIVA codes for engine combustion modeling. In addition, choosing a local ALE method on immersed moving parts represented by overset grid that is 2nd order spatially accurate, allows for easy grid generation from CAD to fluid grid while also provide for robustness in handling any possible moving parts configuration without any code modifications. The combined methods employed produce a minimal amount of computational effort as compared to fully resolved grids at the same accuracy. We demonstrate the solver on benchmark problems for the all flow regimes as follows: (1) 2-D backward-facing step using h-adaption, (2) 2-D driven cavity, (3) 2-D natural convection in a differentially heat cavity with h-adaptation, (4) NACA 0012 airfoil in 2-D, (5) supersonic flows over compression ramps, (6) 2-D natural convection in a differentially heat cavity with hp-adaptation, (7) 3-D natural convection in a differentially heat sphere with hp-adaptation. In addition, we show the new moving parts algorithm for working for a 2-D piston; the immersed moving parts method also for valves and pistons, vanes, etc... The movement is performed using an overset grid method and is 2nd order accurate in space, and never produces a tangle grid, that is, robust system at any resolution and any parts configuration. We also show CHT for the currently distributed KIVA-4mpi software and some fairly automatic grid generation using Sandia's Cubit unstructured grid generator. A new electronic web-based manual for KIVA-4 has been developed as well.

02 PETROLEUM↗