Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock physics”

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 577 records · Page 32

Interstellar Dust: Physical Processes

Dust is formed in stellar environments, and destroyed by sputtering, shattering and vaporization in shock waves due to cloud-cloud collisions and supernova blast waves. Dust is also destroyed during star formation. We review the dust formation and destruction balance. The calculated destruction time-scale is less than or equal to one billion years and the star dust injection time-scale is approx. 2.5 billion years. Hence, the fractions of elemental carbon and silicon locked up in stardust are less than 0.3 and less than 0.15, respectively. An efficient ISM dust formation route is therefore implied. In particular, in dense clouds dust grows; through the processes of coagulation and the accretion of gas phase molecules e.g. H20, CO, CH4. These icy materials may then be photoprocessed to refractory materials in more diffuse regions. The resulting carbonaceous grain mantle may actually be the glue that holds the coagulated grains together.

Jones, A. P.↗

Novel method to control explosive shock sensitivity: A mesoscale study to understand the effect of thermally expandable microsphere (TEM) inclusions in high explosives (HE) microstructure

When the void content and/or void structure of a high explosive (HE) is altered by some means (i.e., bulk heating or mechanical damage), the shock initiation behavior of the material changes. The ability to precisely predict the change in shock sensitivity after an HE has undergone microstructural changes is a crucial capability in multi-scale reactive flow models. Here, we utilize thermally expandable microspheres (TEMs) as a dopant in a polymer bonded explosive (PBX) matrix to alter the shock initiation properties in a controlled fashion. Here, using a mesoscale modeling approach, we evaluated how a single TEM (before and after thermal expansion) behaves under shock compression, as well as how the matrix PBX in the direct vicinity of the TEM is affected. We first examined the effect of an unexpanded TEM in the explosive matrix and found that its presence does not significantly perturb the bulk flow and by extension will not affect bulk sensitivity. Next, we examined the effect of an expanded TEM and found that its presence significantly perturbs the flow via hydrodynamic jetting, which causes a secondary shock wave with a strength that exceeds that of the incident wave. Finally, we showed that this secondary shock interacts with the downstream porosity to ignite a larger fraction of the overall pore volume, commensurate with the secondary shock strength and the affected volume, increasing the global (bulk) shock sensitivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Theory of discrete wave packets in the solar wind.

Discrete wave packets were observed by Ogo 5 and earlier satellites. These waves were believed to be in the whistler mode. Since their group velocities were found to be smaller than the solar-wind speed, these waves could not have been generated in the bow shock and could not have propagated upstream later. The present theory discusses a mechanism similar to that of the echo phenomenon in plasma physics discovered in recent years. The present theory enables us to explain (a) why the wave packets were associated with the bow shock, (b) why the wave packets were characterized by coherent oscillations, and (c) why the wave packets had group velocities smaller than the solar wind and yet could still occur in the solar wind. In short, our theory is able to interpret all the essential features deduced from the observational data.

Wu, C. S.↗

An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa

Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.

97 MATHEMATICS AND COMPUTING↗

Intermittent Electron-only Reconnection at Lunar Mini-magnetospheres

Abstract Lunar crustal magnetic anomalies (LCMA) are sub-ion-gyroradius structures that have been shown to stand off the solar wind (SW) plasma from the Moon’s surface, forming shock-like discontinuities and reflecting incident SW protons. In this Letter, the results of high-resolution, 2D fully kinetic simulations show a bursty electron-only magnetic reconnection in the SW-LCMA interaction region, characterized by the quasiperiodic formation and ejection of magnetic islands and strong parallel electron flows along the X-point separator lines. The islands are observed to modify the magnetic pressure pileup and Hall electric field above the LCMA, leading to sharp increases in reflected protons that drive electromagnetic fluctuations downstream and short distances upstream in the SW.

79 ASTRONOMY AND ASTROPHYSICS↗

Euler equations and the Sod shock tube problem

The Euler equations are a subset of the magnetohydrodynamic (MHD) equations in the infinitely collisional, unmagnetized limit. MHD modeling is central to many areas of plasma physics, ranging from low-temperature glow discharges to inertial confinement fusion. An important aspect of the Euler equations is their ability to describe states with discontinuities, such as shock waves. A standard benchmark test for numerical implementation of the Euler equations is the Sod shock tube. In this test, the system is initialized at rest with a pressure and density discontinuity, which results in a shock wave traveling into the low-pressure region and a rarefaction wave traveling into the high-pressure region. Starting with the presentation of the Euler equations, a numerical algorithm is presented here to solve these equations in one dimension. This is followed by an overview of the Sod shock tube problem that includes the precise initial setup and the analytic solution. Finally, the analytic solution is compared with results from numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling Hydrodynamic Instabilities, Shocks, and Radiation Waves in High Energy Density Experiments [Dissertation]

This thesis presents the computational design, modeling, and analysis of three experiments in high-energy density physics (HEDP), all of them concerning fundamental radiation flows. The first experiment is a laboratory astrophysics experiment to investigate the role of the Kelvin-Helmholtz instability (KHI) in the process of galactic filaments supplying gas to galactic halos. The achieved goal was to provide a first study in which the role of the instability is maximal and predict behavior in future iterations of an experiment accessing a more radiative regime where the role of the instability is stifled. This experiment would help answer how certain galaxies are able to grow so rapidly and produce many stars as the KHI limits this process. The second experiment, COAX, is a radiation flow experiment with a novel spectroscopy diagnostic configuration, designed to spatially measure the temperature of a radiation wave as it travels down a doped foam. A key result of this work was the development of a synthetic spectroscopy application and application of modern spectroscopy comparison techniques to provide our first temperature reconstructions from the experimental data. This experimental platform serves as the launching ground for a number of new experiments that vary the basic premise and thus is foundational to our ongoing research. The final experiment is a full integration of modeling, design, and theoretical development for the Radishock experiment. This experimental platform studies the head-on collision of a radiation wave with a counter-propagating shock, and like COAX, uses spectroscopy to diagnose and detect the interaction. My research analyzes the successful shots, indicating aspects of successful detections and suggests improvements to future iterations of the design.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Molecular Viscosity and Diffusivity Effects in Transitional and Shock-Driven Mixing Flows

Here, this paper investigates the importance of molecular viscosity and diffusivity for the prediction of transitional and shock-driven mixing flows featuring high and low Reynolds and Mach number regions. Two representative problems are computed with implicit large-eddy simulations using the inviscid Euler equations (EE) and viscous Navier-Stokes equations (NSE): the Taylor-Green vortex at Reynolds number Re = 3000 and initial Mach number Ma = 0.28, and an air- SF 6 -air gas curtain subjected to two shock waves at Ma = 1.2. The primary focus is on differences between NSE and EE predictions due to viscous effects. The outcome of the paper illustrates the advantages of utilizing NSE. In contrast to the EE, where the effective viscosity decreases upon grid refinement, NSE predictions can be assessed for simulations of flows with transition to turbulence at prescribed constant Re. The NSE can achieve better agreement between solutions and reference data, and the results converge upon grid refinement. On the other hand, the EE predictions do not converge with grid refinement, and can only exhibit similarities with the NSE results at coarse grid resolutions. We also investigate the effect of viscous effects on the dynamics of the coherent and turbulent fields, as well as on the mechanisms contributing to the production and diffusion of vorticity. The results show that nominally inviscid calculations can exhibit significantly varying flow dynamics driven by changing effective resolution-dependent Reynolds number, and highlight the role of viscous processes affecting the vorticity field. These tendencies become more pronounced upon grid refinement. The discussion of the results concludes with the assessment of the computational cost of inviscid and viscous computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

(abstract) 3D Electromagnetic Plasma Particle Simulations

A 3D electromagnetic plasma particle-in-cell code has been developed using the General Concurrent PIC algorithm. The GCPIC algorithm uses a domain decomposition to divide the computation among the processors. Particles must be exchanged between processors as they move. The efficiencies for 1-, 2-, and 3-dimensional partitions of the three dimensional domain are compared, and the algorithm is found to be very efficient even when a large fraction (e.g., 30%) of the particles must be exchanged at every time step. This PIC code has been used to perform simulations of a variety of space plasma physics problems. Results of three applications will be discussed: 1) plasma disturbances induced by moving conducting bodies in a magnetized plasma; 2) plasma plume interactions; and 3) solar wind termination shock.

electromagnetic plasma particles simulations 3D al↗

The Magnitude and Waveform of Shock Waves Induced by X-ray Lasers in Water

The high energy densities deposited in materials by focused X-ray laser pulses generate shock waves which travel away from the irradiated region, and can generate complex wave patterns or induce phase changes. We determined the time-pressure histories of shocks induced by X-ray laser pulses in liquid water microdrops, by measuring the surface velocity of the microdrops from images recorded during the reflection of the shock at the surface. Measurements were made with ~30 µm diameter droplets using 10 keV X-rays, for X-ray pulse energies that deposited linear energy densities from 3.5 to 120 mJ/m; measurements were also made with ~60 µm diameter drops for a narrower energy range. At a distance of 15 µm from the X-ray beam, the peak shock pressures ranged from 44 to 472 MPa, and the corresponding time-pressure histories of the shocks had a fast quasi-exponential decay with positive pressure durations estimated to range from 2 to 5 ns. Knowledge of the amplitude and waveform of the shock waves enables accurate modeling of shock propagation and experiment designs that either maximize or minimize the effect of shocks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effects of Void Morphology on Detonation Initiation of PBX 9502

Accurate predictive modeling of high explosives in abnormal conditions, such as a fire, is critical to personnel safety. Modeling depends heavily on precise physical characterization of the high explosive in question. The aim of this thesis was to characterize the detonation response to shock and the microstructure at high-temperature (250°C) of PBX 9502 and use these data to attempt to model these results using a new reactive burn model called SURF (Scaled Uniform Reactive Front). In order to characterize the shock response of PBX 9502 at high-temperature a 1D gas gun experiment was designed and executed that provided shock to detonation transition data. These data were then used to do a preliminary calibration of SURF for PBX 9502 at 250°C. Small-angle neutron scattering (SANS) was employed to characterize the void morphology of PBX 9502 as a function of temperature. This was done in order to feed a modified version of SURF called physically informed - SURF (π-SURF) which uses the void size and quantity distribution data as measured by SANS along with the deflagration and detonation characteristics of a high explosive to inform the SURF model. The π-SURF model was then used to model both ambient and high-temperature PBX 9502 SDT response and the results were compared to the physically measured results. It was found that π-SURF works in ambient cases. At high-temperature, however, it was determined that not enough information is known about hot spot formation to use the π-SURF model for high-temperature scenarios.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Direct imaging of shock wave splitting in diamond at Mbar pressure

Understanding the behavior of matter at extreme pressures of the order of a megabar (Mbar) is essential to gain insight into various physical phenomena at macroscales—the formation of planets, young stars, and the cores of super-Earths, and at microscales—damage to ceramic materials and high-pressure plastic transformation and phase transitions in solids. Under dynamic compression of solids up to Mbar pressures, even a solid with high strength exhibits plastic properties, causing the induced shock wave to split in two: an elastic precursor and a plastic shock wave. This phenomenon is described by theoretical models based on indirect measurements of material response. The advent of x-ray free-electron lasers (XFELs) has made it possible to use their ultrashort pulses for direct observations of the propagation of shock waves in solid materials by the method of phase-contrast radiography. However, there is still a lack of comprehensive data for verification of theoretical models of different solids. Here, we present the results of an experiment in which the evolution of the coupled elastic–plastic wave structure in diamond was directly observed and studied with submicrometer spatial resolution, using the unique capabilities of the x-ray free-electron laser (XFEL). The direct measurements allowed, for the first time, the fitting and validation of the 2D failure model for diamond in the range of several Mbar. Our experimental approach opens new possibilities for the direct verification and construction of equations of state of matter in the ultra-high-stress range, which are relevant to solving a variety of problems in high-energy-density physics.

36 MATERIALS SCIENCE↗

Swept shock/boundary-layer interactions: Scaling laws, flowfield structure, and experimental methods

A general review is given of several decades of research on the scaling laws and flowfield structures of swept shock wave/turbulent boundary layer interactions. Attention is further restricted to the experimental study and physical understanding of the steady-state aspects of these flows. The interaction produced by a sharp, upright fin mounted on a flat plate is taken as an archetype. An overall framework of quasiconical symmetry describing such interactions is first developed. Boundary-layer separation, the interaction footprint, Mach number scaling, and Reynolds number scaling are then considered, followed by a discussion of the quasiconical similarity of interactions produced by geometrically-dissimilar shock generators. The detailed structure of these interaction flowfields is next reviewed, and is illustrated by both qualitative visualizations and quantitative flow images in the quasiconical framework. Finally, the experimental techniques used to investigate such flows are reviewed, with emphasis on modern non-intrusive optical flow diagnostics.

Settles, Gary S.↗

Mg II line profiles of the Mira S Carinae

A series of high-dispersion IUE observations obtained to investigate the evolution of the shock structure of the Mira variable S Carinae have produced (despite limited phase coverage) a set of five exceptionally interesting spectra of the Mg II h and k lines. The two primary findings of these observations are (1) there is significant emission from both the h and k lines at positions corresponding to velocities of -150 km/s relative to the stellar photosphere and (2) the k/h flux ratio of the Mg II doublet remains below the theoretically predicted values which lie in the range 2:1 to 1:1 and shows a smooth dependence on the optical phase. Archival studies of other Miras (e.g., R Car) indicate that S Car is not unique in possessing unusual and highly variable Mg II h and k line profiles. Several possible physical processes that may explain these results, including radiative transfer effects in a shocked atmosphere, circumstellar scattering, and circumstellar absorption, are discussed.

Bookbinder, Jay A.↗

Using simultaneous x-ray diffraction and velocity interferometry to determine material strength in shock-compressed diamond

We determine the strength of laser shock-compressed polycrystalline diamond at stresses above the Hugoniot elastic limit using a technique combining x-ray diffraction from the Linac Coherent Light Source with velocity interferometry. X-ray diffraction is used to measure lattice strains, and velocity interferometry is used to infer shock and particle velocities. These measurements, combined with density-dependent elastic constants calculated using density functional theory, enable determination of material strength above the Hugoniot elastic limit. Furthermore, our results indicate that diamond retains approximately 20 GPa of strength at longitudinal stresses of 150–300 GPa under shock compression.

36 MATERIALS SCIENCE↗

ArcjetCV: a new machine learning application for extracting time-resolved recession measurements from arc jet test videos

Arc jet Computer Vision (ArcjetCV) is a software application built to automate analysis of arc jet ground test video footage. This includes tracking material recession, sting arm motion, and the shock-material standoff distance. This provides a new capability to resolve and validate new physics associated with non-linear processes. This is an essential step to reduce testing, modeling, and validation uncertainties for heatshield material performance. ArcjetCV uses several types of machine learning (convolutional neural net: CNN, decision tree: DT, k-means unsupervised clustering: KM) to automate the video processing pipeline. These include inferring the start/stop of time segments of interest (1D CNN), measuring the time-dependent 2D recession of the material samples (2D CNN, DT), measuring the time-dependent shock standoff distance (2D CNN, DT), and post-processing cleaning of the recession data (KM). The software also provides a graphical user interface for ease of use. The results of using this tool on arc jet videos show non-linear time-dependent effects can be important for certain materials.

machine learning↗

ArcjetCV: A New Machine Learning Application for Extracting Time-Resolved Recession Measurements From Arc Jet Test Videos

Arc jet Computer Vision (ArcjetCV) is a software application built to automate analysis of arc jet ground test video footage. This includes tracking material recession, sting arm motion, and the shock-material standoff distance. This provides a new capability to resolve and validate new physics associated with non-linear processes. This is an essential step to reduce testing, modeling, and validation uncertainties for heatshield material performance. ArcjetCV uses several types of machine learning (convolutional neural net: CNN, decision tree: DT, k-means unsupervised clustering: KM) to automate the video processing pipeline. These include inferring the start/stop of time segments of interest (1D CNN), measuring the time-dependent 2D recession of the material samples (2D CNN, DT), measuring the time-dependent shock standoff distance (2D CNN, DT), and post-processing cleaning of the recession data (KM). The software also provides a graphical user interface for ease of use. The results of using this tool on arc jet videos show non-linear time-dependent effects can be important for certain materials.

machine learning↗

ArcjetCV: A New Machine Learning Application for Extracting Time-Resolved Recession Measurements From Arc Jet Test Videos

Arc jet Computer Vision (ArcjetCV) is a software application built to automate analysis of arc jet ground test video footage. This includes tracking material recession and the shock-material standoff distance. This provides a new capability to resolve and validate new physics associated with non-linear processes. This is an essential step to reduce testing, modeling, and validation uncertainties for heatshield material performance. ArcjetCV uses several types of machine learning (convolutional neural net: CNN, decision tree: DT, k-means unsupervised clustering: KM) to automate the video processing pipeline. These include inferring the start/stop of time segments of interest (1D CNN), measuring the time-dependent 2D recession of the material samples (2D CNN, DT), measuring the time-dependent shock standoff distance (2D CNN, DT), and post-processing cleaning of the recession data (KM). The software also provides a graphical user interface for ease of use. The results of using this tool on arc jet videos show non-linear time-dependent effects can be important for certain materials and characterizing certain failure modes.

machine learning↗