Diagnostic development through the decades (mostly a Los Alamos perspective)

Similar documents
Shock Wave Experiments Using PDV: Window Characterization at 1550 nm

One & Two Dimensional Diagnostics for Detonators & Boosters

Extension of the Planar Noh Problem to Aluminum, Iron, Copper, and Tungsten

Push-pull PDV analysis

A METHODOLOGY TO VALIDATE 3D ARBITRARY LAGRANGIAN EULERIAN CODES WITH APPLICATIONS TO ALEGRA

Modified Symmetry Cell Approach for Simulation of Surface Enhancement Over Large Scale Structures

Using a Relativistic Electron Beam to Generate Warm Dense Matter for Equation of State Studies

HIGH ACCURACY ISENTROPIC COMPRESSION STUDIES WITH HIGH EXPLOSIVE PULSED POWER

NAVAL POSTGRADUATE SCHOOL THESIS

An Approach to Model Plastic Deformation of Metallic Plates in Hypervelocity Impact Experiments

Szendrei's model follows the time history of the crater's radial growth, resulting in the following expression for the final crater radius (r f ):

EQUATION OF STATE FOR MODELING THE DETONATION REACTION ZONE

Semiconductor Optoelectronics Prof. M. R. Shenoy Department of Physics Indian Institute of Technology, Delhi. Lecture - 13 Band gap Engineering

Sound velocities in highly-oriented pyrolytic graphite shocked to 18. GPa: orientational order dependence and elastic instability

Use of the Z Accelerator for Isentropic and Shock Compression Studies

High pressure Hugoniot measurements using converging shocks

Theory and Detonation Products Equations of State for a New Generation of Combined Effects Explosives. Dr. Ernest L. Baker

Sensors, Signals and Noise 1 COURSE OUTLINE. Introduction Signals and Noise Filtering Sensors: Strain Gauges. Signal Recovery, 2017/2018 Strain Gauges

Influence of shockwave profile on ejection of micron-scale material from shocked Sn surfaces: An experimental study

Resolving the Velocity Vector in Two Dimensions

Broadband Laser Ranging: Optimizing the Design for Different Experiments

Los Alamos NATIONAL LABORATORY LA-UR Double Shock Initiation of the HMX Based Explosive EDC-37

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department

Available online at ScienceDirect. Procedia Engineering 103 (2015 ) Impact Compaction of a Granular Material

Underwater explosion (non-contact high-intensity and/or near-field) induced shock loading of structures

Advection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II

Good morning everyone, and welcome again to MSLs World Metrology Day celebrations.

The Hugoniot and Strength of Ultem 1000 Polyetherimide

ENZYME SCIENCE AND ENGINEERING PROF. SUBHASH CHAND DEPARTMENT OF BIOCHEMICAL ENGINEERING AND BIOTECHNOLOGY IIT DELHI LECTURE 6

Photonic Displacement Interferometer

CHAPTER 1. Introduction

Refractive-Index Measurements of LiF Ramp Compressed to 800 GPa

Constitutive Modelling of the Shock Behaviour of a Highly Porous Material

Reaching the Extremes of Planet Formation Shockwave Experiments in the Giant Impact Regime

Shock Response of Polymethylmethacrylate

Slide 1. Temperatures Light (Optoelectronics) Magnetic Fields Strain Pressure Displacement and Rotation Acceleration Electronic Sensors

MATERIAL MECHANICS, SE2126 COMPUTER LAB 2 PLASTICITY

Investigation on Mie-Grüneisen type shock Hugoniot equation of state for concrete

arxiv: v2 [cond-mat.mtrl-sci] 24 Nov 2009

Supporting Online Material for

Preliminary investigations of HE performance characterization using SWIFT

Longitudinal and lateral stress measurements in shock-loaded gabbro and granite

Physics 9 Monday, February 13, 2012

D-D FUSION NEUTRONS FROM A STRONG SPHERICAL SHOCK WAVE FOCUSED ON A DEUTERIUM BUBBLE IN WATER. Dr. Michel Laberge General Fusion Inc.

Force and Acceleration Lecture 4

Mechanical Shock Testing for LIGA Materials Characterization

Application of a Computational Glass Model to the Shock Response of Soda-Lime Glass

PHYS 643 Week 4: Compressible fluids Sound waves and shocks

High-temperature and pressure aluminum reactions in carbon dioxide rich post-detonation environments

NATURE OF SHORT HIGH AMPLITUDE PULSES IN A PERIODIC DISSIPATIVE LAMINATE METAMATERIAL. Jolla CA USA USA USA

ISENTROPIC COMPRESSION OF HIGH EXPLOSIVES WITH THE Z ACCELERATOR

Compressive Residual Stress Optimization in Laser Peening of a Curved Geometry

Chamber Development Plan and Chamber Simulation Experiments

Modelling hypervelocity impact in DYNA3D J. Campbell, R. Vignjevic College of Aeronautics, Cranfield University, Cranfield,

Laser Ablation Studies at UCSD and Plans for Time and Space Resolved Ejecta Measurements

Los Alamos National Laboratory Hydrodynamic Methods Applications and Research 1 LA-UR

Big Bang, Black Holes, No Math

SIMULATIONS OF HYPERVELOCITY IMPACTS WITH SMOOTHED PARTICLE HYDRODYNAMICS

A BURN MODEL BASED ON HEATING DUE TO SHEAR FLOW: PROOF OF PRINCIPLE CALCULATIONS. F. J. Zerilli, R. H. Guirguis, and C. S. Coffey

INFLUENCE OF INITIAL DENSITY ON THE REACTION ZONE FOR STEADY-STATE DETONATION OF HIGH EXPLOSIVES

Electricity Electrostatics Types of materials Charging an Object Van de Graaff Generator

VMS-GeoMil. Background

UNVEILING THE ULTIMATE LAWS OF NATURE: DARK MATTER, SUPERSYMMETRY, AND THE LHC. Gordon Kane, Michigan Center for Theoretical Physics Warsaw, June 2009

Engineering Note. Measuring the Polarization Rotation Angle of a Faraday Rotator Mirror with the Optical Vector Analyzer

Strength effects in an imploding cylinder with constant mass-to-explosive loading

Application of Proton Radiography to High Energy Density Research

A Numerical Study on the Effect of Explosive Reactive Armour on a Lightweight Armoured Vehicle s Hull

Note on Posted Slides

This message was made with 100% recycled electrons

THEORETICAL AND EXPERIMENTAL MESOSCALE STUDIES OF IMPACT- LOADED GRANULAR EXPLOSIVE AND SIMULANT MATERIALS

arxiv:cond-mat/ v2 [cond-mat.mtrl-sci] 26 Nov 1999

Mapping Rise Time Information with Down-Shift Analysis

P-22: Hydrodynamic and X-Ray Physics

NUMERICAL SIMULATION OF CRATERING EFFECTS IN ADOBE

An Electromagnetic Gauge Technique for Measuring Particle Velocity of Electrically Conductive Samples in Extreme Conditions

Elastic Properties of Solid Materials. Notes based on those by James Irvine at

Plasma Formation and Self-focusing in Continuum Generation

Force and Motion 20 N. Force: Net Force on 2 kg mass = N. Net Force on 3 kg mass = = N. Motion: Mass Accel. of 2 kg mass = = kg m/s 2.

Electromagnetic Theory Prof. D. K. Ghosh Department of Physics Indian Institute of Technology, Bombay

Physics of Explosions

PHYSICS. Chapter 5 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc.

The Liquid Metal Ion Source A Hot Ion Source. Jon Orloff Professor Emeritus, University of Maryland and FEI Company Hillsboro, Oregon

1. Liquid Wall Ablation 2. FLiBe Properties

We search for the ether. Next time: The ether is missing Conspiracy? I think not!

STRESS, STRAIN AND DEFORMATION OF SOLIDS

The Gram-Schmidt Process

Dynamic Initiator Experiments using X-Ray Phase Contrast Imaging at the Advanced Photon Source

( )( b + c) = ab + ac, but it can also be ( )( a) = ba + ca. Let s use the distributive property on a couple of

The numerical analysis of food processing using shock wave

The University of Melbourne Engineering Mechanics

Year 7 Science. 7C1: The Particle Model. PPA Challenge

Development of a cryogenic compact interferometric displacement sensor

Abstract Submitted for the DPP99 Meeting of The American Physical Society

ASTRO 114 Lecture Okay. What we re going to discuss today are what we call radiation laws. We ve

SINGLE AND DOUBLE SHOCK INITIATION OF EDC37

PS2 Lab 1: Measurement and Uncertainty Supplemental Fall 2013

Contents. I Introduction 1. Preface. xiii

A Brief History of the Use of the Metropolis Method at LANL in the 1950s

Physical Modeling and Simulation Rev. 2

Gravitational-Wave Astronomy - a Long Time Coming Livia Conti, for the Virgo Collaboration Fred Raab, for the LIGO Scientific Collaboration

Transcription:

DOE/NV-5946--114 Diagnostic development through the decades (mostly a Los Alamos perspective) R. S. Hixson Some material taken from: Thunder in the Mountains by John Taylor, SCCM 1983 (Santa Fe) This work was done by National Security Technologies, LLC, under Contract No. DE AC5 06NA5946 with the.s. Department of Energy and supported by the Site-Directed Research and Development Program.

Outline Introduction: Why make time-resolved measurements at interfaces? The early years: Pins and such Middle years: Time-resolved quartz gauges, capacitive gauges, etc. Then there is temperature! More lately: VISAR and PDV MPDV Simulations/ simple integral experiments Conclusions --

Introduction: Why do diagnostic development for shock physics research? At the most fundamental level, we need to measure a couple of velocities; then apply the jump conditions: ρ ( u ) 0 p u p0 (1) = 1 ρ ( u u ) mass s p0 () P P 0 = ρ 0 (u s u p0 )(u p u p0 ) momentum (3) E E 0 = 1 / (P + P 0 )(V 0 V) energy This is a system of 3 equations in 5 unknowns... P, ρ, E, u s, u p Early work focused on measuring u s and inferring u p, then using these equations to get P, E, and ρ While this does not get us all we want, it does get us the principal Hugoniot To measure other physical observables requires more effort -3-

ncertainty analyses essential; errors propagate For symmetric impact, assuming errors are random, these equations hold Take p =0.5 d ;true if density in flyer and target are equal ncertainties in ρ 0 are typically very small; immersion density But not for unsymmetric impact conditions It can be shown that for very compressible materials, error in V (or ρ) can become very large 1 0 0 + + = d d s s P P δ δ ρ δρ δ 1 0 0 0 1) ( + + = d d s s V V δ δ ρ ρ ρ δρ δ ( ) ( ) 1 0 0 + + = P S P P S P S S δ δ ρ δρ ρ δρ

What do we get from time-resolved measurements? From the earliest days, shock wave experimentalists have wanted to make time-resolved measurements: Getting wave structure from pins is hard to do! Time-resolved wave profiles at interfaces contain a wealth of information HEL Shock rise/velocity Particle velocity in peak state (with some assumptions) Evolution of tension Phase transition in compression Reverse phase transition Release isentrope/interface release states -5- Particle Velocity (mm/µs) p (km/s).0 1.6 1. 0.8 0.4 0.34 0.3 0.30 0.8 0.6 0.4 0. 0.0 0.18 0.16 0.14 0.1 0.10 0.08 0.06 0.04 0.0 0.00 Zr 0 (HP) ω à β Zr I (LP) Zr II (VLP) Zirconium Wave Profiles.5e-6 3.0e-6 3.5e-6 ω à β α à ω t (s) Rigg, Knudsen (006) 0.0 0.0 0.1 0. 0.3 0.4 0.5 0.6 0.7 Time (µs) Rigg (006 56-01-33 56-01-34 56-01-36

But, in-situ measurements would be even better Velocimetry techniques are always looking at an interface Free surface; interface between target and vacuum/air: mirror image approximation Target/window interface: Wave interactions! Many efforts to do in-situ: Manganin gauges: Piezo resistance Elecromagnetic gauges: Faraday s law, insulators But, this turns out to be difficult. Embedded gauges perturb the measurement in some ways and, thus, are not widely used So, velocimetry at interfaces has become a well-used tool for looking at wave structure -6-

A bit of history the Manhattan project (Los Alamos), etc. Shock physics research began in the 19th century: Rankine (~1870) and Hugoniot (~1887), as well described by Cheret and Johnson (Classic papers in shock compression, Springer-Verlag) This was mostly theory, but with connections to artillery problems Experimental tool to make quantitative measurements really did not exist Later, the Manhattan project (194 1945) brought the need for dynamic properties of materials into a clear focus As before, mostly theory. Experimental tools were lacking. Interestingly, Bethe and Teller wrote a research paper on shock waves and equation of state (EOS) in 1941, predating the project ( Deviations from thermal equilibrium in shock waves ) Research on making better and more reproducible explosives began early on, ~1943. A lot of effort went into this at LANL s S-site. More on this later By late 1943 there was a recognition that to move metal with explosives required some understanding of the strength of the metal! No help here from the jump conditions or hydrodynamic theory. -7-

History. Hans Bethe (T Division leader) and others worked on the underlying physics theoretically; no experimental program to speak of Data table from Birch up to 0.04 Mbar were used, and TFD theory used (above 50 Mbar?) This lead to some early publications: H. Bethe et al.: Shock hydrodynamics and blast waves (1944); Penney, Energy dissipation in plastic deformation (LA-155) (1944); Hirschfelder, Shock hydrodynamics (LA-165) (1944); Metropolis, metals EOS (LA-08) (1945), etc. By 1944 Teller et al. had published a derivation of the Mie-Gruneisen EOS (LA-144) Computer simulations (such as they were!) started around this time as well. Von Neumann (and others): Solution of hydrodynamic flow equations, concept of artificial viscosity All this should be familiar to current researchers: add in an experimental component, and that how programs are now! -8-

But even most theorists agree: Experiments are needed! Research on making better explosives lead to their use to do metals research This work began in 1944 in G-8 under Froman and Critchfield: Proposed to do experiments with contacting pins First Hugoniot measurements made for the project with in-contact (unsteady) HE drive: Elastic waves observed Early work on making free surface motion measurements using a capacitance gauge by Marshall in 1945 Lab reorganized in 1946 (!!), new experimental division was stood up (M); M-4 took over this work using pin technique Interestingly, first attempts (I m aware of) to make shock temperature measurements was 1946 47-9-

More modern research; experimental technique and diagnostic evolution After 1948, many researchers returned to academia. But the staff who remained and new hires began really pushing on the experimental program: Lab re-organized (!!!) by Bradbury in 1948: GMX division formed; experiments in GMX-6 (GMX-4 doing pin work for other reasons) By 195 the flash gap method (Walsh and Christian) replaced pins for many EOS experiments But there were still drive issues! (in-contact) And, pin techniques were still used for looking at multi-wave (elastic-plastic; phase change) structure (Minshall; iron) In 1954 Fowler/Walsh/Christian did new iron flash gap experiments; discovered an anomaly in compressibility. Minshall repeated with pins. At odds with theory and Bridgeman; turned out to be the α-ε phase change at 13 GPa. Tourmaline pressure gauge records also showed multi-wave structure Bancroft et al. (1956) -10-

Timeline Around 1958 Taylor began looking at measuring shock state temperatures The early capacitive gauge technique was later revived, and it showed that both elastic and plastic waves can have structure and be less than sharp. This demonstrated the need for good time-resolved data, especially at low stress. But, back in 1955, McQueen and Walsh (and others!) began developing the HE flyer technique; allowed much higher stresses than the in-contact method and flat top (steady) waves (but with some flyer bow) When combined with the flash gap method, a data production factory was developed. Result was the current Marsh compendium of data. -11-

HE flyers + flash gaps = lots of data PIs like Walsh and McQueen defined the directions and the experiments Firing site technicians fired the shots, collected data Data analysts reduced data, gave to PIs Most focus was on relatively high pressures where the techniques worked best By about 1956 Sandia had established a shock physics effort somewhat focused on lower pressures: Lundergan, Barker, later Asay, etc. In the early 1970s, Lynn Barker at Sandia began developing an interferometer to make particle velocity measurements: 1 st reference to a laser interferometer is in 1971

The VISAR era: ~197 present Barkers early interferometers (~1965) became mature by 197 and were named VISAR Based on the WAMI concept, following earlier work done at Lockheed This was a revolutionary invention, and soon it led to two landmark research papers: Shockwave studies of PMMA, Fused Silica, and Sapphire, Barker & Hollenbach, 1970 Shockwave study of the phase transition in iron, Barker and Hollenbach, 1974; this work reconfirmed the earlier work of Minshall in considerably more detail; demonstrated how complex shock wave structure can be! In 1979 Will Hemsing (Los Alamos) introduced the push-pull VISAR modification, which made things even better The VISAR was the state-of-the-art tool for time-resolved experiments up to the development of the PDV There were other ways to get time-resolved data: Manganin gauges Quartz gauges Inclined mirrors EM embedded gauges But really, none of these have eve been used as extensively as VISAR since the 1970s -13-

For many years VISAR was the preferred technique Some beautiful time-resolved data have been collected p (km/s) 0.5 0.4 0.3 0. -Nb (6%) 8.4/4.5/55 kbar/100kbar 0.1 0.0-0.1 -e-6 0 e-6 4e-6 6e-6 8e-6 time (S) Velocity (mm/µs) 0.8 0.7 0.6 0.5 0.4 0.3 0. 0.1 0.0-0.1 VISAR Record for Shot prad 0198 0.73-0.5 0.0 0.5 1.0 1.5.0.5 Time (µs) p (km/s) Fe Spall above and below transformation 1 Hixson 0.0e-6 4.0e-6 6.0e-6 8.0e-6 1.0e-5 1.e-5 t (s) -14-

They come in many shapes and sizes Compact VISAR Early are big, hard to align More modern VISAR s have been made by NSTec/STL that really are pretty user friendly! Fast VISAR Variable air-delay VISAR -15-

Line Velocimeter for SNL Status: Cavity design complete Parts fabrication complete Assembly in process Designing array input optics

Any issues with VISAR? Well, OK, nothing is perfect! Some issues that come to mind are: Alignment of optics, which is needed to change fringe constant Issues caused by collecting unshifted light Target surface preparation; finish somewhere between specular and diffuse; loss of light upon emergence of shock wave Data analysis; centering Lissajous figure Every once in a while, we d get a data set that seemed to defy analysis! Time response limited by cavity fill time, detector response, and digitizer response Relatively hard to make multi-channel instrument; laser power Maybe best suited for planar geometry? But when it works, data with good signal to noise can be obtained -17-

Then came PDV (Strand, et. Al RSI 77 #8 083108, 006) Not my area of expertise; but what are pros and cons? Pros (?): Relatively easy to build from telecomunications parts No Lissajous stuff; replaced by sliding Fourier transform method which may be easier (?) Relatively easy to make with many channels! Cons (?): Signal/noise ratio? Can analyze either for good time response or good s/n? Any others? This takes us to MPDV, a revolutionary invention! I think MPDV has the potential to revolutionize fundamental shock physics research methods The ability to flood a target with data points is powerful We already have some examples of this -18-

MPDV: A new tool Rod penetrator example: Lexan rod hitting steel plate A lot of data! How do we interpret it all? Work in progress mpdv -19-

MPDV can also be used on very fundamental experiments Think about being able to look at many, many points of the back surface of a sample; with or without a window Edge effects can be directly observed Variation of response with location on target determined; spall strength for example (see La Lone talk next) Improvements to shock velocity measurements; stepped targets Improvements to sound velocity measurement techniques Many, many possibilities -0-

But time-resolved data have other significant benefits Even simple experiments contain complicated data; how can we use this information? Experiments used to get fundamental data also contain some integral data; both are needed to develop a better numerical modeling capability Hydrodynamic simulations Physics models Comparison of results Integral experiments Small scale experiments -1-

We can use experiments with a range of complexity to see how well our computer simulations are doing: CTH simulations of a symmetric impact copper spall experiment reveal: Johnson-Cook does well for copper EOS gets the plastic shock state right But, spall signature is challenging Right zone sizes needed! Iron wave profile is another story: Here J-C is challenged Phase change stress (pressure?) off both in compression and release Both CTH and WONDY get the rerefaction shock; CTH even gets the timing right (sound speed) We must match the whole wave profile to have confidence in our physics models & codes!

Details on Cu spall shot Simple copper flyer plate spall experiment contains integral information EOS and tensile strength info integrated Compressive strength and EOS info integrated (multiple waves Release info is complex Hydrocode modeling of such an experiment may help to identify strengths/ weaknesses of physics models and computational technique sed CTH; models chosen: MG EOS for copper Johnson-Cook library model for strength Pmin (PFRAC) model for tensile strength ( 13 kbar) is spall strength derived from data Choice of zone size matters! VISAR -3-

Conclusions Shock compression experiments are hard to diagnose! We would like to be able to do in-situ measurements, but this is difficult Diagnostic development has been continuously pursued over the years following the Manhattan project, and arguably the greatest contribution to better understanding how dynamic deformation works has come from velocimetry I think both VISAR and PDV are great tools; my personal choice for planar impact experiments would be both! Don t throw away your VISAR s! MPDV is poised to give us even greater insights! Keep up the good work!