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,675 records · Page 93

High-energy Emission from the Composite Supernova Remnant MSH 15-56

MSH 1556 (G326.3-1.8) is a composite supernova remnant (SNR) that consists of an SNR shell and a displaced pulsar wind nebula (PWN) in the radio. We present XMM-Newton and Chandra X-ray observations of the remnant that reveal a compact source at the tip of the radio PWN and complex structures that provide evidence for mixing of the supernova (SN) ejecta with PWN material following a reverse shock interaction. The X-ray spectra are well fitted by a non-thermal power-law model whose photon index steepens with distance from the presumed pulsar, and a thermal component with an average temperature of 0.55 keV. The enhanced abundances of silicon and sulfur in some regions, and the similar temperature and ionization timescale, suggest that much of the X-ray emission can be attributed to SN ejecta that have either been heated by the reverse shock or swept up by the PWN. We find one region with a lower temperature of 0.3 keV that appears to be in ionization equilibrium.Assuming the Sedov model, we derive a number of SNR properties, including an age of 16,500 yr. Modeling of the gamma-ray emission detected by Fermi shows that the emission may originate from the reverse shock-crushed PWN.

Temim, Tea↗

Generalized seismic analysis

There is a constant need to be able to solve for enforced motion of structures. Spacecraft need to be qualified for acceleration inputs. Truck cargoes need to be safeguarded from road mishaps. Office buildings need to withstand earthquake shocks. Marine machinery needs to be able to withstand hull shocks. All of these kinds of enforced motions are being grouped together under the heading of seismic inputs. Attempts have been made to cope with this problem over the years and they usually have ended up with some limiting or compromise conditions. The crudest approach was to limit the problem to acceleration occurring only at a base of a structure, constrained to be rigid. The analyst would assign arbitrarily outsized masses to base points. He would then calculate the magnitude of force to apply to the base mass (or masses) in order to produce the specified acceleration. He would of necessity have to sacrifice the determination of stresses in the vicinity of the base, because of the artificial nature of the input forces. The author followed the lead of John M. Biggs by using relative coordinates for a rigid base in a 1975 paper, and again in a 1981 paper . This method of relative coordinates was extended and made operational as DMAP ALTER packets to rigid formats 9, 10, 11, and 12 under contract N60921-82-C-0128. This method was presented at the twelfth NASTRAN Colloquium. Another analyst in the field developed a method that computed the forces from enforced motion then applied them as a forcing to the remaining unknowns after the knowns were partitioned off. The method was translated into DMAP ALTER's but was never made operational. All of this activity jelled into the current effort. Much thought was invested in working out ways to unshakle the analysis of enforced motions from the limitations that persisted.

Butler, Thomas G.↗

A study of X-ray and infrared emissions from dusty nonradiative shock waves

We have constructed models that predict the dynamic evolution and infrared (IR) emission of grains behind nonradiative shock waves. We present a self-consistent treatment of the effect of grain destruction and heating on the ionization structure and X-ray emission of the postshock gas. Incorporating thermal sputtering, collisional heating, and deceleration of grains in the postshock flow, we predict the IR and X-ray fluxes from the dusty plasma as a function of swept-up column density. Heavy elements such as C, O, Mg, S, Si and Fe are initially depleted from the gas phase but are gradually returned as the grains are destroyed. The injected neutral atoms require some time to 'catch up' with the ionization state of the ambient gas. The nonequilibrium ionization state and gradient in elemental abundances in the postshock flow produces characteristic X-ray signatures that can be related to the age of the shock and amount of grain destruction. We study the effects of preshock density and shock velocity on the X-ray and IR emission from the shock. We show that the effects of graindestruction on the X-ray spectra of shock waves are substantial. In particular, temperatures derived from X-ray spectra of middle-aged remnants are likely to be overestimated by approximately 15% if cosmic abundances are assumed. Due to the long timescales for grain destruction in X-ray gases over a wide range of temperatures, we suggest that future X-ray spectra studies of supernova remnants be based on depleted abundances instead of cosmic abundances. Our model predictions agree reasonably well with IRAS and Einstein IPC observations of the Cygnus Loop.

Vancura, Olaf↗

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗

Implosion dynamics of triple-nozzle gas-puff z pinches on COBRA

Experiments on the 1-MA, 220-ns COBRA generator at Cornell University were conducted to provide detailed measurements of structured cylindrical gas-puff z pinches. In the experiments, a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet initialize the z-pinch load with various working gases, radial density profiles, and externally applied axial magnetic fields. Planar laser-induced fluorescence provides a measure of the initial neutral gas density of the load, while three-frame laser shearing interferometry and multi-frame extreme ultraviolet (XUV) cameras reveal the formation and propagation of a magneto-Rayleigh–Taylor (MRT) unstable shock layer. Implosion trajectories are compared to simple, experimentally informed models and found to be in good agreement. Differences in the structure of the accelerating plasma sheath and evolution of the MRT instability are observed for different gas species and axial magnetic field strengths, correlating with differences in pinch uniformity and x-ray emission. Here, the average instability growth is compared to linear MRT theory predictions using the instantaneous acceleration of the best-fit implosion models and characteristic instability wavelength, with the effective Atwood number and seed perturbation size as fit parameters. For high density argon center jets, ionization prior to the arrival of the imploding plasma sheath suggests a heating mechanism consistent with photoionization by XUV self-emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Galactic Cosmic Ray Intensity Response to Interplanetary Coronal Mass Ejections/Magnetic Clouds in 1995-2009

We summarize the response of the galactic cosmic ray (CGR) intensity to the passage of the more than 300 interplanetary coronal mass ejections (ICMEs) and their associated shocks that passed the Earth during 1995-2009, a period that encompasses the whole of Solar Cycle 23. In approx.80% of cases, the GCR intensity decreased during the passage of these structures, i.e., a "Forbush decrease" occurred, while in approx.10% there was no significant change. In the remaining cases, the GCR intensity increased. Where there was an intensity decrease, minimum intensity was observed inside the ICME in approx.90% of these events. The observations confirm the role of both post-shock regions and ICMEs in the generation of these decreases, consistent with many previous studies, but contrary to the conclusion of Reames, Kahler, and Tylka (Astrophys. 1. Lett. 700, L199, 2009) who, from examining a subset of ICMEs with flux-rope-like magnetic fields (magnetic clouds) argued that these are "open structures" that allow free access of particles including GCRs to their interior. In fact, we find that magnetic clouds are more likely to participate in the deepest GCR decreases than ICMEs that are not magnetic clouds.

Richardson, I. G.↗

Stream dynamics between 1 AU and 2 AU: A detailed comparison of observations and theory

A radial alignment of three solar wind stream structures observed by IMP-7 and -8 (at 1.0 AU) and Voyager 1 and 2 (in the range 1.4 to 1.8 AU) in late 1977 is presented. It is demonstrated that several important aspects of the observed dynamical evolution can be both qualitatively and quantitatively described with a single-fluid 2-D MHD numerical model of quasi-steady corotating flow, including accurate prediction of: (1) the formation of a corotating shock pair at 1.75 AU in the case of a simple, quasi-steady stream; (2) the coalescence of the thermodynamic and magnetic structures associated with the compression regions of two neighboring, interacting, corotating streams; and (3) the dynamical destruction of a small (i.e., low velocity-amplitude, short spatial-scale) stream by its overtaking of a slower moving, high-density region associated with a preceding transient flow. The evolution of these flow systems is discussed in terms of the concepts of filtering and entrainment.

Burlaga, L. F.↗

Stream dynamics between 1 AU and 2 AU - A comparison of observations and theory

A radial alignment of three solar wind stream structures observed by IMP-7 and -8 (at 1.0 AU) and Voyager 1 and 2 (in the range 1.4 to 1.8 AU) in late 1977 is presented. It is demonstrated that several important aspects of the observed dynamical evolution can be both qualitatively and quantitatively described with a single-fluid 2-D MHD numerical model of quasi-steady corotating flow, including accurate prediction of: (1) the formation of a corotating shock pair at 1.75 AU in the case of a simple, quasi-steady stream; (2) the coalescence of the thermodynamic and magnetic structures associated with the compression regions of two neighboring, interacting, corotating streams; and (3) the dynamical destruction of a small (i.e., low velocity-amplitude, short spatial-scale) stream by its overtaking of a slower moving, high-density region associated with a preceding transient flow. The evolution of these flow systems is discussed in terms of the concepts of filtering and entrainment.

Burlaga, L. F.↗

Quasiconical flowfield structure of the three-dimensional single fin interaction

A series of conical and three-dimensional computations have been performed for the swept oblique shock wave/turbulent boundary-layer interaction generated by a 20-deg sharp fin at Mach 4 and freestream Reynolds number of 2.18 x 10 exp 5 based on the incoming boundary-layer thickness. The Reynolds-averaged compressible Navier-Stokes equations are employed with turbulence incorporated using the Baldwin-Lomax and Jones-Launder models. The computed results are basically similar for both turbulence models and display general agreement with experimental data for surface pressure and surface flow direction, although underestimating the size of the primary vortex. The computed three-dimensional flowfield displays quasiconical behavior of the surface pressure, surface flow direction, and flowfield contours of static pressure, density, and Mach number over the extent of the computational domain except for an inception region near the fin leading edge. Certain features of the flowfield model are not observed in the computations, namely, a 'normal' shock near the attachment line, transonic shocklets in the expansion region, and secondary separation. The absence of these features in the computation is believed to be indirectly attributable to limitations in the turbulence models.

Knight, Doyle D.↗

Quantifying motion blur by imaging shock front propagation with broadband and narrowband X-ray sources

Time-integrated radiography using MeV Bremsstrahlung X-ray sources is the norm for imaging during system-level testing of components and structures under dynamic condition. One source of error in the analysis of the time-integrated radiography data sets stems from motion blur which smears out sharp interfaces to a greater degree with longer exposure times, which become necessary to provide sufficient signal-to-noise with low X-ray penetration of objects of interest. To quantify motion blur, a 1D shock wave through PMMA was investigated experimentally at The Dynamic Compression Sector at The Advanced Photon Source (DCS@APS) with tapered broadband and 25.46 ± 1.06 keV narrowband X-rays. Four cameras with different exposure times were used for each experiment to compare the effect that exposure time has on motion blur. In addition, our methodology to accurately simulate motion blur in terms of transmission and shape is presented and compared to our experimental results and quantified. There is a high level of agreement between the experimental and simulation results across the range of data sets investigated in this study with a percent difference range of 0.29–1.31% for the four shots. The methodology of this work serves as a steppingstone towards a physically validated model that could be used in conjunction with experimental results to deconvolve physical parameters, densities, and interfaces of interest in a way that would not be possible with experimental results alone.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Microstructural Changes and Pb Mobility During the Zircon to Reidite Transformation: Implications for Planetary Impact Chronology

Impact events modify and leave behind a complex history of rock metamorphism on terrestrial planets. Evidence for an impact event may be recorded in physical changes to minerals, such as mineral deformation and formation of high P-T polymorphs, but also in the form of chemical fingerprints, such as enhanced elemental diffusion and isotopic mixing. Here we explore laboratory shock-induced physical and chemical changes to zircon and feldspar, the former of which is of interest because its trace elements abundances and isotope ratios are used extensively in geochemistry and geochronology. To this end, a granular mixture of Bishop Tuff sanidine and Kuehl Lake zircon, both with well characterized Pb isotope compositions, was prepared and then shocked via a flat plate accelerator. The peak pressure of the experiment, as calculated by the impedance matching method, was ~24 GPa although a broader range of P-T conditions is anticipated due to starting sample porosity. Unshocked and shocked materials were characterized via scanning electron microscopy (SEM), electron backscatter diffraction (EBSD), and Raman spectroscopy. These methods show that the starting zircon material had abundant metamict regions, and the conversion of the feldspar to glass in the post-shock material. Analyses of the shocked product also yielded multiple occurrences of the high-pressure ZrSiO 4 polymorph reidite, with some domains up to 300 μm across. The possibility of U-Pb system disturbance was evaluated via laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) and secondary ion mass spectrometry (SIMS). The isotopic data reveal that disturbance of the U-Pb geochronometer in the reidite was minimal (<2% for the main U-Pb geochronometers). To better constrain the P-T conditions during the shock experiment, we complement impedance matching pressure calculations with iSALE2D impact simulations. The simulated results yield a range of P-T conditions experienced during the experiment and show that much of the sample may have reached >30 GPa, which is consistent with formation of reidite. In the recovered shocked material, we identified lamellae of reidite, some of which interlock with zircon lamellae. Reidite {112} twins were identified, which we interpret to have formed to reduce stress between the crystal structure of the host zircon and reidite. These two findings support the interpretation that shear transformation enabled the transition of zircon to reidite. The size and presence of reidite found here indicate that this phase is probably common in impact-shocked crustal rocks that experienced ~25 to ~35 GPa, especially when the target material has porosity. Additionally, shock loading of the zircon and transformation to reidite at these pressures in porous materials is unlikely to significantly disturb the U-Pb system in zircon and that the reidite inherits the primary U and Pb elemental and isotopic ratios from the zircon.

Bishop Tuff↗

Shock-driven amorphization and melting in Fe 2 O 3

We present measurements on Fe 2⁢ O 3 amorphization and melt under laser-driven shock compression up to 209(10) GPa via time-resolved in situ x-ray diffraction. At 122(3) GPa, a diffuse signal is observed indicating the presence of a noncrystalline phase. Structure factors have been extracted up to 182(6) GPa showing the presence of two well-defined peaks. A rapid change in the intensity ratio of the two peaks is identified between 145(12) and 151(12) GPa, indicative of a phase change. The noncrystalline diffuse scattering is consistent with shock amorphization of Fe 2 ⁢O 3 between 122(3) and 145(12) GPa, followed by an amorphous-to-liquid transition above 151(12) GPa. Upon release, a noncrystalline phase is observed alongside crystalline α–Fe 2⁢ O 3 . The extracted structure factor and pair distribution function of this release phase resemble those reported for Fe 2 ⁢O 3 melt at ambient pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Coronal Mass Ejections: a Summary of Recent Results

Coronal mass ejections (CMEs) have been recognized as the most energetic phenomenon in the heliosphere, deriving their energy from the stressed magnetic fields on the Sun. This paper highlights some of the recent results on CMEs obtained from the Solar and Heliospheric Observatory (SOHO) and the Solar Terrestrial Relations Observatory (STEREO) missions. The summary of the talk follows. SOHO observations revealed that the CME rate is almost a factor of two larger than previously thought and varied with the solar activity cycle in a complex way (e.g., high-latitude CMEs occurred in great abundance during the solar maximum years). CMEs were found to interact with other CMEs as well as with other large-scale structures (coronal holes), resulting in deflections and additional particle acceleration. STEREO observations have confirmed the three-dimensional nature of CMEs and the shocks surrounding them. The EUV signatures (flare arcades, corona) dimming, filament eruption, and EUV waves) associated with CMEs have become vital in the identification of solar sources from which CMEs erupt. CMEs with speeds exceeding the characteristic speeds of the corona and the interplanetary medium drive shocks, which produce type II radio bursts. The wavelength range of type II bursts depends on the CME kinetic energy: type II bursts with emission components at all wavelengths (metric to kilometric) are due to CMEs of the highest kinetic energy. Some CMEs, as fast as 1600 km/s do not produce type II bursts, while slow CMEs (400 km/s) occasionally produce type II bursts. These observations can be explained as the variation in the ambient flow speed (solar wind) and the Alfven speed. Not all CME-driven shocks produce type II bursts because either they are subcritical or do not have the appropriate geometry. The same shocks that produce type II bursts also produce solar energetic particles (SEPs), whose release near the Sun seems to be delayed with respect to the onset of type II bursts. This may indicate a subtle difference in the acceleration of the ions and 10 keV electrons needed to produce type II bursts. Surprisingly, some shocks lacking type II bursts are associated with energetic storm particle events (ESPs) pointing to the importance of electron escape from the shock for producing the radio emission. CMEs slow down or accelerate in the interplanetary medium because of the drag force, which modifies the transit time of CMEs and shocks. Halo CMEs that appear to surround the occulting disk were known before the SOHO era as occasional events. During the SOHO era, they became very prominent because of their ability to impact Earth and producing geomagnetic storms. Halo CMEs are generally more energetic than ordinary CMEs, which means they can produce severe impact on Earth's magnetosphere. Their origin close to the disk center of the Sun ensures direct impact on the magnetosphere, although their internal magnetic structure is crucial in causing storms. The solar sources of CMEs that produce SEP events at Earth, on the other hand, are generally in the western hemisphere because of the magnetic connectivity. Thus, CMEs are very interesting from the point of view of plasma physics as well as practical implications because of their space weather impact.

Gopalswamy, Nat↗

Early results from ISEE-A electric field measurements

In the solar wind and in middle latitude regions of the magnetosphere, spacecraft sheath fields obscure the ambient field under low plasma flux conditions such that valid measurements are confined to periods of moderately intense flux. Initial results show: (1) that the DC electric field is enhanced by roughly a factor of two in a narrow region at the front, increasing B, edge of the bow shock, (2) that scale lengths for large changes in E at the subsolar magnetopause are considerably shorter than scale lengths associated with the magnetic structure of the magnetopause, and (3) that the transverse distribution of B-aligned E-fields between the outer magnetosphere and ionospheric levels must be highly complex to account for the random turbulent appearance of the magnetospheric fields and the lack of corresponding time-space variations at ionospheric levels. Spike-like, non-oscillatory, fields lasting less than 0.2 seconds are occasionally seen at the bow shock and at the magnetopause and also intermittently appear in magnetosheath and plasma sheet regions under highly variable field conditions.

Heppner, J. P.↗

Effects of transient thermal shock on the microstructures and corrosion properties of a reduced activation high-entropy alloy

The devices would subject transient thermal shocks (TTS) during operation under extremely harsh conditions of nuclear fusion reactors, which inevitably exert significant impact on the microstructure and performance of structural materials. Here, a reduced activation VCrFeTa 0.2 W 0.2 high-entropy alloy (HEA) was developed by vacuum arc melting. The effects of electron beam induced TTS on its microstructure, microhardness, and corrosion properties were investigated. The results indicate that the weight fraction of each phase changes after TTS, showing a significant decrease in the content of BCC1 phase and an increase in the content of BCC2 and Laves phases. The content of BCC1 phase continues to decrease slightly with increasing the power of TTS. Besides, the microhardness of the alloy increases from ~673 HV to ~714 HV after TTS treatments. In the reduced activation HEA, TTS results in a relatively low corrosion current density of ~0.472 μA/cm 2 in 3.5 wt% NaCl solution, around one-third of the current density observed in the as-cast sample. Furthermore, the VCrFeTa 0.2 W 0.2 HEA after TTS exhibits a pitting potential of ~1.165 V SCE , which is much higher than that of the as-cast sample. The refined composite multiscale entropy method is employed to analyze the influence of TTS on current fluctuation behavior during the corrosion process. The reduced activation VCrFeTa 0.2 W 0.2 HEA exhibits excellent properties in harsh environments after TTS, thereby showing advantageous property in the field of nuclear structural materials. Furthermore, TTS is an efficient and controllable strategy for the improvement of the HEAs’ microstructures and performances.

36 MATERIALS SCIENCE↗

Characteristics of the interplanetary shocks formed by a sudden increase in the velocity of the solar wind from a coronal hole

Coronal holes are the sources of the solar wind and, according to recent YOKOH observations, may undergo rapid changes which are associated with manifestations of explosive solar activity. Rapid changes in a hole's structure will produce rapid changes in the characteristics of the wind emerging from it and, in the particular c se of a sudden increase in wind velocity, this may lead to the formation of an interplanetary shock. We discuss the characteristics of shocks formed in such a way and compare them with interplanetary observations.

Bravo, S.↗

Virtual texture analysis to investigate the deformation mechanisms in metal microstructures at the atomic scale

Understanding the deformation behavior of metallic materials at high strain rates requires the characterization of plasticity contributors such as twins, phase transformed regions, and dislocations. However, predicting the contributions from phase transformation and twinning relies on a complete understanding of the selection of variants for various loading orientations and the evolution of their volume fractions. This manuscript presents a new virtual texture (VirTex) analysis approach to characterize phase transformation and twinning variants in deformed microstructures generated using molecular dynamics (MD) simulations. Furthermore, the VirTex method involves the construction of a rotation matrix to calculate the angle/axis pairs and misorientation angles for each atom in the microstructure. Any changes in the orientation angle from angle/axis pairs and/or structure types are analyzed to determine the nucleation and evolution of variants in the microstructure. The study uses shock deformed single-crystal Fe, Ta, and Cu to analyze the variant selections for phase transformation or twinning or both in BCC and FCC systems. In addition, the VirTex analysis is able to characterize the phase transformation and twinning variants in nanocrystalline Fe and Ta microstructures. Besides characterizing variants, orientation mapping also provides an accelerated and on-the-fly approach for quantifying twin fractions in MD microstructures.

36 MATERIALS SCIENCE↗

Nature of the bonded-to-atomic transition in liquid silica to TPa pressures

First-principles calculations and analysis of the thermodynamic, structural, and electronic properties of liquid SiO 2 characterize the bonded-to-atomic transition at 0.1–1.6 TPa and 10 4 –10 5 K (1–9 eV), the high-energy-density regime relevant to understanding planetary interiors. We find strong ionic bonds that become short-lived due to high kinetics during the transition, with sensitivity of the transition temperature to pressure, and our calculated Hugoniots agree with past experimental data. Furthermore, these results reconcile previous experimental and theoretical findings by clarifying the nature of the bond dissociation process in early Earth and “rocky” (oxide) constituents of large planets.

74 ATOMIC AND MOLECULAR PHYSICS↗