Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shock structure”

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 1,603 records · Page 89

The 3D Structure of Eta Carinae's Nebula: A Definitive Picture from High-Dispersion Near-IR Spectra

High resolution long-slit spectra obtained with the Phoenix spectrograph on Gemini South provide our most accurate probe of the 3D structure of the Homunculus Nebula around Eta Carinae. Emission from molecular hydrogen at 2.122 microns traces a very thin outer skin, which contains the vast majority of the more than 10 solar masses of material in the nebula. This emission, in turn, yields our first definitive picture of the exact shape of the nebula, plus the latitude dependence of the mass-loss rate, speed, kinetic energy, shell thickness, and other properties associated with Eta Car's 19th century explosion. This will be critical for testing any models for the outburst mechanism. A preliminary analysis suggests that explosion from a critically rotating star was the dominant mechanism that gave rise to both the bipolar shape of the nebula and the production of its equatorial disk. [Fe II] emission in the near IR traces a geometrically thicker but less massive shell found on the inner surface of the H2 skin --- this is either a reverse shock that decelerates Eta Car's wind or a warm PDR. [Fe Ill emission also clarifies the structure of an inner "Little Homunculus" seen previously in HST/STlS spectra. Comparing these two tracers of cool molecular gas and warm partially-ionized gas resolves some significant confusion about the complex structure noted in previous studies.

Smith, N.↗

Mechanics of pore array collapse and interaction in shock-compressed polymethyl methacrylate (PMMA)

Recent studies on dynamic pore collapse have revealed significant development of shear localization, which can lead to material failure in porous structures and hot spot generation in energetic materials. These findings have dramatically improved the understanding of failure mechanisms during pore collapse but also prompt further investigation of realistic porous materials. In particular, porous media consist of many pores and porous networks. Even in low-porosity materials, pores can form in close proximity during the manufacturing process, leading to the critical question of pore–pore interaction during collapse under dynamic loading conditions. This study investigates, via plate impact experiments coupled with high-speed internal digital image correlation and shadowgraphy techniques, the collapse of two pores in shock-compressed PMMA at stresses between 0.4 and 1 GPa. The results of these experiments provide new insights into shear localization in pore collapse, in addition to distinct interactions between pores. Shadowgraphy measurements reveal novel, direct visualization of shear band development and crack evolution from pore surfaces. Spacing between adiabatic shear bands is measured over a range of impact stresses and is predicted accurately by the Grady–Kipp model. Pore interactions are found to effect a transition in the impact stress threshold at which different failure mechanisms initiate and are also found to possibly influence preferential sites for shear cracking. Throughout the study, numerical and theoretical models are leveraged to understand shear localization behavior. The role of baroclinicity and wave interactions between the pores is used to elucidate interaction mechanisms between pores.

Lawlor, Barry P. [California Institute of Technolo↗

Molecular dynamics simulations of inelastic x-ray scattering from shocked copper

By taking the spatial and temporal Fourier transforms of the coordinates of the atoms in molecular dynamics simulations conducted using an embedded-atom-method potential, we calculate the inelastic scattering of x rays from copper single crystals shocked along [001] to pressures of up to 70 GPa. Above the Hugoniot elastic limit, we find that the copious stacking faults generated at the shock front introduce strong quasi-elastic scattering (QES) that competes with the inelastic scattering signal, which remains discernible within the first Brillouin zone; for specific directions in reciprocal space outside the first zone, the QES dominates the inelastic signal overwhelmingly. The synthetic scattering spectra we generate from our Fourier transforms suggest that energy resolutions of order 10 meV would be required to distinguish inelastic from quasi-elastic scattering within the first Brillouin zone of shock-loaded copper. We further note that high-resolution inelastic scattering also affords the possibility of directly measuring particle velocities via the Doppler shift. These simulations are of relevance to future planned inelastic scattering experiments at x-ray Free Electron Laser facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Magnetospheric structure.

Magnetospheric research over the period from 1969 to 1971 is surveyed by describing significant results obtained from studies of the bow shock, magnetosheath, magnetopause, polar cusp, convection-driving electric fields, and the composition of trapped particles with charge greater than one. Particle and field observations are related to general questions of particle origin and destination as well as the mechanism by which the entire system remains in equilibrium.

Kavanagh, L. D., Jr.↗

the role of shock waves in modulation of galactic cosmic rays

The understanding of modulation of the galactic cosmic rays has considerably progressed by the exploration by space probes of major heliospheric structures, such as the Corotating Interaction Regions, the neutral sheet, and the compression regions of intense heliospheric magnetic fields. Also relevant in this context were the detections in the outer heliosphere of long lasting Forbush type decreases of cosmic ray intensity. The results of recent theoretical studies on the changes in intensity and energy, at different location from the Sun, induced by the passage of shocks across the heliosphere are presented. In this version of the research, the simplest cases of modulation of uGV and 2GV particles by single or several shocks during periods of positive and negative solar field polarity are reviewed. The results of the theoretical aspects of the search is reported. The comparison of the theoretical predictions with space probe data allows conclusions to be drawn on the role of shocks on the modulation on both the 11 and 22 year galactic cosmic ray cycles in the outer heliosphere and on the plausibility of the models and parameters used.

Gall, R.↗

MESSENGER Observations of Mercury's Bow Shock and Magnetopause

MESSENGER'S 14 January 2008 encounter with Mercury will provide the first new observations of the solar wind interaction with this planet since the Mariner 10 flybys that took place over 30 years ago. The closest approach distance for this first MESSENGER flyby is targeted for an altitude of 200 km as compared with the 707 km and 327 km attained by Mariner 10 on 29 March 1974 and 16 March 1975, respectively. The locations of the bow shock and magnetopause boundaries observed by MESSENGER will be examined and compared against those found in the earlier Mariner 10 measurements and the predictions of theoretical models and numerical simulations. The structure of the magnetopause will be investigated for the presence of flux transfer events or other evidence of magnetic reconnection as will the more general implications of these new MESSENGER bow shock and magnetopause observations for the global solar wind interaction with Mercury.

Slavin, James A.↗

A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries

Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.

97 MATHEMATICS AND COMPUTING↗

Simulation of Relativistic Shocks and Associated Radiation from Turbulent Magnetic Fields

Recent PIC simulations of relativistic electron-positron (electron-ion) jets injected into a stationary medium show that particle acceleration occurs at shocked regions. Simulations show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields and particle acceleration. These magnetic fields contribute to the electron's transverse deflection behind the shock. The jitter'' radiation from deflected electrons in turbulent magnetic fields has different properties than synchrotron radiation, which is calculated in a uniform magnetic field. This jitter radiation may be important for understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets in general, and supernova remnants. We will present detailed spectra for conditions relevant of various astrophysical sites of shock formation via the Weibel instability. In particular we will discuss the application to GRBs and SNRs

Nishikawa, K.-I.↗

Simulation of Relativistic Shocks and Associated Self-Consistent Radiation

Recent PIC simulations of relativistic electron-positron (electron-ion) jets injected into a stationary medium show that particle acceleration occurs at shocked regions. Simulations show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields and particle acceleration. These magnetic fields contribute to the electron's transverse deflection behind the shock. The "jitter" radiation from deflected electrons in turbulent magnetic fields has different properties than synchrotron radiation, which is calculated in a uniform magnetic field. This jitter radiation may be important for understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets in general, and supernova remnants. We will present detailed spectra for conditions relevant of various astrophysical sites of shock formation via the Weibel instability. In particular we will discuss the application to GRBs and SNRs.

Nishikawa, K.-I.↗

Unsteady shock-vortex interaction on a flexible delta wing

Unsteady Navier-Stokes computations have been carried out for simulating transonic flows over a clipped delta wing undergoing oscillatory and ramp motions, including flexibility. The implicit upwind algorithm has been validated by comparing the solutions with experimental data for the oscillatory pitching motion cases. The numerical and experimental results agree well at moderate angles of attack, where a leading-edge vortex develops. The ramp motion cases have demonstrated the effects of unsteadiness of the flow field and structural flexibility on the wing responses. For the 10 deg ramp motion, a vortex breakdown is observed. The interaction with the shock wave plays an essential role in the process of the breakdown observed in the present calculation.

Obayashi, Shigeru↗

Contact and pressure balance structures in two-fluid cosmic-ray hydrodynamics

The role of cosmic-ray-modified contact discontinuities and pressure balance structures in two-fluid cosmic-ray hydrodynamics in one Cartesian space dimension are investigated by means of analytic and numerical solution examples, as well as by weakly nonlinear asymptotics. The fundamental wave modes of the two-fluid cosmic-ray hydrodynamic equations in the long-wavelength limit consist of the backward and forward propagating cosmic-ray-modified sound waves, with sound speed dependent on both the cosmic-ray and thermal gas pressures; the contact discontinuity; and a pressure balance mode in which the sum ofthe cosmic ray and thermal gas pressure perturbations is zero. The pressure balance mode, like the contact discontinuity is advected with the background flow. The interaction of the pressure balance mode with the contact discontinuity is investigated by means of the method of multiple scales. The thermal gas and cosmic-ray pressure perturbations satisfy a linear diffusion equation, and entropy perturbations arising from nonisentropic initial conditions for the thermal gas are frozen into the fluid. The contact discontinuity and pressure balance eigenmodes both admit nonzero perturbations in the thermal gas, whereas the cosmic-ray-modified sound waves are isentropic. The total entropy perturbation is shared between the contact discontinuity and pressure balance eigenmodes, and examples are given in which there is a transfer of entropy between the two modes. In particular, N-wave type density disturbances are obtained which arise as a result of the entropy transfer between the two modes. A weakly nonlinear geometric optics perturbation expansion is used to study the long timescale evolution of the short-wavelength entropy wave and the thermal gas sound waves in a slowly varying, large-scale background flow. The weakly nonlinear geometric optics expansion is also used to generalize previous studies of squeezing instability for short-wavelength sound waves in the two fluid model, by including a weakly nonlinear wave steepening term that leads to shock formation, as well as the effect of long time and space dependence of the background flow. Implications of cosmic-ray-modified pressure balance structures and contact discontinuities in models of the interaction of traveling interplanetary shocks and compression and rarefraction waves with the solar wind termination shock are briefly discussed.

Webb, G. M.↗

AWSD Reactive Burn Model for the HMX‐Based High Explosive LX‐04

An Arrhenius–Wescott–Stewart–Davis (AWSD) reactive burn model is applied to describe shock initiation and detonation properties of the HMX-based high explosive LX-04. The parameters in the model are calibrated to data from multiple sources. The thermodynamic equations of state used in the model are calibrated to a combination of thermochemical calculations for HMX and LX-04 as well as experimentally-measured cylinder expansion results for LX-04. The kinetic parameters are calibrated to velocity data from gas gun experiments performed using EDC-32—a high explosive with the same chemical composition as LX-04 but different structural properties, and scaled rate stick data for PBX 9012. The AWSD model is shown to accurately describe the shock initiation and propagation of LX-04. Very good agreement is observed between the available experimental data and the AWSD model output. The presented results constitute an accurate LX-04 reactive burn model for use in engineering-scale models and simulations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

New Constraints on the Melting Temperature and Phase Stability of Shocked Iron up to 270 GPa Probed by Ultrafast X-Ray Absorption Spectroscopy

Studying the properties and phase diagram of iron at high-pressure and high-temperature conditions has relevant implications for Earth’s inner structure and dynamics and the temperature of the inner core boundary (ICB) at 330 GPa. Also, a hexagonal-closed packed to body-centered cubic (bcc) phase transition has been predicted by many theoretical works but observed only in a few experiments. The recent coupling of high-power laser with advanced x-ray sources from synchrotrons allows for novel approaches to address these issues. Here, we present a study on shock compressed iron up to 270 GPa and 5800 K probed by single-pulse (100 ps FWHM) x-ray absorption spectroscopy (XAS). Based on the analysis of the XAS spectra, we provide structural identification and bulk temperature measurements along the Hugoniot up to the melting. These results rule out the predicted transition to a high-temperature bcc phase and allow one to discriminate among existing equations of state models and melting curves. In particular, we report the first bulk temperature measurement in shock compressed iron on the melting plateau at 240(20) GPa and 5345(600) K. The melting curve resulting from our work extrapolates to a temperature of 6202(514) K at 330 GPa and represents a refined upper bound for the ICB temperature. Published by the American Physical Society 2024

Balugani, S.↗

Distributed Aerodynamic Sensing and Processing Toolbox

A Distributed Aerodynamic Sensing and Processing (DASP) toolbox was designed and fabricated for flight test applications with an Aerostructures Test Wing (ATW) mounted under the fuselage of an F-15B on the Flight Test Fixture (FTF). DASP monitors and processes the aerodynamics with the structural dynamics using nonintrusive, surface-mounted, hot-film sensing. This aerodynamic measurement tool benefits programs devoted to static/dynamic load alleviation, body freedom flutter suppression, buffet control, improvement of aerodynamic efficiency through cruise control, supersonic wave drag reduction through shock control, etc. This DASP toolbox measures local and global unsteady aerodynamic load distribution with distributed sensing. It determines correlation between aerodynamic observables (aero forces) and structural dynamics, and allows control authority increase through aeroelastic shaping and active flow control. It offers improvements in flutter suppression and, in particular, body freedom flutter suppression, as well as aerodynamic performance of wings for increased range/endurance of manned/ unmanned flight vehicles. Other improvements include inlet performance with closed-loop active flow control, and development and validation of advanced analytical and computational tools for unsteady aerodynamics.

Brenner, Martin↗

Cryptoblemes: A New Discovery with Major Economic Implications and Profound Changes to the Geologic Paradigm

Cryptoblemes are subtle impact shock signatures imprinted by cosmic debris on the crustal surfaces of lunar planetary bodes. These signatures constitute a complex cumulative overprinting of topographic, structural geophysical, and tectonic patterns that have a conspicuous radial centric multiringed symmetry. The geometry and distribution of cryptoblemes on Earth is comparable to the size and density of impact features on lunar planetary surfaces. Analysis of satellite imagery, sea-floor sonar, side-looking radar and aerial photographs of specific sites reveals new criteria for the identification and confirmation of impact-shock signatures. These criteria include joint and foliation patterns with asbestiform minerals, ribbon-quartz, spheroidal weathering, domal exfoliation, pencil shale, and shock spheres, which may originate from hydrocavitation of water-saturated sedimentary rocks. Cryptoblemes may also be associated with breccia pipes, sinkholes, buttes, mesas, and bogs, high-Rn anomalies, nodular concentrations, and earthquake epicenters. Major implications of cryptobleme identification include exploratory targeting of hydrocarbon and mineral deposits and the explanation of their origins. Analysis of known mineral deposits, structural traps and sedimentary basins show a direct correlation with cryptobleme patterns. Significant geologic paradigm shifts related to cryptoblemes include mountain building processes, structural orogenies, induced volcanism, earthquake origins, hydrocarbon diagenesis, formation mineral deposits, continental rifting, and plate movements, magnetic overprinting and local regional, and global geologic extinction and speciation patterns. Two figures provide a comparison between a multiring impact overprint in water and multiring cryptobleme in the U.S. basin range. (Additional information is contained in the original document).

Windolph, J., Jr.↗

Aeroacoustics of Turbulent High-Speed Jets

Aeroacoustic noise generation in a supersonic round jet is studied to understand in particular the effect of turbulence structure on the noise without numerically compromising the turbulence itself. This means that direct numerical simulations (DNS's) are needed. In order to use DNS at high enough Reynolds numbers to get sufficient turbulence structure we have decided to solve the temporal jet problem, using periodicity in the direction of the jet axis. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. Therefore in order to answer some questions about the turbulence we will partially compromise the overall structure of the jet. The first section of chapter 1 describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. In the second section we present preliminary work done using a TVD numerical scheme on a CM5. This work is only two-dimensional (plane) but shows very interesting results, including weak shock waves. However this is a nonviscous computation and the method resolves the shocks by adding extra numerical dissipation where the gradients are large. One wonders whether the extra dissipation would influence small turbulent structures like small intense vortices. The second chapter is an extensive discussion of preliminary numerical work using the spectral method to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which are solved in O(N) steps. A very accurate highly resolved DNS of a turbulent jet flow is expected.

Rao, Ram Mohan↗

Deciphering Extreme Mineral Records; Microstructural Phase Heritage of Shocked Materials

High-pressure minerals such as coesite, stishovite, and reidite are often used as indicators of shock metamorphism. However, high-pressure and high-temperature phases are metastable at ambient conditions and, therefore, often revert to a more energetically favorable phase. During solid-state transformations, the crystallographic orientations of the stable polymorph are controlled by the transformation pathway from the parent. Microstructural orientation analysis using electron backscatter diffraction (EBSD) of the stable phase can reveal systematic intercrystalline orientation relationships (OR) diagnostic of solid-state transformations from these high-pressure or high-temperature polymorph phases (Cayron et al., 2006). Knowledge of polymorph stability fields for a given system can thus be used to infer minimum shock P and/or T conditions that rocks have experienced, which are otherwise unavailable using traditional thermobarometers such as element partition coefficients. Using orientation relationships to infer the former presence of mineral phases has been termed microstructural “phase heritage”, and has successfully been used to probe for evidence of extreme pressures and temperatures resulting from hypervelocity impacts. The reconstruction of polymorphs of ZrSiO4 (zircon, reidite) and ZrO2 (baddeleyite) have further elucidated the behavior of these phases under hypervelocity impact conditions (Timms et al., 2017a; White et al., 2018), and OR analysis of shocked monazite has identified a previously unrecognized tetragonal high-P polymorph (Erickson et al., 2019). Zircon will convert to the low tetragonal polymorph reidite at shock stresses above ~21 GPa. However, upon decompression at temperatures above 1200 °C reidite will revert to zircon (Kusaba et al., 1985). Recrystallized zircon neoblasts encased in impact melt often show systematic misorientation relationships of 90° about <110> indicative of reversion from the reidite (Fig. 1), and preserving cryptic evidence of the high-P history that is often erased by the post shock thermal spike (Cavosie et al., 2016). Microstructural OR analysis of a dissociated zircon corona composed of baddeleyite (monoclinic-ZrO2) has identified the former presence of cubic zirconia within impact melt glass from the Mistastin Lake impact structure, Labrador, CA. Based on stability fields in the ZrO2-SiO2 binary phase diagram, cubic zirconia records thermal conditions in excess of 2370 °C at ambient P, indicating that the Mistastin impact melt is the hottest rock identified on Earth’s surface (Timms et al., 2017b). Shock-deformed monazite (monoclinic La,Ce,ThPO4) containing lamellae comprised of interlocking laths has been identified from the Haughton Dome impact structure, Nunavut, Canada and Nördlinger-Ries Crater, southern Germany. The lath-structured lamellae are composed of four systematic orientation variants, which OR analyses indicate originate from a tetragonal parent whereby the [010]monoclinic aligns with either [100] or [010]tetragonal (Erickson et al., 2019). These results highlight utility of microstructural phase heritage analyses for the identification of unstable high-P or high-T polymorphs that uniquely record the extremely transient shock conditions produced by hypervelocity impacts.

Timmons Erickson↗