Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock capturing”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Validation of predictive performance models for supersonic gas-jet nozzles at the Laboratory for Laser Energetics

We present results characterizing the neutral-density distributions produced by the supersonic nozzles used in experiments on the OMEGA-60 and OMEGA-EP laser systems at the University of Rochester’s Laboratory for Laser Energetics (LLE). Axisymmetric Fluent® simulations using LLE nozzle specifications capture the viscous effects, gas expansion, and shock waves that complicate flow predictions for offsets above the nozzle exit. Here, these simulations show good agreement with neutral-density measurements obtained using a four-wave shearing interferometer. An analytical form is given for the plateau length. Fits to simulation data for boundary layer thickness, mean plateau density, and density ramps are given as functions of nozzle offset and nozzle backing pressure for a number of nozzles and gases.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

In Situ X-ray Diffraction of High Density Polyethylene During Dynamic Compression: Polymer Chain Compression and Phase Transformations

Polymers are widely implemented materials often used for safety and security in our daily lives. However, we do not have a thorough understanding of the phase transformations of these materials when exposed to extreme condition, high strain rate, events. To study polymers in these severe environments, we performed plate impact experiments that impart a one-dimensional shock wave, at high pressures and temperatures, to observe the materials response. Generally, velocimetry is the experimental diagnostic used to capture the movement of high velocity interfaces to then calculate the mechanical variables associated with shock compression. However, when a material undergoes a phase transition, such as between two crystalline phases or melting to decomposition during shock compression, velocimetry does not explicitly indicate what structural change occurred. At the Dynamic Compression Sector at the Advanced Photon Source gas gun driven compression has been paired with x-ray diffraction, providing the ability to study both the shock parameters and structural changes concurrently. We explored high density polyethylene (HDPE) both on its reactants Hugoniot (<10.12 GPa) and within its transition region (24.3 GPa) where HDPE is hypothesized to decompose from reactants to products. With x-ray diffraction we observed both polymer chain compression and crystal structure reorientation from orthorhombic to monoclinic above 6 GPa, as well as either a possible phase change or chemical decomposition at 24.3 GPa. With these experiments we were able to determine the crystalline and amorphous Hugoniot loci at pressures below the polymer transition region.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Mathematical modelling of the concave front in the adjacent high explosive detonation problem

This study presents an analysis of the transition-zone in adjacent high explosive (HE) detonation problems which uses a $D, 𝜅, \dot{D}$ relationship, where $D$ is the detonation front-normal velocity, 𝜅 is the detonation front curvature and $\dot{D}$ is the time derivative of detonation front-normal velocity. Our approach extends the traditional $(D, 𝜅)$ model to accurately predict the behaviour of both diverging and converging detonation shock fronts. Our findings affirm that a hyperbolic type of front evolution equation, enhanced with wave acceleration, provides a robust framework for modelling complex shock front dynamics in HE materials. This approach not only captures the natural effects of straightness and boundary slope jumps in the transition-zone but also bridges the gap between mathematical predictions and experimental observations, offering insights into the behaviour of both diverging and converging detonation propagations in a homogeneous HE.

acceleration↗
DOE OSTI2022

Spatial and Temporal Interpolation Analysis Process of Shock Loading

The following document serves to describe the current process that is utilized by W-13 analysts to approximate the pressure-time boundary condition seen by a test object from shock tube loading. The intent of this report is to capture the current capability and to encourage continued growth and development. Information on the required input variables, mapping process, and generated output information is covered in the following sections. Throughout the document, gaps and weaknesses of the mapping process are noted to encourage future development efforts.

42 ENGINEERING↗
DOE OSTI2025

Acoustic Observations of the OSIRIS-REx Sample Return Capsule Re-Entry from Wendover Airport

The Origins, Spectral Interpretation, Resource Identification, and Security‐Regolith Explorer sample return capsule (SRC) re‐entered the Earth’s atmosphere at hypersonic speeds from interplanetary space on 24 September 2023. The current work reports on 18 ground‐based acoustic sensors deployed at Wendover Airport, the same location that the Genesis and Stardust SRC re‐entries were recorded. Four different sensors (Chaparral Physics, Gem, Wilson Engineering Research and Development [WERD], and RedVox) were deployed in close proximity to compare their performance. All the sensors captured an N‐wave signal associated with the SRC re‐entry shock wave followed by a broadband coda. The Chaparral Physics array served as the high‐fidelity reference measurement. The N‐wave signal had a peak‐to‐peak amplitude of 4.07 Pa with a fundamental frequency of 4.98 Hz from 167.5° measured clockwise from north, nearly perpendicular to the SRC trajectory. In addition, high coherence in the coda was shown to be associated reflections from the surrounding mountains. In general, the more economical sensors (Gem, WERD, and RedVox) produced results that were consistent with these observations and sensor specifications. Beamforming with these single sensors arranged as an array showed agreement with the high‐fidelity array to within a couple of degrees. Furthermore, the current high‐fidelity results were compared with the measurements during the Genesis and Stardust SRC re‐entries. All three entries produced a broadband fundamental peak at a frequency that was inversely related to the SRC diameter as well as evidence of reflections from the surrounding topography.

KC, Real J. [Oklahoma State University, Stillwater↗
DOE OSTI2020

Megahertz-rate shock-wave distortion cancellation via phase conjugate digital in-line holography

Abstract Holography is a powerful tool for three-dimensional imaging. However, in explosive, supersonic, hypersonic, cavitating, or ionizing environments, shock-waves and density gradients impart phase distortions that obscure objects in the field-of-view. Capturing time-resolved information in these environments also requires ultra-high-speed acquisition. To reduce phase distortions and increase imaging rates, we introduce an ultra-high-speed phase conjugate digital in-line holography (PCDIH) technique. In this concept, a coherent beam passes through the shock-wave distortion, reflects off a phase conjugate mirror, and propagates back through the shock-wave, thereby minimizing imaging distortions from phase delays. By implementing the method using a pulse-burst laser setup at up to 5 million-frames-per-second, time-resolved holograms of ultra-fast events are now possible. This technique is applied for holographic imaging through laser-spark plasma-generated shock-waves and to enable three-dimensional tracking of explosively generated hypersonic fragments. Simulations further advance our understanding of physical processes and experiments demonstrate ultra-high-speed PCDIH techniques for capturing dynamics.

Science & Technology - Other Topics↗
DOE OSTI2022

Hybrid interferometric diagnostic for the refractive index measurement across nearly discrete shock waves

The measurement of high temperature gas properties is key for characterizing high speed flows. Nearly discrete changes in density across shock waves, in particular, are difficult to resolve through traditional fringe-counting interferometric methods. Existing techniques for estimating large fringe jumps are either resolution limited or require specialized window configurations. In this Letter, we describe a unique hybrid interferometric technique that combines narrowband fringes for high resolution and broadband fringes as an absolute reference to measure changes in refractive index with a resolution of up to 7 × 10 −8 across nearly discrete index changes of up to 1.5 × 10 −4 . By capturing fringes with an ultrahigh-speed camera, the refractive index changes across discrete shock fronts can be estimated inside a shock tube with high accuracy and time resolution. First, a novel hybrid calibration method for tracking finite fringes is discussed. Next, this technique is used to measure the post-initial-shock refractive indices for Mach 2.7 to 4.2 flows (pressures from 90.4 to 228.4 kPa). Results are then compared with theoretical values showing agreement within 2%.

Wang, Gwendolyn T. (ORCID:000000016271893X)↗
DOE OSTI2023

Post-detonation fireball modeling: Validation of freeze out approximations

A numerical simulation study was performed to examine the post-detonation reaction processes produced by the detonation of a 12 mm diameter hemispherical pentaerythritol tetranitrate (PETN) explosive charge. The simulations used a finite rate detailed chemical reaction model consisting of 59 species and 368 reactions to capture post-detonation reaction processes including air dissociation from Mach 19+ shock waves that initially break out of the PETN charge, reactions within the detonation products during expansion, and afterburning when the detonation products mix with the shock heated air. The multi-species and thermodynamically complete Becker-Kistiakowsky-Wilson real-gas equation of state is used for the gaseous phase to allow for the mixing of reactive species. A recent simplified reactive burn model is used to propagate the detonation through the charge and allow for detailed post-detonation reaction processes. The computed blast, shock structures, and mole fractions of species within the detonation products agree well with experimental measurements. A comparison of the simulation results to equilibrium calculations indicates that the assumption of a local equilibrium is fairly accurate until the detonation products rapidly cool to temperatures in the range of 1500-1900 K by expansion waves. Below this range, the computed results show mole fractions that are nearly chemically frozen within the detonation products for a significant portion of expansion. In conclusion, these results are consistent with the freeze out approximation used in the blast modeling community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Observations on the deformation of metal microspheres in shock-driven polymer flows

We report that solid particles can be fragmented by a fast-moving fluid if their velocity difference is great enough, such as during the atmospheric entry of meteoroids or the shock compression of engineered particulate composites. The extent of particle deformation and breakup in such systems is poorly understood because the necessary extreme conditions make observation difficult and data scarce. To meet this need, experiments combining ultrafast synchrotron-based radiography with plate impact loading were performed at the dynamic compression sector at the advanced photon source. Metal microspheres of several densities and strengths (Au, Ta, and W) were placed inside a polymer matrix. A planar shock wave was then produced in the polymer by the impact of a gun-launched flyer plate. X-ray images of the resulting flow were collected at ~150ns intervals. These images document the progression of particle deformation across a range of flow conditions and particle materials. They show that the extent of deformation is sensitive to the ratio of drag stress to particle strength. The deforming particle's shape is determined by the initial shock–particle interaction, fluid stagnation pressure, and vorticity, each acting on its own timescale. A set of scaling relationships is presented to capture these observations and enable comparison with prior hydrodynamic data. The result is a framework for predicting the conditions under which strong particles are severely deformed by a shock-driven flow.

36 MATERIALS SCIENCE↗
DOE OSTI2021

New Design Yields Robust Large-Area Framing Camera

Manuscript to be submitted to Review of Scientific Instruments to determine interest in possible publication. Website: https://aip.scitation.org/rsi/authors/manuscript High-speed imaging is a core diagnostic capability for the Nevada National Security Site (NNSS). Data from experiments that measure the behavior of materials under extreme conditions are essential for the stockpile mission. For example, during subcritical experiments, materials that are shocked and compressed under high-explosive drive change behavior rapidly over timescales measured in millionths or even billionths of a second. Characterizing and understanding this behavior requires radiographic and/or optical imaging to track motion on these timescales. Framing cameras provide a discrete series of images over a short period of time in a manner that mimics a high-speed movie camera, but with individual shutter times that must be extremely short in order to capture freeze-frame images of rapidly evolving phenomena such as high-explosive shock-driven ejecta, dynamic compression of metals, and high-velocity fluid flow.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗
DOE OSTI2022

An Efficient B-Spline Lagrangian/Eulerian Method for Compressible Flow, Shock Waves, and Fracturing Solids

This study presents a new method for modeling the interaction between compressible flow, shock waves, and deformable structures, emphasizing destructive dynamics. Extending advances in time-splitting compressible flow and the Material Point Methods (MPM), we develop a hybrid Eulerian and Lagrangian/Eulerian scheme for monolithic flow-structure interactions. We adopt the second-order WENO scheme to advance the continuity equation. To stably resolve deforming boundaries with sub-cell particles, we propose a blending treatment of reflective and passable boundary conditions inspired by the theory of porous media. The strongly coupled velocity-pressure system is discretized with a new mixed-order finite element formulation employing B-spline shape functions. Shock wave propagation, temperature/density-induced buoyancy effects, and topology changes in solids are unitedly captured.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Kinetic simulation study of magnetized collisionless shock formation on a terawatt laser system

Perpendicular, magnetized, collisionless shocks in hydrogen and neon plasmas are studied with 2-D particle-in-cell simulations for parameters accessible to experiments on OMEGA EP. The simulations are performed with realistic ion-electron mass ratios by which the relative importance of different micro-instabilities can be accurately captured. The dispersion relation of the modified two-stream instability (MTSI), the main dissipation mechanism for these shocks, is used to find suitable parameters for upcoming experiments. Simulations show that magnetized collisionless shocks can be readily formed within a few tenths of an ion gyro-period in both hydrogen and neon gases, with a background magnetic field of 50 T, achievable using the magneto-inertial fusion electrical discharge system (MIFEDS). Here, a portion of the ions are reflected to the upstream region and accelerated in both shock normal and tangential directions, indicating the formation of a supercritical shock. Shock front reformation is seen in longer time 1-D simulations. The results show that the formation time and width of these shocks are determined by MTSI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

New Design Yields Robust Large-Area Framing Camera

Framing cameras provide a discrete series of images over a short period of time in a manner that mimics a high-speed movie camera, but with individual shutter times that must be extremely short in order to capture freeze-frame images of rapidly evolving phenomena such as high-explosive shock-driven ejecta, dynamic compression of metals, and high-velocity fluid flow. The NNSS has designed and fielded a new, large-area, gated framing camera called Kraken. The camera design emphasized manufacturability and flexibility by improving imager yield and creating a camera architecture to shorten design cycles. Design strategies and field data are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗