Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHOCK”

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

Creating Calibration Curves to Determine Shock Pressure in Clinopyroxene

Impact cratering is an important geological process that occurs on every rocky body in the solar system. It alters the texture and mineralogy of rocks via shock metamorphism. The peak shock pressures experienced by a rock are traditionally evaluated using qualitative optical methods however, quantitative methods do exist. One such method was developed by Uchizono et al., who used X-ray Diffraction (XRD) to measure lattice strain (ε) in several artificially shocked olivine grains using XRD peak broadening as a function of tan θ, where θ is the diffraction angle. They plotted the ε values against the known peak shock pressures experienced by the olivine grains. Using this calibration curve, the precise shock pressure experienced by a grain of olivine can be determined using its measured ε value. Another method was developed by McCausland et al. and Izawa et al., who used in situ XRD to measure strain-related mosaicity (SRM) of olivine in several ordinary chondrites and enstatite in enstatite chondrites, respectively. They plotted these results against the shock stage estimates for these meteorites. Using these plots, meteorites can be assigned to shock stage bins by measuring the SRM of olivine and/or enstatite. Both methods are useful for evaluating shock metamorphism, however, they have limitations. Uchizono et al.’s calibration curve has been successfully applied to martian meteorites, however it can only be applied to olivine-bearing rocks. McCausland et al.’s and Izawa et al.’s SRM method is uncalibrated and is limited to binning meteorites by shock stage. This work aims to expand on both methods by creating calibration curves for clinopyroxene (CPX): one for ε, similar to Uchizono et al.’s calibration curve for olivine, and one for SRM. This will extend the application of shock calibration methods to a greater variety of rock types. Preliminary results are presented herein.

Jenkins, Laura E.↗

High-Order Shock Fitting with Finite Element Methods

A moving-grid, shock-tracking, finite element method has been implemented that can achieve high-order accuracy for flow simulations with shocks. In this approach, element edges in the computational mesh are fitted to the shock front and moved with the shock throughout the simulation. The Euler equations are solved on the moving mesh in an arbitrary Lagrangian-Eulerian framework. Three different methods for specifying the shock motion have been tested,and the order of accuracy of the resulting solutions was verified. It was found that one of the methods is applicable only for problems with a constant upstream flow state, whereas the other two methods converge to the expected solution for general shock problems. When using a finite element method with polynomial degree p, the two general shock motion methods converged with the expected order of accuracy of (p+1). Finally, accurate solutions are obtained for supersonic nozzle and blunt body bow shock problems. Based on the results, it is concluded that the shock tracking method can preserve high-order solution accuracy in the vicinity of a shock wave.

Luke D'Aquila↗

Climate Variability and Simultaneous Breadbasket Yield Shocks as Observed in Long-Term Yield Records

That climate variability and change can potentially force multiple simultaneous breadbasket crop yield shocks has been established. But research quantifying the mechanisms behind such simultaneous shocks has been constrained by short records of crop yields. Here we compile a dataset of subnational crop yields in 25 countries dating back to 1900 to study the frequency and trends in multiple breadbasket yield shocks and how large-scale climate anomalies on interannual timescales have affected multiple breadbasket yield shocks over the last century. We find that major simultaneous breadbasket yield shocks have occurred in at least three, four, or five of nine breadbaskets 10.3%, 2.3% and 1.1% of the time for maize and 18.4%, 4.6% and 2.3% of the time for wheat. Furthermore, we find that multiple breadbasket yield shocks decreased in frequency even as those breadbaskets experience increasingly frequent climate-related shocks. For both maize and wheat breadbaskets, there were fewer simultaneous yield shocks during the 1975–2017 time period as compared to 1931–1975. Finally, we find that interannual modes of climate variability - such as the El Niño Southern Oscillation (ENSO), the Indian Ocean Dipole (IOD), and the North Atlantic Oscillation (NAO) - have all affected the relative probability of simultaneous yield shocks in pairs of breadbaskets by up to 20–40% in both maize and wheat breadbaskets. While past literature has focused on the effects of ENSO, we find that at the global scale the NAO affects the overall number of wheat yield shocks most strongly despite only affecting northern hemisphere breadbaskets.

agriculture↗

Species Separation and Other Kinetic Effects at Shock Wave Fronts in Plasma Mixtures Assessed using Molecular Dynamics

The propagation of a shock through a binary plasma mixture of hydrogen isotopes is simulated with molecular dynamics (MD) and a multispecies continuum hydrodynamics model. Separation of the two species due to their mass difference is observed at the shock front in both simulations. We analyze the shock structure in detail to determine which aspects of the wave propagation and shock profile obey single-fluid hydrodynamics and which require multi-component hydrodynamics with Fickian and barodiffusion. In order to reduce possible sources of uncertainty due to transport models used in the continuum simulations, we have performed separate, equilibrium MD simulations which verify the accuracy of the recently proposed Kinetic Molecular Dynamics (KMD) viscosity model, even at the relatively cold pre-shock temperature of 1 eV. We find the post-shock temperature and density are generally consistent between continuum and atomic-scale calculations, but the MD simulations predict a significantly larger thickness of the leading light-ion shock front, with some associated differences in the species separation. We analyze the shape and anisotropy of the MD velocity distributions in detail and compute the electric field at the shock interface from both the MD and continuum simulations. We find that the continuum hydrodynamics calculation captures the major features of the electric fi eld through the electron pressure term in the barodiffusion model and conclude that the differences observed in the shock structures from the two simulations arise largely from non-Maxwellian effects in the longitudinal velocities (the velocities in the direction that the shock propagates) that are not captured in the hydrodynamic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advanced and Exploratory Shock Sensing Mechanisms.

Mechanical component response to shock environments must be predictable in order to ensure reliability and safety. Whether the shock input results from accidental drops during transportation to projectile impact scenarios, the system must irreversibly transition into a safe state that is incapable of triggering the component . With this critical need in mind, the 2017 Nuclear Weapons Summer Product Realization Institute (NW SPRINT) program objective sought the design of a passive shock failsafe with emphasis on additively manufactured (AM) components. Team Advanced and Exploratory (A&E) responded to the challenge by designing and delivering multiple passive shock sensing mech anisms that activate within a prescribed mechanical shock threshold. These AM failsafe designs were tuned and validated using analytical and computational techniques including the shock response spectrum (SRS) and finite element analysis (FEA). After rapid prototyping, the devices experienced physical shock tests conducted on Sandia drop tables to experimentally verify performance. Keywords: Additive manufacturing, dynamic system, failsafe, finite element analysis, mechanical shock, NW SPRINT, shock respon se spectrum

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Kinetic simulations comparing quasi-parallel and quasi-perpendicular piston-driven collisionless shock dynamics in magnetized laboratory plasmas

Magnetized collisionless shocks are common in astrophysical systems, and scaled versions can be created in laboratory experiments by utilizing laser-driven piston plasmas to create these shocks in a magnetized background plasma. A key parameter for these experiments is the angle θB between the shock propagation direction and the background magnetic field. We performed quasi-1D piston-driven shock simulations to explore shock formation, evolution, and key observables relevant to laboratory experiments for a range of shock angles between θB=90° to θB=30°. Our results show that the spatial and temporal scales of shock formation for all angles considered are similar when expressed in terms of the perpendicular component of the magnetic field. In a steady state, ion and electron temperatures become more isotropic, and the electron-to-ion temperature ratio is higher for smaller θB. At θB=30°, ion heating parallel to the magnetic field becomes dominant, associated with more ions being reflected at one discontinuity and subsequently trapped by the next discontinuity due to shock reformation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-order shock-fitted detonation propagation in high explosives

In this paper, a highly accurate numerical shock and material interface fitting scheme composed of fifth-order spatial and third- or fifth-order temporal discretizations is applied to the two-dimensional reactive Euler equations in both slab and axisymmetric geometries. High rates of convergence are not typically possible with shock-capturing methods as the Taylor series analysis breaks down in the vicinity of discontinuities. Furthermore, for typical high explosive (HE) simulations, the effects of material interfaces at the charge boundary can also cause significant computational errors. Fitting a computational boundary to both the shock front and material interface (i.e. streamline) alleviates the computational errors associated with captured shocks and thus opens up the possibility of high rates of convergence for multi-dimensional shock and detonation flows. Several verification tests, including a Sedov blast wave, a Zel'dovich–von Neumann–Döring (ZND) detonation wave, and Taylor–Maccoll supersonic flow over a cone, are utilized to demonstrate high rates of convergence to nontrivial shock and reaction flows. Comparisons to previously published shock-capturing multi-dimensional detonations in a polytropic fluid with a constant adiabatic exponent (PF-CAE) are made, demonstrating significantly lower computational error for the present shock and material interface fitting method. For an error on the order of 10 m / s , which is similar to that observed in experiments, shock-fitting offers a computational savings on the order of 1000. In addition, the behavior of the detonation phase speed is examined for several slab widths to evaluate the detonation performance of PBX 9501 while utilizing the Wescott–Stewart–Davis (WSD) model, which is commonly used in HE modeling. It is found that the thickness effect curve resulting from this equation of state and reaction model using published values is dramatically more steep than observed in recent experiments. Utilizing the present fitting strategy, in conjunction with a nonlinear optimizer, a new set of reaction rate parameters improves the correlation of the model to experimental results. In conclusion, this new model is tested against two dimensional slabs as a validation test.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗

Evaluation of the relative importance of preheat from hohlraum x rays and a radiative shock on a low-density foam

Indirectly driven shock-tube experiments were performed on the Omega Laser Facility to evaluate the relative importance of hohlraum x ray and radiative shock preheat sources on a low-density foam. X rays emitted from the hohlraum and a subsequent shock wave are channeled into a low-density foam sample, which contains a plastic preheat-witness disk placed downstream of the foam. Simultaneous radiographic measurements of the shock speed in the foam and the expansion rate of the witness disk due to preheat allow for the observation of effects from the hohlraum's low-energy and high-energy x-ray spectrum. Here we showed, from simulations, that low-energy x rays from the hohlraum are preferentially absorbed near the ablator surface (where the hohlraum and the shock tube meet), while higher-energy x rays largely pass through the ablator and foam and are volumetrically absorbed by the witness disk. Reproducing the experimentally measured shock speed and expansion of the witness disk simultaneously, we extracted the temperature evolution of preheated foam from the simulation and evaluated the relative importance of preheat sources on a low-density foam from hohlraum x-ray radiation and radiative shock. We found that radiation from the shock front was more effective at preheating the low-density foam than the high-energy x rays from the hohlraum. This shock-tube preheat experiment is important for understanding the results of the MARBLE experiments at the National Ignition Facility because initial conditions of foam-filled MARBLE capsules are sensitive to preheat.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

THE MAXIMUM ENERGY OF ACCELERATED PARTICLES IN RELATIVISTIC COLLISIONLESS SHOCKS

The afterglow emission from gamma-ray bursts (GRBs) is usually interpreted as synchrotron radiation from electrons accelerated at the GRB external shock that propagates with relativistic velocities into the magnetized interstellar medium. By means of multi-dimensional particle-in-cell simulations, we investigate the acceleration performance of weakly magnetized relativistic shocks, in the magnetization range 0 {approx}< {sigma} {approx}< 10{sup -1}. The pre-shock magnetic field is orthogonal to the flow, as generically expected for relativistic shocks. We find that relativistic perpendicular shocks propagating in electron-positron plasmas are efficient particle accelerators if the magnetization is {sigma} {approx}< 10{sup -3}. For electron-ion plasmas, the transition to efficient acceleration occurs for {sigma} {approx}< 3 Multiplication-Sign 10{sup -5}. Here, the acceleration process proceeds similarly for the two species, since the electrons enter the shock nearly in equipartition with the ions, as a result of strong pre-heating in the self-generated upstream turbulence. In both electron-positron and electron-ion shocks, we find that the maximum energy of the accelerated particles scales in time as {epsilon}{sub max}{proportional_to}t {sup 1/2}. This scaling is shallower than the so-called (and commonly assumed) Bohm limit {epsilon}{sub max}{proportional_to}t, and it naturally results from the small-scale nature of the Weibel turbulence generated in the shock layer. In magnetized plasmas, the energy of the accelerated particles increases until it reaches a saturation value {epsilon}{sub sat}/{gamma}{sub 0} m{sub i}c {sup 2} {approx} {sigma}{sup -1/4}, where {gamma}{sub 0} m{sub i}c {sup 2} is the mean energy per particle in the upstream bulk flow. Further energization is prevented by the fact that the self-generated turbulence is confined within a finite region of thickness {proportional_to}{sigma}{sup -1/2} around the shock. Our results can provide physically grounded inputs for models of non-thermal emission from a variety of astrophysical sources, with particular relevance to GRB afterglows.

79 ASTRONOMY AND ASTROPHYSICS↗

New evidence for efficient collisionless heating of electrons at the reverse shock of a young supernova remnant

Although collisionless shocks are ubiquitous in astrophysics, certain key aspects of them are not well understood. In particular, the process known as collisionless electron heating, whereby electrons are rapidly energized at the shock front, is one of the main open issues in shock physics. Here, we present the first clear evidence for efficient collisionless electron heating at the reverse shock of Tycho's supernova remnant (SNR), revealed by Fe K diagnostics using high-quality X-ray data obtained by the Suzaku satellite. We detect Kβ (3p → 1s) fluorescence emission from low-ionization Fe ejecta excited by energetic thermal electrons at the reverse shock front, which peaks at a smaller radius than Fe Kα (2p → 1s) emission dominated by a relatively highly ionized component. Comparisons with our hydrodynamical simulations imply instantaneous electron heating to a temperature 1000 times higher than expected from Coulomb collisions alone. The unique environment of the reverse shock, which is propagating with a high Mach number into rarefied ejecta with a low magnetic field strength, puts strong constraints on the physical mechanism responsible for this heating and favors a cross-shock potential created by charge deflection at the shock front. Our sensitive observation also reveals that the reverse shock radius of this SNR is about 10% smaller than the previous measurement using the Fe Kα morphology from the Chandra observations. Since strong Fe Kβ fluorescence is expected only from low-ionization plasma where Fe ions still have many 3p electrons, this feature is key to diagnosing the plasma state and distribution of the immediate postshock ejecta in a young SNR.

79 ASTRONOMY AND ASTROPHYSICS↗

A merger shock in A2034

We present a 250 ks Chandra observation of the cluster merger A2034 with the aim of understanding the nature of a sharp edge previously characterized as a cold front. The new data reveal that the edge is coherent over a larger opening angle and is significantly more bow-shock-shaped than previously thought. Within ∼27° about the axis of symmetry of the edge, the density, temperature, and pressure drop abruptly by factors of 1.83{sub −0.08}{sup +0.09}, 1.85{sub −0.41}{sup +0.41}, and 3.4{sub −0.7}{sup +0.8}, respectively. This is inconsistent with the pressure equilibrium expected of a cold front and we conclude that the edge is a shock front. We measure a Mach number M=1.59{sub −0.07}{sup +0.06} and corresponding shock velocity v {sub shock} ≅ 2057 km s{sup –1}. Using spectra collected at the MMT with the Hectospec multi-object spectrograph, we identify 328 spectroscopically confirmed cluster members. Significantly, we find a local peak in the projected galaxy density associated with a bright cluster galaxy that is located just ahead of the nose of the shock. The data are consistent with a merger viewed within ∼23° of the plane of the sky. The merging subclusters are now moving apart along a north-south axis approximately 0.3 Gyr after a small impact parameter core passage. The gas core of the secondary subcluster, which was driving the shock, appears to have been disrupted by the merger. Without a driving 'piston,' we speculate that the shock is dying. Finally, we propose that the diffuse radio emission near the shock is due to the revival of pre-existing radio plasma that has been overrun by the shock.

79 ASTRONOMY AND ASTROPHYSICS↗

PROPERTIES OF A CORONAL SHOCK WAVE AS A DRIVER OF EARLY SEP ACCELERATION

Coronal mass ejectmons (CMEs) are thought to drive collisionless shocks in the solar corona, which in turn have been shown to be capable of accelerating solar energetic particles (SEPs) in minutes. It has been notoriously difficult to extract information about energetic particle spectra in the corona, owing to a lack of in situ measurements. It is possible, however, to combine remote observations with data-driven models in order to deduce coronal shock properties relevant to the local acceleration of SEPs and their heliospheric connectivity to near-Earth space. We present such novel analysis applied to the 2011 May 11 CME event on the western solar limb, focusing on the evolution of the eruption-driven, dome-like shock wave observed by the Atmospheric Imaging Assembly (AIA) EUV telescopes on board the Solar Dynamics Observatory spacecraft. We analyze the shock evolution and estimate its strength using emission measure modeling. We apply a new method combining a geometric model of the shock front with a potential field source surface model to estimate time-dependent field-to-shock angles and heliospheric connectivity during shock passage in the low corona. We find that the shock was weak, with an initial speed of ∼450 km s{sup –1}. It was initially mostly quasi-parallel, but a significant portion of it turned quasi-perpendicular later in the event. There was good magnetic connectivity to near-Earth space toward the end of the event as observed by the AIA instrument. The methods used in this analysis hold a significant potential for early characterization of coronal shock waves and forecasting of SEP spectra based on remote observations.

79 ASTRONOMY AND ASTROPHYSICS↗

SHORT ACCELERATION TIMES FROM SUPERDIFFUSIVE SHOCK ACCELERATION IN THE HELIOSPHERE

The analysis of time profiles of particles accelerated at interplanetary shocks allows particle transport properties to be inferred. The frequently observed power-law decay upstream, indeed, implies a superdiffusive particle transport when the level of magnetic field variance does not change as the time interval from the shock front increases. In this context, a superdiffusive shock acceleration (SSA) theory has been developed, allowing us to make predictions of the acceleration times. In this work we estimate for a number of interplanetary shocks, including the solar wind termination shock, the acceleration times for energetic protons in the framework of SSA and we compare the results with the acceleration times predicted by standard diffusive shock acceleration. The acceleration times due to SSA are found to be much shorter than in the classical model, and also shorter than the interplanetary shock lifetimes. This decrease of the acceleration times is due to the scale-free nature of the particle displacements in the framework of superdiffusion. Indeed, very long displacements are possible, increasing the probability for particles far from the front of the shock to return, and short displacements have a high probability of occurrence, increasing the chances for particles close to the front to cross the shock many times.

79 ASTRONOMY AND ASTROPHYSICS↗

The synchrotron maser emission from relativistic shocks in Fast Radio Bursts: 1D PIC simulations of cold pair plasmas

The emission process of Fast Radio Bursts (FRBs) remains unknown. We look into whether the synchrotron maser emission from relativistic shocks in a magnetar wind can explain the observed FRB properties. We perform particle-in-cell (PIC) simulations of perpendicular shocks in cold pair plasmas, checking our results for consistency among three PIC codes. We confirm that a linearly polarized X-mode wave is self-consistently generated by the shock and propagates back upstream as a precursor wave. We find that at magnetizations σ ≳ 1 (i.e. ratio of Poynting flux to particle energy flux of the pre-shock flow) the shock converts a fraction $f_\xi ^{\prime } \approx 7 \times 10^{-4}/\sigma ^2$ of the total incoming energy into the precursor wave, as measured in the shock frame. The wave spectrum is narrow-band (fractional width ≲1-3), with apparent but not dominant line-like features as many resonances concurrently contribute. The peak frequency in the pre-shock (observer) frame is $\omega ^{\prime \prime }_{\rm peak} \approx 3 \gamma _{\rm s | u} \omega _{\rm p}$, where γs|u is the shock Lorentz factor in the upstream frame and ωp the plasma frequency. At σ ≳ 1, where our estimated $\omega ^{\prime \prime }_{\rm peak}$ differs from previous works, the shock structure presents two solitons separated by a cavity, and the peak frequency corresponds to an eigenmode of the cavity. Our results provide physically grounded inputs for FRB emission models within the magnetar scenario.

79 ASTRONOMY AND ASTROPHYSICS↗

SIMULATION OF ENERGETIC PARTICLE TRANSPORT AND ACCELERATION AT SHOCK WAVES IN A FOCUSED TRANSPORT MODEL: IMPLICATIONS FOR MIXED SOLAR PARTICLE EVENTS

We use numerical solutions of the focused transport equation obtained by an implicit stochastic differential equation scheme to study the evolution of the pitch-angle dependent distribution function of protons in the vicinity of shock waves. For a planar stationary parallel shock, the effects of anisotropic distribution functions, pitch-angle dependent spatial diffusion, and first-order Fermi acceleration at the shock are examined, including the timescales on which the energy spectrum approaches the predictions of diffusive shock acceleration theory. We then consider the case that a flare-accelerated population of ions is released close to the Sun simultaneously with a traveling interplanetary shock for which we assume a simplified geometry. We investigate the consequences of adiabatic focusing in the diverging magnetic field on the particle transport at the shock, and of the competing effects of acceleration at the shock and adiabatic energy losses in the expanding solar wind. We analyze the resulting intensities, anisotropies, and energy spectra as a function of time and find that our simulations can naturally reproduce the morphologies of so-called mixed particle events in which sometimes the prompt and sometimes the shock component is more prominent, by assuming parameter values which are typically observed for scattering mean free paths of ions in the inner heliosphere and energy spectra of the flare particles which are injected simultaneously with the release of the shock.

79 ASTRONOMY AND ASTROPHYSICS↗

SUPRATHERMAL ELECTRONS AT SATURN'S BOW SHOCK

The leading explanation for the origin of galactic cosmic rays is particle acceleration at the shocks surrounding young supernova remnants (SNRs), although crucial aspects of the acceleration process are unclear. The similar collisionless plasma shocks frequently encountered by spacecraft in the solar wind are generally far weaker (lower Mach number) than these SNR shocks. However, the Cassini spacecraft has shown that the shock standing in the solar wind sunward of Saturn (Saturn's bow shock) can occasionally reach this high-Mach number astrophysical regime. In this regime Cassini has provided the first in situ evidence for electron acceleration under quasi-parallel upstream magnetic conditions. Here we present the full picture of suprathermal electrons at Saturn's bow shock revealed by Cassini . The downstream thermal electron distribution is resolved in all data taken by the low-energy electron detector (CAPS-ELS, <28 keV) during shock crossings, but the higher energy channels were at (or close to) background. The high-energy electron detector (MIMI-LEMMS, >18 keV) measured a suprathermal electron signature at 31 of 508 crossings, where typically only the lowest energy channels (<100 keV) were above background. We show that these results are consistent with the theory in which the “injection” of thermal electrons into an acceleration process involves interaction with whistler waves at the shock front, and becomes possible for all upstream magnetic field orientations at high Mach numbers like those of the strong shocks around young SNRs. A future dedicated study will analyze the rare crossings with evidence for relativistic electrons (up to ∼1 MeV).

79 ASTRONOMY AND ASTROPHYSICS↗

First detection of a pulsar bow shock nebula in far-UV: PSR J0437-4715

Pulsars traveling at supersonic speeds are often accompanied by cometary bow shocks seen in Hα. We report on the first detection of a pulsar bow shock in the far-ultraviolet (FUV). We detected it in FUV images of the nearest millisecond pulsar J0437−4715 obtained with the Hubble Space Telescope. The images reveal a bow-like structure positionally coincident with part of the previously detected Hα bow shock, with an apex at 10″ ahead of the moving pulsar. Its FUV luminosity, L(1250--2000 A-ring )≈5 × 10{sup 28} erg s{sup −1}, exceeds the Hα luminosity from the same area by a factor of 10. The FUV emission could be produced by the shocked interstellar medium matter or, less likely, by relativistic pulsar wind electrons confined by strong magnetic field fluctuations in the bow shock. In addition, in the FUV images we found a puzzling extended (≃3″ in size) structure overlapping with the limb of the bow shock. If related to the bow shock, it could be produced by an inhomogeneity in the ambient medium or an instability in the bow shock. We also report on a previously undetected X-ray emission extending for about 5″ ahead of the pulsar, possibly a pulsar wind nebula created by shocked pulsar wind, with a luminosity L(0.5–8 keV) ∼ 3 × 10{sup 28} erg s{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗