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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Systematic search for high-energy gamma-ray emission from bow shocks of runaway stars

Context. It has been suggested that the bow shocks of runaway stars are sources of high-energy gamma rays (E > 100 MeV). Theoretical models predicting high-energy gamma-ray emission from these sources were followed by the first detection of non-thermal radio emission from the bow shock of BD+43°3654 and non-thermal X-ray emission from the bow shock of AE Aurigae. Aims. We perform the first systematic search for MeV and GeV emission from 27 bow shocks of runaway stars using data collected by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope (Fermi). Methods. We analysed 57 months of Fermi-LAT data at the positions of 27 bow shocks of runaway stars extracted from the Extensive stellar BOw Shock Survey catalogue (E-BOSS). A likelihood analysis was performed to search for gamma-ray emission that is not compatible with diffuse background or emission from neighbouring sources and that could be associated with the bow shocks. Results. None of the bow shock candidates is detected significantly in the Fermi-LAT energy range. We therefore present upper limits on the high-energy emission in the energy range from 100MeV to 300 GeV for 27 bow shocks of runaway stars in four energy bands. For the three cases where models of the high-energy emission are published we compare our upper limits to the modelled spectra. Our limits exclude the model predictions for ζ Ophiuchi by a factor ≈ 5.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Numerical simulation of shock initiation of Ni/Al multilayered composites

The initiation of chemical reaction in cold-rolled Ni/Al multilayered composites by shock compression is investigated numerically. A simplified approach is adopted that exploits the disparity between the reaction and shock loading timescales. The impact of shock compression is modeled using CTH simulations that yield pressure, strain, and temperature distributions within the composites due to the shock propagation. The resulting temperature distribution is then used as initial condition to simulate the evolution of the subsequent shock-induced mixing and chemical reaction. To this end, a reduced reaction model is used that expresses the local atomic mixing and heat release rates in terms of an evolution equation for a dimensionless time scale reflecting the age of the mixed layer. The computations are used to assess the effect of bilayer thickness on the reaction, as well as the impact of shock velocity and orientation with respect to the layering. Computed results indicate that initiation and evolution of the reaction are substantially affected by both the shock velocity and the bilayer thickness. In particular, at low impact velocity, Ni/Al multilayered composites with thick bilayers react completely in 100 ms while at high impact velocity and thin bilayers, reaction time was less than 100 μs. Quantitative trends for the dependence of the reaction time on the shock velocity are also determined, for different bilayer thickness and shock orientation.

36 MATERIALS SCIENCE↗
DOE OSTI2014

The link between shocks, turbulence, and magnetic reconnection in collisionless plasmas

Global hybrid (electron fluid, kinetic ions) and fully kinetic simulations of the magnetosphere have been used to show surprising interconnection between shocks, turbulence, and magnetic reconnection. In particular, collisionless shocks with their reflected ions that can get upstream before retransmission can generate previously unforeseen phenomena in the post shocked flows: (i) formation of reconnecting current sheets and magnetic islands with sizes up to tens of ion inertial length. (ii) Generation of large scale low frequency electromagnetic waves that are compressed and amplified as they cross the shock. These “wavefronts” maintain their integrity for tens of ion cyclotron times but eventually disrupt and dissipate their energy. (iii) Rippling of the shock front, which can in turn lead to formation of fast collimated jets extending to hundreds of ion inertial lengths downstream of the shock. The jets, which have high dynamical pressure, “stir” the downstream region, creating large scale disturbances such as vortices, sunward flows, and can trigger flux ropes along the magnetopause. This phenomenology closes the loop between shocks, turbulence, and magnetic reconnection in ways previously unrealized. These interconnections appear generic for the collisionless plasmas typical of space and are expected even at planar shocks, although they will also occur at curved shocks as occur at planets or around ejecta.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Dislocation structure produced by an ultrashort shock pulse

We found an ultrashort shock pulse driven by a femtosecond laser pulse on iron generates a different dislocation structure than the shock process which is on the nanosecond timescale. The ultrashort shock pulse produces a highly dense dislocation structure that varies by depth. According to transmission electron microscopy, dislocations away from the surface produce microbands via a network structure similar to a long shock process, but unlike a long shock process dislocations near the surface have limited intersections. Considering the dislocation motion during the shock process, the structure near the surface is attributed to the ultrashort shock duration. This approach using an ultrashort shock pulse will lead to understanding the whole process off shock deformation by clarifying the early stage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Compaction shock dissipation in low density granular explosive

The microstructure of granular explosives can affect dissipative heating within compaction shocks that can trigger combustion and initiate detonation. Because initiation occurs over distances that are much larger than the mean particle size, homogenized (macroscale) theories are often used to describe local thermodynamic states within and behind shocks that are regarded as the average manifestation of thermodynamic fields at the particle scale. In this paper, mesoscale modeling and simulation are used to examine how the initial packing density of granular HMX (C{sub 4}H{sub 8}N{sub 8}O{sub 8}) C{sub 4}H{sub 8}N{sub 8}O{sub 8} having a narrow particle size distribution influences dissipation within resolved, planar compaction shocks. The model tracks the evolution of thermomechanical fields within large ensembles of particles due to pore collapse. Effective shock profiles, obtained by averaging mesoscale fields over space and time, are compared with those given by an independent macroscale compaction theory that predicts the variation in effective thermomechanical fields within shocks due to an imbalance between the solid pressure and a configurational stress. Reducing packing density is shown to reduce the dissipation rate within shocks but increase the integrated dissipated work over shock rise times, which is indicative of enhanced sensitivity. In all cases, dissipated work is related to shock pressure by a density-dependent power law, and shock rise time is related to pressure by a power law having an exponent of negative one.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Mesoscale simulation of granular materials under weak shock compaction–pore size distribution effects

This research established a systematic method to generate various pore-size distributions (PSDs) and studied the effect of PSDs on the shock compaction response of granular materials using two-dimensional mesoscale simulations under identical porosity. Simulations utilized various PSDs for three particle shapes (circle, ellipse, and square). Contacting particle configurations using three PSDs, characterized by spatially uniform distributed pores to heterogeneous distributed pores, and non-contacting particle configurations under a single case of PSD were tested. The PSD of generated particle sets was characterized using coordination number, mean diameter, and bimodality coefficient as statistical metrics. Mesoscale simulations showed that regardless of the conditions of pore distributions, shock compaction of granular materials consistently demonstrates a precursor, shock compaction front, and end. However, the shock compaction velocity of contacting particles was dependent on the PSDs despite the constant initial porosity. The compaction velocity was faster in particle configurations with relatively uniform pore distributions than in heterogeneous pore distributions, which our study demonstrated can be attributed to particle rearrangement during compaction. Circular-shaped particles had high sensitivity in shock compaction response to the various PSDs. Furthermore, a contacting particle configuration tended to propagate the shock compaction wave relatively faster than particles that were in a non-contact configuration. This study established the relative importance of considering PSD as a metric over the coordination number in studies of the shock compaction response of granular materials. Further, insights are provided on the evolving shock substructure to characterize the shock compaction response of granular materials.

36 MATERIALS SCIENCE↗
DOE OSTI2013

A REVERSE SHOCK IN GRB 130427A

We present extensive radio and millimeter observations of the unusually bright GRB 130427A at z = 0.340, spanning 0.67-12 days after the burst. We combine these data with detailed multi-band UV, optical, NIR, and Swift X-ray observations and find that the broadband afterglow emission is composed of distinct reverse shock and forward shock contributions. The reverse shock emission dominates in the radio/millimeter and at ∼< 0.1 days in the UV/optical/NIR, while the forward shock emission dominates in the X-rays and at ∼> 0.1 days in the UV/optical/NIR. We further find that the optical and X-ray data require a wind circumburst environment, pointing to a massive star progenitor. Using the combined forward and reverse shock emission, we find that the parameters of the burst include an isotropic kinetic energy of E{sub K,{sub iso}} ≈ 2 × 10{sup 53} erg, a mass loss rate of M-dot ∼3×10{sup -8} M{sub ☉} yr{sup –1} (for a wind velocity of 1000 km s{sup –1}), and a Lorentz factor at the deceleration time of Γ(200 s) ≈ 130. Due to the low density and large isotropic energy, the absence of a jet break to ≈15 days places only a weak constraint on the opening angle, θ{sub j} ∼> 2.°5, and therefore a total energy of E{sub γ} + E{sub K} ∼> 1.2 × 10{sup 51} erg, similar to other gamma-ray bursts (GRBs). The reverse shock emission is detectable in this burst due to the low circumburst density, which leads to a slow cooling shock. We speculate that this property is required for the detectability of reverse shocks in radio and millimeter bands. Following on GRB 130427A as a benchmark event, observations of future GRBs with the exquisite sensitivity of the Very Large Array and ALMA, coupled with detailed modeling of the reverse and forward shock contributions, will test this hypothesis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

O{sub 2} emission toward Orion H{sub 2} peak 1 and the role of FUV-illuminated C-shocks

Molecular oxygen (O{sub 2}) has been the target of ground-based and space-borne searches for decades. Of the thousands of lines of sight surveyed, only those toward Rho Ophiuchus and Orion H{sub 2} Peak 1 have yielded detections of any statistical significance. The detection of the O{sub 2} N{sub J} = 3{sub 3}–1{sub 2} and 5{sub 4}–3{sub 4} lines at 487.249 GHz and 773.840 GHz, respectively, toward Rho Ophiuchus has been attributed to a short-lived peak in the time-dependent, cold-cloud O{sub 2} abundance, while the detection of the O{sub 2} N{sub J} = 3{sub 3}–1{sub 2}, 5{sub 4}–3{sub 4} lines, plus the 7{sub 6}–5{sub 6} line at 1120.715 GHz, toward Orion has been ascribed to time-dependent preshock physical and chemical evolution and low-velocity (12 km s{sup −1}) non-dissociative C-type shocks, both of which are fully shielded from far-ultraviolet (FUV) radiation, plus a postshock region that is exposed to an FUV field. We report a re-interpretation of the Orion O{sub 2} detection based on new C-type shock models that fully incorporate the significant effects the presence of even a weak FUV field can have on the preshock gas, shock structure, and postshock chemistry. In particular, we show that a family of solutions exists, depending on the FUV intensity, that reproduces both the observed O{sub 2} intensities and O{sub 2} line ratios. The solution in closest agreement with the shock parameters inferred for H{sub 2} Peak 1 from other gas tracers assumes a 23 km s{sup −1} shock impacting gas with a preshock density of 8 × 10{sup 4} cm{sup −3} and G{sub o} = 1, substantially different from that inferred for the fully shielded shock case. As pointed out previously, the similarity between the LSR velocity of all three O{sub 2} lines (≈11 km s{sup −1}) and recently measured H{sub 2}O 5{sub 32}–4{sub 41} maser emission at 620.701 GHz toward H{sub 2} Peak 1 suggests that the O{sub 2} emission arises behind the same shocks responsible for the maser emission, though the O{sub 2} emission is almost certainly more extended than the localized high-density maser spots. Since maser emission arises along lines of sight of low-velocity gradient, indicating shock motion largely perpendicular to our line of sight, we note that this geometry can explain not only the narrow (≲3 km s{sup −1}) observed O{sub 2} line widths despite their excitation behind a shock but also why such O{sub 2} detections are rare.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Plasma physical parameters along CME-driven shocks. II. Observation–simulation comparison

In this work, we compare the spatial distribution of the plasma parameters along the 1999 June 11 coronal mass ejection (CME)-driven shock front with the results obtained from a CME-like event simulated with the FLIPMHD3D code, based on the FLIP-MHD particle-in-cell method. The observational data are retrieved from the combination of white-light coronagraphic data (for the upstream values) and the application of the Rankine–Hugoniot equations (for the downstream values). The comparison shows a higher compression ratio X and Alfvénic Mach number M{sub A} at the shock nose, and a stronger magnetic field deflection d toward the flanks, in agreement with observations. Then, we compare the spatial distribution of M{sub A} with the profiles obtained from the solutions of the shock adiabatic equation relating M{sub A}, X, and θ{sub Bn} (the angle between the upstream magnetic field and the shock front normal) for the special cases of parallel and perpendicular shock, and with a semi-empirical expression for a generically oblique shock. The semi-empirical curve approximates the actual values of M{sub A} very well, if the effects of a non-negligible shock thickness δ{sub sh} and plasma-to magnetic pressure ratio β{sub u} are taken into account throughout the computation. Moreover, the simulated shock turns out to be supercritical at the nose and sub-critical at the flanks. Finally, we develop a new one-dimensional Lagrangian ideal MHD method based on the GrAALE code, to simulate the ion-electron temperature decoupling due to the shock transit. Two models are used, a simple solar wind model and a variable-γ model. Both produce results in agreement with observations, the second one being capable of introducing the physics responsible for the additional electron heating due to secondary effects (collisions, Alfvén waves, etc.).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Structures of strong shocks in low-density helium and neon gases

Strong shocks are essential components in many high-energy-density environments such as inertial confinement fusion implosions. However, the experimental measurements of the spatial structures of such shocks are sparse. In this paper, the soft x-ray emission of a shock front in a helium gas mixture (90% helium, 10% neon) and a pure neon gas was spatially resolved using an imaging spectrometer. We observe that the shock width in the helium mixture gas is about twice as large as in the pure neon gas. Moreover, they exhibit different precursor layers, where electron temperature greatly exceeds ion temperature, extending for more than ~350 µm with the helium gas mixture but less than 30 µm in the pure neon. At the shock front, calculations show that the electrons are strongly collisional with mean-free path two orders of magnitude shorter than the characteristic length of the shock. However, the helium ions can reach a kinetic regime as a consequence of their mean-free path being comparable to the scale of the shock. A radiation-hydrodynamic simulation demonstrates the impact of thermal conduction on the formation of the precursors with charge state, Z, playing a major role in heat flow and the precursor formation in both the helium mixture and the pure neon gases. Particle-in-cell simulations are also performed to study the ion kinetic effects on the formation of the observed precursors. A group of fast-streaming ions is observed leading the shock only in the helium gas mixture. Additionally, both effects explain the longer precursor layer in the helium shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Shock-Driven Hydrodynamic Instability Growth Near Phase Boundaries and Material Property Transitions: Final Report

Activities for this grant included: 1) Development of dynamic impact experiments to probe strength and phase transition influence on dynamic deformation, 2) development of modern strength and phase aware simulation capabilities, 3) and post-processing of experimental data with simulation and closed form analytical techniques. Two different dynamic experiments were developed to probe material strengths in solid metals (largely copper and iron in this effort). In the first experiment a flyer plate impacts a flat target with an opposite rippled surface that is partially supported by a weaker window material. Post mortem analysis of the target sample showed a strong and repeatable residual plastic deformation dependence on grain orientation. Yield strengths for strain rates near 10 5 s -1 and plastic strains near ~50% were estimated to be around 180 to 240 MPa, varying in this range with grain orientation. Unfortunately dynamic real-time measurements were difficult with this setup due to diagnostic laser scattering; hence, an additional experimental setup was developed to complement these results. In the second set of experiments a rippled surface was ablated by a controlled laser pulsed, which launched a rippled shock front to an opposite initially flat diagnostic surface that was monitored in real-time with spatially resolved velocimetry techniques, e.g., line VISAR in addition to Transient Imaging Displacement Interferometry (TIDI) displacement measurements. This setup limited the displacements at the diagnostic surface to a reasonable level for TIDI measurements (~ less than one micrometer). These experiments coupled with analytical and numerical solutions provided evidence that viscous and elastic deviatoric strength affect shock front perturbation evolution in clearly different ways. Particularly, normalized shock front perturbation amplitudes evolve with viscosity (η) and perturbation wavelength (λ) as η/λ, such that increasing viscosity (or decreasing the initial wavelength) delays the perturbation decay. Conversely our experimental data, analysis and simulations show that for materials with elastic yield strength Y the normalized shock perturbation amplitude evolves with Yλ/A 0 , which shows wavelength increases have the opposite effect as in viscous materials and perturbation decay is also dependent on initial amplitude A 0 (viscous materials are independent of this parameter). Materials where strength had clear strain rate dependence, e.g., such as a PTW material law, behaved similarly to materials with only an effective yield stress (elastic-perfectly plastic) in the shock front perturbation studies obeying a Y eff λA 0 relationship where Y eff was a constant (near ~400 MPa for Cu for strain rates around 10 6 s -1 ). Magnitude changes in strain rate would increase Y eff as would be expected from the PTW behavior, but small perturbations (typical of regions behind the shock front) near a mean had little effect. Additional work based on simulations showed that phase transformation kinetics can affect the behavior of the perturbed shock front as well as the evolution of the RM-like instability that develops due to the imprint of the perturbed shock front on the initially flat surface as the shock breaks out.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Propagation of a Strong Shock Over a Random Bed of Spherical Particles

Propagation of a strong shock through a bed of particles results in complex wave dynamics such as a reflected shock, a transmitted shock, and highly unsteady flow inside the particle bed. In this paper we present three-dimensional numerical simulations of shock propagation in air over a random bed of particles. We assume the flow is inviscid and governed by the Euler equations of gas dynamics. Simulations are carried out by varying the volume fraction of the particle bed at a fixed shock Mach number. We compute the unsteady inviscid streamwise and transverse drag coefficients as a function of time for each particle in the random bed as a function of volume fraction. We show that (i) there are significant variations in the peak drag for the particles in the bed, (ii) the mean peak drag as a function of streamwise distance through the bed decreases with a slope that increases as the volume fraction increases, and (iii) the deviation from the mean peak drag does not correlate with local volume fraction. We also present the local Mach number and pressure contours for the different volume fractions to explain the various observed complex physical mechanisms occurring during the shock-particle interactions. Since the shock interaction with the random bed of particles leads to transmitted and reflected waves, we compute the average flow properties to characterize the strength of the transmitted and reflected shock waves and quantify the energy dissipation inside the particle bed. Finally, to better understand the complex wave dynamics in a random bed, we consider a simpler approximation of a planar shock propagating in a duct with a sudden area change. We obtain Riemann solutions to this problem, which are used to compare with fully resolved numerical simulations.

42 ENGINEERING↗
DOE OSTI2020

Formation and Diagnostics of Collisionless Shocks in Laser-Produced Plasmas (Final Report)

Collisions of supersonic flows occur frequently in astrophysics, and the resulting shock waves are responsible for the properties of many astrophysical phenomena, such as supernova remnants, Gamma Ray Bursts and jets from Active Galactic Nuclei. Because of the low density of astrophysical plasmas, the mean free path due to Coulomb collisions is typically very large; therefore, most shock waves in astrophysics are “collisionless” and form due to plasma instabilities and self-generated magnetic fields. Laboratory experiments at DOE high energy density (HED) facilities can achieve the conditions necessary for the formation of collisionless shocks and will provide a unique avenue for studying the nonlinear physics of shock waves. This project investigates the parameters and experimental configurations needed for the creation of collisionless shocks in the laboratory experiments. This study is being performed with state of the art 3D kinetic plasma simulations. Recent advances in ab-initio particle-in-cell (PIC) simulations of astrophysical shocks enable us to simulate shock formation on timescales and in geometries relevant to the laboratory experiments. The main goal of the proposal is to ascertain which signatures of shock physics can be extracted from the realistic experiments, accounting for limitations imposed by finite spatial and time scales in the lab. To aid the interpretation of current and future experimental data, PIC simulation data have been analyzed using the same diagnostic techniques as used in the laboratory, such as proton radiography and neutronics, to unambiguously predict the experimental signatures of relevant collisionless shock phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Shock-Induced Solid-State Reactions in Powders: An Experimentally-Based Reassessment

We have performed a series of ten planar impact experiments on six different materials, including certain reactive powders and two inert materials, using Sandia's 89 - mm powder gun at the STAR facility. Time - resolved particle - velocity histories were determined during each of the experiments from one or more VISAR measurements. We have analyzed the results of these measurements 1) by using jump conditions to determine shock and first reshock states and 2) by comparing measured particle velocity histories to synthetic histories predicted by one - dimensional computational analyses using the CTH shock physics code with various models for inert and reactive materials . These comparisons are consistent with the conclusion for these particular reactive powders, that for the duration of shock loading either 1) there is insignificant reaction or 2) the products of any reaction are indistinguishable from the reactants under the experimental conditions. Shock and reshock states were extracted for shock pressures between 5 and 40 GPa. Densities were at or greater than the theoretical maximum zero - pressure density of the starting mixture. This result would be expected if there were no reaction or negligible reaction for the first two shock states. Two experiments were performed on one reactive powder in a "ring - down " geometry to look for evidence of vapor production on pressure release. In both cases, the measured velocity continued to increase slowly over a period of microseconds for the du ration of the experiment. This observation suggests that vapor is produced along the release path, but information about the mechanism for vapor production cannot be extracted from these data. While it is possible that vapor is produced by a shock - induced reaction involving more than one of the original constituents, a simpler interpretation is that the vapor is made up of products of shock - induced decomposition reactions and/or simple vaporization of the constituents as would be expected to take place under the conditions of these experiments. Other sources of vapor could be water adsorbed on grain surfaces and air originally in the voids. Thus it is not necessary to invoke significant recombination reactions to explain the data. However, in the absence of ring - down control experiments, the possibility remains open. These conclusions are different from those of previous workers, but reassessment of a subset of the earlier data yields results consistent with the present work, i.e., the shock compression data do not provide evidence for strong exothermic reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multi-Axis Resonant Plate Shock Testing Evaluation and Test Specification Development

Resonant plate testing is a shock test method that is frequently used to simulate pyroshock events in the laboratory. Recently, it was discovered that if the unit under test is installed at an off-center location, a tri-axial accelerometer would record a shock response in three directions and the resulting shock response spectra implied that the test may have qualified the component in three directions simultaneously. The purpose of this research project was to evaluate this idea of multi-axis shock testing to determine if it was truly a multi-axis shock environment and if such a test could be used as an equivalent component qualification test. A study was conducted using generic, additively manufactured components tested on a resonant plate, along with an investigation of plate motion to evaluate the component response to off- center plate excitation. The data obtained here along with the analytical simulations performed indicate that off-center resonant plate tests are actually not three-axis shock tests, but rather single axis shocks at an arbitrary angle dictated by the location of the unit under test on the plate. This conclusion is supported by the fact that only one vectored shock input is provided to the component in a resonant plate test. Thus, the output response is a coupled response of the transverse plate vibration and the rotational motion of the component on the plate. Additionally, a multi-axis shock test defined by three single axis test environments always results in a significant component over-test in one direction.

42 ENGINEERING↗
DOE OSTI2025

Numerical Treatment of Shock-induced Nuclear Burning in Double Detonation Type Ia Supernovae

We present a benchmark problem to assess the treatment of shock-induced nuclear burning in the context of double detonation Type Ia supernovae. In a stratified white dwarf model, we implement a shock-detection criterion that suppresses burning in zones characterized by compression and significant pressure gradients, controlled by a tunable parameter, f shock . One-dimensional simulations, using the open-source Castro suite, were conducted across three treatments—burning fully enabled, and burning suppressed with f shock = 2/3 and f shock = 1—across three spatial resolutions (5.0, 2.5, and 0.3125 km). At the finest resolution, the burning-enabled and f shock = 1 models converge, while the f shock = 2/3 front continues to show slight offset behavior. Since most simulations are carried out at much lower resolutions, our tests support the idea that burning in shocks should always be disabled in practice. We also observe that the behavior of lower-resolution simulations remains extremely sensitive to the choice of f shock .

GPU computing↗
NASA NTRS1974

On the local time dependence of the bow shock wave structure

In the first 6 months after its launch, Ogo 3 crossed the earth's bow shock over 500 times. From this group, a set of 494 shock crossings were chosen for analysis. These crossings, as they were recorded by the UCLA/JPL search coil magnetometer, were scanned and classified according to the nature of the plasma waves detected near the shock. More than 85% of the shocks detected fell into a single category showing the predominance of two independent wave trains near the shock, the higher frequency appearing upstream and the lower downstream. The other 15%, which constitute an upper limit, appear to be composed of shocks dominated by a single wave pattern and of chaotic shocks showing no orderly progression of wave frequencies as the shock was penetrated. This division of wave pattern was found to occur at all local times, that is, in all regions where the satellite penetrated the shock.

Olson, J. V.↗
NASA NTRS1977

Pioneer 10, 11 observations of evolving solar wind streams and shocks beyond 1 AU

Observations between 1 and 5 AU by Pioneers 10 and 11 have led to the identification of large numbers of interplanetary shocks. Both forward and reverse shocks, which begin to develop beyond 1.5 AU and which frequently appear as shock pairs, are found to accompany solar wind streams. The number of forward shocks continues to increase out to at least 5 AU. Reverse shocks are seen less often than forward shocks and, in some instances, disappear at larger distances. There is evidence that the shocks are corotating in the solar frame, as anticipated theoretically. The evolution of solar wind streams beyond 1 AU is profoundly affected by the shocks. A thick interaction region, with large enhancements in density, temperature, field strength and fluctuation level, forms in the region originally characterized by a positive velocity gradient. The solar wind and magnetic field properties adjacent to, and within, the interaction regions have been studied to determine their qualitative behavior and characteristic changes with distance. Several interplanetary shocks generated by solar flares have also been identified and analyzed.

Smith, E. J.↗