Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHOCK WAVES”

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 19 records

Introduction to Shock Waves and Shock Wave Research

M-9 and a number of other organizations at LANL and elsewhere study materials in dynamic processes. Often, this is described as “shock wave research,” but in reality is broader than is implied by that term. Most of our work is focused on dynamic compression and associated phenomena, but you will find a wide variety of things we do that, while related, are not simple compression of materials, but involve a much richer variety of phenomena. This tutorial will introduce some of the underlying physics involved in this work, some of the more common types of phenomena we study, and common techniques. However, the list will not be exhaustive by any means.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Chiral Shock Waves

We study the shock waves in relativistic chiral matter. We argue that the conventional Rankine-Hugoinot relations are modified due to the presence of chiral transport phenomena. We show that the entropy discontinuity in a weak shock wave is quadratic in the pressure discontinuity when the effect of chiral transport becomes sufficiently large. We also show that rarefaction shock waves, which do not exist in usual nonchiral fluids, can appear in chiral matter. The direction of shock wave propagation is found to be completely determined by the direction of the vorticity and the chirality of fermions. These features are exemplified by shock propagation in dense neutrino matter in the hydrodynamic regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

3D Burgers Equation in Relativistic Plasma in the Presence of Electron and Negative Ion Trapping: Evolution of Shock Wave

The characteristics of shock waves in a relativistic plasma in the presence of nonisothermal electrons and nonisothermal negative ions is investigated by deriving the evolution equation in terms of a modified 3D Burgers equation, or trapped 3D Burgers equation. The solution of this equation is examined analytically to study the salient characteristics of shock waves in such plasma. The nonlinear coefficient is found to have the lowest (highest) value when the negative ions move toward thermal equilibrium with a dip-shaped electron distribution (when both electrons and negative ions follow a dip-shaped distribution) for a particular value of relativistic factor, and it remains in an intermediate state when both electrons and negative ions follow a flat-topped distribution. On the other hand, the dissipative coefficient is found to decrease (increase) with increasing relativistic parameter (viscous parameter). A profound effect of the trapped state of both electrons and negative ions and the temperature ratio between positive ions and electrons (negative ions and electrons) on the structure of the shock wave is also seen. However, it has been noticed that the trapped parameter of electrons has a dominating control over the shock potential profile than the trapped parameter of negative ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structure of Shock Waves and Inelasticity in Shock-Compressed Cemented Tungsten Carbides

In this study, shock wave experiments are conducted on General Carbide cemented tungsten carbide with 3.7wt.% cobalt binder to determine its shock-induced compression behavior up to 100 GPa. The measured wave profiles indicate the cemented tungsten carbide to undergo elastic-plastic deformation during shock compression. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot Elastic Limit (HEL), an elastic-plastic ramp indicating substantial post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, the shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relation, and the longitudinal stress (σ x ) vs. specific volume (V) curve for the samples. The HEL of the material was determined to lie between 4.41 GPa and 4.58 GPa. The U s - u p relation was determined to be U s = 4.97 + 1. 457u p for particle velocities greater than 0.75 km/s. The measured plastic shock velocities for particle velocities less than 0.7 km/s were found to be larger than those predicted using the linear U s - u p Hugoniot relationship, indicating the cemented WC samples to preserve substantial shear strength in the post-yield deformation region. No phase transformation was observed up to 100 GPa.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantitative x ray phase contrast imaging of oblique shock wave–interface interactions

Oblique shock wave–interface interactions of gases and liquids have been extensively studied in shock tubes using optical methods to measure equation-of-state (EOS) parameters. However, this is difficult with solids due to their opaqueness to visible light. X ray phase contrast imaging (XPCI) has the penetrative strength to probe solids while still being sensitive to mass density and enhancing the visibility of material boundaries. We investigate the accuracy and repeatability of measuring the mean value of the average mass density (areal density divided by thickness) over region S (B S ) and flow deflection angle (θ) from XPCI images of a sample. To that end, a Hough transform-based method for measuring θ is developed. To measure B S , the XPCI image intensity probability density function (PDF) is modeled accounting for the spatial distribution of x ray energy, scintillator response, and pulse-to-pulse variation in the x ray intensity. In addition, a Monte Carlo-based algorithm for computing the B S PDF is developed. Both methods are validated on an impact-generated oblique shock wave interacting at a solid polymer-to-polymer interface. Furthermore this is accomplished through a comparison to hydrodynamic simulations using well-established EOS. Under the modeling framework for the XPCI image intensity, B S is computed with an accuracy of < 0.1% and precision of 3%–5%, while θ has an uncertainty of 0.2°, respectively. This shows that the XPCI-based model that is developed here could be an invaluable tool for high-fidelity testing of hydrodynamic models in shock polar configurations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An integrated approach to optimizing concentration shock wave electrodialysis using 2D multicell simulation and response surface models

Shock wave electrodialysis (SWED) is a highly promising technique for energy-efficient ion separation in the context of a circular economy. This paper presents a approach way of modeling and improving SWED using a two-dimensional multicell model combined with the COMSOL program and response surface methodology. The model integrates the Nernst-Planck equation, Darcy's law, and first-order electroosmosis to examine the local concentration, flux of ionic species, distribution of current, and velocity of flow in SWED cells under various operating conditions. We first illustrate the clear depiction of concentration, velocity, and electric potential distribution through contours which aids in identifying optimal operating conditions and designing scalable SWED systems. The results emphasize the significance of surface charge density and voltage in influencing the features of shock waves for obtaining effective ion separation while optimizing energy consumption and improving current efficiency by controlling the retention time of feed flow. Here, this study defines two crucial characteristics of shock waves, namely the length of the flat depletion zone of a fully developed shock wave (shock wave height) and the distance of shock wave propagation (shock wave length). These properties significantly impact separation performance, as determined by the simulation results. Additionally, the response surface methodology is incorporated with the COMSOL models to develop predictive models and graph responses, enabling a more comprehensive understanding of the interactions between parameters and performance indicators, such as removal ratio, energy consumption, and water recovery. Finally, this work suggests design tactics for expanding SWED processes and outlines potential areas for further research. This research provides valuable insights into the prospective applications, design optimization, and scalability of SWED in the field of electrokinetic separation technologies for green chemistry and a circular economy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A comparison of explosively driven shock wave radius versus time scaling approaches

Abstract Explosively driven shock wave radius versus time profiles are frequently used to document energy release and relative explosive performance. Recently, two universal shock wave radius versus time profiles have been presented in the literature, which demonstrate the ability to represent explosively driven shock wave profiles for all explosive sources in any fluid environment. These two universal shock wave profiles are examined here relative to each other and relative to a commonly used nonlinear shock wave profile, which is fit to experimental data for individual explosive materials. The nonlinear profile, originally developed by Dewey, is examined here, and a universal non-dimensional form of the equation is proposed. The universal shock wave profiles are all found to be relatively similar, but with slight variations in a transition region of non-dimensional radii $$0.15\lesssim R^*\lesssim 2$$ 0.15 ≲ R ∗ ≲ 2 . The variations in this region result in different estimations of energy release or blast strength between the curve fits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

How plane are plane shock waves in solids

The notion of plane shock waves is a macroscopic, very fruitful idealization of near discontinuous disturbance propagating at supersonic speed. Such a picture is comparable to the picture of shorelines seen from a very high altitude. When viewed at the grain scale where the structure of solids is inherently heterogeneous and stochastic, features of shock waves are non-laminar and field variables, such as particle velocity and pressure, fluctuate. This paper reviews select aspects of such fluctuating nonequilibrium features of plane shock waves in solids with focus on grain scale phenomena and raises the need for a paradigm change to achieve a deeper understanding of plane shock waves in solids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Weak-shock wave propagation in polymer-based particulate composites

Shock waves are common in polymer-based particulate composites that are subjected to intermediate to high-velocity impact loading. However, quantitative information on the spatial variation of stress, particle velocities, and energy dissipation during the formation and propagation of weak-shock waves is limited. In this paper, a systematic experimental study is conducted to understand the characteristics of weak-shocks in polymer-bonded particulate composites. Specimens made of polymer-bonded sugar were subjected to a projectile impact loading, at varying velocities, using a modified Hopkinson pressure bar apparatus. Full-field displacement and strains of the deformed samples were obtained with the help of an ultrahigh-speed imaging and digital image correlation technique. Using the full-field displacement data, the shock wave velocity, shock front thickness, and the full-field stress fields are calculated. From the spatial stress field and the strain rate data, the spatial energy dissipation profile is also estimated. The effect of impact velocity on the spatial stress profile, shock wave velocity, and energy dissipation are discussed.

42 ENGINEERING↗

Surface Properties of Double-Fin Generated Shock-Wave/Boundary-Layer Interactions

Shock-Wave/Boundary-Layer Interaction (SBLI) is ubiquitous in high-speed air vehicles. SBLI surface flowfield properties of crossing shocks induced by symmetric Double-Fins (DFs) at fin angles of 8° and 10° at Mach 2 are experimentally investigated; a 10° Single-Fin (SF) is also explored as a comparison. Surface visualization is used to capture distinct flowfield structures such as separation and upstream influence. Steady and unsteady Pressure-Sensitive Paints (PSPs) are used to obtain global surface pressure fields. Characteristics of DF-SBLI surface topology are discussed by comparing the corresponding SF-SBLI. Here, the effect of fin angle on DF-SBLI is characterized, and unsteady dynamics of the pressure field are examined. Spectral Proper Orthogonal Decomposition (SPOD) is performed on unsteady PSP, and traveling surface pressure waves and their frequency-dependent behavior are discussed. Dispersion relation of surface pressure waves is revealed by performing two-dimensional space–time Fourier transform on unsteady PSP, and the results are compared with the phase velocity distributions obtained from SPOD. Accompanying group velocity behavior is examined for distinct SBLI regions wherein wave propagation is found dependent on its frequency. Steady surface pressures along the centerline are compared with previous experimental data. Fin-tip spacing is found to be an essential parameter for scaling of upstream influence length.

33 ADVANCED PROPULSION SYSTEMS↗

Numerical study of variable density turbulence interaction with a normal shock wave

Accurate numerical simulations of shock–turbulence interaction (STI) are conducted with a hybrid monotonicity-preserving–compact-finite-difference scheme for a detailed study of STI in variable density flows. Theoretical and numerical assessments of data confirm that all turbulence scales as well as the STI are well captured by the computational method. Linear interaction approximation (LIA) convergence tests conducted with the shock-capturing simulations exhibit a similar trend of converging to LIA predictions to shock-resolving direct numerical simulations (DNS). The effects of density variations on STI are studied by comparing the results corresponding to an upstream multi-fluid mixture with the single-fluid case. Here the results show that for the current parameter ranges, the turbulence amplification by the normal shock wave is much higher and the reduction in turbulence length scales is more significant when strong density variations exist. Turbulent mixing enhancement by the shock is also increased and stronger mixing asymmetry in the postshock region is observed when there is significant density variation. The turbulence structure is strongly modified by the shock wave, with a differential distribution of turbulent statistics in regions having different densities. The dominant mechanisms behind the variable density STI are identified by analysing the transport equations for the Reynolds stresses, vorticity, normalized mass flux and density specific volume covariance.

42 ENGINEERING↗

Multiscale modeling of shock wave localization in porous energetic material

Shock wave interactions with defects, such as pores, are known to play a key role in the chemical initiation of energetic materials. The shock response of hexanitrostilbene is studied through a combination of large-scale reactive molecular dynamics and mesoscale hydrodynamic simulations. In order to extend our simulation capability at the mesoscale to include weak shock conditions (< 6 GPa), atomistic simulations of pore collapse are used here to define a strain-rate-dependent strength model. Comparing these simulation methods allows us to impose physically reasonable constraints on the mesoscale model parameters. In doing so, we have been able to study shock waves interacting with pores as a function of this viscoplastic material response. Finally, we find that the pore collapse behavior of weak shocks is characteristically different than that of strong shocks.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A numerical study of shock waves generated through laser ablation of explosives

Shock waves resulting from irradiation of energetic materials with a pulsed ultraviolet laser source have been shown to be an effective indicator for explosives detection. In this study, the features of shock wave propagation are explored theoretically. The initial stage of the shock motion is simulated as a one-dimensional process. As the nonlinear wave expands to form a blast wave, a system of conservation equations, simplified to the Euler equations, is employed to model wave propagation. The Euler equations are solved numerically by the 5th order weighted essentially non-oscillatory finite difference scheme with the time integration carried out using the 3rd order total variation diminishing Runge Kutta method. The numerical results for the shock wave evolution are compared with those obtained from experiments with a meltcast 2,6-dinitrotoluene sample. The calculations lay a theoretical foundation for a recently investigated technique for photoacoustically sensing explosives using a vibrometer.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The preplasma effect on the properties of the shock wave driven by a fast electron beam

Strong shock wave generation by a mono-energetic fast electron beam in a plasma with an increasing density profile is studied theoretically. The proposed analytical model describes the shock wave characteristics for a homogeneous plasma preceded by a low density precursor. The shock pressure and the time of shock formation depend on the ratio of the electron stopping length to the preplasma areal density and on the initial energy of injected electrons. The conclusions of theoretical model are confirmed in numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Propagation and dispersion of shock waves in magnetoelastic materials

Previous studies examining the response of magnetoelastic materials to shock waves have predominantly focused on applications involving pulsed power generation, with limited attention given to the actual wave propagation characteristics. This study provides detailed magnetic and mechanical measurements of magnetoelastic shock wave propagation and dispersion. Laser generated rarefacted shock waves exceeding 3 GPa with rise times of 10 ns were introduced to samples of the magnetoelastic material Galfenol. The resulting mechanical measurements reveal the evolution of the shock into a compressive acoustic front with lateral release waves. Importantly, the wave continues to disperse even after it has decayed into an acoustic wave, due in large part to magnetoelastic coupling. The magnetic data reveal predominantly shear wave mediated magnetoelastic coupling, and were also used to noninvasively measure the wave speed. The external magnetic field controlled a 30% increase in wave propagation speed, attributed to a 70% increase in average stiffness. Lastly, magnetic signals propagating along the sample over 20× faster than the mechanical wave were measured, indicating these materials can act as passive antennas that transmit information in response to mechanical stimuli.

36 MATERIALS SCIENCE↗

The research into the propagation law of the shock wave of a gas explosion inside a building

Based on the dissipation rate conservation equations of turbulent kinetic energy in the k‐ ε turbulence model, a complicated three‐dimensional finite element model of a kitchen filled with gas mixture is developed by using the open source field operation and manipulation (OpenFOAM). Two representative kitchens were used to investigate the propagation law of the shock wave of a gas explosion inside a building by considering the key characteristics of the blast shock wave. The influence of some crucial parameters, such as initial conditions and kitchen parameters, on the properties of the blast shock wave is investigated. The basic steps to predict the peak pressure of the blast shock wave are given in consideration of the initial condition and the kitchen whilst the injury effect of the blast shock wave on the humans and animals is evaluated. The research results indicate that the pressure time history and the peak pressure space distribution are greatly influenced by the kitchen design layout. The coupled interaction between the initial temperature and gas volume concentration, especially at the upper and lower explosion limits of the gas, significantly affects the peak pressure. The peak pressure varies significantly with the opening and the buffer; however, it has little relation with the width, length, and height of the kitchen. The proposed method can accurately and effectively predict the peak pressure of the blast shock wave inside buildings. In terms of the peak pressure space distribution of the explosion shock wave, the peak pressure is much higher than the threshold of the killing pressure, which is unsafe for the humans and animals in the building.

Lin, Shu-Chao↗